Skip to content

Phase-Space Conventions

This page fixes the phase-space conventions used in the Dynamics volume. It is a reference page for signs, factors of 2πℏ2\pi\hbar, and product rules. The full explanations live in Weyl Transform, Wigner Function, Star Product, and Moyal Bracket.

The conventions are chosen so that

⟨x∣p⟩=12πℏeipx/ℏ,\langle x|p\rangle = \frac{1}{\sqrt{2\pi\hbar}}e^{ipx/\hbar},

and

x⋆p−p⋆x=iℏ.x\star p-p\star x=i\hbar.

Different books may place signs or factors in different places. Translate definitions before comparing Wigner functions, Weyl symbols, or star products.

The position-momentum convention is

ψ(x)=12πℏ∫−∞∞dp eipx/ℏϕ(p),\psi(x) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty}dp\, e^{ipx/\hbar}\phi(p),

and

ϕ(p)=12πℏ∫−∞∞dx e−ipx/ℏψ(x).\phi(p) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty}dx\, e^{-ipx/\hbar}\psi(x).

Consequently,

⟨x∣p⟩=12πℏeipx/ℏ.\langle x|p\rangle = \frac{1}{\sqrt{2\pi\hbar}}e^{ipx/\hbar}.

This sign choice fixes the sign in the Wigner and Weyl transforms below. The site-wide convention page is Fourier Transform Conventions.

For an operator AA, its Weyl symbol is

AW(x,p)=∫−∞∞dy e−ipy/ℏ⟨x+y2|A|x−y2⟩.A_W(x,p) = \int_{-\infty}^{\infty}dy\, e^{-ipy/\hbar} \left\langle x+\frac{y}{2} \middle| A \middle| x-\frac{y}{2} \right\rangle.

Equivalently, if

A(x1,x2)=⟨x1∣A∣x2⟩,A(x_1,x_2) = \langle x_1|A|x_2\rangle,

then

AW(x,p)=∫dy e−ipy/ℏA(x+y2,x−y2).A_W(x,p) = \int dy\, e^{-ipy/\hbar} A\left( x+\frac{y}{2}, x-\frac{y}{2} \right).

The inverse transform is

⟨x1|A|x2⟩=∫dp2πℏ eip(x1−x2)/ℏAW(x1+x22,p).\left\langle x_1\middle|A\middle|x_2\right\rangle = \int\frac{dp}{2\pi\hbar}\, e^{ip(x_1-x_2)/\hbar} A_W\left( \frac{x_1+x_2}{2}, p \right).

The Weyl symbol is tied to symmetric ordering. For example,

(12(x^p^+p^x^))W=xp.\left( \frac12(\hat x\hat p+\hat p\hat x) \right)_W = xp.

The Wigner function of a density operator ρ\rho is the normalized Weyl symbol

Wρ(x,p)=12πℏρW(x,p).W_\rho(x,p) = \frac{1}{2\pi\hbar}\rho_W(x,p).

Explicitly,

Wρ(x,p)=12πℏ∫−∞∞dy e−ipy/ℏ⟨x+y2|ρ|x−y2⟩.W_\rho(x,p) = \frac{1}{2\pi\hbar} \int_{-\infty}^{\infty}dy\, e^{-ipy/\hbar} \left\langle x+\frac{y}{2} \middle| \rho \middle| x-\frac{y}{2} \right\rangle.

For a pure state,

Wψ(x,p)=12πℏ∫dy e−ipy/ℏψ(x+y2)ψ∗(x−y2).W_\psi(x,p) = \frac{1}{2\pi\hbar} \int dy\, e^{-ipy/\hbar} \psi\left(x+\frac{y}{2}\right) \psi^*\left(x-\frac{y}{2}\right).

With this normalization,

∫dx dp Wρ(x,p)=Tr⁡ρ.\int dx\,dp\,W_\rho(x,p) = \operatorname{Tr}\rho.

For a normalized state, the integral is 11.

The position marginal is

∫dp Wρ(x,p)=⟨x∣ρ∣x⟩.\int dp\,W_\rho(x,p) = \langle x|\rho|x\rangle.

The momentum marginal is

∫dx Wρ(x,p)=⟨p∣ρ∣p⟩.\int dx\,W_\rho(x,p) = \langle p|\rho|p\rangle.

These are ordinary probability densities. The full Wigner function is not an ordinary joint probability density, because it can be negative.

The trace measure for one degree of freedom is

dΓ=dx dp2πℏ.d\Gamma = \frac{dx\,dp}{2\pi\hbar}.

For suitable trace-class operators,

Tr⁡A=∫dΓ AW(x,p).\operatorname{Tr}A = \int d\Gamma\,A_W(x,p).

For products of operators,

Tr⁡(AB)=∫dΓ AW(x,p)BW(x,p).\operatorname{Tr}(AB) = \int d\Gamma\,A_W(x,p)B_W(x,p).

Since Wρ=ρW/(2πℏ)W_\rho=\rho_W/(2\pi\hbar), expectation values are written without an extra 1/(2πℏ)1/(2\pi\hbar):

Tr⁡(ρA)=∫dx dp Wρ(x,p)AW(x,p).\operatorname{Tr}(\rho A) = \int dx\,dp\, W_\rho(x,p)A_W(x,p).

This formula resembles a classical phase-space average, but the Wigner function may be negative and operator products use the star product.

For nn Cartesian degrees of freedom, write

q=(q1,…,qn),p=(p1,…,pn).\mathbf q=(q^1,\ldots,q^n), \qquad \mathbf p=(p_1,\ldots,p_n).

The Wigner function is

Wρ(q,p)=1(2πℏ)n∫dny e−ip⋅y/ℏ×⟨q+y2|ρ|q−y2⟩.\begin{aligned} W_\rho(\mathbf q,\mathbf p) &= \frac{1}{(2\pi\hbar)^n} \int d^n\mathbf y\, e^{-i\mathbf p\cdot\mathbf y/\hbar} \\ &\quad\times \left\langle \mathbf q+\frac{\mathbf y}{2} \middle| \rho \middle| \mathbf q-\frac{\mathbf y}{2} \right\rangle. \end{aligned}

The trace measure is

dΓn=dnq dnp(2πℏ)n.d\Gamma_n = \frac{d^n\mathbf q\,d^n\mathbf p} {(2\pi\hbar)^n}.

Again,

∫dnq dnp Wρ(q,p)=Tr⁡ρ.\int d^n\mathbf q\,d^n\mathbf p\, W_\rho(\mathbf q,\mathbf p) = \operatorname{Tr}\rho.

For non-Cartesian coordinates or constrained phase spaces, this flat measure is not automatically correct.

The Weyl symbol of an operator product is the star product of Weyl symbols:

(AB)W=AW⋆BW.(AB)_W = A_W\star B_W.

For one degree of freedom,

(A⋆B)(x,p)=A(x,p)exp⁡[iℏ2(∂x←∂p→−∂p←∂x→)]B(x,p).(A\star B)(x,p) = A(x,p) \exp\left[ \frac{i\hbar}{2} \left( \overleftarrow{\partial_x}\overrightarrow{\partial_p} - \overleftarrow{\partial_p}\overrightarrow{\partial_x} \right) \right] B(x,p).

The first terms are

A⋆B=AB+iℏ2{A,B}PB+O(ℏ2).A\star B = AB + \frac{i\hbar}{2}\{A,B\}_{\rm PB} + O(\hbar^2).

The basic check is

x⋆p=xp+iℏ2,p⋆x=xp−iℏ2.x\star p = xp+\frac{i\hbar}{2}, \qquad p\star x = xp-\frac{i\hbar}{2}.

Therefore,

x⋆p−p⋆x=iℏ.x\star p-p\star x=i\hbar.

For nn degrees of freedom, set

z=(q1,…,qn,p1,…,pn),\mathbf z=(q^1,\ldots,q^n,p_1,\ldots,p_n),

and

Ω=(0I−I0).\Omega = \begin{pmatrix} 0&I\\ -I&0 \end{pmatrix}.

Then the compact convention is

A⋆B=Aexp⁡[iℏ2∂z← TΩ∂z→]B.A\star B = A \exp\left[ \frac{i\hbar}{2} \overleftarrow{\partial_{\mathbf z}}^{\,T} \Omega \overrightarrow{\partial_{\mathbf z}} \right] B.

The Poisson bracket convention is

{A,B}PB=∂A∂x∂B∂p−∂A∂p∂B∂x.\{A,B\}_{\rm PB} = \frac{\partial A}{\partial x} \frac{\partial B}{\partial p} - \frac{\partial A}{\partial p} \frac{\partial B}{\partial x}.

For many Cartesian degrees of freedom,

{A,B}PB=∑a=1n(∂A∂qa∂B∂pa−∂A∂pa∂B∂qa).\{A,B\}_{\rm PB} = \sum_{a=1}^{n} \left( \frac{\partial A}{\partial q^a} \frac{\partial B}{\partial p_a} - \frac{\partial A}{\partial p_a} \frac{\partial B}{\partial q^a} \right).

This sign convention gives Hamilton’s equations in the form

q˙a={qa,H}PB=∂H∂pa,\dot q^a=\{q^a,H\}_{\rm PB} = \frac{\partial H}{\partial p_a},

and

p˙a={pa,H}PB=−∂H∂qa.\dot p_a=\{p_a,H\}_{\rm PB} = -\frac{\partial H}{\partial q^a}.

The Moyal bracket is

{A,B}M=1iℏ(A⋆B−B⋆A).\{A,B\}_M = \frac{1}{i\hbar} \left( A\star B-B\star A \right).

It obeys

{A,B}M={A,B}PB+O(ℏ2)\{A,B\}_M = \{A,B\}_{\rm PB} +O(\hbar^2)

for smooth functions in the formal small-ℏ\hbar expansion.

For a closed system with density operator ρ\rho,

∂ρ∂t=1iℏ[H,ρ].\frac{\partial\rho}{\partial t} = \frac{1}{i\hbar}[H,\rho].

With

Wρ=12πℏρW,W_\rho = \frac{1}{2\pi\hbar}\rho_W,

the phase-space evolution equation is

∂W∂t={HW,W}M.\frac{\partial W}{\partial t} = \{H_W,W\}_M.

For

HW(x,p)=p22m+V(x),H_W(x,p)=\frac{p^2}{2m}+V(x),

the leading classical part is

∂W∂t=−pm∂W∂x+V′(x)∂W∂p+O(ℏ2).\frac{\partial W}{\partial t} = - \frac{p}{m}\frac{\partial W}{\partial x} + V'(x)\frac{\partial W}{\partial p} + O(\hbar^2).

This equals the classical Liouville equation with the Poisson-bracket sign convention above:

∂W∂t={HW,W}PB+O(ℏ2).\frac{\partial W}{\partial t} = \{H_W,W\}_{\rm PB} + O(\hbar^2).

For quadratic Hamiltonians, the Moyal bracket equals the Poisson bracket exactly.

ItemConvention hereCommon alternative
Wigner prefactor1/(2πℏ)n1/(2\pi\hbar)^n1/hn1/h^n with h=2πℏh=2\pi\hbar
Separation phasee−ip⋅y/ℏe^{-i\mathbf p\cdot\mathbf y/\hbar}opposite sign if $\langle x
Trace measurednq dnp/(2πℏ)nd^nq\,d^np/(2\pi\hbar)^nabsorb measure into symbol definitions
Star productexp⁡[(iℏ/2)Λ]\exp[(i\hbar/2)\Lambda]opposite sign with opposite Fourier convention
Moyal bracket(A⋆B−B⋆A)/(iℏ)(A\star B-B\star A)/(i\hbar)definitions that absorb iℏi\hbar
Wigner functionWρ=ρW/(2πℏ)nW_\rho=\rho_W/(2\pi\hbar)^ncalling ρW\rho_W itself the Wigner function
  • Forgetting the prefactor in the Wigner transform.
  • Mixing a Weyl symbol from one convention with a Wigner function from another.
  • Treating W(x,p)W(x,p) as a positive joint probability distribution.
  • Using ordinary multiplication for symbols of noncommuting operator products.
  • Reversing the Poisson-bracket sign when translating Hamilton’s equations.
  • Dropping the phase-space trace measure 1/(2πℏ)n1/(2\pi\hbar)^n.
  • Assuming flat Cartesian formulas apply unchanged to spin, angle variables, curved spaces, constrained systems, or gauge systems.
  • H. Weyl, The Theory of Groups and Quantum Mechanics, Dover, 1931.
  • E. Wigner, “On the quantum correction for thermodynamic equilibrium,” Physical Review 40, 749-759, 1932.
  • J. E. Moyal, “Quantum mechanics as a statistical theory,” Proceedings of the Cambridge Philosophical Society 45, 99-124, 1949.
  • M. Hillery, R. F. O’Connell, M. O. Scully, and E. P. Wigner, “Distribution functions in physics: Fundamentals,” Physics Reports 106, 121-167, 1984.
  • C. K. Zachos, D. B. Fairlie, and T. L. Curtright, eds., Quantum Mechanics in Phase Space, World Scientific, 2005.
  • W. P. Schleich, Quantum Optics in Phase Space, Wiley-VCH, 2001.
  1. Starting from the Wigner definition, show that ∫dx dp Wρ(x,p)=Tr⁡ρ\int dx\,dp\,W_\rho(x,p)=\operatorname{Tr}\rho.
Solution

Use

Wρ(x,p)=12πℏ∫dy e−ipy/ℏ⟨x+y2|ρ|x−y2⟩.W_\rho(x,p) = \frac{1}{2\pi\hbar} \int dy\, e^{-ipy/\hbar} \left\langle x+\frac{y}{2} \middle| \rho \middle| x-\frac{y}{2} \right\rangle.

Integrating over pp gives

∫dp2πℏe−ipy/ℏ=δ(y).\int\frac{dp}{2\pi\hbar}e^{-ipy/\hbar} = \delta(y).

Therefore

∫dp Wρ(x,p)=⟨x∣ρ∣x⟩.\int dp\,W_\rho(x,p) = \langle x|\rho|x\rangle.

Integrating over xx gives

∫dx ⟨x∣ρ∣x⟩=Tr⁡ρ.\int dx\,\langle x|\rho|x\rangle = \operatorname{Tr}\rho.
  1. Use the star product to compute x⋆p−p⋆xx\star p-p\star x.
Solution

For x⋆px\star p, only the first derivative term contributes:

x⋆p=xp+iℏ2(∂x∂x∂p∂p)=xp+iℏ2.x\star p = xp + \frac{i\hbar}{2} \left( \frac{\partial x}{\partial x} \frac{\partial p}{\partial p} \right) = xp+\frac{i\hbar}{2}.

Similarly,

p⋆x=px−iℏ2=xp−iℏ2.p\star x = px - \frac{i\hbar}{2} = xp-\frac{i\hbar}{2}.

Subtracting gives

x⋆p−p⋆x=iℏ.x\star p-p\star x=i\hbar.
  1. Check the Poisson-bracket sign for H=p2/(2m)+V(x)H=p^2/(2m)+V(x).
Solution

With

{A,B}PB=AxBp−ApBx,\{A,B\}_{\rm PB} = A_xB_p-A_pB_x,

one finds

{x,H}PB=1⋅pm−0=pm,\{x,H\}_{\rm PB} = 1\cdot\frac{p}{m}-0 = \frac{p}{m},

and

{p,H}PB=0−1⋅V′(x)=−V′(x).\{p,H\}_{\rm PB} = 0- 1\cdot V'(x) = -V'(x).

These are Hamilton’s equations for the particle.

  1. A text defines the Wigner function with e+ipy/ℏe^{+ipy/\hbar} while keeping the same 1/(2πℏ)1/(2\pi\hbar) prefactor. What changed?
Solution

The separation Fourier sign changed. This is usually tied to using the conjugate position-momentum convention, or equivalently to replacing pp by −p-p relative to the convention used here. Star-product signs and formulas involving momentum arguments must be translated consistently; one cannot mix the two conventions term by term.