Chapter 3
Chapter 3
JJ II
J I
By L. Salasnich
Page 1 of 43
Go Back
1. Introduction
J I
5.Spontaneous Emission
Page 2 of 43 6. Absorption
Go Back
7. Stimulated Emission
Full Screen
8. Selection Rules
Close
9. Einstein Coefficients
Quit
1. Introduction
Full Screen I Finally, we discuss the lifetime of an atomic state and the line width of
Close an electromagnetic transition.
Quit
2. The Total and Interaction Hamiltonians
~
p~ → p~ − q A, and ~ →E
E ~ − qφ, (1)
~ is the vector potential and φ is the scalar potential. Recall
where A
∂
that the energy operator is E = i~ ∂t and the momentum operator is
p~ = −i~∇.
Home Page • The original problem in the absence of the fields is a Schrodinger
Title Page
problem HΨ = EΨ,
" #
Contents p~2 ∂Ψ(~r, t)
+ V (~r) Ψ(~r, t) = i~ . (2)
2m ∂t
JJ II
where V (~r) given by
J I
N N
X 1X qi q j
Page 5 of 43 V (~r) = qi ϕ i = . (3)
i=1
2 i,j=1 4π0 | ~ri − ~rj |
Go Back
Close
• When the EM fields are now turned on, this changes to
" #
Quit 1 ~ 2 + V (~r) + qφ Ψ(~r, t) = i~ ∂Ψ(~r, t) ,
(~p − q A) (4)
2m ∂t
1 ~ 2 + V (~r) + qφ + Hrad ,
Title Page H = (~p − q A) (6)
2m
Contents p~2 q ~ 2 2
~ + q A + qφ + Hrad .
= + V (~r) − (A · p~ + p~ · A)
2m 2m 2m
JJ II
- The first two terms are the free particle Hamiltonian,
J I
- the next four terms are due to the particle’s interaction with the
Page 6 of 43
EM fields, and
Go Back
- the last term is the free-photon Hamiltonian, whose eigenstates
Full Screen
are states with certain numbers of photons of given wave vec-
Close tors and polarizations.
~ · p~)ψ = A
(A ~ · (−i~∇ψ), and ~ = −i~∇ · (Aψ).
(~p · A)ψ ~
But,
~ = (∇ · A)ψ
∇ · Aψ ~ +A~ · ∇ψ = A
~ · ∇ψ,
JJ II p~2 q ~ q 2 A2
H= + V (~r) − (A · p~) + + Hrad . (7)
J I
2m m 2m
Page 7 of 43 • Note that there are only two interaction terms - one that is linear in
~ and one that is quadratic in A.
A ~
Go Back
Full Screen • The linear terms will involve individual photon creation and annihi-
lation terms, corresponding to single-photon processes that can be
Close
analyzed in first order perturbation theory (PT).
Quit
Contents where Ĥint describes the interaction between the atomic electrons
JJ II and the radiation field.
Contents
At this point, it is worth noting the following facts:
JJ II – In a basic process where the final state has one less photon than
J I the initial state, a photon was absorbed by the matter. On the
other hand, if the final state has one more photon than the ini-
Page 9 of 43
tial state, a photon was emitted by the matter.
Go Back
– In absorption, the absorbed photon could come from any di-
Full Screen
rection dΩ. In emission, the emitted photon can go out in any
Close direction dΩ.
Quit
In either case, some density of states is involved in the calcula-
tion of the rate for the process, consistent with conservation of
energy and other variables.
3. Classical Electrodynamics
Page 10 of 43
where p~2i /(2mi ) = p2i /(2mi ) is the kinetic energy of i-th particle.
Go Back
• The presence of the electromagnetic field is modelled by the Hamil-
Full Screen
tonian
Close H = Hshif t + Hrad , (13)
Quit where
N ~ i )2
X (~pi − qi A
Hshif t = + qi φi , (14)
i=1
2mi
~i = A
with A ~ i (~ri , t), φi = φi (~ri , t), where φ(~r, t) and A(~
~ r, t) are the elec-
tromagnetic potentials and
Z " #
0 ~2 1 ~ 2 (~r, t) .
Ĥrad = d3~r E (~r, t) + B (15)
2 T 2µ0
Home Page is the radiation Hamiltonian.
Title Page
• We observe that by using the Hamilton equations
Contents ∂Hshif t
~r˙i = . (16)
∂~pi
JJ II
J I ∂Hshif t
p~˙i = − . (17)
∂~ri
Page 11 of 43
on the shift Hamiltonian (14) it is straightforward to derive Newton
Go Back equations with the Lorentz force acting on the i-th particle as
Full Screen
h ∂A~i i
¨
mi~ri = qi ∇i φi − ~ ~ i + ~vi × B
+ ~vi × (∇i × Ai ) = qi (E ~ i) (18)
∂t
Close
where ~vi = ~r˙i and, as always, the electric field E
~ is obtained as
Quit
~
~ = −∇φ − ∂ A .
E (19)
∂t
~ can be written as
while the magnetic field B
~ = ∇ × A.
B ~ (20)
~ T which, in the
• In Eq. (15) it appears the transverse electric field E
Coulomb gauge ∇ · A ~ = 0, is related to the total electric field E
~ by
~ L and transverse E
the following decomposition into longitudinal E ~T
fields
Home Page
~ =E
E ~L + E
~T (21)
Title Page
such that
~
Contents ~ L = ∇φ,
E ~ T = − ∂A .
E (22)
∂t
JJ II
~ L gives rise to the instan-
• Remarkably, the longitudinal electric field E
J I
taneous electrostatic interaction between the charges. Indeed, Eq.
Page 12 of 43 (14) can be rewritten as
N ~ i )2
Go Back X (~pi − qi A
Hshif t = + HC , (23)
Full Screen i=1
2mi
where
Close N N
X 1X qi qj
HC = qi φi = . (24)
Quit
i=1
2 i,j=1 4π0 | ~ri − ~rj |
is the Hamiltonian of the Coulomb interaction.
Go Back
• Notice, however, that this Hamiltonian does not take into account
Full Screen
the possible spin of particles.
Close
Quit
4. Quantum Electrodynamics in the Dipole Ap-
proximation
Home Page
• The quantization of electrodynamics is obtained promoting the clas-
sical Hamiltonian of the system to a quantum Hamiltonian.
Title Page
Contents
• For simplicity we consider the hydrogen atom with Hamiltonian
p~2 e2
JJ II Ĥmatt = − , (29)
2m 4π0 | ~r |
J I
where p~ = −i~∇ is the linear momentum operator of the electron in
Page 14 of 43 the state | p~i, in the presence of the radiation field with Hamiltonian
Go Back XX
Ĥrad = ~ωk â†ks âks , (30)
Full Screen s
k
Close where âks and â†ks are the annihilation and creation operators of the
photon in the state | ksi.
Quit
Contents
H = H 0 + HD , (34)
JJ II
where
J I
H = H0 + Hrad , (35)
Page 16 of 43
is the unperturbed Hamiltonian, whose eigenstates are of the form
Go Back
| âi | ...nks ...i =| âi⊗ | ...nks ...i, (36)
Full Screen
Close where | âi is the eigenstate of Ĥmatt with eigenvalue Ea and | ...nks ...i
P
is the eigenstate of Ĥrad with eigenvalue ks ~ωk nks , i.e.,
Quit
Ĥ0 | âi | ...n̂ks ...i = (Ĥmatt + Ĥrad ) | âi | ...nks ...i (37)
X
= Ea + ~ωk nks | âi | ...nks ...i.
ks
The time-dependent perturbing Hamiltonian is instead given by
r
e X ~ h i
ĤD (t) = âks e−iωk t + â†ks eiωk t ûks · p~ˆ. (38)
m ks 20 ωk V
Home Page
Fermi golden rule: Given the initial | ii and final | f i eigenstates
Title Page
of the unperturbed Hamiltonian Ĥ0 under the presence to the time-
Contents dependent perturbing Hamiltonian ĤD , the probability per unit time
JJ II
of the transition from | ii to | f i is given by
2π
J I Wif = | h| f | ĤD (0) | ii |2 δ(Ei − Ef ), (39)
~
Page 17 of 43
with the constraint of energy conservation.
Go Back
• This is the so-called Fermi golden rule, derived in 1926 by Paul Dirac
Full Screen
on the basis of the first order time-dependent perturbation theory,
Close
and named “golden rule” few years later by Enrico Fermi.
Quit
5. Spontaneous Emission
• Let us now apply the Fermi golden rule to the very interesting case
of the hydrogen atom in the state | bi and the radiation field in the
Home Page
vacuum state | 0i.
Title Page
JJ II | ii =| bi | 0i, (40)
J I
• Notice that, because we are considering the hydrogen atom, one has
Page 18 of 43
Go Back
Ĥmatt | bi = Eb | bi, (41)
Full Screen
where
mc2 α2 13.6 eV
Eb = − 2
'− , (42)
Close 2nb n2b
Quit
is the well-known quantization formula of the nonrelativistic hydro-
gen atom with quantum number nb = 1, 2, 3, ... and
e2 1
α= ' . (43)
(4π0 )~c 137
is the fine-structure constant.
Home Page
Title Page
Contents
JJ II
Go Back
• In addition we suppose that the final state is
Full Screen
| f i =| ai | nks i, (44)
Close
Quit
i.e., the final atomic state is | ai and the final photon state is | ksi =|
1ks i = â†ks | 0i.
• Obviously, we have
Ĥmatt | ai = Ea | ai, (45)
where
mc2 α2 13.6 eV
Ea = − 2
'− , (46)
2na n2a
with na = 1, 2, 3, ....
Home Page
• In this process, where Eb > Ea , there is the spontaneous production
Title Page
of a photon from the vacuum: a phenomenon strictly related to the
Contents quantization of the electromagnetic field. From Eqs. (38) and (39)
JJ II one finds
2π e 2 ~
J I spo
Wba = | ûks · ha | p~ˆ | bi |2 δ(Eb − Ea − ~ωk ), (47)
~ m 20 ωk V
Page 20 of 43
because
Go Back
âk0 s0 | ii = âk0 s0 | bi | 0i =| biâk0 s0 | 0i = 0, (48)
Full Screen
while
Close
and consequently
m m
Close ha | p~ˆ | bi = ha | [~r, Ĥmatt ] | bi = ha | [~rĤmatt − Ĥmatt~r] | bi (53)
i~ i~
m
Quit = (Eb − Ea )hb | ~r | ai = −imωba ha | ~r | bi,
i~
where ωba = (Eb − Ea )/~, and consequently
2
spo πωba
Wba = | ûks · ha | e~r | bi |2 δ(~ωba − ~ωk ). (54)
V 0 ωk
• The delta function is eliminated by integrating over the final photon
states. The photon states allowed by the periodic boundary condi-
tion form a continuous energy spectrum as the normalization be-
comes infinite.
Home Page
Title Page
• For a photon emitted at a solid angle dΩ, the number of allowed
states in the energy interval [~ω, ~(ω+dω)] can be written as ρ~ω,dΩ d(~ω),
Contents
where
JJ II V | ~k |2 d | ~k | dΩ V ω2
ρ~ω,dΩ = = dΩ.
(2π)3 d(~ω) 8π 3 ~c3
J I
• Consequently, the probability of spontaneous emission in the solid
Page 22 of 43
angle dΩ becomes
V k 2 (ωba )
Go Back X X sp Z X sp
spo
Wba = Wba,ks = 3)
dΩ Wba,ks , (55)
Full Screen k s
(2π s=1,2
Quit • Because ûk1 , ûk2 , and n̂ = ~k/k form a orthonormal system of vectors,
setting ~rab = ha | ~r | bi one finds
Home Page
namely
X
Title Page | ûks · ha | ~r | bi |2 =| ha | ~r | bi |2 sin2 (θ). (58)
s=1,2
Contents
Contents
JJ II
J I
Page 24 of 43
Go Back
Full Screen
Close
Quit
6. Absorption
Title Page
• Thus we suppose that the initial state is
Contents
| ii =| bi | nks i, (61)
Page 25 of 43 where Ea < Eb . From Eqs. (38) and (39) one finds
2π e 2 ~
Go Back
abs
Wab = nks | ûks ·hb | p~ˆ | ai |2 δ(Ea +~ωk −Eb ), (63)
~ m 20 ωk V
Full Screen because
√
Close hf | p~ˆâk0 s0 | ii = nks hb | p~ˆ | aiδk0 ,k δs0 ,s , hf | p~ˆâ†k0 s0 | ii = 0. (64)
Quit
• Note that with respect to Eq. (47) in Eq. (63) there is also the multi-
plicative term nks . We can follow the procedure of the previous sec-
tion to get
2
abs πωba
Wab = nks | ûks · hb | e~r | ai |2 δ(~ωba − ~ωk ). (65)
V 0 ωk
• Again the delta function can be eliminated by integrating over the
final photon states but here one must choose the functional depen-
dence of nks . We simply set
Title Page and after integration over k and s, from Eq. (65) we get
3
ωba
Contents
abs
Wab = 3
| hb | d~ | ai |2 n(ωba ) = Wba
spo
n(ωba ). (67)
3π0 ~c
JJ II
Full Screen • For a photon emitted at a solid angle dΩ, the number of allowed
states in the energy interval [~ω, ~(ω+dω)] can be written as ρ~ω,dΩ d(~ω),
Close
where
Quit V | ~k |2 d | ~k | dΩ
ρ~ω,dΩ = .
(2π)3 d(~ω)
Equivalently, this can be written in terms of the frequency interval
[ω, ω + dω] to be
V ω 2 dΩ
ρω,dΩ = .
8π 3 c3
Home Page
Title Page
Contents
JJ II
J I
Figure 2: Optical absorption.
Page 27 of 43
so that the energy density (density of allowed states per unit volume,
Go Back
V ) of the states, n(ω), with energy ~ω is obtained by integrating over
Full Screen
the solid angle, dΩ. That is, (taking the spin into account)
ω 2 [~ωn(ω)] ~ω 3
Close
Z
ρ(ω) = 2 × dΩ = n(ω). (68)
Quit
8π 3 c3 π 2 c3
where
1
n(ω) = ,
e~ω/(kB T )
−1
is the Fermi-Dirac distribution function with where kB is the Boltz-
mann constant and T the absolute temperature.
• Note that Eq. (68) represents the thermal distribution of photons,
with ρ(ω) the energy density per unit of angular frequency - the thermal-
equilibrium Planck formula.
J I
Page 28 of 43
Go Back
Full Screen
Close
Quit
7. Stimulated Emission
Title Page
• Thus, we suppose that the initial state is
while
√
â†k0 s0 | ii = â†k0 s0 | bi | nks i =| biâ†k0 s0 | nks i = nks + 1 | bi | nks + 1i,
(74)
Home Page
Title Page
Contents
JJ II
Figure 3: Stimulated emission.
J I
Page 30 of 43
and consequently
Go Back √
hf | p~ˆâk0 s0 | ii = 0, hf | p~ˆâ†k0 s0 | ii = nks + 1ha | p~ˆ | biδk0 ,k δs0 ,s .
Full Screen
(75)
Close
Quit
• Note that with respect to Eq. (64) in Eq. (72) there is the factor nks + 1
instead of nks . It is straightforward to follow the procedure of the
previous section to obtain
3
ωba
sti
Wba = 3
| ha | d~ | bi |2 [n(ωba ) + 1] = Wab
abs spo
+ Wba , (76)
3π0 ~c
sti
which shows that the probability Wba of stimulated emission re-
spo
duces to the spontaneous one Wba when n(ωba ) = 0. It is then useful
to introduce
sti sti spo
W̃ba = Wba − Wba , (77)
Home Page
Home Page
Title Page
Contents
JJ II
J I
Page 32 of 43
Go Back
Full Screen
Close
Quit
8. Selection Rules
Contents
• It is possible to show that the matrix elements are zero for certain
JJ II
pairs of states. If so, the transition is not allowed (at least in the
J I electric dipole approximation), and the results can be summarized
Page 33 of 43
in terms of simple selection rules governing the allowed changes in
quantum numbers in transitions.
Go Back
Full Screen
• Since the electric dipole d~ = −e~r changes sign under parity (~r → −~r),
the matrix element hb || d~ | ai is zero if the states | ai and | bi have
Close
the same parity. Therefore
Quit
– the parity of the wavefunction must change in an electric dipole
transition.
• This means that if ψa (~r) and ψb (~r) are the eigenfunctions of the states
| ai and | bi and, for instance, both eigenfunctions are even, i.e.,
Page 34 of 43
• Let us recall that in our calculations a generic eigenstate | ai =| nlmi
of the matter Hamiltonian Ĥmat is such that
Go Back
Full Screen ψa (~r) = ψnlm (~r) = h~r | ai = h~r | nlmi = Rnl (r)Ylm (θ, φ), (81)
Close
where Ylm (θ, φ) is the spherical harmonic function.
Quit
• The spherical harmonics satisfy the orthonormalization condition
Z π Z 2π
sin θdθ Ylast
0 m0 (θ, φ)Ylm (θ, φ)dφ = δl,l0 δm,m0 . (82)
0 0
∗
where Yl,m (θ, φ) = (−1)m Yl,−m (θ, φ).
• Moreover under parity transformation one finds
ψnlm (−~r) = Rnl (r)Ylm (π−θ, φ+π) = (−1)l Rnl (r)Ylm (θ, φ) = (−1)l ψnlm (~r),
(83)
Home Page
thus the radial part Rnl (r) of the wave function is unchanged and the
Title Page parity of the state is fully determined from the angular part.
Contents
• The generic dipole matrix element is given by
JJ II Z ∞ Z π Z 2π
0 0 0 ~ 3
hn l m | d | nlmi = −e r Rn0 l0 (r)Rnl (r)dr sin θdθ Ylast
0 m0 (θ, φ)
J I 0 0 0
× (cos θ sin φ, sin θ sin φ, cos φ)Ylm (θ, φ)dφ. (84)
Page 35 of 43
Full Screen
• Setting ∆l = l − l0 and ∆m = m − m0 , by using the properties of
Close the spherical harmonics it is possible to prove that the dipole matrix
element hn0 l0 m0 | d~ | nlmi is different from zero if only if ∆l = ±1 and
Quit
∆m = 0, ± 1.
Home Page • This means that in the electric dipole transitions the photon carries
off (or brings in) one unit of angular momentum.
Title Page
Contents • It is important to stress that the above selection rules are obtained
within the dipole approximation. This means that they can be vio-
JJ II
lated by rare electromagnetic transitions involving higher multipo-
J I
larities.
Page 36 of 43
• Concluding this section we observe that in the general case of N par-
Go Back
ticles with positions ~ri and electric charges qi , i = 1, 2, ..., N , the elec-
Full Screen tric dipole momentum is defined
Close N
X
d~ = qi~ri . (86)
Quit i=1
~ˆ D = −d~ · E,
H ~ˆ (87)
which is the familiar expression of the coupling between the electric
~ˆ
dipole d~ and the electric field E.
Home Page
Title Page
Contents
JJ II
J I
Page 37 of 43
Go Back
Full Screen
Close
Quit
9. Einstein Coefficients
J I and consequently
dNa dNb
=− . (89)
Page 38 of 43 dt dt
Go Back • Einstein suggested that, if the atoms are exposed to an electromag-
Full Screen
netic radiation of energy density per unit of frequency ρ(ω), the rate
of change of atoms in the state | ai must be
Close
dNa
Quit = Aba Nb + Bba ρ(ωba )Nb − Bab ρ(ωba )Na , (90)
dt
where the parameters Aba , Bba , and Bab are known as Einstein coeffi-
cients.
- The first two terms in this equation describe the increase of the
number of atoms in | ai due to spontaneous and stimulated
transitions from | bi,
Home Page - while the third term takes into account the decrease of the num-
ber of atoms in | ai due to absorption with consequent transi-
Title Page
tion to | bi.
Contents
• Einstein was able to determine the coefficients Aba , Bba , and Bab by
JJ II
supposing that the two rates in Eqs. (89) and (90) must be equal to
J I zero at thermal equilibrium, i.e.,
Page 39 of 43
dNa dNb
=− = 0. (91)
Go Back dt dt
Nb e−βEb
= −βEa = e−β(Eb −Ea ) = e−β~ωba . (93)
Na e
Einstein got
Aba
ρ(ωba ) = . (94)
Bab eβ~ωba − Bba
• At thermal equilibrium we know that
Home Page
3
~ωba 1
Title Page ρ(ωba ) = 2 3 β~ω
. (95)
π c e ba −1
Contents
• It follows that 3
JJ II ~ωba
Aba = Bba , Bab = Bba . (96)
π 2 c3
J I
• Notice that in this way Einstein obtained the coefficient Aba of spon-
Page 40 of 43
taneous decay by simply calculating the coefficient of stimulated de-
Go Back
cay Bba .
Full Screen
• Actually, by using the results we have obtained in the previous sec-
Close tion, it is clear that
Quit 3
ωba
spo
Aba = Wba = | ha | d~ | bi |2 , (97)
3π0 ~c3
sti
W̃ba π 2 c3
Bba = = 3
Aba , (98)
ρ(ωba ) ~ωba
abs
Wab π 2 c3
Bab = = 3
Aba . (99)
ρ(ωba ) ~ωba
Home Page
Title Page
Contents
JJ II
J I
Page 41 of 43
Go Back
Full Screen
Close
Quit
Main References:
Contents
JJ II
J I
Page 42 of 43
Go Back
Full Screen
Close
Quit
Home Page
Title Page
Contents
JJ II
J I Thank you !
Page 43 of 43
Go Back
Full Screen
Close
Quit