Physics 88 Problem Set 7 Due Wednesday /4/ 7)a)Consider the proton to be a uniformly charged sphere of radius f m Determine the correction to the s ground state energy 4 points) This is a standard problem One can calculate the electric filed from Gauss law as in a freshman course and deduce the potential Φ r) as a function of r The energy of the electron is eφ r) For r > a the potential is e 4πɛ r For r < a only the charge inside contributes to the electric field and so it is er3 /a 3 ) 4πɛ = r 4πɛ and so the potential contains a contribution 4πɛ e r a 3 apart from an additive contribution chosen to make the potential continuous So we have for r < a V r) = e r 4πɛ a 3 + 3 ) a The contribution to the Hamiltonian is eφr) and this differs form the unperturbed Hamiltonian only for r < a The difference is the perturbing Hamiltonian and is given by H = q r a a + 3 a ) r for r < a and vanishes for r > a Now we need to evaluate s H s Since the nucleus is much smaller than the ohr radius it makes eminent sense to approximate ψ s r) by ψ s ) which is given by where a πa 3 is the ohr radius Since the proton is not a spherical distribution with a sharp edge doing the integral carefully is unwarranted We just need Thus we obtain 4π E ) = a q 4 a 5 This is approximately 936 khz ) dr r H = q a πa3 5 a a = πq a 5 ) 3 9 ev 387 9 ev r a 3 This is not a theoretical problem only One can calculate the effect of the nuclear radius on the S state and check how the difference in energy between the S and S states is affected by this term Indeed it has been used to provide an estimate of the charge radius
of the proton! These incredibly precise experiments were conducted by the group of Ted Hänsch who shared the 5 Nobel Prize in Physics for precision spectroscopy The S S transition is especially slow as we will see next quarter b) Consider the nuclear shell model Neutrons and protons are treated as independent particles in a three-dimensional isotropic harmonic oscillator potential Each nucleon has spin-orbit coupling a l s with a > that is a perturbation on the confining harmonic oscillator potential Write down the three lowest energy single-particle states by specifying the correct quantum numbers You are not asked to write down the wave functions Determine the ground state angular momentum and parity if you can) of i) tritium 3 H, ii) 5 5 points) The energy levels of the three-dimensional oscillator can be found on page 35 of Shankar The lowest enrgy state has n l = The first excited state is characterized by n = l = and is three-fold degenerate The next state has n = and l = or and this has degeneracy of 6 Since a > the states in which l and s are parallel have the lowest energy The ground state has zero orbital angular momentum and no spin-orbit effects It can be denoted by S / The next state splits into a P 3/ and a P / level with the P 3/ level being lower in energy H so = a [jj + ) ll + ) ss + )] which has the value a[5/4) 3/4)] = a for P 3/ and a [3/4) 3/4)] = a for P / What is the center of gravity of these six states?) Thus the lowest three levels are S / with energy 3/) hω, P 3/ with energy 5/) hω a with degeneracy 4 and then the P / state with energy 5/) hω + a and degeneracy In tritium the two neutrons occupy the S / state and yield a net angular momentum of zero The one proton occupies the S / state and so J = / for this state The six neutrons completely fill the S / and P 3/ shells with a net angular momentum zero and the 5 protons fill the same two shells leaving a hole behind Thus, the angular momentum is J = 3/ If you are a nucleophile you should figure out the parity of these states Maria Goeppert-Mayer s Nobel lecture is a nice readable introduction to magic numbers in nuclei and the nuclear shell model
c) The dissociation energy of H is 446 ev You are given that the dissociation energy of the D molecule is 454 ev In the simplest approximation would you expect any difference in the dissociation energies for the two molecules Please justify your answer Now consider the vibrational motion of the nuclei and the associated zero-point energy to explain the difference Find the zero-point energy of the H molecule 5 points) If one treats the nuclei as infinitely massive in computing the energy of the molecule in the ground state we compute ɛ + R) the symmetric state for the two electrons) only the charge enters the Hamiltonian and therefore, both H and D yield the same dissociation energy In the adiabatic approximation we can interpret ɛ + R) as the potential energy of interaction between the two nuclei Expanding around the minimum of the potential that is assumed to occur at R we have ɛ + R) ɛ + R ) + k R R ) Thus the spring constant is the same for both the molecules The frequency of the associated oscillator will differ because of the difference in masses What is the energy of the diatomic molecule? The state is specified by n the excitation of the one-dimensional harmonic oscillator corresponding to the vibrational motion and J the rotational quantum number: E = ɛ + R ) + hω n + ) + JJ + ) + higher order terms In the ground state n = and J = The dissociation energy is the difference in energy between the ground state energy of the molecule and that of two non-interacting atoms R ) taken to be zero Thus the dissociation energy is given by E dis = ɛ + R ) hω ɛ + R ) < and so the dissociation energy is positive It is reduced by the zero-point fluctuations Since ɛ + R ) is the same for the hydrogen and deuterium molecules, we have E dis,d E dis,h = hω H hω D ) Since ω = k/µ where k is the spring constant and µ is the reduced mass and k is the same for both the molecules we have ω D ω H =, 3
Thus, given the difference in dissociation energies to be 8 ev we obtain the zero-point energy of hydrogen: hω H ) = 8 ev hω H = 7 ev 4
7 The uneperturbed Hamiltonian is γ S z with eigenvectors + and with energies E + = γ h/ and E = γ h/ respectively using obvious notation The perturbation is H = γσ x / The first-order corrections + σ x + and σ x vanish The second-order corrections are given by It is easy to see that E + = H + E + E E = γ h 4 The corrections to the eigenstates are easily calculated: + = H + E + E = + + H E E + = γ h 4 = = + It is easy to find the exact eigenvalues and eigenvectors by diagonalizing H = γ h ) So the eigenvalues are Λ ± = γ h + Expansion yields γ h ) + The unnormalized eigenvector for the ground state can be written down by inspection: ) ) + This agrees with the perturbation result The other eigenvector yields ) ) + + 74 The Thomas-Reiche-Kuhn sum rule is not only useful in spectroscopy we will discuss oscillator strengths next quarter) but has inspired many sum rules in condensed 5
matter and particle physics One usually starts from a commutator of operators or currents In our case, the usual starting point would be i n [X, P ] n = use completeness and P = im [H, X] to obtain the sum rule We start from the given expression We write down the following obvious steps: E E n ) X n = E E n ) n X X n = [ n XH X n n HX X n ] where we have used n H = n E n etc Completeness of the set { } allows one to simplify the right-hand side: E E n ) X n = n XHX HX n ) Since n HX n = E n n X n = n X H n we can write E E n ) X n = n XHX HX X H n Now we note that XHX HX X H = [X, [H, X]] as is verified by expansion Using [H, X] = [ ] P m, X = i h m P we obtain This yields the celebrated sum rule [X, [H, X]] = i h m [X, P ] = h m The other part is elementary algebra using h X = mω a + a ) and h X n = mω [ n δ,n + n + δ,n+ ] Thus we obtain X n = h mω [n δ,n + n + ) δ,n+ ] 6
For the first term on the right-hand side above we have E E n = hω and from the second term E E n = hωthus we obtain E E n ) X n = h h hω n + n + )) = mω m 77 Caveat: I might have pulled a fast one over you by doing an even number of sign errors or canceling errors since I knew the answer The algebra involved in the problem is trivial once you guess Ω I wanted to make sure you read Method on page 46 and understood this clever technique H = qfx Since we need [Ω, H ] X it is natural to guess Ω = AP and determine A We have [AP, H ] = mω A [P, X ] = i hmω A X qf So we choose A = i to obtain [Ω, H hmω ] = H Now that we have Ω we need the commutators in 739 Since this is a simple harmonic oscillator we employ creation and annihilation operators Using h X = mω a + m hω a ) and P = i a a) we obtain the different matrix elements: n Ω n n P n = n H Ω n = qf mω n a + a)a a) n = qf mω Substituting into Equation 739 the second term vanishes and we obtain the same result as in equation 76 73 The dipole selection rule according to Shankar states that l = l ± and thus in the Hamiltonian for the Stark effect for the n = states of hydrogen all the matrix elements with l = l vanish So only matrix elements between s and p states can be non-vanishing 7
From 7 we also need m = for the matrix elements of Z Thus the only independent, possibly and actually!) non-vanishing matrix element is H This determines its Hermitian conjugate, of course We need eer cos θ = π dφ π dθ sin θ Y cos θ Y drr R r) r R r) where all the functions are known and the integrals can be evaluated The wave functions written as a product of R nl r) and Yl m θ, φ) are given by ψ,, = ψ,, = ) / a 3 r ) a ) / r 4a 3 e r a a e r a 4π ) 3 cos θ 3) 4π Let us assume you wimp out and go to Mathematica The radial integral are easily evaluated in units of a, u r/a H = eea u cos θ): In[8]:= R[u_]:=/Sqrt[]) - u/) Exp[-u/] In[9]:= R[u_]:=Sqrt[/4] u Exp[-u/] In[]:= Integrate[R[u] u^3 R[u], {u,,infinity}] Out[]= -3 Sqrt[3] The angular integral, including a factor of π for the φ integral In[4]:= Y[theta_]:= Sqrt[3/4 Pi)] Cos[theta] In[5]:= Y[theta_]:= /Sqrt[4 Pi] In[8]:= Integrate[ Pi Sin[theta] Y[theta] Cos[theta] Y[theta], {theta,, Pi}] Out[8]= /Sqrt[3] Multiplying the values of the two integrals and restoring the factor eea we obtain eer cos θ = 3 eea 4) The Hamiltonian H within this subspace is given by 3eEa Ĥ 3eEa = 5) 8
So the state and are unaffected while and are mixed; we can diagonalize H in this subspace and find the corrections to the eigenvalue E as E ) = ± 3eEa 6) The eigenvectors are [ ± ] 7) Thus the fourfold degeneracy is partially lifted and there is an energy change linear in E Note that the splitting could be guessed in order of magnitude as eea on dimensional grounds 9