---
title: 'Wigner Current: Quantum Phase-Space Flow'
url: https://www.emergentmind.com/topics/wigner-current
type: topic
---

# Wigner Current: Quantum Phase-Space Flow

Searching arXiv for recent and foundational papers on Wigner current.
Wigner current is the phase-space current associated with the Wigner function \(W\), defined by the exact continuity equation
\[
\partial_t W(x,p,t)+\partial_x J_x(x,p,t)+\partial_p J_p(x,p,t)=0.
\]
In Wigner’s phase-space formulation of quantum mechanics, it replaces point-wise trajectories by a conserved quasi-probability flow. Its vector field, fieldlines, and stagnation sets provide a quantum phase portrait for both closed and open dynamics, while also exposing structures that have no classical analogue, including flow reversal in regions of negative \(W\), discrete topological defects, and nontrivial limits as \(\hbar\to0\) or anharmonicity vanishes [1411.3511, 2307.16510].

## 1. Definition and continuity structure

For a one-dimensional system with density operator \(\hat\rho(t)\), one convenient definition of the Wigner function is
\[
W(x,p,t)=\frac{1}{2\pi\hbar}\int_{-\infty}^{\infty}dy\;e^{ipy/\hbar}\,
\langle x-y/2|\hat\rho(t)|x+y/2\rangle.
\]
If the Hamiltonian is \(H(x,p)=p^2/(2m)+V(x)\), the von Neumann equation maps to the Moyal equation
\[
\partial_t W(x,p,t)=\{\!\!\{H,W\}\!\!\}(x,p,t),
\]
with
\[
\{\!\!\{f,g\}\!\!\}
=
\frac{1}{i\hbar}(f\star g-g\star f)
=
\frac{2}{\hbar}f(x,p)\sin\!\Bigl[\frac{\hbar}{2}
\bigl(\overleftarrow{\partial_x}\overrightarrow{\partial_p}
-\overleftarrow{\partial_p}\overrightarrow{\partial_x}\bigr)\Bigr]g(x,p).
\]
Reorganizing the Moyal expansion yields the continuity equation above, with the closed-system current
\[
J_x(x,p,t)=\frac{p}{m}W(x,p,t),
\]
\[
J_p(x,p,t)=
-\sum_{n=0}^{\infty}\frac{(-1)^n}{(2n+1)!}\Bigl(\frac{\hbar}{2}\Bigr)^{2n}
\partial_x^{2n+1}V(x)\,\partial_p^{2n}W(x,p,t).
\]
An equivalent form isolates the classical force term,
\[
J_p=-\partial_xV\,W
+\sum_{n=1}^{\infty}\frac{(-1)^n}{(2n+1)!}\Bigl(\frac{\hbar}{2}\Bigr)^{2n}
\partial_x^{2n+1}V\,\partial_p^{2n}W.
\]
Physically, \(J\) is the flow of quasi-probability in phase space: quantum dynamics remains conservative in the variables \((W,J)\), even though \(W\) is not a positive probability density [2307.16510].

## 2. Closed-system flow and its relation to classical mechanics

The classical Liouville current for \(H=p^2/(2m)+V(x)\) is
\[
J_x^{\rm cl}=\frac{p}{m}W,\qquad
J_p^{\rm cl}=-V'(x)W.
\]
The quantum current agrees with this at leading order, but the higher even powers of \(\hbar\) in \(J_p\) are purely quantum corrections. These terms encode interference and tunneling effects, and they give rise to non-Liouvillian flow, in the sense that \(\nabla\!\cdot J\neq0\) locally; they also produce regions where \(W<0\) and the current reverses, forming the “whorls” discussed in phase-space analyses [1411.3511, 2307.16510].

The harmonic oscillator is the special case in which the current reduces to
\[
J_x=W\,p,\qquad J_p=-W\,x,
\]
so the fieldlines are circular. For harmonic eigenstates, the Wigner function has alternating positive and negative rings, and the zero-rings of \(W\) coincide with stagnation lines of \(J\). This exact coincidence is exceptional. As soon as anharmonicity is introduced, the zeros of \(J_x\) and \(J_p\) shift differently, and the stagnation structure becomes discrete rather than line-like [1411.3511].

A central result for weakly anharmonic systems is that neither the limit \(\hbar\to0\) nor the limit of vanishing anharmonicity yields pointwise convergence of the Wigner current to the classical flow or to the harmonic-oscillator circular flow. This rejects a common intuition that the current field should collapse smoothly, point by point, onto the classical portrait [1411.3511].

## 3. Stagnation points, winding numbers, and topological constraints

A stagnation point is a phase-space point at which \(J=0\). In the harmonic oscillator, continuous stagnation lines appear because \(J_x\) and \(J_p\) both vanish wherever \(W=0\). In weakly anharmonic systems, those lines split into isolated crossings of the \(J_x=0\) and \(J_p=0\) contours. The resulting discrete defects organize the current fieldline pattern [1411.3511].

Each stagnation point carries an integer Poincaré–Hopf index, or winding number,
\[
\omega=\frac{1}{2\pi}\oint_{\mathcal L} d\phi,
\]
where \(\phi(x,p)=\arg(J_x+iJ_p)\) and \(\mathcal L\) is a small counter-clockwise loop encircling the defect. In the weakly anharmonic analyses, only \(\omega=+1\) and \(\omega=-1\) occur, corresponding respectively to vortices and saddles. Time evolution can create or annihilate pairs of opposite charge, but the total sum of \(\omega\) over the plane is conserved [1411.3511].

This topological viewpoint extends naturally to time-dependent superpositions. For the ground–first-excited superposition
\[
\psi=\cos\theta\,\psi_0+\sin\theta\,e^{-i\Delta E\,t/\hbar}\psi_1,
\]
the zero-circle of \(W\) remains circular while its center spirals, and its stagnation points move rigidly “like markers on a Ferris-wheel.” Under slow Rabi driving, the pattern expands and contracts; opposite-charge pairs may be created or annihilated, but the total index remains invariant [1411.3511].

In an open-system optical realization, the ideal squeezed-vacuum dynamics of a degenerate optical parametric oscillator has a central stagnation point at the origin with topological charge \(Q=-1\). Experimental work showed that this charge remains invariant under both weakly and strongly decohering conditions, provided no pair-creation or pair-annihilation event transfers a defect through the origin. In that sense the central charge is topologically protected [2111.08285].

## 4. Open-system extensions and Lindblad structure

For a single bosonic mode weakly coupled to a thermal bath, the phase-space continuity equation acquires environmental terms,
\[
\partial_tW+\nabla\!\cdot J+\nabla\!\cdot J_{\rm env}=0,
\qquad
J_{\rm env}=J_{\rm damp}+J_{\rm diff}.
\]
In units where \(m=\hbar=\omega_0=1\), the environmental currents are
\[
J_{\rm damp}(x,p)=-\frac{\gamma}{2}\,W(x,p)\,(x,p)^T,
\qquad
J_{\rm diff}(x,p)=
-\frac{\gamma}{2}\Bigl(\bar n+\frac12\Bigr)(\partial_xW,\partial_pW)^T.
\]
The damping part is a friction-type drift back toward the origin; the diffusion part spreads \(W\) out. Together they enforce relaxation to a thermal Gaussian [2307.16510].

In dimensionful variables, the same split can be written as
\[
J_{\rm damp}(q,p)= -\frac{\gamma}{2}(qW,pW),
\]
\[
J_{\rm diff}(q,p)=
-\frac{\gamma\hbar}{2}\Bigl(\bar n+\frac12\Bigr)
\bigl((m\omega_0)^{-1}\partial_qW,\;m\omega_0\,\partial_pW\bigr).
\]
This decomposition makes it possible to analyze decoherence geometrically. In particular, for a region where \(W<0\), Gauss’s theorem yields an exact flux law for the integrated negativity. Diffusion always points inward on a negative region and therefore shrinks it, whereas the system current can point outward on parts of the boundary and can partially counteract decoherence. In the harmonic case negativity always decays; in the driven Duffing case, the leading quantum correction can pump probability so as to enlarge the negative island during part of a drive cycle, and a stabilization criterion can be written in terms of the boundary fluxes [1703.04844].

A further structural result is that requiring all dynamical terms to be expressible as a total divergence in phase space forces the dissipative coupling into Lindblad form. In the star-product representation, only the combinations
\[
2a\star W\star a^\dagger-a^\dagger\star a\star W-W\star a^\dagger\star a
\]
and
\[
2a^\dagger\star W\star a-a\star a^\dagger\star W-W\star a\star a^\dagger
\]
produce the needed divergence terms. The resulting phase-space dissipator is exactly the Wigner-space image of the operator-space Lindblad theorem, and hence preserves trace, hermiticity, and complete positivity [2307.16510].

## 5. Canonical models and generalized phase spaces

Weakly anharmonic one-dimensional bound systems admit a threefold classification. For
\[
V(x)=\frac12x^2+\alpha_\nu x^\nu+\cdots
\]
with small \(\alpha_\nu\), one obtains hard potentials (\(\nu\) even, \(\alpha_\nu>0\)), soft potentials (\(\nu\) even, \(\alpha_\nu<0\)), and odd potentials (\(\nu\) odd, \(\alpha_\nu<0\)). In hard potentials, the zero-circles of \(W\) are squeezed along \(x\) and stretched along \(p\); in soft potentials the deformation is reversed; odd potentials behave as hard on one side and soft on the other, with a cut-and-reconnect structure near \(p=0\). For the first excited state, the hard quartic and soft quartic cases yield nine stagnation points per ring, while the odd cubic case yields seven because part of the axis-crossing structure avoids \(p=0\) through reconnection [1411.3511].

For the Kerr oscillator,
\[
H_\Lambda(x,p)=\frac12(p^2+x^2)+\Lambda^2\Bigl[\frac12(p^2+x^2)\Bigr]^2,
\]
Oliva and Steuernagel derived
\[
J(x,p)=
\begin{pmatrix}p\\-x\end{pmatrix}
\Bigl[1+\Lambda^2\Bigl(x^2+p^2-\frac14\Delta\Bigr)\Bigr]W(x,p),
\qquad
\Delta=\partial_x^2+\partial_p^2.
\]
Because \(J\cdot(x,p)=0\), the current is everywhere tangent to circles \(r=\mathrm{const}\), and the probability on each circle is constant in time. Decomposing \(J=j+J^Q\), with
\[
j=W\,v,\qquad v=(p,-x)(1+\Lambda^2r^2),
\]
\[
J^Q=-\frac{\Lambda^2}{4}(p,-x)\Delta W,
\]
isolates a purely quantum term that suppresses classical shear. The global measure
\[
\Pi(t)=\iint dx\,dp\,W\,\partial_r(-\nabla\times J^Q)
\]
initially grows and then settles to a negative plateau for \(\Lambda^2>0\), signaling effective “viscosity” and the enforcement of the Zurek scale. Deviations of \(\Pi(t)\) from its steady behavior pick out fractional revivals and special nonclassical states [1811.02952].

Wigner current has also been formulated in several nonstandard geometries. For variable beam splitters with bilinear mode-mixing Hamiltonian
\[
\hat H_M=\frac{\pi}{2}(\hat x_a\hat p_b-\hat p_a\hat x_b),
\]
the reduced current of one outgoing mode after tracing out the other is expressed through cross-moments of the two-mode Wigner function. In this setting, coherent–Fock mixing produces circulating currents around toroidal negative regions, the Hong–Ou–Mandel configuration yields strictly radial currents with sink–source structure inside the negative “hole,” and the ratio \(w_x=J_x/W\) becomes singular wherever \(W=0\), signaling volume-nonconserving quantum dynamics [2308.06706].

For spin systems, the SU(2) Wigner function on the sphere satisfies
\[
\partial_tW(\theta,\phi|t)+\nabla\!\cdot J(\theta,\phi|t)=0,
\]
and for Kerr evolution \(\hat H=\chi \hat S_z^2\) one finds
\[
J_\theta=0,\qquad
J_\phi=(\chi/\epsilon)\sin\theta\cos\theta\,\hat\Gamma(\theta,\mathcal L^2)\,W,
\qquad \epsilon=\frac{1}{2S+1}.
\]
In the classical limit \(\epsilon\to0\), \(\hat\Gamma\to1\) and the flow becomes a twist around the \(z\)-axis, \(\dot\phi\propto\cos\theta\). Quantum mechanically, additional stagnation loops, elliptic points, and hyperbolic saddles appear; these are absent in the classical current and arise from the full deformation operator \(\hat\Gamma\) [1812.03097].

In multidimensional quantum billiards, the Wigner current satisfies
\[
\partial_tW+\nabla_x\!\cdot j_W^x+\nabla_p\!\cdot j_W^p=0,
\qquad
j_W^x(x,p)=\frac{p}{m}W(x,p),
\]
while the momentum component becomes a surface integral over the boundary manifold in chord space,
\[
j_W^p(x,p,t)=
-\frac{i}{(2\pi\hbar)^N m}
\oint_{\omega(x,y)}e^{-ip\cdot y/\hbar}
\nabla_x\!\bigl[\psi^*(x-y/2)\psi(x+y/2)\bigr]\,
d^{N-1}y.
\]
An equivalent description uses the boundary potential \(\delta'_S(x)/(2m)\). Here the wall-induced quantum effects are carried entirely by \(j_W^p\); at finite \(\hbar\), the surface integral retains nonlocal contributions even in regions that are otherwise free [2408.14164].

## 6. Experimental reconstruction and interpretive significance

The first experimental reconstruction of Wigner current was reported for squeezed-vacuum states generated by a degenerate optical parametric oscillator. In the rotating frame, the ideal system Hamiltonian is
\[
\hat H=\frac{i\hbar}{2}\chi^{(2)}(\alpha^*\hat a^2-\alpha\hat a^{\dagger2}),
\]
with associated system current
\[
\mathbf J_{\rm sys}(x,p)
=
\chi^{(2)}|\alpha|
\begin{pmatrix}
xW\cos\theta+\frac{p}{\omega_0}W\sin\theta\\
\omega_0xW\sin\theta-pW\cos\theta
\end{pmatrix}.
\]
In the squeezed-vacuum phase \(\theta=\pm\pi/2\), this reduces to hyperbolic flow that pushes the Wigner density along the squeezed and anti-squeezed quadratures. The environment adds the damping and diffusion currents discussed above, yielding the experimentally observed “push-and-pull” picture: damping pulls inward, while diffusion pushes outward and enforces Heisenberg’s uncertainty relations even at zero temperature [2111.08285].

The reported setup used a bow-tie cavity with a PPKTP crystal pumped at \(532\,\mathrm{nm}\) to generate squeezed vacuum at \(1064\,\mathrm{nm}\). Below-threshold operation yielded pure squeezing up to \(\sim15\) dB. State tomography employed balanced homodyne detection at a \(2.5\,\mathrm{MHz}\) sideband with common-mode rejection \(>80\) dB, and a machine-learning-enhanced quantum-state-tomography algorithm reconstructed \(\rho\) for 20 discrete pump powers. Because the squeezing operator implemented by the OPO is formally equivalent to time evolution under \(\hat H\), pump power could be mapped to an effective time \(\tau_{\rm eff}\propto|\alpha|\) [2111.08285].

The current itself was reconstructed numerically from a sequence of Wigner functions \(W_j(x,p)\) at effective times \(\tau_j\). An initial guess \(\mathbf J^{\rm init}\) was formed from \(\mathbf J_{\rm sys}+\mathbf J_{\rm env}\) with fitted \(\gamma\) and \(\bar n\), and the discrete continuity equation was then solved by minimizing \(\lVert(\delta J_x,\delta J_p)\rVert_1\) subject to the continuity constraint, cast as a linear-programming problem. The reconstructed current \(\mathbf J_{\rm exp}=\mathbf J^{\rm init}+\delta\mathbf J\) showed the expected separation into system, damping, and diffusion components, and confirmed that the central stagnation point retained \(Q=-1\) under both weak and strong decoherence [2111.08285].

These developments clarify the interpretive role of Wigner current. It is not merely an auxiliary rewriting of the Moyal equation. In closed systems it resolves flow reversal, interference, and topological defects more finely than a trajectory-based description; in open systems it separates Hamiltonian, damping, and diffusion mechanisms; and in generalized settings such as beam splitters, spin phase spaces, and billiards it provides a common geometric language for nonclassical transport, decoherence, and boundary effects [1411.3511, 2308.06706].

Source: https://www.emergentmind.com/topics/wigner-current