Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Gaussian pure states have nonnegative Wigner function

Proved
Wigner.gaussian_wigner_nonneg

by Lucas · Oct 18, 2026 · Mathlib 0df444a (Lean v4.33.1)

analysisquantum-mechanics

Gaussian pure states have nonnegative Wigner functions: the converse of Hudson's theorem. Let a,b,c∈Ca, b, c \in \mathbb{C}a,b,c∈C with Re(a)>0\mathrm{Re}(a) > 0Re(a)>0, so that

ψ(x)=e−ax2+bx+c\psi(x) = e^{-ax^2 + bx + c}ψ(x)=e−ax2+bx+c

is a square-integrable Gaussian wave packet, and let ℏ>0\hbar > 0ℏ>0. Then the Wigner function of ψ\psiψ is nonnegative at every phase-space point: W(x,p)≥0W(x,p) \ge 0W(x,p)≥0 for all x,px, px,p. The source text notes that "the Wigner function of any such state is a Gaussian distribution on phase space" — coherent and squeezed states are exactly the pure states with classical phase-space pictures.

Formalization Note No normalization of ψ\psiψ is assumed: positivity is invariant under the scalar factor ece^cec. The hypothesis Re(a)>0\mathrm{Re}(a) > 0Re(a)>0 is the source's assumption and is what makes ψ\psiψ a genuine wave packet.

Preamble
import Mathlib
import Definitions.Def_wigner

open MeasureTheory
Formal statement
namespace Wigner

theorem gaussian_wigner_nonneg {a b c : ℂ} {ℏ x p : ℝ} (ha : 0 < Complex.re a)
    (hℏ : 0 < ℏ) :
    0 ≤ wigner (fun y => Complex.exp (-(a * (y : ℂ) ^ 2) + b * (y : ℂ) + c)) ℏ x p := by sorry

end Wigner
Source
Wikipedia, "Wigner quasiprobability distribution" (oldid 1376057304), section "Positivity of the function". https://en.wikipedia.org/wiki/Wigner_quasiprobability_distribution
Read-back

What the Lean code literally says, in plain math · glm-5-latest-short

Let a,b,ca, b, ca,b,c be arbitrary complex numbers and ℏ,x,p\hbar, x, pℏ,x,p arbitrary real numbers, under two hypotheses: (a) Re⁡(a)>0\operatorname{Re}(a) > 0Re(a)>0; (b) 0<ℏ0 < \hbar0<ℏ. The conclusion: 0≤wigner(ψ,ℏ,x,p)0 \le \mathrm{wigner}(\psi, \hbar, x, p)0≤wigner(ψ,ℏ,x,p), where ψ\psiψ is the specific function y↦e−ay2+by+cy \mapsto e^{-a y^2 + b y + c}y↦e−ay2+by+c (with yyy coerced to a complex number in the exponent). Hypothesis (a) makes ψ\psiψ square-integrable, though the statement does not assume this separately; no normalization of ψ\psiψ is assumed.

Human review
  • Endorsed by Shuze Chen · Oct 18, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Lucas · Oct 18, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me