Gaussian States and Wigner Functions

Gaussian states are the tractable core of phase-space quantum mechanics. Their Wigner functions are Gaussian functions on phase space, so the state is controlled by its first moments and covariance matrix. This makes Gaussian states the natural bridge between wave packets, harmonic-oscillator ground states, coherent states, squeezed states, semiclassical dynamics, and continuous-variable quantum information.

The Wigner Function page defines the transform. This page specializes to the Gaussian sector and explains why it is closed under quadratic Hamiltonian evolution.

Continuous-Variable Quantum Computation uses this moment and symplectic framework to audit Gaussian computational transformations, finite-squeezed cluster states, continuous outcomes, and departures from Gaussian closure. The phase-space definitions and covariance derivations remain here.

Continuous-Variable Platforms owns the corresponding optical and microwave hardware, non-Gaussian resource boundary, noise ledger, and experimental evidence.

Phase-Space Variables and Moments

For one degree of freedom, collect canonical variables into the phase-space vector

$$\mathbf z = \begin{pmatrix} x\\ p \end{pmatrix}.$$

The first moment is

$$\mathbf d = \langle\mathbf z\rangle = \begin{pmatrix} \langle x\rangle\\ \langle p\rangle \end{pmatrix}.$$

The covariance matrix is

$$V_{ij} = \frac12 \left\langle \Delta z_i\Delta z_j+\Delta z_j\Delta z_i \right\rangle, \qquad \Delta z_i=z_i-d_i.$$

Explicitly,

$$V = \begin{pmatrix} (\Delta x)^2 & C_{xp}\\ C_{xp} & (\Delta p)^2 \end{pmatrix},$$

where

$$C_{xp} = \frac12 \langle \Delta x\,\Delta p+\Delta p\,\Delta x \rangle.$$

The covariance matrix records widths and phase-space tilt. For a Gaussian Wigner function, no higher moments are independent data.

Gaussian Wigner Functions

With the Wigner normalization used in this volume,

$$\int dx\,dp\,W(x,p)=1,$$

a one-mode Gaussian Wigner function is

$$W(\mathbf z) = \frac{1}{2\pi\sqrt{\det V}} \exp\left[ -\frac12 (\mathbf z-\mathbf d)^T V^{-1} (\mathbf z-\mathbf d) \right].$$

For a diagonal covariance matrix this becomes

$$W(x,p) = \frac{1}{2\pi\sigma_x\sigma_p} \exp\left[ -\frac{(x-x_0)^2}{2\sigma_x^2} -\frac{(p-p_0)^2}{2\sigma_p^2} \right],$$

where $\mathbf d=(x_0,p_0)^T$ and $V=\operatorname{diag}(\sigma_x^2,\sigma_p^2)$.

This formula looks like an ordinary classical Gaussian probability density. It is still a Wigner function: it represents a quantum density operator only when its covariance matrix satisfies the uncertainty constraint.

Uncertainty Constraint

Let

$$\Omega = \begin{pmatrix} 0 & 1\\ -1 & 0 \end{pmatrix}$$

be the symplectic matrix. A physical one-mode Gaussian covariance matrix satisfies

$$V+\frac{i\hbar}{2}\Omega \geq 0.$$

For one degree of freedom this is equivalent to

$$\det V \geq \frac{\hbar^2}{4}.$$

In components,

$$(\Delta x)^2(\Delta p)^2-C_{xp}^2 \geq \frac{\hbar^2}{4}.$$

Pure one-mode Gaussian states saturate this inequality:

$$\det V = \frac{\hbar^2}{4}.$$

Mixed Gaussian states have larger covariance determinant. Thus a broad positive Gaussian on phase space is not automatically a pure minimum-uncertainty state.

Gaussian Wave Packets

A minimum-uncertainty wave packet with no $xp$ correlation has

$$V = \begin{pmatrix} \sigma_x^2 & 0\\ 0 & \sigma_p^2 \end{pmatrix}, \qquad \sigma_x\sigma_p=\frac{\hbar}{2}.$$

Its Wigner function is

$$W(x,p) = \frac{1}{\pi\hbar} \exp\left[ -\frac{(x-x_0)^2}{2\sigma_x^2} -\frac{(p-p_0)^2}{2\sigma_p^2} \right].$$

This is the phase-space version of the Gaussian wave packets developed in Gaussian Wave Packets. The center $(x_0,p_0)$ gives the classical-looking location of the packet, while the covariance matrix gives its quantum spread.

If $C_{xp}\ne0$, the Gaussian ellipse is tilted in phase space. Such correlations appear naturally during free evolution: a wave packet spreads because position becomes correlated with momentum.

Harmonic-Oscillator Ground State

For the harmonic oscillator

$$H = \frac{p^2}{2m} + \frac12m\omega^2x^2,$$

the ground state has covariance matrix

$$V_0 = \begin{pmatrix} \dfrac{\hbar}{2m\omega} & 0\\ 0 & \dfrac{m\hbar\omega}{2} \end{pmatrix}.$$

The Wigner function is

$$W_0(x,p) = \frac{1}{\pi\hbar} \exp\left[ -\frac{m\omega x^2}{\hbar} -\frac{p^2}{m\hbar\omega} \right].$$

Its covariance determinant is

$$\det V_0 = \frac{\hbar^2}{4},$$

so it is a pure minimum-uncertainty Gaussian. Its phase-space ellipse is aligned with the oscillator energy contours after the natural scaling of $x$ and $p$.

Coherent States

A coherent state is a displaced oscillator ground state. In Wigner phase space, displacement changes the first moment but not the covariance:

$$\mathbf d = \begin{pmatrix} x_0\\ p_0 \end{pmatrix}, \qquad V=V_0.$$

Therefore

$$W_{\rm coh}(x,p) = \frac{1}{\pi\hbar} \exp\left[ -\frac{m\omega (x-x_0)^2}{\hbar} -\frac{(p-p_0)^2}{m\hbar\omega} \right].$$

Under harmonic-oscillator evolution, the center follows the classical orbit and the covariance returns to itself. This is the phase-space reason coherent states are classical-like oscillator states. The Hilbert-space construction is canonical in Coherent States; this page records the Wigner-function viewpoint.

Squeezed States Preview

A squeezed state changes the covariance while preserving the uncertainty product for a pure Gaussian. In a simple oscillator-aligned convention,

$$(\Delta x)^2 = \frac{\hbar}{2m\omega}e^{2r}, \qquad (\Delta p)^2 = \frac{m\hbar\omega}{2}e^{-2r},$$

with $C_{xp}=0$. The determinant remains

$$(\Delta x)^2(\Delta p)^2 = \frac{\hbar^2}{4}.$$

For $r\gt0$, the Wigner ellipse is stretched in $x$ and squeezed in $p$. For $r\lt0$, the roles reverse. More general squeezed states also rotate the ellipse, producing nonzero $C_{xp}$.

Squeezing does not violate the uncertainty principle. It redistributes uncertainty between conjugate quadratures while keeping the symplectic area bounded below by $\hbar/2$.

Evolution Under Quadratic Hamiltonians

Gaussian states remain Gaussian under Hamiltonians that are at most quadratic in $x$ and $p$. The reason is that Wigner evolution is exactly classical for quadratic Hamiltonians: higher Moyal corrections vanish.

Let the classical phase-space flow be linear,

$$\mathbf z(t) = S(t)\mathbf z(0),$$

where $S(t)$ is symplectic:

$$S(t)\Omega S(t)^T=\Omega.$$

Then a Gaussian evolves by

$$\mathbf d(t) = S(t)\mathbf d(0), \qquad V(t) = S(t)V(0)S(t)^T.$$

For a free particle,

$$x(t)=x(0)+\frac{t}{m}p(0), \qquad p(t)=p(0).$$

The covariance evolves as

$$\begin{aligned} V_{xx}(t) &= V_{xx}(0) + \frac{2t}{m}V_{xp}(0) + \frac{t^2}{m^2}V_{pp}(0), \\ V_{xp}(t) &= V_{xp}(0) + \frac{t}{m}V_{pp}(0), \\ V_{pp}(t) &= V_{pp}(0). \end{aligned}$$

This is wave-packet spreading in covariance-matrix language. The Wigner function shears in phase space rather than diffusing like a classical noisy distribution.

Why Gaussian States Are Special

Gaussian Wigner functions are special for three reasons:

• They are determined entirely by first and second moments.
• They remain Gaussian under quadratic Hamiltonians and linear canonical transformations.
• Their Wigner functions are nonnegative, so many calculations resemble classical Gaussian probability theory while still obeying quantum uncertainty.

The last point must be handled carefully. Positivity of a Wigner function is not the same as classicality in every operational sense, and mixed non-Gaussian states can also have nonnegative Wigner functions. Gaussian states are tractable because the phase-space calculus closes, not because quantum mechanics has disappeared.

Common Mistakes

• Treating a positive Gaussian Wigner function as proof that the state is classical in all respects.
• Forgetting the covariance constraint $\det V\geq\hbar^2/4$.
• Confusing displacement, which changes first moments, with squeezing, which changes covariance.
• Assuming all positive Wigner functions are Gaussian.
• Applying quadratic-flow covariance formulas to anharmonic Hamiltonians without including Moyal corrections.
• Forgetting that mixed Gaussian states need not saturate the uncertainty bound.

Cross-Links

• Optical Phase-Space Distributions for $P$, Wigner, and $Q$ representations of coherent, thermal, and squeezed light.
• Squeezed Light for optical generation, calibrated quadrature measurements, and realistic loss limits.
• Wigner Function
• Coherent States in Phase Space
• Phase-Space Dynamics
• Weyl Transform
• Moyal Bracket
• Gaussian Wave Packets
• Coherent States
• Quantum Harmonic Oscillator
• Formula Sheet
• Wigner Function Notebook

References

• E. Wigner, “On the Quantum Correction For Thermodynamic Equilibrium,” Physical Review 40, 749-759, 1932.
• M. Hillery, R. F. O’Connell, M. O. Scully, and E. P. Wigner, “Distribution functions in physics: Fundamentals,” Physics Reports 106, 121-167, 1984.
• W. P. Schleich, Quantum Optics in Phase Space, Wiley-VCH, 2001.
• A. Serafini, Quantum Continuous Variables: A Primer of Theoretical Methods, CRC Press, 2017.
• C. Weedbrook et al., “Gaussian quantum information,” Reviews of Modern Physics 84, 621-669, 2012.
• C. K. Zachos, D. B. Fairlie, and T. L. Curtright, eds., Quantum Mechanics in Phase Space, World Scientific, 2005.

Exercises

1. Normalize the one-mode Gaussian Wigner function.

Solution

For positive definite $V$,

$$\int d^2z\, \exp\left[ -\frac12 (\mathbf z-\mathbf d)^T V^{-1} (\mathbf z-\mathbf d) \right] = 2\pi\sqrt{\det V}.$$

Therefore

$$\int d^2z\,W(\mathbf z) = \frac{1}{2\pi\sqrt{\det V}} 2\pi\sqrt{\det V} = 1.$$

2. Check the uncertainty determinant for the harmonic-oscillator ground state.

Solution

For the ground state,

$$V_0 = \begin{pmatrix} \dfrac{\hbar}{2m\omega} & 0\\ 0 & \dfrac{m\hbar\omega}{2} \end{pmatrix}.$$

Thus

$$\det V_0 = \frac{\hbar}{2m\omega} \frac{m\hbar\omega}{2} = \frac{\hbar^2}{4}.$$

The pure Gaussian saturates the one-mode uncertainty bound.

3. Show that squeezing preserves the covariance determinant in the aligned convention.

Solution

With

$$(\Delta x)^2 = \frac{\hbar}{2m\omega}e^{2r}, \qquad (\Delta p)^2 = \frac{m\hbar\omega}{2}e^{-2r},$$

and $C_{xp}=0$, the determinant is

$$\det V = (\Delta x)^2(\Delta p)^2 = \frac{\hbar^2}{4} e^{2r}e^{-2r} = \frac{\hbar^2}{4}.$$

Squeezing changes the ellipse shape while preserving its symplectic area for a pure Gaussian.

4. Derive the free-particle covariance shear.

Solution

For a free particle,

$$\Delta x(t) = \Delta x(0)+\frac{t}{m}\Delta p(0), \qquad \Delta p(t) = \Delta p(0).$$

Then

$$\begin{aligned} V_{xx}(t) &= \langle\Delta x(t)^2\rangle \\ &= V_{xx}(0) + \frac{2t}{m}V_{xp}(0) + \frac{t^2}{m^2}V_{pp}(0), \end{aligned}$$

and

$$V_{xp}(t) = \frac12 \langle \Delta x(t)\Delta p(t) + \Delta p(t)\Delta x(t) \rangle = V_{xp}(0)+\frac{t}{m}V_{pp}(0).$$

Finally, $V_{pp}(t)=V_{pp}(0)$ because $\Delta p$ is constant.