Jump to content

User:Thypingvinen/sandbox/Product operator formalism

From Wikipedia, the free encyclopedia

Motivation for sets of spin-1/2 particles

[edit]

Throughout this section, the reduced Planck constant for convenience.

The product operator formalism is usually applied to sets of spin-1/2 particles, since the fact that the individual operators satisfy , where 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.

The main idea of the formalism is to make it easier to follow the system density operator , which evolves under a Hamiltonian according to the Liouville-von Neumann equation as

For a time-independent Hamiltonian, the density operator inherits its solutions from the Schrödinger time-evolution operator as

Density operator-state duality

[edit]

Suppose a single spin-1/2 is in the state , which is an eigenstate of the z-spin operator , that is . Similarly . Making use of the expansion of a Hermitian operator in terms of projections onto its eigenkets with eigenvalues as , the associated density operator is

where is the identity operator. Similarly, the density operator for the state is

Since the spin operators are all traceless and the expectation value of an operator for a system with density operator is , the terms proportional to the unit operator 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 corresponds to a density operator , while the state corresponds to a density operator . In exactly the same manner, polarisation along the positive x-axis, that is a state , corresponds to a density operator . 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 , the terms in the density operator (ignoring the identity on its own) can be interpreted as representing

  • - longitudinal magnetisation
  • - in-phase transverse magnetisation, which is the observable quantity in NMR.
  • - anti-phase longitudinal magnetisation
  • - longitudinal two-spin order
  • - other coherences, which are more difficult to interpret, but may evolve into other terms

where eg is a shorthand for the Kronecker product , where is the identity operator on the spin, and similarly is a shorthand for .

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 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 , so this is sometimes more convenient.

Evolution of the density operator

[edit]

Consider operators that obey the cyclic commutation relations

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.

Introduce also the commutation superoperator of an operator (in our case, this is more formally related to the adjoint representation of the Lie algebra whose elements are ), which acts as

In particular, for the cyclic operators, we have

and consequently for integer

An identity for two operators is

which can be derived by putting where is a scalar parameter, differentiating both sides with respect to , and noting that both sides satisfy the same differential equation in that parameter, with the same initial condition at . In particular, for some scalar parameter , we have

(1)

where the final equality follows from recognising the Taylor series for sine and cosine. Now suppose that the density operator at time zero is , and it is allowed to freely evolve under the Hamiltonian where is some scalar. Using the results above, the density operator at some later time will be given by

(2)

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 are the generators of rotations about the axes. An example of such a set of generators is just the spin operators themselves.

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

.

With more specific reference to the radiofrequency pulses applied during NMR experiments, a hard pulse with tip angle around a direction is written as above the arrow and corresponds to taking as the rotation generator in Equation 1. 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.

An example: the 180°-refocussing pulse

[edit]

The Hamiltonian for a single spin evolving under a chemical shift of angular frequency is

which means that in an ensemble of many such spins with slightly different chemical shifts, there is a dephasing of the magnetisation in the - plane. Consider the pulse sequence

where is a time interval. Starting in an equilibrium state with all the polarisation along the -axis, the evolution of an individual spin in the ensemble is

Hence this sequence refocusses 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 taking the intial state to the final state as , where explicitly

Essentially we want to find the propagator in the form , 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 for unitary operator , operator and function , as well as the mathematical similarity of the the spin operators with the physical rotation generators, which allow us to write

Hence 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 entagled propogators via the Cayley-Hamilton theorem and eigendecompositions) are used.

An example: the amplitude of a Hahn echo in an inhomogenous magnetic field

[edit]

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

where again is a time interval. Liberally dropping irrelevant terms, the evolution for a single spin with offset up to just after the second pulse is

Now consider an ensemble of spins in a magnetic field that is sufficiently inhomogenous 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 term, in the sense that for an individual spin we have

where we used the identities and .

It is the spins in the new population that has been generated by the second pulse, namely the one with , that will lead to the formation of an echo after evolution for the next interval. Therefore, remembering to include the introudced 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

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 propogator as

and then evaluate a transition probability by considering the result of applying this to a state representing polarisation in the transverse plane.

A more complicated example: DEPT (Distortionless Enhancement by Polarisation Transfer)

[edit]

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, CH2 and CH3 groups. It does this by exploiting the heteronuclear carbon-hydrogen -coupling and varying the tip angle of the final pulse in the sequence. The basic pulse sequence is shown below.

The basic pulse sequence for the DEPT NMR experiment

Under the weak coupling assumption, the chemical shift terms commute with the -coupling term in the Hamiltonian. Hence we can ignore the refocussed chemical shift (see § An example: the 180°-refocussing pulse) in the two intervals containing -pulses, namely and , and additionally refrain from evaluating the chemical shift evolution in the last period . The pulse separation time is adjusted to the coupling strength (with associated Hamiltonian coefficient ) such that it satisfies

,

because then the first term in the evolved density operator in Equation 2 vanishes under the pure coupling evolution between the pulses.

CH

[edit]

Label the hydrogen spin as , and the carbon spin by . For illustrative purposes, we assume that the equilibrium state only has polarisation on the -spin (in reality, there will also be polarisation on the spin, with the relative populations determined by the thermal Boltzmann factors). The -coupling Hamiltonian is

which gives the following evolution

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

and consequently the next cyclic rotation

where we used the 'mixed-product identity' , which relates the matrix and Kronecker products for compatible dimensions of , and also the fact that since the two eigenvalues of any of the spin-1/2 operators are , any of their squares are given by by the Cayley–Hamilton theorem.

Note also that the term is invariant under the -coupling evolution. That is that the term commutes with the Hamiltonian, and in this case, that can be manually confirmed by evaluating the commutator using the matrix representations of the spin operators.

CH2

[edit]

Now label the two hydrogen spins as and the carbon spin by . The -coupling Hamiltonian is now

which gives the following evolution

where 'others' denotes various terms that can safely be ignored because they will not evolve into observable transverse polarisation on the target spin . The required cyclic commutators for dealing with the -coupling evolution are the following three sets (and their versions if needed)

CH3

[edit]

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

Bibliography

[edit]
  • Hennig, Jürgen (1991). "Echoes—how to generate, recognize, use or avoid them in MR-imaging sequences. Part I: Fundamental and not so fundamental properties of spin echoes". Concepts in Magnetic Resonance. 3 (3): 125–143. doi:10.1002/cmr.1820030302.
  • Mateescu, Gheorghe D; Valeriu, Adrian (1993). "2D NMR Density Matrix and Product Operator Treatment". Journal of Chemical Education. 70 (6): A172. Bibcode:1993JChEd..70S.172.. doi:10.1021/ed070pA172.3.
  • Siminovitch, David; Untidt, Thomas; Nielsen, Niels Chr. (2004). "Exact effective Hamiltonian theory. II. Polynomial expansion of matrix functions and entangled unitary exponential operators". The Journal of Chemical Physics. 120 (1): 51–66. Bibcode:2004JChPh.120...51S. doi:10.1063/1.1628216. PMID 15267261.
  • Zhang, Yuning; Han, Fei; Jerschow, Alexej (2013). "Product Operator Formalism". EMag Res. doi:10.1002/9780470034590.emrstm1310. ISBN 978-0470034590.