Problems 6

Part III ✎ Problems Phillips pp. 128–134 Β· ~24 min read

  • zero-point energy
  • raising and lowering operators
  • degeneracy

Nothing in this set solves a differential equation. Every answer comes from an inequality, an integral or a symmetry, which is the habit the chapter was really teaching.

Eleven problems in four groups. 1–3 attack the ground state from three independent directions β€” the uncertainty principle, the Gaussian, and the classically forbidden region β€” and between them explain why the zero-point energy is 12ℏω\frac12\hbar\omega and not something else. 4 and 8 change the potential: a spring that stretches but cannot compress, and the oscillator in two dimensions, which quietly builds angular-momentum eigenstates four chapters early. 5–7 set states moving and then measure a real molecule. 9–11 are marked β€œfor readers who studied Section 6.6” and are that section’s algebra β€” they supply three results the chapter has already used on credit.

Group 1 β€” the zero-point energy, three ways

1 The floor, from uncertainty alone derivation

Use Ξ”x Δpβ‰₯12ℏ\Delta x\,\Delta p \ge \frac12\hbar to derive a lower bound on the energy of a particle of mass mm in a harmonic oscillator of angular frequency Ο‰\omega.

(a) Given ⟨E⟩=⟨p2⟩2m+12mΟ‰2⟨x2⟩\langle E\rangle = \dfrac{\langle p^2\rangle}{2m} + \dfrac12m\omega^2\langle x^2\rangle, show that if ⟨x⟩=⟨p⟩=0\langle x\rangle = \langle p\rangle = 0 then

⟨E⟩β‰₯ℏ28m(Ξ”x)2+12mΟ‰2(Ξ”x)2\langle E\rangle \ge \frac{\hbar^2}{8m(\Delta x)^2} + \frac12m\omega^2(\Delta x)^2

(b) Show that the minimum value of F(Ξ”x)=A2(Ξ”x)2+B2(Ξ”x)2F(\Delta x) = \dfrac{A^2}{(\Delta x)^2} + B^2(\Delta x)^2 is 2AB2AB.

(c) Hence show that ⟨E⟩β‰₯12ℏω\langle E\rangle \ge \frac12\hbar\omega.

βœ“ Worked solution

(a) With ⟨x⟩=0\langle x\rangle = 0 we have ⟨x2⟩=(Ξ”x)2\langle x^2\rangle = (\Delta x)^2 exactly, and likewise ⟨p2⟩=(Ξ”p)2\langle p^2\rangle = (\Delta p)^2. The uncertainty principle then gives Ξ”pβ‰₯ℏ/2Ξ”x\Delta p \ge \hbar/2\Delta x, so

⟨p2⟩=(Ξ”p)2β‰₯ℏ24(Ξ”x)2\langle p^2\rangle = (\Delta p)^2 \ge \frac{\hbar^2}{4(\Delta x)^2}

Substituting both into ⟨E⟩\langle E\rangle produces the stated inequality.

(b) Differentiate: Fβ€²=βˆ’2A2/(Ξ”x)3+2B2Ξ”xF' = -2A^2/(\Delta x)^3 + 2B^2\Delta x, which vanishes when (Ξ”x)4=A2/B2(\Delta x)^4 = A^2/B^2, i.e. Ξ”x=A/B\Delta x = \sqrt{A/B}. At that point F=A2β‹…BA+B2β‹…AB=AB+AB=2ABF = A^2\cdot\frac{B}{A} + B^2\cdot\frac{A}{B} = AB + AB = 2AB.

(Without calculus: F=(AΞ”xβˆ’BΞ”x)2+2ABF = \left(\frac{A}{\Delta x} - B\Delta x\right)^2 + 2AB, and a square is never negative. The minimum is where the square vanishes.)

(c) Match: A2=ℏ2/8mA^2 = \hbar^2/8m and B2=12mΟ‰2B^2 = \frac12m\omega^2, so A=ℏ/8mA = \hbar/\sqrt{8m} and B=Ο‰m/2B = \omega\sqrt{m/2}. Then

2AB=2⋅ℏ8mβ‹…Ο‰m2=2ℏω116=ℏω22AB = 2\cdot\frac{\hbar}{\sqrt{8m}}\cdot\omega\sqrt{\frac m2} = 2\hbar\omega\sqrt{\frac{1}{16}} = \frac{\hbar\omega}{2}

so ⟨E⟩β‰₯12ℏω\langle E\rangle \ge \frac12\hbar\omega.

2 The ground state is a Gaussian, and how wide derivation

Using the properties of the Gaussian distribution from problem 2 of chapter 2, show that the position probability density of the ground state is a Gaussian distribution with standard deviation Οƒ=ℏ/2mΟ‰\sigma = \sqrt{\hbar/2m\omega}.

βœ“ Worked solution

From Table 6.1, ψ0(x)=(aΟ€)βˆ’1/2eβˆ’x2/2a2\psi_0(x) = (a\sqrt\pi)^{-1/2}e^{-x^2/2a^2} with a=ℏ/mΟ‰a = \sqrt{\hbar/m\omega}, so the probability density is

∣ψ0(x)∣2=1aΟ€eβˆ’x2/a2|\psi_0(x)|^2 = \frac{1}{a\sqrt\pi}e^{-x^2/a^2}

The standard Gaussian of standard deviation Οƒ\sigma is 1Οƒ2Ο€eβˆ’x2/2Οƒ2\dfrac{1}{\sigma\sqrt{2\pi}}e^{-x^2/2\sigma^2}. Comparing exponents, 12Οƒ2=1a2\dfrac{1}{2\sigma^2} = \dfrac{1}{a^2}, so

Οƒ=a2=ℏ2mΟ‰\sigma = \frac{a}{\sqrt2} = \sqrt{\frac{\hbar}{2m\omega}}

and the prefactors agree automatically, since both densities are normalized.

Note the factor of 2\sqrt2. The eigenfunction ψ0\psi_0 has β€œwidth” aa in the sense that it falls to 1/e1/\sqrt e there; the density ∣ψ0∣2|\psi_0|^2 is narrower by 2\sqrt2, because squaring a Gaussian halves its variance. That is also why Ξ”x=a/2\Delta x = a/\sqrt2 and not aa β€” the same 2\sqrt2 that made problem 1’s bound come out exactly.

3 How often is the particle where it cannot be? numerical

Find the amplitude AA of a classical particle with the same energy as a quantum particle in the oscillator ground state. Write down an expression for the probability of finding the quantum particle in the classically forbidden region ∣x∣>A|x| > A.

πŸ’‘ Phillips' own hint

Show that if a classical particle of amplitude AA has energy 12ℏω\frac12\hbar\omega then A=a=ℏ/mΟ‰A = a = \sqrt{\hbar/m\omega}. The probability of finding the quantum particle beyond it is then

2aΟ€βˆ«a∞eβˆ’x2/a2 dx\frac{2}{a\sqrt\pi}\int_a^\infty e^{-x^2/a^2}\,\mathrm dx
βœ“ Worked solution

The turning point. A classical oscillator of amplitude AA has E=12mΟ‰2A2E = \frac12m\omega^2A^2, by Eq. (6.5). Setting that equal to 12ℏω\frac12\hbar\omega:

A2=ℏmΟ‰=a2⟹A=aA^2 = \frac{\hbar}{m\omega} = a^2 \qquad\Longrightarrow\qquad A = a

The classical turning point is exactly one oscillator length β€” which is not a coincidence but the reason aa deserves its name.

The probability. Integrating ∣ψ0∣2|\psi_0|^2 over both forbidden tails:

P(∣x∣>a)=2aΟ€βˆ«a∞eβˆ’x2/a2 dx=erfc⁑(1)=0.157P(|x| > a) = \frac{2}{a\sqrt\pi}\int_a^\infty e^{-x^2/a^2}\,\mathrm dx = \operatorname{erfc}(1) = 0.157

About 16% β€” roughly one measurement in six finds the particle somewhere a classical particle of the same energy could not reach. This is not a rare leak like the tunnelling of Β§5.2a, where TT was 10βˆ’510^{-5}; it is routine, and it happens because the ground state’s own width is the same size as the classical turning point.

Group 2 β€” two variations on the potential

4 A spring that stretches but cannot compress conceptual

Consider

V(x)={∞ifΒ x<012mΟ‰2x2ifΒ x>0V(x) = \begin{cases}\infty & \text{if } x < 0\\ \frac12m\omega^2x^2 & \text{if } x > 0\end{cases}

Sketch the ground and first excited eigenfunctions, and give their energies.

πŸ’‘ Phillips' own hint

The potential is identical to the harmonic oscillator for 0<x<∞0 < x < \infty, but presents an infinite barrier for x<0x < 0. The energies are 32ℏω\frac32\hbar\omega, 72ℏω\frac72\hbar\omega, … because the eigenfunctions must satisfy: ψ=0\psi = 0 for x<0x < 0; ψ\psi identical to a harmonic oscillator eigenfunction for x>0x > 0; and ψ\psi continuous at x=0x = 0 β€” which is only possible if the x>0x > 0 piece is an oscillator eigenfunction with n=1,3,5,…n = 1, 3, 5, \dots

βœ“ Worked solution

For x>0x > 0 the equation is unchanged, so any solution must be an oscillator eigenfunction there. The infinite wall forces ψ(0)=0\psi(0) = 0, and only the odd eigenfunctions vanish at the origin β€” ψ1,ψ3,ψ5,…\psi_1, \psi_3, \psi_5, \dots

So the spectrum of the half-oscillator is the odd half of the full oscillator’s:

E=32ℏω,Β 72ℏω,Β 112ℏω, …E = \tfrac32\hbar\omega,\ \tfrac72\hbar\omega,\ \tfrac{11}{2}\hbar\omega,\ \dots

The ground state is 32ℏω\frac32\hbar\omega β€” three times the full oscillator’s zero-point energy β€” and the first excited state is 72ℏω\frac72\hbar\omega. The sketches are the right-hand halves of ψ1\psi_1 and ψ3\psi_3 in Fig. 6.2, with zero drawn to the left of the origin, and renormalized by 2\sqrt2 since only half of each remains.

Note the spacing: adjacent levels are now 2ℏω2\hbar\omega apart, not ℏω\hbar\omega. Removing half the well doubled the level spacing.

8 The two-dimensional oscillator, and a preview of angular momentum derivation

For a particle in a two-dimensional oscillator with

[βˆ’β„22m(βˆ‚2βˆ‚x2+βˆ‚2βˆ‚y2)+12mΟ‰2(x2+y2)]ψnx,ny=Enx,nyψnx,ny\left[-\frac{\hbar^2}{2m}\left(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}\right) + \frac12m\omega^2(x^2+y^2)\right]\psi_{n_x,n_y} = E_{n_x,n_y}\psi_{n_x,n_y}

(a) Verify that ψnx,ny(x,y)=ψnx(x)ψny(y)\psi_{n_x,n_y}(x,y) = \psi_{n_x}(x)\psi_{n_y}(y) is an eigenfunction with Enx,ny=(nx+ny+1)ℏωE_{n_x,n_y} = (n_x + n_y + 1)\hbar\omega.

(b) Draw an energy level diagram and indicate the degeneracies.

(c) By expressing ψ1,0\psi_{1,0} and ψ0,1\psi_{0,1} in plane polar coordinates (r,Ο•)(r,\phi), find ψa\psi_a and ψb\psi_b obeying βˆ’iβ„β€‰βˆ‚Οˆa/βˆ‚Ο•=+β„Οˆa-i\hbar\,\partial\psi_a/\partial\phi = +\hbar\psi_a and βˆ’iβ„β€‰βˆ‚Οˆb/βˆ‚Ο•=βˆ’β„Οˆb-i\hbar\,\partial\psi_b/\partial\phi = -\hbar\psi_b.

πŸ’‘ Phillips' own hint

For (a) and (b), modify the argument of §6.5. For (c), consider linear superpositions of the form ψ1,0±iψ0,1\psi_{1,0} \pm i\psi_{0,1}.

βœ“ Worked solution

(a) Identical to Β§6.5’s separation with one dimension removed: each H^x\hat H_x acts only on xx, so the product passes through and the energies add, giving (nx+12)+(ny+12)=nx+ny+1(n_x+\frac12) + (n_y+\frac12) = n_x + n_y + 1 in units of ℏω\hbar\omega.

(b) Since only N=nx+nyN = n_x + n_y enters, level NN has a degeneracy of N+1N+1 β€” the pairs (0,N),(1,Nβˆ’1),…,(N,0)(0,N), (1,N-1), \dots, (N,0). A linear ladder of degeneracies 1, 2, 3, 4, against the three-dimensional oscillator’s triangular 1, 3, 6, 10.

(c) From Table 6.1, ψ0∝eβˆ’x2/2a2\psi_0 \propto e^{-x^2/2a^2} and ψ1∝x eβˆ’x2/2a2\psi_1 \propto x\,e^{-x^2/2a^2}, so

ψ1,0∝x eβˆ’r2/2a2,ψ0,1∝y eβˆ’r2/2a2\psi_{1,0} \propto x\,e^{-r^2/2a^2},\qquad \psi_{0,1} \propto y\,e^{-r^2/2a^2}

using x2+y2=r2x^2 + y^2 = r^2. Then

ψa=ψ1,0+iψ0,1∝(x+iy)eβˆ’r2/2a2=r e+iΟ•eβˆ’r2/2a2\psi_a = \psi_{1,0} + i\psi_{0,1} \propto (x + iy)e^{-r^2/2a^2} = r\,e^{+i\phi}e^{-r^2/2a^2}

and ψb=ψ1,0βˆ’iψ0,1∝r eβˆ’iΟ•eβˆ’r2/2a2\psi_b = \psi_{1,0} - i\psi_{0,1} \propto r\,e^{-i\phi}e^{-r^2/2a^2}. The Ο•\phi dependence is now a bare eΒ±iΟ•e^{\pm i\phi}, so

βˆ’iβ„βˆ‚Οˆaβˆ‚Ο•=βˆ’iℏ(i)ψa=+β„Οˆa,βˆ’iβ„βˆ‚Οˆbβˆ‚Ο•=βˆ’β„Οˆb-i\hbar\frac{\partial\psi_a}{\partial\phi} = -i\hbar(i)\psi_a = +\hbar\psi_a,\qquad -i\hbar\frac{\partial\psi_b}{\partial\phi} = -\hbar\psi_b

as required.

Problem 8(b) β€” the two-dimensional ladder

energy (Δ§Ο‰)E = Δ§Ο‰1.00Γ—1E = 2Δ§Ο‰2.00Γ—2E = 3Δ§Ο‰3.00Γ—3E = 4Δ§Ο‰4.00Γ—4

Click a gold level for its energy and degeneracy.

No transitions shown in this view.

Degeneracy N + 1, growing linearly. Click a level for the (nβ‚“, n_y) pairs. Compare Fig. 6.5's three-dimensional version, where the counts are the triangular numbers instead.

Group 3 β€” states that move, and a molecule

5 Two levels, and the motion they produce derivation

At t=0t = 0 a particle in V=12mΟ‰2x2V = \frac12m\omega^2x^2 has Ξ¨(x,0)=12[ψ0(x)+ψ1(x)]\Psi(x,0) = \frac{1}{\sqrt2}[\psi_0(x) + \psi_1(x)], with ψ0\psi_0 and ψ1\psi_1 real, normalized and orthogonal.

(a) Write Ξ¨(x,t)\Psi(x,t). (b) Show it is normalized. (c) Show ∣Ψ∣2|\Psi|^2 oscillates at Ο‰\omega. (d) Show ⟨x⟩=Acos⁑ωt\langle x\rangle = A\cos\omega t with A=∫ψ0xψ1 dxA = \int\psi_0x\psi_1\,\mathrm dx.

πŸ’‘ Phillips' own hint

For (a), use Eq. (6.15). For (b), follow the steps that led to Eq. (4.53). For (c), use E1βˆ’E0=ℏωE_1 - E_0 = \hbar\omega.

βœ“ Worked solution

(a) Each term carries its own phase:

Ξ¨(x,t)=12[ψ0(x)eβˆ’iE0t/ℏ+ψ1(x)eβˆ’iE1t/ℏ]\Psi(x,t) = \frac{1}{\sqrt2}\left[\psi_0(x)e^{-iE_0t/\hbar} + \psi_1(x)e^{-iE_1t/\hbar}\right]

(b) ∣Ψ∣2=12[ψ02+ψ12+2ψ0ψ1cos⁑(Ο‰t)]|\Psi|^2 = \frac12[\psi_0^2 + \psi_1^2 + 2\psi_0\psi_1\cos(\omega t)] after using E1βˆ’E0=ℏωE_1 - E_0 = \hbar\omega. Integrating, the first two terms give 1 each and the cross term gives zero by orthogonality, so ∫∣Ψ∣2=12(1+1+0)=1\int|\Psi|^2 = \frac12(1+1+0) = 1 for all tt.

(c) That same expression has a single time-dependent term, ψ0ψ1cos⁑ωt\psi_0\psi_1\cos\omega t β€” so the density sloshes at exactly Ο‰\omega. The frequency is (E1βˆ’E0)/ℏ(E_1-E_0)/\hbar, and for this potential adjacent levels are always ℏω\hbar\omega apart.

(d) Multiplying by xx and integrating, the ψ02\psi_0^2 and ψ12\psi_1^2 terms vanish (both are even, xx is odd), leaving

⟨x⟩=(∫ψ0 xβ€‰Οˆ1 dx)cos⁑ωt\langle x\rangle = \left(\int\psi_0\,x\,\psi_1\,\mathrm dx\right)\cos\omega t

so A=∫ψ0xψ1 dxA = \int\psi_0x\psi_1\,\mathrm dx, which problem 11 evaluates as a/2=0.707aa/\sqrt2 = 0.707a.

6 Three levels, and two frequencies numerical

For Ξ¨(x,t)=13[ψ0eβˆ’iE0t/ℏ+ψ1eβˆ’iE1t/ℏ+ψ2eβˆ’iE2t/ℏ]\Psi(x,t) = \sqrt{\tfrac13}\left[\psi_0e^{-iE_0t/\hbar} + \psi_1e^{-iE_1t/\hbar} + \psi_2e^{-iE_2t/\hbar}\right] with real normalized ψn\psi_n:

(a) the expectation value of the energy; (b) the uncertainty in the energy; (c) show ∣Ψ∣2=A(x)+B(x)cos⁑ωt+C(x)cos⁑2Ο‰t|\Psi|^2 = A(x) + B(x)\cos\omega t + C(x)\cos 2\omega t and find AA, BB, CC.

πŸ’‘ Phillips' own hint

For (a) and (b), use Eq. (4.54). For (c), use eiΞΈ+eβˆ’iΞΈ=2cos⁑θe^{i\theta} + e^{-i\theta} = 2\cos\theta.

βœ“ Worked solution

(a) Each level has probability 13\frac13:

⟨E⟩=13(12+32+52)ℏω=32ℏω\langle E\rangle = \tfrac13\left(\tfrac12 + \tfrac32 + \tfrac52\right)\hbar\omega = \tfrac32\hbar\omega

(b) ⟨E2⟩=13(14+94+254)(ℏω)2=3512(ℏω)2\langle E^2\rangle = \frac13\left(\frac14+\frac94+\frac{25}4\right)(\hbar\omega)^2 = \frac{35}{12}(\hbar\omega)^2, so

Ξ”E=3512βˆ’94 ℏω=23 ℏω=0.816 ℏω\Delta E = \sqrt{\tfrac{35}{12} - \tfrac94}\,\hbar\omega = \sqrt{\tfrac23}\,\hbar\omega = 0.816\,\hbar\omega

(c) Expanding ∣Ψ∣2|\Psi|^2 gives one term per pair of levels, oscillating at the difference of their energies over ℏ\hbar. Pairs (0,1)(0,1) and (1,2)(1,2) are adjacent and beat at Ο‰\omega; the pair (0,2)(0,2) is two rungs apart and beats at 2Ο‰2\omega. So

A(x)=13(ψ02+ψ12+ψ22)A(x) = \tfrac13\left(\psi_0^2+\psi_1^2+\psi_2^2\right)B(x)=23(ψ0ψ1+ψ1ψ2),C(x)=23β€‰Οˆ0ψ2B(x) = \tfrac23\left(\psi_0\psi_1 + \psi_1\psi_2\right),\qquad C(x) = \tfrac23\,\psi_0\psi_2

BB has two terms and CC only one β€” because two of the three pairs are adjacent. And note that ⟨x⟩\langle x\rangle would keep only the BB terms: problem 11 shows ∫ψ0xψ2=0\int\psi_0x\psi_2 = 0, so the 2Ο‰2\omega component is visible in the density but not in the average position, exactly as Β§6.3 claimed.

7 The nitric oxide bond numerical

Transitions between adjacent vibrational levels of NO give infrared radiation of wavelength Ξ»=5.33Β ΞΌm\lambda = 5.33\ \mu\mathrm{m}. Find the elastic constant kk of the bond. The reduced mass of NO is 7.467.46 amu.

πŸ’‘ Phillips' own hint

Use Eq. (6.20) to show that k=1548Β N mβˆ’1k = 1548\ \mathrm{N\,m^{-1}}.

βœ“ Worked solution

Invert Eq. (6.20), exactly as Β§6.4 did for CO:

Ξ»=2Ο€cΞΌk⟹k=ΞΌ(2Ο€cΞ»)2\lambda = 2\pi c\sqrt{\frac{\mu}{k}} \qquad\Longrightarrow\qquad k = \mu\left(\frac{2\pi c}{\lambda}\right)^2

With ΞΌ=7.46Γ—1.6605Γ—10βˆ’27=1.239Γ—10βˆ’26Β kg\mu = 7.46 \times 1.6605\times10^{-27} = 1.239\times10^{-26}\ \mathrm{kg} and 2Ο€c/Ξ»=3.534Γ—1014Β sβˆ’12\pi c/\lambda = 3.534\times10^{14}\ \mathrm{s^{-1}}:

k=1547Β N mβˆ’1k = 1547\ \mathrm{N\,m^{-1}}

matching the hint’s 15481548.

Worth comparing. Β§6.4 found k=1907Β N mβˆ’1k = 1907\ \mathrm{N\,m^{-1}} for CO by the same method. NO’s bond is about 19% softer β€” and NO has one fewer bonding electron pair than CO, which is exactly the direction chemistry predicts. Two infrared wavelengths, two bond stiffnesses, and a chemical comparison that never required touching either molecule.

Group 4 β€” the ladder algebra

These three are marked β€œfor readers who studied Section 6.6”, and they supply three results the chapter used on credit about its raising and lowering operators : that the ladder goes down as well as up, what the step constants are, and the selection rule.

Apply the operator, land on the next rung β€” Eqs. (6.35)–(6.37)

-4-2024-1.0-0.500.51.0q = x / aψ (normalized)
  • ψ1 β€” what went in
  • ψ2 β€” what it should be
  • the operator’s output
polynomial in2qΞ΅ = 1.5
polynomial out4qΒ² βˆ’ 2Ξ΅ = 2.5
|aβ‚™|Β² = 2(n+1)4.0000expected 4

The solid curve is what the operator produced; the thick pale one is ψ2 drawn independently. They coincide, which is Eq. (6.37): applying [q βˆ’ d/dq] to a solution of energy Ξ΅ gives a solution of energy Ξ΅ + 1, so one state and one operator generate the whole spectrum.

Watch the polynomial row. Starting from 1 and raising repeatedly gives 2q, then 4qΒ² βˆ’ 2, then 8qΒ³ βˆ’ 12q β€” the Hermite polynomials, appearing on their own rather than being looked up. Those are exactly the brackets in Table 6.1. The operators here act on the polynomial alone β€” raising is p ↦ 2qΒ·p βˆ’ pβ€² and lowering is p ↦ pβ€² β€” so every coefficient shown is exact integer arithmetic, with no numerical differentiation anywhere.

Natural units Δ§ = m = Ο‰ = 1, so q = x/a and Ξ΅β‚™ = n + Β½. Curves are normalized for display; the |aβ‚™|Β² row uses the unnormalized norms, which is what problem 10 asks for.

9 The lowering operator derivation

Show that when [q+d/dq][q + \mathrm d/\mathrm dq] acts on ψn\psi_n with eigenvalue Ο΅n\epsilon_n, it gives ψnβˆ’1\psi_{n-1} with eigenvalue Ο΅nβˆ’1=Ο΅nβˆ’1\epsilon_{n-1} = \epsilon_n - 1.

βœ“ Worked solution

The mirror image of Β§6.6’s raising argument, with the two forms swapped. Start from Eq. (6.31) applied to the nn-th state:

[qβˆ’ddq][q+ddq]ψn=(2Ο΅nβˆ’1)ψn\left[q - \frac{\mathrm d}{\mathrm dq}\right]\left[q + \frac{\mathrm d}{\mathrm dq}\right]\psi_n = (2\epsilon_n - 1)\psi_n

Apply [q+d/dq][q + \mathrm d/\mathrm dq] to both sides:

[q+ddq][qβˆ’ddq]([q+ddq]ψn)=(2Ο΅nβˆ’1)([q+ddq]ψn)\left[q + \frac{\mathrm d}{\mathrm dq}\right]\left[q - \frac{\mathrm d}{\mathrm dq}\right]\left(\left[q + \frac{\mathrm d}{\mathrm dq}\right]\psi_n\right) = (2\epsilon_n - 1)\left(\left[q + \frac{\mathrm d}{\mathrm dq}\right]\psi_n\right)

Compare with Eq. (6.30), which reads [… ][… ]ψ=(2Ο΅+1)ψ[\dots][\dots]\psi = (2\epsilon + 1)\psi. Matching constants, 2Ο΅+1=2Ο΅nβˆ’12\epsilon + 1 = 2\epsilon_n - 1, so Ο΅=Ο΅nβˆ’1\epsilon = \epsilon_n - 1.

So [q+d/dq][q + \mathrm d/\mathrm dq] steps down one rung β€” except at n=0n = 0, where Eq. (6.40) gives zero and the ladder stops. Try it in the widget above.

10 The step constants derivation

If [qβˆ’d/dq]ψn=anψn+1[q - \mathrm d/\mathrm dq]\psi_n = a_n\psi_{n+1}, show that ψn\psi_n and ψn+1\psi_{n+1} have the same normalization if ∣an∣2=2(n+1)|a_n|^2 = 2(n+1). Similarly, if [q+d/dq]ψn=bnψnβˆ’1[q + \mathrm d/\mathrm dq]\psi_n = b_n\psi_{n-1}, show ∣bn∣2=2n|b_n|^2 = 2n.

πŸ’‘ Phillips' own hint

Write down the normalization integrals; note that integration by parts gives

βˆ«βˆ’βˆžβˆždfdqg dq=βˆ’βˆ«βˆ’βˆžβˆžfdgdq dq\int_{-\infty}^{\infty}\frac{\mathrm df}{\mathrm dq}g\,\mathrm dq = -\int_{-\infty}^{\infty}f\frac{\mathrm dg}{\mathrm dq}\,\mathrm dq

when fg→0fg \to 0 at infinity; and use Eq. (6.30).

βœ“ Worked solution

Take the norm of the raised state:

∣an∣2=∫([qβˆ’ddq]ψn)2dq|a_n|^2 = \int\left(\left[q-\tfrac{\mathrm d}{\mathrm dq}\right]\psi_n\right)^2\mathrm dq

Integration by parts moves one operator across, turning [qβˆ’d/dq][q - \mathrm d/\mathrm dq] into [q+d/dq][q + \mathrm d/\mathrm dq] β€” the sign flips because d/dq\mathrm d/\mathrm dq is anti-Hermitian while qq is Hermitian. So

∣an∣2=∫ψn[q+ddq][qβˆ’ddq]ψn dq|a_n|^2 = \int\psi_n\left[q+\tfrac{\mathrm d}{\mathrm dq}\right]\left[q-\tfrac{\mathrm d}{\mathrm dq}\right]\psi_n\,\mathrm dq

and the operator in the middle is exactly Eq. (6.30)β€˜s, which returns (2Ο΅n+1)ψn(2\epsilon_n + 1)\psi_n. With Ο΅n=n+12\epsilon_n = n + \frac12:

∣an∣2=2ϡn+1=2n+2=2(n+1)|a_n|^2 = 2\epsilon_n + 1 = 2n + 2 = 2(n+1)

The lowering case is identical using Eq. (6.31), giving ∣bn∣2=2Ο΅nβˆ’1=2n|b_n|^2 = 2\epsilon_n - 1 = 2n β€” and note it vanishes at n=0n = 0, which is Eq. (6.40) once more.

11 Matrix elements, and the selection rule derivation

(a) Show that ∫ψmβˆ—qβ€‰Οˆn dq=12(anΞ΄m,n+1+bnΞ΄m,nβˆ’1)\int\psi_m^*q\,\psi_n\,\mathrm dq = \frac12(a_n\delta_{m,n+1} + b_n\delta_{m,n-1}).

(b) Show ∫ψnβˆ—q2ψn dq=n+12\int\psi_n^*q^2\psi_n\,\mathrm dq = n + \frac12 and ∫ψnβˆ—(βˆ’d2/dq2)ψn dq=n+12\int\psi_n^*\left(-\mathrm d^2/\mathrm dq^2\right)\psi_n\,\mathrm dq = n + \frac12.

(c) Hence verify ⟨x⟩=0\langle x\rangle = 0, ⟨x2⟩=(n+12)a2\langle x^2\rangle = (n+\frac12)a^2, ⟨p⟩=0\langle p\rangle = 0 and ⟨p2⟩=(n+12)ℏ2/a2\langle p^2\rangle = (n+\frac12)\hbar^2/a^2.

βœ“ Worked solution

(a) The hint is the key identity: solve the two operator definitions for qq and for d/dq\mathrm d/\mathrm dq separately β€”

q=12[q+ddq]+12[qβˆ’ddq]q = \tfrac12\left[q+\tfrac{\mathrm d}{\mathrm dq}\right] + \tfrac12\left[q-\tfrac{\mathrm d}{\mathrm dq}\right]

So qψn=12(bnψnβˆ’1+anψn+1)q\psi_n = \frac12(b_n\psi_{n-1} + a_n\psi_{n+1}): multiplying by qq always produces a mixture of the two neighbouring states and nothing else. Taking the inner product with ψm\psi_m and using orthonormality, Eq. (6.41), leaves exactly the stated Kronecker deltas.

This is the selection rule. ∫ψmβˆ—qψn=0\int\psi_m^*q\psi_n = 0 unless ∣mβˆ’n∣=1|m - n| = 1 β€” not because of any special property of these functions, but because qq is a sum of one raising and one lowering operator.

(b) Apply qq twice. q2ψnq^2\psi_n reaches ψnβˆ’2\psi_{n-2}, ψn\psi_n and ψn+2\psi_{n+2}; only the middle survives the inner product with ψn\psi_n, and its coefficient is 14(anβˆ’1bn+bn+1an)=14(2n+2(n+1))=n+12\frac14(a_{n-1}b_n + b_{n+1}a_n) = \frac14(2n + 2(n+1)) = n + \frac12. The kinetic term follows the same way, or immediately from Eq. (6.29): since [βˆ’d2/dq2+q2][-\mathrm d^2/\mathrm dq^2 + q^2] has expectation 2Ο΅n=2n+12\epsilon_n = 2n+1 and the q2q^2 half is n+12n + \frac12, the other half must be n+12n + \frac12 too.

That equality is Β§6.3’s even split β€” half potential, half kinetic β€” now derived rather than quoted.

(c) Restore units with x=qax = qa and p=βˆ’iℏ d/dxp = -i\hbar\,\mathrm d/\mathrm dx. ⟨x⟩∝∫ψnqψn=0\langle x\rangle \propto \int\psi_nq\psi_n = 0 by part (a) with m=nm = n; ⟨x2⟩=(n+12)a2\langle x^2\rangle = (n+\frac12)a^2; and ⟨p2⟩=(ℏ2/a2)(n+12)\langle p^2\rangle = (\hbar^2/a^2)(n+\frac12). Hence

Ξ”x Δp=(n+12)ℏ\Delta x\,\Delta p = \left(n+\tfrac12\right)\hbar

which is the result Β§6.3 stated without proof.

Check yourself

0 / 6 answered

  1. 1.Problem 1 gets the exact ground-state energy from the uncertainty principle alone. The same argument applied to a box in chapter 4 was 39.5Γ— loose. What is different here?

  2. 2.Problem 3 finds the ground-state particle beyond the classical turning point about 16% of the time. Why is that so much larger than chapter 5's tunnelling probabilities?

  3. 3.Problem 4's half-oscillator β€” a spring that stretches but cannot compress β€” has ground-state energy . Where does that come from?

  4. 4.Problem 8(c) builds and finds states of definite angular momentum . What does that reveal about degenerate levels?

  5. 5.Problem 11(a) shows unless . What makes that true?

  6. 6.Problem 10 gives , and Β§6.6's bridge wrote the annihilation matrix with on its superdiagonal. How are those connected?