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.

Tensor operators

So far we have mostly discussed the action of transformations on quantum states. But transformations are also defined by their action on operators:

AUAU=AA \to U A U^{\dagger} = A'

In the case that the transformation is infintesimal, U=1iϵK+O(ϵ2)U = {\bf 1} - i \eps K + \cO(\eps^2) with KK Hermitian, then to O(ϵ)\cO(\eps), we have

A=Aiϵ[K,A]A' = A - i \eps [K,A]

We cna compare this to the change in classical observables under a transformation generated by an obserbale K(p,q)K(p,q): for this %\delta A = - {K, A}where where { , }$ is the Poisson bracket.

Now operators also form a vector space under addition. So we can expect them to form irreps of the transformation group. Here we will discuss such operators in the case of the rotation group.

Example: position operators

To get our heads around the concept of the transformation of operators. To see how these transform, consider the matrix element

χ(x^i)ψ=χU(R)x^iU(R)ψ=d3xχ(R1x)xiψ(R1x)\begin{align} \bra{\chi} ({\hat x}^i)' \ket{\psi} & = \bra{\chi} U(R) {\hat x}^i U^{\dagger}(R) \ket{\psi}\\ & = \int d^3 x \chi^*(R^{-1} x) x^i \psi(R^{-1} x) \end{align}

Now, we can change variables as x(Rx){\vec x} \to {\vec (Rx)}. Because RR has determinant 1, the measure factor does not change and we have```{math}

χ(x^i)ψ=d3xχ(x)Rjixjψ(x)\bra{\chi} ({\hat x}^i)' \ket{\psi} = \int d^3 x \chi^*(x) R^i_j x^j \psi(x)

This is true for any pair of states, so we have (x^i)=Rjix^j({\hat x}^i)' = R^i_j {\hat x}^j; the operator transforms the same way as a vector in physical space.

In the case of infinitesimal transformations, we can expand out RR. Instead, let us consider the simple case of an infinitesimal rotation about the zz axis. In this case, U1iϕJzU \sim {\bf 1} - \frac{i\phi}{\hbar} J_z, and so (x^i)=x^iiϕ[Jz,x^i]({\hat x}^i)' = {\hat x}^i - \frac{i\phi}{\hbar} [J_z, {\hat x}^i]. Thus

x^=x^+ϕy^y^=y^ϕy^z^=z^\begin{align} {\hat x}' & = {\hat x} + \phi {\hat y}\\ {\hat y}' & = {\hat y} - \phi {\hat y}\\ {\hat z}' & = {\hat z} \end{align}

This is precisely how a vector transforms under infinitesimal rotations about the zz axis. The first two transformations can be reassembled like so:

(x^±iy^)=(x^±iy^)(1±iϕ)({\hat x} \pm i {\hat y})' = ({\hat x} \pm i {\hat y})(1 \pm i \phi)

Under finite rotations, A=eiϕJz/AeiϕJz/A' = e^{-i \phi J_z/\hbar} A e^{i\phi J_z/\hbar}. To see how this changes with ϕ\phi, we compute

ϕA=ϕ(eiϕJz/AeiϕJz/)=i(eiϕJz/[Jz,A]eiϕJz/)\frac{\del}{\del \phi} A' = \frac{\del}{\del \phi} \left(e^{-i \phi J_z/\hbar} A e^{i\phi J_z/\hbar}\right) = \frac{- i}{\hbar} \left(e^{-i \phi J_z/\hbar} [J_z,A] e^{i\phi J_z/\hbar}\right)

Now in the present case we have

ϕ(x+iy)=yix=i(x+iy)ϕ(xiy)=i(xiy)ϕz=0\begin{align} \frac{\del}{\del\phi} (x + i y)' & = y' - i x' = i (x + i y)'\\ \frac{\del}{\del\phi} (x - i y)' & = - i (x - i y)'\\ \frac{\del}{\del\phi} z' = 0 \end{align}

Integraing these equations from ϕ=0\phi = 0, we find

(x+iy)=(x+iy)eiϕ(xiy)=(x+iy)eiϕz=z\begin{align} (x + i y)' & = (x + i y) e^{i\phi}\\ (x - i y)' & = (x + i y) e^{-i\phi}\\ z' & = z \end{align}

In other words, under rotations, x±iyx \pm i y transform as the m=±1m = \pm 1 components of a spin-1 irrep, and zz as the m=0m = 0 component. We can show that the commutators of J±J_{\pm} with x±iyx \pm i y and zz raise and lower mm as expected.

Spherical tensor operators

More generally, we denote the irreducible spherical tensor operators Omj\cO^j_m for fixed jj as operators which transform as a spin-jj irrep under commutators with J{\vec J}. In particular, this means

[Jz,Omj]=mOmj[J±,Omj]=j(j+1)m(m±1)Om±1j\begin{align} [J_z, \cO^j_m] & = \hbar m \cO^j_m\\ [J_{\pm}, \cO^j_m] & = \hbar \sqrt{j(j+1) - m(m\pm 1)} \cO^j_{m\pm 1} \end{align}

Under finite rotations, we have

U(R)OmjU(R)=dm,mjOmjU(R) \cO^j_m U(R)^{\dagger} = d^j_{m,m'} \cO^j_{m'}

so the operators transform just as the j,m\ket{j,m} states do.

The simples tensor operators have specific names:

  1. An operator O00\cO^0_0 is called a scalar operator; it does not change under rotations. Thee include operators such as p^2{\hat {\vec p}}^2, J2{\vec J}^2, and so on.

  2. The operators O±121/2\cO^{1/2}_{\pm \half} are called spinor operators. Off the tiop of my head the examples I can think of appear in quantum field theory.

  3. The operators O±11\cO^1_{\pm 1}, O01\cO^1_0 are called vector operators. We can find linear combinations which transform the same way as a vector in R3\CR^3. Examples as x^{\hat {\vec x}}, p^{\hat {\vec p}},
    J^{\hat {\vec J}}.

Relationship to cartesian tensors

In general when people say “tensors” they have in mind objects with lots of indices that run form 1,d1,\ldots d (where dd is the dimension fo spacetime). One example is the totally antisymmetric tensor ϵijk\eps_{ijk}.

In general, a tensor in R3\CR^3 is an object of the form

Ti1,,imj1,,jn ,T^{i_1,\ldots,i_m}{}_{j_1,\ldots,j_n}\ ,

where ik,jk{1,,d}i_k, j_k \in \{1,\ldots,d\}, that transforms as follows under rotations RR:

Ti1,,imj1,,jnRi1i1Eimim(R1)j1j1(R1)jmjmTi1,,imj1,,jnT^{i_1,\ldots,i_m}{}_{j_1,\ldots,j_n} \to R^{i_1}_{i'_1}\ldots E^{i_m}_{i'_m} (R^{-1})^{j'_1}_{j_1}\ldots (R^{-1})^{j'_m}_{j_m}T^{i'_1,\ldots,i'_m}{}_{j'_1,\ldots,j'_n}

The upper indices are called contravariant indices and the lower indices covariant indices. The reason for this nomenclature is based on how the transform relative to the gradient operator xi\frac{\del}{\del x^i}; the chain rule shows that

xi(R1)ikxk\frac{\del}{\del x^i} \to (R^{-1})^k_i \frac{\del}{\del x^k}

Thus the upper indices transform oppositely to this (thus contravariant), and the lower indices transform in the same way (thus covariant).

Note that the coordinates xix^i are thus contravariant vectors. As it happens, the momenta pip_i are covariant vectors. This is compatible with their representation in quantum mechanics p^i=ixi{\hat p}_i = \frac{\hbar}{i} \frac{\del}{\del x^i}.

Such tensors form reducible representations of the rotation group. We can see they behave like tensor products of the spin-1 irrep. From these, we can try to work out combinations that transform as irreps. The general theory is complicated; we will work out a simple example.

Consider a 2-index contravariant tensor TijT^{ij}. This has 9 total components. Since under the rules for tensor products of irreps of the rotation group, D1D1=D2D1D0D_1\otimes D_1 = D_2 \oplus D_1 \oplus D_0, we expect TijT^{ij} to be written in terms of a scalatr operator, a vector operator, and a j=2j = 2 tensor operator. To see this note that the total angular omentum in D1D2D_1\otimes D_2 is invariant under exchange of these two factors. For the tensor, this corresponds to the exchange of the two indices. Therefore, we expect a given irrep to correspond to tensors which come back to themselves up to ±1\pm 1 under the exchange of indices.

First, we can form T = \frac{1]{3} \delta_{ij} T^{ij}. Orthogonality of the rotation matrices guarantees that this is invariant under rotations; so this is a scalar operator.

Next, we can consider Tasij=12(TijTji)T_{as}^{ij} = \half \left(T^{ij} - T^{ji}\right). This has three components. The following object

Ti=12ϵijkTjkT_i = \frac{1}{2} \epsilon_{ijk} T^{jk}

depends on the same three components; it transforms as a (covariant) vector. So the antisymmetric part of TT is a vector opwerator.

Finally, this leaves Tsij=12(Tij+Tji)δijTT_s^{ij} = \half(T^{ij} + T^{ji}) - \delta_{ij} T. This is traceless and symmetric. By removing the trace, the resulting symmetric tensor has 5 independent parameters. Thus the components transform as operators of the form Om2\cO^2_m, m{2,1,0,1,2}m \in \{-2,1,0,1,2\}.

Action on states and the Wigner-Eckart theorem

Let us consider a spherical temsor operator Om1j1\cO^{j_1}_{m_1} acting on a state α,j2m2\ket{\alpha, j_2 m_2}. where α\alpha denote sadditional quantum numbers. Under finite rotations, this transforms as

Om1j1α,j2m2U(R)Om1j1α,j2m2=U(R)Om1j1U(R)U(R)α,j2m2=dm1,m1j1Om1j1dm2,m2j2α,j2m2\begin{align} \cO^{j_1}_{m_1}\ket{\alpha, j_2 m_2} & \to U(R) \cO^{j_1}_{m_1}\ket{\alpha, j_2 m_2}\\ & = U(R) \cO^{j_1}_{m_1}U(R)^{\dagger} U(R) \ket{\alpha, j_2 m_2}\\ & = d^{j_1}_{m_1,m_1'} \cO^{j_1}_{m_1'} d^{j_2}_{m_2,m_2'}\ket{\alpha, j_2 m_2'} \end{align}

But this behaves exactly like s state in the tensor product Dj1Dj2D_{j_1}\otimes D_{j_2}. Furthermore, if we act on this state with the operator JzJ_z, using [Jz,Omj]=mOmj[J_z,\cO^j_m] = \hbar m \cO^j_m, we have

JzOm1j1α,j2m2=(m1+m2)Om1j1α,j2m2J_z \cO^{j_1}_{m_1}\ket{\alpha, j_2 m_2} = \hbar(m_1 + m_2) \cO^{j_1}_{m_1}\ket{\alpha, j_2 m_2}

Thus,

Om1j1j2m2,α=j=j1j2j=j1+j2cjββ,j,m1+m2\cO^{j_1}_{m_1}\ket{j_2 m_2,\alpha} = \oplus_{j = |j_1 - j_2|}^{j = j_1 + j_2} c_j^{\beta} \ket{\beta, j, m_1 + m_2}

This leads to the following important theorem, the Wigner-Eckart theorem. It comes from taking matrix elements of the commutators of Jz,J±J_z, J_{\pm} with O\cO. We will not prove it here, but simply state it:

α1,j1,m1Om2j2α3,j3,m3=α1,j1Oj2α)3,j3j1,m2;j2,m2j3,m3\bra{\alpha_1, j_1, m_1} \cO^{j_2}_{m_2} \ket{\alpha_3, j_3, m_3} = \langle \alpha_1, j_1 || \cO^{j_2} || \alpha)3, j_3 \rangle \brket{j_1,m_2; j_2,m_2}{j_3,m_3}

Here j1,m2;j2,m2j3,m3\brket{j_1,m_2; j_2,m_2}{j_3,m_3} is a Clebsch-Gordon coefficient for the expansion of Dj1Dj2D_{j_1}\otimes D_{j_2} into irreps with angular momentum j3j_3; while α1,j1Oj2α3,j3 \langle \alpha_1, j_1 || \cO^{j_2} || \alpha_3, j_3 \rangle is a reduced matrix element. The essential point is that it is independent of mkm_k; all of the mkm_k dependence, and the information about which j3j_3 appear at all, are contained in the Clebsch-Gordon coefficients. This is kinematic information, goverened entirely by the properties under rotation. The dynamical informatio, dependent on the detailed nature of the system, is contained in the reduced matrix element. To calculate this, we need merely compute the matrix element for specific set of mkm_ks. We can then divide by the Clebsch-Gordon coefficient for that set, to extract the reduced matrix element; all other matrix elements as above are then determined by the reduced matrix element and the Clebsch-Gordon coefficients which are known (and can be looked up).