<<

4.1 Matrices

The of vector is the realm of . All rotations can be described by the multiplication of matrices. Although we will not use matrix multiplications very often in describing the rotations of NMR, the methods will serve as a useful model to introduce and visualize some rather bizarre ideas that we will encounter in coupled systems. Matrix methods are also very easy to program on a computer and can be of great assistance in analyzing complicated pulse sequences.

( 10 0 a Ix N & N ( 0 cos sin b Iy 0 sinN cosN ( (1) c Iz NIIˆ Yx I Vector

This is the matrix representation for a rotation. The one shown is for a rotation of an N about the x axis. The matrix with 3 rows and 3 columns (3 X 3) is the rotation matrix and it "operates" on the 3 X 1 vector matrix which represents the magnetization vector.

To multiply a matrix and a vector, first the top row of the matrix is multiplied element by element with the column vector, then the sum of the products becomes the top element in the resultant vector. The next row times the column vector gives the middle element of the resultant and likewise for the third.

( % ( % ( abcIx aIx bIy cIz ' ( % ( %( def Iy dIx eIy fIz (2) gh i ( % ( %( Iz gIx hIy iIz

E B If a vector, I, has the components [0*Ixyz ,0*I ,1*I ] or Iz, then the a rotation of 90 Ix ( /2 radian) be represented as: 10 0 0 ? 0 cos(90E) &sin(90E) 0 ' ? 0 sin(90E) cos(90E) 1 ? (3)

90EIˆ ' Y x Iz which since cos(90E) = 0 and sin(90E) = 1 gives:

10 0 0 0 00&1 0'&1 (4) 01 0 1 0

The resulting vector is [0*Ixyz,-1*I ,0*I ] or -Iy. In standard product notation, B Iz = /2Îx=> -Iy The rotation matrices for rotations of a three dimensional vector around the three coordinate axes are: 10 0 cosN 0 sinN cosT &sinT 0 0 cosN &sinN 010 sinT cosT 0 0 sinN cosN &sinN 0 cosN 001 (5)

NIIˆ NIIˆ TIIˆ Yx Yy Yz

In most sequences there are more than one rotation. As an example, consider a E T sequence that rotates Iz by 90 about the X axis followed by a rotation of t about the Z axis. In standard notation:

B T Iz = /2Îx=> = IztÎ => ?

In matrix notation:

cosTIt &sinTIt 0 10 0 0 ? sinTIt cosTIt 0 (0 cos90E &sin90E ( 0 ' ? 0010 sin90E cos90E 1 ? (6)

TIIˆ NIIˆ Yz Yx Iz

Note that the order of time sequential matrices is from right to left. The first rotation is next to the vector and the next rotation is placed to the left. The order of rotations is very important. Once the physical order of matrices is established, the order of multiplication is irrelevant. The rightmost matrix can multiply the vector yielding a resultant vector which is then multiplied by the next matrix to give the final vector. Alternatively, the two matrices can be multiplied together first and then the resultant matrix can multiply the vector. Both results will be the same. However, transposing the matricies will, in general, completely change the outcome. To multiply two matrices together:

abc ABC aA%bD%cG aB%bE%cH aC%bF%cI def(DEF'dA%eD%fG dB%eE%fH dC%eF%fI (7) gh i GH I gA%hD%iG gB%hE%iH gC%hF%iI Using Equation 7 we can calculate the overall rotation of the sequence.

T & T cos It sin It 0 10 0 0 ? T T ( E & E ( ' sin It cos It 0 0 cos90 sin90 0 ? 0010 sin90E cos90E 1 ?

T & T cos It sin It 0 10 0 0 ? T T (& ( ' sin It cos It 0 00 1 0 ? (8) 00101 0 1 ?

T T T cos It 0 sin It 0 sin It T & T ('& T sin It 0 cos It 0 cos It 01 0 1 0

T T The result is the vector [sin( t)*Ix, -cos( t)*Iy, 0*Iz].

Using operator formalism,

B B B Iz = /2Îx=> Iz cos /2 - Iy sin /2, which is equal to

-Iy.

Since cos B/2 = 0 and sin B/2 = 1. For the second rotation:

T T T -Iy = IztÎ => -Iy cos It + Ix sin It

Which is (whew!) the same as calculated by matricies. Obviously, if the order of the rotations is reversed:

T B Iz = IztÎ => = /2Îx=>

The answer is completely different.

T B Iz = IztÎ => Iz = /2Îx=> -Iy The operator notation is essentially an expanded version of more compact matrix method. When the sequences become very complicated (reality) and there is no possibility for the luxury of simplification, then the matrix method is a very nice tool to use especially if a computer is available.

So as to not create too much confusion for those who know something about density matricies, the rotation matricies are not the . In fact, the magnetization vector is identical to the reduced spin density matrix. The rotation matricies presented here are related to unitary matricies called superoperators.

As an example how one might use the rotation matrices to simplify calculations consider a single isolated spin I subjected to the following pulse sequence: 90x 180x I $ t $$ t Ž

Calculating this sequence by product operators:

B T T T Iz = /2Îx=> -Iy = IztÎ => -Iy cos It + Ix sin It B T T = Îx=> Iy cos It + Ix sin It T T T T = IztÎ => ( Iy cos It - Ix sin It ) cos It T T T +( Ix cos It + Iy sin It ) sin It / T T Iy (cos² It + sin² It) / Iy With such a simple result, it would seem there should be a simpler way to calculate the sequence. Figure 4.1. Vector representation of the refocusing effect of a 180E pulse placed in the center of free precession period.

If we look only at the rotations:

B T B T = /2Îx=> = IztÎ => = Îx=> = IztÎ =>

We can simplify this sequence by modifying it in the following manner. By introducing a rotation followed by the inverse of the same rotation we do no cause any net of the spin system. For example:

B B = Îx=> =- Îx=> rotates the spin system by 180E then by -180E. The overall rotation is zero this is 100 010 001

IE

equivalent to multiplication by the IE. This can easily be demonstrated by multiplying two inverse matricies.

We can formally introduce these rotations anywhere in a pulse sequence. Using the above sequence:

B B B T B T = /2Îx=> = Îx=> =- Îx=> = IztÎ => = Îx=> = IztÎ =>

We have not accomplished anything. However, if we now look at the series of rotations underlined in the next sequence,

B B B T B T = /2Îx=> = Îx=> =- Îx=> = IztÎ => = Îx=> = IztÎ => , we can simplify the sequence by means of matrix algebra.

10 0 cosTt &sinTt 0 10 0 0&10 sinTt cosTt 0 0&10 00&1 00100&1 &B T B Ix tIz Ix

cosTt &sinTt 0 10 0 &sinTt &cosTt 0 0&10 00&100&1

cosTt sinTt 0 &sinTt cosTt 0 001 &T tIz The shows that the sequence:

B T B =- Îx=> = IztÎ => = Îx=> is equivalent to the single rotation:

T =- IztÎ => .

Substituting this rotation for the equivalent three rotations in the sequence:

B B T T = /2Îx=> = Îx=> =- IztÎ => = IztÎ =>

Now two inverse rotations are adjacent:

B B T T = /2Îx=> = Îx=> =- IztÎ => = IztÎ =>

And their product is equal to IE. These rotations cancel, yielding:

B B = /2Îx=> = Îx=> or simply

B =3 /2Îx=>

T We have eliminated the = IztÎ => term completely and made the calculation much easier.

B Iz =3 /2Îx=> Iy

This sequence is referred to as a spin echo sequence, and is used to eliminate the T chemical shift term, = IztÎ =>, during a period of a pulse sequence.