Problems 5

Part III ✎ Problems Phillips pp. 103–108 Β· ~34 min read

  • parity
  • probability current
  • resonance
  • alpha decay

Three of these are printed wrong, and every correction was found by running the arithmetic rather than by consulting anything.

Eight problems in four pairs. 1–2 rebuild Β§5.1’s bound-state count on two new potentials and introduce parity , which chapter 8 onward cannot do without. 3–4 run the machinery backwards: one extracts a phase shift from a given well, the other extracts a potential from a given wave function β€” and turns out to be the hydrogen atom four chapters early. 5–6 supply the two things Β§5.2a asserted: that R+T=1R + T = 1, and the exact transmission above a barrier. 7–8 are Coulomb barriers, in both directions β€” protons tunnelling in to fuse, and alpha particles tunnelling out to decay.

Group 1 β€” parity, and counting bound states again

1 The symmetric well, and what parity is derivation

A particle of mass mm moves in

V(x)={0ifΒ βˆ’βˆž<x<βˆ’aβˆ’V0ifΒ βˆ’a<x<a0ifΒ a<x<+∞.V(x) = \begin{cases} 0 & \text{if } -\infty < x < -a\\ -V_0 & \text{if } -a < x < a\\ 0 & \text{if } a < x < +\infty.\end{cases}

Because VV is symmetric about x=0x = 0, the eigenfunctions split into two families: symmetric ones with ψ(x)=+ψ(βˆ’x)\psi(x) = +\psi(-x), said to have positive parity, and antisymmetric ones with ψ(x)=βˆ’Οˆ(βˆ’x)\psi(x) = -\psi(-x), negative parity.

(a) Show that a positive-parity eigenfunction with E=βˆ’β„2Ξ±2/2mE = -\hbar^2\alpha^2/2m has the form

ψ(x)={Ae+Ξ±xifΒ βˆ’βˆž<x<βˆ’aCcos⁑k0xifΒ βˆ’a<x<+aAeβˆ’Ξ±xifΒ +a<x<+∞,\psi(x) = \begin{cases} Ae^{+\alpha x} & \text{if } -\infty < x < -a\\ C\cos k_0x & \text{if } -a < x < +a\\ Ae^{-\alpha x} & \text{if } +a < x < +\infty,\end{cases}

with k0=2m(E+V0)/ℏ2k_0 = \sqrt{2m(E+V_0)/\hbar^2}.

(b) Show that continuity of ψ\psi and dψ/dx\mathrm d\psi/\mathrm dx at the well edges implies α=k0tan⁑k0a\alpha = k_0\tan k_0a.

(c) By solving Ξ±=k0tan⁑k0a\alpha = k_0\tan k_0a together with Ξ±2+k02=w2\alpha^2 + k_0^2 = w^2 graphically, where w=2mV0/ℏ2w = \sqrt{2mV_0/\hbar^2}, count the bound states.

(d) Confirm that the negative-parity bound-state energies are identical to those of the well in Fig. 5.1.

πŸ’‘ Phillips' own hint

For (a), the eigenvalue equation is Eq. (5.5) inside the well and Eq. (5.7) outside it. For (b) and (c), modify the algebra that led to Eq. (5.11) and the graphical construction of Fig. 5.2.

βœ“ Worked solution

(a) Outside the well E<V=0E < V = 0, so d2ψ/dx2=Ξ±2ψ\mathrm d^2\psi/\mathrm dx^2 = \alpha^2\psi and the solutions are eΒ±Ξ±xe^{\pm\alpha x}. Finiteness picks e+Ξ±xe^{+\alpha x} on the left and eβˆ’Ξ±xe^{-\alpha x} on the right. Inside, E+V0>0E + V_0 > 0 gives d2ψ/dx2=βˆ’k02ψ\mathrm d^2\psi/\mathrm dx^2 = -k_0^2\psi, so ψ\psi is a combination of cos⁑k0x\cos k_0x and sin⁑k0x\sin k_0x β€” and positive parity kills the sine, because sin⁑\sin is odd. The same coefficient AA appears on both outside pieces for the same reason.

That is the whole content of parity here: it is a symmetry of the Hamiltonian, and it halves the work by forbidding one of the two solutions in each region before any boundary condition is applied.

(b) Match at x=ax = a:

Ccos⁑k0a=Aeβˆ’Ξ±a,βˆ’Ck0sin⁑k0a=βˆ’Ξ±Aeβˆ’Ξ±aC\cos k_0a = Ae^{-\alpha a},\qquad -Ck_0\sin k_0a = -\alpha Ae^{-\alpha a}

Dividing the second by the first removes AA and CC together:

k0tan⁑k0a=αk_0\tan k_0a = \alpha

Matching at x=βˆ’ax = -a gives the identical condition β€” which is the payoff of parity. One interface does the work of two.

(c) See the correction below.

(d) For negative parity the inside solution is Csin⁑k0xC\sin k_0x, and the same division gives k0cot⁑k0a=βˆ’Ξ±k_0\cot k_0a = -\alpha β€” which is Eq. (5.11) exactly. The reason is worth stating: an antisymmetric function must vanish at the origin, and ψ(0)=0\psi(0) = 0 is precisely the condition the infinite wall imposed in Fig. 5.1. A symmetric well of half-width aa therefore contains a perfect copy of Fig. 5.1’s spectrum, hidden in its odd states.

2 The three-dimensional spherical well derivation

A particle moves in three dimensions in

V(r)={βˆ’V0ifΒ r<R0ifΒ r>R,V(r) = \begin{cases} -V_0 & \text{if } r < R\\ 0 & \text{if } r > R,\end{cases}

with eigenfunctions governed by βˆ’β„22mβˆ‡2ψ+V(r)ψ=Eψ-\frac{\hbar^2}{2m}\nabla^2\psi + V(r)\psi = E\psi. Consider spherically symmetric eigenfunctions ψ(r)=u(r)/r\psi(r) = u(r)/r.

(a) Show that βˆ’β„22md2udr2+V(r)u=Eu-\frac{\hbar^2}{2m}\frac{\mathrm d^2u}{\mathrm dr^2} + V(r)u = Eu.

(b) Solving for uu inside and outside and imposing the boundary conditions, show that a bound state with E=βˆ’β„2Ξ±2/2mE = -\hbar^2\alpha^2/2m has

ψ(r)={Csin⁑(k0r)/rifΒ r<RAeβˆ’Ξ±r/rifΒ r>R\psi(r) = \begin{cases} C\sin(k_0r)/r & \text{if } r < R\\ Ae^{-\alpha r}/r & \text{if } r > R\end{cases}

(c) Show that there is one such bound state if ℏ2Ο€28mR2<V0<9ℏ2Ο€28mR2\dfrac{\hbar^2\pi^2}{8mR^2} < V_0 < \dfrac{9\hbar^2\pi^2}{8mR^2}.

βœ“ Worked solution

(a) For a function of rr alone the radial Laplacian is βˆ‡2ψ=1rd2(rψ)dr2\nabla^2\psi = \frac{1}{r}\frac{\mathrm d^2(r\psi)}{\mathrm dr^2}. With ψ=u/r\psi = u/r we have rψ=ur\psi = u, so βˆ‡2ψ=1rd2udr2\nabla^2\psi = \frac{1}{r}\frac{\mathrm d^2u}{\mathrm dr^2}. Substituting and multiplying through by rr gives the stated equation β€” a one-dimensional eigenvalue problem in rr.

This is the same substitution Β§5.2b used for the two-proton problem, and the reason it is worth doing once carefully: it converts every spherically symmetric 3-D problem in this book into a 1-D one.

(b) Inside, uβ€²β€²=βˆ’k02uu'' = -k_0^2u so u=Csin⁑k0r+Cβ€²cos⁑k0ru = C\sin k_0r + C'\cos k_0r. Here is the new ingredient: ψ=u/r\psi = u/r must stay finite at the origin, so u(0)=0u(0) = 0, which kills the cosine. Outside, uβ€²β€²=Ξ±2uu'' = \alpha^2u and finiteness at infinity leaves u=Aeβˆ’Ξ±ru = Ae^{-\alpha r}. Dividing by rr gives the stated ψ\psi.

(c) Matching uu and uβ€²u' at r=Rr = R and dividing:

k0cot⁑k0R=βˆ’Ξ±k_0\cot k_0R = -\alpha

which is Eq. (5.11) with aβ†’Ra \to R β€” identical to Fig. 5.1’s condition, and for exactly the reason problem 1(d) gave: u(0)=0u(0) = 0 is an infinite wall at the origin. So the counting carries over unchanged: the first bound state needs w>Ο€/2Rw > \pi/2R, the second w>3Ο€/2Rw > 3\pi/2R, where w=2mV0/ℏ2w = \sqrt{2mV_0/\hbar^2}. Squaring those thresholds,

ℏ2Ο€28mR2<V0<9ℏ2Ο€28mR2\frac{\hbar^2\pi^2}{8mR^2} < V_0 < \frac{9\hbar^2\pi^2}{8mR^2}

as required.

Group 2 β€” running the machinery backwards

3 A phase shift, from a well and an energy numerical

A particle with energy E=ℏ2Ο€2/2ma2E = \hbar^2\pi^2/2ma^2 is scattered by the well of Fig. 5.1 with depth V0=2ℏ2Ο€2/ma2V_0 = 2\hbar^2\pi^2/ma^2. Use Eq. (5.19) to show that the phase shift can be taken as Ξ΄(E)=3.53\delta(E) = 3.53 radians.

(Any integer multiple of Ο€\pi may be added to a phase shift satisfying Eq. (5.19).)

βœ“ Worked solution

With ℏ=m=a=1\hbar = m = a = 1 the numbers are clean. From E=Ο€2/2E = \pi^2/2 and V0=2Ο€2V_0 = 2\pi^2:

k=2E=Ο€,k0=2(E+V0)=5 πk = \sqrt{2E} = \pi,\qquad k_0 = \sqrt{2(E+V_0)} = \sqrt{5}\,\pi

Eq. (5.19) reads k0cot⁑k0a=kcot⁑(ka+Ξ΄)k_0\cot k_0a = k\cot(ka+\delta). Since cot⁑\cot has period Ο€\pi and ka=Ο€ka = \pi, the left side collapses:

cot⁑δ=k0kcot⁑k0a=5 cot⁑(5 π)=2.4431\cot\delta = \frac{k_0}{k}\cot k_0a = \sqrt5\,\cot(\sqrt5\,\pi) = 2.4431

so Ξ΄=arccot⁑(2.4431)=0.3888\delta = \operatorname{arccot}(2.4431) = 0.3888, and the family of valid values is …,Β βˆ’2.7528,Β 0.3888,Β 3.5304, …\dots,\ -2.7528,\ 0.3888,\ 3.5304,\ \dots spaced by Ο€\pi. The book quotes the third of these.

Why the ambiguity is real and not sloppiness. Eq. (5.19) constrains sin⁑(kx+Ξ΄)\sin(kx+\delta) only through a cotangent, and shifting Ξ΄\delta by Ο€\pi flips the sign of the whole outside wave function β€” which is not an observable change, since ψ\psi and βˆ’Οˆ-\psi describe the same state. Only differences of Ξ΄\delta between energies matter, which is why Β§5.1b’s time delay Ο„=2ℏ dΞ΄/dE\tau = 2\hbar\,\mathrm d\delta/\mathrm dE is well defined even though Ξ΄\delta itself is not.

4 Given the wave function, find the potential verification

(a) A particle has definite energy E=βˆ’β„2Ξ±2/2mE = -\hbar^2\alpha^2/2m and eigenfunction

ψ(x)={Nxeβˆ’Ξ±xifΒ 0≀x<∞0elsewhere,\psi(x) = \begin{cases} Nxe^{-\alpha x} & \text{if } 0 \le x < \infty\\ 0 & \text{elsewhere,}\end{cases}

with NN and Ξ±\alpha positive real constants. Verify that the potential is

V(x)={βˆ’Ξ±β„2/mxifΒ 0≀x<∞∞elsewhere.V(x) = \begin{cases} -\alpha\hbar^2/mx & \text{if } 0 \le x < \infty\\ \infty & \text{elsewhere.}\end{cases}

(b) Given that this ψ\psi is the ground state, sketch the first excited state.

πŸ’‘ Phillips' own hint

For (a), substitute ψ\psi and EE into βˆ’β„22md2ψdx2+V(x)ψ=Eψ-\frac{\hbar^2}{2m}\frac{\mathrm d^2\psi}{\mathrm dx^2} + V(x)\psi = E\psi and find the V(x)V(x) that makes it hold. For (b), the first excited state has a node between x=0x = 0 and x=∞x = \infty.

βœ“ Worked solution

(a) This problem runs the eigenvalue equation backwards. Rearranged, it is not a differential equation at all but a formula:

V(x)=E+ℏ22mΟˆβ€²β€²ΟˆV(x) = E + \frac{\hbar^2}{2m}\frac{\psi''}{\psi}

Differentiate twice:

Οˆβ€²=N(1βˆ’Ξ±x)eβˆ’Ξ±x,Οˆβ€²β€²=N(Ξ±2xβˆ’2Ξ±)eβˆ’Ξ±x\psi' = N(1-\alpha x)e^{-\alpha x},\qquad \psi'' = N(\alpha^2x - 2\alpha)e^{-\alpha x}

so Οˆβ€²β€²/ψ=Ξ±2βˆ’2Ξ±/x\psi''/\psi = \alpha^2 - 2\alpha/x. Substituting, with E=βˆ’β„2Ξ±2/2mE = -\hbar^2\alpha^2/2m:

V(x)=βˆ’β„2Ξ±22m+ℏ22m(Ξ±2βˆ’2Ξ±x)=βˆ’Ξ±β„2mxV(x) = -\frac{\hbar^2\alpha^2}{2m} + \frac{\hbar^2}{2m}\left(\alpha^2 - \frac{2\alpha}{x}\right) = -\frac{\alpha\hbar^2}{mx}

The Ξ±2\alpha^2 terms cancel exactly β€” which is what makes this ψ\psi an eigenfunction rather than merely a function. And V=∞V = \infty for x<0x < 0 is what forces ψ(0)=0\psi(0) = 0, which the factor of xx in Nxeβˆ’Ξ±xNxe^{-\alpha x} delivers.

(b) The first excited state has one node at some x>0x > 0: it rises, comes back through zero, dips to a negative minimum, and decays to zero from below. It is also more spread out, because it is less tightly bound.

Group 3 β€” the two things Β§5.2a asserted

5 R + T = 1, from the probability current derivation

For a stationary state Ξ¨(x,t)=ψ(x)eβˆ’iEt/ℏ\Psi(x,t) = \psi(x)e^{-iEt/\hbar} incident on the barrier of Fig. 5.4, with ψE=AIeikx+AReβˆ’ikx\psi_E = A_Ie^{ikx} + A_Re^{-ikx} on the left and ψE=ATeikx\psi_E = A_Te^{ikx} on the right:

(a) Show that on the left j=∣AI∣2ℏkmβˆ’βˆ£AR∣2ℏkmj = |A_I|^2\dfrac{\hbar k}{m} - |A_R|^2\dfrac{\hbar k}{m}, and on the right j=∣AT∣2ℏkmj = |A_T|^2\dfrac{\hbar k}{m}.

(b) By noting that ∣Ψ∣2|\Psi|^2 is constant in time for a stationary state, show that ∣AI∣2=∣AR∣2+∣AT∣2|A_I|^2 = |A_R|^2 + |A_T|^2, i.e. R+T=1R + T = 1.

βœ“ Worked solution

(a) The time factor eβˆ’iEt/ℏe^{-iEt/\hbar} cancels between Ξ¨\Psi and Ξ¨βˆ—\Psi^*, so jj may be computed from ψE\psi_E alone. On the right this is immediate: with ψ=ATeikx\psi = A_Te^{ikx}, βˆ‚Οˆ/βˆ‚x=ikATeikx\partial\psi/\partial x = ikA_Te^{ikx}, and

j=iℏ2m[ATβˆ—eβˆ’ikx(βˆ’ik)ATeikxβ‹…(βˆ’1)+… ]=∣AT∣2ℏkmj = \frac{i\hbar}{2m}\left[A_T^*e^{-ikx}(-ik)A_Te^{ikx}\cdot(-1) + \dots\right] = |A_T|^2\frac{\hbar k}{m}

On the left the algebra produces four terms. Two of them are the incident and reflected currents, +∣AI∣2ℏk/m+|A_I|^2\hbar k/m and βˆ’βˆ£AR∣2ℏk/m-|A_R|^2\hbar k/m β€” note the sign, which is the reflected wave carrying probability the other way. The two cross terms, proportional to AIβˆ—AReβˆ’2ikxA_I^*A_R e^{-2ikx} and its conjugate, cancel exactly, because they enter as a quantity minus its own complex conjugate multiplied by ii β€” leaving something real that turns out to be zero. That cancellation is the crux: without it jj would depend on xx, and no conservation statement would be possible.

(b) For a stationary state ∣Ψ(x,t)∣2=∣ψE(x)∣2|\Psi(x,t)|^2 = |\psi_E(x)|^2 has no time dependence at all, so for any x1,x2x_1, x_2,

j(x2,t)βˆ’j(x1,t)=ddt[∫x1x2∣Ψ∣2 dx]=0j(x_2,t) - j(x_1,t) = \frac{\mathrm d}{\mathrm dt}\left[\int_{x_1}^{x_2}|\Psi|^2\,\mathrm dx\right] = 0

The current is therefore the same everywhere. Equating its value on the left to its value on the right:

∣AI∣2ℏkmβˆ’βˆ£AR∣2ℏkm=∣AT∣2ℏkm|A_I|^2\frac{\hbar k}{m} - |A_R|^2\frac{\hbar k}{m} = |A_T|^2\frac{\hbar k}{m}

The common factor ℏk/m\hbar k/m cancels β€” and it can only cancel because kk is the same on both sides, which is true here because V=0V = 0 on both sides and would not be true for a step. Dividing by ∣AI∣2|A_I|^2 gives 1βˆ’R=T1 - R = T.

6 Above the barrier: exact T, and resonances derivation

For a particle with energy EE above the barrier of Fig. 5.4, with kBk_B and kk defined by E=ℏ2kB2/2m+VB=ℏ2k2/2mE = \hbar^2k_B^2/2m + V_B = \hbar^2k^2/2m, show that

T=∣AT∣2∣AI∣2=11+(Ssin⁑kBa)2,R=(Ssin⁑kBa)21+(Ssin⁑kBa)2T = \frac{|A_T|^2}{|A_I|^2} = \frac{1}{1 + (S\sin k_Ba)^2},\qquad R = \frac{(S\sin k_Ba)^2}{1 + (S\sin k_Ba)^2}

where S=k2βˆ’kB22kkBS = \dfrac{k^2 - k_B^2}{2kk_B}. Show that the barrier is completely transparent for certain energies.

πŸ’‘ Phillips' own hint

ψ\psi is Eq. (5.29) for x<0x < 0, Eq. (5.31) for 0<x<a0 < x < a and Eq. (5.34) for x>ax > a. Impose continuity of ψ\psi and dψ/dx\mathrm d\psi/\mathrm dx at x=0x = 0 and x=ax = a.

βœ“ Worked solution

Write q≑kBq \equiv k_B and apply the four matching conditions. With AI=1A_I = 1 the result is the same closed form Β§5.2a derived, since the algebra never cared whether qq was real or imaginary:

AT=4kq eβˆ’ika(k+q)2eβˆ’iqaβˆ’(kβˆ’q)2eiqaA_T = \frac{4kq\,e^{-ika}}{(k+q)^2e^{-iqa} - (k-q)^2e^{iqa}}

Taking ∣AT∣2|A_T|^2 and simplifying the denominator with ∣e±iqa∣=1|e^{\pm iqa}|=1 gives

T=4k2q24k2q2+(k2βˆ’q2)2sin⁑2qa=11+(k2βˆ’q22kq)2sin⁑2qaT = \frac{4k^2q^2}{4k^2q^2 + (k^2-q^2)^2\sin^2 qa} = \frac{1}{1 + \left(\frac{k^2-q^2}{2kq}\right)^2\sin^2 qa}

which is the stated result with S=(k2βˆ’kB2)/2kkBS = (k^2-k_B^2)/2kk_B. Then R=1βˆ’TR = 1 - T gives the second formula immediately β€” so problem 5’s conservation law does half of problem 6’s work.

Complete transparency. T=1T = 1 requires sin⁑kBa=0\sin k_Ba = 0, i.e.

kBa=nΟ€,n=1,2,3,…k_Ba = n\pi,\qquad n = 1, 2, 3,\dots

The barrier width is then a whole number of half wavelengths, so the waves reflected from its two faces are exactly out of phase and cancel. Nothing is absorbed and nothing is reflected: the barrier becomes invisible.

Problem 6's result, drawn

00.51.01.52.02.53.000.20.40.60.81.0energy E / V_Bprobabilitycurrent E
  • T exact
  • T from Eq. (5.40)
  • R exact

Three things are worth waiting for on this plot. In the hatched region (E < V_B) T is not zero β€” that is tunnelling, and on the log scale you can watch it fall through ten orders of magnitude as the barrier widens. The dashed curve is Eq. (5.40): it tracks the exact answer beautifully deep in the hatched region and then fails near E = V_B, where it can even exceed 1, because eβˆ’2Ξ²a β†’ 1 there and the approximation’s one assumption has failed. And above the barrier, T returns to 1 at a discrete set of resonances, k_B a = nΟ€ β€” a whole number of half wavelengths across the barrier, which is the same interference condition that makes an anti-reflection coating work.

E / V_B 1.500
k_B (inside) 2.0000
R 0.216060
T 0.783940

R + T = 1.000000000000 β€” Eq. (5.36), and it holds to the last digit the arithmetic carries, for every setting of every slider. The energy is above the barrier, where classical physics promises certain transmission β€” yet R = 0.2161 of the particles come back. Reflection from a step down in wavelength is a wave effect with no classical counterpart at all.

Above the barrier the exact T oscillates and touches 1 exactly at k_B a = nΟ€. The dips between are where S sin k_B a is largest. Below the barrier, the same curve is Β§5.2a's tunnelling result β€” one formula covers both, because q = k_B above and q = iΞ² below.

Group 4 β€” Coulomb barriers, in and out

7 Two protons, then two helium nuclei numerical

Find the classical distance of closest approach for two protons with an energy of approach of 2Β keV2\ \mathrm{keV}. Estimate the probability that they penetrate the Coulomb barrier. Compare with the corresponding probability for two 4He^4\mathrm{He} nuclei at the same energy of approach.

πŸ’‘ Phillips' own hint

Use Eqs. (5.44), (5.50) and (5.51).

βœ“ Worked solution

The turning point. Eq. (5.44) with e2/4πϡ0=1.44Β MeV fme^2/4\pi\epsilon_0 = 1.44\ \mathrm{MeV\,fm}:

rC=e24πϡ0E=1.44Β MeV fm0.002Β MeV=720Β fmr_C = \frac{e^2}{4\pi\epsilon_0 E} = \frac{1.44\ \mathrm{MeV\,fm}}{0.002\ \mathrm{MeV}} = 720\ \mathrm{fm}

360 times the 2Β fm2\ \mathrm{fm} range of the nuclear force.

Two protons. EG=493Β keVE_G = 493\ \mathrm{keV} from Β§5.2b, so

EG/E=493.1/2=15.70,T≃eβˆ’15.70=1.5Γ—10βˆ’7\sqrt{E_G/E} = \sqrt{493.1/2} = 15.70,\qquad T \simeq e^{-15.70} = 1.5\times10^{-7}

Two helium nuclei. Eq. (5.51) must be generalized, because it was written for unit charges. Restoring them, e2β†’Z1Z2e2e^2 \to Z_1Z_2e^2:

EG=(Z1Z2e24πϡ0ℏc)22Ο€2ΞΌc2E_G = \left(\frac{Z_1Z_2e^2}{4\pi\epsilon_0\hbar c}\right)^2 2\pi^2\mu c^2

For 4He^4\mathrm{He}, Z1Z2=4Z_1Z_2 = 4 and ΞΌ=mHe/2β‰ˆ2mp\mu = m_{\mathrm{He}}/2 \approx 2m_p, so EGE_G grows by 42=164^2 = 16 from the charge and by about 4 from the mass β€” a factor of 64, giving EG=31.3Β MeVE_G = 31.3\ \mathrm{MeV}. The exponent grows by 64=8\sqrt{64} = 8:

EG/E=125.2,T≃eβˆ’125.2=4.3Γ—10βˆ’55\sqrt{E_G/E} = 125.2,\qquad T \simeq e^{-125.2} = 4.3\times10^{-55}

48 orders of magnitude below the proton case, at the same energy.

8 Alpha decay, and the Geiger–Nuttall law numerical

Tunnelling through a Coulomb barrier also governs alpha decay. In the simplest model the alpha particle is preformed inside the nucleus and trapped by a potential like Fig. 5.5; the decay rate Ξ»\lambda is the frequency Ξ½\nu with which it strikes the barrier, times the penetration probability of Eq. (5.50).

Write an approximate expression for Ξ»\lambda in terms of Ξ½\nu, EGE_G and the energy released EE. Given that 235U^{235}\mathrm{U} has a half-life of 7.1Γ—1087.1\times10^8 years with E=4.68Β MeVE = 4.68\ \mathrm{MeV}, estimate the half-life of 239Pu^{239}\mathrm{Pu}, for which E=5.24Β MeVE = 5.24\ \mathrm{MeV}.

πŸ’‘ Phillips' own hint

The half-life is inversely proportional to the decay rate.

βœ“ Worked solution

The expression.

λ≃νexp⁑[βˆ’EGE],t1/2=ln⁑2λ∝1Ξ½exp⁑[+EGE]\lambda \simeq \nu\exp\left[-\sqrt{\frac{E_G}{E}}\right],\qquad t_{1/2} = \frac{\ln 2}{\lambda} \propto \frac{1}{\nu}\exp\left[+\sqrt{\frac{E_G}{E}}\right]

The estimate. Both decays are an alpha particle leaving a heavy nucleus of almost the same charge (Z=90Z = 90 and 9292), so EGE_G and Ξ½\nu are very nearly common to the two and cancel in a ratio:

t1/2(Pu)t1/2(U)=exp⁑[EGEPuβˆ’EGEU]\frac{t_{1/2}(\mathrm{Pu})}{t_{1/2}(\mathrm{U})} = \exp\left[\sqrt{\frac{E_G}{E_{\mathrm{Pu}}}} - \sqrt{\frac{E_G}{E_{\mathrm U}}}\right]

With Z1Z2=2Γ—90Z_1Z_2 = 2\times90 and ΞΌβ‰ˆ3.93 u\mu \approx 3.93\,u, Eq. (5.51) gives EGβ‰ˆ1.25Γ—105Β MeVE_G \approx 1.25\times10^5\ \mathrm{MeV}, so EG/E\sqrt{E_G/E} falls from 163.3163.3 to 154.3154.3 β€” a drop of 8.978.97 β€” and

t1/2(Pu)β‰ˆ7.1Γ—108Γ—eβˆ’8.97β‰ˆ9Γ—104Β yearst_{1/2}(\mathrm{Pu}) \approx 7.1\times10^8 \times e^{-8.97} \approx 9\times10^4\ \text{years}

The measured value is 2.4Γ—1042.4\times10^4 years. A factor of four, from a model with one adjustable input, across a range of four and a half orders of magnitude in half-life.

What the answer is really showing. A 12%12\% increase in the energy released shortens the half-life by a factor of nearly 80008000. This is the Geiger–Nuttall law, discovered empirically in 1911 and unexplained until Gamow derived exactly this exponent in 1928 β€” the first successful application of quantum mechanics to the nucleus. It is why measured alpha half-lives run from microseconds to longer than the age of the universe while the energies behind them vary by less than a factor of three.

Check yourself

0 / 6 answered

  1. 1.In problem 1 the well is symmetric, and matching at alone is enough β€” the condition from is automatically the same. Why?

  2. The book states there are two bound states for .

    2.How many bound states does the symmetric well actually hold at ?

  3. 3.Problem 4 asks you to verify that belongs to . What is this potential, really?

  4. 4.In problem 5, the factor cancels from every term and leaves . What does that cancellation quietly depend on?

  5. 5.The book prints problem 6's answer with where belongs. Why is this misprint so hard to catch?

  6. 6.Problem 8: the energy released rises 12% from U to Pu, and the estimated half-life falls by a factor of about 8000. Where does that leverage come from?