Rotation Operators

In order to emphasize various things I have coloured operators, eigenfunctions etc.:

eigenfunctions

operators

eigenvalues

Rotation operators are so named because when they operate on a state vector they produce a 'rotation' in the state space. As nmr spectroscopists we are generally interested in what is happening to our spin systems as time passes and, along with the density operator, rotation operators are very helpful to us. They are very easy to visualize for an object in Cartesian coordinate space. For simpliciy, let's look at the rotation of a 2D vector:

This is actually a rotation of coordinate axes rather than a vector rotation but it is entirely equivalent. Our problem is, given starting values of x0 and y0 in the coordinate axes x and y, what will the new values x1 and y1 be in the new coordinate axes x' and y'. Looking at the diagram and the red line and using a little basic trigonometry you should be able to see that:

A = cos(f)x0 and B = sin(f)y0

and since A+B=x1 we can write:

x1 = cos(f)x0 + sin(f)y0

Similarly, for y1 we find:

y1 = -sin(f)x0 + cos(f)y0

We can write this in matrix notation:

( x1 ) = ( cos(f) sin(f) ) ( x0)
( y1 )      ( -sin(f) cos(f) ) ( x0)

or even more compactly:

V' = AV

Yet another way to look at this is to characterize the vector with a complex number:

Algebraically we write:

r = x + iy

Simple trigonometry gives us:

x = |r|cos(f)
y = |r|sin(f)

and substituting into the first equation:

r = |r|cos(f) + i|r|sin(f)
= |r|(cos(f) + isin(f))

Using the Euler identity:

e(if) = cos(f) + isin(f)

we rewrite:

r = |r|e(if)

So we now have an exponential representation of our vector. What about rotation? Let's multiply |r|e(if) by e(w):

r = |r|e(if)e(w)
= |r|e(if + w)

This must represent a new vector that has phase angle f + w:

So, our original vector, at phase angle f has been rotated by w degrees counterclockwise. e(w) can be considered to have operated on the original vector, r = |r|e(if), to produce a rotation of that vector. Obviously, to produce a clockwise rotation we would use -w instead.

Rotation operators for angular momentum operators are slightly more complicated but their form is exponential as in the previous example. We start with the time-dependent Schrodinger equation:

d|Y(t)>/dt = -iH|Y(t)>

where H is the hamiltonian operator.

This is just a simple 1st order differential equation with a simple solution:

|Y(t)> = exp(-iHt)|Y(0)>

where |Y(0)> is the state of the system at time 0 and exp(-iHt) is the exponential form of the hamiltonian operator.

Now let's see what happens if we calculate the expectation value for the operator, Ix under the influence of a magnetic field. Our hamiltonian operator becomes -ghB0Iz or, in frequency units, -wIz:

<Ix> = <Y(t)|Ix|Y(t)>

= <(exp(iwIzt)Y(0))*|Ix|exp(iwIzt)Y(0)>

= <Y(0)|exp(-iwIzt)Ixexp(iwIzt)|Y(0)>

= <Y(0)|I'x|Y(0)>

where I have defined:

I'x = exp(-iwIzt)Ixexp(iwIz)

I'x is related to Ix by a rotation about the z-axis by wt=q radians. To show this we start with the following definition:

g(q) = exp(-iqIz)Ixexp(iqIz)

The first derivative of g is:

g'(q) = exp(-iqIz)(-iIzIx + iIxIz)exp(iqIz)

Remembering our angular momentum commutators:

[Iz,Ix] = iIy

we substitute:

g'(w) = exp(-iqIz)Iyexp(iqIzq)

Next, the second derivative:

g''(w) = exp(-iqIz)(-iIzIy + iIyIz)exp(iqIz)

= -exp(-iqIz)Ixexp(iqIz)

= -g(q)

or, in standard differential equation form:

g''(q) + g(q) = 0

The solution to this equation is:

g(q) = Acos(q) + Bsin(q)

A equals g(0) which is just Ix and B equals g'(0) which is Iy so:

g(q) = Ixcos(q) + Iysin(q)

So:

I'x = exp(-iqIz)Ixexp(iqIz) = Ixcos(q) + Iysin(q)

Similarly:

I'y = exp(-iqIz)Iyexp(iqIz) = -Ixsin(q) + Iycos(q)

Does this look familiar? Look back up this page to the very beginning where vectors in Cartesian space were being discussed. This equation has exactly the same form as the simple rotation equations. Evidently, an expression such as:

exp(-iqIz)Ixexp(iqIz)

corresponds to a rotation about the z-axis by q radians. Sometimes you will see this equation shortened to:

R(q)IxR(-q)