2.7: Solution of the Dirac equation for a free particle
- Page ID
- 20886
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\dsum}{\displaystyle\sum\limits} \)
\( \newcommand{\dint}{\displaystyle\int\limits} \)
\( \newcommand{\dlim}{\displaystyle\lim\limits} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\(\newcommand{\longvect}{\overrightarrow}\)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\ket}[1]{\left| #1 \right>}\)
\(\newcommand{\bra}[1]{\left< #1 \right|}\)
\(\newcommand{\braket}[2]{\left< #1 \vphantom{#2} \right| \left. #2 \vphantom{#1} \right>}\)
\(\newcommand{\braopket}[3]{\left< #1 \vphantom{#2}\vphantom{#3} \right| #2 \vphantom{#1}\vphantom{#3} \left| #3 \vphantom{#1}\vphantom{#2} \right>}\)
\(\newcommand{\qmvec}[1]{\mathbf{\vec{#1}}}\)
\(\newcommand{\op}[1]{\hat{\mathbf{#1}}}\)
\(\newcommand{\expect}[1]{\langle #1 \rangle}\)
\(\newcommand{\dfn}[1]{\emph{\textbf{#1}}}\)
The Dirac Hamiltonian takes the form
where

Using \({\bf P}= (\hbar/i)\nabla\), in the coordinate basis, the Dirac equation for a free particle reads
Since the operator on the left side is a 4\(\times\)4 matrix, the wave function \(\Psi({\bf r},t)\) is actually a four-component vector of functions of \({\bf r}\) and \(t\):

which is called a four-component Dirac spinor. In order to generate an eigenvalue problem, we look for a solution of the form
which, when substituted into the Dirac equation gives the eigenvalue equation
Note that, since \(H\) is only a function of \({\bf P}\), then \([{\bf P},H]=0\) so that the eigenvalues \({\bf p}\) of \({\bf P}\) can be used to characterize the states. In particular, we look for free-particle (plane-wave) solutions of the form:
where \(u_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/p[8]/span, line 1, column 1
Since the matrix on the left is expressible in terms of 2\(\times\)2 blocks, we look for \(u_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/p[10]/span, line 1, column 1

Therefore, writing the equation in matrix form, we find

or

which yields two equations
| \(\displaystyle \left(E-mc^2\right)\phi_ (click for details) - c\stackrel{\rightarrow}{\sigma}\cdot{\bf p}\chi_Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[12]/table/tbody/tr[1]/td[1]/span[1], line 1, column 1
(click for details) \)Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[12]/table/tbody/tr[1]/td[1]/span[2], line 1, column 1
| \(\textstyle =\) | \(\displaystyle 0\) | |
| \(\displaystyle - c\stackrel{\rightarrow}{\sigma}\cdot{\bf p}\phi_ (click for details) + \left(E+mc^2\right)\chi_Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[12]/table/tbody/tr[2]/td[1]/span[1], line 1, column 1
(click for details) \)Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[12]/table/tbody/tr[2]/td[1]/span[2], line 1, column 1
| \(\textstyle =\) | \(\displaystyle 0\) |
From the second equation:
Note, one could also solve the first for \(\phi_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/p[15]/span, line 1, column 1
Using the first of these, then a single equation for \(\phi_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/p[16]/span, line 1, column 1
However,
Hence, we have the condition
Since \(\phi_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/p[19]/span, line 1, column 1
We see that the eigenvalues can be positive or negative. A plot of the energy levels is shown below:
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/span[1], line 1, column 1
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/span[2], line 1, column 1
We will show that for \(E>0\), an appropriate solution is to take
If this is the case, then
However,

so that

so that the full solution \(u_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[5]/span, line 1, column 1

Note that when \({\bf p}=0\), the third and fourth components of \(u_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[6]/span, line 1, column 1

which are both forward propagating solutions. These correspond to particle solutions, in particular, a spin-1/2 particle propagating forward in time with an energy equal to the rest mass energy.
When \(E<0\), we take
so that
By the same reasoning, the solution for \(u_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[10]/span, line 1, column 1

so that in the limit \({\bf p}=0\), and \(E_0=-mc^2\),

which describes particles moving backward in times. Thus, the interpretation is that the negative energy solutions correspond to anti-particles, the the components, \(\phi_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[12]/span[1], line 1, column 1
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[12]/span[2], line 1, column 1
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[12]/span[3], line 1, column 1
In the non-relativistic limit, for \(E>0\), we have
so that
since \(mc^2 \gg {\bf p}^2/2m\), it follows that
Neglecting it, and recalling that for \(E>0\),
the eigenfunctions reduce to

The lower component has become reduntant, and the eigenfunctions just correspond to those of a free particle with an attached spin eigenfunction \({\left(\matrix{1 \cr 0}\right)}\) or \({\left(\matrix{0 \cr 1}\right)}\) for \(m_s=\hbar/2\) or \(-\hbar/2\), respectively. For \(E>0\), the lower component, \(\chi_
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[18]/span[1], line 1, column 1
Callstack:
at (Courses/New_York_University/G25.2666:_Quantum_Chemistry_and_Dynamics/2:_General_theory_of_spin/2.7:_Solution_of_the_Dirac_equation_for_a_free_particle), /content/body/div[19]/div/p[18]/span[2], line 1, column 1


