Continuity Conditions of the Wave Function

In non-relativistic quantum mechanics, the wave function $\psi(x)$ serves as the fundamental description of a particle's state. To ensure that $\psi(x)$ yields physically meaningful results—such as well-defined probability densities and expectation values—it must satisfy several mathematical criteria: it must be square-integrable (normalizable), single-valued, and, crucially, continuous.

The requirement for continuity is not an arbitrary rule but is deeply rooted in the mathematical structure of the time-independent Schrödinger equation:

$$-\frac{\hbar^{2}}{2m}\frac{d^{2}\psi(x)}{dx^{2}}+V(x)\psi(x)=E\psi(x)$$

Because this is a second-order linear differential equation, the behavior of the wave function and its derivatives at any point $x_0$ is dictated by the nature of the potential $V(x)$. If the potential is well-behaved, the wave function must also be well-behaved to prevent the equation from becoming mathematically undefined.

The Mathematical Necessity of Continuity

To understand why $\psi(x)$ and its derivative $\psi'(x)$ must be continuous, we examine the relationship between the kinetic and potential energy terms.

If the potential $V(x)$ is a finite function at a point $x_0$, then the term $V(x)\psi(x)$ remains finite. For the equality in the Schrödinger equation to hold, the kinetic energy term—specifically the second derivative $\psi''(x)$—must also remain finite. If $\psi''(x)$ were to become infinite (a singularity), the equation would break down.

Deriving the Matching Conditions

We can rigorously demonstrate the continuity requirements by integrating the Schrödinger equation over a vanishingly small interval $[x_0 - \epsilon, x_0 + \epsilon]$ centered at $x_0$:

$$\int_{x_0-\epsilon}^{x_0+\epsilon} \left[ -\frac{\hbar^{2}}{2m}\frac{d^{2}\psi(x)}{dx^{2}} + V(x)\psi(x) \right] dx = E \int_{x_0-\epsilon}^{x_0+\epsilon} \psi(x) dx$$

As we take the limit $\epsilon \to 0$:

1. The integral on the right-hand side (the energy term) vanishes because the interval width goes to zero.
2. The integral of the potential term $V(x)\psi(x)$ also vanishes, provided $V(x)$ does not contain a singularity like a Dirac delta function.

This leaves us with the integral of the second derivative:

$$-\frac{\hbar^{2}}{2m} \int_{x_0-\epsilon}^{x_0+\epsilon} \frac{d^2\psi}{dx^2} dx = -\frac{\hbar^{2}}{2m} \left[ \psi'(x_0^+) - \psi'(x_0^-) \right] = 0$$

This result implies that $\psi'(x_0^+) = \psi'(x_0^-)$, meaning the first derivative must be continuous. Furthermore, if $\psi(x)$ itself were discontinuous at $x_0$, the first derivative would behave like a $\delta$-function, and the second derivative would behave like the derivative of a $\delta$-function, making the Schrödinger equation impossible to satisfy for any finite $V(x)$.

Classification of Continuity Requirements

The specific constraints on $\psi(x)$ and $\psi'(x)$ depend entirely on the type of potential encountered at the boundary.

1. Finite Potentials

For any potential that does not diverge to infinity or involve a delta function, both the wave function and its first derivative must be continuous. This allows for a smooth "matching" of solutions across different regions of space.

2. Infinite Potential Barriers

When $V(x) \to \infty$, the particle is strictly forbidden from entering that region. To keep the energy finite, the wave function must vanish ($\psi = 0$) within the infinite barrier. While $\psi(x)$ remains continuous at the boundary (approaching zero from both sides), the slope $\psi'(x)$ can change abruptly, as there is no mathematical requirement for the derivative to be smooth when the potential is infinite.

3. Delta Function Potentials

In the case of a Dirac delta potential, $V(x) = \lambda\delta(x - x_0)$, the potential is infinitely sharp. While $\psi(x)$ remains continuous to ensure the probability density is well-defined, the first derivative undergoes a specific jump discontinuity proportional to the strength of the delta function:

$$\psi'(x_0^+) - \psi'(x_0^-) = \frac{2m\lambda}{\hbar^2}\psi(x_0)$$

Case Studies in Quantum Mechanics

The Infinite Square Well

In an infinite square well of width $L$, the potential is zero inside and infinite outside. The continuity condition at the boundaries ($x=0$ and $x=L$) forces the wave function to be zero: $\psi(0) = \psi(L) = 0$. This leads to the quantized energy levels and the sinusoidal standing wave solutions:
$$\psi_n(x) = A \sin\left(\frac{n\pi x}{L}\right)$$
Here, the continuity of $\psi$ is satisfied, but $\psi'$ is clearly discontinuous at the walls.

The Finite Square Barrier

Consider a particle encountering a barrier of height $V_0$ and width $a$. Unlike the infinite well, the particle has a non-zero probability of tunneling through the barrier. To solve this, we must apply matching conditions at the boundaries $x=0$ and $x=a$:

• $\psi_{\text{left}}(0) = \psi_{\text{middle}}(0)$ and $\psi_{\text{left}}'(0) = \psi_{\text{middle}}'(0)$
• $\psi_{\text{middle}}(a) = \psi_{\text{right}}(a)$ and $\psi_{\text{middle}}'(a) = \psi_{\text{right}}'(a)$

Solving this system of equations allows us to calculate the transmission and reflection coefficients, which are central to understanding quantum tunneling.

The Delta Function Potential

The delta function is a powerful tool for modeling short-range forces. By applying the jump condition derived earlier, we can find bound states for an attractive delta potential ($\lambda < 0$). The discontinuity in $\psi'$ is the mathematical manifestation of the "kick" the particle receives from the singular potential.

Practical Considerations in Computational Physics

When solving the Schrödinger equation numerically (e.g., using the Finite Difference Method), the continuity conditions are not just theoretical constraints—they are essential for stability.

• Grid Discretization: If the wave function is not forced to be continuous across the discrete grid points, the numerical solver may produce "ghost" solutions or diverge due to artificial singularities.
• Boundary Conditions: In simulations of open systems, researchers often use Absorbing Boundary Conditions (ABC) to prevent waves from reflecting off the edges of the computational box. Even in these complex setups, the internal continuity of $\psi$ must be strictly maintained.
• Handling Discontinuities: When a simulation encounters a sharp potential jump, the mesh must be sufficiently fine, or specialized matching algorithms must be used to ensure the derivative jump is captured accurately without introducing numerical noise.

Summary

The continuity conditions of the wave function are the bridge between the abstract differential form of the Schrödinger equation and the physical reality of quantum states.

• For smooth, finite potentials, the wave function and its derivative must be continuous.
• For infinite barriers, the wave function must be continuous (and typically zero), but the derivative may jump.
• For singular delta potentials, the wave function remains continuous, but the derivative experiences a calculated jump.

Mastering these matching conditions is fundamental to solving nearly every standard problem in quantum mechanics, from the simplest bound states to complex scattering phenomena.