Skip to content

Moyal Bracket

The Moyal bracket is the phase-space version of the quantum commutator. In Hilbert-space language, dynamics is generated by

1iℏ[H,ρ].\frac{1}{i\hbar}[H,\rho].

In Wigner phase space, the same closed-system dynamics is generated by

{HW,Wρ}M,\{H_W,W_\rho\}_M,

where HWH_W is the Weyl symbol of the Hamiltonian, WρW_\rho is the Wigner function, and {⋅,⋅}M\{\cdot,\cdot\}_M is the Moyal bracket. It reduces to the classical Poisson bracket when the relevant ℏ\hbar corrections are negligible.

The Weyl transform maps operators to phase-space functions. Operator multiplication does not become ordinary multiplication. Instead, it becomes the star product. For one degree of freedom, a common convention is

(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 arrows specify which factor the derivative acts on. Equivalently, define the bidifferential operator Λ\Lambda by

AΛB=∂A∂x∂B∂p−∂A∂p∂B∂x.A\Lambda B = \frac{\partial A}{\partial x} \frac{\partial B}{\partial p} - \frac{\partial A}{\partial p} \frac{\partial B}{\partial x}.

Then the star product can be written schematically as

A⋆B=Aexp⁡(iℏ2Λ)B.A\star B = A\exp\left( \frac{i\hbar}{2}\Lambda \right)B.

The lowest 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),

where

{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}

is the classical Poisson bracket.

The Moyal bracket is the antisymmetric part of the star product:

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

Using Λ\Lambda, this may be written as

{A,B}M=2ℏAsin⁡(ℏ2Λ)B.\{A,B\}_M = \frac{2}{\hbar} A\sin\left( \frac{\hbar}{2}\Lambda \right)B.

Expanding the sine gives

{A,B}M=AΛB−ℏ224AΛ3B+O(ℏ4).\{A,B\}_M = A\Lambda B - \frac{\hbar^2}{24} A\Lambda^3B + O(\hbar^4).

The first term is the Poisson bracket. The remaining terms encode quantum corrections due to noncommutativity.

Let AWA_W and BWB_W be Weyl symbols of operators AA and BB. Then

(1iℏ[A,B])W={AW,BW}M.\left( \frac{1}{i\hbar}[A,B] \right)_W = \{A_W,B_W\}_M.

This is the conceptual core of the Moyal bracket. It is not a new law of motion; it is the commutator written in phase-space language.

For canonical variables, the Moyal bracket reproduces the basic commutator:

{x,p}M=1,\{x,p\}_M=1,

matching

1iℏ[x^,p^]=1.\frac{1}{i\hbar}[\hat x,\hat p]=1.

Because xx and pp are linear functions, there are no higher-derivative corrections in this example.

The closed-system density operator obeys the Liouville–von Neumann equation

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

Taking the Wigner transform gives

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

For

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

with the standard Weyl symbol

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

the equation becomes

∂W∂t=−pm∂W∂x+∑ℓ=0∞(−1)ℓ(2ℓ+1)!(ℏ2)2ℓ×d2ℓ+1Vdx2ℓ+1∂2ℓ+1W∂p2ℓ+1.\begin{aligned} \frac{\partial W}{\partial t} &= - \frac{p}{m}\frac{\partial W}{\partial x} + \sum_{\ell=0}^{\infty} \frac{(-1)^\ell}{(2\ell+1)!} \left( \frac{\hbar}{2} \right)^{2\ell} \\ &\quad \times \frac{d^{2\ell+1}V}{dx^{2\ell+1}} \frac{\partial^{2\ell+1}W}{\partial p^{2\ell+1}}. \end{aligned}

The ℓ=0\ell=0 term is

dVdx∂W∂p.\frac{dV}{dx}\frac{\partial W}{\partial p}.

Together with the kinetic term, it gives the classical Liouville flow:

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

The Moyal bracket has the formal classical limit

lim⁡ℏ→0{A,B}M={A,B}PB,\lim_{\hbar\to0}\{A,B\}_M = \{A,B\}_{\rm PB},

provided the functions are smooth enough on the relevant phase-space scale. More physically, the Poisson bracket dominates when the Wigner function and Hamiltonian vary slowly compared with the ℏ\hbar-scale oscillations where quantum interference is important.

For quadratic Hamiltonians, all third and higher derivatives vanish, so the Moyal bracket equals the Poisson bracket exactly. This is why free particles, harmonic oscillators, and general quadratic systems move Wigner functions by classical phase-space flows even though the states remain quantum.

For anharmonic systems, higher derivatives of V(x)V(x) generate quantum corrections. These corrections are the phase-space version of noncommutative operator dynamics.

For

H=p22m,H=\frac{p^2}{2m},

the Moyal equation reduces to

∂W∂t=−pm∂W∂x.\frac{\partial W}{\partial t} = - \frac{p}{m}\frac{\partial W}{\partial x}.

Equivalently,

∂W∂t+pm∂W∂x=0.\frac{\partial W}{\partial t} + \frac{p}{m}\frac{\partial W}{\partial x} =0.

The solution is transport along classical straight-line trajectories:

W(x,p,t)=W(x−ptm,p,0).W(x,p,t) = W\left( x-\frac{pt}{m}, p, 0 \right).

The Wigner function may still have quantum features, such as negative regions or interference fringes, but the free-particle flow translates them classically in phase space.

For

H=p22m+12mω2x2,H=\frac{p^2}{2m}+\frac12m\omega^2x^2,

the Moyal bracket again equals the Poisson bracket exactly:

∂W∂t=−pm∂W∂x+mω2x∂W∂p.\frac{\partial W}{\partial t} = - \frac{p}{m}\frac{\partial W}{\partial x} + m\omega^2x\frac{\partial W}{\partial p}.

This is rigid rotation in oscillator phase space. Coherent-state Wigner functions therefore rotate like classical oscillator ellipses without changing shape.

Energy eigenstate Wigner functions are stationary because their density operators commute with HH. In phase-space language,

{H,Wn}M=0.\{H,W_n\}_M=0.

For an anharmonic potential

V(x)=λ4x4,V(x)=\frac{\lambda}{4}x^4,

the first quantum correction comes from V′′′(x)=6λxV'''(x)=6\lambda x. The Wigner evolution contains

∂W∂t=−pm∂W∂x+λx3∂W∂p−ℏ2λx4∂3W∂p3+⋯ .\frac{\partial W}{\partial t} = - \frac{p}{m}\frac{\partial W}{\partial x} + \lambda x^3\frac{\partial W}{\partial p} - \frac{\hbar^2\lambda x}{4} \frac{\partial^3 W}{\partial p^3} +\cdots.

The first two terms are the classical Liouville flow. The third term is the leading Moyal correction.

  • Replacing the star product by ordinary multiplication in commutators.
  • Assuming the Moyal bracket is always just the Poisson bracket.
  • Forgetting that the Hamiltonian must be represented by the correct Weyl symbol.
  • Ignoring sign and normalization conventions in the Wigner transform.
  • Treating exact classical-looking flow for quadratic Hamiltonians as proof that the state is classical.
  • Applying the ℏ\hbar expansion where the Wigner function has sharp interference features on the ℏ\hbar scale.
  • H. J. Groenewold, “On the Principles of elementary quantum mechanics,” Physica 12, 405-460, 1946.
  • J. E. Moyal, “Quantum mechanics as a statistical theory,” Mathematical 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. B. Case, “Wigner functions and Weyl transforms for pedestrians,” American Journal of Physics 76, 937-946, 2008.
  1. Show that {x,p}M=1\{x,p\}_M=1.
Solution

Because xx and pp are linear functions, the star product terminates after the first derivative term:

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

and

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

Therefore

{x,p}M=1iℏ(x⋆p−p⋆x)=1.\{x,p\}_M = \frac{1}{i\hbar} \left( x\star p-p\star x \right) =1.
  1. Derive the free-particle Wigner equation from the Moyal bracket.
Solution

For

H=p22m,H=\frac{p^2}{2m},

the Hamiltonian depends only on pp. Since it is quadratic, the Moyal bracket equals the Poisson bracket:

{H,W}M={H,W}PB.\{H,W\}_M = \{H,W\}_{\rm PB}.

Compute

{H,W}PB=∂H∂x∂W∂p−∂H∂p∂W∂x=−pm∂W∂x.\{H,W\}_{\rm PB} = \frac{\partial H}{\partial x} \frac{\partial W}{\partial p} - \frac{\partial H}{\partial p} \frac{\partial W}{\partial x} = - \frac{p}{m}\frac{\partial W}{\partial x}.

Thus

∂W∂t=−pm∂W∂x.\frac{\partial W}{\partial t} = - \frac{p}{m}\frac{\partial W}{\partial x}.
  1. Why do quadratic Hamiltonians have exact classical phase-space flow in the Wigner representation?
Solution

The Moyal bracket expansion is

{A,B}M={A,B}PB−ℏ224AΛ3B+O(ℏ4).\{A,B\}_M = \{A,B\}_{\rm PB} - \frac{\hbar^2}{24}A\Lambda^3B + O(\hbar^4).

If the Hamiltonian is at most quadratic in xx and pp, all third and higher derivatives of the Hamiltonian vanish. Therefore the higher Moyal terms vanish, and

{H,W}M={H,W}PB.\{H,W\}_M=\{H,W\}_{\rm PB}.

The state may still be nonclassical; the statement is about the form of the phase-space flow.

  1. Find the leading Moyal correction for V(x)=λx4/4V(x)=\lambda x^4/4.
Solution

For

V(x)=λ4x4,V(x)=\frac{\lambda}{4}x^4,

one has

V′(x)=λx3,V′′′(x)=6λx.V'(x)=\lambda x^3, \qquad V'''(x)=6\lambda x.

The potential part of the Wigner equation begins

V′(x)∂W∂p−ℏ224V′′′(x)∂3W∂p3+⋯ .V'(x)\frac{\partial W}{\partial p} - \frac{\hbar^2}{24} V'''(x) \frac{\partial^3W}{\partial p^3} +\cdots.

Substituting V′′′(x)=6λxV'''(x)=6\lambda x gives

λx3∂W∂p−ℏ2λx4∂3W∂p3+⋯ .\lambda x^3\frac{\partial W}{\partial p} - \frac{\hbar^2\lambda x}{4} \frac{\partial^3W}{\partial p^3} +\cdots.