<<

Lecture #2: NMR Operators, Populations, and Coherences • Topics – description of NMR – The Density – NMR in Liouville Space – Populations and coherences • Handouts and Reading assignments – Kowalewski, Chapter 1 – van de Ven, appendices A-C

1 Why Quantum for NMR?

• Classical picture of NMR (i.e. Bloch equations) is valid for collections of non-interacting spins.

• Individual spins don’t have an associated T1 or T2. • NMR relaxation is described statistically and is driven by interactions among spins, e.g. – Nuclear dipolar and J couplings – Chemical exchange – -nuclear interactions: shielding, dipolar coupling, J coupling • Such interactions are best described using semi-classical or fully quantum mechanical descriptions.

2 Wave Functions • For the classical concept of a trajectory (succession in time of the state of a classical particle), we substitute the concept of the ! of a particle characterized by a , ψ ( r , t ) . ! • ψ ( r , t ) - contains all info possible to obtain about the particle € - interpreted as a amplitude of the particle’s presence € with the probability density given by:

! ! 2 3! d Ρ(r ,t) = Cψ(r ,t) d r , C constant.

! 1 ! 2 ! ∫ dΡ(r ,t) =1 = ∫ ψ(r ,t) d 3r << ∞ C € square-integrable! ! 2 ! - wave functions typically normalized, i.e. ∫ ψ(r ,t) d 3r =1 3 € €

€ QM Math: Kets and Bras • Important to be familiar with the concept of treating functions as elements in a . • Example: consider f(n) for 1 ≤ n ≤ N, n integer. " f % f 1 3 ! $ ' f2 fN f2 f1 f f = f(n) 4 $ ' ... $ " ' n f # N & ! • Each wave function ψ ( r ) is associated with an element of a linear vector space F and is denoted as ψ (called a “ket”). Dirac ! € notation ψ (r ) ⇔ ψ ∈ F € “ ” • Elements of the dual space of€ F are denoted as χ (called a bra ).

• A linear vector€ space with a defined metric, in this case the scalar product χ ψ , is called a . € 4 QM Math: Bases and Operators n • We can express ψ as a of other kets: ψ = ∑ci fi i=1 N

• If it follows that: ∑ci f i = 0 implies ci = 0 for all i , then the kets i=1 fi are called linearly independent. • In an N-dimensional vector space, N linearly independent kets € constitute a , and every element ψ in F can be written as n € ψ = ψ f ∑ i i where ψi is the coefficient of the ket f i . i=1 $" 1 if i = j f f • We typically deal with orthonormal€ basis: i j = # %$ 0 if i ≠ j € € • An operator, O ˆ , generates one ket from another: Oˆ ψ = ξ Example of a linear operator: Oˆ ( λ ψ + µ ξ ) = λOˆ ψ + µOˆ ξ scalars € € 5 € QM Math: Vectors and Matrices • When a ket, ψ , is expressed as a linear combination of basis kets, f ! i , then a corresponding column vector, ψ , can be constructed. # ψ1 & N % ( € ! ψ2 ! † ∗ ∗ ∗ ψ = ψ f ψ = % ( ψ ψ = ψ ,ψ ,…,ψ ∑ i i ( 1 2 n ) € i 1 % " ( € = % ( $ψ n ' • Similarly, operators have representations. ! ! € € € x Oˆ y where O is a n x n matrix with elements: ˆ = x = Oy Oij = f i O f j € z # x0 & # 0 & z ! % ( ! % 2 2 ( Remember: vector and 0 y x y ψ ψ = % 0 ( ψ = % 0 + 0 ( €matrix representations€ % ( % ( € depends on the basis set! € $ z0 ' $ z0 ' € y x 0 x " 0 y basis:{ x , y , z } { x " , y " , z }

x y "€ € 6 € € € € € € QM Math: Eigenkets and Eigenvalues • If Oˆ ψ = λ ψ

scalar then ψ is called an eigenket of O ˆ with eigenvalue λ. € • The of a matrix is defined at the sum of the diagonal elements: € € Tr(O) = ∑Oii i • Since the trace of a matrix is invariant under a change of basis, we can also talk about the trace of the corresponding operator. €

Tr Oˆ = λ ˆ ( ) ∑ i where λi are the eigenvalues of O i

7 € € QM Math: Liouville Space

• Operators defined on an n-dimensional Hilbert space, are themselves elements of an n2-dimensional vector space known as Liouville space. • The trace is the metric in Liouville space: ( Aˆ | Bˆ ) = Tr ( Aˆ † B ˆ ) , and in contrast to Hilbert space, the product of two Liouville Space elements is defined. a Liouville space Aˆ , Bˆ ∈ L Aˆ B ˆ = Cˆ ∈ L In general, AˆBˆ ≠ BˆAˆ. • Operators that work on elements of a Liouville space are called superoperators (denoted with a double hat). € ˆ ˆ ˆ € AO 1 = O2 , and we will deal with just one superoperator: the ˆ : Aˆ B ˆ ≡ Aˆ , Bˆ = Aˆ B ˆ − Bˆ Aˆ Do superoperators have [ ] associated eigenoperators?

€ What about supermatrices? 8 € Summary of Vector Spaces

Hilbert Space H NMR can be described {ψi },i = 1,..., n in Hilbert space ...

Liouville Space L …or in Liouville (operator algebra) space. ˆ 2 Oi ,i = 1,..., n { } Key Concept NMR is easiest to Superoperator Algebra S understand in ! ˆ $ 4 Liouville space! " Sˆ % ,i = 1,..., n # i &

See R. Ernst, G. Bodenhausen, and A. Wokaun, Principles of Nuclear Magnetic Resonance in One and Two Dimensions, Oxford University Press, 1990. 9 Postulates of QM • The state of a physical system is defined by a ket ψ ( t ) that belongs to a linear vector space known as a Hilbert Space. • Every measurable quantity A is described by an Hermitian operator A ˆ whose eigenkets form a basis in state space.

• The only possible result of the measurement of A is an eigenvalue ˆ of the A , the probability of obtaining eigenvalue an is 2 ( ) = un ψ where u is the associated eigenket, and the state P an n of the system immediately after the measurement is ψ = un .

• The of the ψ ( t ) is governed by the Schrödinger € € equation: ∂ ψ(t) = −iHˆ (t)ψ(t) ∂t where H ˆ ( t ) , known as the Hamiltonian, is the operator for the H(t) associated with the total energy of the system. € 10 , Angular Momentum, and Magnetic Moment • Spin, angular momentum, and magnetic moment operators are linearly related. ˆ ˆ µˆ p = γL p = γ!I p , p = {x,y,z} magnetic moment angular momentum spin with the following : Iˆ ,Iˆ iIˆ Iˆ ,Iˆ iIˆ Iˆ ,Iˆ iIˆ [ x y ] = z [ y z ] = x [ z x ] = y € 1 1 • For a spin ½ particle: Iˆ + = + + and Iˆ − = − − z 2 z 2 ˆ eigenkets of I z € € € • Matrix representation in { + , − } basis: € 1" 0 1% 1# 0 −i& 1# 1 0 & I x, = $ ' I y, = % ( I z, = % ( 2# 1 0& 2$ i 0 ' 2$ 0 −1'

The

11 € € € Isolated Spin in a Magnetic Field • Goal: Find the appropriate wavefunction ψ ( t ) that describes a system consisting of a nucleus (spin = ½) in a uniform magnetic field. • Procedure: € ∂ – Given: Schrödinger’s Equation: ψ(t) = −iHˆ (t)ψ(t) ∂t – Find Hˆ (t). – Solve for ψ . € – Compute quantities of interest: e.g. components of

magnetic moment µˆx , µˆy , and µˆz . € Previously showed these correspond to the

familiar quantities Mx, My, and Mz.

12 A Collection of Spins

• In a typical experiment, the number of nuclear spins, N, can be very big, e.g. 1024, and the complete quantum state is described by the wavefunction (or state vector): Rather unwieldly N for large N! ψ = ∑cn ψn n=1

• To make matters worse, we rarely, if ever, know ψ precisely, e.g. cn may be different for each spin. At best, we have only a statistical model for the state of the system.

€ 13 Density Operator for a Pure State • Consider a spin system in a pure state with normalized state vector

ψ(t) = ∑ci (t) ui where { u i } form an . i • We previously showed that we could define an operator:

σˆψ (t) = ψ(t) ψ(t) for which € € - Expectation of observable Aˆ : Aˆ (t) u ˆ (t) u u Aˆ u Tr ˆ (t)Aˆ = ∑ j σ ψ i i j = {σ ψ } i, j - Time evolution: ∂ ˆ ∂ ! # σ ψ (t) = " ψ(t) ψ(t) $ ∂t ∂t = −i[Hˆ , ψ(t) ψ(t) ] ˆ = −iHσˆψ 14 € Density Operator: Statistical Mixture

• Consider a system consisting of a statistical mixture of states ψ n ˆ with associated pn, and ai be an eigenvalue of A with associated eigenket ui . • While there is no “average state vector”, it turns out there is an “average operator” known as the density operator. ui ui a = p a p Tr ˆ Pˆ … in general: P( i ) ∑ nPn ( i )= ∑ n {σ n i } n n " € % ˆ ˆ ˆ ˆ = Tr#∑ pnσ nPi & = Tr{σ Pi } $ n ' ˆ ˆ where σ = ∑ p n σ n is, by definition, the density operator for the system. n ∂ ˆ • It is easy to show that: Aˆ = Tr σˆ Aˆ and σˆ = −iHˆσˆ { } ∂t Liouville-von Neumann ensemble average equation 15 System of independent spin= ½ nuclei ˆ • Most convenient basis set is the eigenkets of H 0 : { + , − } + − $ 2 * ' c c c + " % L - “longitudinal magnetization” & + + − ) L I σ = $ ' I - “transverse magnetization” & * 2 ) € −€ c c c # I L & “vector” in a 4D % + − − ( € Liouville space called € “ ” • σ for a one-spin system can be expressed as: space € " L1 % $ ' # 1 0& # 0 1& # 0 0& # 0 0& I $ 1 ' σ = L1% ( + I1% €( + I2% ( + L2% ( σˆ € $ 0 0' $ 0 0' $ 1 0' $ 0 1' $ I2 ' € ˆ ˆ ˆ ˆ $ ' In operator form: σˆ = L1T11 + I1T12 + I2T21 + L2T22 # L2 &

1 i 1 # 1 0& # 0 2& # 0 − 2& € # 2 0 & • or … σ = a1% ( + a2% ( + a3% ( + a4 % ( € 1 i 1 Iˆ $ 0 1' $ 2 0' $ 2 0 ' $ 0 − 2' z € € Iˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ y σˆ = a1E + a2I x + a3I y + a4 I z E ,I x,I y,I z { } A 3D subspace of σˆ “product operator” basis set coherence space Iˆ 16 € x

€ € Coupled Two-Spin System ˆ • Eigenkets of H0 : { ++ , + − , −+ , − − }

++ + − −+ − − Off-diagonal elements called “coherences” ++ ! P1 C1,2 C1,3 C1,4 $ + − # C P C C & σ 2,1 2 2,3 2,4 −+ # C3,1 C3,2 P3 C3,4 & − − # & " C4,1 C4,2 C4,3 P4 % Diagonal elements called “populations” € • As with the single-spin case, the product operators form a convenient orthonormal basis set (16 in total). 1 ˆ ˆ ˆ ˆ ˆ ˆ ˆ In-phase single quantum coherences 2 E , I x, I y, I z, S x, S y, S z ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ 2I xS z, 2I yS z, 2I zS x, 2I zS y Anti-phase single quantum coherences ˆ ˆ 2I zS z Longitudinal two-spin order € 2Iˆ Sˆ , 2Iˆ Sˆ , 2Iˆ Sˆ , 2Iˆ Sˆ Linear combinations of double and x x y y x y y x zero quantum coherences € 17 € € Phase Coherence “Coherent superposition of quantum states”

Arrows represents spin axis

18 Product Operators Pictorial Examples (note, we are now considering pairs of spins):

Sˆ Iˆ z I x S

€ Net tendency for S spins to be +z€ No net tendency for S spins in any direction No net tendency for I spins in any direction Net tendency for I spins to be +x

2Iˆ Sˆ z z ? ?

€ No net tendency for I or S spins to be ±z If I or S is ±z, increased probability paired spin is ±z

19 Observations about σ ˆ ˆ ˆ ˆ • Let σ = ∑ b n B n where B n are an orthonormal set of Hermitian basis operators,n i.e. Bˆ | Bˆ = Tr Bˆ †Bˆ = Tr Bˆ Bˆ = δ . ( n m ) ( n m ) ( n m ) nm Bˆ Tr ˆ Bˆ = b Tr Bˆ Bˆ n = (σ n ) ∑ m ( m n ) = bn m ˆ ˆ Coefficients of the expansion€ are the ensemble ∴σˆ = € Bn Bn € ∑ € averages of the expected values of the respective n operators! € € € ˆ 1 ˆ€ ˆ ˆ ˆ ˆ ˆ ˆ • Example: Bn ∈ { 2 E, Ix, Sx, Iy, Sy,…, 2IzSz }, the 16 two-spin product operators Signals observable with Rf coils: € M Tr ˆIˆ Iˆ ˆ ˆ ˆ x = γ! {σ x } = γ! x M y = γ!Tr{σ Iy } = γ! Iy Some signals not directly observable with Rf coils: ˆ ˆ ˆ ˆ ˆ M z = γ! Iz Cxz = γ! 2IxSz Czz = γ! 2IzSz “antiphase x” “longitudinal two-spin order”

20 Liouville-von Neumann Equation • Time evolution of σ ˆ ( t ) : ∂ ˆ σˆ = −i[Hˆ ,σˆ ] = −iHˆ σ ˆ ∂t • Hamiltonian€ Hˆ - operator corresponding to energy of the system

Iˆ € € z σˆ (0) ˆ ˆ σ (t) I y No longer the case that σ ˆ ( t ) just € undergoes simple rotations in € operator space. ˆ I x

€ 21 Thermal Equilibrium

• The spin expressed in the eigenkets of the ˆ Hamiltonian H 0 has the following form: Off-diagonal elements called $ P1 C1,2 ! C1,n ' “coherences” (= 0 at thermal equilibrium). σ = & C2,1 P2 " # ) & # " " Cn−1,n ) % Cn,1 ! Cn,n−1 Pn ( Diagonal elements called “populations” • Using the Boltzmann distribution one can verify: 1 Hˆ /kT ˆ σˆ = e−! 0 where Z = Tr e−!H0 /kT . € B Z ( ) remember Hˆ defined as Eˆ / ! • Using the high temperature approximation:

1 σˆ ≈ Eˆ − ! Hˆ (spin 1/2 particles) B 2 ( kT 0 )

22 A little foreshadowing… The • As always, finding the Hamiltonian is the key step. • Our ultimate strategy will be to write the Hamiltonian as the sum of a large static component plus a small time-varying perturbation. ˆ ˆ ˆ H = H0 + H1 (t)

Rad226a/BioE326a Rad226b/BioE326b • Then look for an equation of the form: Relaxation superoperator ∂ ˆ ˆ ˆ ˆ ˆ ˆ σ = −iH0σ − Γ(σ −σ B ) ∂t Rotations Relaxation This equation is known as the Master Equation* of NMR. ! Note the similarity to dM ! Mx xˆ + My yˆ M − M zˆ = γM × B zˆ − − ( z 0 ) Bloch’s equations: dt 0 T T 2 1 Rotations Relaxation terms

*I lied in Rad226a, as what I called the “Master Equation” was one that only applied to branch diagrams. 23

€ Next lecture: The Nuclear Spin Hamiltonian

24