Lecture #8 Redfield theory of NMR relaxation

Similar documents
Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Linear Response Theory: The connection between QFT and experiments

ψ(t) = V x (0)V x (t)

Math 334 Fall 2011 Homework 11 Solutions

Two Coupled Oscillators / Normal Modes

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Chapter 7: Solving Trig Equations

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Chapter 2. First Order Scalar Equations

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

) were both constant and we brought them from under the integral.

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

Section 4.4 Logarithmic Properties

Lecture #9 Redfield theory of NMR relaxation

8. Basic RL and RC Circuits

Continuous Time Linear Time Invariant (LTI) Systems. Dr. Ali Hussein Muqaibel. Introduction

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Product Operators. Fundamentals of MR Alec Ricciuti 3 March 2011

Vehicle Arrival Models : Headway

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Class Meeting # 10: Introduction to the Wave Equation

10. LINEAR RESPONSE THEORYEquation Chapter 10 Section 1

HW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)

Math 333 Problem Set #2 Solution 14 February 2003

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

13.3 Term structure models

Numerical Dispersion

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

Waves are naturally found in plasmas and have to be dealt with. This includes instabilities, fluctuations, waveinduced

ψ ( t) = c n ( t) t n ( )ψ( ) t ku t,t 0 ψ I V kn

Book Corrections for Optimal Estimation of Dynamic Systems, 2 nd Edition

Non-equilibrium Green functions I

Notes for Lecture 17-18

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

CHAPTER 12 DIRECT CURRENT CIRCUITS

6.302 Feedback Systems Recitation 4: Complex Variables and the s-plane Prof. Joel L. Dawson

Final Spring 2007

Direct Current Circuits. February 19, 2014 Physics for Scientists & Engineers 2, Chapter 26 1

Section 4.4 Logarithmic Properties

k B 2 Radiofrequency pulses and hardware

Chapter 3 Boundary Value Problem

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Estimation of Poses with Particle Filters

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

4. The multiple use forestry maximum principle. This principle will be derived as in Heaps (1984) by considering perturbations

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

10. State Space Methods

System of Linear Differential Equations

Section 7.4 Modeling Changing Amplitude and Midline

Today in Physics 218: radiation reaction

Math 2214 Solution Test 1A Spring 2016

Relaxation. T1 Values. Longitudinal Relaxation. dm z dt. = " M z T 1. (1" e "t /T 1 ) M z. (t) = M 0

2001 November 15 Exam III Physics 191

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

KINEMATICS IN ONE DIMENSION

Chapter 6. Systems of First Order Linear Differential Equations

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Basic Circuit Elements Professor J R Lucas November 2001

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

6.003: Signal Processing

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

ψ ( t) = c n ( t ) n

Maxwell s Equations and Electromagnetic Waves

Linear Time-invariant systems, Convolution, and Cross-correlation

( ) is the stretch factor, and x the

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Sub Module 2.6. Measurement of transient temperature

LAPLACE TRANSFORM AND TRANSFER FUNCTION

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:

Circuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in

Y 0.4Y 0.45Y Y to a proper ARMA specification.

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

Some Basic Information about M-S-D Systems

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

From Complex Fourier Series to Fourier Transforms

A Bayesian Approach to Spectral Analysis

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Lecture Notes 2. The Hilbert Space Approach to Time Series

NMR Spectroscopy: Principles and Applications. Nagarajan Murali 1D - Methods Lecture 5

CHAPTER 2 Signals And Spectra

Chapter 8 The Complete Response of RL and RC Circuits

ln 2 1 ln y x c y C x

Transcription:

Lecure #8 Redfield heory of NMR relaxaion Topics The ineracion frame of reference Perurbaion heory The Maser Equaion Handous and Reading assignmens van de Ven, Chapers 6.2. Kowalewski, Chaper 4. Abragam Chaper VIII.C, pp 272-284, 1955. 1

The Maser Equaion When he Hamilonian can be wrien as he sum of a large saic componen plus a small ime-varying perurbaion Ĥ = Ĥ + ( Ĥ1 ) our goal is o find an equaion of he form: σ = iĥ σ Γ σ ( σ B ) Relaxaion Relaxaion superoperaor from which we can calculae any direc and cross relaxaion erms of ineres. In general, Γ differs across differen relaxaion mechanisms. Reminder: Bloch s d M! equaions are: = γ M! B z M xx + M y y ( M M z ) z d T 2 Roaions Roaions Relaxaion erms T 1 2

Redfield Theory Redfield heory, also known as Wangsness, Bloch, and Redfield (WBR) heory, is more general han ha derived by Solomon. Boh Solomon and Redfield heory rely on 2 nd -order perurbaion heory (can be a limiaion, bu usually no for liquids) Like Solomon s approach, Redfield heory is semi-classical, using he same Bolzmann correcion for hermal equilibrium values. However, raher han direcly dealing wih energy level populaions, he heory is derived in erms of he densiy operaor. This allows for a more general descripion of relaxaion and allows he derivaion of relaxaion raes from muliple mechanisms, e.g. dipolar coupling, CSA, scalar relaxaion of he 1 s and 2 nd kind, ec. 3

Liouville-von Neumann equaion We ve seen where we can express he spin Hamilonian as he sum of a large saic componen plus a small ime-dependen spin laice ineracion erm. Ĥ() = Ĥ + Ĥ1() perurbaion Saring wih he Liouville-von Neumann equaion for he densiy operaor dσ d = iĥ σ I is helpful o swich o a roaing frame of reference where he spin-laice ineracion erm is isolaed from he saic par of Ĥ. 4

The Ineracion Frame Le σ! = e iĥ σ H! ( ) = e iĥ Ĥ Roaion abou he Ĥ axis in Liouville space H1! = e iĥ Ĥ 1, hen d σ d! = d d eiĥ σ which simplifies o: = iĥ e iĥ σ + e iĥ d σ d, d σ d! = i H1! σ! (see homework) Hence in his frame of reference, known as he ineracion frame, where he ime dependence of he densiy operaor depends only on H 1!. 5

Redfield heory Sep 1. Le s sar by formally inegraing d σ d! = i H1! σ!. σ! = σ! = σ! = σ! i H1 " ( " ) σ! ( " )d! i H1 " ( " ) σ! d " H1 " ( " ) H1! (!! ) σ! (!! )d!! d! i H1! (! ) σ!!!!! d! H1!! H1!! σ d!! d! + i!!! This is saring o ge ugly, so we ll jus keep he firs 3 erms (o be jusified laer). 6

Redfield heory (con.) Sep 2 is o ake he ensemble average (echnically his is an ensemble of ensembles), and noe.. H1! (! ) σ! d " = because H1! This leads o. σ! σ!! = = H1! (! ) H1! (!! ) σ! H1! and we ve assumed H1! and σ! are uncorrelaed. d!! d! If, we can always incorporae he non-zero par ino Ĥ. Technically, from here on we should be using σ! jus going o go wih he simpler noaion of σ!. bu we re 7

Redfield heory (con.) Sep 3. Choose = Δ very small, such ha σ! Δσ " = σ! ( Δ) σ" Δ! = H1 " ( " ) H1! (!! ) σ! Sep 4. Define a new variable, τ =!!! d!! d! σ! This assumpion is he one ha allows us o drop hose higher order erms. In essence, σ! is assume o vary slowly in ime as compared o H1!. Δσ " = σ! ( Δ) σ" Δ! = H1 " ( " ) H1! (! τ ) σ! dτd # 8

Redfield heory (con.) Sep 5. Inroduce a correlaion superoperaor: G τ = H1! (! ) H1! (! τ ) Sep 6. Make some more assumpions. Δσ" Δ Δ = Ĝ τ σ" dτ Assume Δ sufficienly small ha Δ σ" Δ d σ" d, bu ha Δ is sufficienly large ha we can exend he inegraion o. dσ! d = Ĝ τ σ! dτ Can we acually find such a Δ?!? 9

Redfield heory assumpions Le s look a some numbers Consider Δ = 1-6 s. Typically NMR relaxaion imes are on he order of milliseconds o seconds, hence Δ σ" d σ" Δ d. Tissue waer correlaion imes due o molecular umbling are on he order of 1-9 s. Hence G Δ for Δ = 1-6 s. In general, Redfield heory is valid for relaxaion imes several orders of magniude longer ha he correlaions imes driving he relaxaion processes, e.g. liquids. 1

Redfield heory (con.) Sep 7. We need o find an explici expression for Ĥ 1 H1! = e iĥ Ĥ 1. The bes choice is o express eigenoperaors of H. Ĥ 1 as a linear combinaion of eigenvalue Le  q be an eigenoperaors of Ĥ : Ĥ  q = e q  q Examples. Le Ĥ = ω I Î z ω S Ŝ z hen Î + = Î x + iî y H Î + = Ĥ ( Î x + iî ) y = ω I iî ω y IÎ x = ω IÎ+ Î + Ŝ + H Î + Ŝ + = ω I +ω S (See Problem Se 1 Î+Ŝ+ for more examples) 11

Redfield heory (con.) Hence, le Ĥ1 = F q q  q Eigenoperaor Random funcions of ime (ypically dependen on molecular orienaion) We firs noe, ha for Ĥ 1 o represen a physical process, Ĥ 1 mus be Hermiian. However, he F q Âqs can be complex. Thus, for every F q Âq he sum mus also conain a erm F * q Âq *, which, we will denoe as F q  q. 12

Redfield heory (con.) Noing Ĥ  q = e q  q e iĥ  q = e ie q  q We have H1! = e iĥ Ĥ 1 q = F q  q e ie q We can now wrie he full expression for he correlaion superoperaor as G τ = H1! (! ) H1! (! τ ) = F p! F q (! τ )  p  q e ie p! p q e ie q! τ Using a secular approximaion, one can show ha only he erms for which e p = -e q need o be kep. The oher erms average ou as hey oscillae fas as compared o he relaxaion rae of σ!. 13

Thus Ĝ τ Redfield heory (con.) = F q " q F q ( " τ )  q  q e ie qτ Defining a se of correlaion funcions: G q ( τ ) = F q ( " )F q " τ yields d σ! d = G q ( τ )e ie qτ  q  q σ! dτ. q Sep 9. Define a se of specral densiy funcions: J ( q e ) q = G q τ e ie qτ dτ, wih G q ( τ ) = G q e τ τ c Correlaion ime characerisic of he perurbaion. 14

Redfield heory (con.) The specral densiy funcions are: J ( q e ) q = G q τ e ie qτ dτ, As defined, J q s, are complex, however, in pracice he real par is much larger han he imaginary componen *. J q ( ω) = G q ( τ )e iωτ dτ G q ( τ )cosω q τ = G q wih G q ( τ ) = G q e τ τ c A more formal reamen shows he imaginary componens basically cancel due o he combinaion of he erms: A q  q and  q  q hus we ll use.. Correlaion ime characerisic of he perurbaion. τ c 1+ω 2 τ c 2. * If you re curious abou hese imaginary erms, look up dynamic frequency shifs in he NMR lieraure. 15

The Relaxaion Superoperaor Puing i all ogeher, we can now define he relaxaion superoperaor as: Γ = J ( q e ) q  q  q q wih J ( q e ) q = F q " F q ( " τ )e ie qτ dτ Subsiuing: dσ! d = Γ σ! We can now add he Bolzmann correcion and finally swich back o he laboraory frame * (see homework) o yield: dσ d = iĥ σ Γ σ ( σ B ) The Maser Equaion of NMR * Noe he relaxaion superoperaor is he same is boh frames of reference, as befis relaxaion imes. 16

Nex Lecure: Redfield heory II 17