Chapter 1 • Theory & Derivations

Unit 1: Foundations & Historical Development of Quantum Mechanics

Comprehensive foundational treatment of the breakdown of classical Newtonian and Maxwellian physics at microscopic scales: blackbody spectral radiance and ultraviolet divergence, Planck's quantized energy hypothesis, the photoelectric effect and Einstein's photon momentum, Compton scattering kinematics, de Broglie matter waves, the Heisenberg uncertainty principle, Schrödinger wave mechanics, Hermitian operator algebra, and the foundational postulates of quantum theory.

1.1Failure of Classical Mechanics & The Ultraviolet Catastrophe

By the late nineteenth century, classical mechanics (Newtonian dynamics) and classical electrodynamics (Maxwell's field equations) successfully described macroscopic phenomena. However, applying classical equipartition of energy to electromagnetic cavity radiation led to a profound mathematical catastrophe.

The Classical Rayleigh-Jeans Radiation Law

Consider an isothermal cavity (hohlraum) at absolute temperature \(T\) containing standing electromagnetic waves with perfectly conducting walls of volume \(V\).

In classical physics, each spatial standing wave mode possesses two independent transverse polarization states. The number of electromagnetic cavity modes per unit volume in the frequency interval \([\nu, \nu + d\nu]\) is given by the three-dimensional density of states:

g(\nu) d\nu = \frac{8\pi \nu^2}{c^3} d\nu

According to the classical equipartition theorem of statistical mechanics, every quadratic degree of freedom in thermal equilibrium carries an average thermal energy of \(\frac{1}{2} k_B T\). Because each standing electromagnetic mode corresponds to an equivalent one-dimensional harmonic oscillator (possessing kinetic and potential energy quadratic terms), the classical average energy per mode is:

\langle \epsilon \rangle_{\text{class}} = k_B T

Multiplying the mode density by the classical average energy yields the Rayleigh-Jeans spectral energy density:

\rho_{\text{RJ}}(\nu, T) d\nu = g(\nu) \langle \epsilon \rangle_{\text{class}} d\nu = \frac{8\pi \nu^2}{c^3} k_B T d\nu

Expressed in terms of wavelength \(\lambda = c / \nu\):

\rho_{\text{RJ}}(\lambda, T) d\lambda = \frac{8\pi k_B T}{\lambda^4} d\lambda

The Ultraviolet Catastrophe

While the Rayleigh-Jeans law agreed with infrared measurements at long wavelengths (\(\lambda \rightarrow \infty\)), it diverged catastrophically as \(\lambda \rightarrow 0\) (\(\nu \rightarrow \infty\)):

\lim_{\nu \rightarrow \infty} \rho_{\text{RJ}}(\nu, T) = \infty

Integrating the total radiant energy density across all frequencies yielded an impossible physical infinity:

U_{\text{total}} = \int_0^\infty \rho_{\text{RJ}}(\nu, T) d\nu = \frac{8\pi k_B T}{c^3} \int_0^\infty \nu^2 d\nu = \infty

This divergence, termed the ultraviolet catastrophe by Paul Ehrenfest, meant that any heated cavity should instantaneously radiate infinite energy into high-frequency ultraviolet, X-ray, and gamma-ray modes—a physical impossibility.

Planck's Revolutionary Quantization Hypothesis

In December 1900, Max Planck resolved the catastrophe by proposing that the atomic resonators constituting the cavity walls cannot absorb or emit energy continuously. Instead, energy exchange occurs strictly in discrete, indivisible packets called quanta:

\epsilon_n = n h \nu, \quad n \in \{0, 1, 2, 3, \dots\}

where \(h = 6.62607015 \times 10^{-34}\text{ J}\cdot\text{s}\) is Planck's constant.

Applying Maxwell-Boltzmann statistics to the discrete energy states, the quantum statistical average energy of a cavity oscillator is:

\langle \epsilon \rangle_{\text{quant}} = \frac{\sum_{n=0}^\infty (n h \nu) e^{-n h \nu / k_B T}}{\sum_{n=0}^\infty e^{-n h \nu / k_B T}}

Letting \(x = e^{-h\nu / k_B T}\):

\sum_{n=0}^\infty x^n = \frac{1}{1 - x}, \quad \sum_{n=0}^\infty n x^n = x \frac{d}{dx} \left( \frac{1}{1 - x} \right) = \frac{x}{(1 - x)^2}

Hence:

\langle \epsilon \rangle_{\text{quant}} = h\nu \frac{x}{1 - x} = \frac{h\nu}{e^{h\nu / k_B T} - 1}

Multiplying by the mode density \(g(\nu)\) yields Planck's radiation law:

\rho(\nu, T) d\nu = \frac{8\pi h \nu^3}{c^3} \frac{1}{e^{h\nu / k_B T} - 1} d\nu

In wavelength coordinates:

\rho(\lambda, T) d\lambda = \frac{8\pi h c}{\lambda^5} \frac{1}{e^{h c / \lambda k_B T} - 1} d\lambda

Asymptotic Limits and Wien's Displacement Law

**Low Frequency / Long Wavelength Limit (\(h\nu \ll k_B T\)):**

Expanding the exponential \(e^{h\nu / k_B T} \approx 1 + \frac{h\nu}{k_B T}\):

\langle \epsilon \rangle \approx \frac{h\nu}{1 + \frac{h\nu}{k_B T} - 1} = k_B T

reproducing the classical Rayleigh-Jeans expression.

**High Frequency / Short Wavelength Limit (\(h\nu \gg k_B T\)):**

\(e^{h\nu / k_B T} \gg 1\):

\rho(\nu, T) \approx \frac{8\pi h \nu^3}{c^3} e^{-h\nu / k_B T}

reproducing Wien's empirical exponential distribution and preventing the ultraviolet catastrophe.

**Wien's Displacement Law:**

Differentiating \(\rho(\lambda, T)\) with respect to \(\lambda\) and setting \(\frac{\partial \rho}{\partial \lambda} = 0\) leads to the transcendental equation \(5(1 - e^{-y}) - y = 0\) where \(y = \frac{h c}{\lambda_{\text{max}} k_B T}\). Numerical solution gives \(y \approx 4.965114\), yielding:

\lambda_{\text{max}} T = \frac{h c}{4.965114 k_B} = b \approx 2.89777 \times 10^{-3}\text{ m}\cdot\text{K}

Planck Blackbody Radiance & Compton Wavepacket Dispersion Unit 1: Quantum Foundations

Temperature T: 5500 K

Packet Width σ₀: 1.2 Å

Wavenumber k₀: 1.5 Å⁻¹

1.2Photoelectric & Compton Effects: Particle Nature of Light

While Planck viewed quantization as a mathematical property of atomic cavity resonators, Albert Einstein (1905) established that electromagnetic radiation itself propagates as localized, corpuscular quanta—photons—each carrying discrete energy and momentum.

The Photoelectric Effect

When monochromatic ultraviolet or visible light strikes a clean metallic surface, electrons (photoelectrons) are ejected. Classical wave theory predicted that:

Electron kinetic energy should increase with increasing light wave intensity (electric field amplitude squared).
Electron emission should occur after a measurable time lag at low light intensities while the continuous wavefront accumulates sufficient energy.
Photoelectrons should be ejected at any frequency provided the intensity is sufficiently large.

Experimental observations by Heinrich Hertz, Philipp Lenard, and Robert Millikan directly contradicted classical theory:

Electron emission occurs instantaneously (time lag \(< 10^{-9}\text{ s}\)) even under ultra-weak illumination.
A well-defined threshold frequency \(\nu_0\) exists below which zero photoelectrons are emitted regardless of beam intensity.
The maximum kinetic energy \(K_{\text{max}}\) depends linearly on light frequency \(\nu\) and is completely independent of light intensity.
Increasing light intensity increases the photocurrent (number of ejected electrons per second) without changing \(K_{\text{max}}\).

Einstein's Photoelectric Equation

Einstein proposed that a single photon of energy \(h\nu\) collides with a single bound electron in the metal. Overcoming the surface electrostatic potential barrier requires a characteristic binding energy termed the work function \(\Phi\). The conservation of energy yields:

h\nu = \Phi + K_{\text{max}} = \Phi + \frac{1}{2} m_e v_{\text{max}}^2

Expressing \(K_{\text{max}}\) in terms of the stopping potential \(V_s\) required to reduce the photocurrent to zero (\(K_{\text{max}} = e V_s\)):

e V_s = h\nu - \Phi \implies V_s = \left(\frac{h}{e}\right)\nu - \frac{\Phi}{e}

The slope of \(V_s\) versus frequency \(\nu\) is universally equal to \(h/e\), providing an independent measurement of Planck's constant.

The Compton Effect

In 1923, Arthur H. Compton directed monochromatic X-rays (\(\lambda \sim 0.07\text{ nm}\)) at a graphite target and observed that the scattered radiation contained both the incident wavelength \(\lambda\) and an unexpected shifted, longer wavelength \(\lambda'\).

Treating the collision between an X-ray photon and an initially stationary free electron (\(m_0\)) using relativistic mechanics:

**Incident Photon:** Energy \(E = h\nu\), momentum \(\mathbf{p} = \frac{h\nu}{c} \hat{\mathbf{i}}\).
**Target Electron (at rest):** Energy \(E_e = m_0 c^2\), momentum \(\mathbf{p}_e = 0\).
**Scattered Photon (at angle \(\theta\)):** Energy \(E' = h\nu'\), momentum \(\mathbf{p}' = \frac{h\nu'}{c}\).
**Recoil Electron (at angle \(\phi\)):** Relativistic energy \(E_e' = \sqrt{p_e^2 c^2 + m_0^2 c^4}\), momentum \(\mathbf{p}_e\).

From conservation of relativistic momentum:

\mathbf{p} = \mathbf{p}' + \mathbf{p}_e \implies \mathbf{p}_e = \mathbf{p} - \mathbf{p}'

Squaring both sides:

p_e^2 = p^2 + p'^2 - 2 p p' \cos\theta = \left(\frac{h\nu}{c}\right)^2 + \left(\frac{h\nu'}{c}\right)^2 - 2 \left(\frac{h\nu}{c}\right)\left(\frac{h\nu'}{c}\right) \cos\theta

From conservation of total relativistic energy:

h\nu + m_0 c^2 = h\nu' + E_e' \implies E_e' = h(\nu - \nu') + m_0 c^2

Equating \((E_e')^2 = p_e^2 c^2 + m_0^2 c^4\):

[h(\nu - \nu') + m_0 c^2]^2 = p_e^2 c^2 + m_0^2 c^4

Expanding and substituting \(p_e^2\):

h^2(\nu - \nu')^2 + 2 h m_0 c^2 (\nu - \nu') + m_0^2 c^4 = h^2 \nu^2 + h^2 \nu'^2 - 2 h^2 \nu \nu' \cos\theta + m_0^2 c^4

Subtracting common terms yields:

-2 h^2 \nu \nu' + 2 h m_0 c^2 (\nu - \nu') = -2 h^2 \nu \nu' \cos\theta

Dividing throughout by \(2 h m_0 c^2 \nu \nu'\):

\frac{1}{\nu'} - \frac{1}{\nu} = \frac{h}{m_0 c^2} (1 - \cos\theta)

Multiplying by \(c\) (\(\lambda = c/\nu\)) yields the Compton shift equation:

\Delta\lambda = \lambda' - \lambda = \frac{h}{m_0 c} (1 - \cos\theta) = \lambda_C (1 - \cos\theta)

where \(\lambda_C = \frac{h}{m_0 c} \approx 2.42631 \times 10^{-12}\text{ m} = 0.02426\text{ Å}\) is the Compton wavelength of the electron.

1.3de Broglie Hypothesis & Matter Wave Diffraction

In 1924, Louis de Broglie extended the wave-particle duality of light to all material particles. If electromagnetic waves possess particle-like characteristics (photons with \(p = h/\lambda\)), then material particles (electrons, protons, neutrons) must possess intrinsic wave-like properties.

The de Broglie Relations

For any particle with relativistic energy \(E\) and linear momentum \(p = m v\), de Broglie associated a matter wave of frequency \(\nu\) and wavelength \(\lambda\):

\lambda = \frac{h}{p} = \frac{h}{m v}

\nu = \frac{E}{h}

Expressed in terms of the reduced Planck constant \(\hbar = \frac{h}{2\pi}\), wavevector \(k = \frac{2\pi}{\lambda}\), and angular frequency \(\omega = 2\pi\nu\):

\mathbf{p} = \hbar \mathbf{k}, \quad E = \hbar \omega

Thermal and Non-Relativistic Kinetic Scaling

For a non-relativistic particle of mass \(m\) accelerated through an electrostatic potential difference \(V\):

K = \frac{p^2}{2m} = e V \implies p = \sqrt{2 m e V}

The corresponding de Broglie wavelength is:

\lambda = \frac{h}{\sqrt{2 m e V}}

For an electron (\(m_e = 9.109 \times 10^{-31}\text{ kg}\), \(e = 1.602 \times 10^{-19}\text{ C}\)):

\lambda_e = \frac{1.226}{\sqrt{V}}\text{ nm} = \sqrt{\frac{150}{V}}\text{ Å}

For an accelerating voltage \(V = 100\text{ V}\), \(\lambda_e \approx 0.123\text{ nm}\) (1.23 Å), which is of the exact order of interatomic lattice spacings in crystalline solids (\(d \sim 1 - 3\text{ Å}\)).

Experimental Verification: Davisson-Germer & Thomson Experiments

In 1927, Clinton Davisson and Lester Germer scattered slow electrons (\(54\text{ eV}\)) from a target nickel single crystal. They observed a pronounced peak in scattered electron intensity at an angle \(\phi = 50^\circ\).

Using Bragg's diffraction law for atomic lattice planes with spacing \(D = 0.215\text{ nm}\):

n \lambda = 2 d \sin\theta = 2 (D \sin(\phi/2)) \sin\theta = D \sin\phi = 0.215\text{ nm} \times \sin(50^\circ) = 0.165\text{ nm}

The theoretical de Broglie wavelength for a \(54\text{ eV}\) electron is:

\lambda = \frac{h}{\sqrt{2 m_e (54\text{ eV})}} = 0.167\text{ nm}

The extraordinary agreement (\(0.165\text{ nm}\) vs \(0.167\text{ nm}\)) definitively proved that electrons travel as coherent wave fields capable of constructive and destructive interference.

1.4Heisenberg Uncertainty Principle & Wavepacket Dispersion

The wave nature of matter implies that a microscopic particle cannot be localized to an infinitesimal geometric point while possessing a sharply defined momentum. In 1927, Werner Heisenberg formulated this fundamental limitation.

Mathematical Origin from Fourier Transform Pairs

A localized quantum particle is mathematically represented as a wavepacket formed by the continuous superposition of monochromatic plane waves:

\psi(x, 0) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^\infty \phi(k) e^{i k x} dk

where \(\phi(k)\) is the momentum-space probability amplitude obtained via inverse Fourier transform:

\phi(k) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^\infty \psi(x, 0) e^{-i k x} dx

For a standard Gaussian wavepacket:

\psi(x, 0) = \left( \frac{1}{2\pi \sigma_x^2} \right)^{1/4} \exp\left( -\frac{x^2}{4\sigma_x^2} + i k_0 x \right)

Its Fourier transform in wavenumber space is also Gaussian:

\phi(k) = \left( \frac{2\sigma_x^2}{\pi} \right)^{1/4} \exp\left( -\sigma_x^2 (k - k_0)^2 \right)

The standard deviation in spatial position is \(\Delta x = \sigma_x\), and the standard deviation in wavenumber is \(\Delta k = \frac{1}{2\sigma_x}\).

Using de Broglie's relation \(p = \hbar k \implies \Delta p = \hbar \Delta k\):

\Delta x \cdot \Delta p = \sigma_x \cdot \left( \frac{\hbar}{2\sigma_x} \right) = \frac{\hbar}{2}

For any arbitrary (non-Gaussian) normalized wavefunction, the product of standard deviations satisfies the general Heisenberg uncertainty inequality:

\Delta x \cdot \Delta p_x \ge \frac{\hbar}{2}

Generalized Robertson-Schrödinger Uncertainty Relation

In operator quantum mechanics, for any two Hermitian operators \(\hat{A}\) and \(\hat{B}\):

\Delta A = \sqrt{\langle \hat{A}^2 \rangle - \langle \hat{A} \rangle^2}, \quad \Delta B = \sqrt{\langle \hat{B}^2 \rangle - \langle \hat{B} \rangle^2}

The general uncertainty product is bounded by their commutator \([\hat{A}, \hat{B}] = \hat{A}\hat{B} - \hat{B}\hat{A}\):

\Delta A \cdot \Delta B \ge \frac{1}{2} |\langle [\hat{A}, \hat{B}] \rangle|

For position \(\hat{x} = x\) and momentum \(\hat{p}_x = -i\hbar \frac{\partial}{\partial x}\):

[\hat{x}, \hat{p}_x] f(x) = x \left(-i\hbar \frac{\partial f}{\partial x}\right) - \left(-i\hbar \frac{\partial (x f)}{\partial x}\right) = i\hbar f(x) \implies [\hat{x}, \hat{p}_x] = i\hbar \hat{I}

Substituting into Robertson's theorem yields \(\Delta x \cdot \Delta p_x \ge \frac{1}{2} |i\hbar| = \frac{\hbar}{2}\).

Time-Energy Uncertainty Relation

A complementary uncertainty relation connects the lifetime \(\Delta t\) of an excited quantum state to the uncertainty in its energy \(\Delta E\):

\Delta E \cdot \Delta t \ge \frac{\hbar}{2}

If an excited atomic or molecular state has a finite radiative lifetime \(\tau\), the emitted photon energy cannot be monochromatic. It possesses an intrinsic natural linewidth \(\Delta \nu = \frac{1}{2\pi \tau}\).

1.5Schrödinger Wave Equation: Time-Dependent & Stationary States

In 1926, Erwin Schrödinger developed wave mechanics, replacing classical trajectories \(\mathbf{r}(t)\) with a continuous, complex-valued probability amplitude wavefunction \(\Psi(\mathbf{r}, t)\).

The Time-Dependent Schrödinger Equation (TDSE)

For a single non-relativistic particle of mass \(m\) moving in an external scalar potential field \(V(\mathbf{r}, t)\), the state vector evolves according to:

i\hbar \frac{\partial \Psi(\mathbf{r}, t)}{\partial t} = \hat{H} \Psi(\mathbf{r}, t)

where \(\hat{H}\) is the Hamiltonian operator:

\hat{H} = \hat{T} + \hat{V} = -\frac{\hbar^2}{2m} \nabla^2 + V(\mathbf{r}, t)

In one Cartesian dimension:

i\hbar \frac{\partial \Psi(x, t)}{\partial t} = -\frac{\hbar^2}{2m} \frac{\partial^2 \Psi(x, t)}{\partial x^2} + V(x, t) \Psi(x, t)

Separation of Variables & The Time-Independent Schrödinger Equation

When the potential energy is static (\(V(\mathbf{r}, t) = V(\mathbf{r})\)), we apply the separation of variables ansatz:

\Psi(\mathbf{r}, t) = \psi(\mathbf{r}) \phi(t)

Substituting into the TDSE:

i\hbar \psi(\mathbf{r}) \frac{d\phi(t)}{dt} = \left[ -\frac{\hbar^2}{2m} \nabla^2 \psi(\mathbf{r}) + V(\mathbf{r}) \psi(\mathbf{r}) \right] \phi(t)

Dividing throughout by \(\psi(\mathbf{r}) \phi(t)\):

\frac{i\hbar}{\phi(t)} \frac{d\phi(t)}{dt} = \frac{1}{\psi(\mathbf{r})} \left[ -\frac{\hbar^2}{2m} \nabla^2 \psi(\mathbf{r}) + V(\mathbf{r}) \psi(\mathbf{r}) \right] = E

Because the left side depends purely on time \(t\) while the right side depends purely on position \(\mathbf{r}\), both sides must equal a spatial-temporal constant \(E\) (the total energy).

Solving the temporal differential equation:

\frac{d\phi}{dt} = -\frac{i E}{\hbar} \phi \implies \phi(t) = \exp\left(-\frac{i E t}{\hbar}\right)

The spatial component satisfies the Time-Independent Schrödinger Equation (TISE):

\hat{H} \psi(\mathbf{r}) = E \psi(\mathbf{r})

\left[ -\frac{\hbar^2}{2m} \nabla^2 + V(\mathbf{r}) \right] \psi(\mathbf{r}) = E \psi(\mathbf{r})

Properties of Stationary States

A state described by a single energy eigenfunction \(\psi_n(\mathbf{r})\) with eigenvalue \(E_n\) has total wavefunction:

\Psi_n(\mathbf{r}, t) = \psi_n(\mathbf{r}) e^{-i E_n t / \hbar}

**Time-Invariant Probability Density:**

|\Psi_n(\mathbf{r}, t)|^2 = \Psi_n^(\mathbf{r}, t) \Psi_n(\mathbf{r}, t) = \psi_n^(\mathbf{r}) e^{+i E_n t / \hbar} \psi_n(\mathbf{r}) e^{-i E_n t / \hbar} = |\psi_n(\mathbf{r})|^2

The probability distribution is strictly static over time (hence stationary state).

**Stationary Expectation Values:**

For any time-independent operator \(\hat{A}\):

\langle \hat{A} \rangle(t) = \int \Psi_n^ \hat{A} \Psi_n d\tau = \int \psi_n^ \hat{A} \psi_n d\tau = \text{constant}

**Superposition and Non-Stationary Dynamics:**

If a system is prepared in a linear combination of two or more distinct stationary states:

\Psi(\mathbf{r}, t) = c_1 \psi_1(\mathbf{r}) e^{-i E_1 t / \hbar} + c_2 \psi_2(\mathbf{r}) e^{-i E_2 t / \hbar}

the probability density beats at the Bohr transition frequency \(\omega_{21} = \frac{E_2 - E_1}{\hbar}\):

|\Psi(\mathbf{r}, t)|^2 = |c_1|^2 |\psi_1|^2 + |c_2|^2 |\psi_2|^2 + 2 \text{Re}\left[ c_1^ c_2 \psi_1^ \psi_2 e^{-i (E_2 - E_1) t / \hbar} \right]

1.6Operators, Eigenvalues, Hermitian Adjoints & Commutation Rules

In quantum mechanics, physical observables are represented by linear operators acting on vectors in a complex Hilbert space \(\mathcal{H}\).

Operator Algebra and Linearity

An operator \(\hat{A}\) is linear if for all wavefunctions \(\psi_1, \psi_2\) and complex scalars \(c_1, c_2\):

\hat{A} (c_1 \psi_1 + c_2 \psi_2) = c_1 \hat{A}\psi_1 + c_2 \hat{A}\psi_2

The eigenvalue equation for an operator \(\hat{A}\) is:

\hat{A} \psi_n = a_n \psi_n

where \(\psi_n\) is the eigenfunction and \(a_n\) is the corresponding scalar eigenvalue.

Hermitian (Self-Adjoint) Operators

Because physical measurements yield real numbers, all physical observables must correspond to Hermitian operators.

The Hermitian adjoint (conjugate transpose) \(\hat{A}^\dagger\) of an operator \(\hat{A}\) is defined by the inner product relation:

\langle \phi | \hat{A} \psi \rangle = \int \phi^ (\hat{A} \psi) d\tau = \int (\hat{A}^\dagger \phi)^ \psi d\tau = \langle \hat{A}^\dagger \phi | \psi \rangle

An operator is Hermitian if \(\hat{A}^\dagger = \hat{A}\):

\int \phi^ (\hat{A} \psi) d\tau = \int (\hat{A} \phi)^ \psi d\tau

Theorem 1: Real Eigenvalues

Let \(\hat{A} \psi = a \psi\). Taking the inner product with \(\psi\):

\langle \psi | \hat{A} \psi \rangle = a \langle \psi | \psi \rangle

Since \(\hat{A}\) is Hermitian:

\langle \psi | \hat{A} \psi \rangle = \langle \hat{A} \psi | \psi \rangle = a^* \langle \psi | \psi \rangle

Therefore, \((a - a^*) \langle \psi | \psi \rangle = 0\). Since \(\langle \psi | \psi \rangle > 0\):

a = a^* \implies a \in \mathbb{R}

Theorem 2: Orthogonality of Eigenfunctions

Let \(\hat{A} \psi_1 = a_1 \psi_1\) and \(\hat{A} \psi_2 = a_2 \psi_2\) with distinct eigenvalues \(a_1 \neq a_2\):

\langle \psi_1 | \hat{A} \psi_2 \rangle = a_2 \langle \psi_1 | \psi_2 \rangle

\langle \hat{A} \psi_1 | \psi_2 \rangle = a_1^* \langle \psi_1 | \psi_2 \rangle = a_1 \langle \psi_1 | \psi_2 \rangle

Subtracting the two equations:

(a_2 - a_1) \langle \psi_1 | \psi_2 \rangle = 0

Since \(a_1 \neq a_2\):

\langle \psi_1 | \psi_2 \rangle = \int \psi_1^* \psi_2 d\tau = 0

Eigenfunctions belonging to distinct eigenvalues of a Hermitian operator are strictly orthogonal.

Commutators and Simultaneous Observables

The commutator of two operators \(\hat{A}\) and \(\hat{B}\) is defined as:

[\hat{A}, \hat{B}] = \hat{A}\hat{B} - \hat{B}\hat{A}

If \([\hat{A}, \hat{B}] = 0\), the operators **commute**. There exists a complete set of mutual simultaneous eigenfunctions \(\psi_{a,b}\) such that:

\hat{A} \psi_{a,b} = a \psi_{a,b}, \quad \hat{B} \psi_{a,b} = b \psi_{a,b}

Both observables can be measured simultaneously to arbitrary precision (\(\Delta A \cdot \Delta B = 0\)).

If \([\hat{A}, \hat{B}] \neq 0\), the operators are **incompatible**. Measuring one observable inherently disrupts the state and creates fundamental uncertainty in the other.

1.7Postulates of Quantum Mechanics & Dirac Bracket Formalism

The axiomatic mathematical foundation of modern non-relativistic quantum theory is codified in Dirac's bra-ket formalism.

The Six Foundational Postulates

Postulate 1: The State Space

The physical state of a quantum system is completely specified by a state vector (ket) \(|\Psi(t)\rangle\) residing in a complex Hilbert space \(\mathcal{H}\). The state is normalized such that:

\langle \Psi | \Psi \rangle = \int |\Psi(\mathbf{r}, t)|^2 d\tau = 1

Max Born's statistical interpretation dictates that \(P(\mathbf{r}) d\tau = |\Psi(\mathbf{r}, t)|^2 d\tau\) represents the probability of finding the particle in differential volume \(d\tau\).

Postulate 2: Physical Observables

To every physically measurable dynamical observable \(A\) in classical mechanics, there corresponds a linear Hermitian operator \(\hat{A}\) acting in Hilbert space. Position and momentum operators obey the canonical commutation relation:

[\hat{x}_j, \hat{p}_k] = i\hbar \delta_{jk} \hat{I}

Postulate 3: Measurement Eigenvalues

The only possible outcomes of a single precise measurement of an observable \(A\) are the eigenvalues \(a_n\) of its associated Hermitian operator \(\hat{A}\):

\hat{A} |a_n\rangle = a_n |a_n\rangle

Postulate 4: Measurement Probabilities & State Collapse

Any normalized state \(|\Psi\rangle\) can be expanded in terms of the complete orthonormal eigenbasis \(\{|a_n\rangle\}\) of \(\hat{A}\):

|\Psi\rangle = \sum_n c_n |a_n\rangle, \quad c_n = \langle a_n | \Psi\rangle

The probability \(P(a_n)\) of measuring eigenvalue \(a_n\) is:

P(a_n) = |\langle a_n | \Psi \rangle|^2 = |c_n|^2

Immediately upon measuring the value \(a_n\), the state vector instantaneously collapses to the corresponding eigenstate \(|a_n\rangle\).

Postulate 5: Time Evolution

Between non-demolition measurements, the time evolution of the state vector is governed by the time-dependent Schrödinger equation:

i\hbar \frac{d}{dt} |\Psi(t)\rangle = \hat{H} |\Psi(t)\rangle

Postulate 6: Pauli Antisymmetry (Symmetrization Postulate)

The total state vector of a system of identical indistinguishable particles must be symmetric under particle permutation for bosons (integer spin \(S = 0, 1, 2\)) and antisymmetric for fermions (half-integer spin \(S = 1/2, 3/2\)):

\hat{P}_{12} |\Psi(1, 2)\rangle = (-1)^{2S} |\Psi(1, 2)\rangle

1.8Postulates of Quantum Mechanics & Dirac's Bra-Ket Formalism

The axiomatic mathematical foundation of quantum mechanics is formalized through six rigorous postulates, unified by Paul Dirac's elegant bra-ket notation in Hilbert space \(\mathcal{H}\).

Postulate 1: The State of a Quantum System

The physical state of a microscopic system at time \(t\) is completely described by a state vector (ket) \(|\Psi(t)\rangle\) residing in a complex Hilbert space \(\mathcal{H}\). The inner product \(\langle \Psi | \Psi \rangle\) is normalized to unity:

\langle \Psi | \Psi \rangle = \int_{-\infty}^\infty |\Psi(\mathbf{r}, t)|^2 d\tau = 1

The ket contains all knowable physical information about the system.

Postulate 2: Physical Observables and Hermitian Operators

To every physically measurable dynamical observable \(\mathcal{A}\) (position, momentum, kinetic energy, angular momentum) there corresponds a linear Hermitian operator \(\hat{A}\) acting in \(\mathcal{H}\):

\hat{A} = \hat{A}^\dagger \iff \langle \phi | \hat{A} \psi \rangle = \langle \hat{A} \phi | \psi \rangle

Hermiticity guarantees two foundational physical properties:

All eigenvalues \(a_n\) of \(\hat{A}\) are strictly **real numbers**: \(a_n \in \mathbb{R}\).
Eigenvectors \(|\psi_n\rangle\) belonging to distinct eigenvalues are mutually **orthogonal**: \(\langle \psi_m | \psi_n \rangle = \delta_{m n}\), forming a complete basis of \(\mathcal{H}\).

Postulate 3: Measurement Outcomes

The only possible numerical values obtained from an ideal measurement of physical observable \(\mathcal{A}\) are the eigenvalues \(a_n\) of the corresponding operator \(\hat{A}\):

\hat{A} |a_n\rangle = a_n |a_n\rangle

Postulate 4: Probabilistic Interpretation (The Born Rule)

When observable \(\mathcal{A}\) is measured on a system in normalized state \(|\Psi\rangle\), the probability \(P(a_n)\) of obtaining eigenvalue \(a_n\) (assuming non-degeneracy) is the absolute square of the projection amplitude:

P(a_n) = |\langle a_n | \Psi \rangle|^2 = \langle \Psi | \hat{P}_n | \Psi \rangle

where \(\hat{P}_n = |a_n\rangle\langle a_n|\) is the projection operator.

The quantum expectation value is:

\langle A \rangle = \langle \Psi | \hat{A} | \Psi \rangle = \sum_n a_n P(a_n)

Postulate 5: Wavefunction Collapse (State Reduction)

Immediately after a measurement of \(\mathcal{A}\) that yields the eigenvalue \(a_n\), the state vector collapses discontinuously into the corresponding eigenstate:

|\Psi_{\text{after}}\rangle = \frac{\hat{P}_n |\Psi\rangle}{\sqrt{\langle \Psi | \hat{P}_n | \Psi \rangle}} = |a_n\rangle

Postulate 6: Time Evolution (Schrödinger Equation)

Between measurements, the time evolution of the state vector is continuous, deterministic, and unitary, governed by the time-dependent Schrödinger equation:

i\hbar \frac{\partial}{\partial t} |\Psi(t)\rangle = \hat{H}(t) |\Psi(t)\rangle

For a time-independent Hamiltonian \(\hat{H}\), the formal solution is generated by the unitary evolution operator \(\hat{U}(t, 0) = e^{-i \hat{H} t / \hbar}\):

|\Psi(t)\rangle = \hat{U}(t, 0) |\Psi(0)\rangle = e^{-i \hat{H} t / \hbar} |\Psi(0)\rangle

Advanced Mathematical Supplement: Spectral Theorem & Rigged Hilbert Spaces

Mathematical Rigor: The Spectral Theorem for Unbounded Self-Adjoint Operators

In standard elementary quantum mechanics, physical operators such as position \(\hat{x}\) and momentum \(\hat{p}\) are often treated as simple linear matrices. However, in functional analysis, \(\hat{x}\) and \(\hat{p}\) are unbounded self-adjoint operators acting on the infinite-dimensional Hilbert space \(\mathcal{L}^2(\mathbb{R})\).

Domain Issues and Self-Adjointness vs Hermiticity

A linear operator \(\hat{A}\) with domain \(D(\hat{A}) \subset \mathcal{H}\) is symmetric (Hermitian) if:

\langle \phi | \hat{A} \psi \rangle = \langle \hat{A} \phi | \psi \rangle \quad \forall \phi, \psi \in D(\hat{A})

However, the adjoint operator \(\hat{A}^\dagger\) has domain:

D(\hat{A}^\dagger) = \{ \phi \in \mathcal{H} : \exists \eta \in \mathcal{H} \text{ such that } \langle \phi | \hat{A} \psi \rangle = \langle \eta | \psi \rangle \; \forall \psi \in D(\hat{A}) \}

An operator is truly self-adjoint if and only if:

\hat{A} = \hat{A}^\dagger \quad \text{and} \quad D(\hat{A}) = D(\hat{A}^\dagger)

By Stone's theorem, only strictly self-adjoint operators can generate strongly continuous one-parameter unitary groups representing physical time evolution:

\hat{U}(t) = \exp\left( -\frac{i \hat{H} t}{\hbar} \right)

The Spectral Decomposition

According to the Spectral Theorem for unbounded self-adjoint operators, there exists a unique projection-valued measure \(E(\lambda)\) on the Borel \(\sigma\)-algebra of \(\mathbb{R}\) such that:

\hat{A} = \int_{\sigma(\hat{A})} \lambda \, dE(\lambda)

where the spectrum \(\sigma(\hat{A}) = \sigma_{\text{point}}(\hat{A}) \cup \sigma_{\text{cont}}(\hat{A})\) decomposes into:

**Point Spectrum \(\sigma_{\text{point}}\)**: Discrete eigenvalues with normalizable eigenfunctions in \(\mathcal{L}^2(\mathbb{R})\) (bound states).
**Continuous Spectrum \(\sigma_{\text{cont}}\)**: Continuous eigenvalues whose generalized eigenfunctions lie outside \(\mathcal{L}^2(\mathbb{R})\) (scattering states).

The Rigged Hilbert Space (Gel'fand Triplet)

To place plane waves \(\psi_p(x) = \frac{1}{\sqrt{2\pi\hbar}} e^{i p x / \hbar}\) and Dirac delta functions \(\delta(x - x_0)\) on rigorous mathematical footing, quantum mechanics employs the Gel'fand Triplet:

\Phi \subset \mathcal{H} \subset \Phi^\times

\(\Phi\): The nuclear space of rapidly decreasing test functions (Schwartz space \(\mathcal{S}(\mathbb{R})\)).
\(\mathcal{H}\): The conventional Hilbert space of square-integrable functions \(\mathcal{L}^2(\mathbb{R})\).
\(\Phi^\times\): The dual space of tempered distributions (\(\mathcal{S}'(\mathbb{R})\)), containing plane waves and delta distributions.

In this rigged Hilbert space, Dirac's bra-ket formalism is mathematically rigorous, and the completeness relation holds identically:

\hat{I} = \int_{-\infty}^\infty |p\rangle\langle p| \, dp

Research Monograph: Quantum Decoherence & The Emergence of Classical Reality

Introduction & The Measurement Problem

In standard Copenhagen quantum mechanics, the linear unitary Schrödinger equation describes continuous deterministic wave evolution, while measurement introduces an ad hoc, discontinuous "collapse" of the wavefunction. For decades, the boundary between the quantum micro-world and classical macro-world remained an unresolved philosophical dilemma.

Modern quantum physics resolves this paradox through the theory of Quantum Decoherence, pioneered by H. Dieter Zeh, Wojciech Zurek, and Erich Joos.

Decoherence via Environmental Entanglement

No realistic macroscopic or mesoscopic quantum system is truly isolated. A central quantum system \(\mathcal{S}\) constantly scatters and interacts with an immense environmental bath \(\mathcal{E}\) of photons, air molecules, and phonon modes.

Let the system be in a spatial superposition:

|\Psi_{\mathcal{S}}\rangle = c_1 |x_1\rangle + c_2 |x_2\rangle

As the environment scatters off the system, the combined state becomes entangled:

|\Psi_{\mathcal{S}+\mathcal{E}}(t)\rangle = c_1 |x_1\rangle |E_1(t)\rangle + c_2 |x_2\rangle |E_2(t)\rangle

where \(|E_1(t)\rangle\) and \(|E_2(t)\rangle\) are normalized environmental states that have interacted with the particle at positions \(x_1\) and \(x_2\).

The Reduced Density Matrix & Off-Diagonal Decay

Because an observer cannot monitor the \(10^{20}\) environmental degrees of freedom, all observable properties of \(\mathcal{S}\) are obtained by tracing out the environment:

\hat{\rho}_{\mathcal{S}}(t) = \text{Tr}_{\mathcal{E}} \left( |\Psi_{\mathcal{S}+\mathcal{E}}(t)\rangle\langle \Psi_{\mathcal{S}+\mathcal{E}}(t)| \right) = |c_1|^2 |x_1\rangle\langle x_1| + |c_2|^2 |x_2\rangle\langle x_2| + c_1 c_2^ \langle E_2(t) | E_1(t) \rangle |x_1\rangle\langle x_2| + c_1^ c_2 \langle E_1(t) | E_2(t) \rangle |x_2\rangle\langle x_1|

Because the environmental states rapidly scatter into mutually orthogonal configurations:

\langle E_2(t) | E_1(t) \rangle \approx e^{-\Lambda_{\text{dec}} (x_1 - x_2)^2 t}

The off-diagonal coherence terms vanish with an astonishingly rapid decoherence timescale:

\tau_{\text{dec}} \sim \frac{\hbar^2}{m k_B T (x_1 - x_2)^2 \gamma}

For a macroscopic dust grain (\(10\text{ \mu m}\)) in air separated by \(1\text{ \mu m}\), \(\tau_{\text{dec}} \sim 10^{-31}\text{ seconds}\)—many orders of magnitude faster than any thermalization or dynamical timescale.

Decoherence dynamically destroys quantum phase interference, dynamically selecting the familiar localized classical position eigenstates (Einselection / Pointer States) without requiring wavefunction collapse.

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 1.1: Problem 1: Blackbody Peak Radiance & Solar Surface Temperature from Wien's Law

The solar emission spectrum measured outside Earth's atmosphere peaks at a maximum spectral wavelength \(\lambda_{\text{max}} = 502\text{ nm}\).

Calculate the effective surface blackbody temperature \(T_{\odot}\) of the Sun using Wien's displacement law constant \(b = 2.89777 \times 10^{-3}\text{ m}\cdot\text{K}\).
Derive the ratio of spectral energy densities predicted by the quantum Planck law versus the classical Rayleigh-Jeans formula at this peak wavelength \(\lambda_{\text{max}}\).

Foundational Example 1.2: Problem 2: Photoelectric Stopping Potential & Maximum Kinetic Energy of Potassium

Potassium metal has a work function \(\Phi = 2.29\text{ eV}\). Ultraviolet radiation of wavelength \(\lambda = 320\text{ nm}\) irradiates a clean potassium photocathode in high vacuum.

Determine the threshold frequency \(\nu_0\) and threshold wavelength \(\lambda_0\) for potassium.
Calculate the maximum kinetic energy \(K_{\text{max}}\) (in \(\text{eV}\) and Joules) of the emitted photoelectrons.
Calculate the stopping potential \(V_s\) required to bring the photocurrent to zero.
Calculate the maximum velocity \(v_{\text{max}}\) of the ejected electrons.

Intermediate Example 1.3: Problem 3: Compton Scattering Angle, Electron Recoil Energy & Relativistic Kinematics

A beam of monochromatic X-rays with incident wavelength \(\lambda_0 = 0.0500\text{ nm}\) is scattered through an angle \(\theta = 120^\circ\) by quasi-free electrons in a beryllium target.

Calculate the wavelength \(\lambda'\) of the Compton scattered X-rays.
Determine the kinetic energy \(K_e\) transferred to the recoil electron (in \(\text{eV}\)).
Calculate the recoil angle \(\phi\) of the scattered electron relative to the incident photon trajectory.

Comprehensive Multi-Step Solution:

Step 1: Wavelength of Scattered X-ray

The Compton formula is:

\lambda' - \lambda_0 = \lambda_C (1 - \cos\theta)

where the electron Compton wavelength is \(\lambda_C = \frac{h}{m_e c} = 2.42631 \times 10^{-12}\text{ m} = 0.0024263\text{ nm}\).

For \(\theta = 120^\circ\):

\cos(120^\circ) = -0.5 \implies 1 - \cos\theta = 1 - (-0.5) = 1.5

The Compton shift is:

\Delta\lambda = 0.0024263\text{ nm} \times 1.5 = 0.0036395\text{ nm}

The scattered wavelength is:

\lambda' = \lambda_0 + \Delta\lambda = 0.0500\text{ nm} + 0.0036395\text{ nm} = 0.05364\text{ nm} = 5.364 \times 10^{-11}\text{ m}

Step 2: Recoil Electron Kinetic Energy

From conservation of energy:

K_e = E_0 - E' = h c \left( \frac{1}{\lambda_0} - \frac{1}{\lambda'} \right) = h c \frac{\lambda' - \lambda_0}{\lambda_0 \lambda'}

Using \(h c = 1239.84\text{ eV}\cdot\text{nm}\):

E_0 = \frac{1239.84\text{ eV}\cdot\text{nm}}{0.0500\text{ nm}} = 24796.8\text{ eV} = 24.797\text{ keV}

E' = \frac{1239.84\text{ eV}\cdot\text{nm}}{0.0536395\text{ nm}} = 23114.3\text{ eV} = 23.114\text{ keV}

The kinetic energy imparted to the electron is:

K_e = 24796.8\text{ eV} - 23114.3\text{ eV} = 1682.5\text{ eV} \approx 1.683\text{ keV}

Step 3: Recoil Angle of the Electron

From momentum conservation along the transversal axis (perpendicular to incident photon):

0 = p' \sin\theta - p_e \sin\phi \implies p_e \sin\phi = \frac{h}{\lambda'} \sin\theta

Along the longitudinal axis:

p_0 = p' \cos\theta + p_e \cos\phi \implies p_e \cos\phi = \frac{h}{\lambda_0} - \frac{h}{\lambda'} \cos\theta

Dividing the transverse by longitudinal equations:

\tan\phi = \frac{\frac{1}{\lambda'} \sin\theta}{\frac{1}{\lambda_0} - \frac{1}{\lambda'} \cos\theta} = \frac{\sin\theta}{\left(\frac{\lambda'}{\lambda_0}\right) - \cos\theta}

Substitute \(\sin(120^\circ) = \frac{\sqrt{3}}{2} \approx 0.86603\), \(\cos(120^\circ) = -0.5\), and \(\frac{\lambda'}{\lambda_0} = \frac{0.0536395}{0.0500} = 1.07279\):

\tan\phi = \frac{0.86603}{1.07279 - (-0.5)} = \frac{0.86603}{1.57279} \approx 0.55063

Taking the arctangent:

\phi = \arctan(0.55063) \approx 28.84^\circ

Intermediate Example 1.4: Problem 4: de Broglie Wavelength of Thermal Neutrons & Macromolecular Crystallography

Thermal neutrons in a research nuclear reactor emerge in thermal equilibrium with heavy water at \(T = 300\text{ K}\).

Calculate the root-mean-square momentum \(p_{\text{rms}}\) and the corresponding thermal de Broglie wavelength \(\lambda_{\text{th}}\) of the neutrons (\(m_n = 1.67493 \times 10^{-27}\text{ kg}\)).
Compare this wavelength to the covalent bond length of a carbon-carbon single bond (\(d_{\text{C-C}} = 1.54\text{ Å}\)) and explain why thermal neutron diffraction is uniquely suited for locating hydrogen atoms in protein crystallography.

Comprehensive Multi-Step Solution:

Step 1: Thermal Root-Mean-Square Momentum & de Broglie Wavelength

From the equipartition theorem for a three-dimensional monoatomic gas:

\langle K \rangle = \frac{p_{\text{rms}}^2}{2 m_n} = \frac{3}{2} k_B T \implies p_{\text{rms}} = \sqrt{3 m_n k_B T}

Substitute \(m_n = 1.67493 \times 10^{-27}\text{ kg}\), \(k_B = 1.38065 \times 10^{-23}\text{ J/K}\), and \(T = 300\text{ K}\):

p_{\text{rms}} = \sqrt{3 \times (1.67493 \times 10^{-27}\text{ kg}) \times (1.38065 \times 10^{-23}\text{ J/K}) \times (300\text{ K})}

p_{\text{rms}} = \sqrt{2.0825 \times 10^{-47}} \approx 4.5634 \times 10^{-24}\text{ kg}\cdot\text{m/s}

The corresponding de Broglie wavelength is:

\lambda_{\text{th}} = \frac{h}{p_{\text{rms}}} = \frac{6.62607 \times 10^{-34}\text{ J}\cdot\text{s}}{4.5634 \times 10^{-24}\text{ kg}\cdot\text{m/s}} \approx 1.452 \times 10^{-10}\text{ m} = 1.452\text{ Å} = 0.1452\text{ nm}

Step 2: Comparison with Covalent Bond Length and Structural Significance

The thermal neutron wavelength (\(1.452\text{ Å}\)) is remarkably close to the canonical \(\text{C-C}\) single bond length (\(1.54\text{ Å}\)).

In X-ray crystallography, X-rays scatter off the electron cloud, meaning scattering power scales with atomic number squared (\(Z^2\)). Hydrogen (\(Z = 1\)) scatters X-rays exceptionally weakly compared to carbon (\(Z=6\)), nitrogen (\(Z=7\)), or oxygen (\(Z=8\)), making hydrogen positions notoriously difficult to resolve.
In neutron diffraction, neutrons scatter off atomic nuclei via the strong nuclear force. The coherent neutron scattering length for hydrogen (\(b_H = -3.74\text{ fm}\)) and deuterium (\(b_D = +6.67\text{ fm}\)) is of the same order of magnitude as carbon (\(b_C = 6.65\text{ fm}\)) and oxygen (\(b_O = 5.80\text{ fm}\)).
Combined with sub-angstrom thermal wavelengths that avoid radiation damage, thermal neutron diffraction enables precise three-dimensional localization of protonation states in enzyme active sites.

Advanced Example 1.5: Problem 5: Minimum Uncertainty Product for a Gaussian Wavepacket

A particle of mass \(m\) is described by the real normalized ground-state Gaussian wavepacket:

\psi(x) = \left( \frac{2\alpha}{\pi} \right)^{1/4} \exp(-\alpha x^2)

where \(\alpha > 0\) is a real constant.

Verify that \(\psi(x)\) is properly normalized on the interval \((-\infty, \infty)\).
Calculate the expectation values \(\langle x \rangle\), \(\langle x^2 \rangle\), and evaluate the position uncertainty \(\Delta x\).
Calculate the expectation values \(\langle p_x \rangle\), \(\langle p_x^2 \rangle\), and evaluate the momentum uncertainty \(\Delta p_x\).
Prove that the uncertainty product \(\Delta x \cdot \Delta p_x\) strictly achieves the minimum theoretical lower bound allowed by the Heisenberg uncertainty principle.

Comprehensive Multi-Step Solution:

Step 1: Verification of Normalization

Standard Gaussian integral identity:

\int_{-\infty}^\infty e^{-a x^2} dx = \sqrt{\frac{\pi}{a}}

Evaluate \(\langle \psi | \psi \rangle\):

\int_{-\infty}^\infty |\psi(x)|^2 dx = \left( \frac{2\alpha}{\pi} \right)^{1/2} \int_{-\infty}^\infty e^{-2\alpha x^2} dx = \sqrt{\frac{2\alpha}{\pi}} \cdot \sqrt{\frac{\pi}{2\alpha}} = 1

Hence \(\psi(x)\) is normalized.

Step 2: Evaluation of \(\Delta x\)

By symmetry, because \(|\psi(x)|^2\) is an even function of \(x\):

\langle x \rangle = \int_{-\infty}^\infty x |\psi(x)|^2 dx = 0

Using the standard Gaussian moment identity \(\int_{-\infty}^\infty x^2 e^{-a x^2} dx = \frac{1}{2a} \sqrt{\frac{\pi}{a}}\) with \(a = 2\alpha\):

\langle x^2 \rangle = \sqrt{\frac{2\alpha}{\pi}} \int_{-\infty}^\infty x^2 e^{-2\alpha x^2} dx = \sqrt{\frac{2\alpha}{\pi}} \cdot \left( \frac{1}{4\alpha} \sqrt{\frac{\pi}{2\alpha}} \right) = \frac{1}{4\alpha}

The spatial standard deviation is:

\Delta x = \sqrt{\langle x^2 \rangle - \langle x \rangle^2} = \sqrt{\frac{1}{4\alpha}} = \frac{1}{2\sqrt{\alpha}}

Step 3: Evaluation of \(\Delta p_x\)

The momentum operator is \(\hat{p}_x = -i\hbar \frac{d}{dx}\):

\frac{d\psi}{dx} = -2\alpha x \psi(x)

\langle p_x \rangle = -i\hbar \int_{-\infty}^\infty \psi(x) \left( -2\alpha x \psi(x) \right) dx = 0 \quad (\text{integrand is odd})

For \(\langle p_x^2 \rangle = -\hbar^2 \int_{-\infty}^\infty \psi(x) \frac{d^2\psi}{dx^2} dx\), integrating by parts:

\langle p_x^2 \rangle = \hbar^2 \int_{-\infty}^\infty \left| \frac{d\psi}{dx} \right|^2 dx = \hbar^2 \int_{-\infty}^\infty (-2\alpha x \psi)^2 dx = 4 \alpha^2 \hbar^2 \langle x^2 \rangle

Substituting \(\langle x^2 \rangle = \frac{1}{4\alpha}\):

\langle p_x^2 \rangle = 4 \alpha^2 \hbar^2 \left( \frac{1}{4\alpha} \right) = \alpha \hbar^2

The momentum standard deviation is:

\Delta p_x = \sqrt{\langle p_x^2 \rangle - \langle p_x \rangle^2} = \sqrt{\alpha \hbar^2} = \hbar \sqrt{\alpha}

Step 4: Heisenberg Uncertainty Product

Multiplying the two uncertainties:

\Delta x \cdot \Delta p_x = \left( \frac{1}{2\sqrt{\alpha}} \right) \cdot (\hbar \sqrt{\alpha}) = \frac{\hbar}{2}

The result exactly equals \(\frac{\hbar}{2}\). This proves that the Gaussian wavepacket is the unique state of minimum uncertainty.

Advanced Example 1.6: Problem 6: Proof of Hermiticity and Commutator Relations for Position and Kinetic Energy

Prove rigorously that the one-dimensional momentum operator \(\hat{p} = -i\hbar \frac{d}{dx}\) is Hermitian for square-integrable wavefunctions vanishing at infinity (\(\psi(\pm\infty) = 0\)).
Evaluate the commutator \([\hat{x}, \hat{T}]\) where \(\hat{T} = \frac{\hat{p}^2}{2m} = -\frac{\hbar^2}{2m} \frac{d^2}{dx^2}\) is the kinetic energy operator.
Use your result to derive the Ehrenfest theorem expression for the time evolution of the position expectation value \(\frac{d\langle x \rangle}{dt}\).

Comprehensive Multi-Step Solution:

Step 1: Proof of Hermiticity of \(\hat{p}\)

To be Hermitian, \(\hat{p}\) must satisfy:

\langle \phi | \hat{p} \psi \rangle = \langle \hat{p} \phi | \psi \rangle \iff \int_{-\infty}^\infty \phi^ \left(-i\hbar \frac{d\psi}{dx}\right) dx = \int_{-\infty}^\infty \left(-i\hbar \frac{d\phi}{dx}\right)^ \psi dx

Evaluate the left-hand integral using integration by parts:

\int_{-\infty}^\infty \phi^ \left(-i\hbar \frac{d\psi}{dx}\right) dx = -i\hbar \left[ \phi^ \psi \right]_{-\infty}^\infty - (-i\hbar) \int_{-\infty}^\infty \frac{d\phi^*}{dx} \psi dx

Because \(\phi, \psi \in L^2(\mathbb{R})\), the boundary term vanishes: \(\left[ \phi^* \psi \right]_{-\infty}^\infty = 0\).

Now rewrite the remaining integral:

i\hbar \int_{-\infty}^\infty \frac{d\phi^}{dx} \psi dx = \int_{-\infty}^\infty \left( -i\hbar \frac{d\phi}{dx} \right)^ \psi dx = \langle \hat{p}\phi | \psi \rangle

Hence \(\hat{p}^\dagger = \hat{p}\), proving that the momentum operator is Hermitian.

Step 2: Evaluation of \([\hat{x}, \hat{T}]\)

The kinetic energy operator is \(\hat{T} = \frac{\hat{p}^2}{2m}\).

Using the commutator identity \([\hat{A}, \hat{B}\hat{C}] = [\hat{A}, \hat{B}]\hat{C} + \hat{B}[\hat{A}, \hat{C}]\):

[\hat{x}, \hat{p}^2] = [\hat{x}, \hat{p}]\hat{p} + \hat{p}[\hat{x}, \hat{p}]

Since the canonical commutator is \([\hat{x}, \hat{p}] = i\hbar \hat{I}\):

[\hat{x}, \hat{p}^2] = (i\hbar)\hat{p} + \hat{p}(i\hbar) = 2i\hbar \hat{p}

Therefore:

[\hat{x}, \hat{T}] = \frac{1}{2m} [\hat{x}, \hat{p}^2] = \frac{1}{2m} (2i\hbar \hat{p}) = \frac{i\hbar}{m} \hat{p}

Step 3: Ehrenfest's Theorem for \(\frac{d\langle x \rangle}{dt}\)

For any time-independent operator \(\hat{A}\), the general equation of motion for expectation values is:

\frac{d\langle \hat{A} \rangle}{dt} = \frac{i}{\hbar} \langle [\hat{H}, \hat{A}] \rangle

For position \(\hat{A} = \hat{x}\) in a local potential \(\hat{H} = \hat{T} + V(\hat{x})\):

Since \([\hat{x}, V(\hat{x})] = 0\):

[\hat{H}, \hat{x}] = [\hat{T}, \hat{x}] = -[\hat{x}, \hat{T}] = -\frac{i\hbar}{m} \hat{p}

Substitute into Ehrenfest's theorem:

\frac{d\langle x \rangle}{dt} = \frac{i}{\hbar} \left\langle -\frac{i\hbar}{m} \hat{p} \right\rangle = \frac{\langle p \rangle}{m}

This reproduces the classical kinematic velocity relation \(\frac{dx}{dt} = \frac{p}{m}\).

Advanced Example 1.7: Problem 7: Non-Stationary State Superposition Dynamics & Bohr Beat Frequency

A quantum particle in a one-dimensional system is prepared at \(t = 0\) in a normalized non-stationary state consisting of an equal linear superposition of the ground state \(\psi_1(x)\) (energy \(E_1\)) and first excited state \(\psi_2(x)\) (energy \(E_2\)):

\Psi(x, 0) = \frac{1}{\sqrt{2}} \psi_1(x) + \frac{1}{\sqrt{2}} \psi_2(x)

Write the full time-dependent state \(\Psi(x, t)\).
Derive the time-dependent probability density \(P(x, t) = |\Psi(x, t)|^2\) and identify the oscillation frequency \(\omega_{21}\).
If \(E_1 = 2.0\text{ eV}\) and \(E_2 = 5.5\text{ eV}\), calculate the beat frequency \(\nu_{21}\) (in THz) and oscillation period \(T_{\text{osc}}\) (in fs).
Evaluate the expectation value of the Hamiltonian \(\langle \hat{H} \rangle\) as a function of time.

Comprehensive Multi-Step Solution:

Step 1: Time-Dependent Wavefunction

Each stationary eigenstate evolves with its phase factor \(\exp(-i E_n t / \hbar)\):

\Psi(x, t) = \frac{1}{\sqrt{2}} \psi_1(x) e^{-i E_1 t / \hbar} + \frac{1}{\sqrt{2}} \psi_2(x) e^{-i E_2 t / \hbar}

Step 2: Probability Density & Beat Frequency

Taking the squared modulus:

P(x, t) = \Psi^*(x, t) \Psi(x, t)

P(x, t) = \left( \frac{1}{\sqrt{2}} \psi_1^ e^{i E_1 t / \hbar} + \frac{1}{\sqrt{2}} \psi_2^ e^{i E_2 t / \hbar} \right) \left( \frac{1}{\sqrt{2}} \psi_1 e^{-i E_1 t / \hbar} + \frac{1}{\sqrt{2}} \psi_2 e^{-i E_2 t / \hbar} \right)

Assuming real spatial eigenfunctions \(\psi_1, \psi_2\):

P(x, t) = \frac{1}{2} \psi_1^2(x) + \frac{1}{2} \psi_2^2(x) + \frac{1}{2} \psi_1(x) \psi_2(x) \left( e^{-i (E_2 - E_1) t / \hbar} + e^{i (E_2 - E_1) t / \hbar} \right)

Using Euler's identity \(e^{i\theta} + e^{-i\theta} = 2\cos\theta\):

P(x, t) = \frac{1}{2} \psi_1^2(x) + \frac{1}{2} \psi_2^2(x) + \psi_1(x) \psi_2(x) \cos(\omega_{21} t)

where the Bohr transition angular frequency is:

\omega_{21} = \frac{E_2 - E_1}{\hbar}

Step 3: Numerical Beat Frequency and Period

Given:

\Delta E = E_2 - E_1 = 5.5\text{ eV} - 2.0\text{ eV} = 3.5\text{ eV} = 3.5 \times 1.60218 \times 10^{-19}\text{ J} = 5.6076 \times 10^{-19}\text{ J}

The linear oscillation frequency is:

\nu_{21} = \frac{\Delta E}{h} = \frac{5.6076 \times 10^{-19}\text{ J}}{6.62607 \times 10^{-34}\text{ J}\cdot\text{s}} \approx 8.463 \times 10^{14}\text{ Hz} = 846.3\text{ THz}

The oscillation period is:

T_{\text{osc}} = \frac{1}{\nu_{21}} = \frac{1}{8.463 \times 10^{14}\text{ s}^{-1}} \approx 1.182 \times 10^{-15}\text{ s} = 1.182\text{ fs}

Step 4: Expectation Value of Energy

\langle \hat{H} \rangle = \int \Psi^*(x, t) \hat{H} \Psi(x, t) dx

Since \(\hat{H} \psi_n = E_n \psi_n\) and \(\langle \psi_m | \psi_n \rangle = \delta_{mn}\):

\langle \hat{H} \rangle = |c_1|^2 E_1 + |c_2|^2 E_2 = \left(\frac{1}{\sqrt{2}}\right)^2 E_1 + \left(\frac{1}{\sqrt{2}}\right)^2 E_2 = \frac{E_1 + E_2}{2}

\langle \hat{H} \rangle = \frac{2.0\text{ eV} + 5.5\text{ eV}}{2} = 3.75\text{ eV}

The energy expectation value is completely independent of time, demonstrating energy conservation in isolated quantum systems.

Advanced Example 1.8: Commutator Algebra, Ehrenfest Theorem & Heisenberg Uncertainty

For the position operator \(\hat{x}\) and momentum operator \(\hat{p}_x = -i\hbar \frac{\partial}{\partial x}\), evaluate the fundamental commutator \([\hat{x}, \hat{p}_x]\) and the commutator \([\hat{x}^2, \hat{p}_x]\).
Using the generalized Robertson-Schrödinger uncertainty relation \(\sigma_A \sigma_B \ge \frac{1}{2} |\langle [\hat{A}, \hat{B}] \rangle|\), prove that \(\Delta x \Delta p_x \ge \frac{\hbar}{2}\).
State Ehrenfest's theorem for the expectation value of momentum \(\frac{d\langle p_x \rangle}{dt}\) and show how it recovers Newton's second law of motion in the classical limit.

Comprehensive Multi-Step Solution:

Step 1: Commutator Evaluations

**Fundamental commutator \([\hat{x}, \hat{p}_x]\):**

Apply the commutator to an arbitrary test function \(f(x)\):

[\hat{x}, \hat{p}_x] f(x) = \hat{x} (\hat{p}_x f) - \hat{p}_x (\hat{x} f) = x \left( -i\hbar \frac{\partial f}{\partial x} \right) - \left( -i\hbar \frac{\partial}{\partial x}(x f) \right)

Using the product rule \(\frac{\partial}{\partial x}(x f) = f + x \frac{\partial f}{\partial x}\):

[\hat{x}, \hat{p}_x] f(x) = -i\hbar x \frac{\partial f}{\partial x} + i\hbar \left( f + x \frac{\partial f}{\partial x} \right) = i\hbar f(x)

Since this holds for any well-behaved function \(f(x)\):

[\hat{x}, \hat{p}_x] = i\hbar \hat{I}

**Commutator \([\hat{x}^2, \hat{p}_x]\):**

Using the operator identity \([\hat{A}\hat{B}, \hat{C}] = \hat{A}[\hat{B}, \hat{C}] + [\hat{A}, \hat{C}]\hat{B}\) with \(\hat{A} = \hat{B} = \hat{x}\) and \(\hat{C} = \hat{p}_x\):

[\hat{x}^2, \hat{p}_x] = \hat{x} [\hat{x}, \hat{p}_x] + [\hat{x}, \hat{p}_x] \hat{x} = \hat{x} (i\hbar) + (i\hbar) \hat{x} = 2i\hbar \hat{x}

Step 2: Proof of Heisenberg Uncertainty Relation

The Robertson-Schrödinger uncertainty theorem for any two Hermitian operators \(\hat{A}\) and \(\hat{B}\) states:

\sigma_A^2 \sigma_B^2 \ge \left( \frac{1}{2i} \langle [\hat{A}, \hat{B}] \rangle \right)^2 + \left( \frac{1}{2} \langle \{ \Delta\hat{A}, \Delta\hat{B} \} \rangle \right)^2 \ge \frac{1}{4} |\langle [\hat{A}, \hat{B}] \rangle|^2

Substituting \(\hat{A} = \hat{x}\) and \(\hat{B} = \hat{p}_x\):

[\hat{x}, \hat{p}_x] = i\hbar \implies |\langle [\hat{x}, \hat{p}_x] \rangle| = |i\hbar| = \hbar

Therefore:

\sigma_x^2 \sigma_p^2 \ge \frac{1}{4} \hbar^2 \implies \sigma_x \sigma_p \ge \frac{\hbar}{2}

Writing in standard notation \(\Delta x \equiv \sigma_x\) and \(\Delta p_x \equiv \sigma_p\):

\Delta x \Delta p_x \ge \frac{\hbar}{2}

The lower bound \(\hbar/2\) is attained strictly by minimum-uncertainty Gaussian wavepackets.

Step 3: Ehrenfest Theorem and Classical Limit

For any time-independent operator \(\hat{A}\), the Heisenberg equation of motion yields:

\frac{d\langle A \rangle}{dt} = \frac{i}{\hbar} \langle [\hat{H}, \hat{A}] \rangle

For the momentum operator \(\hat{p}_x\) with Hamiltonian \(\hat{H} = \frac{\hat{p}_x^2}{2m} + V(\hat{x})\):

[\hat{H}, \hat{p}_x] = [V(\hat{x}), \hat{p}_x] = i\hbar \frac{\partial V}{\partial x}

Therefore:

\frac{d\langle p_x \rangle}{dt} = \frac{i}{\hbar} \langle i\hbar \frac{\partial V}{\partial x} \rangle = -\left\langle \frac{\partial V}{\partial x} \right\rangle = \langle F_x \rangle

In the classical macroscopic limit, the wavepacket is tightly localized (\(\Delta x \rightarrow 0\)), so \(\langle \frac{\partial V}{\partial x} \rangle \approx \frac{\partial V}{\partial x}\Big|_{\langle x \rangle}\).

The equation reduces to:

\frac{d\langle p_x \rangle}{dt} = -\frac{dV}{dx}\Big|_{\langle x \rangle} = F_x(\langle x \rangle)

This precisely recovers Newton's second law \(\mathbf{F} = \frac{d\mathbf{p}}{dt}\), proving that classical mechanics emerges as the macroscopic expectation value of quantum mechanics.

Advanced Example 1.9: Wavepacket Dispersion and Group Velocity of Free Matter Waves

A free electron of mass \(m\) is described at \(t = 0\) by a normalized one-dimensional Gaussian wavepacket of initial spatial width \(\sigma_0\):

\psi(x, 0) = \left( \frac{1}{2\pi \sigma_0^2} \right)^{1/4} \exp\left( -\frac{x^2}{4\sigma_0^2} + \frac{i p_0 x}{\hbar} \right)

Calculate the initial momentum expectation value \(\langle p \rangle\), momentum uncertainty \(\sigma_p\), and verify that \(\sigma_x \sigma_p = \hbar / 2\).
By solving the time-dependent Schrödinger equation, express the width of the wavepacket \(\sigma(t)\) as a function of time \(t\).
If an electron has initial localization \(\sigma_0 = 1.00\text{ \AA}\) (\(10^{-10}\text{ m}\)), calculate the time \(\tau\) required for the wavepacket's spatial width to double.

Comprehensive Multi-Step Solution:

Step 1: Initial Expectation Values and Minimum Uncertainty

The initial probability density is:

|\psi(x, 0)|^2 = \left( \frac{1}{2\pi \sigma_0^2} \right)^{1/2} \exp\left( -\frac{x^2}{2\sigma_0^2} \right)

By symmetry:

\langle x \rangle = 0, \quad \langle x^2 \rangle = \sigma_0^2 \implies \sigma_x = \sigma_0

The momentum expectation value is:

\langle p \rangle = \int_{-\infty}^\infty \psi^* \left( -i\hbar \frac{\partial}{\partial x} \right) \psi \, dx = p_0

For a Gaussian wavefunction, the momentum variance is:

\langle p^2 \rangle = p_0^2 + \frac{\hbar^2}{4\sigma_0^2} \implies \sigma_p^2 = \langle p^2 \rangle - \langle p \rangle^2 = \frac{\hbar^2}{4\sigma_0^2} \implies \sigma_p = \frac{\hbar}{2\sigma_0}

The uncertainty product is:

\sigma_x \sigma_p = \sigma_0 \left( \frac{\hbar}{2\sigma_0} \right) = \frac{\hbar}{2}

The Gaussian wavepacket is an exact minimum-uncertainty state.

Step 2: Time Evolution and Wavepacket Dispersion

The free-particle propagator is \(K(x, t; x', 0) = \sqrt{\frac{m}{2\pi i \hbar t}} \exp\left[ \frac{i m (x - x')^2}{2\hbar t} \right]\).

Performing the Gaussian convolution integral yields the time-dependent wavefunction:

\psi(x, t) = \left( \frac{1}{2\pi \sigma_0^2} \right)^{1/4} \frac{1}{\sqrt{1 + i \hbar t / 2 m \sigma_0^2}} \exp\left[ -\frac{(x - p_0 t / m)^2}{4\sigma_0^2 (1 + i \hbar t / 2 m \sigma_0^2)} + \frac{i p_0 x}{\hbar} - \frac{i p_0^2 t}{2 m \hbar} \right]

The probability density remains Gaussian, with its center moving at the classical group velocity \(v_g = p_0 / m\):

|\psi(x, t)|^2 = \frac{1}{\sqrt{2\pi [\sigma(t)]^2}} \exp\left[ -\frac{(x - v_g t)^2}{2 [\sigma(t)]^2} \right]

where the time-dependent spatial width is:

\sigma(t) = \sigma_0 \sqrt{1 + \left( \frac{\hbar t}{2 m \sigma_0^2} \right)^2}

Step 3: Doubling Time Calculation for an Electron

We require \(\sigma(\tau) = 2 \sigma_0\):

\sqrt{1 + \left( \frac{\hbar \tau}{2 m \sigma_0^2} \right)^2} = 2 \implies 1 + \left( \frac{\hbar \tau}{2 m \sigma_0^2} \right)^2 = 4 \implies \left( \frac{\hbar \tau}{2 m \sigma_0^2} \right)^2 = 3

\frac{\hbar \tau}{2 m \sigma_0^2} = \sqrt{3} \implies \tau = \frac{2\sqrt{3} m \sigma_0^2}{\hbar}

Substitute \(m = 9.10938 \times 10^{-31}\text{ kg}\), \(\sigma_0 = 1.00 \times 10^{-10}\text{ m}\), and \(\hbar = 1.05457 \times 10^{-34}\text{ J}\cdot\text{s}\):

\tau = \frac{2 \sqrt{3} (9.10938 \times 10^{-31}\text{ kg})(1.00 \times 10^{-20}\text{ m}^2)}{1.05457 \times 10^{-34}\text{ J}\cdot\text{s}}

\tau = \frac{3.1556 \times 10^{-50}}{1.05457 \times 10^{-34}} \approx 2.99 \times 10^{-16}\text{ seconds} \approx 0.30\text{ femtoseconds}

An electron localized within an atomic dimension (\(1\text{ \AA}\)) doubles its spatial spread in just 0.3 femtoseconds, illustrating how extreme quantum dispersion is for microscopic particles.

The Classical Rayleigh-Jeans Radiation Law

The Ultraviolet Catastrophe

Planck's Revolutionary Quantization Hypothesis

Asymptotic Limits and Wien's Displacement Law

The Photoelectric Effect

Einstein's Photoelectric Equation

The Compton Effect

The de Broglie Relations

Thermal and Non-Relativistic Kinetic Scaling

Experimental Verification: Davisson-Germer & Thomson Experiments

Mathematical Origin from Fourier Transform Pairs

Generalized Robertson-Schrödinger Uncertainty Relation

Time-Energy Uncertainty Relation

The Time-Dependent Schrödinger Equation (TDSE)

Separation of Variables & The Time-Independent Schrödinger Equation

Properties of Stationary States

Operator Algebra and Linearity

Hermitian (Self-Adjoint) Operators

Commutators and Simultaneous Observables

The Six Foundational Postulates

Postulate 1: The State Space

Postulate 2: Physical Observables

Postulate 3: Measurement Eigenvalues

Postulate 4: Measurement Probabilities & State Collapse

Postulate 5: Time Evolution

Postulate 6: Pauli Antisymmetry (Symmetrization Postulate)

Postulate 1: The State of a Quantum System

Postulate 2: Physical Observables and Hermitian Operators

Postulate 3: Measurement Outcomes

Postulate 4: Probabilistic Interpretation (The Born Rule)

Postulate 5: Wavefunction Collapse (State Reduction)

Postulate 6: Time Evolution (Schrödinger Equation)

Advanced Mathematical Supplement: Spectral Theorem & Rigged Hilbert Spaces

Mathematical Rigor: The Spectral Theorem for Unbounded Self-Adjoint Operators

Domain Issues and Self-Adjointness vs Hermiticity

The Spectral Decomposition

The Rigged Hilbert Space (Gel'fand Triplet)

Research Monograph: Quantum Decoherence & The Emergence of Classical Reality

Introduction & The Measurement Problem

Decoherence via Environmental Entanglement

The Reduced Density Matrix & Off-Diagonal Decay

Worked Problems & Step-by-Step Quantum Derivations

Comprehensive Multi-Step Solution:

Step 1: Effective Surface Temperature Calculation

Step 2: Planck to Rayleigh-Jeans Ratio at \(\lambda_{\text{max}}\)

Comprehensive Multi-Step Solution:

Step 1: Threshold Frequency and Wavelength

Step 2: Incident Photon Energy and Maximum Kinetic Energy

Step 3: Stopping Potential

Step 4: Maximum Electron Velocity

Comprehensive Multi-Step Solution:

Step 1: Wavelength of Scattered X-ray

Step 2: Recoil Electron Kinetic Energy

Step 3: Recoil Angle of the Electron

Comprehensive Multi-Step Solution:

Step 1: Thermal Root-Mean-Square Momentum & de Broglie Wavelength

Step 2: Comparison with Covalent Bond Length and Structural Significance

Comprehensive Multi-Step Solution:

Step 1: Verification of Normalization

Step 2: Evaluation of \(\Delta x\)

Step 3: Evaluation of \(\Delta p_x\)

Step 4: Heisenberg Uncertainty Product

Comprehensive Multi-Step Solution:

Step 1: Proof of Hermiticity of \(\hat{p}\)

Step 2: Evaluation of \([\hat{x}, \hat{T}]\)

Step 3: Ehrenfest's Theorem for \(\frac{d\langle x \rangle}{dt}\)

Comprehensive Multi-Step Solution:

Step 1: Time-Dependent Wavefunction

Step 2: Probability Density & Beat Frequency

Step 3: Numerical Beat Frequency and Period

Step 4: Expectation Value of Energy

Comprehensive Multi-Step Solution:

Step 1: Commutator Evaluations

Step 2: Proof of Heisenberg Uncertainty Relation

Step 3: Ehrenfest Theorem and Classical Limit

Comprehensive Multi-Step Solution:

Step 1: Initial Expectation Values and Minimum Uncertainty

Step 2: Time Evolution and Wavepacket Dispersion

Step 3: Doubling Time Calculation for an Electron