A particle moves non-relativistically
in a three-dimensional harmonic oscillator potential. The potential energy
function is U(r) = αr^{2}. Use
spherical coordinates.

(a) Use the trial function
ψ = A exp(-br^{2}) in the
variational method to find the ground-state energy and the normalized wave
function.

(b) Comment on the quality of the following trial functions. Do not
do the calculations, just comment on the functions' good and bad points.

ψ = A exp(-br)

ψ = A sech(-br) [sech(x) = 1/cosh(x) =
2/(e^{x} + e^{-x})]

ψ = A (1 + b^{2}r^{2})^{-1
}ψ = A sin(br)/(br)

Note: [∫_{0}^{∞} r^{4}exp(-r^{2})dr = 3π^{½}/8.
∫_{0}^{∞}
r^{2}exp(-r^{2})dr =
π^{½}/4.]

Solution:

- Concepts:

The variational method - Reasoning:

The variational method often yields a very good estimate for the ground state energy of a system. - Details of the calculation:

(a) The ground state has zero angular momentum. Therefore we have

ψ(r) = [u(r)/r]Y_{00}(θ,φ).

((-ħ^{2}/2m)(∂^{2}/∂r^{2}) + ar^{2}))u_{k0}(r) = E_{k0}u_{k0}(r).

Choose u_{k0}(r) = B r exp(-br^{2}). Here b is an adjustable parameter.

u(r) must be zero at r = 0 and at r = infinity and should not have any nodes in the region 0 < r < infinity.

B^{2}∫_{0}^{∞}r^{2}exp(-2br^{2})dr = 1, B^{2}= 8*2^{½}b^{3/2}/π^{½}.

Then

<H> = B^{2}∫_{0}^{∞}r exp(-br^{2})((-ħ^{2}/2 m)(∂^{2}/∂r^{2}) + αr^{2}) r exp(-br^{2})dr

= B^{2}(-ħ^{2}/2m)∫_{0}^{∞}[-6br^{2}+ 4r^{4}b^{2}]exp(-2br^{2})dr + B^{2}α∫_{0}^{∞}r^{4}exp(-2br^{2})dr

= B^{2}[(3ħ^{2}/m)(b^{-½}2^{-3/2})∫_{0}^{∞}r^{2}exp(-r^{2})dr - B^{2}[(ħ^{2}/2m)(b^{-½}2^{-½})∫_{0}^{∞}r^{4}exp(-r^{2})dr

+ B^{2}α(b^{-5/2}2^{-5/2})∫_{0}^{∞}r^{4}exp(-r^{2})dr

= (ħ^{2}/2m)3b + 3α/(4b).

[∫_{0}^{∞}r^{4}exp(-r^{2})dr = 3π^{½}/8. ∫_{0}^{∞}r^{2}exp(-r^{2})dr = π^{½}/4.]

d<H>/db = 0, (3ħ^{2}/2m) = 3α/(4b^{2}) , b^{2}= αm/(2ħ^{2}).

<H> = 3ħα^{½}/(2m)^{½}. E(upper limit) = 3ħα^{½}/(2m)^{½}.

Setting α = ½mω^{2}, we see that E(upper limit) = (3/2)(ħω), the exact result.

If ψ = A exp(-br^{2}), then the normalization constant A = B/(4π)^{½}.

(b) The wave function must satisfy the boundary conditions and be normalizable. rψ (r) must be zero at r = 0 and at r = infinity. In the region 0 < r < infinity rψ and its derivative with respect to r must be continuous. In the region 0 < r < infinity rψ for the ground state has no nodes.

ψ = A exp(-br) is an acceptable trial function

ψ = A sech(-br) is an acceptable trial function

ψ = A (1 + b^{2}r^{2})^{-1}is an acceptable trial function.^{ }ψ = A sin(br)/(br) is not an acceptable trial function. It has nodes and is not normalizable.

Assume that you are asked to use the variational method to estimate the
ground state energy of a particle with mass m, given a one-dimensional potential
energy function U(x). Assume U(x) = -Bx, x < 0, U(x) = Fx, x > 0,
where B and F are positive constants with the appropriate units.

(a) Assume B ≠ F. Choose a trial wave function ψ(x) with two
adjustable parameters and explain why it is an acceptable wave function.
(You do not have to find the ground-state enery.)

(b) Assume B = F. Choose a trial wave function ψ(x) with one
adjustable parameter and explain why it is an acceptable wave function.
(You do not have to find the ground-state enery.)

(c) Let B --> ∞. Choose a trial wave function ψ(x)
with one adjustable parameter and estimate the ground state energy

Solution:

- Concepts:

The variational method - Reasoning:

The trial wave functions ψ(x) must satisfy the boundary conditions at infinity, ψ(±∞) = 0, or ψ(x) = 0 in regions where U(x) = ∞, the derivative dψ/dx must be continuous at all x except where U(x) --> ∞, and the ground state wave function has no nodes except at x = ±∞ or where U(x) --> ∞. Any function that satisfies these conditions can be written as a linear combination of eigenfunctions of H and is an acceptable trial function giving <H> > E_{0}. - Details of the calculation:

(a) The trial wave function must be asymmetric. For example, we can use the a linear combination of exp(-c^{2}x^{2}) and x*exp(-c^{2}x^{2}).

ψ(x) = N[exp(-cx^{2}) + b*x*exp(-cx^{2})]. This wave function satisfies the conditions stated above. Here b and c are the adjustable parameters. We minimize <H> with respect to b and c to find an upper limit for the ground state energy.

(b) The trial wave function must be symmetric. For example, ψ(x) = N exp(-c^{2}x^{2}). This wave function satisfies the conditions stated above. Here c is the trial parameter.

(c) Choose the trial function ψ(x) = N x exp(-cx) to make the integrals easy to evaluate. This wave function satisfies the conditions stated above. Here c is the adjustable parameter.

Normalize the wave function: N^{2}∫_{0}^{∞}x^{2}exp(-2cx) dx = N^{2}/(4c^{3}) = 1, N^{2}= 4c^{3}.

<H> = N^{2}∫_{0}^{∞}x exp(-cx)[(-ħ^{2}/2m)(∂^{2}/∂x^{2}) + Fx]x exp(-cx) dx

= N^{2}(-ħ^{2}/2m)∫_{0}^{∞}(-2cx + c^{2}x^{2}) exp(-2cx) dx + N^{2}F∫_{0}^{∞}x^{3}exp(-2cx) dx

= N^{2}(ħ^{2}/4mc - ħ^{2}/8mc + 6F/(2c)^{4}).

d<H>/dc = 0 --> c^{3}= 3Fm/(2ħ^{2}).

Insert c back into the equation for <H>.

E = 1.96 F^{⅔}ħ^{⅔}/m^{⅓}.

Assume the potential for the deuterium is given by

U(r) = -U_{0}exp(-r^{2}/a^{2}), with U_{0} = 37 MeV and a = 2.2*10^{-15 }m.

Estimate the ground state energy of the deuterium.

Solution:

- Concepts:

The variational method - Reasoning:

The variational method often yields a very good estimate for the ground state energy of a system. -
Details of the calculation:

The relative motion of the proton and neutron is described in the same way as the motion of a fictitious particle of reduced mass m in a central potential. The ground state has zero angular momentum. Therefore we have

Φ_{k00}(r) = [u_{k0}(r)/r]Y_{00}(θ,φ).

((-ħ^{2}/2μ)(∂^{2}/∂r^{2}) - U_{0}exp(-r^{2}/a^{2}))u_{k0}(r) = E_{k0}u_{k0}(r).

Choose u_{k0}(r) = A r exp(-r^{2}/b^{2}). Here b is an adjustable parameter.

u_{k0}(r) must be zero at r = 0.

A^{2}∫_{0}^{∞}r^{2}exp(-2r^{2}/b^{2})dr = 1, A^{2}= 8*2^{½}/[π^{½}b^{3}].

Then

<H> = A^{2}∫_{0}^{∞}r exp(-r^{2}/b^{2})((-ħ^{2}/2μ)(∂^{2}/∂r^{2}) - U_{0}exp(-r^{2}/a^{2})) r exp(-r^{2}/b^{2})dr

= A^{2}(-ħ^{2}/2μ)∫_{0}^{∞}[-6r^{2}/b^{2}+ 4r^{4}/b^{4}]exp(-2r^{2}/b^{2})dr - A^{2}U_{0}∫_{0}^{∞}r^{2}exp(-r^{2}(2/b^{2}+ 1/a^{2}))dr

= A^{2}[(3ħ^{2}/μ)(b/2^{3/2})∫_{0}^{∞}r^{2}exp(-r^{2})dr - A^{2}[(ħ^{2}/2μ)(b/2^{½})∫_{0}^{∞}r^{4}exp(-r^{2})dr

+ A^{2}U_{0}[a^{2}b^{2}/(2a^{2}-b^{2})]^{3/2}∫_{0}^{∞}r^{2}exp(-r^{2})dr

= A^{2}[(3ħ^{2}/μ)(b/2^{3/2})(π^{½}/4) - A^{2}[(ħ^{2}/2μ)(b/2^{½})(3π^{½}/8)

+ A^{2}U_{0}[a^{2}b^{2}/(2a^{2}+b^{2})]^{3/2}(π^{½}/4).

[∫_{0}^{∞}r^{4}exp(-r^{2})dr = 3π^{½}/8. ∫_{0}^{∞}r^{2}exp(-r^{2})dr = π^{½}/4.]

<H> = 3ħ^{2}/(2μb^{2}) - U_{0}(2a^{2}/(2a^{2}+b^{2}))^{3/2}.

Let x = 2a^{2}/b^{2}.

<H> = 3ħ^{2}x/(4μa^{2}) - U_{0}/(1 + 1/x)^{3/2}.

d<H>/dx = 0 --> 3ħ^{2}/(4μa^{2}) - (3U_{0}/2x^{2})(1 + 1/x)^{-5/2}= 0.

ħ = 6.58*10^{-22}MeV-s, μ = 938/2 MeV/c^{2}_{, }a = 2.2*10^{-15 }m

ħ^{2}/(μa^{2}) = 17.2 MeV,

0.348 = (3/2)(1/x^{2})(1 + 1/x)^{-5/2}. 0.232 = (1/x^{2})(1 + 1/x)^{-5/2},

x = 0.64.

b^{2}= 3.125 a^{2}.

<H>/U_{0}= 0.223 - 0.244 = -0.0213

E(upper limit) = -0.786 MeV

The system is barely bound.

Consider the half space linear potential in one dimension

U(x) = ∞, x < 0, U(x) = Fx, x > 0.

Use the variational method to estimate the ground state energy. Choose
your own variational wave function.

Solution:

- Concepts:

The variational method - Reasoning:

We are asked to use the variational method to estimate the ground state energy. -
Details of the calculation:

H = p^{2}/(2m) + Fx, x > 0.

The eigenfunctions of H must be zero at x = 0 and x = ∞. The derivative dψ/dx must be continuous at all x except x = 0. Any function that satisfies these boundary conditions can be written as a linear combination of eigenfunctions of H and is an acceptable trial function giving <H> > E_{0}.

We may choose as the trial function for x > 0 the eigenfunction of the first excited state of the harmonic oscillator. It satisfies the boundary conditions. We choose

ψ(x) = (√a/π^{¼}) 2ax exp(-½a^{2}x^{2}), with a^{2}= mω/ħ, for x >0,

ψ(x) = 0 for x < 0.

Here a is the adjustable parameter.

Note the change in the normalization constant, since we are only looking at x > 0.

<H> = <p^{2}/(2m)> + F<x>, x > 0.

<p^{2}/(2m)> = ½<T>_{harm. osc}= ¼<E>_{harm. osc}= (3/8)ħω = (3/8)ħ^{2}a^{2}/m,

since <T> = <U>.

<x> = (a/√π)∫_{0}^{∞}4a^{2}x^{3}exp(-a^{2}x^{2}) dx = (4/(a√π))∫_{0}^{∞}x^{3}exp(-x^{2}) dx

= (2/(a√π))∫_{0}^{∞}y exp(-y) dy = 2/(a√π).

<H> = (3/8)ħ^{2}a^{2}/m + 2F/(a√π).

d<H>/da = 0 --> a^{3}= (8/3)Fm/(ħ^{2}√π).

<H> = (1/a)3F/√π = (F^{⅔}ħ^{⅔}/m^{⅓})[(π^{⅓}/√π)(3/8)^{⅓}3]

= 1.86 (F^{⅔}ħ^{⅔}/m^{⅓}).

Exact solutions:

E = 1.85 F^{⅔}ħ^{⅔}/m^{⅓}.

Let us choose a different trial function.

ψ(x) = N x exp(-cx), N^{2}∫_{0}^{∞}x^{2}exp(-2cx) dx = N^{2}/(4c^{3}) = 1, N^{2}= 4c^{3}.

Here c is the adjustable parameter.

<H> = N^{2}∫_{0}^{∞}x exp(-cx)[(-ħ^{2}/2m)(∂^{2}/∂x^{2}) + Fx]x exp(-cx) dx

= N^{2}(-ħ^{2}/2m)∫_{0}^{∞}(-2cx + c^{2}x^{2}) exp(-2cx) dx + N^{2}F∫_{0}^{∞}x^{3}exp(-2cx) dx

= N^{2}(ħ^{2}/4mc - ħ^{2}/8mc + 6F/(2c)^{4}).

d<H>/dc = 0 --> c^{3}= 3Fm/(2ħ^{2}).

Insert c back into the equation for <H>.

E = 1.96 F^{⅔}ħ^{⅔}/m^{⅓}.

A particle is bound to a site at the origin by a
linear “string” potential, U(r) = br, where b is the string tension.
In 3D the Hamiltonian is H = a∇^{2 }+ br, where
a = -ħ^{2}/2m.

(a) The energy eigenvalues are all proportional to a^{q}b^{p}.
Use dimensional analysis to determine both p and q.

(b) Assume a trial wave function of the form ψ(r) = Bexp(-d^{2}r^{2}/2).

(1) Determine the normalization
constant B.

(2) Calculate the expected energy and find the
optimum variational parameter d.

Solution:

- Concepts:

Dimensional analysis, the variational method - Reasoning:

We are instructed to apply those methods -
Details of the calculation:

(a) The eigenvalue equation for H is Hψ(r) = (a∇^{2 }+ br)ψ(r) = Ca^{q}b^{p}ψ(r).

Here C is a dimensionless constant.

The units of a^{q}b^{p}are the units of energy, in SI units a^{q}b^{p}has units of J.

a = -ħ^{2}/2m has units of Jm^{2}and b must have units of J/m.

For the units we therefore have (Jm^{2})^{q}(J/m)^{p}=J^{q+p}m^{2q-p}= J.

q + p = 1, 2q - p = 0, q = ⅓, p = ⅔.

(b) Let ψ(r) = Bexp(-d^{2}r^{2}/2).

Normalization: 4πB^{2}∫_{0}^{∞}r^{2}exp(-d^{2}r^{2}/)dr = 1, B^{2}= (d/π^{½})^{3}.

Using the given wave function we find <H> = <ψ|H|ψ>.

<H> = 4πaB^{2}∫_{0}^{∞}exp(-d^{2}r^{2}/2)(1/r^{2})(d/dr)[(r^{2}(d/dr)exp(-d^{2}r^{2}/2)]r^{2}dr

+ 4πbB^{2}∫_{0}^{∞}exp(-d^{2}r^{2})r^{3}dr

= (4πaB^{2}/d)∫_{0}^{∞}(x^{4}-3x^{2})exp(-x^{2})dx + (4πbB^{2}/d^{4})∫_{0}^{∞}x^{3}exp(-x^{2})dx

= (4πaB^{2}/d)[3π^{½}/8 - 3π^{½}/4] + (2πbB^{2}/d^{4})∫_{0}^{∞}yexp(-y)dy

= (4πaB^{2}/d)[3π^{½}/8 - 3π^{½}/4] + (2πbB^{2}/d^{4})

= -(4d^{2})[3a/8 - (b/(2d^{3}π^{½})].

We now minimize <H> with respect to d.

d<H>/d(d) = 0. -8d[3a/8 - (b/(2d^{3}π^{½})] - (4d^{2})(3b/(2d^{4}π^{½}) = 0.

d^{3}= -2b/(3aπ^{½}) = 2b/(3|a|π^{½}).

d = [2b/(3|a|π^{½})]^{⅓}is the optimum variational parameter.

<H> = ([2b/(3|a|π^{½})]^{⅔})[9|a|/2] is an upper limit for the ground state energy.

It
is known that the stable H^{-} exists (two electrons bound to a
proton). Estimate the
ground state energy of H^{-} using the variational method.
Assume that each electron moves in a 1s orbit.
Neglect spin.

Useful
information:

ψ_{1s}(**r**) = (1/πa_{0}^{3})^{½}exp(-r/a_{0})

< ψ_{1s} |(1/r)| ψ_{1s} > = 1/a_{0}

< ψ_{1s} |∇^{2}| ψ_{1s} > = -1/a_{0}^{2}

∫∫d^{3}r d^{3}r'|ψ_{1s}(**r**)|^{2}|ψ_{1s}(**r'**)|^{2}(1/|**r **-** r**'|) = 5/(8a_{0})

Solution:

- Concepts:

The variational method, the H^{-}negative ion - Reasoning:

Neglecting spin, the Hamiltonian for the H^{-}negative ion is

H = P_{1}^{2}/(2m) + P_{1}^{2}/(2m) - e^{2}/r_{1}- e^{2}/r_{2}+ e^{2}/|**r**_{1 }-**r**_{2}|.

e^{2}= q_{e}^{2}/(4πε_{0}) in SI units.

We are asked to find the ground state energy using the variational method. - Details of the calculation:

Let us use as a trial function

ψ(**r**_{1},**r**_{2}) = (1/πa^{3})exp(-(r_{1}+ r_{2})/a)

(We are guided by the ground state wave function of hydrogen )

ψ(**r**) = (1/πa_{0}^{3})^{½}exp(-r/a_{0})

The adjustable parameter in our trial function is a. The wave function is normalized.

ψ(**r**_{1},**r**_{2}) is a product function, ψ(**r**_{1},**r**_{2}) = Φ_{1}(**r**_{1})Φ_{2}(**r**_{2}).

Φ_{1}(**r**_{1}) = (1/πa^{3})^{½}exp(-r_{1}/a), Φ_{2}(**r**_{2}) = (1/πa^{3})^{½}exp(-r_{2}/a).

<ψ|p^{2}_{1}/(2m)|ψ> = <Φ_{1}|p^{2}_{1}/(2m)|Φ_{1}> = -ħ^{2}/(2m)< Φ_{1}|∇_{(1)}^{2}|Φ_{1}> = ħ^{2}/(2ma^{2}).

<ψ|p^{2}_{2}/(2m)|ψ> = <Φ_{2}|p^{2}_{2}/(2m)|Φ_{2}> = -ħ^{2}/(2m)< Φ_{2}|∇_{(2)}^{2}|Φ_{2}> = ħ^{2}/(2ma^{2}).

<ψ|e^{2}/r_{1}|ψ> = e^{2}<Φ_{1}|1/r_{1}|Φ_{1}> = e^{2}/a.

<ψ|e^{2}/r_{2}|ψ> = e^{2}<Φ_{2}|1/r_{2}|Φ_{2}> = e^{2}/a.

<ψ|e^{2}/|**r**_{1 }-**r**_{2}||ψ> = ∫d^{3}r_{1}d^{3}r_{2}|Φ_{1}(**r**_{1})|^{2}|Φ_{2}(**r**_{2})|^{2}(1/|**r**_{1 }-**r**_{2}|) = 5e^{2}/(8a).

Therefore <H> = ħ^{2}/(ma^{2}) - 2e^{2}/a + 5e^{2}/(8a) = ħ^{2}/(ma^{2}) - 11e^{2}/(8a).

d<H>/da = -2ħ^{2}/(ma^{3}) + 11e^{2}/(8a^{2}) = 0, a = 16ħ^{2}/(11me^{2})

E_{ground state}(upper bound) = -(12½56)me^{4}/ħ^{2}= -0.945 * 13.6 eV

This calculation would indicate that the H^{-}ion is not stable, since the ground state energy is less negative than that of a free electron and a neutral hydrogen atom.

(H^{-}is a stable negative ion.)