Chapter 2 • Theory & Derivations

Unit 2: Exact Quantum Models I: Particle in a Box & Ring

Rigorous quantum mechanical formulation of bounded particle models: the one-dimensional infinite potential well, Dirichlet boundary conditions, energy quantization, zero-point energy, symmetry-enforced parity, extension to 2D and 3D cuboidal and cubic boxes, spatial degeneracy, density of states, electron on a ring under periodic boundary conditions, angular momentum quantization, and the Free Electron Molecular Orbital (FEMO) model applied to conjugated polyenes and aromatic systems.

§2.1 One-Dimensional Infinite Potential Well: Boundary Conditions & Quantization

The particle in a one-dimensional box represents the simplest exact solvable quantum system exhibiting boundary-condition-induced energy quantization.

Mathematical Formulation

A point particle of mass \(m\) is constrained to move along the \(x\)-axis between rigid boundaries at \(x = 0\) and \(x = L\). The potential energy function is:

V(x) = \begin{cases} 0 & 0 < x < L \\ \infty & x \le 0 \text{ or } x \ge L \end{cases}

Inside the box (\(0 < x < L\)), the TISE is:

-\frac{\hbar^2}{2m} \frac{d^2\psi(x)}{dx^2} = E \psi(x) \implies \frac{d^2\psi(x)}{dx^2} + k^2 \psi(x) = 0

where the wavenumber is \(k = \frac{\sqrt{2mE}}{\hbar}\). The general solution is:

\psi(x) = A \sin(kx) + B \cos(kx)

Application of Dirichlet Boundary Conditions

To prevent infinite potential energy, the wavefunction must vanish identically outside the well: \(\psi(x) = 0\) for \(x \le 0\) and \(x \ge L\). Continuity requires:

1. At \(x = 0\):

\psi(0) = A \sin(0) + B \cos(0) = B = 0

Thus, \(\psi(x) = A \sin(kx)\).

2. At \(x = L\):

\psi(L) = A \sin(k L) = 0

To avoid the trivial solution \(\psi(x) \equiv 0\) (no particle exists), we require:

k L = n \pi, \quad n \in \{1, 2, 3, \dots\}

The quantum number \(n = 0\) is disallowed because it yields \(\psi(x) \equiv 0\), violating normalization. Negative integers \(n < 0\) yield identical physical states up to an unobservable global phase factor (\(\sin(-x) = -\sin(x)\)).

Quantized Energy Eigenvalues

Equating \(k = \frac{n\pi}{L} = \frac{\sqrt{2mE}}{\hbar}\) yields the quantized energy spectrum:

E_n = \frac{n^2 \pi^2 \hbar^2}{2 m L^2} = \frac{n^2 h^2}{8 m L^2}, \quad n = 1, 2, 3, \dots

Normalization of Eigenfunctions

Integrating the probability density over the box:

\int_0^L |\psi_n(x)|^2 dx = |A|^2 \int_0^L \sin^2\left(\frac{n\pi x}{L}\right) dx = |A|^2 \frac{L}{2} = 1 \implies A = \sqrt{\frac{2}{L}}

Thus, the normalized stationary eigenfunctions are:

\psi_n(x) = \sqrt{\frac{2}{L}} \sin\left(\frac{n\pi x}{L}\right)

Physical Hallmarks

Zero-Point Energy: The ground state (\(n = 1\)) has non-zero energy \(E_1 = \frac{h^2}{8 m L^2} > 0\). A particle confined to a finite domain cannot be at rest (\(E = 0\)), as having zero energy would imply \(\Delta p = 0\) while \(\Delta x \le L\), violating \(\Delta x \cdot \Delta p \ge \hbar/2\).
Nodal Structure: State \(\psi_n(x)\) has exactly \(n - 1\) internal nodes (zeros inside the box). Higher kinetic energy directly correlates with higher spatial curvature (\(d^2\psi/dx^2\)) and more nodes.
Classical Correspondence: As \(n \rightarrow \infty\), the spacing between successive nodes becomes microscopic. The rapid spatial oscillation averages out such that \(\langle |\psi_n(x)|^2 \rangle \approx 1/L\), recovering the classical uniform probability distribution (Bohr's correspondence principle).

§2.2 Symmetry, Parity & The Symmetric Infinite Well

Choosing coordinate origins that exploit physical symmetries greatly simplifies matrix element evaluations and reveals conservation laws.

The Center-Symmetric Box

Let the box be centered at the origin, spanning the interval \([-a, +a]\) with total width \(L = 2a\). The potential is:

V(x) = \begin{cases} 0 & -a < x < a \\ \infty & |x| \ge a \end{cases}

Notice that the Hamiltonian is invariant under spatial inversion (parity operator \(\hat{\Pi} x = -x\)):

\hat{H}(-x) = \hat{H}(x) \implies [\hat{H}, \hat{\Pi}] = 0

Because \(\hat{H}\) and \(\hat{\Pi}\) commute, the non-degenerate stationary states must be simultaneous eigenfunctions of parity with eigenvalues \(\pm 1\):

\psi(-x) = \pm \psi(x)

Solving the Symmetric Well

The general solution inside the well is:

\psi(x) = A \cos(kx) + B \sin(kx)

1. Even Parity Solutions (\(\psi(-x) = +\psi(x) \implies B = 0\)):

Boundary condition at \(x = \pm a\):

\cos(k a) = 0 \implies k a = \frac{n \pi}{2}, \quad n = 1, 3, 5, \dots

\psi_n(x) = \frac{1}{\sqrt{a}} \cos\left(\frac{n\pi x}{2a}\right), \quad n = 1, 3, 5, \dots

2. Odd Parity Solutions (\(\psi(-x) = -\psi(x) \implies A = 0\)):

Boundary condition at \(x = \pm a\):

\sin(k a) = 0 \implies k a = \frac{n \pi}{2}, \quad n = 2, 4, 6, \dots

\psi_n(x) = \frac{1}{\sqrt{a}} \sin\left(\frac{n\pi x}{2a}\right), \quad n = 2, 4, 6, \dots

Combining both branches with \(L = 2a\):

E_n = \frac{n^2 \pi^2 \hbar^2}{2 m (2a)^2} = \frac{n^2 h^2}{8 m L^2}, \quad n = 1, 2, 3, \dots

The energy spectrum is identical, but the alternating even/odd parity character drastically simplifies dipole transition selection rules: \(\langle \psi_m | \hat{x} | \psi_n \rangle = 0\) whenever \(m\) and \(n\) possess identical parity.

§2.3 Two- and Three-Dimensional Boxes: Separation of Variables & Degeneracy

Extending confinement to higher spatial dimensions introduces the concept of quantum degeneracy, where multiple distinct wavefunctions share identical energy eigenvalues.

The Three-Dimensional Rectangular Box

Consider a particle of mass \(m\) confined within a 3D cuboidal cavity of dimensions \(L_x, L_y, L_z\):

V(x, y, z) = \begin{cases} 0 & 0 < x < L_x, \ 0 < y < L_y, \ 0 < z < L_z \\ \infty & \text{otherwise} \end{cases}

Inside the box, the TISE is:

-\frac{\hbar^2}{2m} \left( \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2} \right) \psi(x, y, z) = E \psi(x, y, z)

Using separation of variables \(\psi(x, y, z) = X(x) Y(y) Z(z)\):

-\frac{\hbar^2}{2m} \left( \frac{X''}{X} + \frac{Y''}{Y} + \frac{Z''}{Z} \right) = E

Each independent coordinate satisfies an identical 1D equation:

-\frac{\hbar^2}{2m} X''(x) = E_x X(x), \quad -\frac{\hbar^2}{2m} Y''(y) = E_y Y(y), \quad -\frac{\hbar^2}{2m} Z''(z) = E_z Z(z)

with \(E = E_x + E_y + E_z\).

Quantized Eigenvalues and Wavefunctions

Applying Dirichlet boundary conditions on each axis:

\psi_{n_x, n_y, n_z}(x, y, z) = \sqrt{\frac{8}{L_x L_y L_z}} \sin\left(\frac{n_x \pi x}{L_x}\right) \sin\left(\frac{n_y \pi y}{L_y}\right) \sin\left(\frac{n_z \pi z}{L_z}\right)

E_{n_x, n_y, n_z} = \frac{h^2}{8m} \left( \frac{n_x^2}{L_x^2} + \frac{n_y^2}{L_y^2} + \frac{n_z^2}{L_z^2} \right), \quad n_x, n_y, n_z \in \{1, 2, 3, \dots\}

Symmetry and Systematic Degeneracy in a Cubic Box

When the box possesses cubic symmetry (\(L_x = L_y = L_z = L\)):

E_{n_x, n_y, n_z} = \frac{h^2}{8 m L^2} (n_x^2 + n_y^2 + n_z^2)

1. Ground State (1, 1, 1):

\(E_{1,1,1} = 3 \frac{h^2}{8mL^2}\). Non-degenerate (degeneracy \(g = 1\)).

2. First Excited State (2, 1, 1), (1, 2, 1), (1, 1, 2):

\(E = (4 + 1 + 1) \frac{h^2}{8mL^2} = 6 \frac{h^2}{8mL^2}\). Three distinct states share identical energy; triply degenerate (\(g = 3\)).

3. Second Excited State (2, 2, 1), (2, 1, 2), (1, 2, 2):

\(E = (4 + 4 + 1) \frac{h^2}{8mL^2} = 9 \frac{h^2}{8mL^2}\). Triply degenerate (\(g = 3\)).

4. Accidental Degeneracy:

Higher energy levels may exhibit degeneracies not dictated by cubic spatial symmetry alone. For example, \((3, 3, 3)\) with \(n_x^2 + n_y^2 + n_z^2 = 27\) can mix with different sums of three squares, termed accidental degeneracy.

§2.4 Electron in a Ring: Periodic Boundary Conditions & Aromaticity

When a particle is constrained to move along a circular path of radius \(R\), the spatial domain is closed and continuous, replacing Dirichlet boundary conditions with periodic boundary conditions.

The Schrödinger Equation in Polar Coordinates

Consider an electron of mass \(m\) moving on a planar circular ring of radius \(R\) (circumference \(L = 2\pi R\)) in the \(xy\)-plane with \(V(\phi) = 0\). The arc length coordinate is \(s = R \phi\). The kinetic energy operator in angular coordinate \(\phi\) is:

\hat{H} = -\frac{\hbar^2}{2m R^2} \frac{d^2}{d\phi^2} = -\frac{\hbar^2}{2 I} \frac{d^2}{d\phi^2}

where \(I = m R^2\) is the moment of inertia. The TISE is:

-\frac{\hbar^2}{2 I} \frac{d^2\psi(\phi)}{d\phi^2} = E \psi(\phi) \implies \frac{d^2\psi(\phi)}{d\phi^2} + m_l^2 \psi(\phi) = 0

where \(m_l = \frac{\sqrt{2 I E}}{\hbar}\). The general solution is:

\psi(\phi) = A e^{i m_l \phi} + B e^{-i m_l \phi}

Periodic Boundary Conditions

Because \(\phi\) and \(\phi + 2\pi\) represent the identical physical point in space:

\psi(\phi + 2\pi) = \psi(\phi) \implies e^{i m_l (\phi + 2\pi)} = e^{i m_l \phi} \implies e^{i 2\pi m_l} = 1

Euler's identity requires:

m_l \in \{0, \pm 1, \pm 2, \pm 3, \dots\}

Notice that \(m_l = 0\) is physically allowed here (unlike the 1D box) because \(e^0 = 1\) is non-zero and normalizable.

Quantized Energy Levels & Degeneracy

The energy eigenvalues are:

E_{m_l} = \frac{m_l^2 \hbar^2}{2 I} = \frac{m_l^2 \hbar^2}{2 m R^2}

Ground state (\(m_l = 0\)): \(E_0 = 0\). The electron has zero angular momentum and zero zero-point energy (non-degenerate, \(g = 1\)).
Excited states (\(|m_l| \ge 1\)): Every energy level is doubly degenerate (\(g = 2\)) corresponding to clockwise (\(+m_l\)) and counter-clockwise (\(-m_l\)) circulating currents.

Normalization

\int_0^{2\pi} |\psi(\phi)|^2 d\phi = |A|^2 \int_0^{2\pi} d\phi = 2\pi |A|^2 = 1 \implies A = \frac{1}{\sqrt{2\pi}}

Thus:

\psi_{m_l}(\phi) = \frac{1}{\sqrt{2\pi}} e^{i m_l \phi}

These functions are simultaneous eigenfunctions of the \(z\)-component of orbital angular momentum:

\hat{L}_z = -i\hbar \frac{\partial}{\partial \phi} \implies \hat{L}_z \psi_{m_l}(\phi) = m_l \hbar \psi_{m_l}(\phi)

§2.5 Free Electron Molecular Orbital (FEMO) Model for Conjugated Polyenes

The Free Electron Molecular Orbital (FEMO) model, developed by Hans Kuhn (1949), treats \(\pi\)-electrons in conjugated polyenes as independent particles in an effective one-dimensional box.

Physical Model Assumptions

In a conjugated polyene \(\text{H}_2\text{C}=\text{CH}-(\text{CH}=\text{CH})_k-\text{CH}=\text{CH}_2\) containing \(N\) conjugated carbon atoms:

The \(\sigma\)-bonding skeleton provides a uniform, flat potential well of length \(L\).
The \(N\) \(\pi\)-electrons move freely along the carbon backbone without inter-electronic repulsion.
The box length \(L\) is approximated by the sum of carbon-carbon bond lengths plus boundary extensions \(\delta\) at the terminal carbons (accounting for \(\pi\)-cloud spill-over):

L = (N - 1) d_{\text{C-C}} + 2 \delta \approx N d_{\text{avg}}

where \(d_{\text{avg}} \approx 1.39 - 1.40\text{ Å}\) is the average resonant \(\text{C-C}\) bond distance.

Aufbau Principle & HOMO-LUMO Transition

According to the Pauli exclusion principle, each spatial orbital \(\psi_n\) accommodates at most two electrons with antiparallel spins (\(m_s = \pm 1/2\)). For an even number of \(\pi\)-electrons \(N\):

The Highest Occupied Molecular Orbital (HOMO) corresponds to quantum number:

n_{\text{HOMO}} = \frac{N}{2}

The Lowest Unoccupied Molecular Orbital (LUMO) corresponds to quantum number:

n_{\text{LUMO}} = \frac{N}{2} + 1

The lowest-energy electronic absorption transition is the \(\text{HOMO} \rightarrow \text{LUMO}\) excitation:

\Delta E = E_{\text{LUMO}} - E_{\text{HOMO}} = \frac{h^2}{8 m_e L^2} \left[ \left(\frac{N}{2} + 1\right)^2 - \left(\frac{N}{2}\right)^2 \right] = \frac{h^2}{8 m_e L^2} (N + 1)

Expressing the absorption wavelength \(\lambda\):

\Delta E = \frac{h c}{\lambda} \implies \lambda = \frac{8 m_e c L^2}{h (N + 1)}

Substituting \(L \approx N d\):

\lambda \approx \frac{8 m_e c d^2}{h} \frac{N^2}{N + 1}

As conjugation length \(N\) increases, \(\lambda\) shifts bathochromically into the visible spectrum, explaining the intense colors of carotenoids (\(\beta\)-carotene, \(N=22\), orange) and cyanine dyes.

§2.6 Quantum Mechanical Tunneling: Finite Potential Wells & Barriers

When the potential walls have finite height \(V_0 < \infty\), quantum wavefunctions penetrate into classically forbidden regions where \(E < V(x)\), giving rise to quantum mechanical tunneling.

The Finite One-Dimensional Square Well

Consider a particle of mass \(m\) in a finite well:

V(x) = \begin{cases} 0 & -a \le x \le a \\ V_0 & |x| > a \end{cases}

For bound states with \(0 < E < V_0\):

Inside the Well (\(|x| < a\)):

\frac{d^2\psi}{dx^2} + k^2 \psi = 0, \quad k = \frac{\sqrt{2mE}}{\hbar}

Solutions are oscillatory: \(\cos(kx)\) (even) or \(\sin(kx)\) (odd).

Outside the Well (\(|x| > a\)):

\frac{d^2\psi}{dx^2} - \kappa^2 \psi = 0, \quad \kappa = \frac{\sqrt{2m(V_0 - E)}}{\hbar}

Solutions are exponentially decaying evanescent waves:

\psi(x) = C e^{-\kappa x} \quad (x > a), \quad \psi(x) = C e^{+\kappa x} \quad (x < -a)

Continuity Conditions and Penetration Depth

Matching \(\psi(x)\) and its first derivative \(\psi'(x)\) at the boundary \(x = a\) for even states yields:

k \tan(k a) = \kappa

The characteristic penetration depth \(\delta\) where the wavefunction decays to \(1/e\) of its boundary value is:

\delta = \frac{1}{\kappa} = \frac{\hbar}{\sqrt{2m(V_0 - E)}}

Because \(\psi(x) \neq 0\) outside the well, there is a finite probability of finding the particle in the classically forbidden zone.

Rectangular Potential Barrier & Transmission Coefficient

For a barrier of width \(L\) and height \(V_0 > E\), the quantum transmission coefficient (tunneling probability) is:

T = \frac{1}{1 + \frac{V_0^2}{4 E (V_0 - E)} \sinh^2(\kappa L)}

In the thick/high barrier limit (\(\kappa L \gg 1\)):

T \approx 16 \frac{E}{V_0} \left(1 - \frac{E}{V_0}\right) \exp(-2\kappa L)

Tunneling probability decays exponentially with barrier width \(L\) and the square root of particle mass \(\sqrt{m}\). This explains why proton tunneling occurs readily in enzyme active sites and ammonia inversion, whereas deuteron tunneling is suppressed by an order of magnitude.

§2.7 Density of States & Transition to the Classical Continuum

In macroscopic condensed matter and statistical mechanics, quantum systems contain billions of closely spaced energy levels, requiring a continuous density of states representation.

Derivation of the Density of States in 3D

For a particle in a 3D cubic box of volume \(V = L^3\), the energy levels satisfy:

E = \frac{h^2}{8 m L^2} (n_x^2 + n_y^2 + n_z^2) = \frac{\hbar^2 \pi^2}{2 m L^2} n^2

where \(n^2 = n_x^2 + n_y^2 + n_z^2\). In the three-dimensional quantum number space \((n_x, n_y, n_z)\), each quantum state occupies a unit volume of \(1 \times 1 \times 1 = 1\). Because \(n_x, n_y, n_z > 0\), the states are restricted to the positive octant (\(1/8\) of a sphere). The total number of states with radius less than \(n\) is:

N(E) = \frac{1}{8} \left( \frac{4}{3} \pi n^3 \right) = \frac{\pi}{6} n^3

Expressing \(n\) in terms of energy \(E\):

n = \frac{L}{\pi \hbar} \sqrt{2 m E} \implies N(E) = \frac{\pi}{6} \left( \frac{L}{\pi \hbar} \sqrt{2m E} \right)^3 = \frac{V}{6\pi^2 \hbar^3} (2m E)^{3/2}

The density of states \(g(E) = \frac{dN}{dE}\) is the number of available quantum states per unit energy:

g(E) = \frac{dN}{dE} = \frac{V}{4\pi^2 \hbar^3} (2m)^{3/2} E^{1/2} = \frac{2\pi V}{h^3} (2m)^{3/2} \sqrt{E}

Dimensionality and Density of States Scaling

The energy dependence of the density of states depends fundamentally on the spatial dimensionality \(d\) of confinement:

3D Bulk: \(g(E) \propto E^{1/2}\) (continuous parabolic increase)
2D Quantum Well: \(g(E) \propto E^0 = \text{constant}\) (step-like staircase)
1D Quantum Wire: \(g(E) \propto E^{-1/2}\) (van Hove singularities)
0D Quantum Dot: \(g(E) = \sum_i \delta(E - E_i)\) (discrete delta peaks)

This dimensional scaling underpins semiconductor quantum well lasers, carbon nanotubes, and colloidal nanocrystal quantum dots.

§2.8 Symmetric Double-Well Potentials & Quantum Inversion Tunneling

Double-well potential energy surfaces represent the archetypal model for chemical isomerization, proton transfer in hydrogen bonds, and pyramidal inversion in molecules such as ammonia (\(\text{NH}_3\)).

The Double-Well Hamiltonian

Consider a symmetric one-dimensional potential with two degenerate classical minima at \(x = \pm x_0\) separated by a central potential barrier of height \(V_0\) at \(x = 0\):

V(x) = \frac{\lambda}{4} (x^2 - x_0^2)^2 = -a x^2 + b x^4 + V_0

Classically, a particle with energy \(E < V_0\) is trapped forever in either the left well (\(x < 0\)) or right well (\(x > 0\)).

Parity and Energy Level Splitting

Because the Hamiltonian is symmetric under spatial inversion \(\hat{\Pi}: x \rightarrow -x\), \([\hat{H}, \hat{\Pi}] = 0\). All non-degenerate eigenstates must possess definite parity:

1. Symmetric (Even parity) Ground State: \(\psi_S(-x) = +\psi_S(x)\)

2. Antisymmetric (Odd parity) First Excited State: \(\psi_A(-x) = -\psi_A(x)\)

Let \(|L\rangle\) and \(|R\rangle\) denote localized quasi-ground states in the left and right wells. The true stationary eigenstates are quantum superpositions:

|S\rangle = \frac{1}{\sqrt{2}} (|L\rangle + |R\rangle), \quad E_S = E_0 - \frac{\Delta E}{2}

|A\rangle = \frac{1}{\sqrt{2}} (|L\rangle - |R\rangle), \quad E_A = E_0 + \frac{\Delta E}{2}

The energy difference \(\Delta E = E_A - E_S\) is the tunnel splitting.

Tunneling Dynamics and Inversion Frequency

If the system is prepared at \(t = 0\) localized entirely in the left well:

|\Psi(0)\rangle = |L\rangle = \frac{1}{\sqrt{2}} (|S\rangle + |A\rangle)

Its time evolution is:

|\Psi(t)\rangle = \frac{1}{\sqrt{2}} \left[ |S\rangle e^{-i E_S t / \hbar} + |A\rangle e^{-i E_A t / \hbar} \right] = e^{-i E_0 t / \hbar} \left[ \cos\left(\frac{\Delta E \, t}{2\hbar}\right) |L\rangle - i \sin\left(\frac{\Delta E \, t}{2\hbar}\right) |R\rangle \right]

The probability of finding the particle in the right well oscillates periodically:

P_R(t) = |\langle R | \Psi(t) \rangle|^2 = \sin^2\left( \frac{\Delta E \, t}{2\hbar} \right)

The period of complete quantum tunneling through the barrier is \(\tau = \frac{\pi \hbar}{\Delta E}\), corresponding to the inversion tunneling frequency:

\nu_{\text{inv}} = \frac{\Delta E}{h}

For ammonia (\(\text{NH}_3\)), the nitrogen atom tunnels through the plane of three hydrogen atoms with \(\Delta E \approx 0.7935\text{ cm}^{-1}\), producing the celebrated microwave inversion frequency of \(\nu_{\text{inv}} \approx 23.79\text{ GHz}\), which was utilized in 1954 to build the world's first maser.

## Advanced Mathematical Supplement: Feynman Path Integral Formulation

The Feynman Path Integral Representation of Solvable Wells

In 1948, Richard Feynman reformulated quantum mechanics by replacing the operator differential equation with a sum over all possible classical trajectories connecting spacetime points \((x_a, t_a)\) and \((x_b, t_b)\).

The Quantum Propagator

The transition amplitude (quantum propagator) is:

K(x_b, t_b; x_a, t_a) = \langle x_b | e^{-i \hat{H} (t_b - t_a) / \hbar} | x_a \rangle

In Feynman's functional path integral formulation:

K(x_b, t_b; x_a, t_a) = \int \mathcal{D}[x(t)] \exp\left( \frac{i}{\hbar} S[x(t)] \right)

where the classical action functional is:

S[x(t)] = \int_{t_a}^{t_b} \left[ \frac{1}{2} m \dot{x}(t)^2 - V(x(t)) \right] dt

Method of Images for the Infinite Potential Well

For a particle in an infinite potential well of width \(L\) (\(V(x) = 0\) for \(0 < x < L\)), every path must remain confined strictly within the boundaries \(0 < x < L\). Using the method of images, the path integral over the box equals the sum over all classical trajectories of a free particle in an infinite periodic lattice with alternating mirror reflections at \(x = 0\) and \(x = L\):

K(x_b, t_b; x_a, 0) = \sum_{n=-\infty}^\infty \left[ K_{\text{free}}(x_b - x_a + 2nL, t) - K_{\text{free}}(-x_b - x_a + 2nL, t) \right]

where \(K_{\text{free}}(x, t) = \sqrt{\frac{m}{2\pi i \hbar t}} \exp\left( \frac{i m x^2}{2\hbar t} \right)\). By applying the Poisson summation formula \(\sum_{n=-\infty}^\infty e^{2\pi i n k} = \sum_{m=-\infty}^\infty \delta(k - m)\), this sum over topological winding numbers converts identically into the eigenfunction expansion:

K(x_b, t; x_a, 0) = \sum_{n=1}^\infty \frac{2}{L} \sin\left(\frac{n\pi x_b}{L}\right) \sin\left(\frac{n\pi x_a}{L}\right) \exp\left( -\frac{i E_n t}{\hbar} \right)

with \(E_n = \frac{n^2 \pi^2 \hbar^2}{2 m L^2}\). The path integral establishes that the discrete quantum energy spectrum arises from the constructive interference of an infinite family of virtual classical trajectories reflecting between the boundaries.

## Research Monograph: Quantum Chaos & Wigner-Dyson Random Matrix Spectral Statistics

Semiclassical Mechanics in Non-Integrable Systems

In classical mechanics, systems are either integrable (possessing as many independent constants of motion in involution as degrees of freedom, such as the circular billiard) or chaotic (such as the Sinai or Bunimovich stadium billiard, where adjacent trajectories diverge exponentially with positive Lyapunov exponent \(\lambda_L > 0\)). In quantum mechanics, the Schrödinger equation is strictly linear, meaning quantum wavepackets cannot diverge exponentially indefinitely. The signature of classical chaos in quantum mechanics is termed Quantum Chaos.

The Bohigas-Giannoni-Schmit (BGS) Conjecture

In 1984, Oriol Bohigas, Marie-Joya Giannoni, and Charles Schmit conjectured a profound universality:

1. Integrable Systems (Berry-Tabor Conjecture):

The energy levels are uncorrelated. The probability distribution \(P(s)\) of normalized adjacent energy level spacings \(s = (E_{n+1} - E_n) / \langle \Delta E \rangle\) follows a Poisson distribution:

P_{\text{Poisson}}(s) = e^{-s}

Notice that \(P(0) = 1\): level clustering is favored; states can be arbitrarily close together without interacting.

2. Chaotic Systems (BGS Conjecture):

The quantum energy spectra exhibit universal level repulsion (\(P(0) = 0\)). The spectral statistics match the eigenvalue distributions of Random Matrix Theory (RMT):

For systems with time-reversal symmetry (Gaussian Orthogonal Ensemble, GOE):

P_{\text{GOE}}(s) = \frac{\pi}{2} s \exp\left( -\frac{\pi}{4} s^2 \right) \quad (\text{Wigner surmise with linear repulsion } P(s) \propto s)

For systems with broken time-reversal symmetry (Gaussian Unitary Ensemble, GUE):

P_{\text{GUE}}(s) = \frac{32}{\pi^2} s^2 \exp\left( -\frac{4}{\pi} s^2 \right) \quad (\text{Quadratic repulsion } P(s) \propto s^2)

Wavefunction Scarring

Martin Gutzwiller's trace formula connects the quantum density of states directly to the periodic classical orbits of the system:

\rho(E) = \bar{\rho}(E) + \frac{1}{\pi\hbar} \sum_{\gamma} \sum_{k=1}^\infty \frac{T_\gamma}{\sqrt{|\det(\mathbf{M}_\gamma^k - \mathbf{I})|}} \cos\left( \frac{k S_\gamma}{\hbar} - \frac{k \mu_\gamma \pi}{2} \right)

Eric Heller discovered that high-energy quantum eigenstates of classically chaotic billiards are not uniformly ergodic; instead, they exhibit intense, localized standing-wave ridges along the shortest unstable classical periodic orbits—a quantum interference phenomenon known as wavefunction scarring.

Worked Problems & Step-by-Step Quantum Derivations

Multi-step solved problems covering Planck distribution, photoelectric kinetics, Compton shift, de Broglie wavelengths, uncertainty relations, and Hermitian operator commutation algebra.

Foundational Example 2.1: Problem 1: Ground State & Excitation Energies for an Electron in a Nanoscale Well

An electron is trapped in a one-dimensional infinite potential well of width \(L = 0.500\text{ nm}\).

Calculate the ground-state zero-point energy \(E_1\) in Joules and electron-volts (\(\text{eV}\)).
Determine the energy of the first two excited states (\(E_2\) and \(E_3\)).
Calculate the wavelength \(\lambda\) of the photon absorbed when the electron undergoes a transition from \(n = 1 \rightarrow n = 2\).

Comprehensive Multi-Step Solution:

Step 1: Ground-State Zero-Point Energy

The energy levels for a 1D infinite well are:

E_n = \frac{n^2 h^2}{8 m_e L^2}

Substitute \(n = 1\), \(h = 6.62607 \times 10^{-34}\text{ J}\cdot\text{s}\), \(m_e = 9.10938 \times 10^{-31}\text{ kg}\), and \(L = 0.500 \times 10^{-9}\text{ m}\):

E_1 = \frac{(1)^2 (6.62607 \times 10^{-34}\text{ J}\cdot\text{s})^2}{8 (9.10938 \times 10^{-31}\text{ kg})(0.500 \times 10^{-9}\text{ m})^2}

E_1 = \frac{4.39048 \times 10^{-67}}{1.82188 \times 10^{-48}} = 2.40986 \times 10^{-19}\text{ J}

In electron-volts:

E_1 = \frac{2.40986 \times 10^{-19}\text{ J}}{1.60218 \times 10^{-19}\text{ J/eV}} \approx 1.5041\text{ eV} \approx 1.504\text{ eV}

Step 2: Excited State Energies

Since \(E_n = n^2 E_1\):

For \(n = 2\):

E_2 = 2^2 E_1 = 4 \times 1.5041\text{ eV} = 6.0164\text{ eV} \approx 6.016\text{ eV}

For \(n = 3\):

E_3 = 3^2 E_1 = 9 \times 1.5041\text{ eV} = 13.5369\text{ eV} \approx 13.537\text{ eV}

Step 3: Transition Wavelength for \(n = 1 \rightarrow 2\)

The transition energy is:

\Delta E = E_2 - E_1 = 3 E_1 = 3 \times 2.40986 \times 10^{-19}\text{ J} = 7.22958 \times 10^{-19}\text{ J} = 4.5123\text{ eV}

The absorbed photon wavelength is:

\lambda = \frac{h c}{\Delta E} = \frac{(6.62607 \times 10^{-34}\text{ J}\cdot\text{s})(2.99792 \times 10^8\text{ m/s})}{7.22958 \times 10^{-19}\text{ J}} \approx 2.7476 \times 10^{-7}\text{ m} = 274.8\text{ nm}

This transition lies in the near-ultraviolet spectrum.

Foundational Example 2.2: Problem 2: Spatial Probability Integration & Central Region Locality

For a particle in a one-dimensional box of length \(L\) in its ground state \(\psi_1(x) = \sqrt{\frac{2}{L}} \sin\left(\frac{\pi x}{L}\right)\):

Calculate the probability of finding the particle in the central third of the box, \(x \in [L/3, 2L/3]\).
Compare this result to the classical probability for a uniform distribution and explain the physical difference.

Comprehensive Multi-Step Solution:

Step 1: Quantum Probability Calculation

The probability is given by integrating \(|\psi_1(x)|^2\):

P = \int_{L/3}^{2L/3} |\psi_1(x)|^2 dx = \frac{2}{L} \int_{L/3}^{2L/3} \sin^2\left(\frac{\pi x}{L}\right) dx

Using the trigonometric identity \(\sin^2\theta = \frac{1 - \cos(2\theta)}{2}\):

P = \frac{2}{L} \int_{L/3}^{2L/3} \frac{1 - \cos\left(\frac{2\pi x}{L}\right)}{2} dx = \frac{1}{L} \left[ x - \frac{L}{2\pi} \sin\left(\frac{2\pi x}{L}\right) \right]_{L/3}^{2L/3}

Evaluate at the upper limit \(x = 2L/3\):

\text{Upper} = \frac{2L}{3} - \frac{L}{2\pi} \sin\left(\frac{4\pi}{3}\right) = \frac{2L}{3} - \frac{L}{2\pi} \left(-\frac{\sqrt{3}}{2}\right) = \frac{2L}{3} + \frac{\sqrt{3} L}{4\pi}

Evaluate at the lower limit \(x = L/3\):

\text{Lower} = \frac{L}{3} - \frac{L}{2\pi} \sin\left(\frac{2\pi}{3}\right) = \frac{L}{3} - \frac{L}{2\pi} \left(\frac{\sqrt{3}}{2}\right) = \frac{L}{3} - \frac{\sqrt{3} L}{4\pi}

Subtracting the lower limit from the upper limit:

P = \frac{1}{L} \left[ \left(\frac{2L}{3} - \frac{L}{3}\right) + \frac{\sqrt{3} L}{2\pi} \right] = \frac{1}{L} \left[ \frac{L}{3} + \frac{\sqrt{3} L}{2\pi} \right] = \frac{1}{3} + \frac{\sqrt{3}}{2\pi}

Evaluating numerically:

P \approx 0.33333 + \frac{1.73205}{6.28318} = 0.33333 + 0.27566 = 0.60899 \approx 0.609\text{ (or } 60.9\%\text{)}

Step 2: Classical Comparison

For a classical particle bouncing back and forth at constant speed:

P_{\text{class}} = \frac{\Delta x}{L} = \frac{2L/3 - L/3}{L} = \frac{1}{3} \approx 0.333\text{ (or } 33.3\%\text{)}

The quantum probability (\(60.9\%\)) is nearly double the classical expectation because the ground-state wavefunction forms a standing half-wave with maximum probability amplitude at the center of the well (\(x = L/2\)).

Intermediate Example 2.3: Problem 3: Degeneracy and Level Spectrum in a 2D Square Box

A particle of mass \(m\) is confined within a two-dimensional square box of length \(L\) along both \(x\) and \(y\).

Write the expression for the energy levels \(E_{n_x, n_y}\) in units of \(E_0 = \frac{h^2}{8 m L^2}\).
Determine the energies and degeneracies of the lowest six energy levels.
If the square well is deformed into a rectangle with \(L_x = L\) and \(L_y = 2L\), calculate the new energies of the states that were originally degenerate with \(n_x^2 + n_y^2 = 5\).

Comprehensive Multi-Step Solution:

Step 1: Energy Level Expression for 2D Square Box

E_{n_x, n_y} = \frac{h^2}{8 m L^2} (n_x^2 + n_y^2) = E_0 (n_x^2 + n_y^2), \quad n_x, n_y \in \{1, 2, 3, \dots\}

Step 2: Lowest Six Energy Levels and Degeneracies

Listing states by increasing values of \(n_x^2 + n_y^2\):

1. Level 1 (\(n_x=1, n_y=1\)):

\(n_x^2 + n_y^2 = 1 + 1 = 2 \implies E = 2 E_0\). Degeneracy \(g = 1\) (non-degenerate).

2. Level 2 (\(n_x=1, n_y=2\) and \(n_x=2, n_y=1\)):

\(n_x^2 + n_y^2 = 1 + 4 = 5 \implies E = 5 E_0\). Degeneracy \(g = 2\) (doubly degenerate).

3. Level 3 (\(n_x=2, n_y=2\)):

\(n_x^2 + n_y^2 = 4 + 4 = 8 \implies E = 8 E_0\). Degeneracy \(g = 1\).

4. Level 4 (\(n_x=1, n_y=3\) and \(n_x=3, n_y=1\)):

\(n_x^2 + n_y^2 = 1 + 9 = 10 \implies E = 10 E_0\). Degeneracy \(g = 2\).

5. Level 5 (\(n_x=2, n_y=3\) and \(n_x=3, n_y=2\)):

\(n_x^2 + n_y^2 = 4 + 9 = 13 \implies E = 13 E_0\). Degeneracy \(g = 2\).

6. Level 6 (\(n_x=1, n_y=4\) and \(n_x=4, n_y=1\)):

\(n_x^2 + n_y^2 = 1 + 16 = 17 \implies E = 17 E_0\). Degeneracy \(g = 2\).

Step 3: Symmetry Breaking in Rectangular Box (\(L_x = L, L_y = 2L\))

For a rectangular box:

E_{n_x, n_y} = \frac{h^2}{8m} \left( \frac{n_x^2}{L^2} + \frac{n_y^2}{(2L)^2} \right) = E_0 \left( n_x^2 + \frac{n_y^2}{4} \right)

Evaluating the two states that were previously degenerate at \(5 E_0\):

State \((1, 2)\):

E_{1, 2} = E_0 \left( 1^2 + \frac{2^2}{4} \right) = E_0 (1 + 1) = 2.00 E_0

State \((2, 1)\):

E_{2, 1} = E_0 \left( 2^2 + \frac{1^2}{4} \right) = E_0 (4 + 0.25) = 4.25 E_0

Deforming the square into a rectangle breaks the spatial reflection symmetry across the diagonal (\(x \leftrightarrow y\)), completely lifting the degeneracy.

Intermediate Example 2.4: Problem 4: FEMO Model of Octatetraene Absorption Wavelength

1,3,5,7-Octatetraene has \(N = 8\) conjugated \(\pi\)-electrons.

Assuming an average \(\text{C-C}\) bond distance \(d_{\text{avg}} = 0.140\text{ nm}\) and a box length \(L = (N - 1) d_{\text{avg}} + 2 \delta\) with terminal extension \(\delta = 0.070\text{ nm}\) (total \(L = 8 \times 0.140\text{ nm} = 1.120\text{ nm}\)):

Identify the quantum numbers of the HOMO and LUMO.
Calculate the \(\text{HOMO} \rightarrow \text{LUMO}\) transition energy \(\Delta E\) in \(\text{eV}\).
Calculate the predicted absorption wavelength \(\lambda_{\text{abs}}\) (in \(\text{nm}\)).

Compare this prediction with the experimental maximum absorption of octatetraene (\(\lambda_{\text{exp}} \approx 304\text{ nm}\)).

Comprehensive Multi-Step Solution:

Step 1: Identification of Frontier Orbitals

With \(N = 8\) \(\pi\)-electrons, each spatial level holds 2 electrons:

n_{\text{HOMO}} = \frac{N}{2} = \frac{8}{2} = 4

n_{\text{LUMO}} = n_{\text{HOMO}} + 1 = 5

Step 2: Transition Energy Calculation

The energy difference is:

\Delta E = E_5 - E_4 = \frac{h^2}{8 m_e L^2} (5^2 - 4^2) = \frac{h^2}{8 m_e L^2} (25 - 16) = 9 \frac{h^2}{8 m_e L^2}

Given \(L = 1.120\text{ nm} = 1.120 \times 10^{-9}\text{ m}\):

\frac{h^2}{8 m_e L^2} = \frac{(6.62607 \times 10^{-34}\text{ J}\cdot\text{s})^2}{8 (9.10938 \times 10^{-31}\text{ kg})(1.120 \times 10^{-9}\text{ m})^2} = \frac{4.39048 \times 10^{-67}}{9.1415 \times 10^{-48}} = 4.8028 \times 10^{-20}\text{ J}

In electron-volts:

\frac{h^2}{8 m_e L^2} = \frac{4.8028 \times 10^{-20}\text{ J}}{1.60218 \times 10^{-19}\text{ J/eV}} \approx 0.29977\text{ eV}

Therefore:

\Delta E = 9 \times 0.29977\text{ eV} \approx 2.6979\text{ eV} \approx 2.698\text{ eV}

In Joules:

\Delta E = 9 \times 4.8028 \times 10^{-20}\text{ J} = 4.3225 \times 10^{-19}\text{ J}

Step 3: Absorption Wavelength & Comparison

The theoretical absorption wavelength is:

\lambda_{\text{abs}} = \frac{h c}{\Delta E} = \frac{(6.62607 \times 10^{-34}\text{ J}\cdot\text{s})(2.99792 \times 10^8\text{ m/s})}{4.3225 \times 10^{-19}\text{ J}} \approx 4.595 \times 10^{-7}\text{ m} = 459.5\text{ nm}

Comparing with the experimental value \(\lambda_{\text{exp}} \approx 304\text{ nm}\): The simple FEMO model overestimates the wavelength because it assumes completely zero electron-electron Coulomb repulsion and equal bond lengths. In real octatetraene, bond length alternation (alternating double and single bonds) creates a periodic potential modulation that opens a wider bandgap, shifting the absorption to \(304\text{ nm}\).

Advanced Example 2.5: Problem 5: Electron in a Benzene Ring & Aromatic Delocalization Energy

Model the six \(\pi\)-electrons of a benzene ring (\(\text{C}_6\text{H}_6\)) as independent particles on a circular ring of radius \(R = 0.139\text{ nm}\).

Write the energy formula \(E_{m_l}\) in terms of \(\hbar, m_e, R\).
Populate the six \(\pi\)-electrons into the ring energy levels according to the Pauli principle and calculate the total \(\pi\)-electron ground-state energy \(E_{\text{total}}\) in \(\text{eV}\).
Determine the lowest electronic excitation energy \(\Delta E\) and predicted absorption wavelength \(\lambda\).
Explain how this circular model naturally rationalizes Hückel's \((4n + 2)\) aromatic stability rule.

Comprehensive Multi-Step Solution:

Step 1: Energy Expression on a Ring

The energy levels for an electron on a ring are:

E_{m_l} = \frac{m_l^2 \hbar^2}{2 m_e R^2}, \quad m_l = 0, \pm 1, \pm 2, \dots

Calculate the base energy unit \(\epsilon_0 = \frac{\hbar^2}{2 m_e R^2}\): With \(\hbar = 1.05457 \times 10^{-34}\text{ J}\cdot\text{s}\), \(m_e = 9.10938 \times 10^{-31}\text{ kg}\), and \(R = 0.139 \times 10^{-9}\text{ m}\):

\epsilon_0 = \frac{(1.05457 \times 10^{-34}\text{ J}\cdot\text{s})^2}{2 (9.10938 \times 10^{-31}\text{ kg})(0.139 \times 10^{-9}\text{ m})^2} = \frac{1.1121 \times 10^{-68}}{3.5209 \times 10^{-49}} = 3.1586 \times 10^{-20}\text{ J} \approx 0.19714\text{ eV}

Step 2: Ground-State Electronic Configuration

Populating 6 \(\pi\)-electrons:

Ground level \(m_l = 0\) (degeneracy \(g = 1\)): Holds 2 electrons (spin up, spin down).

E(m_l = 0) = 2 \times (0) = 0

First excited level \(m_l = \pm 1\) (degeneracy \(g = 2\)): Holds 4 electrons (2 in \(+1\), 2 in \(-1\)).

E(m_l = \pm 1) = 4 \times (1^2 \epsilon_0) = 4 \epsilon_0

Total ground-state energy:

E_{\text{total}} = 0 + 4 \epsilon_0 = 4 \times 0.19714\text{ eV} \approx 0.7886\text{ eV} \approx 1.263 \times 10^{-19}\text{ J}

Step 3: Lowest Electronic Transition

The HOMO is \(m_l = \pm 1\), and the LUMO is \(m_l = \pm 2\). The transition energy is:

\Delta E = E(m_l = 2) - E(m_l = 1) = (2^2 - 1^2) \epsilon_0 = 3 \epsilon_0

\Delta E = 3 \times 0.19714\text{ eV} \approx 0.5914\text{ eV} = 9.476 \times 10^{-20}\text{ J}

The transition wavelength is:

\lambda = \frac{h c}{\Delta E} = \frac{(6.62607 \times 10^{-34}\text{ J}\cdot\text{s})(2.99792 \times 10^8\text{ m/s})}{9.476 \times 10^{-20}\text{ J}} \approx 2.096 \times 10^{-6}\text{ m} = 2096\text{ nm}

Step 4: Physical Rationale of Hückel's \((4n+2)\) Rule

In any circular potential:

The lowest level (\(m_l = 0\)) is single-fold degenerate, accommodating exactly 2 electrons.
All subsequent energy levels (\(|m_l| = 1, 2, 3, \dots\)) are doubly degenerate (\(g = 2\)), each holding exactly 4 electrons (2 pairs).
To achieve a closed-shell electronic configuration with all occupied degenerate shells completely filled, the total number of \(\pi\)-electrons must be:

N_\pi = 2 + 4n = 4n + 2, \quad n \in \{0, 1, 2, \dots\}

This provides a direct physical derivation of Hückel's aromaticity rule.

Advanced Example 2.6: Problem 6: Quantum Tunneling Transmission Coefficient Through a Rectangular Barrier

An electron with kinetic energy \(E = 2.00\text{ eV}\) approaches a rectangular potential energy barrier of height \(V_0 = 5.00\text{ eV}\) and width \(L = 0.200\text{ nm}\).

Calculate the decay constant \(\kappa\) inside the barrier (in \(\text{m}^{-1}\)) and the characteristic penetration depth \(\delta\).
Calculate the exact quantum transmission probability \(T\) through the barrier.
Determine how much the transmission probability drops if the barrier width is doubled to \(L = 0.400\text{ nm}\).

Comprehensive Multi-Step Solution:

Step 1: Decay Constant and Penetration Depth

The barrier height is \(V_0 - E = 5.00\text{ eV} - 2.00\text{ eV} = 3.00\text{ eV}\). In Joules:

V_0 - E = 3.00 \times 1.60218 \times 10^{-19}\text{ J} = 4.80654 \times 10^{-19}\text{ J}

The decay wavevector \(\kappa\) is:

\kappa = \frac{\sqrt{2 m_e (V_0 - E)}}{\hbar} = \frac{\sqrt{2 (9.10938 \times 10^{-31}\text{ kg})(4.80654 \times 10^{-19}\text{ J})}}{1.05457 \times 10^{-34}\text{ J}\cdot\text{s}}

\kappa = \frac{\sqrt{8.7569 \times 10^{-49}}}{1.05457 \times 10^{-34}} = \frac{9.3578 \times 10^{-25}}{1.05457 \times 10^{-34}} \approx 8.8736 \times 10^9\text{ m}^{-1}

The penetration depth is:

\delta = \frac{1}{\kappa} = \frac{1}{8.8736 \times 10^9\text{ m}^{-1}} \approx 1.127 \times 10^{-10}\text{ m} = 0.1127\text{ nm} = 1.127\text{ Å}

Step 2: Transmission Probability for \(L = 0.200\text{ nm}\)

Evaluate \(\kappa L\):

\kappa L = (8.8736 \times 10^9\text{ m}^{-1})(0.200 \times 10^{-9}\text{ m}) = 1.7747

Evaluate \(\sinh(\kappa L)\):

\sinh(1.7747) = \frac{e^{1.7747} - e^{-1.7747}}{2} = \frac{5.8985 - 0.1695}{2} = 2.8645

\sinh^2(\kappa L) = (2.8645)^2 \approx 8.2054

Now evaluate the prefactor:

\frac{V_0^2}{4 E (V_0 - E)} = \frac{(5.00)^2}{4 (2.00)(3.00)} = \frac{25.0}{24.0} \approx 1.04167

The exact transmission coefficient is:

T = \frac{1}{1 + \frac{V_0^2}{4 E (V_0 - E)} \sinh^2(\kappa L)} = \frac{1}{1 + (1.04167)(8.2054)} = \frac{1}{1 + 8.5473} = \frac{1}{9.5473} \approx 0.10474\text{ (or } 10.47\%\text{)}

Step 3: Doubling the Barrier Width to \(L = 0.400\text{ nm}\)

For \(L = 0.400\text{ nm}\):

\kappa L = 2 \times 1.7747 = 3.5494

\sinh(3.5494) = \frac{e^{3.5494} - e^{-3.5494}}{2} \approx \frac{34.793 - 0.0287}{2} \approx 17.382

\sinh^2(\kappa L) = (17.382)^2 \approx 302.13

T_{0.4} = \frac{1}{1 + (1.04167)(302.13)} = \frac{1}{1 + 314.72} = \frac{1}{315.72} \approx 3.167 \times 10^{-3}\text{ (or } 0.317\%\text{)}

The transmission probability drops by a factor of:

\frac{T_{0.2}}{T_{0.4}} = \frac{0.10474}{0.003167} \approx 33.1\text{ times}

Doubling the barrier width attenuates the tunneling current by over 97%.

Advanced Example 2.7: Problem 7: Complete Orthonormality Matrix and Momentum Expectation Value

For a particle in a 1D box of length \(L\):

Prove by direct integration that the eigenfunctions \(\psi_m(x)\) and \(\psi_n(x)\) satisfy the orthonormality condition:

\int_0^L \psi_m^*(x) \psi_n(x) dx = \delta_{mn}

Calculate the expectation value of linear momentum \(\langle p \rangle\) in any arbitrary stationary state \(\psi_n(x)\).
Calculate the expectation value of kinetic energy \(\langle T \rangle\) and verify that \(\langle T \rangle = E_n\).

Comprehensive Multi-Step Solution:

Step 1: Proof of Orthonormality

The eigenfunctions are \(\psi_n(x) = \sqrt{\frac{2}{L}} \sin\left(\frac{n\pi x}{L}\right)\). Evaluate the inner product:

I_{mn} = \int_0^L \psi_m^*(x) \psi_n(x) dx = \frac{2}{L} \int_0^L \sin\left(\frac{m\pi x}{L}\right) \sin\left(\frac{n\pi x}{L}\right) dx

Using the product-to-sum identity \(\sin A \sin B = \frac{1}{2}[\cos(A - B) - \cos(A + B)]\):

I_{mn} = \frac{1}{L} \int_0^L \left[ \cos\left(\frac{(m - n)\pi x}{L}\right) - \cos\left(\frac{(m + n)\pi x}{L}\right) \right] dx

Case 1: \(m = n\):

I_{nn} = \frac{1}{L} \int_0^L \left[ \cos(0) - \cos\left(\frac{2n\pi x}{L}\right) \right] dx = \frac{1}{L} \left[ x - \frac{L}{2n\pi} \sin\left(\frac{2n\pi x}{L}\right) \right]_0^L = \frac{1}{L} [L - 0] = 1

Case 2: \(m \neq n\):

I_{mn} = \frac{1}{L} \left[ \frac{L}{(m - n)\pi} \sin\left(\frac{(m - n)\pi x}{L}\right) - \frac{L}{(m + n)\pi} \sin\left(\frac{(m + n)\pi x}{L}\right) \right]_0^L

Because \(\sin(k\pi) = 0\) for all integers \(k\), both sine terms evaluate to zero at both limits:

I_{mn} = 0 \quad (\text{for } m \neq n)

Combining both cases proves \(\langle \psi_m | \psi_n \rangle = \delta_{mn}\).

Step 2: Momentum Expectation Value \(\langle p \rangle\)

\langle p \rangle = \int_0^L \psi_n^*(x) \left(-i\hbar \frac{d}{dx}\right) \psi_n(x) dx

Substitute \(\psi_n(x)\):

\frac{d\psi_n}{dx} = \sqrt{\frac{2}{L}} \left(\frac{n\pi}{L}\right) \cos\left(\frac{n\pi x}{L}\right)

\langle p \rangle = -i\hbar \left(\frac{2}{L}\right) \left(\frac{n\pi}{L}\right) \int_0^L \sin\left(\frac{n\pi x}{L}\right) \cos\left(\frac{n\pi x}{L}\right) dx

Using \(\sin\theta \cos\theta = \frac{1}{2} \sin(2\theta)\):

\int_0^L \sin\left(\frac{2n\pi x}{L}\right) dx = \left[ -\frac{L}{2n\pi} \cos\left(\frac{2n\pi x}{L}\right) \right]_0^L = -\frac{L}{2n\pi} [\cos(2n\pi) - \cos(0)] = -\frac{L}{2n\pi} [1 - 1] = 0

Thus:

\langle p \rangle = 0

The average momentum is zero because the particle moves with equal probability in the \(+x\) and \(-x\) directions, forming a standing wave.

Step 3: Kinetic Energy Expectation Value

The kinetic energy operator is \(\hat{T} = -\frac{\hbar^2}{2m} \frac{d^2}{dx^2}\):

\frac{d^2\psi_n}{dx^2} = -\left(\frac{n\pi}{L}\right)^2 \psi_n(x)

Therefore:

\hat{T} \psi_n(x) = \frac{\hbar^2}{2m} \left(\frac{n\pi}{L}\right)^2 \psi_n(x) = E_n \psi_n(x)

Evaluating the expectation value:

\langle T \rangle = \int_0^L \psi_n^* \hat{T} \psi_n dx = E_n \int_0^L |\psi_n|^2 dx = E_n \cdot 1 = E_n

Since \(V(x) = 0\) inside the box, the total energy is purely kinetic: \(E_n = \langle T \rangle\).

Advanced Example 2.8: Inversion Tunneling Splitting and Inversion Time in Ammonia

In the ammonia molecule (\(\text{NH}_3\)), the nitrogen atom tunnels through the plane of three hydrogen atoms between two equivalent pyramidal configurations. The tunnel splitting between the symmetric ground state and antisymmetric first excited state is measured to be \(\Delta\tilde{\nu} = 0.7935\text{ cm}^{-1}\).

Calculate the energy splitting \(\Delta E\) in Joules and in electron volts.
Determine the microwave inversion tunneling frequency \(\nu_{\text{inv}}\) in GHz.
If the nitrogen atom is prepared initially localized in the left pyramidal configuration at \(t = 0\), calculate the time required for complete inversion to the right configuration.

Comprehensive Multi-Step Solution:

Step 1: Energy Splitting \(\Delta E\)

Given \(\Delta\tilde{\nu} = 0.7935\text{ cm}^{-1}\):

\Delta E = h c \Delta\tilde{\nu}

Using \(h = 6.62607 \times 10^{-34}\text{ J}\cdot\text{s}\) and \(c = 2.99792 \times 10^{10}\text{ cm/s}\):

h c = (6.62607 \times 10^{-34})(2.99792 \times 10^{10}) \approx 1.98645 \times 10^{-23}\text{ J}\cdot\text{cm}

\Delta E = (1.98645 \times 10^{-23}\text{ J}\cdot\text{cm}) \times 0.7935\text{ cm}^{-1} \approx 1.5762 \times 10^{-23}\text{ J}

Converting to electron volts (\(1\text{ eV} = 1.60218 \times 10^{-19}\text{ J}\)):

\Delta E = \frac{1.5762 \times 10^{-23}\text{ J}}{1.60218 \times 10^{-19}\text{ J/eV}} \approx 9.838 \times 10^{-5}\text{ eV} = 98.38\text{ \mu eV}

Step 2: Inversion Frequency \(\nu_{\text{inv}}\)

\nu_{\text{inv}} = c \Delta\tilde{\nu} = (2.99792 \times 10^{10}\text{ cm/s}) \times 0.7935\text{ cm}^{-1} \approx 2.3789 \times 10^{10}\text{ Hz} = 23.789\text{ GHz}

The inversion transition occurs at \(23.79\text{ GHz}\) in the microwave K-band (\(\lambda \approx 1.26\text{ cm}\)).

Step 3: Complete Inversion Time \(\tau_{\text{inv}}\)

The wavefunction evolves as:

P_R(t) = \sin^2\left( \frac{\Delta E \, t}{2\hbar} \right) = \sin^2(\pi \nu_{\text{inv}} t)

Complete inversion (\(P_R = 1\)) occurs when the phase reaches \(\pi/2\):

\pi \nu_{\text{inv}} \tau_{\text{inv}} = \frac{\pi}{2} \implies \tau_{\text{inv}} = \frac{1}{2 \nu_{\text{inv}}}

Substitute \(\nu_{\text{inv}} = 2.3789 \times 10^{10}\text{ s}^{-1}\):

\tau_{\text{inv}} = \frac{1}{2 \times 2.3789 \times 10^{10}\text{ s}^{-1}} \approx 2.102 \times 10^{-11}\text{ s} = 21.02\text{ picoseconds}

The nitrogen atom tunnels back and forth through the barrier approximately 24 billion times every second, inverting every 21 picoseconds.

Advanced Example 2.9: Two-Dimensional Circular Quantum Billiard and Bessel Functions

Consider a particle of mass \(m\) confined within a 2D circular billiard (quantum disk) of radius \(R\) with hard walls: \(V(r) = 0\) for \(r < R\) and \(\infty\) for \(r \ge R\).

Set up the Schrödinger equation in polar coordinates \((r, \phi)\) and separate variables \(\psi(r, \phi) = R(r) \Phi(\phi)\).
Show that the radial equation is Bessel's differential equation and state the boundary conditions.
Express the quantized energy eigenvalues in terms of the zeros of Bessel functions \(\alpha_{m, n}\), and calculate the ground-state energy for an electron in a quantum corral of radius \(R = 5.0\text{ nm}\).

Comprehensive Multi-Step Solution:

Step 1: Separation of Variables in Polar Coordinates

In polar coordinates \((r, \phi)\), the Laplacian is:

\nabla^2 = \frac{\partial^2}{\partial r^2} + \frac{1}{r} \frac{\partial}{\partial r} + \frac{1}{r^2} \frac{\partial^2}{\partial \phi^2}

The time-independent Schrödinger equation for \(r < R\) is:

-\frac{\hbar^2}{2m} \left( \frac{\partial^2 \psi}{\partial r^2} + \frac{1}{r} \frac{\partial \psi}{\partial r} + \frac{1}{r^2} \frac{\partial^2 \psi}{\partial \phi^2} \right) = E \psi

Let \(k = \sqrt{2m E}/\hbar\) and substitute \(\psi(r, \phi) = R(r) \Phi(\phi)\):

\frac{r^2}{R(r)} \left( \frac{d^2 R}{dr^2} + \frac{1}{r} \frac{dR}{dr} \right) + k^2 r^2 = -\frac{1}{\Phi(\phi)} \frac{d^2 \Phi}{d\phi^2} = m_l^2

The angular equation \(\frac{d^2 \Phi}{d\phi^2} = -m_l^2 \Phi\) with cyclic boundary condition \(\Phi(\phi + 2\pi) = \Phi(\phi)\) requires:

m_l \in \{0, \pm 1, \pm 2, \dots\}, \quad \Phi(\phi) = \frac{1}{\sqrt{2\pi}} e^{i m_l \phi}

Step 2: Radial Equation and Bessel Functions

The radial equation becomes:

r^2 \frac{d^2 R}{dr^2} + r \frac{dR}{dr} + (k^2 r^2 - m_l^2) R = 0

Let \(\rho = k r\):

\rho^2 \frac{d^2 R}{d\rho^2} + \rho \frac{dR}{d\rho} + (\rho^2 - m_l^2) R = 0

This is Bessel's differential equation of order \(m_l\). The general solution is a linear combination of Bessel functions of the first kind \(J_{m_l}(\rho)\) and Neumann functions (Bessel functions of the second kind) \(Y_{m_l}(\rho)\):

R(r) = C_1 J_{m_l}(k r) + C_2 Y_{m_l}(k r)

Because \(Y_{m_l}(k r) \rightarrow -\infty\) as \(r \rightarrow 0\), physical wavefunctions must have \(C_2 = 0\):

R(r) = C J_{m_l}(k r)

Boundary condition at the wall \(r = R\):

R(R) = 0 \implies J_{m_l}(k R) = 0

Step 3: Quantized Energy Spectrum and Numerical Calculation

Let \(\alpha_{m_l, n}\) denote the \(n\)-th positive root (zero) of the Bessel function \(J_{m_l}(x)\):

k R = \alpha_{m_l, n} \implies k = \frac{\alpha_{m_l, n}}{R}

The quantized energy eigenvalues are:

E_{m_l, n} = \frac{\hbar^2 k^2}{2m} = \frac{\hbar^2 \alpha_{m_l, n}^2}{2 m R^2} \quad (m_l = 0, 1, 2, \dots; \; n = 1, 2, 3, \dots)

The absolute ground state occurs for \(m_l = 0\) (zero angular momentum) and the first root \(n = 1\). The first zero of \(J_0(x)\) is:

\alpha_{0, 1} \approx 2.4048

For an electron (\(m_e = 9.10938 \times 10^{-31}\text{ kg}\)) in a corral of radius \(R = 5.0\text{ nm} = 5.0 \times 10^{-9}\text{ m}\):

E_{0, 1} = \frac{(1.05457 \times 10^{-34}\text{ J}\cdot\text{s})^2 (2.4048)^2}{2(9.10938 \times 10^{-31}\text{ kg})(5.0 \times 10^{-9}\text{ m})^2} = \frac{(1.1121 \times 10^{-68})(5.7832)}{(1.8219 \times 10^{-30})(2.5 \times 10^{-17})}

E_{0, 1} = \frac{6.4315 \times 10^{-68}}{4.5547 \times 10^{-47}} \approx 1.412 \times 10^{-21}\text{ J}

Converting to electron volts:

E_{0, 1} = \frac{1.412 \times 10^{-21}\text{ J}}{1.60218 \times 10^{-19}\text{ J/eV}} \approx 0.00881\text{ eV} = 8.81\text{ meV}

This describes the standing de Broglie electron wave patterns observed with scanning tunneling microscopes (STM) in artificial atomic quantum corrals.