# Tensor operator

"Spherical tensor operator" redirects to here. For the closely related concept see spherical basis.

In pure and applied mathematics, particularly quantum mechanics and computer graphics and applications therefrom, a tensor operator generalizes the notion of operators which are scalars and vectors. A special class of these are spherical tensor operators which apply the notion of the spherical basis and spherical harmonics. The spherical basis closely relates to the description of angular momentum in quantum mechanics and spherical harmonic functions.

## Rotations of quantum states

### Quantum rotation operator

The rotation operator about the unit vector n (defining the axis of rotation) through angle θ is

$U[R(\theta, \hat{\mathbf{n}})] = \exp\left(-\frac{i\theta}{\hbar}\hat{\mathbf{n}}\cdot\widehat{\mathbf{J}}\right)$

where J = (Jx, Jy, Jz) are the rotation generators (also the angular momentum matrices):

$J_x = \begin{pmatrix} 0 & 0 & 0 \\ 0 & 0 & -i \\ 0 & i & 0 \end{pmatrix}\,\quad J_y = \begin{pmatrix} 0 & 0 & i \\ 0 & 0 & 0 \\ -i & 0 & 0 \end{pmatrix}\,\quad J_z = \begin{pmatrix} 0 & -i & 0 \\ i & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}$

and let $\widehat{R} = \widehat{R}(\theta,\hat{\mathbf{n}})$ is a rotation matrix. Taylor expanding to first order in θ gives the infinitesimal rotation operator:

$U[R(\theta, \hat{\mathbf{n}})] = 1 - \frac{i\theta}{\hbar}\hat{\mathbf{n}}\cdot\widehat{\mathbf{J}}$

An operator $\widehat{\Omega}$ is invariant under a unitary transformation U if:

$\widehat{\Omega} = {U}^\dagger \widehat{\Omega} U$

in this case for the rotation $\widehat{U}(R)$:

$\widehat{\Omega} = {U(R)}^\dagger \widehat{\Omega} U(R) = \exp\left(\frac{i\theta}{\hbar}\hat{\mathbf{n}}\cdot\widehat{\mathbf{J}}\right) \widehat{\Omega} \exp\left(-\frac{i\theta}{\hbar}\hat{\mathbf{n}}\cdot\widehat{\mathbf{J}}\right)$

### Angular momentum eigenkets

The orthonormal basis set for total angular momentum is |j, m, where j is the total angular momentum quantum number and m is the magnetic angular momentum quantum number, which takes values −j, −j + 1, ..., j − 1, j. A general state

$| \psi \rangle = \sum_m |j,m\rangle$

in the space rotates to a new state |ψ by:

$|\bar{\psi} \rangle = U(R)|\psi \rangle$

Using the completeness condition:

$I = \sum_m |j , m' \rangle \langle j, m' |$

we have

$| \bar{\psi} \rangle = I U(R)|\psi \rangle = \sum_{mm'} | j , m' \rangle \langle j, m' | U(R) |j , m \rangle$

Introducing the Wigner D matrix elements:

${D(R)}^{(j)}_{m'm} = \langle j, m' | U(R) |j,m \rangle$

gives the matrix multiplication:

$| \bar{\psi} \rangle = \sum_{mm'} D^{(j)}_{m'm} | j , m \rangle \quad \Rightarrow \quad | \bar{\psi} \rangle = D^{(j)} | \psi \rangle$

For one basis ket:

$| \overline{j , m} \rangle = \sum_{mm'} {D(R)}^{(j)}_{m'm} | j , m \rangle$

For the case of orbital angular momentum, the eigenstates |l, m of the orbital angular momentum operator L and solutions of Laplace's equation on a 3d sphere are spherical harmonics:

$Y_\ell^m( \theta , \varphi ) = \langle \theta,\phi | \ell , m \rangle = \sqrt{{(2\ell+1)\over 4\pi}{(\ell-m)!\over (\ell+m)!}} \, P_\ell^m ( \cos{\theta} ) \, e^{i m \phi }$

where Pm is an associated Legendre polynomial, is the orbital angular momentum quantum number, and m is the orbital magnetic quantum number which takes the values −, − + 1, ... − 1, The formalism of spherical harmonics have wide applications in applied mathematics, and are closely related to the formalism of spherical tensors, as shown below.

Spherical harmonics are functions of the polar and azimuthal angles, θ and ϕ respectively, which can be conveniently by collected into a unit vector n(θ, ϕ) pointing in the direction of those angles, in the Cartesian basis it is:

$\hat{\mathbf{n}}(\theta,\phi) = \cos\theta \sin\phi \mathbf{e}_x + \sin\theta \sin\phi \mathbf{e}_y + \cos\phi \mathbf{e}_z$

So a spherical harmonic can also be written Ym = n|, m. Spherical harmonic states |m, rotate according to the inverse rotation matrix U(R−1), while |n rotates by the initial rotation matrix $\widehat{U}(R)$.

$| \overline{\ell , m} \rangle = \sum_{m'} D_{m'm}^{(\ell)} U(R^{-1}) | \ell , m' \rangle\,,\quad | \overline{\hat{\mathbf{n}}} \rangle = U(R) | \hat{\mathbf{n}}\rangle$

## Rotations of tensorial operators

### Scalar operators

A scalar operator is invariant under rotations:[1]

$U(R)^\dagger \widehat{S} U(R) = \widehat{S}$

and we have a simple result, that the scalar operator commutes with the rotation generators:

$\left[ \widehat{S}_a , \widehat{J}_b \right] = 0$

Examples of a scalar operators include

$\widehat{E} \psi = i\hbar\frac{\partial}{\partial t} \psi$
$\widehat{V}(\mathbf{r},t)\psi = V(\mathbf{r},t)\psi$
$\widehat{T}(\mathbf{r},t)\psi = -\frac{\hbar^2}{2m} \nabla^2 (\mathbf{r},t)$
$\widehat{\mathbf{L}}\cdot\widehat{\mathbf{S}} = \widehat{L}_x \widehat{S}_x + \widehat{L}_y \widehat{S}_y + \widehat{L}_z \widehat{S}_z \,.$

### Vector operators

Vector operators (as well as pseudovector operators) can be rotated according to:[1]

${U(R)}^\dagger \widehat{V}_i U(R) = \sum_j R_{ij} \widehat{V}_j$

from this and the infinitesimal rotation operator and its Hermitian conjugate, one can derive the commutation relation with the rotation generator:

$\left[ \widehat{V}_a , \widehat{J}_b \right] = i \hbar \varepsilon_{abc} \widehat{V}_c$

where εijk is the Levi-Civita symbol, which all vector operators must satisfy, by construction. As the symbol εijk is a pseudotensor, pseudovector operators are invariant up to a sign: +1 for proper rotations and −1 for improper rotations.

Vector operators include

$\widehat{\mathbf{r}} \psi = \mathbf{r} \psi$
$\widehat{\mathbf{p}} \psi = -i\hbar \nabla \psi$

and peusodovector operators include

$\widehat{\mathbf{L}} \psi = -i\hbar \mathbf{r} \times \nabla \psi$
• as well the spin operator S, and hence the total angular momentum
$\widehat{\mathbf{J}} = \widehat{\mathbf{L}}+\widehat{\mathbf{S}}\,.$

In Dirac notation:

$\langle\bar{\psi}|\widehat{V}_a|\bar{\psi}\rangle = \langle \psi | {U(R)}^\dagger \widehat{V}_a U(R) | \psi \rangle = \sum_b R_{ab} \langle \psi | \widehat{V}_b | \psi \rangle$

and since |ψ is any quantum state, the same result follows:

$\widehat{V}_a = {U(R)}^\dagger \widehat{V}_a U(R) = \sum_b R_{ab} \widehat{V}_b$

Note that here, the term "vector" is used two different ways: kets such as |ψ are elements of abstract Hilbert spaces, while the vector operator is defined as a quantity whose components transform in a certain way under rotations.

### Spherical vector operators

A vector operator in the spherical basis is V = (V+1, V0, V−1) where the components are:[1]

$V_{+1}=-\frac{1}{\sqrt{2}}(V_x + i V_y)\,\quad V_{-1}=\frac{1}{\sqrt{2}}(+ V_x - i V_y)\,,\quad V_0 = V_z \,,$

and the commutators with the rotation generators are:

$\left[J_z, V_q\right] = q V_q$
$\left[J_\pm, V_0\right] = \sqrt{2} V_\pm$
$\left[J_\pm, V_\mp\right] = \sqrt{2} V_0$
$\left[J_\pm, V_\pm\right] = 0$

where q is a placeholder for the spherical basis labels (+1, 0, −1), and:

$J_\pm = J_x \pm i J_y \,,$

(some authors may place a factor of 1/2 on the left hand side of the equation) and raise (J+) or lower (J) the total magnetic quantum number m by one unit. In the spherical basis the generators are:

$J_{\pm 1} = \mp \frac{1}{\sqrt{2}} J_\pm \,,\quad J_0 = J_z$

The rotation transformation in the spherical basis (originally written in the Cartesian basis) is then:

${U(R)}^\dagger \widehat{V}_q U(R) = \sum_{q'} {D(R)}^{(1)}_{qq'} \widehat{V}_{q'}$

One can generalize the vector operator concept easily to tensorial operators, shown next.

### Tensor operators and their reducible and irreducible representations

A tensor operator can be rotated according to:[1]

$U(R)^\dagger \widehat{T}_{pqr\cdots} U(R) = R_{pi}R_{qj}R_{rk}\cdots \widehat{T}_{ijk\cdots}$

Consider a dyadic tensor with components T = aibj, this rotates infinitesimally according to:

$U(R)^\dagger \widehat{T}_{pq} U(R) = R_{pi} R_{qj} \widehat{T}_{ij} = R_{pi} \widehat{a}_i R_{qj} \widehat{b}_j$

Cartesian dyadic tensors of the form

$\hat{\mathbf{T}} = \mathbf{e}_i \widehat{a}_i \otimes \mathbf{e}_j \widehat{b}_j = \mathbf{e}_i \otimes \mathbf{e}_j \widehat{a}_i \widehat{b}_j$

where a and b are two vector operators:

$\hat{\mathbf{a}} = \mathbf{e}_i \widehat{a}_i \,,\quad \hat{\mathbf{b}} = \mathbf{e}_j \widehat{b}_j$

are reducible, which means they can be re-expressed in terms of a and b as a rank 0 tensor (scalar), plus a rank 1 tensor (an antisymmetric tensor), plus a rank 2 tensor (a symmetric tensor with zero trace):

$\mathbf{T} = \mathbf{T}^{(1)} + \mathbf{T}^{(2)} + \mathbf{T}^{(3)}$

where the first term

$\widehat{T}^{(1)}_{ij} = \frac{\widehat{a}_k \widehat{b}_k}{3}\delta_{ij}$

includes just one component, a scalar equivalently written (a·b)/3, the second

$\widehat{T}^{(2)}_{ij} = \frac{1}{2}[\widehat{a}_i \widehat{b}_j - \widehat{a}_j \widehat{b}_i] = \widehat{a}_{[i} \widehat{b}_{j]}$

includes three independent components, equivalently the components of (a×b)/2, and the third

$\widehat{T}^{(3)}_{ij} = \frac{1}{2}(\widehat{a}_i \widehat{b}_j + \widehat{a}_j \widehat{b}_i) - \frac{\widehat{a}_k \widehat{b}_k}{3}\delta_{ij} = \widehat{a}_{(i} \widehat{b}_{j)} - T^{(1)}_{ij}$

includes five independent components. Throughout, δij is the Kronecker delta, the components of the identity matrix. The number in the superscripted brackets denotes the tensor rank. These three terms are irreducible, which means they cannot be decomposed further and still be tensors satisfying the defining transformation laws under which they must be invariant. These also correspond to the number of spherical harmonic functions 2 + 1 for = 0, 1, 2, the same as the ranks for each tensor. Each of the irreducible representations T(1), T(2) ... transform like angular momentum eigenstates according to the number of independent components.

### Spherical tensor operators

Continuing the previous example of the second order dyadic tensor T = ab, casting each of a and b into the spherical basis and substituting into T gives the spherical tensor operators of the second order, which are:

$\widehat{T}^{(2)}_{\pm 2} = \widehat{a}_{\pm 1} \widehat{b}_{\pm 1}$
$\widehat{T}^{(2)}_{\pm 1} = \frac{1}{\sqrt{2}}\left( \widehat{a}_{\pm 1} \widehat{b}_0 + \widehat{a}_0 \widehat{b}_{\pm 1} \right)$
$\widehat{T}^{(2)}_{0} = \frac{1}{\sqrt{6}}\left( \widehat{a}_{+1} \widehat{b}_{-1} + \widehat{a}_{-1} \widehat{b}_{+1} + 2 \widehat{a}_0 \widehat{b}_0 \right)$

Using the infinitesimal rotation operator and its Hermitian conjugate, one can derive the commutation relation in the spherical basis:

$\left[J_a, \widehat{T}^{(2)}_{q} \right] = \sum_{q'} {D(J)}^{(2)}_{qq'} \widehat{T}_{q'}^{(2)}$

and the finite rotation transformation in the spherical basis is:

${U(R)}^\dagger \widehat{T}^{(2)}_q U(R) = \sum_{q'} {D(R)}^{(2)}_{qq'} \widehat{T}_{q'}^{(2)}$

In general, tensor operators can be constructed from two perspectives.[2]

One way is to specify how spherical tensors transform under a physical rotation - a group theoretical definition. A rotated angular momentum eigenstate can be decomposed into a linear combination of the initial eigenstates: the coefficients in the linear combination consist of Wigner rotation matrix entries. Spherical tensor operators are sometimes defined as the set of operators that transform just like the eigenkets under a rotation.

A spherical tensor Tq(k) of rank k is defined to rotate into Tq(k) according to:

${U(R)}^\dagger \widehat{T}_{q'}^{(k)} U(R) = \sum_{q} D^{(k)}_{qq'} \widehat{T}_{q}^{(k)}$

where q = k, k − 1, ..., −k + 1, −k. For spherical tensors, k and q are analogous labels to and m respectively, for spherical harmonics. Some authors write Tkq instead of Tq(k), with or without the brackets enclosing the rank number k.

Another related procedure requires that the spherical tensors satisfy certain commutation relations with respect to the rotation generators Jx, Jy, Jz - an algebraic definition.

The commutation relations of the angular momentum components with the tensor operators are:

$\left[ J_\pm , \widehat{T}_{k}^{q} \right] = \hbar \sqrt{( k \mp q )(k \pm q + 1)} \widehat{T}_{k}^{q \pm 1}$
$\left[ J_z , \widehat{T}_{k}^{q} \right] = \hbar q \widehat{T}_{k}^{q}$

Consider any 3d vector, not just a unit vector, and not just the position vector:

$\mathbf{a} = a_x \mathbf{e}_x + a_y \mathbf{e}_y + a_z \mathbf{e}_z$

then a spherical tensor is a spherical harmonic as a function of this vector a, and in Dirac notation:

$T_q^{(k)} = Y_{\ell = k }^{m = q} (\mathbf{a}) = \langle \mathbf{a}|k,q\rangle$

(the super and subscripts switch places for the corresponding labels k and mq which spherical tensors and spherical harmonics use).

Spherical harmonic states and spherical tensors can also be constructed out of the Clebsch–Gordan coefficients. Irreducible spherical tensors can build higher rank spherical tensors; if Aq1(k1) and Bq2(k2) are two spherical tensors of ranks k1 and k2 respectively, then:

$T_{q}^{(k)} = \sum_q \langle k_1 , k_2 , q_1 , q_2 | k , q , q_1 , q_2 \rangle A_{q_1}^{(k_1)} B_{q_2}^{(k_2)}$

is a spherical tensor of rank k.

## Angular momentum and spherical harmonics

### Orbital angular momentum and spherical harmonics

Orbital angular momentum operators have the ladder operators:

$L_\pm = L_x \pm i L_y$

which raise or lower the orbital magnetic quantum number m by one unit. This has almost exactly the same form as the spherical basis, aside from constant multiplicative factors.

### Spherical tensor operators and quantum spin

Spherical tensors can also be formed from algebraic combinations of the spin operators Sx, Sy, Sz, as matrices, for a spin system with total quantum number j = + s (and = 0). Spin operators have the ladder operators:

$S_\pm = S_x \pm i S_y$

which raise or lower the spin magnetic quantum number ms by one unit.

## Applications

Spherical bases have broad applications in pure and applied mathematics and physical sciences where spherical geometries occur.

### Magnetic resonance

The spherical tensor formalism provides a common platform for treating coherence and relaxation in nuclear magnetic resonance. In NMR and EPR, spherical tensor operators are employed to express the quantum dynamics of particle spin, by means of an equation of motion for the density matrix entries, or to formulate dynamics in terms of an equation of motion in Liouville space. The Liouville space equation of motion governs the observable averages of spin variables. When relaxation is formulated using a spherical tensor basis in Liouville space, insight is gained because the relaxation matrix exhibits the cross-relaxation of spin observables directly.[2]

## References

### Notes

1. ^ a b c d E. Abers (2004). "5". Quantum Mechanics. Addison Wesley. ISBN 978-0-13-146100-0.
2. ^ a b R.D. Nielsen, B.H. Robinson (2006). "The Spherical Tensor Formalism Applied to Relaxation in Magnetic Resonance". p. 270-271.