Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Derivation of L, X and P Operators from Commutation Relations

This document provides a bottom-up derivation starting with the commutation relations of the operators for position XiX_i, momentum PiP_i, and angular momentum LiL_i. We compute their form, demonstrate their infinite-dimensional nature, and prove their uniqueness in the Hilbert space L2(R3)L^2(\mathbb{R}^3).

Commutation Relations

The four fundamental commutation relations are:

  1. [Li,Lj]=iℏϵijkLk[L_i, L_j] = i \hbar \epsilon_{ijk} L_k

  2. [Li,Xj]=iℏϵijkXk[L_i, X_j] = i \hbar \epsilon_{ijk} X_k

  3. [Li,Pj]=iℏϵijkPk[L_i, P_j] = i \hbar \epsilon_{ijk} P_k

  4. [Xi,Pj]=iℏδij[X_i, P_j] = i \hbar \delta_{ij}

All generators are Hermitian. The relation 1. describes the angular momentum Lie algebra, relations 2. and 3. indicate that XiX_i and PiP_i are vector operators with repsect to LiL_i. The relation 4. is the canonical commutation relation for position and momentum.

Here, ϵijk\epsilon_{ijk} is the Levi-Civita symbol, δij\delta_{ij} is the Kronecker delta, and ℏ\hbar is the reduced Planck constant.

Infinite-Dimensional Nature of XiX_i and PiP_i

The operators XiX_i and PiP_i must be infinite-dimensional due to the commutation relation [Xi,Pj]=iℏδij[X_i, P_j] = i \hbar \delta_{ij}. In a finite-dimensional Hilbert space of dimension nn, the trace of a commutator [A,B][A, B] is zero:

Tr([A,B])=Tr(AB−BA)=Tr(AB)−Tr(BA)=Tr(AB)−Tr(AB)=0 .\text{Tr}([A, B]) =\text{Tr}(AB-BA) =\text{Tr}(AB)-\text{Tr}(BA) =\text{Tr}(AB)-\text{Tr}(AB) =0\,.

Thus:

This contradiction implies that nn cannot be finite; hence, the Hilbert space must be infinite-dimensional, such as L2(R3)L^2(\mathbb{R}^3).

Both XiX_i and PiP_i must be infinite-dimensional, and consequently also LiL_i, otherwise the commutators like [Li,Xj][L_i, X_j] would not make sense. All the generators must be from the same space. When infinite-dimensional space is needed, the L2(R3)L^2(\mathbb{R}^3) is usually used, so we will use it as well.

Form of XiX_i and PiP_i

The Stone-von Neumann theorem says that XiX_i and PiP_i are uniquely determined as below, except a unitary transformation.

In the position representation on L2(R3)L^2(\mathbb{R}^3):

To verify:

[Xi,Pj]ψ(x)=XiPjψ(x)−PjXiψ(x)=xi(−iℏ∂ψ∂xj)−(−iℏ∂∂xj(xiψ))[X_i, P_j] \psi(x) = X_i P_j \psi(x) - P_j X_i \psi(x) = x_i \left(-i \hbar \frac{\partial \psi}{\partial x_j}\right) - \left(-i \hbar \frac{\partial}{\partial x_j} (x_i \psi)\right)

Using the product rule:

=−iℏxi∂ψ∂xj+iℏ(∂xi∂xjψ+xi∂ψ∂xj)=−iℏxi∂ψ∂xj+iℏδijψ+iℏxi∂ψ∂xj=iℏδijψ= -i \hbar x_i \frac{\partial \psi}{\partial x_j} + i \hbar \left( \frac{\partial x_i}{\partial x_j} \psi + x_i \frac{\partial \psi}{\partial x_j} \right) = -i \hbar x_i \frac{\partial \psi}{\partial x_j} + i \hbar \delta_{ij} \psi + i \hbar x_i \frac{\partial \psi}{\partial x_j} = i \hbar \delta_{ij} \psi

Thus, [Xi,Pj]=iℏδij[X_i, P_j] = i \hbar \delta_{ij}, confirming the forms.

Unitary Rotation Operator UU

In a previous chapter we derived for any vector operator (we will use the vector operator XiX_i):

UXiU−1=RjiXj .U X_i U^{-1} = R^j{}_i X_j\,.

Now we can write, using the known properties of XiX_i derived above:

Xi∣x⟩=xi∣x⟩ ,X_i |x\rangle = x_i |x\rangle\,,
UXi∣x⟩=Uxi∣x⟩ ,U X_i |x\rangle = U x_i |x\rangle\,,
UXi∣x⟩=xiU∣x⟩ ,U X_i |x\rangle = x_i U |x\rangle\,,
RjiXjU∣x⟩=xiU∣x⟩ ,R^j{}_i X_j U |x\rangle = x_i U |x\rangle\,,
XjU∣x⟩=(R−1)jixiU∣x⟩ .X_j U |x\rangle = (R^{-1})_j{}^i x_i U |x\rangle\,.

Comparing with Xi∣x′⟩=xi′∣x′⟩X_i |x'\rangle = x_i' |x'\rangle we get U∣x⟩=∣x′⟩U |x\rangle=|x'\rangle and xi′=(R−1)jixix_i' = (R^{-1})_j{}^i x_i, so:

U∣x⟩=∣x′⟩=∣R−1x⟩ .U |x\rangle=|x'\rangle = | R^{-1} x\rangle\,.

Then we can compute:

⟨x∣ψ′⟩=ψ′(x)=Uψ(x)=⟨x∣U∣ψ′⟩=⟨R−1x∣ψ′⟩=ψ(R−1x) .\langle x | \psi' \rangle =\psi'(x) = U\psi(x) = \langle x | U | \psi' \rangle = \langle R^{-1} x | \psi' \rangle = \psi(R^{-1} x)\,.

So we derived:

Uψ(x)=ψ(R−1x) .U\psi(x) = \psi(R^{-1} x)\,.

For a rotation by angle θ\theta around axis n\mathbf{n}, the unitary operator UU is:

U=exp⁡(−iℏθnkLk)U = \exp\left(-\frac{i}{\hbar} \theta n^k L_k\right)

Acting on a wave function:

Uψ(x)=ψ(R−1x)U \psi(x) = \psi(R^{-1} x)

where RR is the rotation matrix, and R−1R^{-1} is its inverse.

Derivation of LiL_i

Consider an infinitesimal rotation by δθ\delta \theta around n\mathbf{n}. The rotation matrix approximates as Rij≈δij+δθϵijknkR_{ij} \approx \delta_{ij} + \delta \theta \epsilon_{ijk} n_k, and:

R−1xi≈xi−δθ(n×x)i,(n×x)i=ϵijknjxkR^{-1} x_i \approx x_i - \delta \theta (\mathbf{n} \times \mathbf{x})_i, \quad (\mathbf{n} \times \mathbf{x})_i = \epsilon_{ijk} n_j x_k

Expand ψ(R−1x)\psi(R^{-1} x):

ψ(R−1x)≈ψ(x)−δθϵijknjxk∂ψ∂xi\psi(R^{-1} x) \approx \psi(x) - \delta \theta \epsilon_{ijk} n_j x_k \frac{\partial \psi}{\partial x_i}

For UU:

U≈1−iℏδθnmLm,Uψ(x)≈ψ(x)−iℏδθnmLmψ(x)U \approx 1 - \frac{i}{\hbar} \delta \theta n_m L_m, \quad U \psi(x) \approx \psi(x) - \frac{i}{\hbar} \delta \theta n_m L_m \psi(x)

Equate:

ψ(x)−iℏδθnmLmψ=ψ(x)−δθϵijknjxk∂iψ\psi(x) - \frac{i}{\hbar} \delta \theta n_m L_m \psi = \psi(x) - \delta \theta \epsilon_{ijk} n_j x_k \partial_i \psi

Subtract ψ(x)\psi(x) and simplify:

iℏnmLmψ=ϵijknjxk∂iψ\frac{i}{\hbar} n_m L_m \psi = \epsilon_{ijk} n_j x_k \partial_i \psi

Since this holds for all n\mathbf{n}, equate coefficients:

iℏLiψ=ϵijkxk∂jψ\frac{i}{\hbar} L_i \psi = \epsilon_{ijk} x_k \partial_j \psi

Substitute Pj=−iℏ∂jP_j = -i \hbar \partial_j, so ∂jψ=1−iℏPjψ=iℏPjψ\partial_j \psi = \frac{1}{-i \hbar} P_j \psi = \frac{i}{\hbar} P_j \psi:

iℏLiψ=ϵijkxk(iℏPjψ)=iℏϵijkXkPjψ\frac{i}{\hbar} L_i \psi = \epsilon_{ijk} x_k \left( \frac{i}{\hbar} P_j \psi \right) = \frac{i}{\hbar} \epsilon_{ijk} X_k P_j \psi

Multiply by ℏi\frac{\hbar}{i}:

Liψ=ϵijkXjPkψL_i \psi = \epsilon_{ijk} X_j P_k \psi

Thus, Li=ϵijkXjPkL_i = \epsilon_{ijk} X_j P_k.

Uniqueness of LiL_i

Suppose Li′L_i' also satisfies the commutation relations. Define Di=Li−Li′D_i = L_i - L_i':

DiD_i commutes with all XjX_j and PjP_j. In the irreducible representation on L2(R3)L^2(\mathbb{R}^3), Schur’s lemma implies Di=ciID_i = c_i I. Check the algebra:

[Li,Lj]=[Li′+Di,Lj′+Dj]=[Li′,Lj′]=iℏϵijkLk′,[Li,Lj]=iℏϵijk(Lk′+ckI)[L_i, L_j] = [L_i' + D_i, L_j' + D_j] = [L_i', L_j'] = i \hbar \epsilon_{ijk} L_k', \quad [L_i, L_j] = i \hbar \epsilon_{ijk} (L_k' + c_k I)

Thus, ck=0c_k = 0, so Di=0D_i = 0, and Li=Li′L_i = L_i'.

Summary

The four commutation relations uniquely determine the generators:

No other solution is possible. One can then extend this Hilbert space for example to also include spin, but in the smallest space L2(R3)L^2(\mathbb{R}^3) the above solution is unique.

Group

To determine the group elements, we can define:

T(a)=eia⋅P .T(a) = e^{i {\bf a}\cdot{\bf P}}\,.

That’s a translation operator: T(a)ψ(x)=ψ(x+a)T(a)\psi(x) = \psi(x+a). We derived that P{\bf P} is a translation generator.

Next we compute from BCH:

eia⋅XPje−ia⋅X=Pj−aje^{i {\bf a}\cdot{\bf X}} P_j e^{-i {\bf a}\cdot{\bf X}} =P_j - a_j

This is a momentum translation element, generated by a momentum translation generator X{\bf X}.

Finally:

U=exp⁡(−iℏθnkLk)U = \exp\left(-\frac{i}{\hbar} \theta n^k L_k\right)

is an operator of rotation, due to Uψ(x)=ψ(R−1x)U\psi(x) = \psi(R^{-1} x).

Explicit Example: Rotation Operator Around the Z-Axis

For a rotation by angle θ\theta around the z-axis, with n=(0,0,1)\mathbf{n} = (0, 0, 1):

Angular Momentum Operator LzL_z

Lz=−iℏ(x∂∂y−y∂∂x)L_z = -i \hbar \left( x \frac{\partial}{\partial y} - y \frac{\partial}{\partial x} \right)

In cylindrical coordinates: Lz=−iℏ∂∂ϕL_z = -i \hbar \frac{\partial}{\partial \phi}

Rotation Operator UU

U=exp⁡(−iℏθLz)U = \exp\left(-\frac{i}{\hbar} \theta L_z\right)

Explicitly:

Action on Wave Function

(Uψ)(x,y,z)=ψ(xcos⁡θ+ysin⁡θ,−xsin⁡θ+ycos⁡θ,z)(U \psi)(x, y, z) = \psi(x \cos \theta + y \sin \theta, -x \sin \theta + y \cos \theta, z)
(Uψ)(ρ,ϕ,z)=ψ(ρ,ϕ−θ,z)(U \psi)(\rho, \phi, z) = \psi(\rho, \phi - \theta, z)

Matches ψ(R−1x)\psi(R^{-1} \mathbf{x}), where:

R=(cos⁡θ−sin⁡θ0sin⁡θcos⁡θ0001)R = \begin{pmatrix} \cos \theta & -\sin \theta & 0 \\ \sin \theta & \cos \theta & 0 \\ 0 & 0 & 1 \end{pmatrix}