Canonical commutation relation

In quantum mechanics, the canonical commutation relation is the fundamental relation between canonical conjugate quantities (quantities which are related by definition such that one is the Fourier transform of another). For example, $${\displaystyle [{\hat {x}},{\hat {p}}_{x}]=i\hbar \mathbb {I} }$$

between the position operator x and momentum operator p_x in the x direction of a point particle in one dimension, where [x , p_x] = x p_x − p_x x is the commutator of x and p_x , i is the imaginary unit, and ℏ is the reduced Planck constant h/2π, and ${\displaystyle \mathbb {I} }$ is the unit operator. In general, position and momentum are vectors of operators and their commutation relation between different components of position and momentum can be expressed as $${\displaystyle [{\hat {x}}_{i},{\hat {p}}_{j}]=i\hbar \delta _{ij},}$$ where ${\displaystyle \delta _{ij}}$ is the Kronecker delta.

This relation is attributed to Werner Heisenberg, Max Born and Pascual Jordan (1925),^[1][2] who called it a "quantum condition" serving as a postulate of the theory; it was noted by E. Kennard (1927)^[3] to imply the Heisenberg uncertainty principle. The Stone–von Neumann theorem gives a uniqueness result for operators satisfying (an exponentiated form of) the canonical commutation relation.

By contrast, in classical physics, all observables commute and the commutator would be zero. However, an analogous relation exists, which is obtained by replacing the commutator with the Poisson bracket multiplied by iℏ, $${\displaystyle \{x,p\}=1\,.}$$

This observation led Dirac to propose that the quantum counterparts ${\displaystyle {\hat {f}}}$, ĝ of classical observables f, g satisfy $${\displaystyle [{\hat {f}},{\hat {g}}]=i\hbar {\widehat {\{f,g\}}}\,.}$$

In 1946, Hip Groenewold demonstrated that a general systematic correspondence between quantum commutators and Poisson brackets could not hold consistently.^[4][5]

However, he further appreciated that such a systematic correspondence does, in fact, exist between the quantum commutator and a deformation of the Poisson bracket, today called the Moyal bracket, and, in general, quantum operators and classical observables and distributions in phase space. He thus finally elucidated the consistent correspondence mechanism, the Wigner–Weyl transform, that underlies an alternate equivalent mathematical representation of quantum mechanics known as deformation quantization.^[4][6]

According to the correspondence principle, in certain limits the quantum equations of states must approach Hamilton's equations of motion. The latter state the following relation between the generalized coordinate q (e.g. position) and the generalized momentum p: $${\displaystyle {\begin{cases}{\dot {q}}={\frac {\partial H}{\partial p}}=\{q,H\};\\{\dot {p}}=-{\frac {\partial H}{\partial q}}=\{p,H\}.\end{cases}}}$$

In quantum mechanics the Hamiltonian ${\displaystyle {\hat {H}}}$, (generalized) coordinate ${\displaystyle {\hat {Q}}}$ and (generalized) momentum ${\displaystyle {\hat {P}}}$ are all linear operators.

The time derivative of a quantum state is - ${\displaystyle i{\hat {H}}/\hbar }$ (by Schrödinger equation). Equivalently, since in the Schrödinger picture the operators are not explicitly time-dependent, the operators can be seen to be evolving in time (for a contrary perspective where the operators are time dependent, see Heisenberg picture) according to their commutation relation with the Hamiltonian: $${\displaystyle {\frac {d{\hat {Q}}}{dt}}={\frac {i}{\hbar }}[{\hat {H}},{\hat {Q}}]}$$ $${\displaystyle {\frac {d{\hat {P}}}{dt}}={\frac {i}{\hbar }}[{\hat {H}},{\hat {P}}]\,\,.}$$

In order for that to reconcile in the classical limit with Hamilton's equations of motion, ${\displaystyle [{\hat {H}},{\hat {Q}}]}$ must depend entirely on the appearance of ${\displaystyle {\hat {P}}}$ in the Hamiltonian and ${\displaystyle [{\hat {H}},{\hat {P}}]}$ must depend entirely on the appearance of ${\displaystyle {\hat {Q}}}$ in the Hamiltonian. Further, since the Hamiltonian operator depends on the (generalized) coordinate and momentum operators, it can be viewed as a functional, and we may write (using functional derivatives): $${\displaystyle [{\hat {H}},{\hat {Q}}]={\frac {\delta {\hat {H}}}{\delta {\hat {P}}}}\cdot [{\hat {P}},{\hat {Q}}]}$$ $${\displaystyle [{\hat {H}},{\hat {P}}]={\frac {\delta {\hat {H}}}{\delta {\hat {Q}}}}\cdot [{\hat {Q}},{\hat {P}}]\,.}$$

In order to obtain the classical limit we must then have $${\displaystyle [{\hat {Q}},{\hat {P}}]=i\hbar ~\mathbb {I} .}$$

The group ${\displaystyle H_{3}(\mathbb {R} )}$ generated by exponentiation of the 3-dimensional Lie algebra determined by the commutation relation ${\displaystyle [{\hat {x}},{\hat {p}}]=i\hbar }$ is called the Heisenberg group. This group can be realized as the group of ${\displaystyle 3\times 3}$ upper triangular matrices with ones on the diagonal.^[7]

According to the standard mathematical formulation of quantum mechanics, quantum observables such as ${\displaystyle {\hat {x}}}$ and ${\displaystyle {\hat {p}}}$ should be represented as self-adjoint operators on some Hilbert space. It is relatively easy to see that two operators satisfying the above canonical commutation relations cannot both be bounded. Certainly, if ${\displaystyle {\hat {x}}}$ and ${\displaystyle {\hat {p}}}$ were trace class operators, the relation ${\displaystyle \operatorname {Tr} (AB)=\operatorname {Tr} (BA)}$ gives a nonzero number on the right and zero on the left.

Alternately, if ${\displaystyle {\hat {x}}}$ and ${\displaystyle {\hat {p}}}$ were bounded operators, note that ${\displaystyle [{\hat {x}}^{n},{\hat {p}}]=i\hbar n{\hat {x}}^{n-1}}$, hence the operator norms would satisfy $${\displaystyle 2\left\|{\hat {p}}\right\|\left\|{\hat {x}}^{n-1}\right\|\left\|{\hat {x}}\right\|\geq n\hbar \left\|{\hat {x}}^{n-1}\right\|,}$$ so that, for any n, $${\displaystyle 2\left\|{\hat {p}}\right\|\left\|{\hat {x}}\right\|\geq n\hbar }$$ However, n can be arbitrarily large, so at least one operator cannot be bounded, and the dimension of the underlying Hilbert space cannot be finite. If the operators satisfy the Weyl relations (an exponentiated version of the canonical commutation relations, described below) then as a consequence of the Stone–von Neumann theorem, both operators must be unbounded.

Still, these canonical commutation relations can be rendered somewhat "tamer" by writing them in terms of the (bounded) unitary operators ${\displaystyle \exp(it{\hat {x}})}$ and ${\displaystyle \exp(is{\hat {p}})}$. The resulting braiding relations for these operators are the so-called Weyl relations $${\displaystyle \exp(it{\hat {x}})\exp(is{\hat {p}})=\exp(-ist/\hbar )\exp(is{\hat {p}})\exp(it{\hat {x}}).}$$ These relations may be thought of as an exponentiated version of the canonical commutation relations; they reflect that translations in position and translations in momentum do not commute. One can easily reformulate the Weyl relations in terms of the representations of the Heisenberg group.

The uniqueness of the canonical commutation relations—in the form of the Weyl relations—is then guaranteed by the Stone–von Neumann theorem.

For technical reasons, the Weyl relations are not strictly equivalent to the canonical commutation relation ${\displaystyle [{\hat {x}},{\hat {p}}]=i\hbar }$. If ${\displaystyle {\hat {x}}}$ and ${\displaystyle {\hat {p}}}$ were bounded operators, then a special case of the Baker–Campbell–Hausdorff formula would allow one to "exponentiate" the canonical commutation relations to the Weyl relations.^[8] Since, as we have noted, any operators satisfying the canonical commutation relations must be unbounded, the Baker–Campbell–Hausdorff formula does not apply without additional domain assumptions. Indeed, counterexamples exist satisfying the canonical commutation relations but not the Weyl relations.^[9] (These same operators give a counterexample to the naive form of the uncertainty principle.) These technical issues are the reason that the Stone–von Neumann theorem is formulated in terms of the Weyl relations.

A discrete version of the Weyl relations, in which the parameters s and t range over ${\displaystyle \mathbb {Z} /n}$, can be realized on a finite-dimensional Hilbert space by means of the clock and shift matrices.

The simple formula $${\displaystyle [x,p]=i\hbar \,\mathbb {I} ~,}$$ valid for the quantization of the simplest classical system, can be generalized to the case of an arbitrary Lagrangian ${\displaystyle {\mathcal {L}}}$.^[10] We identify canonical coordinates (such as x in the example above, or a field Φ(x) in the case of quantum field theory) and canonical momenta π_x (in the example above it is p, or more generally, some functions involving the derivatives of the canonical coordinates with respect to time): $${\displaystyle \pi _{i}\ {\stackrel {\mathrm {def} }{=}}\ {\frac {\partial {\mathcal {L}}}{\partial (\partial x_{i}/\partial t)}}.}$$

This definition of the canonical momentum ensures that one of the Euler–Lagrange equations has the form $${\displaystyle {\frac {\partial }{\partial t}}\pi _{i}={\frac {\partial {\mathcal {L}}}{\partial x_{i}}}.}$$

The canonical commutation relations then amount to $${\displaystyle [x_{i},\pi _{j}]=i\hbar \delta _{ij}\,}$$ where δ_ij is the Kronecker delta.

Further, it can be shown that $${\displaystyle [F({\vec {x}}),p_{i}]=i\hbar {\frac {\partial F({\vec {x}})}{\partial x_{i}}};\qquad [x_{i},F({\vec {p}})]=i\hbar {\frac {\partial F({\vec {p}})}{\partial p_{i}}}.}$$

Using ${\displaystyle C_{n+1}^{k}=C_{n}^{k}+C_{n}^{k-1}}$, it can be shown that by mathematical induction $${\displaystyle \left[{\hat {x}}^{n},{\hat {p}}^{m}\right]=\sum _{k=1}^{\min \left(m,n\right)}{{\frac {-\left(-i\hbar \right)^{k}n!m!}{k!\left(n-k\right)!\left(m-k\right)!}}{\hat {x}}^{n-k}{\hat {p}}^{m-k}}=\sum _{k=1}^{\min \left(m,n\right)}{{\frac {\left(i\hbar \right)^{k}n!m!}{k!\left(n-k\right)!\left(m-k\right)!}}{\hat {p}}^{m-k}{\hat {x}}^{n-k}},}$$ generally known as McCoy's formula.^[11]

Canonical quantization is applied, by definition, on canonical coordinates. However, in the presence of an electromagnetic field, the canonical momentum p is not gauge invariant. The correct gauge-invariant momentum (or "kinetic momentum") is

${\displaystyle p_{\text{kin}}=p-qA\,\!}$ (SI units)      ${\displaystyle p_{\text{kin}}=p-{\frac {qA}{c}}\,\!}$ (cgs units),

where q is the particle's electric charge, A is the vector potential, and c is the speed of light. Although the quantity p_kin is the "physical momentum", in that it is the quantity to be identified with momentum in laboratory experiments, it does not satisfy the canonical commutation relations; only the canonical momentum does that. This can be seen as follows.

The non-relativistic Hamiltonian for a quantized charged particle of mass m in a classical electromagnetic field is (in cgs units) $${\displaystyle H={\frac {1}{2m}}\left(p-{\frac {qA}{c}}\right)^{2}+q\phi }$$ where A is the three-vector potential and φ is the scalar potential. This form of the Hamiltonian, as well as the Schrödinger equation Hψ = iħ∂ψ/∂t, the Maxwell equations and the Lorentz force law are invariant under the gauge transformation $${\displaystyle A\to A'=A+\nabla \Lambda }$$ $${\displaystyle \phi \to \phi '=\phi -{\frac {1}{c}}{\frac {\partial \Lambda }{\partial t}}}$$ $${\displaystyle \psi \to \psi '=U\psi }$$ $${\displaystyle H\to H'=UHU^{\dagger },}$$ where $${\displaystyle U=\exp \left({\frac {iq\Lambda }{\hbar c}}\right)}$$ and Λ = Λ(x,t) is the gauge function.

The angular momentum operator is $${\displaystyle L=r\times p\,\!}$$ and obeys the canonical quantization relations $${\displaystyle [L_{i},L_{j}]=i\hbar {\epsilon _{ijk}}L_{k}}$$ defining the Lie algebra for so(3), where ${\displaystyle \epsilon _{ijk}}$ is the Levi-Civita symbol. Under gauge transformations, the angular momentum transforms as $${\displaystyle \langle \psi \vert L\vert \psi \rangle \to \langle \psi ^{\prime }\vert L^{\prime }\vert \psi ^{\prime }\rangle =\langle \psi \vert L\vert \psi \rangle +{\frac {q}{\hbar c}}\langle \psi \vert r\times \nabla \Lambda \vert \psi \rangle \,.}$$

The gauge-invariant angular momentum (or "kinetic angular momentum") is given by $${\displaystyle K=r\times \left(p-{\frac {qA}{c}}\right),}$$ which has the commutation relations $${\displaystyle [K_{i},K_{j}]=i\hbar {\epsilon _{ij}}^{\,k}\left(K_{k}+{\frac {q\hbar }{c}}x_{k}\left(x\cdot B\right)\right)}$$ where $${\displaystyle B=\nabla \times A}$$ is the magnetic field. The inequivalence of these two formulations shows up in the Zeeman effect and the Aharonov–Bohm effect.

All such nontrivial commutation relations for pairs of operators lead to corresponding uncertainty relations,^[12] involving positive semi-definite expectation contributions by their respective commutators and anticommutators. In general, for two Hermitian operators A and B, consider expectation values in a system in the state ψ, the variances around the corresponding expectation values being (ΔA)² ≡ ⟨(A − ⟨A⟩)²⟩, etc.

Then $${\displaystyle \Delta A\,\Delta B\geq {\frac {1}{2}}{\sqrt {\left|\left\langle \left[{A},{B}\right]\right\rangle \right|^{2}+\left|\left\langle \left\{A-\langle A\rangle ,B-\langle B\rangle \right\}\right\rangle \right|^{2}}},}$$ where [A, B] ≡ A B − B A is the commutator of A and B, and {A, B} ≡ A B + B A is the anticommutator.

This follows through use of the Cauchy–Schwarz inequality, since |⟨A²⟩| |⟨B²⟩| ≥ |⟨A B⟩|², and A B = ([A, B] + {A, B})/2 ; and similarly for the shifted operators A − ⟨A⟩ and B − ⟨B⟩. (Cf. uncertainty principle derivations.)

Substituting for A and B (and taking care with the analysis) yield Heisenberg's familiar uncertainty relation for x and p, as usual.

For the angular momentum operators L_x = y p_z − z p_y, etc., one has that $${\displaystyle [{L_{x}},{L_{y}}]=i\hbar \epsilon _{xyz}{L_{z}},}$$ where ${\displaystyle \epsilon _{xyz}}$ is the Levi-Civita symbol and simply reverses the sign of the answer under pairwise interchange of the indices. An analogous relation holds for the spin operators.

Here, for L_x and L_y ,^[12] in angular momentum multiplets ψ = |ℓ,m⟩, one has, for the transverse components of the Casimir invariant L_x² + L_y²+ L_z², the z-symmetric relations

⟨L_x²⟩ = ⟨L_y²⟩ = (ℓ (ℓ + 1) − m²) ℏ²/2 ,

as well as ⟨L_x⟩ = ⟨L_y⟩ = 0 .

Consequently, the above inequality applied to this commutation relation specifies $${\displaystyle \Delta L_{x}\Delta L_{y}\geq {\frac {1}{2}}{\sqrt {\hbar ^{2}|\langle L_{z}\rangle |^{2}}}~,}$$ hence $${\displaystyle {\sqrt {|\langle L_{x}^{2}\rangle \langle L_{y}^{2}\rangle |}}\geq {\frac {\hbar ^{2}}{2}}\vert m\vert }$$ and therefore $${\displaystyle \ell (\ell +1)-m^{2}\geq |m|~,}$$ so, then, it yields useful constraints such as a lower bound on the Casimir invariant: ℓ (ℓ + 1) ≥ |m| (|m| + 1), and hence ℓ ≥ |m|, among others.

• Canonical quantization
• CCR and CAR algebras
• Conformastatic spacetimes
• Lie derivative
• Moyal bracket
• Stone–von Neumann theorem

• Hall, Brian C. (2013), Quantum Theory for Mathematicians, Graduate Texts in Mathematics, vol. 267, Springer.
• Hall, Brian C. (2015), Lie Groups, Lie Algebras and Representations, An Elementary Introduction, Graduate Texts in Mathematics, vol. 222 (2nd ed.), Springer.