Density Matrix Description of NMR BCMB/CHEM 8190

Similar documents
Density Matrix Description of NMR BCMB/CHEM 8190

Let s treat the problem of the response of a system to an applied external force. Again,

( ) () we define the interaction representation by the unitary transformation () = ()

Linear Response Theory: The connection between QFT and experiments

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

Notes on the stability of dynamic systems and the use of Eigen Values.

Mechanics Physics 151

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

Mechanics Physics 151

Born Oppenheimer Approximation and Beyond

Lecture 6: Learning for Control (Generalised Linear Regression)

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

Mechanics Physics 151

Scattering at an Interface: Oblique Incidence

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

FI 3103 Quantum Physics

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

FTCS Solution to the Heat Equation

Lecture VI Regression

Solution in semi infinite diffusion couples (error function analysis)

4. THE DENSITY MATRIX

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Density Matrix Second Order Spectra BCMB/CHEM 8190

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Lecture 9: Dynamic Properties

( ) [ ] MAP Decision Rule

EP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

Chapter 6: AC Circuits

Bandlimited channel. Intersymbol interference (ISI) This non-ideal communication channel is also called dispersive channel

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

Advanced Machine Learning & Perception

On One Analytic Method of. Constructing Program Controls

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Linear Circuit Elements

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

Transcription: Messenger RNA, mrna, is produced and transported to Ribosomes

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

H = d d q 1 d d q N d d p 1 d d p N exp

PHYS 705: Classical Mechanics. Canonical Transformation

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Chapter Lagrangian Interpolation

Chapters 2 Kinematics. Position, Distance, Displacement

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

Fall 2010 Graduate Course on Dynamic Learning

Lecture 2 M/G/1 queues. M/G/1-queue

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

Math 128b Project. Jude Yuen

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

Time-interval analysis of β decay. V. Horvat and J. C. Hardy

Relative controllability of nonlinear systems with delays in control

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...

Fall 2009 Social Sciences 7418 University of Wisconsin-Madison. Problem Set 2 Answers (4) (6) di = D (10)

2.1 Constitutive Theory

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

Material Science Simulations using PWmat

CS286.2 Lecture 14: Quantum de Finetti Theorems II

Comparison of Differences between Power Means 1

Displacement, Velocity, and Acceleration. (WHERE and WHEN?)

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence

2/20/2013. EE 101 Midterm 2 Review

Volatility Interpolation

Polymerization Technology Laboratory Course

Pendulum Dynamics. = Ft tangential direction (2) radial direction (1)

C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology. September 1996

6.003: Signals and Systems. Relations among Fourier Representations

Normal Random Variable and its discriminant functions

QUANTUM TO CLASSICAL TRANSITION IN THE THEORY OF OPEN SYSTEMS

Lecture 11 SVM cont

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds

CHAPTER 10: LINEAR DISCRIMINATION

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

Motion in Two Dimensions

PHYS 1443 Section 001 Lecture #4

Robustness Experiments with Two Variance Components

6.003 Homework #8 Solutions

MCTDH Approach to Strong Field Dynamics

Advanced Macroeconomics II: Exchange economy

Graduate Macroeconomics 2 Problem set 5. - Solutions

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

Guest Lecturer Friday! Symbolic reasoning. Symbolic reasoning. Practice Problem day A. 2 B. 3 C. 4 D. 8 E. 16 Q25. Will Armentrout.

Li An-Ping. Beijing , P.R.China

Chapter Q1. We need to understand Classical wave first. 3/28/2004 H133 Spring

Predator - Prey Model Trajectories and the nonlinear conservation law

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

( ) = 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.

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

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

Chapter 2. First Order Scalar Equations

6.003: Signals and Systems

Transcription:

Densy Marx Descrpon of NMR BCMBCHEM 89

Operaors n Marx Noaon Alernae approach o second order specra: ask abou x magnezaon nsead of energes and ranson probables. If we say wh one bass se, properes vary only because of changes n he coeffcens of each bass se funcon x = h< Ix > - hs s how we calculae observables = c + c + c3 + c4 = j c j j AB example: =. +. +. +. =. +.6 +.8 +. 3 =. +.8 -.6 +. 4 =. +. +. +. < Ix > = j,k c j * c k < j Ix k > We need o calculae < j Ix k > only once f we say wh hs bass se hese can be pu n a n x n marx. Marx equvalen: < Ix > = c, c, * c Ix c

Specal Case: Paul Spn Marces Ix = ½ ½ Noe: < Ix > = ½ < > = < Ix > = ½ < > = ½ Iy = -½ ½ Iz = ½ ½ How do hey work? Try somehng we know: Ix > = ½ ½ ½ = = ½ = ½ ½ Operaors are a marx of numbers, Spn funcons a vecor of numbers

Larger Collecons of Spn ½ Nucle =??????? I Ax =???????? < I Ax > = < ½ > = ½? I Xx = Easer way: drec producs: E I Ax wh X marces

Densy Marces and Observables for an Ensemble of Spns Expecaon values peran o sngle spn properes; we observe ne behavor of an ensemble of spns. For a parcular sysem wo spn ½ nucle operaor represenaons are he same; all varaons are n bass se coeffcens. We need o average over hese coeffcens. < Ix > = c, c, * c Ix c Producs of averaged coeffcens can also be colleced and used n marx form. Ths s called a densy marx,. < Ix > = Tr c j c k * Ix = Tr Ix = j,k c j * c k Ix jk For a spn sysem: c*,c* Ix Ix c = c*, c* Ix c + Ix c Ix Ix c Ix c + Ix c = c*c Ix + c*c Ix + c*c Ix + c*c Ix For a spn sysem: < Ix > = Ix + Ix + Ix + Ix

Solvng for he Tme Dependence of Our observables are me dependen magnezaon precesses. All me dependence can also be pu n bass se coeffcens, or n coeffcen producs of he densy marx. Schrodnger s me dependen equaon H = - h dd Allows us o solve for dc j d, dd c j c k * = c j dc k *d + dc j d c k * Resul n erms of a densy marx s Lovlle von Neuman eq. dd = h { H - H } = h [, H] Implcaons: If we know a any me equlbrum a me and know he Hamlonan H, we can solve for a any me, and calculae any observable as he Tr O O s any operaor.

Tes: Does X Magnezaon Precess n B? Wha elemens of dcae x magnezaon? <I x > = Tr [] [I x ] - examne spn ½ case = Therefore, consder = Wha s H?, H = -hb I z = -hb d d B Tr

Consder evoluon of one elemen dd = B -½ ½ = Soluon: = exp- or: Same happens for = cos sn hence elemens correspondng o x magnezaon precess = =

Densy Marx a Equlbrum Have a roue o an observable: < Ix > = Tr Ix Have a roue o : dd = h [, H ] Need a place o sar: Dagonal elemens: consder c n c n * a probably c n c n * = exp-e n ktz Z - E n ktz; small E n nn = Z - nn ; Z # of saes Noe: ex uses for densy marx really devaon marx

Off-dagonal elemens a equlbrum mus sasfy me dependen Schrodnger Eq. H = h = E; n lm of H = H Soluon s: n = n exphe n Bu, he members of he ensemble do no have a common orgn n me. Hence add a random phase facor: nm = n exphe n + m = n c mn c mn = cosn,m, + snn,m, Densy marx elemen = c c j * = cos cos j + sn cos j +cos I sn j +sn I sn j f = j, Pcos + sn P If I j, cos cos j ec =, and all off dagonals =

Examples for a spn ½ nucleus H = -hb I Z, E = +- h4b eq exp B kt exp B kt eq B kt No me varaon as expeced: dd = h[, H ] = [ H] [H ] = produc of dagonal marces B kt

Wha abou M Z and M X? Cures law for suscepbly kt B kt B kt B Tr Tr M Z Z 4 4 4 ] ][ [ ] ];[ ][ [ Tr M I Tr M X X X X X Expec no X magnezaon a equlbrum

Wha elemens do gve M X? M X Tr We have argued ha hese are relaed o ranson probables We have shown ha hese precess n a magnec feld

Roaon operaors a more general descrpon of roaon dd = h [, H] General soluon: = exp-hhexphh = RR - If bass se elemens are egen funcons of H, R s dagonal and easy o evaluae Chemcal shf evoluon: R Z = exp-b cs I Z R Z exp exp ;

M X and M X Precess a sn sn cos cos exp exp exp exp exp exp exp exp exp exp M X M Y

Roaon operaors for an RF pulse along X axs n he roang frame Hamlonan n he roang frame: H = -B hi X Formal soluon a end of pulse me s gven as: = exp-hh exphh = R X R X - Canno smply nser I X n exponenal operaor o evaluae marx elemens; I X mxes spn saes Elemens of do no jus oscllae, bu conver one elemen no anoher for example - o + somehng we assocae wh z magnezaon o wha we assocae wh y magnezaon n a snb dependen way.

X pulse wh =, ; ; ; ;,,,, R R R R X X X X cos sn sn cos R X For one spn Smlar argumens lead o Y pulse operaons

The One Pulse FID = R Z R X, R - X, R Z - 9 X eq M p p Z

Evoluon of he FID Frs pon n he FID: = R Z d w p R Z - d w ; d w = dwell me; = exp exp exp exp W W W W d d d d = R Z dw R Z - dw; second pon n FID Calculang he FID: FID = Tr{[ ] [M X + M Y ]} Programs lke GAMMA, SPINEVOLUTION work hs way

Densy Marx Smulaon of nd Order Specra For wo spn bass se s:,,, 3... 3....... + 34 + 43 3 + 3 4 + 4,, are assocaed wh lnes n an AX specrum. : for frs order specrum Don have o have a : assocaon f s no an egen funcon of H. c I + c I may evolve coherenly as one lne e, 3, are mxed Calculaon of Mx sll works for a second order sysem