/*! This file is auto-generated */ .wp-block-button__link{color:#fff;background-color:#32373c;border-radius:9999px;box-shadow:none;text-decoration:none;padding:calc(.667em + 2px) calc(1.333em + 2px);font-size:1.125em}.wp-block-file__button{background:#32373c;color:#fff;text-decoration:none} Q18P In Problem 7.7 we found that t... [FREE SOLUTION] | 91Ó°ÊÓ

91Ó°ÊÓ

InProblem 7.7we found that the trial wave function with shielding (Equation 7.27), which worked well for helium, is inadequate to confirm the existence of a bound state for the negative hydrogen ion.

ψ1(r1,r2)≡z3πa3e-z(r1+r2)/a (7.27)

Chandrasekhar used a trial wave function of the form

ψ(r1,r2)≡A[ψ1(r1)ψ2(r2)+ψ2(r1)ψ1(r2)] (7.62).

Where

ψ1(r)≡z13Ï€²¹3e-z1/a,ψ1(r)≡z23Ï€a3e-z2/a,(7.63)

In effect, he allowed two different shielding factors, suggesting that one electron is relatively close to the nucleus, and the other is farther out. (Because electrons are identical particles, the spatial wave function must be symmetries with respect to interchange. The spin state—which is irrelevant to the calculation—is evidently anti symmetric.) Show that by astute choice of the adjustable parameters Z1and Z2you can get<H>less than -13.6.

Answer: :<H>=E1x6+y6(-x8+2x7+12x6y2-12x5y2-18x3y4+118xy6-12y8).

Wherex≡z1+z2.y≡2z1z2Chandrasekhar usedZ1=1.039

(Since this is larger than 1, the motivating interpretation as an effective nuclear charge cannot be sustained, but never mindit’s still an acceptable trial wave function) andZ2=0.283.

Short Answer

Expert verified

<H>=E1x6+y6(-x8+2x7+12x6y2-12x5y2-18x3y4+118xy6-12y8).

Step by step solution

01

Concept.

Total wave function is:

ψ(r→1,r→2)=A[ψ1(r1)ψ2(r2)ψ2(r1)ψ1(r2)],

Where, ψ1(r)≡z13Ï€²¹3e-z1/a,ψ2(r)≡z23Ï€²¹3e-z2r/a

The total wave function when residual forces are neglected is simply the product of the one quasi-particle proton and the one-quasi-particle neutron wave function.

02

Showing that the astute choice of adjustable parameters

Here the astute choice of adjustable parameters

A2(1+2S2+1)=1

where

S=∫ψ1(r)ψ2(r)d3r=Z1Z23Ï€²¹3∫e-(z1+z2)r/a4Ï€°ù2dr=4a3y232a3(z1+z2)3=yx3.A2=121+(y/x)6.

H=-h22m∇12+∇22-e24πo˙01r1+1r2+e24πo˙01r1+r2,Hψ=A-h22m∇12+∇22-e24πo˙0z1r1+z2r2ψ1(r1)ψ2(r2)+-h22m∇12+∇22-e24πo˙0z1r1+z2r2ψ2(r1)ψ1(r2)+Ae24πo˙0Z1-1r1+Z2-1r2ψ1(r1)ψ2(r2)+Z1-1r1+Z2-1r2ψ2(r1)ψ1(r2)+Veeψ.

whereVee≡e24πo˙01r1-r2.

03

The first term:

The term in first curly brackets is z12+z22E1ψ1(r1)ψ2(r2)+z22+z12E1ψ2(r1)ψ1(r2), so

Hψ=z12+z22+Ae24πo˙0z1-1r1+z2-1r2ψ1(r1)ψ2(r2)+z1-1r1+z2-1r2ψ2(r1)ψ1(r2),+Veeψ.H=z12+z22E1+Vee+A2e24πo˙0×ψ1(r1)ψ2(r2)+ψ2(r1)ψ1(r2)|z1-1r1+z2-1r2|ψ1(r1)ψ2(r2)>+z1-1r1+z2-1r2ψ2(r1)ψ1(r2)

Last term is equal to:

localid="1658396571051" =(z1-1)ψ1(r1)1r1ψ1(r1)+(z2-1)ψ2(r2)1r2ψ2(r2)+(z2-1)ψ1(r1)1r1ψ1(r1)ψ2(r2)|ψ1(r2)+(z1-1)ψ1(r1)|ψ2(r1)ψ2(r2)1r2ψ1(r2)+(z2-1)ψ2(r1)1r1ψ1(r1)ψ1(r2)|ψ2(r2)+(z2-1)ψ2(r1)|ψ1(r1)ψ1(r2)1r2ψ2(r2)++(z2-1)ψ2(r1)1r1ψ2(r1)+(z2-1)ψ1(r2)1r2ψ1(r2)

Instead r1,r2,we put r and last expression will simplify:

localid="1658404037903" =2z1-11r1+2z1-11r2+2z1-1ψ1|ψ2ψ11rψ2+2z1-1ψ1|ψ2ψ11rψ2But1r1=ψ1(r)1rψ1(r)=Zla;1r2=Z2a,SoH=Z12+Z22E1+A2e24πo˙021a(Z1-1)Z1+1a(Z2-1)Z2+(Z1+Z2-2)ψ1|ψ2ψ11rψ2+Vee.Andψ1|ψ2=S=(y/x)3,so

ψ11rψ2=Z1Z23Ï€²¹34π∫e-(z1+z2)r/ardr=y32a3aZ1+Z22=y32ax2.H=x2-12y2E1+A2e24Ï€oË™02az12+z22-(Z1+Z2)(x-2)yx3y32x2+Vee=x2-12y2E1+4E1A2x2-12y2-x+12(x-2)y6x5+Vee.

04

The calculation of the values:

Now we need to calculate

Vee=e24πo˙0ψ1r1-r2ψ.=e24πo˙0A2ψ1(r1)ψ2(r2)+ψ2(r1)ψ1(r2)1r1-r2ψ1(r1)ψ2(r2)+ψ2(r1)ψ1(r2)=e24πo˙0A22ψ1(r1)ψ2(r2)1r1-r2ψ1(r1)ψ2(r2)+2ψ1(r1)ψ2(r2)1r1-r2ψ2(r1)ψ1(r2)

First and last term are the same, so we have :

=2e24πo˙0A2(B+C),

localid="1658488118584" whereψ1(r1)ψ2(r2)1r1-r2ψ1(r1)ψ2(r2);C=ψ1(r1)ψ2(r2)1r1-r2ψ2(r1)ψ1(r2).B=Z13Z23Ï€²¹32∫e-2Z1r1/ae-2Z2r2/a1r1-r2d3r1d3r2

ther2integralis∫e-2Z2r2/a1r12+r22-2r1r2cosθ2d3r2

=Ï€²¹3Z23r11-1+Z2r1ae-2Z2r1/aEq.8.25,butwitha→2Z2aB=Z13Z23Ï€²¹32(Ï€²¹3)Z234π∫0∞e-2Z1r1/a1r11-1+Z2r1ae-2(Z1+Z2)r1/ar12dr1

=4Z13a3∫0∞r1e-2Z1r1/a-r1e-2(Z1+Z2)r1/a-Z2ar12e-2(Z1+Z2)r1/adr1.=4Z13a3a2Z12-a2(Z1+Z22-Z2a2a2(Z1+Z23=Z13a1Z12-1(Z1+Z2)2-Z2(Z1+Z2)3

=Z1Z2a(Z1+Z2)1+Z1Z2(Z1+Z2)=y24ax1+y24x2C=Z1Z2Ï€²¹32∫e-Z1r1/ae-Z2r2/ae-Z2r1/ae-Z1r2/a1r1-r2d3r1d3r2=Z1Z23Ï€²¹32∫e-z1+z2r1+r2/a1r1-r2d3r1d3r2

The integral is the same as, only with a →4Z1+Z2a.

Comparing Eqs. 7.20 and 7.25, we see that the integral itself was

Vee=e24Ï€oË™08Ï€²¹32∫e-4(r1+r2)/ar1-r2d3r1d3r2Vee=54ae24Ï€oË™0=-52E1=34eV54aÏ€²¹382=5256Ï€2a5.soC=Z1Z23Ï€²¹325Ï€225645a5Z1+Z25=20aZ1Z23Z1+Z25=516ay6x5

Vee=2Vee=e24πo˙0A2y24ax1+y24x2+516y6x5=2A2(-2E1)y24x1+y24x2+5y44x4.H=E1x2-12y2-21+(y/x)6x2-12y2-x+12(x-2)y6x5-21+(y/x)6y24x1+y24x2+5y44x4.

localid="1658492734438" =E1x6+y6x2-12y2x6+y6-2x2x2-12y2-x+12y6x4-y6x5+y24x+y416x3+5y616x5.=E1x6+y6x8+x2y6-12x6y2-12y8-2x8+x6y2+2x7-x2y6+2xy6--12x5y2-18x3y4-58xy6.=E1x6+y6-x8+2x7+12x6y2-12x5y2-18x3y4+118xy6-12y8.

Mathematical finds the minimum

ofHatx=1.32245,y=1.08505,correspondingtoZ1=1.0392,Z2=0.2832,.

At this point,Hmin=1.0266E,=-13.962eV, which is less than -13.6 eV but not by much.

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with 91Ó°ÊÓ!

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

Suppose we used a minus sign in our trial wave function (Equation 7.37):ψ=A[ψ0(r1)-ψ0(r2)]): Without doing any new integrals, find F(x) (the analog to Equation ) for this case, and construct the graph. Show that there is no evidence of bonding. (Since the variational principle only gives an upper bound, this doesn't prove that bonding cannot occur for such a state, but it certainly doesn't look promising). Comment: Actually, any function of the form ψ=A[ψ0(r1)+eiϕψ0(r2)]has the desired property that the electron is equally likely to be associated with either proton. However, since the Hamiltonian (Equation 7.35 ) is invariant under the interchange P:r1↔r2 , its eigen functions can be chosen to be simultaneously eigen functions of P . The plus sign (Equation 7.37) goes with the eigenvalue +1 , and the minus sign (Equation 7.52 ) with the eigenvalue -1 ; nothing is to be gained by considering the ostensibly more general case (Equation 7.53), though you're welcome to try it, if you're interested.

Use a gaussian trial function (Equation 7.2) to obtain the lowest upper bound you can on the ground state energy of (a) the linear potential V(x)=α|x| (b) the quartic potential:V(x)=αx4

As an explicit example of the method developed inProblem 7.15, consider an electron at rest in a uniform magnetic field B=B2kÁåœfor which the Hamiltonian is (Equation 4.158):

H=-γµþ (4.158).

H0=eBzmSz (7.57).

The eigenspinors, xaarelocalid="1655969802629" xb,andthecorrespondingenergies,EaandEb,aregiven in Equation 4.161. Now we turn on a perturbation, in the form of a uniform field in the x direction:

{x+,withenergyE+=-γB0ħ/2x-,withenergyE-=+-γB0ħ/2 (4.161).

H'=eBxmSx (7.58).

(a) Find the matrix elements of H′, and confirm that they have the structure of Equation 7.55. What is h?

(b) Using your result inProblem 7.15(b), find the new ground state energy, in second-order perturbation theory.

(c) Using your result inProblem 7.15(c), find the variation principle bound on the ground state energy.

a) Use the variational principle to prove that first-order non-degenerate perturbation theory always overestimates (or at any rate never underestimates) the ground state energy.

(b) In view of (a), you would expect that the second-order correction to the ground state is always negative. Confirm that this is indeed the case, by examining Equation 6.15.

(a) Prove the following corollary to the variational principle: If ⟨ψ∣ψgs⟩=0thená±áñ≥Efe whereEfe is the energy of the first excited state. Comment: If we can find a trial function that is orthogonal to the exact ground state, we can get an upper bound on the first excited state. In general, it's difficult to be sure that is orthogonal toψgsi since (presumably) we don't know the latter. However, if the potentialV(x) is an even function of x, then the ground state is likewise even, and hence any odd trial function will automatically meet the condition for the corollary.

(b) Find the best bound on the first excited state of the one-dimensional harmonic oscillator using the trial functionψ(x)=Axe-bx2

See all solutions

Recommended explanations on Physics Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.