In pure and applied mathematics, quantum mechanics and computer graphics, 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. The coordinate-free generalization of a tensor operator is known as a representation operator.[1]
In quantum mechanics, physical observables that are scalars, vectors, and tensors, must be represented by scalar, vector, and tensor operators, respectively. Whether something is a scalar, vector, or tensor depends on how it is viewed by two observers whose coordinate frames are related to each other by a rotation. Alternatively, one may ask how, for a single observer, a physical quantity transforms if the state of the system is rotated. Consider, for example, a system consisting of a molecule of mass
M
p{\hatz
z
90\circ
y
p{\hatx
x
p2/2M
p{\hatz
p{\hatx
In the same way, tensor quantities must be represented by tensor operators. An example of a tensor quantity (of rank two) is the electrical quadrupole moment of the above molecule. Likewise, the octupole and hexadecapole moments would be tensors of rank three and four, respectively.
Other examples of scalar operators are the total energy operator (more commonly called the Hamiltonian), the potential energy, and the dipole-dipole interaction energy of two atoms. Examples of vector operators are the momentum, the position, the orbital angular momentum,
{L}
{S}
Scalar, vector and tensor operators can also be formed by products of operators. For example, the scalar product
{L} ⋅ {S}
{L}
{S}
Here, the indices
i
j
x
y
z
\alpha
q\alpha
\alpha
r\alpha,
i
r\alpha,r\alpha,
{r}\alpha
The rotation operator about the unit vector n (defining the axis of rotation) through angle θ is
where are the rotation generators (also the angular momentum matrices):
and let
\widehat{R}=\widehat{R}(\theta,\hat{n
An operator
\widehat{\Omega}
\widehat{U}(R)
The orthonormal basis set for total angular momentum is
|j,m\rangle
rotates to a new state by:
Using the completeness condition:
we have
Introducing the Wigner D matrix elements:
gives the matrix multiplication:
For one basis ket:
For the case of orbital angular momentum, the eigenstates
|\ell,m\rangle
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 collected into a unit vector n(θ, ϕ) pointing in the direction of those angles, in the Cartesian basis it is:
So a spherical harmonic can also be written
m | |
Y | |
\ell |
=\langlen|\ellm\rangle
|m,\ell\rangle
U(R-1)
|\ell,m\rangle
\widehat{U}(R)
We define the Rotation of an operator by requiring that the expectation value of the original operator
\widehat{A
Now as,
we have,
since,
|\psi\rangle
A scalar operator is invariant under rotations:[2]
This is equivalent to saying a scalar operator commutes with the rotation generators:
Examples of scalar operators include
Vector operators (as well as pseudovector operators) are a set of 3 operators that can be rotated according to:[2]
Any observable vector quantity of a quantum mechanical system should be invariant of the choice of frame of reference. The transformation of expectation value vector which applies for any wavefunction, ensures the above equality. In Dirac notation:where the RHS is due to the rotation transformation acting on the vector formed by expectation values. Since is any quantum state, the same result follows: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.
From the above relation for infinitesimal rotations and the Baker Hausdorff lemma, by equating coefficients of order
\delta\theta
where εijk is the Levi-Civita symbol, which all vector operators must satisfy, by construction. The above commutator rule can also be used as an alternative definition for vector operators which can be shown by using the Baker Hausdorff lemma. As the symbol εijk is a pseudotensor, pseudovector operators are invariant up to a sign: +1 for proper rotations and −1 for improper rotations.
Since operators can be shown to form a vector operator by their commutation relation with angular momentum components (which are generators of rotation), its examples include:
and peusodovector operators include
If
\vec{V}
\vec{W}
\vec{V} ⋅ \vec{W}=
3 | |
\sum | |
i=1 |
\hat{Vi}\hat{Wi}
Under rotation of coordinates, the newly defined operator transforms as:Rearranging terms and using transpose of rotation matrix as its inverse property:Where the RHS is the
\vec{V} ⋅ \vec{W}
A vector operator in the spherical basis is where the components are:[2]
using the various commutators with the rotation generators and ladder operators are:
which are of similar form of
In the spherical basis, the generators of rotation are:
From the transformation of operators and Baker Hausdorff lemma:
{U(R)}\dagger\widehat{V}qU(R)=\widehat{V}q+i
\theta | |
\hbar |
\left[\hat{n} ⋅ \vec{J},\widehat{V}q\right]+\sum
infty | |
k=2 |
| ||||||
⋅ |
\vec{J},
k}{k!}\widehat{V} | |
.]\right) | |
q |
=exp\left({i
\theta | |
\hbar |
\hat{n} ⋅ Ad\vec{J
compared to
U(R)|j,k\rangle=|j,k\rangle-i
\theta | |
\hbar |
\hat{n} ⋅ \vec{J}|j,k
infty | |
\rangle+\sum | |
k=2 |
| |||||
⋅ |
\vec{J}\right)k}{k!}|j,k\rangle=exp\left({-i
\theta | |
\hbar |
\hat{n} ⋅ \vec{J}}\right)|j,k\rangle
it can be argued that the commutator with operator replaces the action of operator on state for transformations of operators as compared with that of states:
U(R)|j,k\rangle=\exp\left({-i
\theta | |
\hbar |
\hat{n} ⋅ \vec{J}}\right)|j,k\rangle=\sumj',k'|j',k'\rangle\langlej',k'|\exp\left({-i
\theta | |
\hbar |
\hat{n} ⋅ \vec{J}}\right)|j,k\rangle=\sumk'
(j) | |
D | |
k'k |
(R)|j,k'\rangle
The rotation transformation in the spherical basis (originally written in the Cartesian basis) is then, due to similarity of commutation and operator shown above:
One can generalize the vector operator concept easily to tensorial operators, shown next.
In general, a tensor operator is one that transforms according to a tensor:where the basis are transformed by
R-1
R
In the subsequent discussion surrounding tensor operators, the index notation regarding covariant/contravariant behavior is ignored entirely. Instead, contravariant components is implied by context. Hence for an n times contravariant tensor:
Note: In general, a tensor operator cannot be written as the tensor product of other tensor operators as given in the above example.
If
\vec{V}
\vec{W}
\hat{T}ij=\hat{Vi}\hat{Wj}
R-1
R
We observe that the subspace spanned by linear combinations of the rank two tensor components form an invariant subspace, ie. the subspace does not change under rotation since the transformed components itself is a linear combination of the tensor components. However, this subspace is not irreducible ie. it can be further divided into invariant subspaces under rotation. Otherwise, the subspace is called reducible. In other words, there exists specific sets of different linear combinations of the components such that they transforms into a linear combination of the same set under rotation.[3] In the above example, we will show that the 9 independent tensor components can be divided into a set of 1, 3 and 5 combination of operators that each form irreducible invariant subspaces.
The subspace spanned by
\{\hat{T}\}
\{\hat{A}\}
\{\hat{S}\}
\hat{A}=
1 | |
2 |
(\hat{T}-\hat{T}ji)
\hat{S}=
1 | |
2 |
(\hat{T}+\hat{T}ji)
\{\hat{T}\}
\{\hat{A}\}
\{\hat{S}\}
\{\hat{A}\}
\{\hat{S}\}
The six independent symmetric component set can be divided into five independent traceless symmetric component and the invariant trace can be its own subspace.
Hence, the invariant subspaces of
\{\hat{T}\}
\hat{t}=
3\hat{T} | |
\sum | |
kk |
\hat{A}=
1 | |
2 |
(\hat{T}-\hat{T}ji)
\hat{S}=
1 | |
2 |
(\hat{T}+\hat{T}ji)-
1 | |
3 |
\hat{t}\deltai
If
\hat{T}ij=\hat{Vi}\hat{Wj}
\{\hat{T}\}
\vec{V} ⋅ \vec{W}
1 | |
2 |
(\hat{V}i\hat{W}j-\hat{V}j\hat{W}i)
1 | |
2 |
(\hat{V}i\hat{W}j+\hat{V}j\hat{W}i)-
1 | |
3 |
(\vec{V} ⋅ \vec{W})\deltai
From the above examples, the nine component
\{\hat{T}\}
\{\hat{T}\}
\hat{T}=
1 | |
3 |
\hat{t}\deltaij+\hat{A}ij+\hat{S}ij
or
\hat{T}=
1 | |
3 |
(\vec{V} ⋅ \vec{W})\deltaij+\left(
1 | |
2 |
(\hat{V}i\hat{W}j-\hat{V}j\hat{W}i)\right)+\left(
1 | |
2 |
(\hat{V}i\hat{W}j+\hat{V}j\hat{W}i)-
1 | |
3 |
(\vec{V} ⋅ \vec{W})\deltai\right)=T(0)+T(1)+T(2)
where:
In general cartesian tensors of rank greater than 1 are reducible. In quantum mechanics, this particular example bears resemblance to the addition of two spin one particles where both are 3 dimensional, hence the total space being 9 dimensional, can be formed by spin 0, spin 1 and spin 2 systems each having 1 dimensional, 3 dimensional and 5 dimensional space respectively. 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. Each of the irreducible representations T(0), T(1), T(2) ... transform like angular momentum eigenstates according to the number of independent components.
It is possible that a given tensor may have one or more of these components vanish. For example, the quadrupole moment tensor is already symmetric and traceless, and hence has only 5 independent components to begin with.
See also: Wigner–Eckart theorem. Spherical tensor operators are generally defined as operators with the following transformation rule, under rotation of coordinate system:
(j) | |
\widehat{T} | |
m |
→ U(R)\dagger
(j) | |
\widehat{T} | |
m |
U(R)=\summ'
(j) | |
D | |
m'm |
(R-1
(j) | |
)\widehat{T} | |
m' |
The commutation relations can be found by expanding LHS and RHS as:
U(R)\dagger
(j) | |
\widehat{T} | |
m |
U(R)=\left(1+
i\epsilon\hat{n | |
⋅ |
\vec{J}}{\hbar}
(j) | ||
+l{O}(\epsilon | \left(1- | |
m |
i\epsilon\hat{n | |
⋅ |
\vec{J}}{\hbar}+l{O}(\epsilon2)\right)=\summ'\langlej,m'|\left(1+
i\epsilon\hat{n | |
⋅ |
\vec{J}}{\hbar}+l{O}(\epsilon2)\right)|j,m\rangle
(j) | |
\widehat{T} | |
m' |
Simplifying and applying limits to select only first order terms, we get:
{[\hat{n} ⋅
(j) | |
{\vec{J}}},\widehat{T} | |
m |
]=\summ'
(j) | |
\widehat{T} | |
m' |
\langlej,m'|\vec{J} ⋅ \hat{n}|j,m\rangle
For choices of
\hat{n}=\hat{x}\pmi\hat{y}
\hat{n}=\hat{z}
Jx
Jy
J\pm
If, only the commutation relations hold, using the following relation,
|j,m\rangle → U(R)|j,m\rangle=exp\left(-{i
\theta | |
\hbar |
\hat{n} ⋅ \vec{J}}\right)|j,m\rangle=\summ'
(j) | |
D | |
m'm |
(R)|j,m'\rangle
we find due to similarity of actions of
J
|j,m\rangle
(j) | |
\widehat{T} | |
m |
(j) | |
\widehat{T} | |
m |
→ U(R)\dagger
(j) | |
\widehat{T} | |
m |
U(R)=exp\left({i
\theta | |
\hbar |
\hat{n} ⋅ ad\vec{J
where the exponential form is given by Baker–Hausdorff lemma. Hence, the above commutation relations and the transformation property are equivalent definitions of spherical tensor operators. It can also be shown that
\{ad\hat{Ji}\}
In the following section, construction of spherical tensors will be discussed. For example, since example of spherical vector operators is shown, it can be used to construct higher order spherical tensor operators. In general, spherical tensor operators can be constructed from two perspectives. 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. Or by continuing the previous example of the second order dyadic tensor T = a ⊗ b, casting each of a and b into the spherical basis and substituting into T gives the spherical tensor operators of the second order.
Combination of two spherical tensors
(k1) | |
A | |
q1 |
(k2) | |
B | |
q2 |
This equation can be used to construct higher order spherical tensor operators, for example, second order spherical tensor operators using two first order spherical tensor operators, say A and B, discussed previously:
Using the infinitesimal rotation operator and its Hermitian conjugate, one can derive the commutation relation in the spherical basis:and the finite rotation transformation in the spherical basis can be verified:
Define an operator by its spectrum:Since for spherical harmonics under rotation:It can also been shown that:Then
m | |
\Upsilon | |
l(\vec{V}) |
\vec{V}
l | |
\Upsilon | |
l(\vec{r})=r |
m | |
Y | |
l(x,y,z) |
(l) | |
\hat{T} | |
m |
The Hermitian adjoint of a spherical tensor may be defined asThere is some arbitrariness in the choice of the phase factor: any factor containing will satisfy the commutation relations.[5] The above choice of phase has the advantages of being real and that the tensor product of two commuting Hermitian operators is still Hermitian.[6] Some authors define it with a different sign on, without the, or use only the floor of .[7]
Orbital angular momentum operators have the ladder operators:
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.
See also: Spin-weighted spherical harmonics and spinor spherical harmonics.
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:
which raise or lower the spin magnetic quantum number ms by one unit.
Spherical bases have broad applications in pure and applied mathematics and physical sciences where spherical geometries occur.
The transition amplitude is proportional to matrix elements of the dipole operator between the initial and final states. We use an electrostatic, spinless model for the atom and we consider the transition from the initial energy level Enℓ to final level En′ℓ′. These levels are degenerate, since the energy does not depend on the magnetic quantum number m or m′. The wave functions have the form,
The dipole operator is proportional to the position operator of the electron, so we must evaluate matrix elements of the form,
where, the initial state is on the right and the final one on the left. The position operator r has three components, and the initial and final levels consist of 2ℓ + 1 and 2ℓ′ + 1 degenerate states, respectively. Therefore if we wish to evaluate the intensity of a spectral line as it would be observed, we really have to evaluate 3(2ℓ′+ 1)(2ℓ+ 1) matrix elements, for example, 3×3×5 = 45 in a 3d → 2p transition. This is actually an exaggeration, as we shall see, because many of the matrix elements vanish, but there are still many non-vanishing matrix elements to be calculated.
A great simplification can be achieved by expressing the components of r, not with respect to the Cartesian basis, but with respect to the spherical basis. First we define,
Next, by inspecting a table of the Yℓm′s, we find that for ℓ = 1 we have,
where, we have multiplied each Y1m by the radius r. On the right hand side we see the spherical components rq of the position vector r. The results can be summarized by,
for q = 1, 0, −1, where q appears explicitly as a magnetic quantum number. This equation reveals a relationship between vector operators and the angular momentum value ℓ = 1, something we will have more to say about presently. Now the matrix elements become a product of a radial integral times an angular integral,
We see that all the dependence on the three magnetic quantum numbers (m′,q,m) is contained in the angular part of the integral. Moreover, the angular integral can be evaluated by the three-Yℓm formula, whereupon it becomes proportional to the Clebsch-Gordan coefficient,
The radial integral is independent of the three magnetic quantum numbers (m′, q, m), and the trick we have just used does not help us to evaluate it. But it is only one integral, and after it has been done, all the other integrals can be evaluated just by computing or looking up Clebsch–Gordan coefficients.
The selection rule m′ = q + m in the Clebsch–Gordan coefficient means that many of the integrals vanish, so we have exaggerated the total number of integrals that need to be done. But had we worked with the Cartesian components ri of r, this selection rule might not have been obvious. In any case, even with the selection rule, there may still be many nonzero integrals to be done (nine, in the case 3d → 2p).The example we have just given of simplifying the calculation of matrix elements for a dipole transition is really an application of the Wigner–Eckart theorem, which we take up later in these notes.
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.[8]