$$\require{cancel}$$

# 9.6 The Photoelectric Effect in Hydrogen

### Introduction

In the photoelectric effect, incoming light causes an atom to eject an electron. We consider the simplest possible scenario: that the atom is hydrogen in its ground state.  The interesting question is: for an ingoing light wave of definite frequency and amplitude, what is the probability of ionization of a hydrogen atom in a given time?  In other words, assuming we can use time-dependent perturbation theory, what is the ionization rate?

Formally, we know what to do.  We must find the interaction Hamiltonian H1, then use Fermi’s Golden Rule for the transition rate with a periodic perturbation:

Ri→f=2πℏ∣∣⟨f|H1|i⟩∣∣2δ(Ef−Ei−ℏω)

But it’s not that easy!  For one thing, the outgoing electron will be in some kind of plane wave state, so whatever convention we adopt for normalizing such states appears in the rate.  But also the δ function is tricky for excitation into the continuum: just how many of these plane wave states satisfy Ef=Ei+ℏω?  We shall discover that with a consistent formalism, these two difficulties cancel each other.

### The Interaction Hamiltonian

Taking the incoming wave to be an electromagnetic field having vector potential

A→(r→,t)=A→ 0cos(k→⋅r→−ωt)

The interaction Hamiltonian is given by replacing the electron kinetic energy term p→2/2m with(p→−qA→/c)2/2m.  The relevant new term is

−(1/2m)(q/c(p→⋅A→+A→⋅p→))=(e/mc)A→⋅p→

since q=−e and ∇→⋅A→=0 in our gauge.

Therefore

H1=(emc)cos(k→⋅r→−ωt)A→0⋅p→=(e2mc)(ei(k→⋅r→−ωt)+e−i(k→⋅r→−ωt))A→0⋅p→.

The two different terms in this expression, having time dependences e−iωt and eiωt will give δ functions δ(Ef−Ei−ℏω) and δ(Ef−Ei+ℏω) respectively in the transition rate.  The e−iωt term therefore corresponds to absorption of a photon, since we are looking at a process in which the electron gains energy, Ef>Ei.  The  eiωt term is for the process where an atom in an excited state emits a photon into the beam and drops in energy.

So the relevant interaction Hamiltonian is

H1(t)=H1e−iωt  where  H1=(e2mc)eik→⋅r→A→0⋅p→.

### Plane Waves: Density of States

We make the assumption that the final state is a plane wave state ∣∣∣k→f⟩∝eik→f⋅r→.

The most straightforward way of handling the plane wave states is to confine the whole system to an extremely large cubical box of side L, and impose periodic boundary conditions (so that plane traveling wave states are allowed).

The big box has volume V=L 3, so the appropriately normalized plane wave states are

∣∣∣k→⟩=1L3/2eik→⋅r→=1V−−√eik→⋅r→

As will become apparent, we need to count how thickly these states are distributed, both in momentum space (or k- space) and in energy.  We’ll begin by reviewing the one-dimensional problem—the three-dimensional case is a simple generalization.

Recall that for particles in one dimension confined to a line of length L with periodic boundary conditions, the allowed values of wave number k were given by eikL=1, so k=2nπ/L withn  an integer.  Thus considering only intervals Δk≫2π/L, the “density of states” in k is L/2π: an interval of length Δk contains (L/2π)Δk=ρ(k)Δk states, where here ρ(k)=L/2π. The density of states in energy, ρ(E), follows from differentiating E=ℏ2k2/2m.  Writing ΔE=(ℏ2k/m)Δk, gives the incremental change ΔE in E for a given incremental change Δk in k, so the two intervals ΔE and Δk must contain the same number of states, that is, ρ(E)ΔE=ρ(k)Δk. It then follows from ρ(k)=L/2π that the one-dimensional density of states in energy

ρ1D(E)=(L/2π)(m/ℏ2k)=(L/2πℏ)m/2E−−−−−√.

Note this one-dimensional density of states goes to infinity as E goes to zero.

In three dimensions, with a cube of side L and periodic boundary conditions, the density of states in k- space is (L/2π)3.  The allowed states can be visualized as the points of a cubic lattice, (kx,ky,kz)=2πL(nx,ny,nz), the n ’s being integers, so each allowed state has associated with it the volume of a small cube (2π/L)3.

To find the three-dimensional density of states in energy, using E=ℏ2k→2/2m, again ΔE=(ℏ2k/m)Δk but now to find the number of states in a small energy range we must multiply by 4πk2, since the states in the energy range lie between two close concentric spheres in k-space.  This gives

ρ(E)=(L/2π)34πk2(m/ℏ2k)=(L/2π)34πk(m/ℏ2)=(V/2π2)(m/ℏ3)2mE−−−−√.

Notice that in contrast to the one-dimensional case, the three-dimensional density of states goes to zero at zero energy.  (Exercise: What happens in two dimensions?)

(Of course, if we are detecting the ejected electron with apparatus restricted to a solid angle dΩ, the 4π is replaced by dΩ.)

The orthogonality condition between the plane wave states is

⟨k→∣∣∣k′→⟩=δk→,k′→

the ordinary Kronecker delta function—not Dirac’s—since the k ’s are an enumerated set,

(kx,ky,kz)=2πL(nx,ny,nz)

the n ’s being integers.

### Finding the Matrix Element

The ground state wave function for hydrogen is

|100⟩=1πa30−−−−√e−r/a0

.

The matrix element entering Fermi’s Golden Rule Ri→f=2πℏ∣∣⟨f|H1|i⟩∣∣2δ(Ef−Ei−ℏω) is therefore:

⟨k→f∣∣∣(e2mc)eik→⋅r→A→0⋅p→|100⟩=∫d3r(1/L)3/2e−ik→f⋅r→(e2mc)eik→⋅r→A→0⋅(−iℏ∇→)1πa30−−−−√e−r/a0

Actually the eik→⋅r→ term is not very important—the wavelength of incoming photons for the usual photoelectric effect is far greater than the size of the hydrogen atom in its ground state (which our integral is limited to) so eik→⋅r→≅1, and we can drop that term.

One point we’ve overlooked is that the electromagnetic wave has a magnetic field just as strong as the electric field, so what about the interaction of this magnetic field with the electron’s magnetic moment?  This turns out to be much weaker than the (e2mc)A→0⋅p→ term: the magnetic interaction

μ→B⋅B→=(e2mc)S→⋅B→

,

and the ratio of this magnetic contribution to the electric one is

(e2mc)S→⋅B→(e2mc)A→0⋅p→≃ℏσ→⋅∇→×A→A→⋅p→≃ℏkp

,

with p∼ℏ/a0, so this ratio is of order a0/λ,  λ being the wavelength of the incoming light, around 100 nm to ionize hydrogen.   So, we can safely ignore the magnetic interaction.

This interaction Hamiltonian H1=(e2mc)A→0⋅p→e−iωt is called the dipole approximation, because it can also be written in terms of the atom’s dipole moment er→.  To see how this comes about, taking |i⟩, |f⟩ to be eigenstates of H=p→2/2m+V(r→), and using [r→,p→]=iℏ,  [r→,H]=(iℏ/m)p→, we find the matrix elements of p→,r→ for this ionization transition are simply related

⟨f|p→|i⟩=(m/iℏ)⟨f|r→H−Hr→|i⟩=(m/iℏ)(Ei−Ef)⟨f|r→|i⟩=imω⟨f|r→|i⟩.

Therefore

⟨f|H1(t)|i⟩=(e2mc)A→0e−iωt⋅⟨f|p→|i⟩=(e2mc)imωA→0e−iωt⋅⟨f|r→|i⟩,

and E→=−(1/c)∂A→/∂t=(iω/2c)A→0e−iωt, from which ⟨f|H1(t)|i⟩=⟨f|−μ→⋅E→(t)|i⟩, with μ→=−er→, the atom’s electric dipole moment.

Nevertheless, for the particular interaction we are considering here, the p→=−iℏ∇→ representation proves more convenient.  (We’ll use the μ→=−er→ representation in later work.)

We must evaluate:

(1/L)3/2(e2mc)1πa30−−−−√∫d3re−ik→f⋅r→A→0⋅(−iℏ∇→)e−r/a0

.

Integration by parts gives the gradient operator acting on the plane wave state,

∫d3re−ik→f⋅r→A→0⋅(−iℏ∇→)e−r/a0=−(A→0⋅p→f)∫d3re−ik→f⋅r→e−r/a0.

The integral is now a Fourier transform of the hydrogen ground state wave function, and is straightforward: choose the z- axis in the direction of k→f, the φ -integration gives 2π, the θ -integration has sinθdθ=−d(cosθ), etc.  The result is (8π/a0)/(a−20+k2f)2.

Finally, we can put this into Fermi’s Golden Rule:

Ri→f=2πℏ∣∣∣∣∣(1/L)3/2(e2mc)1πa30−−−−√(A→0⋅p→f)⎛⎝⎜⎜8π/a0(a−20+k2f)2⎞⎠⎟⎟∣∣∣∣∣2δ(Ef−Ei−ℏω).

To detect the ejected electron, we will have a detector sensitive to some small solid angle, dΩ, not to some precise value of p→f.  There will also be some tiny uncertainty in ∣∣p→f∣∣, equivalent to an energy uncertainty, because for one thing the ejection takes place after a finite time.  This means the δ -function actually has finite width, and by taking our normalizing box big enough, there will be many states within this width—so, effectively, the δ -function is measuring the density of possible outgoing states (see discussion at end).  Recall the density of states in energy for outgoing solid angle dΩ is

ρ(E,dΩ)=(L/2π)3k2(m/ℏ2k)dΩ=(L/2π)3k(m/ℏ2)dΩ,

giving

Ri→f=2πℏ∣∣∣∣∣(1/L)3/2(e2mc)1πa30−−−−√(A→0⋅p→f)⎛⎝⎜⎜8π/a0(a−20+k2f)2⎞⎠⎟⎟∣∣∣∣∣2(L/2π)3kf(m/ℏ2)dΩ.

Notice first that the L3 terms cancel, reassuringly, our result cannot depend on the size of the box chosen for the plane wave states.  Writing pf=ℏkf, and of course p2f/2m=Ei+ℏω, we find

Ri→f=4mpfπa50ℏ4(emc)2(A→0⋅p→f)2(1a−20+(pf/ℏ)2)4dΩ.

Note that the rate is angle-dependent, since (A→0⋅p→f)2=A20p2fcos2θ: ejection is most likely parallel to the electric field.  The total ionization rate is given by integrating the rate over all angles, and on the unit sphere cos2θ¯¯¯¯¯¯¯¯=z2¯¯¯=1/3, so in the above, (A→0⋅p→f)2dΩ→4πA20p2f/3.

### The Photoelectric Cross Section

Imagine now sending this radiation into a gas of hydrogen atoms, many of them, but not enough to shade each other from the radiation significantly.  Energy will be absorbed from the beam as atoms ionize.  What is the rate at which the beam loses energy?  A convenient way of visualizing this rate of loss of energy is to replace each atom by a tiny perfectly absorbent disc oriented with its normal parallel to the beam, the size of these discs such that the beam loses energy at the same rate as it would by ionization.  The area of the disc equivalent to one atom is called the photoelectric cross section.

The energy density in the beam of radiation is

18π(∣∣∣E→∣∣∣2+∣∣∣B→∣∣∣2)=18π(2ω2c2A→20cos2(k→⋅r→−ωt))

Denoting the photoelectric cross section by σ,

energy absorbed per second=σ×c×energy density,

and averaging cos2, this gives the energy absorption rate per atom to be A20ω2σ/8πc.

However, if the rate of ionization of one atom is Ri→f, and that ionization takes energy ℏω from the beam, the rate of energy absorption is just ℏωRi→f, so the ionization cross section is given by

A20ω2σ/8πc=ℏωRi→f

This gives

σ=8πcA20ω2ℏωRi→f=8πcA20ω2ℏω4mpfπa50ℏ4(emc)24πA20p2f3(1a−20+(pf/ℏ)2)4=128ωe2a50ℏ3πp3f3mc(1a−20+(pf/ℏ)2)4.

### Appendix: the Golden Rule Delta Function and the Density of States

For our big box model, the states are infinite in number, but can be counted by going outwards from the origin in k- space, and adopting some convention for ordering those of equal energy.  We can label the states with n→=(nx,ny,nz), a vector with integer components placing the state in k- space, and denote its energy En→.  The contribution of this state to the density of states is a Dirac δ -function δ(E−En→), that is to say, this state contributes 1 to the density of states at the point En→ on the energy axis.  Therefore the density of states in energy ρ(E)=∑n→δ(E−En→), well approximated by the smooth function we derived above.

Now consider the integral over final plane wave states needed in the evaluation of the Golden Rule formula.  That δ -function has finite width (from the time-energy uncertainty principle), so by taking our big box big enough we can have many plane wave states within the width of the Golden Rule δ -function: to picture this, let’s represent it by a function equal to Δ over an interval 1/Δ, zero otherwise.  Then integrating this Golden Rule δ -function with ρ(E) will give a contribution Δ from each state inside the interval of width 1/Δ.  If the states were uniformly distributed in energy, this would give the total number of states in an interval of unit energy—and that is the definition of the density of states.