Product operator formalism explained

In NMR spectroscopy, the product operator formalism is a method used to determine the outcome of pulse sequences in a rigorous but straightforward way. With this method it is possible to predict how the bulk magnetization evolves with time under the action of pulses applied in different directions. It is a net improvement from the semi-classical vector model which is not able to predict many of the results in NMR spectroscopy and is a simplification of the complete density matrix formalism.

In this model, for a single spin, four base operators exist:

I _{x} , I _{y} , I _{z} and E/ 2 which represent respectively polarization (population difference between the two spin states), single quantum coherence (magnetization on the xy plane) and the unit operator. Many other, non-classical operators exist for coupled systems. Using this approach, the evolution of the magnetization under free precession is represented by I _{z} and corresponds to a rotation about the z-axis with a phase angle proportional to the chemical shift of the spin in question:I _{x} \xrightarrow{\omega \tau I _{z}} \cos (\omega \tau)I _{x} - \sin (\omega \tau)I _{y}

Pulses about the x and y axis can be represented by

I _{x} and I _{y} respectively; these allow to interconvert the magnetization between planes and ultimately to observe it at the end of a sequence. Since every spin will evolve differently depending on its shift, with this formalism it is possible to calculate exactly where the magnetization will end up and hence devise pulse sequences to measure the desired signal while excluding others.

The product operator formalism is particularly useful in describing experiments in two-dimensions like COSY, HSQC and HMBC.

Motivation for sets of spin-1/2 particles

\hbar = 1 for convenience.

The product operator formalism is usually applied to sets of spin-1/2 particles, since the fact that the individual operators satisfy

L ^{2} _{x} = L ^{2} _{y} = L ^{2} _{z} \propto 1 , where 1 is the identity operator, makes the commutation relations of product operators particularly simple. In principle the formalism could be extended to higher spins, but in practice the general irreducible spherical tensor treatment is more often used. As such, we consider only the spin-1/2 case below.\rho , which evolves under a Hamiltonian H according to the Liouville-von Neumann equation as(\partial \rho)/(\partial t) = - i [H, \rho]. U(t) = \exp( - i H t) as\rho(t) = U(t)\rho( 0 )U ^{-1} (t) = \exp( - i H t) \rho( 0 ) \exp( + i H t)

Density operator-state duality

Suppose a single spin-1/2

L is in the state |\uparrow\rangle , which is an eigenstate of the z-spin operator L _{z} , that is L _{z|\uparrow\rangle} = (1)/(2) |\uparrow\rangle . Similarly L _{z|\downarrow\rangle} =- (1)/(2) |\downarrow\rangle . Making use of the expansion of a Hermitian operator A in terms of projections onto its eigenkets |a\rangle with eigenvalues a as A = \sum a|a\rangle\langle a| , the associated density operator is\begin{align} \rho _{|\uparrow\rangle} & = |\uparrow\rangle\langle\uparrow|\\ & = (1)/(2) (|\uparrow\rangle\langle\uparrow| + |\downarrow\rangle\langle\downarrow|) + (1)/(2) (|\uparrow\rangle\langle\uparrow| - |\downarrow\rangle\langle\downarrow|)\\ & = (1)/(2) 1 + L _{z, \end{align}} where 1 is the identity operator. Similarly, the density operator for the state |\downarrow\rangle is\rho _{|\downarrow\rangle} = (1)/(2) 1 - L _{z} Since the spin operators L _{x,} L _{y,} L _{z} are all traceless and the expectation value of an operator A for a system with density operator \rho is \langle A\rangle = \operatorname{tr}(\rho A) , the terms proportional to the unit operator 1 do not affect the expectations of the spin operators. Additionally those parts do not evolve in time, since they trivially commute with the Hamiltonian. Therefore those terms can be ignored, and the state |\uparrow\rangle corresponds to a density operator + L _{z} , while the state |\downarrow\rangle corresponds to a density operator - L _{z} . In exactly the same manner, polarisation along the positive x-axis, that is a state |\uparrow _{x\rangle} , corresponds to a density operator + L _{x} . This idea extends naturally to multiple spins, where the states and operators are direct products of single-spin states and operators. Hence operator terms in the density operator have a direct duality with states.

In the case of two spins

L, S , the terms in the density operator (ignoring the identity on its own) can be interpreted as representingL _{z,} S _{z} - longitudinal magnetisationL _{x,} L _{y,} S _{x,} S _{y} - in-phase transverse magnetisation, which is the observable quantity in NMR.2 L _{x} S _{z,} 2 L _{y} S _{z,} 2 L _{z} S _{x,} 2 L _{z} S _{y} - anti-phase longitudinal magnetisation2 L _{z} S _{z} - longitudinal two-spin order2 L _{x} S _{x,} 2 L _{x} S _{y,} 2 L _{y} S _{x,} 2 L _{y} S _{y} - other coherences, which are more difficult to interpret, but may evolve into other termswhere eg L _{z} is a shorthand for the Kronecker product L _{z} ⊗ 1 _{S} , where 1 _{S} is the identity operator on the S spin, and similarly L _{x} S _{z} is a shorthand for L _{x} ⊗ S _{z} .

The factors of two in the 'true' two-spin operators are to allow for convenient commutation relations in this specific spin-1/2 case - see below. Note also that we could instead choose to expand the density operator in the basis

L _{z,} L _{\pm} = L _{x} \pm i L _{y} etc, where the transverse operators have been replaced with raising and lowering operators. With quadrature detection, the observable associated with an individual spin is effectively the non-Hermitian L _{\pm} , so this is sometimes more convenient.

Evolution of the density operator

Consider operators

A, B, C that obey the cyclic commutation relations\begin{align} \left.[A, B]\right. & = i C,\\ \left.[B, C]\right. & = i A,\\ \left.[C, A]\right. & = i B. \end{align} In fact only the first two relations are necessary for the following derivation, but since we are usually working with operators associated with Cartesian directions, such as the individual angular momentum operators, the third commutator follows by a symmetry argument.\hat F of an operator F (in our case, this is more formally related to the adjoint representation of the Lie algebra whose elements are A, B, C ), which acts as\hat F \bullet = [F, \bullet] In particular, for the cyclic operators, we have\begin{align} \hat B A & = [B, A] = - i C,\\ {\hat B} ^{2} A & = [B, [B, A]] = [B, - i C] = - i [B, C] = A, \end{align} and consequently for integer n \ge 0 \begin{align} {\hat B} ^{2n} A & = A,\\ {\hat B} ^{2n+1} A & = - i C. \end{align} An identity for two operators F, G is\exp(F)G\exp( - F) = \exp(\hat F) G, which can be derived by putting F \to t F where t is a scalar parameter, differentiating both sides with respect to t , and noting that both sides satisfy the same differential equation in that parameter, with the same initial condition at t =0 . In particular, for some scalar parameter \theta , we have

where the final equality follows from recognising the Taylor series for sine and cosine. Now suppose that the density operator at time zero is

\rho( 0 ) = A , and it is allowed to freely evolve under the Hamiltonian H = \alpha B where \alpha is some scalar. Using the results above, the density operator at some later time t will be given by

The interpretation of this is that although nuclear spin angular momentum itself is not connected to rotations in three-dimensional space in the same way that angular momentum is, the evolution of the density operator can be viewed as rotations in an abstract space, in which the operators

A, B, C are the generators of rotations about the axes. An example of such a set of generators is just the spin operators L _{x,} L _{y,} L _{z} themselves.

We now also introduce the 'arrow notation' typically used in NMR, which writes the general evolution given above as the shorthand

A \xrightarrow{\alpha B t} \cos\alpha t A - \sin\alpha t C .With more specific reference to the radiofrequency pulses applied during NMR experiments, a hard pulse with tip angle \theta around a direction q is written as (\theta) _{q} above the arrow and corresponds to taking B = L _{q} as the rotation generator in Equation . When there is no ambiguity, the arrow label may be omitted, or be eg text instead.

Note that a more complicated calculation has now been reduced to a simpler procedure that requires no knowledge of the underlying quantum mechanics, especially since the subspaces of cyclic operators can be tabulated in advance.

Examples

The 180°-refocussing pulse

The Hamiltonian for a single spin

L evolving under a chemical shift of angular frequency \omega isH = \omega L _{z,} which means that in an ensemble of many such spins with slightly different chemical shifts, there is a dephasing of the magnetisation in the x -y plane. Consider the pulse sequence\left( (\pi)/(2) \right) _{x} — \tau — (\pi) _{x} — \tau, where \tau is a time interval. Starting in an equilibrium state with all the polarisation along the z -axis, the evolution of an individual spin in the ensemble is\begin{align} L _{z} \xrightarrow{\left( (\pi)/(2) \right) _{x}} - L _{y} \xrightarrow{\omega L _{z} t} - L _{y\cos\omega} t + L _{x\sin\omega} t \xrightarrow{(\pi) _{x}} L _{y\cos\omega} t + L _{x\sin\omega} t \xrightarrow{\omega L _{z} t} &(L _{y\cos\omega} t - L _{x\sin\omega} t)\cos\omega t + (L _{x} \cos\omega t + L _{y} \sin\omega t)\sin\omega t\\ & = L ^{2 \omega} _{y(\cos} t + \sin ^{2\omega} t) + L _{x} ( - \sin\omega t \cos\omega t + \cos\omega t \sin\omega t)\\ & = L _{y. \end{align}} Hence this sequence refocuses the transverse magnetisation produced by the first pulse, independent of the value of the chemical shift.

As an indication of the utility of the formalism, suppose instead that we tried to reach the same result using states only and therefore the Schrödinger time evolution operators. This amounts to trying to simplify the unitary propagator

U taking the initial state |\psi _{0\rangle} to the final state |\psi\rangle as |\psi\rangle = U|\psi _{0\rangle} , where explicitlyU = \exp( - i \omega L _{z} t)\exp( - i \pi L _{x)\exp(-i\omega} L _{z} t). Essentially we want to find the propagator in the form U = \exp C , that is as a single exponential of a combination of operators, because that gives the effective Hamiltonian acting during the sequence. Since the arguments of the exponentials in the original form of the propagator do not commute, this amounts to solving a specific example of the Baker–Campbell–Hausdorff (BCH) problem. In this relatively simple case we can solve the BCH problem using the fact that U f(A) U ^{\dagger} = f(U A U ^{\dagger)} for unitary operator U , operator A and function f , as well as the mathematical similarity of the spin operators with the physical rotation generators, which allow us to write\exp( i \pi L _{x)\exp(-i\omega} L _{z} t)\exp( - i \pi L _{x)} = \exp\left[ - i \omega t\exp( i \pi L _{x)L} _{z\exp(-i\pi} L _{x)\right]} = \exp( i \omega L _{z} t) \implies \exp( - i \omega L _{z} t)\exp( - i \pi L _{x)} = \exp( - i \pi L _{x)\exp(i\omega} L _{z} t). Hence U = \exp( - i \pi L _{x)} and only the effect of the 180° pulse remains, which agrees with the product operator treatment. For larger sequences of pulses this state treatment quickly becomes even more unwieldy, unless more advanced methods such as exact effective Hamiltonian theory (which gives closed-form expressions for the entangled propagators via the Cayley–Hamilton theorem and eigendecompositions) are used.

The amplitude of a Hahn echo in an inhomogeneous magnetic field

As an extension of the refocussing pulse treated above, consider a set of two pulses with arbitrary flip angles

\alpha _{1} and \alpha _{2} , that is sequence\left(\alpha _{1\right)} _{x} — \tau — \left(\alpha _{2\right)} _{x} — \tau, where again \tau is a time interval. Liberally dropping irrelevant terms, the evolution for a single spin with offset \omega up to just after the second pulse isL _{z} \xrightarrow{\left(\alpha _{1\right)} _{x}} - L _{y} \sin\alpha _{1} + … \xrightarrow{\omega L _{z} t} - L _{y} \sin\alpha _{1} \cos\omega t + L _{x} \sin\alpha _{1} \sin\omega t + … \xrightarrow{\left(\alpha _{2\right)} _{x}} L _{x} \sin\alpha _{1} \sin\omega t - L _{y} \sin\alpha _{1} \cos\alpha _{2} \cos\omega t + … . Now consider an ensemble of spins in a magnetic field that is sufficiently inhomogeneous to completely dephase the spins in the interval between the pulses. After the second pulse, we can decompose the remaining terms into a sum of two spin populations differing only in the sign of the L _{y} term, in the sense that for an individual spin we haveL _{x} \sin\omega t - \cos\alpha _{2} L _{y} \cos\omega t = \cos ^{2 (\alpha _{2})/(2)} (L _{x\sin\omega} t - L _{y\cos\omega} t) + \sin ^{2 (\alpha _{2})/(2)} (L _{x\sin\omega} t + L _{y\cos\omega} t), where we used the identities \cos ^{2\theta} + \sin ^{2\theta} = 1 and \cos ^{2\theta} - \sin ^{2\theta} = \cos 2 \theta .

It is the spins in the new population that has been generated by the second pulse, namely the one with

+ L _{y} , that will lead to the formation of an echo after evolution for the next \tau interval. Therefore, remembering to include the \sin\alpha _{1} introduced by the first pulse, the amplitude of the resulting Hahn echo relative to that produced by an ideal 90°—180° refocussing pulse sequence is roughly\sin\alpha ^{2 (\alpha _{2})/(2)} _{1 \sin} . Note that this is not an exact result, because it considers only the refocussing of polarisation that was transverse immediately before the second pulse. In reality there will be further transverse components originating from the tipping of the longitudinal magnetisation that remained after the first pulse. However, for many tip angles, this is a good rule of thumb.

To instead arrive at this result using the state formalism, we would have had to non-trivially evaluate the rotation propagator as

U(\alpha) = \exp( - i \alpha L _{x)} = \cos (\alpha)/(2) 1 -2i \sin (\alpha)/(2) L _{x,} and then evaluate a transition probability by considering the result of applying this to a state representing polarisation in the transverse plane.

DEPT (Distortionless Enhancement by Polarisation Transfer)

DEPT (Distortionless Enhancement by Polarisation Transfer) is a pulse sequence used to distinguish between the multiplicity of hydrogen bonded to carbon, that is it can separate C, CH, CH₂ and CH₃ groups. It does this by exploiting the heteronuclear carbon-hydrogen

J -coupling and varying the tip angle of the final pulse in the sequence. The basic pulse sequence is shown below.

Under the weak coupling assumption, the chemical shift terms commute with the

J -coupling term in the Hamiltonian. Hence we can ignore the refocussed chemical shift (see) in the two intervals containing \pi -pulses, namely ( 1 )\to( 4 ) and ( 3 )\to( 6 ) , and additionally refrain from evaluating the chemical shift evolution in the last (\tau)/(2) period ( 5 )\to( 6 ) . The pulse separation time (\tau)/(2) is adjusted to the coupling strength J (with associated Hamiltonian coefficient \alpha = \pi J ) such that it satisfies(\tau)/(2) ⋅ \pi J = (\pi)/(2) \implies \tau = (1)/(J) ,because then the first term in the evolved density operator in Equation vanishes under the pure coupling evolution between the pulses.

CH

Label the hydrogen spin as

L , and the carbon spin by S . For illustrative purposes, we assume that the equilibrium state only has polarisation on the L -spin (in reality, there will also be polarisation on the S spin, with the relative populations determined by the thermal Boltzmann factors). The J -coupling Hamiltonian isH = \pi J 2 L _{z} S _{z} which gives the following evolution\begin{align} ( 0 )&: L _{z\\ (0)\to(1)&: L} _{z} \xrightarrow{\left( (\pi)/(2) \right) _{xL} } -L_y\\(1)\to(2)&:\ \xrightarrow 2 L_x S_z\\(2)\to(3)&:\ 2 L_x S_z \xrightarrow -2 L_x S_y \xrightarrow -2 L_x S_y\\(3)\to(4)&:\ \xrightarrow \\(4)\to(5)&:\ \xrightarrow 2 L_x S_y \xrightarrow 2 L_x S_y \cos\theta - 2 L_z S_y \sin\theta\\(5)\to(6)&:\ \xrightarrow 2 L_x S_y \cos\theta + S_x\sin\theta\end

The non-trivial commutators used to identify the cyclic subspace for

( 1 ) \to ( 2 ) are\begin{align} \left.[L _{y,} 2 L _{z} S _{z]\right.} & = 2 (L _{y} ⊗ 1 _{L)(L} _{z} ⊗ S _{z)} - 2 (L _{z} ⊗ S _{z)(L} _{y} ⊗ 1 _{L)\\ &=} 2 L _{y} L _{z} ⊗ S _{z} - 2 L _{z} L _{y} ⊗ S _{z\\ &=} 2 [L _{y,} L _{z]} ⊗ S _{z\\ &=} i 2 L _{x} S _{z, \end{align}} and consequently the next cyclic rotation\begin{align} \left.[ 2 L _{z} S _{z,} 2 L _{x} S _{z]\right.} & = 4 [L _{z,} L _{x] ⊗} S ^{2 \\ & =} _{z} 4 i L _{y} ⊗ (1)/(4) 1 _{S\\ &=} i L _{y, \end{align}} where we used the 'mixed-product identity' (A ⊗ B)(C ⊗ D) = AC ⊗ BD , which relates the matrix and Kronecker products for compatible dimensions of A,B,C,D , and also the fact that since the two eigenvalues of any of the spin-1/2 operators S _{x,} S _{y,} S _{z} are s = \pm (1)/(2) , any of their squares are given by s ^{2} 1 _{S} by the Cayley–Hamilton theorem.

Note also that the

2 L _{x} L _{y} term is invariant under the J -coupling evolution. That is that the term commutes with the Hamiltonian, and in this case, that can be manually confirmed by evaluating the commutator [ 2 L _{x} L _{y,} 2 L _{z} S _{z]} = 0 using the matrix representations of the spin operators.

CH₂

Now label the two hydrogen spins as

L, L' and the carbon spin by S . The J -coupling Hamiltonian is nowH = \pi J ( 2 L _{z} S _{z} + 2 L _{z'} S _{z)} which gives the following evolution\begin{align} ( 0 )&: L _{z} + L _{z'\\ (0)\to(1)&: L} _{z} + L _{z'} \xrightarrow{\left( (\pi)/(2) \right) _{xL} } -L_y + L_z' \xrightarrow -L_y - L_y'\\(1)\to(2)&:\ \xrightarrow 2 L_x S_z + 2 L_x' S_z\\(2)\to(3)&:\ 2 L_x S_z + 2 L_x' S_z\xrightarrow -2 L_x S_y - 2 L_x' S_y \xrightarrow -2 L_x S_y - 2 L_x' S_y \xrightarrow -2 L_x S_y - 2 L_x' S_y\\(3)\to(4)&:\ -2 L_x S_y - 2 L_x' S_y \xrightarrow 4 L_x L_z' S_x + 4 L_z L_x' S_x\\(4)\to(5)&:\ 4 L_x L_z' S_x + 4 L_z L_x' S_x \xrightarrow 4 L_x L_z' S_x + 4 L_z L_x' S_x \xrightarrow 4 L_z L_x' S_x \cos\theta - 4 L_z L_z' S_x \sin\theta + \text \xrightarrow -8 L_zL_z' S_x \cos\theta\sin\theta + \text\\(5)\to(6)&:\ -8 L_zL_z' S_x \cos\theta\sin\theta + \text \xrightarrow -4 L_z' S_y \cos\theta\sin\theta + \text \xrightarrow 2 S_x \cos\theta\sin\theta + \text\end

where 'others' denotes various terms that can safely be ignored because they will not evolve into observable transverse polarisation on the target spin

S . The required cyclic commutators for dealing with the J -coupling evolution are the following three sets (and their L \leftrightarrow L' versions if needed)\begin{align} \left.[L _{x} S _{y,} 2 L _{z'} S _{z]\right.} & = 2 L _{x} ⊗ L _{z'} ⊗ [S _{y,} S _{z]\\ &=} i 2 L _{x} L _{z'} S _{x \end{align}} \begin{align} \left.[ 2 L _{z'} S _{z,} 2 L _{x} L _{z'} S _{x]\right.} & = 4 L _{x} ⊗ {L ^{2} _{z'}} ⊗ [S _{z,} S _{x]\\ &=} 4 L _{x} ⊗ (1)/(4) 1 _{L'} ⊗ i S _{y\\ &=} i L _{x} S _{y, \end{align}} \begin{align} \left.[ 2 L _{z} L _{z'} S _{x,} 2 L _{z} S _{z]\right.} & = 4 L ^{2} _{z} ⊗ L _{z'} ⊗ [S _{x,} S _{z]\\ &=} 4 ⋅ (1)/(4) 1 _{L} ⊗ L _{z'} ⊗ - i S _{y\\ &=} i ⋅ - L _{z'} S _{y \end{align}} \begin{align} \left.[ 2 L _{z} L _{z,} - L _{z'} S _{y]\right.} & = - L _{z} ⊗ L _{z'} ⊗ [S _{z,} S _{y]\\ &=} i 2 L _{z} L _{z'} S _{x, \end{align}} \begin{align} \left.[ 2 L _{z'} S _{y,} 2 L _{z'} S _{z]\right.} & = 4 1 _{L} ⊗ L ^{2} _{z'} ⊗ [S _{y,} S _{z]\\ &=} 4 1 _{L} ⊗ (1)/(4) 1 _{L'} ⊗ i S _{x\\ &=} i S _{x \end{align}} \begin{align} \left.[ 2 L _{z'} S _{z,} S _{x]\right.} & = 2 1 _{L} ⊗ L _{z'} ⊗ [S _{z,} S _{x]\\ &=} i 2 L _{z'} ⊗ S _{y. \end{align}}

CH₃

A similar (but more lengthy) treatment gives the final observable term as

3 S _{x} \cos ^{2\theta\sin\theta} .

APT (Attached Proton Test)

Refer to for the notation used in this example.

APT is similar to DEPT in that it detects carbon multiplicity. However, it has additional degeneracies: it gives identical positive signals for C and CH₂, and identical negative signals for CH and CH₃. One variation on the basic pulse sequence is shown below.

The key observation is that since we can again ignore the refocussed chemical shift, the only relevant dynamics occur in the interval with no hydrogen decoupling, where we can consider solely the

J -coupling. By using an interval twice as long as in the DEPT case, we ensure that a density operator of L _{y} at the start of the interval just has its sign inverted following the coupling (since this corresponds to \alpha t = \pi in the general treatment, and \cos\pi = -1 , \sin\pi = 0 ). The Hamiltonians for the couplings to each of the n separate neighbouring hydrogen atoms commute, so the overall effect is to multiply by a factor ( -1 ) ^{n} . This motivates the alternating sign of the signal mentioned above.

References

• Web site: . Understanding NMR Spectroscopy . James Keeler . . 2012-08-05.
• Web site: University of Texas. A Pictorial Representation of Product Operator Formalism. David Donne. David Gorenstein . https://web.archive.org/web/20130514034653/http://www.biophysics.org/Portals/1/PDFs/Education/donne.pdf. 2012-08-05. 2013-05-14.
• 10.1002/cmr.1820030302. Echoes—how to generate, recognize, use or avoid them in MR-imaging sequences. Part I: Fundamental and not so fundamental properties of spin echoes. 1991. Hennig. Jürgen. Concepts in Magnetic Resonance. 3. 3. 125–143.
• 10.1021/ed070pA172.3. 2D NMR Density Matrix and Product Operator Treatment. Journal of Chemical Education. 1993. 70. 6. A172. 1993JChEd..70S.172.. Mateescu. Gheorghe D. Valeriu. Adrian.
• 10.1063/1.1628216. Exact effective Hamiltonian theory. II. Polynomial expansion of matrix functions and entangled unitary exponential operators. 2004. Siminovitch. David. Untidt. Thomas. Nielsen. Niels Chr.. The Journal of Chemical Physics. 120. 1. 51–66. 15267261. 2004JChPh.120...51S.
• Book: 10.1002/9780470034590.emrstm1310. Product Operator Formalism. EMag Res. 2013. Zhang. Yuning. Han. Fei. Jerschow. Alexej. 978-0470034590.