Problems 9

Part V ✎ Problems Phillips pp. 205–212 · ~31 min read

  • reduced mass
  • fine structure

Change one constant and the whole atom rescales. Four of these problems do nothing but that, and between them they cover positronium, muonic hydrogen and every hydrogen-like ion.

Fifteen problems, the most in the book, and they fall into four groups. 1–3 find the atom’s scale by three independent routes. 4–7 are the bookkeeping — normalization, degeneracy, mean radii, the virial theorem. 8–10 compute things the chapter only asserted: the momentum-space wave function, a dipole matrix element, and the spin–orbit energy. And 11–15 are the payoff: change one mass or one charge and the whole atom rescales.

Finding the scale

1 A variational calculation that finds the Bohr radius numerical

An electron in a Coulomb potential has the trial wave function ψ(r)=Ner/a\psi(r) = Ne^{-r/a} with aa a constant.

(a) What is its orbital angular momentum? (b) Show that V=e2/4πϵ0a\langle V\rangle = -e^2/4\pi\epsilon_0a and T=2/2mea2\langle T\rangle = \hbar^2/2m_ea^2. (c) Show that E\langle E\rangle is minimized when a=a0a = a_0, and find that minimum.

💡 Phillips' own hint

(a) A spherically symmetric wave function implies zero orbital angular momentum. (c) Find the minimum of E=T+V\langle E\rangle = \langle T\rangle + \langle V\rangle by setting dE/da=0\mathrm d\langle E\rangle/\mathrm da = 0. The integral 0rkeαrdr=k!/αk+1\int_0^\infty r^ke^{-\alpha r}\,\mathrm dr = k!/\alpha^{k+1} is useful here and in problems 3, 4 and 10.

✓ Worked solution

(a) ψ\psi depends only on rr, so §8.3a gives L^ψ=0\hat{\mathbf L}\psi = 0: zero, in every component and in magnitude.

(b) Both are standard integrals with the hint’s formula. The kinetic one needs the radial Laplacian; both give the stated results.

(c) Minimizing

E=22mea2e24πϵ0a\langle E\rangle = \frac{\hbar^2}{2m_ea^2} - \frac{e^2}{4\pi\epsilon_0a}

gives dE/da=2/mea3+e2/4πϵ0a2=0\mathrm d\langle E\rangle/\mathrm da = -\hbar^2/m_ea^3 + e^2/4\pi\epsilon_0a^2 = 0, so

a=4πϵ02mee2=a0andEmin=e28πϵ0a0=ER=13.6 eVa = \frac{4\pi\epsilon_0\hbar^2}{m_ee^2} = a_0 \qquad\text{and}\qquad \langle E\rangle_{\min} = -\frac{e^2}{8\pi\epsilon_0a_0} = -E_R = -13.6\ \mathrm{eV}

This is the whole ground state, obtained without solving anything. The trial function happens to be exact here — but the method works even when it is not, and it always gives an upper bound on the true ground-state energy.

2 The classical circular orbit derivation

For a classical particle in a Coulomb potential with effective potential Ve(r)V_e(r) from Eq. (9.3):

(a) Show that VeV_e has a minimum at rc=(L2/m)(4πϵ0/e2)r_c = (L^2/m)(4\pi\epsilon_0/e^2) with value L2/2mrc2-L^2/2mr_c^2, and hence that a particle of angular momentum LL has its lowest energy in a circular orbit of that radius. (b) Show that for energy E=e2/8πϵ0aE = -e^2/8\pi\epsilon_0a the orbit ranges between aa2arca - \sqrt{a^2 - ar_c} and a+a2arca + \sqrt{a^2 - ar_c}.

💡 Phillips' own hint

(a) Use dVe/dr=L2/mr3+e2/4πϵ0r2=0\mathrm dV_e/\mathrm dr = -L^2/mr^3 + e^2/4\pi\epsilon_0r^2 = 0. (b) The turning points, where pr=0p_r = 0, satisfy e2/8πϵ0a=L2/2mr2e2/4πϵ0re^2/8\pi\epsilon_0a = L^2/2mr^2 - e^2/4\pi\epsilon_0r.

✓ Worked solution

(a) Setting the derivative to zero gives rcr_c directly. At that radius the particle sits at the bottom of the effective potential with pr=0p_r = 0 and no way to move radially — a circular orbit, and the lowest energy Eq. (9.2) allows for that LL.

(b) Setting pr=0p_r = 0 in Eq. (9.2) gives a quadratic in rr whose roots are the two stated turning points — perihelion and aphelion. When a=rca = r_c the discriminant vanishes and the two coincide: the circular orbit of part (a).

Compare with quantum mechanics. rc=l(l+1)a0r_c = l(l+1)a_0 is Eq. (9.18)‘s minimum, and problem 3(e) shows the quantum state approaches it as ll \to \infty. But the quantum ground state at l=0l = 0 has no classical counterpart at all — a classical electron with zero angular momentum falls straight in.

3 The lowest state for each l, and its classical limit derivation

For u0,l(r)=Nrl+1er/(l+1)a0u_{0,l}(r) = Nr^{l+1}e^{-r/(l+1)a_0}:

(a) Show N2=[2(l+1)a0]2l+31(2l+2)!N^2 = \left[\dfrac{2}{(l+1)a_0}\right]^{2l+3}\dfrac{1}{(2l+2)!}. (b) Show the most probable radius is (l+1)2a0(l+1)^2a_0. (c) Show r=12(2l+3)(l+1)a0\langle r\rangle = \tfrac12(2l+3)(l+1)a_0. (d) Show r2=14(2l+4)(2l+3)(l+1)2a02\langle r^2\rangle = \tfrac14(2l+4)(2l+3)(l+1)^2a_0^2, and that Δr\Delta r becomes small compared with r\langle r\rangle as ll \to \infty. (e) Show that as ll \to \infty both the most probable and the mean radius tend to rcr_c from problem 2.

💡 Phillips' own hint

Use problem 1’s integral throughout. (b) Maximize r2l+2e2r/(l+1)a0r^{2l+2}e^{-2r/(l+1)a_0}. (e) For l1l \gg 1 the angular momentum tends to ll\hbar, and both radii tend to l2a0=L2a0/2l^2a_0 = L^2a_0/\hbar^2.

✓ Worked solution

Parts (a)–(d) are the hint’s integral applied four times.

(e) is the point of the problem. As ll grows:

  • the most probable radius (l+1)2a0l2a0(l+1)^2a_0 \to l^2a_0;
  • the mean radius 12(2l+3)(l+1)a0l2a0\tfrac12(2l+3)(l+1)a_0 \to l^2a_0;
  • and problem 2’s classical rc=l(l+1)a0l2a0r_c = l(l+1)a_0 \to l^2a_0.

All three converge. Meanwhile part (d) gives Δr/r0\Delta r/\langle r\rangle \to 0: the state stops being spread out and becomes a sharp shell at a definite radius.

That is the correspondence principle made concrete. A state of large ll with no radial nodes is a quantum object that behaves like a classical circular orbit — sharp radius, definite energy, and the classical radius to boot. §8.1’s cone picture said the same thing about direction; this says it about position.

Bookkeeping

4 Normalizing the ground state, and an orthogonality check derivation

The hydrogen ground state is ψ1(r)=N1er/a0\psi_1(r) = N_1e^{-r/a_0}.

(a) Normalize it using 0ψ124πr2dr=1\int_0^\infty|\psi_1|^24\pi r^2\,\mathrm dr = 1. (b) Show that 0ψ2(r)ψ1(r)r2dr=0\int_0^\infty\psi_2^*(r)\psi_1(r)r^2\,\mathrm dr = 0 for ψ2=N2(1+λr)er/2a0\psi_2 = N_2(1 + \lambda r)e^{-r/2a_0} if λ=1/2a0\lambda = -1/2a_0.

💡 Phillips' own hint

(a) The eigenfunction is normalized if N1=1/πa03N_1 = 1/\sqrt{\pi a_0^3}. (b) Show the integral vanishes if λ=1/2a0\lambda = -1/2a_0.

✓ Worked solution

(a) 0e2r/a04πr2dr=4π2(2/a0)3=πa03\int_0^\infty e^{-2r/a_0}4\pi r^2\mathrm dr = 4\pi\cdot\dfrac{2}{(2/a_0)^3} = \pi a_0^3, so N1=1/πa03N_1 = 1/\sqrt{\pi a_0^3}.

(b) The integral is a sum of two of problem 1’s integrals with α=3/2a0\alpha = 3/2a_0; requiring it to vanish fixes λ=1/2a0\lambda = -1/2a_0, which is exactly Table 9.1’s 2s state.

What this shows is not arithmetic. Eigenfunctions of a Hermitian operator with different eigenvalues are automatically orthogonal — §7.1 proved it. So this integral had to vanish, and the calculation is really determining λ\lambda: the 2s state is forced into its form by orthogonality to 1s, not by anything about the potential.

5 Proving the degeneracy is n² by induction derivation

Verify that gn=n2g_n = n^2 counts the independent hydrogen states of energy En=ER/n2E_n = -E_R/n^2, given that ll runs 0n10 \dots n-1 and each ll carries 2l+12l+1 values of mlm_l.

(a) Verify gn=n2g_n = n^2 for n=1n = 1 and 22. (b) Show that if it holds for n=kn = k then it holds for n=k+1n = k+1.

✓ Worked solution

(a) n=1n = 1: only l=0l = 0, giving 1 state =12= 1^2. n=2n = 2: l=0l = 0 gives 1 and l=1l = 1 gives 3, total 4=224 = 2^2.

(b) Going from kk to k+1k+1 adds exactly one new value of ll, namely l=kl = k, contributing 2k+12k+1 states. So

gk+1=gk+(2k+1)=k2+2k+1=(k+1)2g_{k+1} = g_k + (2k+1) = k^2 + 2k + 1 = (k+1)^2

The induction works because (k+1)2k2=2k+1(k+1)^2 - k^2 = 2k+1 is precisely the size of the multiplet being added — an identity about squares doing the work of a physical argument.

Doubling for spin gives 2n22n^2, which is the shell capacity 2, 8, 18, 32 that chapter 11 builds the periodic table from.

6 Mean radii in terms of n numerical

Rewrite Eq. (9.25) in terms of nn and ll, and find r\langle r\rangle for the 1s, 2s, 2p, 3s, 3p and 3d states.

✓ Worked solution

Substituting nr=nl1n_r = n - l - 1 into Eq. (9.25) collapses it to

r=a02[3n2l(l+1)]\langle r\rangle = \frac{a_0}{2}\left[3n^2 - l(l+1)\right]

which §9.3 verified is algebraically identical to the book’s form.

state$n$$l$$\langle r\rangle / a_0$
1s10
2s20
2p21
3s30
3p31
3d32

Mean radii from the collapsed formula. Click a cell for detail.

The pattern is 3n23n^2 dominating, so size grows quadratically with nn while ll makes a smaller correction downward.

7 The virial theorem in hydrogen derivation

Use the virial theorem from ch07’s problem 8 to find T\langle T\rangle and V\langle V\rangle for a hydrogen state of principal quantum number nn.

✓ Worked solution

Chapter 7’s problem 8 gave 2T=rdV/dr2\langle T\rangle = \langle r\,\mathrm dV/\mathrm dr\rangle, and for V=k/rV = -k/r, rdV/dr=+k/r=Vr\,\mathrm dV/\mathrm dr = +k/r = -V. So 2T=V2\langle T\rangle = -\langle V\rangle. Combined with En=T+VE_n = \langle T\rangle + \langle V\rangle:

T=En=+ERn2,V=2En=2ERn2\langle T\rangle = -E_n = +\frac{E_R}{n^2},\qquad \langle V\rangle = 2E_n = -\frac{2E_R}{n^2}

For the ground state: T=+13.6\langle T\rangle = +13.6 eV and V=27.2\langle V\rangle = -27.2 eV.

The consequence is worth carrying: E=TE = -\langle T\rangle, so a hydrogen atom that loses energy has its electron speed up. That is the same algebra that gives self-gravitating systems a negative heat capacity, and ch07’s problem 8 drew the same conclusion for satellites and stars.

Computing what the chapter asserted

8 The ground state in momentum space derivation

Using the 3-D Fourier transform, show that the ground state ψ(r)=er/a0/πa03\psi(r) = e^{-r/a_0}/\sqrt{\pi a_0^3} has momentum amplitude

ψ~(p)=22p05/2π(p2+p02)2,p0=/a0\tilde\psi(p) = \frac{2\sqrt2\,p_0^{5/2}}{\pi(p^2+p_0^2)^2},\qquad p_0 = \hbar/a_0

Hence show the most probable p|p| is p0/3p_0/\sqrt3 and the mean is 8p0/3π8p_0/3\pi.

💡 Phillips' own hint

Choose k\mathbf k along zz so eipr/=eiprcosθ/e^{-i\mathbf p\cdot\mathbf r/\hbar} = e^{-ipr\cos\theta/\hbar}, and write d3r=r2drd(cosθ)dϕ\mathrm d^3\mathbf r = r^2\,\mathrm dr\,\mathrm d(\cos\theta)\,\mathrm d\phi. Integrate ϕ\phi from 0 to 2π2\pi, cosθ\cos\theta from 1-1 to 1, and rr from 0 to \infty. For the most probable momentum locate the maximum of 4πp2ψ~(p)24\pi p^2|\tilde\psi(p)|^2; for the mean evaluate p=ψ~(p)pψ~(p)4πp2dp\langle p\rangle = \int\tilde\psi^*(p)\,p\,\tilde\psi(p)\,4\pi p^2\,\mathrm dp.

✓ Worked solution

The angular integrals give 4πsin(pr/)/pr4\pi\hbar\sin(pr/\hbar)/pr, leaving a radial integral that is elementary, and the stated amplitude follows.

The distribution has no sharp cutoff. ψ~2|\tilde\psi|^2 falls off as p8p^{-8}, so arbitrarily large momenta occur with small probability — the electron is localized to a0\sim a_0, so its momentum must be spread. That is the uncertainty principle appearing as a fact about the shape of a function rather than as an inequality.

And the peak is not the mean again: 1/3=0.5771/\sqrt3 = 0.577 against 8/3π=0.8498/3\pi = 0.849. The momentum distribution is skewed for exactly the reason the radial one is — §9.3’s caution applies verbatim.

9 Two dipole matrix elements, one of which vanishes derivation

Show that the integral of Eq. (9.27) is zero for a 2s1s2s \to 1s transition but not for 2p1s2p \to 1s.

💡 Phillips' own hint

The simplest way is to write out the integral in Cartesian coordinates.

✓ Worked solution

2s1s2s \to 1s. Both states are spherically symmetric, so the integrand is (even) × (odd) × (even) = odd, and integrating an odd function over all space gives zero. This is §9.4’s parity argument on its simplest possible case.

2p1s2p \to 1s. Now ψ2pz\psi_{2p} \propto z (for ml=0m_l = 0), and the operator contributes another zz, so the integrand goes as z2z^2even, and the integral survives.

The two cases differ only in parity, which is the entire content of the selection rule. And the consequence is the metastability of 2s that §9.4 computed: 0.14 s against 1.6 ns.

10 Spin–orbit energy — and the book contradicting itself numerical

Show that the spin–orbit energy of Eq. (9.35) has expectation values

Emag2p3/2=196α4mec2,Emag2p1/2=296α4mec2\langle E_{\mathrm{mag}}\rangle_{2p_{3/2}} = \frac{1}{96}\alpha^4m_ec^2, \qquad \langle E_{\mathrm{mag}}\rangle_{2p_{1/2}} = -\frac{2}{96}\alpha^4m_ec^2

Find their difference, and show the wavelengths of 2p3/21s1/22p_{3/2}\to1s_{1/2} and 2p1/21s1/22p_{1/2}\to1s_{1/2} differ by 5.4×1045.4\times10^{-4} nm.

💡 Phillips' own hint

Use Table 9.1’s eigenfunctions to find 1/r3\langle 1/r^3\rangle, and note that J2L2S2=j(j+1)l(l+1)s(s+1)\langle J^2 - L^2 - S^2\rangle = j(j+1) - l(l+1) - s(s+1) in units of 2\hbar^2.

✓ Worked solution

With l=1l = 1, s=12s = \tfrac12, the bracket gives +1+1 for j=32j = \tfrac32 and 2-2 for j=12j = \tfrac12the 1 : −2 ratio in the stated answers. Multiplying by 1/r3\langle 1/r^3\rangle for the 2p state gives the coefficients.

The difference is

196(296)=396=132α4mec2=45.3 μeV\frac{1}{96} - \left(-\frac{2}{96}\right) = \frac{3}{96} = \frac{1}{32}\,\alpha^4m_ec^2 = 45.3\ \mu\mathrm{eV}

and dividing that into the 10.2 eV transition energy shifts the 121.5 nm line by 5.4×1045.4\times10^{-4} nm — the book’s answer, confirmed.

Change one number, rescale the atom

11 Positronium derivation

Positronium is an electron bound to a positron of equal mass. (a) What is the reduced mass? (b) Write down its energy levels. (c) Given that hydrogen’s ground state has r=32a0\langle r\rangle = \tfrac32a_0, what is positronium’s?

💡 Phillips' own hint

See §9.5.

✓ Worked solution

(a) μ=mememe+me=me2\mu = \dfrac{m_em_e}{m_e+m_e} = \dfrac{m_e}{2} — exactly half, the most extreme case possible.

(b) En=μme(ERn2)=ER2n2E_n = \dfrac{\mu}{m_e}\left(-\dfrac{E_R}{n^2}\right) = -\dfrac{E_R}{2n^2}, so the ground state is 6.8-6.8 eV rather than 13.6-13.6 eV.

(c) Since a0=(me/μ)a0=2a0a_0' = (m_e/\mu)a_0 = 2a_0, the mean radius is 32×2a0=3a0\tfrac32 \times 2a_0 = \mathbf{3a_0}.

Half the binding, twice the size. And note what positronium is not: there is no heavy centre, so “the electron orbits the positron” is exactly as wrong as the reverse. The reduced mass handles that automatically.

12 Muonic hydrogen numerical

A muon (207me207m_e) bound to a proton. (a) Reduced mass? (b) Energy levels? (c) Mean ground-state radius?

✓ Worked solution

(a) μ=207memp207me+mp=186me\mu = \dfrac{207m_e\,m_p}{207m_e + m_p} = 186\,m_e — not 207, because the proton is only about 9 times heavier than the muon and recoils appreciably.

(b) En=186ER/n2E_n = -186\,E_R/n^2, so the ground state is 2530-2530 eV — an X-ray transition rather than an ultraviolet one.

(c) r=32a0/186=0.0081a0=2.8×1013\langle r\rangle = \tfrac32a_0/186 = 0.0081\,a_0 = 2.8\times10^{-13} m.

A muonic atom is 186 times smaller and 186 times more tightly bound. The muon orbits so close that the proton is no longer a good point charge — which is what problem 13 is about, and it is why muonic hydrogen is used to measure the proton’s radius.

13 How often is the electron inside the proton? numerical

Show that for Ra0R \ll a_0 a 1s electron has probability about 4(R/a0)34(R/a_0)^3 of being found within RR of the nucleus. Estimate that probability for a proton of radius 2×10152\times10^{-15} m, do the same for a muon in muonic hydrogen, and explain why muonic energy levels are more sensitive to the proton’s size.

✓ Worked solution

Using the corrected 43(R/a0)3\tfrac43(R/a_0)^3:

systemeffective a0a_0P(r<Rproton)P(r < R_{\text{proton}})
ordinary hydrogen5.29×10115.29\times10^{-11} m7×10147\times10^{-14}
muonic hydrogen2.85×10132.85\times10^{-13} m5×1075\times10^{-7}

A factor of about 6×1066\times10^6, because the probability goes as a03a_0^{-3} and the muonic radius is 186 times smaller: 18636×106186^3 \approx 6\times10^6.

That is why muonic atoms measure nuclear size. The energy shift from the proton not being a point charge is proportional to how much of the wave function sits inside it. For an electron that is a part in 101310^{13} — unmeasurable. For a muon it is a part in 10610^6, which is a large effect in modern spectroscopy, and muonic hydrogen is how the proton radius is measured today.

14 Hydrogen-like ions numerical

For a point charge ZeZe: write down EnE_n, write down the 1s, 2s and 2p radial eigenfunctions by inspection of Table 9.1, give the ionization energies of He⁺ and Li²⁺, and explain why relativistic corrections matter more for them.

💡 Phillips' own hint

When an electron is in the Coulomb potential of a point charge ZeZe, the size of a bound state is proportional to 1/Z1/Z and the binding energy to Z2Z^2.

✓ Worked solution

Energies and eigenfunctions. Everywhere a0a_0 appears, replace it by a0/Za_0/Z; everywhere ERE_R appears, replace it by Z2ERZ^2E_R:

En=Z2ERn2,u0,0(r)=2a0/Z(Zra0)eZr/a0E_n = -\frac{Z^2E_R}{n^2},\qquad u_{0,0}(r) = \frac{2}{\sqrt{a_0/Z}}\left(\frac{Zr}{a_0}\right)e^{-Zr/a_0}

and similarly for 2s and 2p.

Ionization energies. He⁺ has Z=2Z = 2: 4×13.6=54.44 \times 13.6 = \mathbf{54.4} eV. Li²⁺ has Z=3Z = 3: 9×13.6=122.59 \times 13.6 = \mathbf{122.5} eV.

Why relativity matters more. The electron’s speed scales as ZαcZ\alpha c, so the relativistic correction — which goes as (v/c)4(v/c)^4 — scales as Z4α4Z^4\alpha^4. Going from hydrogen to Li²⁺ multiplies it by 34=813^4 = 81.

For a heavy element like uranium (Z=92Z = 92), Zα0.67Z\alpha \approx 0.67: the inner electrons are genuinely relativistic, and no non-relativistic treatment works at all. This is why gold is yellow and mercury is liquid — both are relativistic effects on the outer electrons of heavy atoms.

15 Tritium decays — does the electron notice? numerical

A tritium atom’s nucleus beta-decays to ³He, so the electron suddenly finds itself in a Z=2Z = 2 potential. Starting in tritium’s ground state, show that

P=128a06[0r2e3r/a0dr]2P = \frac{128}{a_0^6}\left[\int_0^\infty r^2e^{-3r/a_0}\,\mathrm dr\right]^2

is the probability of ending in the He⁺ ground state, and verify P=0.702P = 0.702.

✓ Worked solution

The physics is ch04’s problem 8. The decay is far faster than the electron can respond, so its wave function is unchanged at that instant — but the Hamiltonian is not. The probability of finding it in a particular new eigenstate is the modulus squared of the overlap:

P=ψ1sHe+(r)ψ1sH(r)d3r2P = \left|\int\psi^{\mathrm{He}^+}_{1s}(\mathbf r)^*\,\psi^{\mathrm H}_{1s}(\mathbf r)\,\mathrm d^3\mathbf r\right|^2

Using ψ1sHer/a0\psi^{\mathrm H}_{1s} \propto e^{-r/a_0} and ψ1sHe+e2r/a0\psi^{\mathrm{He}^+}_{1s} \propto e^{-2r/a_0} (problem 14), the exponentials combine to e3r/a0e^{-3r/a_0} and the constants give the stated prefactor.

The integral is 2/(3/a0)3=227a032/(3/a_0)^3 = \tfrac{2}{27}a_0^3, so

P=128(227)2=512729=0.702P = 128\left(\frac{2}{27}\right)^2 = \frac{512}{729} = \mathbf{0.702}

So 70% of the time the electron lands in the new ground state, and 30% of the time it does not — ending up excited, or ionized outright. Nothing pushed it; the Hamiltonian simply changed underneath a wave function that had no time to adjust.

Angular shapes: |Yl,m(θ,φ)|² for l ≤ 3
peak 9.0a₀⟨r⟩ 10.5a₀r (Bohr radii)|u(r)|²
state 3d, n_r = 0
⟨r⟩ = 10.50 a₀, peak at 9.02 a₀
nodes: 2 angular, 0 radial

m = 0 gives the most nodal circles2 of them, the maximum for this l. With no angular momentum about z, the density is free to pile up at the poles. Note there is no dependence on φ, whatever l and m — because |eimφ|² = 1. A definite Lz means the wave function is completely smeared around the z axis.

Check yourself

0 / 6 answered

  1. 1.Problem 1 minimizes over a trial width and lands exactly on and eV. Why is that not just luck?

  2. 2.Problem 3(e) shows that as the mean radius, the most probable radius, and the classical circular-orbit radius all converge. What is that an instance of?

  3. 3.Problem 10 computes the splitting as . Why does that matter beyond the problem itself?

  4. 4.Problems 11–14 change one mass or one charge and rewrite the whole atom. What makes that legitimate?

  5. 5.Problem 13 finds a muon is about times more likely than an electron to be found inside the proton. Why?

  6. 6.In problem 15 a tritium nucleus decays and 70% of the time the electron ends up in the He⁺ ground state. What kind of calculation is that?