In atomic physics, the fine structure describes the splitting of the spectral lines of atoms due to electron spin and relativistic corrections to the nonrelativistic Schrödinger equation.
The gross structure of line spectra is the line spectra predicted by the quantum mechanics of nonrelativistic electrons with no spin. For a hydrogenic atom, the gross structure energy levels only depend on the principal quantum number n. However, a more accurate model takes into account relativistic and spin effects, which break the degeneracy of the energy levels and split the spectral lines. The scale of the fine structure splitting relative to the gross structure energies is on the order of (Zα)^{2}, where Z is the atomic number and α is the finestructure constant, a dimensionless number equal to approximately
1
/
137
.
The fine structure energy corrections can be obtained by using perturbation theory. To do this one adds three corrective terms to the Hamiltonian: the leading order relativistic correction to the kinetic energy, the correction due to the spinorbit coupling, and the Darwinian term. The full Hamiltonian is given by
H
=
H
0
+
H
k
i
n
e
t
i
c
+
H
s
o
+
H
D
a
r
w
i
n
i
a
n
,
where
H
0
is the Hamiltonian from the Coulomb interaction. These corrections can also be obtained from the nonrelativistic limit of the Dirac equation, since Dirac's theory naturally incorporates relativity and spin interactions.
Classically, the kinetic energy term of the Hamiltonian is
T
=
p
2
2
m
where
p
is the momentum and
m
is the mass of the electron.
However, when considering a more accurate theory of nature via. special relativity, we must use a relativistic form of the kinetic energy,
T
=
p
2
c
2
+
m
2
c
4
−
m
c
2
where the first term is the total relativistic energy, and the second term is the rest energy of the electron. (
c
is the speed of light) Expanding this in a Taylor series ( specifically a binomial series ), we find
T
=
p
2
2
m
−
p
4
8
m
3
c
2
+
⋯
Then, the first order correction to the Hamiltonian is
H
k
i
n
e
t
i
c
=
−
p
4
8
m
3
c
2
Using this as a perturbation, we can calculate the first order energy corrections due to relativistic effects.
E
n
(
1
)
=
⟨
ψ
0

H
′

ψ
0
⟩
=
−
1
8
m
3
c
2
⟨
ψ
0

p
4

ψ
0
⟩
=
−
1
8
m
3
c
2
⟨
ψ
0

p
2
p
2

ψ
0
⟩
where
ψ
0
is the unperturbed wave function. Recalling the unperturbed Hamiltonian, we see
H
0

ψ
0
⟩
=
E
n

ψ
0
⟩
(
p
2
2
m
+
V
)

ψ
0
⟩
=
E
n

ψ
0
⟩
p
2

ψ
0
⟩
=
2
m
(
E
n
−
V
)

ψ
0
⟩
We can use this result to further calculate the relativistic correction:
E
n
(
1
)
=
−
1
8
m
3
c
2
⟨
ψ
0

p
2
p
2

ψ
0
⟩
E
n
(
1
)
=
−
1
8
m
3
c
2
⟨
ψ
0

(
2
m
)
2
(
E
n
−
V
)
2

ψ
0
⟩
E
n
(
1
)
=
−
1
2
m
c
2
(
E
n
2
−
2
E
n
⟨
V
⟩
+
⟨
V
2
⟩
)
For the hydrogen atom,
V
(
r
)
=
−
e
2
4
π
ε
0
r
,
⟨
1
r
⟩
=
1
a
0
n
2
, and
⟨
1
r
2
⟩
=
1
(
l
+
1
/
2
)
n
3
a
0
2
where
a
0
is the Bohr Radius,
n
is the principal quantum number and
l
is the azimuthal quantum number. Therefore, the first order relativistic correction for the hydrogen atom is
E
n
(
1
)
=
−
1
2
m
c
2
(
E
n
2
+
2
E
n
e
2
4
π
ε
0
1
a
0
n
2
+
1
16
π
2
ε
0
2
e
4
(
l
+
1
2
)
n
3
a
0
2
)
=
−
E
n
2
2
m
c
2
(
4
n
l
+
1
2
−
3
)
where we have used:
E
n
=
−
e
2
8
π
ε
0
a
0
n
2
On final calculation, the order of magnitude for the relativistic correction to the ground state is
−
9.056
×
10
−
4
eV
.
For a hydrogenlike atom with
Z
protons, orbital momentum
L
→
and electron spin
S
→
, the spinorbit term is given by:
H
s
o
=
1
2
(
Z
e
2
4
π
ε
0
)
(
g
s
2
m
e
2
c
2
)
L
→
⋅
S
→
r
3
m
e
is the electron mass,
ε
0
is the vacuum permittivity and
g
s
is the spin gfactor.
r
is the distance of the electron from the nucleus.
The spinorbit correction can be understood by shifting from the standard frame of reference (where the electron orbits the nucleus) into one where the electron is stationary and the nucleus instead orbits it. In this case the orbiting nucleus functions as an effective current loop, which in turn will generate a magnetic field. However, the electron itself has a magnetic moment due to its intrinsic angular momentum. The two magnetic vectors,
B
→
and
μ
→
s
couple together so that there is a certain energy cost depending on their relative orientation. This gives rise to the energy correction of the form
Δ
E
S
O
=
ξ
(
r
)
L
→
⋅
S
→
Notice that there is a factor of 2, called the Thomas precession, which comes from the relativistic calculation that changes back to the electron's frame from the nucleus frame.
Since
⟨
1
r
3
⟩
=
Z
3
n
3
a
0
3
1
l
(
l
+
1
2
)
(
l
+
1
)
⟨
L
→
⋅
S
→
⟩
=
ℏ
2
2
[
j
(
j
+
1
)
−
l
(
l
+
1
)
−
s
(
s
+
1
)
]
the expectation value for the Hamiltonian is:
⟨
H
S
O
⟩
=
E
n
2
m
e
c
2
n
j
(
j
+
1
)
−
l
(
l
+
1
)
−
3
4
l
(
l
+
1
2
)
(
l
+
1
)
Thus the order of magnitude for the spinorbital coupling is
Z
4
n
3
(
j
+
1
/
2
)
10
−
5
eV
.
Remark: On the (n,l,s)=(n,0,1/2) and (n,l,s)=(n,1,1/2) energy level, which the fine structure said their level are the same. If we take the gfactor to be 2.0031904622, then, the calculated energy level will be different by using 2 as gfactor. Only using 2 as the gfactor, we can match the energy level in the 1st order approximation of the relativistic correction. When using the higher order approximation for the relativistic term, the 2.0031904622 gfactor may agree with each other. However, if we use the gfactor as 2.0031904622, the result does not agree with the formula, which included every effect.
There is one last term in the nonrelativistic expansion of the Dirac equation. Because it was first derived by Charles Galton Darwin it is referred to as the Darwin term, and is given by:
H
D
a
r
w
i
n
i
a
n
=
ℏ
2
8
m
e
2
c
2
4
π
(
Z
e
2
4
π
ε
0
)
δ
3
(
r
→
)
⟨
H
D
a
r
w
i
n
i
a
n
⟩
=
ℏ
2
8
m
e
2
c
2
4
π
(
Z
e
2
4
π
ε
0
)

ψ
(
0
)

2
ψ
(
0
)
=
0
for
l
>
0
ψ
(
0
)
=
1
4
π
2
(
Z
n
a
0
)
3
2
for
l
=
0
H
D
a
r
w
i
n
i
a
n
=
2
n
m
e
c
2
E
n
2
The Darwin term affects only the sorbit. This is because the wave function of an electron with
l
>
0
vanishes at the origin, hence the delta function has no effect. For example, it gives the 2sorbit the same energy as the 2porbit by raising the 2sstate by 6977145109124427590♠9.057×10^{−5} eV.
The Darwin term changes the effective potential at the nucleus. It can be interpreted as a smearing out of the electrostatic interaction between the electron and nucleus due to zitterbewegung, or rapid quantum oscillations, of the electron. This can be demonstrated by a short calculation
Quantum fluctuations allow for the creation of virtual electronpositron pairs with a lifetime estimated by the uncertainty principle
Δ
t
≈
ℏ
/
Δ
E
≈
ℏ
/
m
c
2
. The distance the particles can move during this time is
ξ
≈
c
Δ
t
≈
ℏ
/
m
c
=
λ
c
, the Compton wavelength. The electrons of the atom interact with those pairs. This yields a fluctuating electron position
r
→
+
ξ
→
. Using a Taylor expansion, the effect on the potential
U
can be estimated:
U
(
r
→
+
ξ
→
)
≈
U
(
r
→
)
+
ξ
⋅
∇
U
(
r
→
)
+
1
2
∑
i
j
ξ
i
ξ
j
∂
i
∂
j
U
(
r
→
)
Averaging over the fluctuations
ξ
→
ξ
¯
=
0
,
ξ
i
ξ
j
¯
=
1
3
ξ
→
2
¯
δ
i
j
,
gives the average potential
U
(
r
→
+
ξ
→
)
¯
=
U
(
r
→
)
+
1
6
ξ
→
2
¯
∇
2
U
(
r
→
)
.
Approximating
ξ
→
2
¯
≈
λ
c
2
, this yields the perturbation of the potential due to fluctuations:
δ
U
≈
1
6
λ
c
2
∇
2
U
=
ℏ
2
6
m
2
c
2
∇
2
U
To compare with the expression above, plug in the Coulomb potential:
∇
2
U
=
−
∇
2
Z
e
2
4
π
ε
0
r
=
4
π
(
Z
e
2
4
π
ε
0
)
δ
(
r
→
)
⇒
δ
U
≈
ℏ
2
6
m
2
c
2
4
π
(
Z
e
2
4
π
ε
0
)
δ
(
r
→
)
This is only slightly different.
Another mechanism that affects only the sstate is the Lamb shift, a further, smaller correction that arises in quantum electrodynamics that should not be confused with the Darwin term. The Darwin term gives the sstate and pstate the same energy, but the Lamb shift makes the sstate higher in energy than the pstate.
The total effect, obtained by summing the three components up, is given by the following expression:
Δ
E
=
E
n
(
Z
α
)
2
n
(
1
j
+
1
2
−
3
4
n
)
,
where
j
is the total angular momentum (
j
=
1
/
2
if
l
=
0
and
j
=
l
±
1
/
2
otherwise). It is worth noting that this expression was first obtained by A. Sommerfeld based on the old Bohr theory; i.e., before the modern quantum mechanics was formulated.
The total effect can also be obtained by using the Dirac equation. In this case, the electron is treated as nonrelativistic. The exact energies are given by
E
j
n
=
−
m
e
c
2
[
1
−
(
1
+
[
α
n
−
j
−
1
2
+
(
j
+
1
2
)
2
−
α
2
]
2
)
−
1
/
2
]
.
This expression, which contains all higher order terms that were left out in the other calculations, expands to first order to give the energy corrections derived from perturbation theory. However, this equation does not contain the hyperfine structure corrections, which are due to interactions with the nuclear spin. Other corrections from quantum field theory such as the Lamb shift and the anomalous magnetic dipole moment of the electron are not included.