Wigner function

In classical physics, the state of a particle is perfectly described by its position and momentum, a single point in what is known as phase space. However, the quantum world, governed by Heisenberg's uncertainty principle, forbids the simultaneous precise knowledge of both properties, seemingly rendering the intuitive phase-space picture obsolete. This article addresses this fundamental conflict by introducing the Wigner function, a brilliant mathematical construction that acts as a bridge between the classical and quantum realms. Across the following chapters, we will explore the core principles and mechanisms of this "quasi-probability" distribution, uncovering how it visualizes quantum states and their dynamics. We will then journey through its surprisingly diverse applications, demonstrating how the Wigner function provides a unified phase-space language for phenomena ranging from classical optics and signal processing to materials science and astrophysics.

Imagine you want to describe a car moving down a road. At any instant, you could say where it is (its position) and how fast it's going (its momentum). You could plot this as a single point on a 2D map, with one axis for position and one for momentum. This map is what physicists call ​​phase space​​. As the car moves, this point traces a path, a trajectory, on the map. This is the simple, comfortable world of classical mechanics.

But in the quantum world, things are not so simple. Werner Heisenberg's uncertainty principle tells us that we cannot know both the precise position and the precise momentum of a particle at the same time. If you try to create a quantum state that is perfectly localized in position, its momentum becomes wildly uncertain, and vice versa. So, how can we even begin to draw a picture in phase space for a quantum particle? Does the very idea of a phase-space map become meaningless?

Not quite. In 1932, the brilliant physicist Eugene Wigner devised a remarkable way to navigate this quantum fog. He invented a mathematical tool, now called the ​​Wigner function​​, $W(x,p)$, which is the closest thing quantum mechanics has to a classical phase-space distribution. It's a function that lives on the classical phase-space map, and it tries to tell you the "likelihood" of finding a particle at position $x$ and momentum $p$. The way it's constructed is beautifully clever. Instead of looking at the wavefunction $\psi(x)$ at a single point, the Wigner function looks at a kind of correlation between the wavefunction at two points, $\psi(x+y)$ and $\psi(x-y)$, symmetric about the position $x$. It then uses a mathematical procedure called a Fourier transform to relate the separation $y$ to the momentum $p$. The formal definition is:

You don't need to digest the details of this integral. The key idea is this: the Wigner function is a bridge between the abstract wavefunction of quantum mechanics and the intuitive phase-space picture of classical mechanics. But it's a rickety bridge, with a fascinating twist. Wigner called it a "quasi-probability" distribution, and that "quasi" is where all the quantum magic lies.

To get a feel for this new tool, let's apply it to the simplest, most well-behaved quantum systems we can find. Consider the ground state—the state of lowest possible energy—of a quantum harmonic oscillator. You can think of this as the quantum version of a mass on a spring at rest, or a single mode of light in an optical cavity in its vacuum state. It's not truly "at rest," because the uncertainty principle requires it to have some minimal jiggle, but it's as calm as a quantum system can be.

If we calculate the Wigner function for this state, we get something wonderfully simple and elegant. It's a perfectly symmetric, positive "mound" centered at the origin of phase space $(x=0, p=0)$:

This is a two-dimensional Gaussian function. It looks exactly like what you might expect for a classical probability distribution of a particle that is most likely to be found at its equilibrium position with zero momentum, but with some random fluctuations around that point. The width of this Gaussian mound is determined by Planck's constant, $\hbar$, which sets the fundamental scale of this quantum jiggle.

Now, what if we "kick" this oscillator? In quantum mechanics, this can create what's known as a ​​coherent state​​. This is a special state that behaves in a very classical manner. If we calculate its Wigner function, we find that it is the exact same Gaussian mound as the ground state, but simply shifted in phase space to be centered on the classical position $x_0$ and momentum $p_0$. As time goes on, the center of this mound follows the exact elliptical path that a classical pendulum would trace in phase space. The Wigner function provides a stunningly literal picture of a localized quantum wave packet behaving like a classical particle.

These first examples are encouragingly classical. Let's push this further. In what other ways does the Wigner function mimic a true probability distribution?

First, and most importantly, it gives the correct "shadows." If you have a 3D object, its shadow on one wall tells you its silhouette from that angle, and its shadow on another wall tells you another silhouette. The Wigner function behaves similarly. If you don't care about momentum and just want to know the probability of finding the particle at position $x$, you can simply add up (integrate) the values of $W(x,p)$ over all possible momenta $p$. The result is, miraculously, $|\psi(x)|^2$, the correct quantum mechanical probability density for position. Likewise, integrating over all positions $x$ gives the correct momentum probability density $|\phi(p)|^2$. Even though $W(x,p)$ itself might be strange, its projections onto the fundamental axes of position and momentum are perfectly well-behaved and physically correct.

Second, it transforms under changes of perspective just as you'd expect. If you re-analyze the system from the viewpoint of someone in a moving car (a Galilean boost with velocity $v$), the Wigner function simply shifts in momentum: the new distribution is just the old one evaluated at $p-mv$. This is precisely what classical intuition would tell you. If you "run the movie backwards" (apply a time-reversal operation), the position remains the same, but the momentum flips its sign: $W_T(x,p) = W(x,-p)$. Again, this is in perfect accord with classical physics.

Finally, for many important systems, its time evolution is classical. For a free particle or a harmonic oscillator—systems whose energy is a quadratic function of position and momentum—the Wigner function evolves in time exactly as a cloud of classical particles would. Each point of the distribution flows along the classical trajectory. For a free particle, this leads to a "shearing" effect in phase space, as faster parts of the distribution move further. For a harmonic oscillator, the entire distribution rotates rigidly around the origin like a spinning platter, with an angular frequency equal to the oscillator's own frequency $\omega$.

Up to this point, you might be wondering what the fuss is about. The Wigner function seems to be a perfectly good, and remarkably intuitive, phase-space distribution. But now we come to the "quasi" part, the feature that makes the Wigner function a profound window into the soul of quantum mechanics.

​​The Wigner function can be negative.​​

A negative probability is, of course, a physical absurdity. You can't have a -10% chance of finding your keys. The fact that the Wigner function can dip below zero is the unambiguous, undeniable signature that it is not a true probability distribution. It is a more exotic object, and these regions of negativity are the smoking gun of quantum interference.

To see this, we don't have to look far. Let's move from the ground state ($n=0$) of the harmonic oscillator to its first excited state ($n=1$). Calculating the Wigner function for this state reveals something striking. It's no longer a simple mound. Instead, it looks like a doughnut, with a positive ring surrounding a central region. But the center isn't just empty—it has a negative value!. Specifically, this negative region corresponds to points in phase space where the classical energy, $H = \frac{p^2}{2m} + \frac{1}{2}m\omega^2 x^2$, is less than $\frac{\hbar\omega}{4}$.

This expression shows it all: whenever the classical energy $H$ is less than $\hbar\omega/4$, the term in the parenthesis is negative, making the whole function negative. This negativity is not a mathematical quirk; it's a fundamental feature of quantum states that are not simple Gaussians. We can even calculate the total "volume" of this negative part, giving a concrete number that quantifies the "non-classicality" of the state.

This raises a deep question. If the quantum world is filled with these bizarre negative quasi-probabilities, why is our everyday macroscopic world so stubbornly classical and positive? The Wigner function provides two beautiful answers to this.

The first answer is related to the act of measurement. Any real-world measurement has finite precision; we can never measure position and momentum with perfect sharpness. This is not just a technical limitation, but a fundamental consequence of the uncertainty principle. What happens if we "blur" our vision of the Wigner function, mimicking a fuzzy measurement? We can do this mathematically by "smoothing" the function, averaging it with a small Gaussian blur.

The result is astounding. If the blur is very fine, the negative regions persist. But if we make the blur just big enough—specifically, with a characteristic size in phase space on the order of Planck's constant $\hbar$—the negativity is guaranteed to wash away, leaving a perfectly positive, classical-like probability distribution. The minimum amount of smoothing needed corresponds to the limit set by the uncertainty principle itself. In a sense, the uncertainty principle acts as a cosmic censor, hiding the most jarring quantum features from any coarse-grained classical viewpoint.

The second answer comes from looking at high-energy states. What does the Wigner function look like for our harmonic oscillator in a state with a very large quantum number, $n \gg 1$? As described in, the function develops extremely rapid oscillations between large positive and large negative values. These quantum interference fringes are concentrated in a narrow band around the classical ellipse corresponding to the energy $E_n = (n+1/2)\hbar\omega$. The "wavelength" of these wiggles in phase space is of order $\hbar$.

To a macroscopic observer, whose instruments are far too crude to resolve such microscopic details, these rapid oscillations simply average out to zero. The only thing that survives this averaging is the net positive ridge located right on the classical trajectory. Thus, in the macroscopic limit, the oscillating quantum Wigner function effectively converges to the classical microcanonical distribution: a distribution that is zero everywhere except on the path that a classical particle with that energy would take. This is the correspondence principle, visualized with breathtaking clarity. The classical world of definite trajectories emerges from the quantum world of interference, not by ignoring quantum mechanics, but by averaging over its intricate, oscillatory structure.

Having acquainted ourselves with the curious nature of the Wigner function—this strange beast that lives in phase space, part probability, part wave—we might be tempted to file it away as a mathematical oddity of quantum theory. But to do so would be to miss the point entirely. The true power and beauty of a physical idea are measured not by its abstract elegance, but by its ability to connect disparate phenomena, to simplify the complex, and to grant us a new and more powerful way of seeing the world. The Wigner function is a tool of precisely this caliber. It is a kind of universal translator, a Rosetta Stone for wave physics, allowing us to read the story of a wave—be it a quantum electron, a classical light beam, or a seismic tremor in a star—in the unified language of phase space.

In this chapter, we will embark on a journey through the vast landscape of its applications. We will see how this single concept illuminates problems in quantum mechanics, revolutionizes our understanding of optical systems, provides a new language for signal processing, and even gives us tools to probe the heart of a material or the interior of the sun. The journey will reveal that the strange rules of phase-space dynamics are not an esoteric feature of quantum mechanics alone, but a deep and unifying principle of all wave phenomena.

Our journey begins in the Wigner function's native land: quantum mechanics. Here, its ability to represent a quantum state in phase space is not just a novelty; it is a profound interpretative and analytical tool. Quantum mechanics is famous for its weirdness, for states that are somehow here and there at the same time. How can we visualize such a thing? The Wigner function gives us a picture.

Consider, for instance, a state that is a superposition of two distinct quantum states, like the vacuum and a coherent state—a so-called "Schrödinger cat" state. If we were to plot the Wigner function for this superposition, we would see two distinct, positive "blobs" in phase space, corresponding to the two individual states. But between them, something marvelous appears: a series of rapidly oscillating fringes, which dip into negative values. These negative regions are the smoking gun of quantum mechanics. They are a direct visualization of quantum interference, a feature with no classical analogue. The Wigner function doesn't just tell us that interference is happening; it shows us where in phase space the quantum "weirdness" is concentrated.

This ability to dissect quantum phenomena extends to dynamic processes as well. Imagine a quantum particle scattering from a potential barrier. The wavefunction in the reflection region is a superposition of the incoming wave and the reflected wave. The Wigner function for this situation beautifully decomposes the state. It shows a distribution for the incident particle moving with momentum $+p$, another for the reflected particle moving with momentum $-p$, and a third, purely quantum feature: an interference term. This interference term is localized at zero momentum ($p=0$) but oscillates in space. It is the phase-space signature of the standing wave pattern created by the colliding incident and reflected waves. The Wigner function allows us to cleanly separate the classical-like parts of the motion from the purely quantum interference effects.

Perhaps the most stunning testament to the Wigner function's power is its seamless applicability to classical wave optics. It turns out that the mathematics governing a quantum particle's wavefunction and a paraxial light beam's complex amplitude are formally identical. This means we can import our entire phase-space toolkit from quantum mechanics and apply it to light. In this new context, position $x$ remains position, but momentum $p$ is replaced by the propagation angle or transverse spatial frequency $k_x$. The Wigner function now describes how the light's energy is distributed across positions and angles simultaneously.

The true magic happens when we consider how light propagates. In the traditional wave picture, calculating how a beam changes as it passes through a series of lenses, apertures, and free space requires wrestling with cumbersome diffraction integrals. The phase-space picture offers a breathtakingly simple alternative. For a vast class of optical systems—so-called first-order or ABCD systems—the complex integral transforms are replaced by simple geometric transformations in phase space.

For example, the propagation of a beam through empty space, a process described by the Fresnel diffraction integral, corresponds to a simple "shear" of its Wigner function in phase space. A ray at a given position $x$ with an angle (momentum) $k_x$ simply moves to a new position $x + L(k_x/k)$ after propagating a distance $L$, while its angle remains the same—exactly as a classical light ray would! Similarly, passing a beam through an ideal thin lens corresponds to a different shear, this time in the momentum coordinate. The entire propagation through a complex optical system of lenses and spaces becomes a sequence of simple matrix multiplications acting on the phase-space coordinates. The formidable physics of wave diffraction is tamed into the simple elegance of ray tracing in phase space.

This framework is not just elegant; it is immensely practical. It allows us to characterize the "quality" of a laser beam. A partially coherent beam, for instance, can be described by a Gaussian Schell-model, and its Wigner function neatly encapsulates both its spatial size and its angular spread. From the second moments of this distribution, one can directly calculate the beam parameter product and the all-important beam quality factor, $M^2$, a standard metric in the laser industry. The Wigner function also provides a powerful diagnostic tool. If a lens suffers from imperfections like spherical aberration, the WDF of a beam passing through it will be distorted from its ideal shape in a characteristic way, providing a direct visual and quantitative measure of the aberration's impact.

The concept of phase space is not limited to position and momentum. Any pair of conjugate variables, connected by a Fourier transform, can form a phase space. A powerful example is the relationship between time and frequency. For any time-dependent signal, like a radio wave or an ultrashort laser pulse, we can define a time-frequency Wigner function. This function, $W(t, \omega)$, tells us the frequency content of the signal at every instant in time.

Consider a linearly chirped Gaussian pulse—a pulse whose frequency changes linearly with time, like a sliding musical note. If we calculate its Wigner function, we find that it forms an elliptical blob in the time-frequency plane. Crucially, this ellipse is tilted. The tilt angle is directly proportional to the chirp rate $\beta$. This provides a beautiful, intuitive picture: as time $t$ increases, the center of the frequency distribution $\omega$ shifts. The Wigner function makes the abstract concept of a chirp visible as a simple geometric feature. This tool is indispensable in fields like radar, sonar, and ultrafast optics, where understanding the time-varying frequency content of signals is paramount.

The reach of the Wigner function extends even further, into the experimental and observational sciences, providing new ways to see the invisible and understand the immense.

In materials science, a technique called 4D Scanning Transmission Electron Microscopy (4D-STEM) allows scientists to map the fields inside materials at the atomic scale. A focused beam of electrons is scanned across a sample, and for each position, a full diffraction pattern is recorded. How can we make sense of this colossal dataset? The Wigner function provides the key. It turns out that the center of mass of the measured diffraction pattern is directly proportional to the average momentum transfer the electrons experience. Using the Wigner function formalism, one can show a remarkably direct relationship: the shift in this center of mass is proportional to the electric field at that point in the sample. By simply tracking how the diffraction pattern shifts as the probe is scanned, scientists can create a direct map of the electric fields that hold atoms together. We are literally watching the force fields of nature at work.

From the infinitesimally small to the astronomically large, the Wigner function remains a faithful guide. In helioseismology, scientists study the waves that constantly ripple through the Sun's interior to learn about its structure. The Sun is an inhomogeneous medium, where the wave speed changes with depth. A full wave simulation is incredibly complex. However, using the WKB approximation (a method for finding approximate solutions to wave equations in slowly varying media), we can analyze the situation with the Wigner function. The result is profound: in a weakly varying medium, the Wigner function of the wave becomes sharply peaked along the path that a classical ray would take. It shows us, in phase space, that the wave's energy is concentrated along the geometric ray path. It provides a rigorous connection between the full, complicated wave picture and the more intuitive ray-tracing picture, allowing astrophysicists to interpret solar oscillations in terms of seismic rays bending through the Sun's fiery plasma.

From quantum cat states to laser beams, from chirped pulses to atomic fields and stellar vibrations, the Wigner function provides a common thread. It teaches us that to truly understand a wave, we must look at it in phase space. There, hidden complexities often resolve into simple geometric pictures, and deep connections between seemingly unrelated fields of physics are laid bare. It is a testament to the unifying power of physical law, and a beautiful example of how a shift in perspective can transform our understanding of the universe.