3.E: The Schrödinger Equation and a Particle in a Box (Exercises)

These are homework exercises to accompany Chapter 3 of McQuarrie and Simon's "Physical Chemistry: A Molecular Approach" Textmap.

Q3.2

Determine from the following operators which are linear and nonlinear:

1. $$\hat{A}f(x)= f(x)^2$$ [square f(x)]
2. $$\hat{A}f(x)= f^*(x)$$ [form the complex conjugate of f(x)]
3. $$\hat{A}f(x)= 0$$ [multiply f(x) by zero]
4. $$\hat{A}f(x)= [f(x)]^{-1}$$ [take the reciprocal of f(x)]
5. $$\hat{A}f(x)= f(0)$$ [evaluate f(x) at x=0]
6. $$\hat{A}f(x)= \ln f(x)$$ [take the log of f(x)]

S3.2

It is important to note that an operator $$\hat{A}$$ is linear if

$\underbrace{\hat{A}[c_1f(x)+c_2f_2(x)]}_{\text{left side}}= \underbrace{c_1\hat{A}f_1(x)+c_2\hat{A}f_2(x) }_{\text{right side}}$

and the operator is nonlinear if

$\underbrace{ \hat{A}[c_1f_1(x)+c_2f_2(x)]}_{\text{left side}} \neq \underbrace{ c_1\hat{A}f_1(x)+c_2\hat{A}f_2(x) }_{\text{right side}}$

a)

Evaluate the left side

\begin{align} \hat{A}[c_1f(x)+c_2f_2(x)] &= [c_1f_1(x)+c_2f_2(x)]^2 \\ &= c_1^2 f_1(x)^2+2c_1f_1(x) c_2f_2(x)+c_2^2f_2(x)^2 \end{align}

Evaluate the right side

$c_1 \hat{A} f_1(x)+c_2\hat{A}f_2(x)=c_1[f_1(x)]^2+c_2[f_2(x)]^2 \neq \hat{A}[c_1f_1(x)+c_2f_2(x)]$

This operator is nonlinear

b)

Evaluate the left side

$\hat{A}[c_1f_1(x)+c_2f_2(x)] = c_1^*f_1^*(x) + c_2^*f_2(x)$

Evaluate the right side

$c_1\hat{A}f_1(x) + c_2\hat{A}f_2(x) = c_1f_1^*(x) + c_2f_2^*(x)$

$= \hat{A}[c_1f_1(x) + c_2f_2(x)]$

This operator is linear

c)

Evaluate the left side

$\hat{A}[c_1f_1(x)+c_2f_2(x)] = 0$

Evaluate the right side

$c_1\hat{A}f_1(x) + c_2\hat{A}f_2(x) = c_1f_1(x) + c_2f_2(x) = 0$

$= \hat{A}[c_1f_1(x) + c_2f_2(x)]$

This operator is linear

d)

Evaluate the left side

$\hat{A}[c_1f_1(x)+c_2f_2(x)] = \dfrac{1}{c_1f_1(x) + c_2f_2(x)}$

Evaluate the right side

$c_1\hat{A}f_1(x) + c_2\hat{A}f_2(x) = \dfrac{c_1}{f_1(x)} + \dfrac{c_2}{f_2(x)}$

$\neq \hat{A}[c_1f_1(x) + c_2f_2(x)]$

This operator is nonlinear

e)

Evaluate the left side

$\hat{A}[c_1f_1(x)+c_2f_2(x)] = c_1f_1(0) + c_2f_2(0)$

Evaluate the right side

$= c_1\hat{A}f_1(x) + c_2\hat{A}f_2(x)$

This operator is linear

f)

Evaluate the left side

$\hat{A}[c_1f_1(x)+c_2f_2(x)] = \ln [c_1f_1(x) + c_2f_2(x)]$

Evaluate the right side

$c_1\hat{A}f_1(x) + c_2\hat{A}f_2(x) = c1 \ln f_1(x) + c_2 \ln f_2(x)$

$\neq \hat{A}[c_1f_1(x) + c_2f_2(x)]$

This operator is nonlinear

Q3.8

Show that for a particle in a box with length a with state $$n=3$$ that there are 3 locations along the x axis where the probability density is at a maximum.

S3.8

The probability density for a particle in a box for state $$n=3$$ is

${\psi}^*{\psi}=\dfrac{2}{a}{\sin}^2\dfrac{3px}{a}$

To maximize the probability density, take its derivative and set it equal to zero and solve for $$x$$ .

$\dfrac{d}{dx}\left[\dfrac{2}{a}{{\sin}^2 \dfrac{3{\pi}x}{a}\ }\right]=\dfrac{2}{a}\cdot 2\cdot {\sin \dfrac{3{\pi}x}{a}\ }\cdot {\cos \dfrac{3{\pi}x}{a}\ }\cdot \dfrac{3\pi}{a}=0$

${\sin \dfrac{3{\pi}x}{a}\ }{\cos \dfrac{3{\pi}x}{a}\ }=0$

We want to not choose values of x that make $${\sin \dfrac{3{\pi}x}{a}\ }=0$$ , as that means that the probability density will be zero. We will only choose the zeros of $${\cos \dfrac{3{\pi}x}{a}\ }$$. So the possible values for x which make

${\cos \dfrac{3{\pi}x}{a}\ } =0$

are

$\dfrac{3px}{a}=\dfrac{2m+1}{2}p\ \ \ \ \ \ m=0,1,2,\dots$

$x=\dfrac{\left(2m+1\right)a}{6}$

We only choose $$m=0,1,2$$ and not 3 because $$m=3$$ would give $$x=\dfrac{7a}{6}$$ which is outside the box. So the locations are

$x=\dfrac{a}{6}$

$x=\dfrac{a}{2}$

$x=\dfrac{5a}{6}$

Q3.13

What range for L is possible for $$\sigma_x$$ given:

$\sigma_x = \sqrt{\langle x^2\rangle - \langle x \rangle^2}$

where L is the length of the box? Hint: Remember that $$\sigma_x$$ is the uncertainty in the position of a particle in a box.

S3.13

For a particle in a box:

$$\langle x \rangle = \dfrac{\text{L}}{2}$$

and

$$\langle x^2 \rangle = \dfrac{L^2}{3}-\dfrac{L^2}{2n^2\pi^2}$$

$$\sigma_x = \sqrt{\dfrac{L^2}{3}-\dfrac{L^2}{2n^2\pi^2} - (\dfrac{\text{L}}{2})^2}$$

By inspection, only values of $$\sigma_x$$ less than $$L$$ will make this statement true.

Q3.14

Using the trigonometric identity

$\cos(2x)=2\cos^{2}x-1$

show that

$\int_0^a 2 \cos^2 \dfrac{n\pi x}{a} -1 dx = 0$

$\int_0^a \cos \dfrac{2n\pi x}{a} dx = 0$

$\dfrac{a}{2n\pi} \sin(2n\pi) = 0$

Q3.18

Is the wavefunction $$\phi_n= \sqrt{\dfrac{2}{L}} \sin{(\dfrac{n\pi x}{L})}$$ orthonormal over $$0 \leq{x} \leq{L}$$. Explain your reasoning.

S3.18

For a wavefunction to be orthonormal, it has to satisfy these conditions 1.) it has to be orthogonal and 2.) it has to be normalized.

To show that it is orthogonal:

$\int_0^L \phi_m \phi\ast_n dx = 0$

when $$m\neq$$ n

To show that that the wavefunction is normalized it must follow that

$\int_0^L \phi_n \phi\ast_n dx = 1$

when m=n

Because our wavefunction satisfies both conditions, it is an orthonormal function.

We find that

$\langle \psi_3 |\psi_3 \rangle = \int \limits _{-\infty}^{\infty} \dfrac{2}{L}(\sin\dfrac{3\pi x}{L})^2 dx = 1$

and

$\langle \psi_4 |\psi_4 \rangle = \int \limits _{-\infty}^{\infty} \dfrac{2}{L}(\sin\dfrac{3\pi x}{L})(\sin\dfrac{4\pi x}{L}) dx = 0$

From orthogonality, we can learn that if n is not equal to m, our dot product will always be zero. But if n = m our dot product will equal 1.

Q3.22

What is the Heisenberg Uncertainty Principle? Do position and momentum follow the uncertainty principle; why or why not? If they do, what is the minimum uncertainty in the velocity of an electron if it is known to be within 1.5nm of a nucleus?

S3.22

The Heisenberg Uncertainty Principle states that two properties that follow cannot be simultaneously measured to arbitrary precision. Position and momentum follow the principle. If one were to try and commute these two operators, one would not get zero and therefore the properties do not commute. If they do not commute then they cannot be measured to arbitrary precision.

We know that

$\Delta x \Delta p \geq \dfrac{\hbar}{2}$

And that p = mv. This gives

$m\Delta x \Delta v \geq \dfrac{\hbar}{2}$

The mass of an electron is known is $$m_e \approx 9.1 \times 10^{-31}\;kg$$. The problem also gives $$\Delta x$$ to be 1.5 nm. From here, it becomes a plug and chug to solve for $$\Delta v$$.

$\Delta v = 3.86 \times 10^4$

Q3.23

Describe the degeneracies of a two-dimensional box whose two sides have different lengths.

S3.23

The energies of a two-dimensional box is given by,

$E = \dfrac{h^2}{8m} \left(\dfrac{n^{2}_{x}}{a^2}+\dfrac{n^{2}_{y}}{b^2}\right)$

We can see that even if a $$\ne$$ b the energy levels will not necessarily be degenerate.

Q3.26

How many degenerate states do the first three energy levels for a three-dimensional particle in a box have if a=b=c?

Q3.27

Metal porphyrin molecules are commonly in many proteins and it has the general structure. This molecule is planar, so we can approximate π electrons as being confined inside a square. What are energy levels and corresponding degeneracies of a particle in a square of side m? Porphyrin molecules have 18 π electrons. If the length of the molecule is 850 pm, what is the lowest energy absorption of the porphyrin molecules? (The experimental value ≈ 17 000 cm－1

S3.26

The first energy level is E(1,1,1) which has no degeneracy.

The second energy state is E(2,1,1)=E(1,2,1)=E(1,1,2), therefore it has three degenerate states.

The third energy state is E(2,2,1)=E(2,1,2)=E(1,2,2), therfeore it has three degenerate states.

Q3.26

For a two dimensional box of width $$w$$ and height $$h=\sqrt{a}w$$, calculate all possible energy combinations between $$E_{11}$$ and $$E_{33}$$ note any degeneracy.

S3.26

The energy of a two dimensional particle in the box has the form,

$E = \dfrac{h^2}{8m}\Bigg(\dfrac{n_x^2}{w^2}+\dfrac{n_y^2}{h^2}\Bigg)$

In this specific case $$h=\sqrt{a}w$$ so we can simplify the problem to,

$E = \dfrac{h^2}{8m}\Bigg(\dfrac{n_x^2}{w^2}+\dfrac{2n_y^2}{w^2}\Bigg)$

Now we can tabulate the energy level indicating degeneracy.

$$E_{xy}$$ Degeneracy $$\dfrac{E8mw^2}{h^2}$$
$$E_{11}$$ 1 3
$$E_{12}$$. $$E_{31}$$ 2 5
$$E_{21}$$ 1 4
$$E_{22}$$ 1 6
$$E_{13}$$. $$E_{32}$$ 2 7
$$E_{23}$$ 1 8
$$E_{33}$$ 1 9

Q3.32

In this problem, we will explore the quantum-mechanical problem of a free particle that is not restricted to a finite region. Remember quantized energies of a particle in a box is a direct result from the boundary conditions set by the confines of the box.

When the potential energy $$V(x)$$ is equal to zero and the Schrodinger equation become

$\dfrac{d^2\psi}{dx^2} + \dfrac{2mE}{\hbar^2}\psi(x) = 0$

The two solutions to this Schrodinger equation are

$\psi_1(x) = A_1e^{ikx}$

$\psi_2(x) = A_2e^{-ikx}$

where $k = \dfrac{(2mE)^(1/2)}{\hbar}$

Show that $$\psi_1(x)$$ and $$\psi_2(x)$$ are solution to the Schrodinger equation where the potential energy $$V(x)$$ is equal to zero

S3.32

In order to prove that $$\psi_1(x)$$ and $$\psi_2(x)$$ are solutions we need to mention a few values

$p = \hbar k → k = p/\hbar →k = \dfrac{(2mE)^{1/2}}{\hbar}$

Now we have

$\dfrac{d^2Ae^{+-ikx}}{dx^2} + \dfrac{2mE}{\hbar^2}Ae^{+-ikx}= 0$

$A(+-ik)^2e^{+-ikx} + \dfrac{2mE}{\hbar^2}Ae^{+-ikx} = 0 → -k^2 + \dfrac{2mE}{\hbar^2} = 0$

Cancel the like terms

Thus, $$k = \dfrac{(2mE)^{1/2}}{\hbar}$$, which equals the original $$k$$ value

Q3.32

Show that $$E$$ had to be a positive value, since when $$E$$ is negative the wave function become unbounded for large $$x$$ values

S3.32

If $$E$$ < 0 then k becomes imaginary, $$k = ik$$

$$\psi = Ae^{+-ikx} = Ae^{+-i(ik)x} Ae^{+-kx}$$

For $$\psi_1(x) = A_1e^{-kx}$$ this will blow up for x → - $$\infty$$

For $$\psi_2(x) = A_2e^{kx}$$ this will blow up for x → $$\infty$$

Q3.32

With $$\hat{P}$$ $$\psi_1(x)$$ and $$\hat{P}$$ $$\psi_2(x)$$ as eigenvalue equations, show that

$\hat{P}\psi_1(x) = -i\hbar \dfrac{d\psi_1}{dx} = \hbar k\psi_1$

and

$\hat{P}\psi_2(x) = -i\hbar \dfrac{d\psi_2}{dx} = -\hbar k\psi_2$

S3.32

$\hat{P}\psi_1(x) = -i\hbar \dfrac{d\psi_1}{dx} = -i\hbar \dfrac{d}{dx}A_1e^{+ikx} = -i\hbar A_1(ik)e^{ikx} = +\hbar kA_1e^{ikx} = + \hbar k\psi_1$

$\hat{P}\psi_2(x) = -i\hbar \dfrac{d\psi_2}{dx} = -i\hbar \dfrac{d}{dx}A_2e^{-ikx} = -i\hbar A_2(-ik)e^{ikx} = -\hbar kA_1e^{-ikx} = -\hbar k\psi_2$

Now we can show that

$E = \dfrac{p^2}{2m} = \dfrac{+-(\hbar)^2}{2m} = \dfrac{\hbar^2 k^2}{2m}$

Q3.32

Show that $$\psi_1^{*}\psi_1(x) = A_1^{*}A_1 = \left | A_1\right |^2$$ and that $$\psi_2^{*}\psi_2(x) = A_2^{*}A_2 = \left | A_2\right |^2$$

S3.32

$\psi_1^{*}\psi_1(x) = (A_1e^{ikx})^*A_1e^{ikx} = A_1^{*}A_1e^{-ikx}e^{ikx} = A_1^{*}A_1e^{-ikx+ikx} = A_1^{*}A_1e^{0} = A_1^{*}A_1$

$\psi_2^{*}\psi_2(x) = (A_2e^{-ikx})^*A_2e^{-ikx} = A_2^{*}A_2e^{ikx -ikx} = A_2^{*}A_2e^{0} = A_2^{*}A_2$

$$\psi$$ has equal probability to be everywhere when $$\Delta x = \infty$$ and $$\Delta p = 0$$

Q3.33

Assuming that a particle is characterized by a standing de Broglie wave, come up with an equation for the allowed energies of a particle in a one-dimensional box.

S3.33

$\lambda = \dfrac{h}{p}$

Because the waves are standing waves, an integral number of half wave-lengths will fit in the box or:

$a = \dfrac{n\lambda}{2}$

and

$a = \dfrac{nh}{2p}$

Solving for $$p$$ yields

$p = \dfrac{nh}{2a}$

and the corresponding energy is

$E = \dfrac {mv^2}{2} = \dfrac{p^2}{2m} =\dfrac{1}{2m} \dfrac {n^2h^2}{4a^2} = \dfrac{n^2h^2}{8ma^2}$

Q3.33

Derive the lowest allowed velocity for a proton in a box of length 10-14m (approximate size of nucleus), assuming the particle is described by a standing deBroglie wave.

S3.33

The deBroglie wavelength is

$λ = \dfrac{h}{p} = \dfrac{h}{m_pv}$

For a one dimensional wave that has nodes on both ends of a box, an integer number of half wavelengths can fit, so

$n \left(\dfrac{λ}{2} \right) = L$

Substituting this wavelength in the deBroglie relationship, one gets

v = hn/2mpL

lowest allowed velocity will have n = 1

v = (6.626*10-34)(1) / 2*(1.67*10-27)(10-14)

= 19.8 *106 m/s

Q3.33

If a particle in a one-dimensional box is described by standing de Broglie waves within the box, derive an equation for the allowed energies. Then use that equation to find the transition energy from n=1 to n=2 given the length of the box is 350 pm and the mass of an electron is $$9.109 \times 10^{-31} kg$$.

S3.33

The de Broglie formula is

$\lambda=\dfrac{h}{p}$

An integral number of half-wavelengths will fit in the box because the waves are standing waves so

$\dfrac{n\lambda}{2}=a$

$\dfrac{nh}{2p}=a$

Then solving for p

$p=\dfrac{nh}{2a}$

Therefore the energy equation is

$E=\dfrac{mv^2}{2}=\dfrac{p^2}{2,m}=\dfrac{1}{2m} \dfrac{n^2h^2}{4a^2}=\dfrac{n^2h^2}{8ma^2}$

Just plug into the equation to find the transition energy

$\Delta E=\dfrac{h^2}{8m_ea^2}(2^2-1^2)$

$\Delta E=\dfrac{(6.626 \times 10^{-34} J \centerdot s)^2 (3) }{8(9.109 \times 10^{-31} kg)(350 \times 10^{-12}m)^2}$

$\Delta E=1.47 \times 10^{-18}J$

Q3.35

Consider the two wavefunctions

$$\hat{\psi}_n(x) = sin\dfrac{n\pi x}{a}$$with even n numbers and

$$\hat{\psi}_n(x) = cos\dfrac{n\pi x}{a}$$with odd n numbers, prove that the wavefunctions can be symmetric and antisymmetric by using the operation x to -x, a is a constant.

Given that the Schrodinger equation has the expression:

$$\hat{H}(x)\psi_n(x) = E_n\psi_n(x)$$

Through the operation x to -x, the equation now becomes:

$$\hat{H}(-x)\psi_n(-x) = E_n\psi_n(-x)$$

Show that $$\hat{H}(x) = \hat{H}(-x)$$ is true to prove the Schrodinger equation.

S3.35

Substituting x by -x, for odd n numbers,

$$\hat{\psi}_n(-x) = cos\dfrac{-n\pi x}{a} = cos\dfrac{n\pi x}{a} = \hat{\psi}_n(x)$$

For even n numbers,

$$\hat{\psi}_n(-x) = sin\dfrac{-n\pi x}{a} = -sin\dfrac{n\pi x}{a} = -\hat{\psi}_n(x)$$

Thus, the wavefunction for odd n number is symmetric and even n numbers is antisymmetric.

And,

$$\hat{H}(x) = -\dfrac{\hbar^2}{2m}\dfrac{d^2}{dx^2} = \hat{H}(x)$$

$$\hat{H}(-x) = -\dfrac{\hbar^2}{2m}\dfrac{d^2}{d(-x)^2} = \hat{H}(x)$$

Thus,

$$\hat{H}(x) = \hat{H}(-x)$$ and $$\hat{H}(x)$$

is an even function of x.

Q3.35

Show that the Hamiltonian for a Rigid Rotor Model is odd.

S3.35

$\hat{H}(x) = \hat{H}(-x)$
so

$\hat{H}= -\dfrac{h^2}{4\pi \mu}\nabla^2$

$\nabla^2 = \dfrac{d^2y}{dx^2} \dfrac{d^2y}{dy^2}\dfrac{d^2y}{dz^2}$

so

$\dfrac{d^2y}{dx^2} (x) = 0$

and

$\dfrac{d^2y}{dx^2}(-x) = 0$

so

$\hat{H}(x) = 0$

and

$\hat{H}(-x) = 0$

so

$\hat{H}(x)= \hat{H}(-x)$