Product operator (NMR)

In the various fields of nuclear magnetic resonance, the product operator mathematical formalism is often used to simplify both the design and the interpretation of often very complex sequences of radio frequency electromagnetic pulses applied to samples under study. Basically, it is a short hand mathematical construct, a set of equations, that is used in place of more complex, although equivalent, matrix multiplication. The formalism uses a rotating frame of reference, in which the central irradiation frequency, say 800 MHz, is fixed on the X- or Y-axis, and the magnetic field, by convention, points towards the postive Z-axis. By convention, I and S indicate magnetic vectors associated with protons or heteroatom, respectively. Subscripts are used to indicate the axial orientation of the magnetic vector. At equilibrium, the net proton magnetic vector is thus I_z. Although even 2 pulse experiments on only 2 protons can lead to equations with 32 parts in the final equation due to chemical shifts, pulses and various forms of coupling, the use and knowledge of phase-cycling techniques or gradients to wipe out most of the terms leads to simplified final equations once accounted for.

Single Pulses (rotations)

Arbitrary pulses (rotations)

I_x -->($\theta$_x) --> I_x

I_x -->($\theta$_y) --> I_xcos($\theta$) + I_zsin($\theta$)

I_x -->($\theta$_z) --> I_xcos($\theta$) - I_ysin($\theta$)

I_y -->($\theta$_x) --> I_ycos($\theta$) - I_zsin($\theta$)

I_y -->($\theta$_y) --> I_y

I_y -->($\theta$_z) --> I_ycos($\theta$) + I_xsin($\theta$)

I_z -->($\theta$_x) --> I_zcos($\theta$) + I_ysin($\theta$)

I_z -->($\theta$_y) --> I_zcos($\theta$) - I_xsin($\theta$)

I_z -->($\theta$_z) --> I_z

90 degree pulses

So called 90 degree ($\pi$/2) pulses, in which magnetization is rotated from one axis to another, are the most widely used single pulses in NMR spectroscopy and the above equations simplify to the following for such pulses.

I_x -->(90_x) --> I_x

I_x -->(90_y) --> I_z

I_x -->(90_z) --> -I_y

I_y -->(90_y) --> I_y

I_y -->(90_z) --> I_x

I_y -->(90_x) --> -I_z

I_z -->(90_z) --> I_z

I_z -->(90_x) --> -I_y

I_z -->(90_y) --> -I_x

Chemical Shift Operators

Nuclei rotate around the XY plane at different frequencies. For example, assuming an 800 MHz central proton frequency, some protons will rotate 800 Hertz, or 1 part-per-million (ppm), faster, while others will rotate about the field more slowly. This difference from the central frequency, expressed in ppm, is called a chemical shift, which is sybolized as $\delta$. The actual frequency, is sybolized as $\omega$ in radians/second or $\nu$ if expressed in Hertz. $\omega$= 2$\pi \nu$. Remembering that the central frequency is fixed on the X-axis, the chemical shifts of each proton will cause them to rotate away from the X-axis towards the Y-axis for faster frequencies and towards the minus Y-axis for slower frequencies. The total angle of the rotation is time dependent, so that during time delay $\tau$, the angle extended is

$\Omega$ = $\omega \tau$

Note that chemical shifts only evolved in the XY plane.

I_x -->{$\delta (\tau )$} --> I_x cos($\omega \tau$) + I_y sin($\omega \tau$)

I_y -->{$\delta (\tau )$} --> I_y cos($\omega \tau$) - I_y sin($\omega \tau$)

I_z -->{$\delta (\tau )$} --> I_z

The one pulse experiment

Ignoring relaxation, coupling and other effects then, for a simple single proton starting at equilibrium, excited with a -90_y pulse, the time-dependent signal observed is:

I_z --> I_x --> I_x cos($\omega \tau$) + I_y sin($\omega \tau$)

Relaxation Operators

Although the previous equations imply that the NMR signal would ring out indefinitely, a number of relaxation processes cause the excited NMR states to relax back to equilibrium. Although these effects can be separated into rates R₁ and R₂, longitudinal (Z-axis) and transverse (XY-plane) effects, respectively, for simplicity one can consider the effect relaxation rate, R, and its characteristic relaxation time T = 1/R. Relaxation effects cause an exponential decay of the observable signal. The effects of relaxation then are expressed as:

I_x -->{$\delta (\tau )$} --> I_x cos($\omega \tau$) + I_y sin($\omega \tau$) --> {R($\tau$)} --> I_x cos($\omega \tau$)e^{-$\tau$/T} + I_y sin($\omega \tau$)e^{-$\tau$/T}