Schrödinger Wave Equation and Model Systems

1. The Schrödinger Wave Equation

The Schrödinger wave equation is a fundamental equation in quantum mechanics that describes how the quantum state of a physical system changes over time. It is the quantum mechanical analogue of Newton's second law of motion in classical mechanics. The equation allows us to calculate the wave function of a system, which contains information about the probability of finding a particle in a certain location and its energy.

There are two main forms of the Schrödinger equation:

1.1. Time-Dependent Schrödinger Equation (TDSE)

The time-dependent Schrödinger equation describes the evolution of a quantum system over time. It is given by:

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

Where:

• $i$ is the imaginary unit ($\sqrt{-1}$).
• $\hbar$ (h-bar) is the reduced Planck constant, equal to $h/2\pi$, where $h$ is Planck's constant.
• $\frac{\partial}{\partial t}$ is the partial derivative with respect to time.
• $\Psi(\mathbf{r}, t)$ is the wave function of the system, which depends on position $\mathbf{r}$ and time $t$. It is a complex-valued function.
• $\hat{H}$ is the Hamiltonian operator, which represents the total energy of the system (kinetic energy + potential energy).

The square of the magnitude of the wave function, $|\Psi(\mathbf{r}, t)|^2$, represents the probability density of finding the particle at position $\mathbf{r}$ at time $t$. The TDSE is used when the Hamiltonian operator changes with time.

1.2. Time-Independent Schrödinger Equation (TISE)

The time-independent Schrödinger equation is used for systems where the potential energy does not change with time. In such cases, the wave function can be separated into spatial and temporal parts, leading to a simpler equation that determines the stationary states (states with definite energy) of the system. It is given by:

$\hat{H} \Psi(\mathbf{r}) = E \Psi(\mathbf{r})$

Where:

• $\hat{H}$ is the time-independent Hamiltonian operator.
• $\Psi(\mathbf{r})$ is the time-independent wave function, depending only on position.
• $E$ is the energy of the system, which is a constant for a stationary state.

This equation is an eigenvalue equation, where the Hamiltonian operator acts on the wave function to yield the same wave function multiplied by a scalar, which is the energy eigenvalue. Solving this equation yields the possible energy levels of the system and the corresponding wave functions.

The Hamiltonian operator for a single particle of mass $m$ in three dimensions is typically given by:

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

Where:

• $\hat{T}$ is the kinetic energy operator: $-\frac{\hbar^2}{2m} \nabla^2$.
• $\nabla^2$ is the Laplacian operator ($\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2}$).
• $\hat{V}$ is the potential energy operator, $V(\mathbf{r})$.

For the time-independent Schrödinger equation, this becomes:

$-\frac{\hbar^2}{2m} \nabla^2 \Psi(\mathbf{r}) + V(\mathbf{r}) \Psi(\mathbf{r}) = E \Psi(\mathbf{r})$

2. Model Systems

To understand the application of the Schrödinger equation, we often study simple, idealized systems called model systems. These systems have simplified potential energy functions, making them solvable. The solutions provide insights into the quantization of energy and the behavior of particles at the atomic and molecular level.

2.1. Particle in a One-Dimensional Box

This is one of the simplest and most important model systems. It considers a particle of mass $m$ confined to move along a line (say, the x-axis) between two impenetrable walls at $x=0$ and $x=L$. Inside the box ($0 < x < L$), the potential energy is zero ($V(x) = 0$). Outside the box ($x \le 0$ or $x \ge L$), the potential energy is infinite ($V(x) = \infty$), meaning the particle cannot exist there.

The time-independent Schrödinger equation inside the box ($0 < x < L$) is:

$-\frac{\hbar^2}{2m} \frac{d^2\Psi(x)}{dx^2} + 0 \cdot \Psi(x) = E \Psi(x)$

Which simplifies to:

$\frac{d^2\Psi(x)}{dx^2} = -\frac{2mE}{\hbar^2} \Psi(x)$

Let $k^2 = \frac{2mE}{\hbar^2}$. The equation becomes:

$\frac{d^2\Psi(x)}{dx^2} = -k^2 \Psi(x)$

The general solution to this second-order differential equation is:

$\Psi(x) = A \sin(kx) + B \cos(kx)$

Where $A$ and $B$ are constants determined by boundary conditions.

The boundary conditions are derived from the fact that the wave function must be continuous. Since the particle cannot be outside the box, the probability of finding it there must be zero. This means $\Psi(x) = 0$ at the boundaries.

• At $x=0$: $\Psi(0) = A \sin(0) + B \cos(0) = A \cdot 0 + B \cdot 1 = B$. Since $\Psi(0) = 0$, we must have $B=0$.
• At $x=L$: $\Psi(L) = A \sin(kL)$. Since $\Psi(L) = 0$ and $A$ cannot be zero (otherwise the wave function would be zero everywhere, meaning no particle), we must have $\sin(kL) = 0$.

The condition $\sin(kL) = 0$ implies that $kL$ must be an integer multiple of $\pi$.

$kL = n\pi$, where $n = 1, 2, 3, \dots$ (n cannot be 0, as it would lead to $k=0$ and $\Psi(x)=0$).

So, $k = \frac{n\pi}{L}$.

Now substitute this back into the expression for $k^2$ to find the allowed energies ($E$):

$k^2 = \frac{2mE}{\hbar^2} \implies \left(\frac{n\pi}{L}\right)^2 = \frac{2mE_n}{\hbar^2}$

$E_n = \frac{n^2\pi^2\hbar^2}{2mL^2}$

Since $\hbar = h/2\pi$, we can also write:

$E_n = \frac{n^2h^2}{8mL^2}$, where $n = 1, 2, 3, \dots$

This is a crucial result: the energy of the particle in a box is quantized, meaning it can only take on specific discrete values determined by the quantum number $n$. The lowest possible energy (ground state) occurs when $n=1$, and it is non-zero ($E_1 = \frac{h^2}{8mL^2}$), which is a consequence of the Heisenberg uncertainty principle.

The corresponding wave functions are:

$\Psi_n(x) = A \sin\left(\frac{n\pi x}{L}\right)$

To find the normalization constant $A$, we use the normalization condition: the total probability of finding the particle in the box must be 1.

$\int_{0}^{L} |\Psi_n(x)|^2 dx = 1$

$\int_{0}^{L} A^2 \sin^2\left(\frac{n\pi x}{L}\right) dx = 1$

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

$A^2 \int_{0}^{L} \frac{1 - \cos\left(\frac{2n\pi x}{L}\right)}{2} dx = 1$

$\frac{A^2}{2} \left[ x - \frac{L}{2n\pi} \sin\left(\frac{2n\pi x}{L}\right) \right]_{0}^{L} = 1$

$\frac{A^2}{2} \left[ (L - 0) - \left(\frac{L}{2n\pi} \sin(2n\pi) - \frac{L}{2n\pi} \sin(0)\right) \right] = 1$

Since $\sin(2n\pi) = 0$ and $\sin(0) = 0$ for integer $n$:

$\frac{A^2}{2} [L] = 1 \implies A^2 = \frac{2}{L} \implies A = \sqrt{\frac{2}{L}}$

So, the normalized wave functions are:

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

2.2. Rigid Rotator

The rigid rotator model describes a system of two masses, $m_1$ and $m_2$, connected by a rigid rod of fixed length $r$. This model is often used to approximate the rotation of diatomic molecules. The rotation occurs in three dimensions, and the system has rotational kinetic energy but no potential energy ($V=0$).

In three dimensions, the Schrödinger equation involves the Laplacian operator $\nabla^2$. For a system with no potential energy, the time-independent Schrödinger equation is:

$-\frac{\hbar^2}{2\mu} \nabla^2 \Psi(\theta, \phi) = E \Psi(\theta, \phi)$

Where $\mu$ is the reduced mass of the system, given by $\mu = \frac{m_1 m_2}{m_1 + m_2}$, and $\Psi$ depends on the spherical coordinates $\theta$ and $\phi$ (since the distance $r$ is fixed).

The solution to this equation leads to quantized energy levels. The energy of a rigid rotator depends on two quantum numbers, $J$ (the rotational quantum number) and $M_J$ (the magnetic quantum number).

The allowed energy levels are given by:

$E_J = \frac{\hbar^2}{2I} J(J+1)$

Where:

• $I$ is the moment of inertia of the rotator, $I = \mu r^2$.
• $J$ is the rotational quantum number, which can take integer values $J = 0, 1, 2, 3, \dots$.
• $M_J$ is the magnetic quantum number, which can take values from $-J$ to $+J$ in integer steps, i.e., $M_J = -J, -J+1, \dots, 0, \dots, J-1, J$. There are $2J+1$ possible values of $M_J$ for each $J$.

The energy $E_J$ depends only on $J$, not on $M_J$. This means that for a given $J$, there are $2J+1$ degenerate states (states with the same energy) corresponding to the different $M_J$ values.

The wave functions for the rigid rotator are spherical harmonics, denoted as $Y_{J, M_J}(\theta, \phi)$. These functions describe the angular distribution of the particle.

2.3. Harmonic Oscillator

The harmonic oscillator model describes a system that experiences a restoring force proportional to its displacement from an equilibrium position. A common example is a mass attached to a spring obeying Hooke's Law. In quantum mechanics, this model is used to describe vibrations of atoms in molecules and solids.

The potential energy for a one-dimensional harmonic oscillator is given by:

$V(x) = \frac{1}{2} kx^2$

Where $k$ is the force constant (related to the stiffness of the spring).

The time-independent Schrödinger equation for a particle of mass $m$ in this potential is:

$-\frac{\hbar^2}{2m} \frac{d^2\Psi(x)}{dx^2} + \frac{1}{2} kx^2 \Psi(x) = E \Psi(x)$

Solving this equation yields quantized energy levels. The solution is non-trivial and typically involves using creation and annihilation operators or solving the differential equation directly, which leads to Hermite polynomials.

The allowed energy levels for the harmonic oscillator are:

$E_v = \left(v + \frac{1}{2}\right) \hbar \omega$

Where:

• $v$ is the vibrational quantum number, which can take integer values $v = 0, 1, 2, 3, \dots$.
• $\omega$ is the classical angular frequency of the oscillator, given by $\omega = \sqrt{\frac{k}{m}}$.
• $\hbar$ is the reduced Planck constant.

The energy levels are equally spaced, with the spacing being $\hbar \omega$.

The ground state energy ($v=0$) is $E_0 = \frac{1}{2}\hbar\omega$. This non-zero minimum energy is known as the zero-point energy, a direct consequence of the uncertainty principle. It means that even at absolute zero temperature, the oscillator still possesses some vibrational energy and is never truly at rest.

The wave functions for the harmonic oscillator are given by:

$\Psi_v(x) = N_v H_v(\alpha x) e^{-\alpha^2 x^2 / 2}$

Where:

• $N_v$ is a normalization constant.
• $H_v(y)$ are the Hermite polynomials ($H_0(y)=1$, $H_1(y)=2y$, $H_2(y)=4y^2-2$, etc.).
• $\alpha = \left(\frac{m\omega}{\hbar}\right)^{1/2}$.

3. Significance and Applications

The Schrödinger equation and these model systems form the bedrock of quantum chemistry and physics. They allow us to understand and predict phenomena at the atomic and molecular level, such as:

• Quantization of energy in atoms and molecules.
• The nature of chemical bonds.
• Spectroscopic properties (absorption and emission of light), which are used to identify and study molecules.
• The behavior of electrons in materials.

For instance, the particle-in-a-box model, though simple, helps explain the delocalization of electrons in conjugated systems (like those in dyes and pigments) and the electronic spectra of such molecules. The rigid rotator is directly applicable to understanding microwave spectroscopy, which probes molecular rotations. The harmonic oscillator is fundamental to interpreting infrared spectroscopy, which reveals information about the vibrational modes of molecules.