In the semiclassical regime, where the action $S$ is large compared to $\hbar$, we can find approximate solutions to the time-independent Schrödinger equation. The Wentzel-Kramers-Brillouin (WKB) approximation provides a powerful method for this. We start by positing a solution of the form $\psi(x) = e^{iS(x)/\hbar}$ and expand the action $S(x)$ as a power series in $\hbar$: $S(x) = S_0(x) + \hbar S_1(x) + \dots$.

Substituting this into the Schrödinger equation and collecting terms of the same order in $\hbar$ allows us to solve for the components of the action. To lowest order, we find that the momentum is purely real in classically allowed regions ($E > V(x)$) and purely imaginary in classically forbidden regions ($E < V(x)$).

The WKB approximation fails at classical turning points, where $E = V(x)$ and thus $p(x) = 0$. The wave function diverges at these points. To connect the solutions across a turning point, we linearize the potential in its vicinity: $V(x) \approx E + V'(x_t)(x - x_t)$.

Substituting this linear potential into the Schrödinger equation yields a standard differential equation.

The two linearly independent solutions are the Airy functions $Ai(z)$ and $Bi(z)$. For a potential barrier on the right, we require a decaying solution, so we choose $Ai(z)$. By examining the asymptotic forms of $Ai(z)$ for large $|z|$, we can match them to the WKB solutions. This procedure yields the connection formulae, which relate the amplitudes and phases of the wavefunctions across the turning points. This is critical for problems like calculating tunneling rates or finding the quantized energy levels of bound states.

While the wavefunction lives in configuration space, semiclassical mechanics finds a natural home in phase space. The Wigner-Weyl transform maps quantum operators into c-number functions on phase space. Applying this transform to the density operator $\hat{\rho}$ yields the Wigner quasi-probability distribution, $P(x,p)$.

Crucially, the Wigner function is not a true probability distribution, as it can take on negative values. These negative regions are a signature of quantum interference and have no classical analogue. The time evolution of the Wigner function is governed by an equation involving the Moyal bracket of the Wigner-transformed Hamiltonian and the Wigner function itself.

These advanced semiclassical methods are indispensable for analyzing complex quantum systems, from calculating tunneling rates in molecules to understanding decoherence in quantum computing.

Time to check your understanding of these advanced concepts.

The interplay between classical intuition and quantum reality, as described by these methods, continues to be a fertile ground for research.