4.1. Interaction Light-Matter 2-7-08
4.1. Interaction Light-Matter 2-7-08
4.1. Interaction Light-Matter 2-7-08
(4.1)
The Hamiltonian for the matter H M is generally (although not necessarily) time independent,
whereas the electromagnetic field H L and its interaction with the matter H LM are time-dependent.
A quantum mechanical treatment of the light would describe the light in terms of photons for
different modes of electromagnetic radiation, which we will describe later.
We will start with a common semiclassical treatment of the problem. For this approach we
treat the matter quantum mechanically, and treat the field classically. For the field we assume that
the light only presents a time-dependent interaction potential that acts on the matter, but the matter
doesnt influence the light. (Quantum mechanical energy conservation says that we expect that the
change in the matter to raise the quantum state of the system and annihilate a photon from the
field. We wont deal with this right now). We are just interested in the effect that the light has on
the matter. In that case, we can really ignore H L , and we have a Hamiltonian that can be solved in
the interaction picture representation:
H H M + H LM ( t )
= H0 + V (t )
(4.2)
Here, well derive the Hamiltonian for the light-matter interaction, the Electric Dipole
Hamiltonian. It is obtained by starting with the force experienced by a charged particle in an
electromagnetic field, developing a classical Hamiltonian for this system, and then substituting
quantum operators for the matter:
4-2
p i
x x
(4.3)
In order to get the classical Hamiltonian, we need to work through two steps: (1) We need
to describe electromagnetic fields, specifically in terms of a vector potential, and (2) we need to
describe how the electromagnetic field interacts with charged particles.
Brief summary of electrodynamics
Lets summarize the description of electromagnetic fields that we will use. A derivation of the
plane wave solutions to the electric and magnetic fields and vector potential is described in the
appendix. Also, it is helpful to review this material in Jackson1 or Cohen-Tannoudji, et al.2
> Maxwells Equations describe electric and magnetic fields ( E , B ) .
> To construct a Hamiltonian, we must describe the time-dependent interaction potential (rather
than a field).
> To construct the potential representation of E and B , you need a vector potential A ( r , t ) and
a scalar potential ( r , t ) . For electrostatics we normally think of the field being related to the
electrostatic potential through E = , but for a field that varies in time and in space, the
electrodynamic potential must be expressed in terms of both A and .
> In general an electromagnetic wave written in terms of the electric and magnetic fields requires
6 variables (the x,y, and z components of E and B). This is an overdetermined problem;
Maxwells equations constrain these. The potential representation has four variables
( Ax , Ay , Az and ), but these are still not uniquely determined. We choose a constraint a
representation or guage that allows us to uniquely describe the wave. Choosing a gauge such
that = 0 (Coulomb gauge) leads to a plane-wave description of E and B :
2 A ( r , t ) +
2
1 A(r ,t )
=0
c2
t 2
A = 0
1
2
Jackson, J. D. Classical Electrodynamics (John Wiley and Sons, New York, 1975).
Cohen-Tannoudji, C., Diu, B. & Lale, F. Quantum Mechanics (Wiley-Interscience, Paris, 1977), Appendix III.
(4.4)
(4.5)
4-3
This wave equation allows the vector potential to be written as a set of plane waves:
A ( r , t ) = A0 e
i k r t
) + A* ei( k r t ) .
0
(4.6)
This describes the wave oscillating in time at an angular frequency and propagating in space
in the direction along the wavevector k , with a spatial period = 2 k . The wave has an
amplitude A0 which is directed along the polarization unit vector . Since A = 0 , we see
that k = 0 or k . From the vector potential we can obtain E and B
E=
A
t
i ( k r t )
i ( k r t )
= i A0 e
e
(4.7)
B = A
i ( k r t )
i ( k r t )
= i ( k ) A0 e
e
(4.8)
magnetic
field
polarization
are
mutually
orthogonal k b .
Also, by comparing eq. (4.6) and (4.7) we see that
the vector potential oscillates as cos(t), whereas the
field oscillates as sin(t). If we define
then,
Note, E0 B0 = k = c .
1
E0 = i A0
2
(4.9)
1
B0 = i k A0
2
(4.10)
E ( r , t ) = E0 sin ( k r t )
(4.11)
B ( r , t ) = B0 b sin ( k r t ) .
(4.12)
4-4
Classical Hamiltonian for radiation field interacting with charged particle
Now, lets find a classical Hamiltonian that describes charged particles in a field in terms of the
vector potential. Start with Lorentz force3 on a particle with charge q:
F = q(E +v B) .
(4.13)
Here v is the velocity of the particle. Writing this for one direction (x) in terms of the Cartesian
components of E , v and B , we have:
Fx = q ( Ex + v y Bz vz By ) .
(4.14)
In Lagrangian mechanics, this force can be expressed in terms of the total potential energy U as
U d U
+
x dt vx
Fx =
(4.15)
Using the relationships that describe E and B in terms of A and , inserting into eq. (4.14), and
working it into the form of eq. (4.15), we can show that:
U = q qv A
(4.16)
This is derived in CTDL,4 and you can confirm by replacing it into eq. (4.15).
Now we can write a Lagrangian in terms of the kinetic and potential energy of the particle
L = T U
L=
1
mv 2 + qv A q
2
(4.17)
(4.18)
L
= mv + qA
v
(4.20)
1 p qA .
v=m
(
)
(4.21)
p=
3
4
(4.19)
4-5
H=
1
m
p ( p qA )
H=
1
2m
( p qA )
q
m
( p qA ) A + q
2
1
p qA ( r , t ) + q ( r , t )
2m
(4.22)
(4.23)
This is the classical Hamiltonian for a particle in an electromagnetic field. In the Coulomb
gauge ( = 0 ) , the last term is dropped.
We can write a Hamiltonian for a collection of particles in the absence of a external field
p2
H 0 = i + V0 ( ri ) .
i 2mi
(4.24)
2
1
H =
pi qi A ( ri ) ) + V0 ( ri ) .
(
i 2mi
(4.25)
Expanding:
H = H0
i
qi
( pi A + A pi ) +
2mi
qi
2m
i
(4.26)
Generally the last term is considered small compared to the cross term. This term should be
considered for extremely high field strength, which is nonperturbative and significantly distorts the
potential binding molecules together. One can estimate that this would start to play a role at
intensity levels >1015 W/cm2, which may be observed for very high energy and tightly focused
pulsed femtosecond lasers. So, for weak fields we have an expression that maps directly onto
solutions we can formulate in the interaction picture:
H = H0 + V (t )
V (t ) =
i
qi
( pi A + A pi ) .
2mi
(4.27)
(4.28)
Now we are in a position to substitute the quantum mechanical momentum for the classical. Here
the vector potential remains classical, and only modulates the interaction strength.
4-6
p = i
V (t ) =
i
(4.29)
i
qi ( i A + A i )
2mi
(4.30)
i qi
A i
mi
(4.31)
q
= i A pi
i mi
V (t ) =
(4.32)
Under most circumstances, we can neglect the wavevector dependence of the interaction
potential.
If
the
wavelength
of
the
field
is
much
larger
than
the
molecular
(4.33)
For interactions, with UV, visible, and infrared (but not X-ray) radiation, k ri r0 << 1 , and setting
r0 = 0 means that eik r 1 . We retain the second term for quadrupole transitions: charge
distribution interacting with gradient of electric field and magnetic dipole.
Now, using A0 = iE0 2 , we write (4.32) as
V (t ) =
iqE0
p eit p e+ it
2m
(4.34)
4-7
qE0
( p ) sin t
m
q
=
( E (t ) p )
m
V (t ) =
(4.35)
E
q
V ( t ) = i ( pi ) 0 sin t
i mi
(4.36)
We are seeking to use this Hamiltonian to evaluate the transition rates induced by V(t) from our
first-order perturbation theory expression. For a perturbation V ( t ) = V0 sin t the rate of
transitions induced by field is
wk =
Vk
( Ek E ) + ( Ek E + )
(4.37)
Now we evaluate the matrix elements of the EDH in the eigenstates for H0:
Vk = k V0
We can evaluate the matrix element k p
qE0
k p
m
(4.38)
Hamiltonian:
[r , H0 ] =
i p
.
m
(4.39)
m
k r H0 H0 r
i
m
=
k r E Ek k r
i
= im k k r .
=
(4.40)
So we have
Vk = iqE0
k
k r
(4.41)
4-8
or for a collection of particles
Vk = iE0
k qi ri
k
k
= iE0 k kl
= iE0
(4.42)
is the dipole operator, and kl is the transition dipole matrix element. We can see that it is the
quantum analog of the classical dipole moment, which describes the distribution of charge density
in the molecule:
= dr r ( r ) .
(4.43)
These expressions allow us to write in simplified form the well known interaction potential
for a dipole in a field:
V ( t ) = E (t )
(4.44)
Then the rate of transitions between quantum states induced by the electric field is
wk =
=
E0
E0
k2
2
kl ( Ek E ) + ( Ek E + )
2
(4.45)
kl (k ) + (k + )
Equation (4.45) is an expression for the absorption spectrum since the rate of transitions
can be related to the power absorbed from the field. More generally we would express the
absorption spectrum in terms of a sum over all initial and final states, the eigenstates of H0:
w fi =
i, f
The
strength
of
interaction
element fi f i .
E0
fi ( fi ) + ( fi + )
between
light
and
f i
matter
is
(4.46)
given
by
the
matrix
distribution between f and i to get effective absorption. This matrix element is the basis of
selection rules based on the symmetry of the states. The vector part says that the light field must
project onto the dipole moment. This allows information to be obtained on the orientation of
molecules, and forms the basis of rotational transitions.
4-9
Relaxation Leads to Line-broadening
Lets combine the results from the last two lectures, and describe absorption to a state that is
coupled to a continuum. What happens to the probability of absorption if the excited state decays
exponentially?
k relaxes exponentially
... for instance by coupling to continuum
Pk exp [ wnk t ]
t0
d k V t
(4.47)
i
bk = eik t Vk ( t )
t
or equivalently
(4.48)
We can add irreversible relaxation to the description of bk , following our early approach:
w
i
bk = eik t Vk ( t ) nk bk
2
t
(4.49)
i
bk = eik t sin t Vk nk bk ( t )
2
t
(4.50)
=
E0 k i(k
e
2i
+ )
i (k )t
wnk
bk ( t )
2
y ( t ) = A e at +
b ei t
.
a + i
(4.51)
(4.52)
4-10
bk ( t ) = A e wnk t /2 +
E0 k
e i ( k + ) t
ei(k )t
2i wnk / 2 + i (k + ) wnk / 2 + i (k )
(4.53)
i t
E0 k
e( k )
2 k iwnk / 2
(4.54)
E02 k
=
4 2
( k )
+ wnk 2 / 4
(4.55)
The linewidth is related to the relaxation rate from k into the continuum n. Also the linewidth is
related to the system rather than the manner in which we introduced the perturbation.