, given by. , I y. and I z. , are self adjoint, meaning that the adjoint of the operator is equal to the operator. This follows as A.

Size: px
Start display at page:

Download ", given by. , I y. and I z. , are self adjoint, meaning that the adjoint of the operator is equal to the operator. This follows as A."

Transcription

1 Further relaxation 6. ntrodution As resented so far the theory is aable of rediting the rate of transitions between energy levels i.e. it is onerned with oulations. The theory is thus erfetly aetable for rediting the rate onstants for relaxation of longitudinal magnetization ("T ") but is inaable of dealing with transverse magnetization ("T "). To do this it is neessary to onsider the evolution of oherenes under the random hamiltonian. A onvenient way of roeeding is to look at the relaxation behaviour of oerators. This fits in well with the Solomon equations an aroah already used for disussing the dynamis of z-magnetization and oulations. The oerator method is less well suited to disussing the relaxation behaviour of individual lines but it will be seen that with some modifiations is an be adated to this use. This Setion starts with a brief revision of the roerties of oerators. 6. Proerties of oerators 6.. Matrix reresentations Oerators were introdued in Setion. as the way in whih observable quantities are reresented in quantum mehanis. t is sometimes onvenient to think of oerators as matries and the matrix reresentation of an oerator an be formed in a artiular basis set of wavefuntions (Setions.. and.7.). The ijth element (meaning the element in row i and olumn j) of the matrix reresentation of an oerator A is A ij given by Aij = φ * i A φ j d τ φ i A φ j where φ i is a basis funtion. The same oerator will have different matrix reresentations in different sets of basis funtions. f a matrix has the roerty that A ji = A * the matrix is said to be hermetian. Oerators whose matrix reresentations are hermetian are alled hermetian oerators. Examles of suh oerators are the familiar angular momentum oerators x y and z. The adjoint of an oerator A A is defined in the following way * ( A ) = A ij ji Hermetian oerators suh as x are self adjoint meaning that the adjoint of the oerator is equal to the oerator. This follows as A * ji = Aij for a hermetian oerator and so A = A. ij ij The raising and lowering oerators + and are not hermetian a fat that an readily be areiated from the matrix reresentations: ij 6

2 = + i = + x y x iy 0 = = These oerators are the adjoints of one another: + = = + Again this is readily areiated from the matrix reresentations. The trae of an oerator is equal to the sum of its diagonal elements { }= Tr A 6.. Basis sets of oerators The density oerator σ an be exanded in a basis of oerators B i : bb bb bb i A ii σ = + + K = where the b i are numerial oeffiients. A basis set of wavefuntions an be desribed as being orthonormal (Setion..) whih means that they have the roerties φi φj = 0 if i j; and = if i = j For oerators orthogonality is defined using the trae. Two oerators B i and B j are orthogonal if Tr{ BB i j }= 0 To omute the trae it would be neessary to find the matrix reresentations of B i and B j multily these together to give another matrix and then add u the diagonal elements of the resulting matrix. The oerators x y and z are all orthogonal to one another. Normalization also involves omutation of the trae of the adjoint of an oerator with itself Tr{ BB i i }= β i Sometimes a basis set of oerators are hosen so that β i = and sometimes the oerators are hosen so that β i is the same for eah oerator but not neessarily =. For examle: { }= { } { }=. Tr Tr = Tr x x y y z z The set of rodut oerators for two sins ( x x z et.) all have β i = ; the fator of in the oerator roduts with two sin oerators is needed to kee the β i the same for all the oerators. 6.. Commutators The ommutator of two oerators is written [AB] and defined as follows [ A B]= AB BA f A and B are numbers or simle funtions then the order in whih they aear does not matter. For examle = and x y = y x. Suh simle funtions are said to ommute; the ommutator between them is zero i i i 6

3 [ xy ]= xy yx= 0 However not all oerators are simle funtions and so they do not neessarily ommute with one another. n artiular the angular momentum oerators obey the following ommutation relations x y and z. [ ]= [ ]= [ ]= i i i x y z y z x z x y The seond ommutator an be found from the first by yli ermutation of the indies: x y y z z x; likewise the third an be found from the seond by further yli ermutation. Using these relationshis and the definitions of + and the following an be derived [ ]= [ ]= [ ]= + z + z + z Oerators for different sins always ommute with one another and an oerator always ommutes with itself. For examle [ x z]= [ x z z x z z]= Commutator relations The following roerties are often useful [ AB ]= [ BA ] [ PA B]= P[ A B] if P ommutes with A and B A+ B C A C B C [ ]= [ ]+ [ ] 6... Oerator roduts The ommutator relations given in Setion 6.. for x y and z aly for any sin. For sin only the following relationshis also aly = i = i and for yli ermutations of x y and z x y z y x z = = = E where Eis the unit oerator or matrix x y z There are also similar seial relationshis for the raising and lowering oerators. = ( ) + = = 0 + = + z + Sine [ x y]= iz and x y= iz it follows that = [ ] = [ ] x y x y y x y x These relationshis are also valid for yli ermutations of x y and z. f these relationshis aly (i.e. for sin half) then it follows that [ AP BQ]= 0 if A and B both ommute with P and Q An examle of the latter is when A and B are oerators of sin and P and Q are oerators of sin [ xy yz]= 0 Another seial ase whih will be of use is z 6

4 AP AQ A P Q P Q A or rovided A ommutes with P and Q. 6 [ ]= [ ]= [ ] = x y z 6. Oerator equations The time evolution of oherenes is onveniently handled using the density matrix introdued in Setion.7. Further it is useful to make an oerator exansion of this density matrix as this fits in well with the rodut oerator aroah (Setion ) used to desribed many multile ulse exeriments. n addition this oerator aroah is losely related to the Solomon equations introdued in Setion 5... The starting oint is the equation of motion of the density matrix (or oerator) dσ () t = ih [ σ () t] where the square braket indiates the ommutator [ A B]= AB BA As before the hamiltonian onsists of a stati art H 0 and a random time deendent art H (t) H = H0 + H() t A transformed density oerator σ T and transformed time deendent hamiltonian H T () t are defined as T T σ ()= t ex( iht) σ() t ex( iht) H ()= t ex( ihth ) () t ex( iht) t an be shown that the equation of motion of σ T is () = () () T dσ t T T ih [ t σ t] [] This transformed reresentation is alled the interation reresentation. t is useful as σ T just evolves under the ation of the random hamiltonian; the evolution due to the large stati hamiltonian has been fatored out. t will be seen that the influene of H 0 aears as a simle hase evolution of H T (). t Equation [] an be solved to seond order in the erturbation reresented by the random hamiltonian to give () = () [ ] T dσ t T T T H t [ H ( t τ ) σ ( t ) ] dτ [] 0 There are various aroximations involved in this derivation all of whih rely on the fat that the random hamiltonian reresents a small erturbation and that the theory is to be used for times longer than the orrelation time τ. Equation [] is the master equation for the motion of the density oerator. However this equation has one defet: it redits that at long times the density oerator goes to zero. This is inorret at long times the density oerator must go to its equilibrium value whih is not zero but orresonds to the equilibrium oulation distribution. The reason for this defet is

5 idential to that disussed in Setion 5.. it omes from the failure to treat the lattie as a quantum objet. The solution to this roblem is either to treat that lattie roerly (whih will be rather omlex) or simly to relae σ T by (σ T σ eq ) where σ eq is the equilibrium density oerator. n all that follows it will be assumed that this relaement is made at the end of the alulation. 6.. Oerator exansion To develo Eq. [] into a useful form H (t) is exanded aording to q H t F t A q ( q) = () where the F (q) (t) are satial funtions and the A (q) are sin oerators. This exansion allows for the ossibility that there are several different terms with different satial fators in the random hamiltonian. Sine ollisions ause the moleule to randomly hange its orientation the satial funtions beome random funtions of time. The oerators A (q) are hosen so that they evolve under the stati hamiltonian aording to = 0 0 q q q ex ih t A ex ih t ex iω t A [] where the frequeny ω (q) is harateristi of the oerator A (q). For examle in the ase of a single sin for whih H 0 = ω 0 z the oerator + evolves in the way given by Eq. []: ex( ih t) ex( ih t)= ex( iω t) with ω (q) = ω 0. Likewise and z evolve at frequenies +ω 0 and 0 resetively. n more omlex ases it may be that A (q) is a sum of oerators all of whih have the same satial funtion F (q) (t) but all of whih do not evolve at the same frequeny under the transformation reresented by Eq. []. For examle in the diolar interation there is a set of oerators S + z z S + S whih all have the same satial deendene but whih under the stati hamiltonian H 0 = ω z + ω S S z evolve at frequenies 0 ω +ω S and ω ω S resetively. To oe with situation suh as these Eq. [] is modified to ( q) ( ex ih t A ex ih t ex i q ) ( q) ( 0 ) 0 ω t A [] = where the sum over allows for the ossibility that the A (q) is a sum of ( q oerators A ) ( q eah with its own assoiated frequeny ω ). The hamiltonian must be hermetian so it follows that for eah oerator with a ositive index q there must be one with a negative index with the two oerators being related by ( q) ( q) [ A ] = A [5] where the dagger ( ) indiates the adjoint. For examle 6 5

6 6 6 [ ] = [ z] = + Likewise for eah F (q) there is an assoiated F ( q) related by F ( q) = F (q)*. t ( q) ( q) follows from these definitions that ω = ω. Substituting the exansion of Eq. [] into the master equation (Eq. []) gives after onsiderable maniulation T dσ t where J q () = () q' q ' q' ( q) ( q' ) T ( A A σ t Jq q ω q ' ) ( [ [ ]] ( ) ( i( ω q ) ( + ω q ' ) ' ' ' ' ) t) ex ω is the Fourier transform of the orrelation funtion G q = ( ) J ω G τ iωτ dτ q q' q q' ex and where the orrelation funtion is defined as q q' = () z q' [6] ( τ ) q q' G τ F t F t τ [7] 6.. Simlifiation by negleting ross orrelation and non-seular terms The relaxation behaviour redited by Eq. [6] deends on the orrelation funtion between all ossible interations i.e. all ossible airs qq'. However as has been seen above it is often the ase or an be assumed to be the ase that there is no orrelation between the flutuations of different terms in the random hamiltonian; in other words ross orrelation an be ignored. f ross orrelation is indeed ignored then the only non-zero orrelation funtions are ones with the same index q for both funtions i.e. () ( τ ) q q * F t F t By definition F ( q) = F (q)* this non zero orrelation funtion an be written () ( τ ) ( ) q q F t F t whih from Eq. [7] is G q q' with q' = q i.e. G qq n Eq. [6] there is a omlex exonential whih auses a hase osillation q q' ( ( + ω ' ) t) ex i ω q q f ross orrelation is ignored q' = q and as ω = ω this term beomes q q q q ex i ω ( ω ) ( ( + ' ) t)= ex( i( ω ω ' ) t) ( q) ( q) f ω and ω ' differ signifiantly then this terms auses a raid hase osillation so the term will not ontribute to relaxation simly as it is onstantly hanging sign. The terms whih will ontribute are those with ω( q) ( q) = ω ' for whih there is no hase osillation. This ondition is only met if = '. Terms whih satisfy this restrition are termed seular ontributions (as oosed to the remainder whih are non-seular). Thus retaining only the seular terms and ignoring ross orrelation simlifies Eq. [6] to

7 T dσ t q q T [ A [ A ( σ t ]] ) Jqq( ω ( q) ' ) q () = () [8] 6.. Relaxation of individual oerators As was desribed in Setion 6.. it is ommon to exand the density oerator is exanded in terms of a basis of oerators B i σ()= t b() t B i where the oeffiients b i arry the time deendene. Evolution under ulses and delays auses the oerators to transform into one another. t would be useful to have a similar way of handling the effet of relaxation on the oerators. n suh an aroah eah oerator would relax in a harateristi way erhas being transferred to other oerators. The theory would redit these rate onstants for self relaxation and for transfer to other oerators. Suh a way of treating the motion of oerators under the influene of relaxation has already been introdued as the Solomon equations. For examle the equation = 0 d z z 0 0 R ( z z ) σ S( Sz Sz ) says that z relaxes with a rate onstant R and that S z is transferred to z with a rate onstant σ S. The aim is to use the master equation Eq. [8] to rodue similar equations for any oerator B i. Suose that the density oerator just ontains a single oerator B with oeffiient b (t). To use Eq. [8] the density oerator has to be transformed to the interation reresentation σ T ()= t ex( iht 0 ) σ() t ex( iht 0 ) = ex( ih0t) b() t Bex( ih0t) = b() t B ex( iω t) where it has been assumed that under this transformation the oerator B aquires a hase modulation at a harateristi frequeny ω. The exonential term an be merged into the oeffiient by defining so that ()= () b' t b t ex iω t ()= () σ T t b' t B The right-hand side of Eq. [8] beomes i [ [ ]] ( ( q A A b t B J ) ' qq ω ) ( q ) q '() q The effet of the double ommutator will in general be to transform the oerator B into other oerators; there will be some oeffiient for this transfer deending on the details of the oerators and the setral densities. n general the result an be written i 6 7

8 () + () + () + K b ' tb b' tb b ' tb where j is the oeffiient for the transfer of oerator to oerator j. So Eq. [8] beomes d ( b' () t B)= b' () t B + b' () t B + b' () t B + K = b' () t B j j j However in general the density oerator will not start out as just onsisting of a single oerator B but will be a sum of many oerators σ T ()= t b' () tb + b' () tb + b' () tbk = b' () t B i i i Eah of these oerators an be transformed into any other oerator by the ation of the double ommutators and if as above the oeffiient for transfer of oerator B j into oerator B i is ji the result is d bi' () t Bi b' t B b' t B b' t B b' tb b' tb b' tb + b ' () tb + b' () tb + b' tb = b' () tb i = () + () + () + + () + () + () + i j ji j i () + K K K This begins to look like a Solomon equation but is rather omlex as it has a sum of derivatives on the left. To ik out just one of these say that of the oerator B k both left and right hand sides are multilied by the oerator B k and the trae taken d j bj' () t Tr{ Bk Bj} jibj' t Tr Bk Bi = () { } The oerators are assumed to be orthogonal and to have a onstant normalization fator β: Tr { B Bq}= 0 if q = β if = q The β anel and so the differential equation beomes d ( bk' () t )= jkbj' () t [9] j This is a Solomon-tye equation. What it says is that if the rate of hange of the "amount" of oerator B k b k (t) deends on the amounts of all other oerators B j b j (t) resent. The rate onstant for transfer of oerator B j into oerator B k is given jk. The rate onstant jk is simly found by omuting i j [ ] q q q jk = Tr ( Bk ) A [ A ' Bj ] Jqq ω β ( ) [0] q 6 8

9 n words what this means that starting with a single oerator B j it is seen how this oerator is transformed into other oerators by the double ommutators; from the resulting set of oerators just the omonent of B k is iked out. A knowledge of the jk rate onstants gives a omlete set of Solomon equations and so a omlete desrition of the relaxation behaviour of the set of oerators. n ratie therefore the theory is used to alulate these jk rather than the relaxation of an arbitrary density oerator Time deendene Equation [9] is written in the interation reresentation with the oeffiients b i '(t) having a hase osillation due to the evolution of the oerators under the stati hamiltonian. ()= () b' t b t ex iω t [] i i i The aim is now to searate the evolution due to relaxation from that due to evolution under the stati hamiltonian. Starting from Eq. [9] the oeffiients b i '(t) are relaed by those from Eq. [] d ( bk' () t )= jkbj' () t j d ( bk() t ex( iωkt) )= jkbj() t ex( iω jt) j dbk () t ex( iωkt)+ bk() t iωkex( iωkt)= jkbj() t ex( iω jt) j dbk () t = bk () tiωk + jkbj() tex ( i( ω j ωk ) t) [] j Equation [] an be interreted in the following way. The first term on the right desribes the hase evolution of the oerator B k due to the stati hamiltonian; it has nothing to do with relaxation. The sum of terms on the right gives the rate of relaxation-indued transfer from oerator B j to B k with rate onstant jk. However eah term on the right inludes a hase osillation at (ω j ω k ); suh terms will not be effetive at transferring B j to B k unless (ω j ω k ) = 0 so that this hase osillation goes away. Unless this is the ase the sign of the relaxation terms will hange raidly (muh more raidly than the rather slow rate of relaxation) and so the net effet will be zero. For examle suose that the two oerators involved are + and S + ; the harateristi frequenies are therefore ω and ω S (the two Larmor frequenies) and as these are not the same any relaxation indued transfer between the two oerators whih in rinile might take lae will atually have no net effet. n ontrast the two oerators z and S z whih both have harateristi frequenies of zero an have relaxation indued transfer between them (just as has been seen in the Solomon equations). So Eq. [] an be written 6 9

10 dbk t bk t iω k j' kbj' t [] j' where the sum over j' is taken to inlude only those terms for whih ω j' = ω k. Conentrating on just the relaxation arts of the equation it is ossible to write Eq. [] in terms of the oerators rather than their oeffiients as was done in the ase of the Solomon equations dbk = jk ' Bj' 6 0 () = () + () j' As before it is imortant to remember that this equation is just a short hand it is not the oerators that are hanging with time but their ontributions to the density oerator. 6.. Sueroerators Often it is aetable to think of relaxation and evolution under the stati hamiltonian as being searate even though in reality they are learly taking lae at the same time. However there are oasions on whih it is essential to onsider the two roesses ating at the same time. Essentially this is what a omlete set of equations like Eq. [] (for all oerators B k ) will do. A onvenient way of exressing suh equations is to use the idea of sueroerators. The starting oint is the equation of motion of the density oerator: dσ () t = ih [ σ () t] For suh an equation it is usual to think of σ and H as being matries exressed in some basis of wavefuntions. However if σ is exressed as a linear ombination of oerators as in the rodut oerator formalism it is then natural to think of σ as a (olumn) vetor of the oeffiients the individual oerators. The equation of motion is then written dσ () t ˆ = ihσ () t [] where Ĥ is the hamiltonian sueroerator; this an be thought of as a matrix. Formally the sueroerator is defined as a ommutator: ˆ Hσ = [ H σ]= Hσ σh f H is time indeendent the solution to Eq. [] is σ()= t ˆ ex iht σ ( 0) Relaxation an be added to Eq. [] by introduing a relaxation sueroerator ˆΓ dσ() t ih ˆ σ t ˆΓ = () { σ () t σ } eq where ˆΓ is defined via

11 [ [ ' ]] q qq ( ) q ˆ T T Γσ t A ( q ) A q σ t J ω ()= () n the ase that σ is written as a vetor of oerators ˆΓ will be a matrix with the elements jk given by Eq. [0]. 6. Examle alulations 6.. One sin with random field The simlest ase to onsider is a single sin exeriening a random field in the x y and z diretions given by the hamiltonian H t ω t ω t ω t = () + () + () x x y y z z t will be onvenient to re-exress the x and y oerators in terms of the raising and lowering oerators to give + + H t ω t ω t ω t ()= () [ + ]+ () [ ]+ () x y i z z [ x y ]+ [ x()+ y() ]+ z() z = () () + ω t iω t ω t iω t ω t From this three searate terms are identified q A (q) F (q) (t) ω (q) 0 z ω z () t 0 + ωx() t iωy() t ω ωx()+ t iωy() t +ω where the stati hamiltonian H 0 has been taken as ω z so that ex( ih0t) ± ex( ih0t)= ex( miω t) ± ex( ih0t) zex( ih0t)= z For eah value of q there is just one term so the index in Eq. [8] is not needed. The orrelation funtion G (τ) is () ( ) G( τ)= F () t F ( t + τ) = ω () t iω () t ω ( t τ)+ iω ( t + τ) [ ][ + ] x y x y [ x x+ x() y( + ) y() x( + )+ y() y( + τ )] = ω () t ω t + τ iω t ω t τ iω t ω t τ ω t ω t f ross orrelation between the x and y omonents is ignored then the ensemble averages of the seond and third terms is zero. The first and seond terms will be written So that ω () t ω ( t + τ)= ω g( τ) ω ( t) ω ( t + τ)= ω g( τ) x x x y y y The orresonding setral density is G J [ x y] ( τ)= ω + ω g τ [ x y] ( ω)= ω + ω j ω 6

12 The orrelation funtion G (τ) is the same as G (τ). G 00 (τ) is similarly given by ( 0) ( 0) G00( τ)= F () t F ( t + τ) = ωz() t ωz( t + τ) = ωz g( τ) and so J00( ω)= ωz j( ω). The aroah will be to start with Eq. [0] [ ] q q q jk = Tr ( Bk ) A [ A ' Bj ] Jqq ω β ( ) [0] q and to ut eah oerator B j in turn into this relationshi. To start with just the art in the round brakets will be evaluated Motion of z The term to alulate is 6 [ ] ( ( q ) qq ω ) ( q) q [ ' z] A A J q where q = 0 + ; the index is not used. The following terms are found = [ [ ]] = [ ] = = [ [ ]] [ ] = [ ] ω = [ ] = q 0 : z z z J00 0 z 0 J q : J ( ω )= J ( ω )= J ω z z + [ ] ( ) J ω zj ω + q : z J The double ommutators generate no new oerators so all that haens is that z relaxes on its own. t is not therefore neessary to go through the formal alulation of the trae. Rather the oeffiient 00 an be iked out as: { + } 00 = J ω J ω = J( ω) J( ω) = [ ωx + ωy] j( ω) [ ωx + ωy] j( ω) = [ ωx + ωy] j( ω) where is has been assumed that the setral densities at ±ω 0 are the same. The Solomon equation is therefore dz = 00z = R z z where R z the longitudinal relaxation rate onstant is Rz = [ ω x + ω y] j( ω ) As exeted the rate of longitudinal relaxation deends on the size of transverse fields and the setral density at the Larmor frequeny.

13 Motion of + The term to alulate this time is ( q) [ ] qq ( ω ) ( q) [ ( q ) ' + ] A A J q where q = 0 + ; the index is not used. The following terms are found [ z z + ] 00=[ z + ] 00= 00 + q 0 : J 0 J 0 J 0 = [ ] = [ ] = [ ] [ ] = [ ]= q : J ω [ + + ] ( )= [ + z] = ( ω ) q : J ω J ω J + The double ommutators generate no new oerators so as before all that haens is that + relaxes on its own. The oeffiient that is the rate onstant for this roess an be extrated by insetion as [ 00 ω ] J 0 J = + = ω j( 0) ω + ω j ω [ ] z x y The Solomon-tye equation is therefore d+ = + = R t + where R t the transverse relaxation rate onstant is [ ] Rt = ωz j( 0)+ ωx + ωy j ω 6... nterretation t is interesting to omare the longitudinal and transverse relaxation rate onstants Rz = [ ω x + ω y] j( ω ) 0 [ ] = + Rt = ω j( 0)+ ω + ω j ω0 ω j 0 R z x y z z Longitudinal relaxation only deends on the setral density at the Larmor frequeny but transverse relaxation deends on the setral density both at the Larmor frequeny and at zero frequeny. The two ontributions to transverse relaxation are quite distint. Longitudinal random fields give a ontribution whih deends on the setral density at zero frequeny; this is alled the seular ontribution as no energy hange is involved. Transverse random fields give a ontribution whih deends on the setral density at the Larmor frequeny; this is alled the non-seular ontribution as hanges of energy are involved as the sins fli. Furthermore longitudinal relaxation is only brought about by transverse magneti fields but both transverse and longitudinal magneti fields give rise to transverse relaxation. The longitudinal fields are assoiated with the setral density at zero frequeny. The seular art of transverse relaxation 6

14 Reall that transverse magnetization is due to the resene of a oherene in the system; this oherene has a hase osillation at the Larmor frequeny and this in turn is what gives rise to the measured reession of transverse magnetization at the Larmor frequeny. Eah sin exerienes the alied magneti field along the z-diretion lus any loal random field ω z (t). As a result the frequeny of the hase osillation of the ontribution to the oherene from eah sin has a sread of values aross the samle. As time roeeds therefore these individual ontributions get out of hase with one another and so the oherene (and hene the transverse magnetization) deays. For suh a mehanism to be effetive at dehasing oherene it is not neessary for there to be any time variation in the fields along the z- diretion all that is required is that there is a distribution of fields aross the samle. However the dehasing aused by a omletely stati distribution of fields is not lassed as relaxation sine the dehasing ould be reversed by the aliation of a 80 refousing ulse. ndeed this kind of dehasing is exatly that whih ours when a field gradient ulse is alied. Relaxation is a dissiative irreversible roess whose effets annot be undone with ulses. So for a distribution of fields in the z-diretion to be effetive at ausing transverse relaxation these fields must be time deendent so that their effet annot be reversed by a refousing ulse. The time deendene has a different role in transverse and longitudinal relaxation. n the latter the time deendene is needed to ause transitions in the former it is needed to inhibit refousing. The seular art of transverse relaxation deends on j(0) whih = τ. So as the orrelation time gets shorter the rate of transverse relaxation dereases. To understand why this is it is useful to imagine first that the sins are frozen and not moving. The distribution of fields aross the samle will give rise to a range of different Larmor frequenies and as a result the line observed in the setrum will be broad. Now suose that the sins start to move; they will jum from osition to osition exeriening a different field eah time (just like hemial exhange but between very many sites). As the juming rate beomes omarable with the original linewih the line will start to narrow (again just like exhange narrowing). What is haening is that there are now so many jums that all the sins are beginning to see the same average frequeny. As the jums beome muh faster than the original linewih the line is narrowed drastially; however there is still a residual wih whih deends on the rate of juming and the original linewih for the stati arrangement. For tyial NMR samles the linewih of the frozen samle would be rather large (tens of khz or more). However tyial orrelation times for liquids are in the s to ns range; these times are an indiation of the time between jums. With suh fast motion omared to the frozen linewih the line observed in the liquid is very muh narrower than in the frozen samle. Nevertheless the extent of narrowing deends on the rate of 6

15 juming and this is why the relaxation rate onstant is roortional to τ. The smaller τ the faster the juming and the narrower the line. The non-seular art of transverse relaxation Longitudinal relaxation is assoiated with hanges in oulations and hene transitions between sin states. Suh roesses involve a transfer of energy between the sins and the lattie. t is also lear that a transition from one sin state to another will ause an interrution in the hase of the ontribution of a artiular sin to the overall oherenes. Thus suh transitions also lead to transverse relaxation. The non-seular ontribution to the transverse rate onstant is R z. The half arises beause one sin fliing auses the oulation differene and hene the z-magnetization to hange by two units. Sin flis are thus twie as effetive at ausing longitudinal relaxation as they are at ausing transverse relaxation. 6.. Variation with orrelation times For simliity it will be assumed that n whih ase ω = ω = ω = ω R x y z Rz ωr j ω Rt ωr j 0 ωr j ω The simlest model is to assume that the orrelation funtion is exonential g(t) = ex( τ/τ ); this gives a lorentzian setral density = = + τ j( ω)= + ωτ The relaxation rate onstants are therefore τω R τ Rz = Rt = ωr τ + + ωτ + ωτ Several imortant oints an be drawn from these relationshis. The first is that the rate of transverse relaxation is always greater than or equal to that of longitudinal relaxation. This follows as the differene is learly always ositive. R Rt R z = τωωτ [5] + ωτ 6... Extreme narrowing The extreme narrowing or fast motion limit is when the orrelation time is so short that ω 0 τ <<. f this is the ase then + ωτ 0 and the setral density is just τ at all frequenies. The rate onstants for longitudinal and transverse relaxation are equal in this limit: R t ex. narrow = = τω ex. narrow R z As the orrelation time inreases out of the extreme narrowing limit the transverse relaxation rate onstant beomes greater than the longitudinal R 6 5

16 rate onstant as an be seen from Eq. [5]. For a 500 MHz setrometer the extreme narrowing limit imlies a orrelation time of muh less than 00 s. A small moleule in a nonvisous solvent might have a orrelation time of 0 s or less whih would lae the moleule easily in the extreme narrowing limit Sin diffusion The sin diffusion or slow motion limit is when ω 0 τ >> so that the setral density at zero frequeny is muh larger that that at the Larmor frequeny. n suh a ase the transverse relaxation rate onstant is dominated by the term whih deends on j(0) R t sin diff. = ωτ R The longitudinal rate onstant goes on getting smaller and smaller as the orrelation time inreases in rinile beoming vanishingly small in the limit ω 0 τ >>. There is thus a strong ontrast in the behaviour of the longitudinal and transverse rates. This differene omes about beause longitudinal relaxation requires motion at the Larmor frequeny but transverse relaxation an be aused either by motion at zero frequeny or at the Larmor frequeny. As was noted in Setion 5.. for a given frequeny the setral density is a maximum when ω 0 τ. Therefore the longitudinal relaxation rate onstant is a maximum when the orrelation time is suh that ω 0 τ. R t R z ( ωτ 0 ) The longitudinal and transverse relaxation rate onstants as a funtion of the orrelation time for fixed Larmor frequeny. For short orrelation times the two rate onstants are equal but as the orrelation time inreases the longitudinal rate onstant goes through a maximum and then falls away. n ontrast the transverse rate onstant ontinues to inrease. 6.. More than one sin with random field Suose that there are two sins eah exeriening a random field so that the hamiltonian is H t ω t ω t ω t ω ts ω ts ω ts = () + () + () + () + () + () whih an be written x x y y z z Sx x Sy y Sz z + () ( ) z S S S z + () 0 0 H ()= t F + F + F + S F + S F + F S 6 6

17 where [ ] ( ± ) S z S x y 0 S S S F ω t F ω t m iω t = () = () () The sin oerators in the hamiltonian are the only ones whih an affet sins oerators in the density oerator and likewise for S sin oerators. So the set of double ommutators [ ] ( ( q ) qq ) ( q) q [ ' ] A A α Sβ J ω q where α and S β are any and S sin oerators an be searated into two arts ( S A q ) A q β α J ' q q [ ] qq( ω ( ) ) [ ] Sq [ ] Sqq( ω ( ) ) q ( A Sq ) A S q α ' Sβ J + [ ] [6] The sub and suersrit labels and S on the oerators setral densities and frequenies distinguish those for the sins from those for the S sins. As an examle of using this searation suose that the relaxation behaviour of the oerator rodut + S z is required. From Eq. [6] [ ] q qq( ) + ( ) ω S A q A q ( ) z ' J q A Sq A S q + ( ) ( ) ( + [ [ Sz ]] J Sq ) ' Sqq ( ω ) q The first braket has already been alulated in the revious setion when the relaxation of + was onsidered; likewise the seond braket was also alulated above when the relaxation of S z was onsidered. So the relaxation rate onstants for + S z is just the sum of these two terms RS ( + z) = R ( + ) + RS ( z) ( ) ( ) ( S) ( S) = ( z ) j+ ( x ) + ( y ) j x y j S + + ω 0 ω ω ω ω ω ω n fat this aroah will work for any oerator with any number of terms in it. The overall relaxation rate onstant will be the sum of the rate onstants for eah individual oerator. Note that with an oerator suh as + S z the result is a mixture of longitudinal and transverse relaxation rates. A seond examle is a double quantum oherene: + S + RS ( + + ) = R ( + ) + RS ( + ) ( ) ( ) = ( z ) j+ ( x ) + ( y ) j ω 0 ω ω ω ( S ) ( S) ( S) + ( z ) j+ ( x ) + ( y ) j S ω 0 ω ω ω [ ] The relaxation does not deend on the setral density at the sum of the Larmor frequenies of and S. 6 7

18 x φ z The diole interation deends on the angles θ and φ whih are made by the vetor joining the two nulei. The stati magneti field is assumed to be along the z-axis. θ y 6.. Two sins with diolar interation The diole relaxation mehanism was desribed in Setion 5..; essentially one sin gives rise to a magneti field at a seond sin due to the magneti diole that the former ossesses. The interation deends on the distane between the two sins r and the angles that the vetor between them makes to the diretion of the alied magneti field. These angles are the usual θ and φ used in sherial olar oordinates and illustrated oosite. The diole hamiltonian an be searated into five terms with index q = 0 ± ± q A (q) F (q) [ sinθosθex m φ ] b[ sin θex miφ ] 0 S z z S + + S + b os θ ± S S ± S z ± + ± z b i ± ± where S b = µγγ 0 h πr in whih γ and γ S are the gyromagneti ratios. The angles and the distane an in rinile all hange randomly with time; it will be assumed that r is fixed and that the angles hange due to moleular tumbling. t is lear from the resene of oerators suh as z and + that the hamiltonian will give rise to longitudinal and transverse relaxation. However as will be seen the really imortant thing about the diole interation is that it gives rise to ross relaxation as was desribed in Setion 5... The details of how this omes about will be seen in this setion. The oerators A (q) evolve under the influene of the stati hamiltonian at harateristi frequenies as given in Eq. []. ( q) ( ex ih t A ex ih t ex i q ) ( q) ( 0 ) 0 ω t A [] = Taking the stati hamiltonian to be ω z + ω S S z the oerators A (q) searate out ( q into the A ) ( q eah with its assoiated frequeny ω ) aording to q ( q) A ( q) ω 0 0 S z z 0 +0 S + ω + ω S 0 S + ω ω S ± S z ± ±ω S ± S z ± ±ω 6 8

19 ± S ± ± m (ω + ω S ) The notation q = ±0 is used to indiate that the oerator with q = +0 is the adjoint of the one with q = 0 as required by Eq. [5]. The oerators with q = evolve at only one frequeny so the index is not required for these; likewise for q = Correlation funtions f the alulation is restrited to seular ontributions (Setion 6..) the required orrelation funtions are ( ) = () ( ) q q Gqq τ F t F t τ and these will be written as = () () = q q q Gqq τ F t F t g τ F g τ where g(τ) is the redued orrelation funtion. The average of the square of F (q) is omuted by averaging over the ensemble. f it is assumed that all angles are equally likely i.e. all orientations are equally robable then this averaging is equivalent to integrating over all angles θ and φ π π ( q) ( q) ( q) F = F F sinθdθdφ π θ = 0 φ = 0 where sinθdθd φ is the volume element in sherial olar o-ordinates. The division by π is needed for normalization as For q = 0 the integral is π π sinθdθdφ = π θ = 0 φ = 0 ( 0) F = b( ) b( ) π os θ os θ sinθdθdφ = 5 b π θ = 0 φ = 0 π The integral is tedious to evaluate by hand Mathematia or some similar rogram makes short work of it though. Similar alulations give the following results for the orrelation funtions and setral densities. q F ( q) J qq ω 0 5 b 5 ± 0 ± 0 b 0 b 0 b j ω b j( ω) b j( ω) From now on for brevity the setral densities J qq J q ( ω). will be written as ω 6 9

20 6... Relaxation of z The relaxation behaviour of the z oerator is determined by evaluating [ ] ( ( q ) q ω ) ( q) q [ ' z] A A J q where the restrition to seular terms has been assumed. The setral densities J q ( ω) have been omuted in the revious setion so the roblem redues to evaluating the double ommutator for all the relevant oerators. n evaluating the ommutators the results given in Setion 6.. will be useful. n addition the following ommutators will be needed: [ ]= [ ]= + S S S S + + z z + + S z Sz The following table shows the evaluation of all the double ommutators q 0 0 S S +0 6 S + S + z 0 6 S S z S S ( q) q [ A [ A ( ' z] ] A q ) [ Q result ] [ z z [ z z z] ] S z z 0 [ ] 0 [ [ ]] 6 [ S + S + ] 6 ( z Sz) [ + [ + ]] 6 [ S + S + ] 6 ( z + Sz) z + [ z z] [ S z + 0] 0 + z [ z z] [ S + z S z] z z [ z + z] [ S z 0] 0 z [ + z z] [ S z S + z] z + + [ ] [ S + + S ] ( z + Sz) [ + + ] [ S S + + ] ( z + Sz) [ ] [ ] [ ] [ ] [ ] [ ] S S S S S S S S z S S z The final result is therefore 6 ( z Sz) J ( ω + ωs)+ 6 ( z Sz) J 0 0( + ω ωs) + J z ( ω)+ J z ( ω) + ( z + Sz) J ( ω + ωs)+ ( z + Sz) J ( ω ω S ) where the aroriate setral densities have been inserted. Assuming that J ( ω)= J ( ± ω)= J ( ± ω) gives q q q q q { ( + )} ( ω ω ) ( ω ω ) J ω ω J ω J ω ω z S S { } Sz J + S 6 J0 S nserting the exliit exressions for the setral densities this gives 6 0

21 { ( + )} { 0 ( ω ω ) 6 5 ( ω ω )} z{ 0 ( ω ωs)+ 0 ( ω)+ 0 ( ω + ωs) } j( ω ω ) j( ω ω ) b j ω ω j ω j ω ω z S S S j + j z S S = j j j b z{ 0 S 0 S } [7] S + An exatly analogous alulation for the oerator S z gives { ( + )} ( ω ω ) ( ω ω ) S j ω ω j ω j ω ω b z S S S { } z 0 j + S 0 j S b whih is easily obtained by swaing the indies and S. [8] From Eq. [7] it is lear that the term multilying z is the relaxation rate onstant for z and the term multilying S z is the rate onstant for transfer from S z to z. Likewise in Eq. [8] the first term is the relaxation rate onstant for S z and the seond is the term for transfer of z to S z. These equations imly that z and S z both relax on their own and are interonverted by relaxation aording to the air of differential equations dz dsz = R z σssz = RS S z σsz where { + + ( + )} + + ( + ) R = j ω ω j ω j ω ω b S 0 S 0 S 0 S { } R = 0 j ω ωs 0 j ω 0 j ω ωs b σs = { 0 j( ω + ωs) 0 j( ω ωs) } b These are reisely the Solomon equations found in Setion 5... The differene is that it now it has been ossible to identify the ontributions to the self and ross relaxation rate onstants. Alternatively these equations an be written in terms of the J q (ω) as { + + ( + )} + + ( + ) R = J ω ω J ω J ω ω 6 0 S S { } { } RS = 6 J0 ω ωs J ωs J ω ωs σ = J ( ω + ω ) J ( ω ω ) S S 6 0 S [9] nterretation Referring to Eq. [9] it is seen that the setral densities J (ω S ) J 0 (ω ω S ) and J (ω + ω S ) all ontribute to the self relaxation of the and S sin. However only the latter two setral densities ontribute to the ross relaxation rate onstant. These two setral densities are assoiated with the terms in the hamiltonian S ± m alled the fli-flo terms as one sin flis one way and the other flis the other and the terms S ± ± alled the fli-fli and flo-flo terms. Diole relaxation is almost unique in having these terms resent and hene giving rise to ross relaxation. This is the reason why the resene of a nulear Overhauser effet whih requires there to be ross relaxation an almost always be assoiated with diolar relaxation. 6

22 Writing in the simlest setral density funtion gives the following exression for σ S σ S τ = 0 + ( ω + ω ) τ 0 + τ ω ω τ S S µ 0h π S γ γ 6 r f and S are both the same nulear seies say rotons it is safe to assume that (ω + ω S ) = ω and (ω ω S ) = 0 giving the simler exression σ S τ µ = τ 0h 0 +ωτ 0 π S γ γ 6 r The first thing to notie about this is that the ross-relaxation rate onstant is roortional to /r 6. Transient NOE exeriments an be used to measure the ross-relaxation rate onstant (see Setion 5..) and so rovided the orrelation time is known it is ossible to determine the distane. Tyially the orrelation time is found from indeendent exeriments suh as the measurement of the arbon- relaxation times and NOE enhanements. Alternatively relative distanes in the same moleule an simly be determined from the ratio of ross-relaxation rate onstants 6 σ r S PQ = 6 σ PQ rs This will work rovided it is valid to assume that the motion of both sin airs S and P Q an be modelled by a single orrelation time (i.e. a rigid moleule). Note that as only diolar relaxation ontributes to the rossrelaxation rate onstant it does not matter if other mehanisms are ating. These will alter the self relaxation rate onstants but not the rossrelaxation rate onstants. n the extreme narrowing limit ω S τ << and so σ S simlifies to ext. narrow µ γ γ S µ γ γ S σs = τ τ = τ 0 0h 0h 0 π 6 r π 6 r This is ositive resulting in a ositive NOE enhanement i.e. the sin reeiving the enhanement will have its magnetization inreased above the equilibrium value. This is tyial of the NOE in small moleules. n the sin diffusion limit the setral density at ω beomes negligible and so σ sin diff. S = µ τ 0h 0 π S γ γ 6 r This is negative giving rise to negative NOE enhanements i.e. the sin reeiving the enhanement will have its magnetization redued when omared to the equilibrium value. The self relaxation rate is R = τ τ µ 0h τ 0 + ωτ 0 ( ω + ωs) τ + π whih simlifies in the two limits to: S γ γ 6 r 6

23 R extr. narrow = { τ } µ 0h π S γ γ 6 r sin diff. S R = µ 0h γ γ τ π 6 0 r n Setion 5... it was shown that the steady-state NOE enhanement is σ η = S SS R f the relaxation is urely due to the diole interation between and S then this enhanement is + in the extreme narrowing limit and in the sin diffusion limit. Somewhere in between the NOE enhanement goes to zero when the two terms in the ross-relaxation rate onstant and equal i.e. τ µ γ γ S σ S = τ 0h 0 ωτ π r 6τ 5 whih is zero when = τ ie ωτ =.. + ωτ At a Larmor frequeny of 500 MHz this zero rossing orresonds to a orrelation time of about 0.5 ns. Suh a orrelation time is harateristi of a medium-sized moleule suh as a short etide or oligo-saharide dissolved in water Relaxation of z S z As in the revious setion the roess starts by evaluating [ ] ( ( q A A S J ) q ω ) ( q) q [ ' z z] q Following through the same roedure as before gives { } S J( ω )+ J( ω ) z z S there is no transfer to z or S z. The differential equation is thus d ( S z z) = RS ( zsz ) where { } R J ω J = + ( ω ) S S = { 0 j( ωs)+ 0 j( ω) } b τ τ = + µ 0h 0 + ωτ S 0 ωτ π + Note that in the sin diffusion limit this rate tends to zero. S γ γ 6 r 6... Relaxation of + Evaluating 6

24 gives 6 ( q) [ ] q( ω ) ( q) [ ( q ) ' + ] A A J q +{ 8 0 S 0( S)+ ( )+ ( + S) } J ( 0)+ J ( ω )+ J ω ω J ω J ω ω There are no transfer terms to other oerators. The differential equation is thus d+ t = R + where t R = { 8 J0( 0)+ J( ωs)+ J0( ω ωs)+ J( ω)+ J( ω + ωs) } = { 8 J0( 0)+ J( ωs) }+ { 6 J0( ω ωs)+ J( ω)+ J( ω + ωs) } = { 8 J0( 0) + J( ω S )}+ R Just as in the random field ase the transverse relaxation rate onstant searates into two arts. The seond art is just half of the longitudinal rate onstant R. The first art has a term whih deends on the setral density at zero frequeny and a term whih deends on the setral density at the Larmor frequeny of sin S. This latter term is a little unusual. t says that the relaxation of transverse magnetization of sin is affeted by the setral density at the Larmor frequeny of sin S. The interretation of this is that this latter setral density is assoiated with roesses whih involve the fliing of sin S (terms suh as z S ± ). As is (diolar) ouled to S when the sin state of S hanges the field seen by the sin hanges and so the reession frequeny of the sin oherene hanges. These sin flis thus give rise to a onstantly varying loal field at sin whih auses loss of hase oherene amongst the sins whih are ontributing to the oherene. The term is seular from the oint of view of the sin as no energy is transferred to or from it. n terms of the redue setral densities the relaxation rate onstant is { + + ( + )} t R = j( 0)+ j( ω )+ j ω ω j ω j ω ω b 0 0 S 0 S 0 0 S Transverse ross relaxation n ontrast to the relaxation behaviour of z and S z there is no ross relaxation transfer between + and S + ; i.e. no transverse ross relaxation. Suh a transfer would ome about through the following double ommutators q= 0 : [ S z z [ S + + ]]= [ S z z S z + ]= 8 S+ [ S + [ S z z + ]]= [ S z z S + z]= 8 S+ q =± : [ Sz [ zs+ + ]]= [ Sz S + + ]= S+ [ S z + [ S z + ]]=[ S z + S z z]= S+ However none of these terms ontribute as they are non-seular. For

25 examle in the first ommutator the frequeny assoiated with z S z is 0 whereas that assoiated with S + is ( ω + ω S ). The resulting term is therefore raidly osillating and thus ineffetive at ausing relaxation. There are two ases in whih these terms beome seular. Firstly when ω = ω S i.e. when the two sins are degenerate. Seondly when it aears that the two sins have the same Larmor frequeny. The first ase is not of muh interest as it would be imossible to detet the transfer of magnetization between degenerate sins. The latter ase is imortant as it is ossible to make it aear that two sins have the same Larmor frequeny by sin loking their magnetization. The simlest sin loking exeriment involves a non-seletive 90 ulse about the y-axis followed by a eriod during whih a strong radiofrequeny field is alied along the x-axis. The field is suffiiently strong that for both sins (whih have different offsets) the effetive field lies very lose to the x-axis. The magnetization whih is along the axis of the alied radiofrequeny field (here x) exerienes no rotation from that field; viewed in the laboratory frame the magnetization is reessing at the same rate as the rotating field. n the rotating frame both aear stati. The magnetization is said to be sin loked. Sine the magnetization from the two sins does not diverge while the sin-loking field is on it aears as if they have the same offset or Larmor frequeny. Thus transverse ross-relaxation an take lae. The results of the ross relaxation an be observed by removing the sin-loking field and then observing the FD in the usual way. The exeriments used to observe suh transverse NOE enhanements as they are alled are generally referred to as ROESY exeriments (or sometime CAMELSPN). Both one- and two-dimensional versions of the exeriment are available. Using the above double ommutators the ross-relaxation rate onstant is omuted as { } tr σs = J( ω )+ J ( ) or in terms of the redued setral densities tr σs = { 0 j( ω )+ 0 j( 0) } b The imortant oint about this rate onstant is that for all values of the orrelation time it is always ositive. Therefore the NOE enhanements are always ositive and there is no range of orrelation times for whih the enhanement rosses zero. ROESY has therefore roved to be a useful tehnique for obtaining distane information on moleules for whih beause of their orrelation times the onventional NOE is lose to zero. (a) (b) S (y) 90 (y) τ m x τ m x One-dimensional ROESY exeriment. The ombination of the seletive 80 ulse and the 90 ulse in (a) generates a state in whih the sin is along x and the S sin is along x. During the eriod of sin loking transverse ross relaxation an take lae. Sequene (b) generates the referene setrum without ross relaxation; it is subtrated from (a) to give the usual differene setrum. τ t t mix Two-dimensional ROESY exeriment whih uses a eriod of sin loking during the mixing time. Off-resonane ROESY n ratie it is diffiult to use a suffiiently strong sin loking field that the sin loking axis is in the transverse lane. Rather the axis will by tilted u into the xz-lane. The tile angle of the effetive field θ deends on the offset Ω and the radio frequeny field strength ω. Ω θ ω Ω eff 6 5

26 ω ω Ω tanθ = sinθ = os θ = Ωeff = ω + Ω Ω Ω Ω eff When the sin loking field is alied omonents of the magnetization whih are not arallel to the field are quikly dehased and lost. Only the art along the sin lok axis need be onsidered. Magnetization whih is sin loked about an axis tilted at angle θ will have a transverse omonent roortional to sinθ this relaxes at the aroriate rate for transverse magnetization and a longitudinal omonent roortional to osθ this relaxes at the longitudinal rate. The rate onstant for the relaxation of sin loked magnetization thus deends on the rate onstants for both transverse and longitudinal magnetization and the angle θ. n an ROE exeriment the two sins whih are ross-relaxing one another will have different offsets and thus have different tilt angles θ and θ S. t is the transverse omonents whih undergo transverse ross relaxation and the longitudinal omonents whih undergo normal ross relaxation. To work out the effetive ross relaxation rate of the sin loked magnetization it is thus neessary to onsider the two omonents searately. z (a) eff (b) S S b θ d a θ S x Geometri onstrutions for the ase where the sin loking axes of and S are tilted to different extents. Case (a) is for transverse and (b) for longitudinal relaxation. Referring to diagram (a) the transverse omonent of labelled a is roortional to sinθ ; this omonent ross relaxes with S with a rate onstant σ t. Any S magnetization generated in this roess aears along the x-axis but only the omonent along the sin loking axis of S will survive. This omonent labelled b is roortional to sinθ S. The effetive transverse ross relaxation rate is thus sinθ sinθ S σ t. Diagram (b) shows the longitudinal omonents. is the longitudinal omonent of the magnetization and is roortional to osθ. This omonent ross relaxes with rate onstant σ l to give S magnetization. The omonent of this magnetization along the S sin lok axis d is roortional to osθ S. The effetive longitudinal ross relaxation rate onstant is thus osθ osθ S σ l. Overall therefore the ross relaxation rate onstant of the sin loked magnetization is given by σ tilted = sinθ sinθ S σ t + osθ osθ S σ l Using a similar geometri onstrution it is ossible to show that the self 6 6

27 relaxation rate onstant is and likewise for S. R tilted = sin θ R t + os θ R l Relaxation of + S z Evaluating gives ( q) [ A A S ] Jq( ω ) ( q) [ ( q ) ' + z] q { ( + )} S J( 0)+ J ω ω J ω J ω ω + z S S there is no transfer to other terms. The differential equation is thus where ds + z = R + Sz S { + + ( + )} R S = J J S J J + z 8 0( 0)+ 0 ω ω ω ω ωs = j( 0)+ j ω ω j ω j ω ω b { 0 0 ( S)+ 0 ( )+ 0 ( + S) } t is interesting to omare this with the relaxation time onstant for + { + + ( + )} t R = J ( 0)+ J ( ω )+ J ω ω J ω J ω ω 8 0 S 0 S S t is immediately lear that the antihase term relaxes more slowly than the inhase term; the differene is the term J( ω S) whih does not ontribute to the relaxation of the antihase term. The effet of this term was ommented on in the revious setion. There is another imortant feature of the differene between the relaxation of inhase and antihase terms. For the antihase term any additional longitudinal relaxation of the S sin R addn. S ontributes diretly to the relaxation of the oerator rodut + S z ; the reasons for this were disussed in Setion 6... R = diole R + addn. R + Sz + Sz S Suh longitudinal relaxation of sin S does not ontribute to the relaxation of +. As suh extra relaxation of S is often found in ratie it is usually the ase that antihase terms relax more quikly than inhase terms Relaxation of multile quantum terms + z Double quantum Evaluating gives ( q) [ A A S ] Jq( ω ) ( q) [ ( q ) ' + + ] q S J( ω )+ J( ω )+ J( ω + ω ) + +{ S S } 6 7

Subject: Modeling of Thermal Rocket Engines; Nozzle flow; Control of mass flow. p c. Thrust Chamber mixing and combustion

Subject: Modeling of Thermal Rocket Engines; Nozzle flow; Control of mass flow. p c. Thrust Chamber mixing and combustion 16.50 Leture 6 Subjet: Modeling of Thermal Roket Engines; Nozzle flow; Control of mass flow Though onetually simle, a roket engine is in fat hysially a very omlex devie and diffiult to reresent quantitatively

More information

Panos Kouvelis Olin School of Business Washington University

Panos Kouvelis Olin School of Business Washington University Quality-Based Cometition, Profitability, and Variable Costs Chester Chambers Co Shool of Business Dallas, TX 7575 hamber@mailosmuedu -768-35 Panos Kouvelis Olin Shool of Business Washington University

More information

Proportional-Integral-Derivative PID Controls

Proportional-Integral-Derivative PID Controls Proortional-Integral-Derivative PID Controls Dr M.J. Willis Det. of Chemial and Proess Engineering University of Newastle e-mail: mark.willis@nl.a.uk Written: 7 th November, 998 Udated: 6 th Otober, 999

More information

Risk Analysis in Water Quality Problems. Souza, Raimundo 1 Chagas, Patrícia 2 1,2 Departamento de Engenharia Hidráulica e Ambiental

Risk Analysis in Water Quality Problems. Souza, Raimundo 1 Chagas, Patrícia 2 1,2 Departamento de Engenharia Hidráulica e Ambiental Risk Analysis in Water Quality Problems. Downloaded from aselibrary.org by Uf - Universidade Federal Do Ceara on 1/29/14. Coyright ASCE. For ersonal use only; all rights reserved. Souza, Raimundo 1 Chagas,

More information

Solutions of burnt-bridge models for molecular motor transport

Solutions of burnt-bridge models for molecular motor transport PHYSICAL REVIEW E 75, 0390 2007 Solutions of burnt-bridge models for moleular motor transort Alexander Yu. Morozov, Ekaterina Pronina, and Anatoly B. Kolomeisky Deartment of Chemistry, Rie University,

More information

CHAPTER 16. Basic Concepts. Basic Concepts. The Equilibrium Constant. Reaction Quotient & Equilibrium Constant. Chemical Equilibrium

CHAPTER 16. Basic Concepts. Basic Concepts. The Equilibrium Constant. Reaction Quotient & Equilibrium Constant. Chemical Equilibrium Proerties of an Equilibrium System CHAPTER 6 Chemial Equilibrium Equilibrium systems are DYNAMIC (in onstant motion) REVERSIBLE an be aroahed from either diretion Pink to blue Co(H O) 6 Cl ---> > Co(H

More information

CONSTRUCTION OF MIXED SAMPLING PLAN WITH DOUBLE SAMPLING PLAN AS ATTRIBUTE PLAN INDEXED THROUGH (MAPD, MAAOQ) AND (MAPD, AOQL)

CONSTRUCTION OF MIXED SAMPLING PLAN WITH DOUBLE SAMPLING PLAN AS ATTRIBUTE PLAN INDEXED THROUGH (MAPD, MAAOQ) AND (MAPD, AOQL) Global J. of Arts & Mgmt., 0: () Researh Paer: Samath kumar et al., 0: P.07- CONSTRUCTION OF MIXED SAMPLING PLAN WITH DOUBLE SAMPLING PLAN AS ATTRIBUTE PLAN INDEXED THROUGH (MAPD, MAAOQ) AND (MAPD, AOQL)

More information

EXTENDED MATRIX CUBE THEOREMS WITH APPLICATIONS TO -THEORY IN CONTROL

EXTENDED MATRIX CUBE THEOREMS WITH APPLICATIONS TO -THEORY IN CONTROL MATHEMATICS OF OPERATIONS RESEARCH Vol. 28, No. 3, August 2003,. 497 523 Printed in U.S.A. EXTENDED MATRIX CUBE THEOREMS WITH APPLICATIONS TO -THEORY IN CONTROL AHARON BEN-TAL, ARKADI NEMIROVSKI, and CORNELIS

More information

CSC321: 2011 Introduction to Neural Networks and Machine Learning. Lecture 11: Bayesian learning continued. Geoffrey Hinton

CSC321: 2011 Introduction to Neural Networks and Machine Learning. Lecture 11: Bayesian learning continued. Geoffrey Hinton CSC31: 011 Introdution to Neural Networks and Mahine Learning Leture 11: Bayesian learning ontinued Geoffrey Hinton Bayes Theorem, Prior robability of weight vetor Posterior robability of weight vetor

More information

SOLVED QUESTIONS 1 / 2. in a closed container at equilibrium. What would be the effect of addition of CaCO 3 on the equilibrium concentration of CO 2?

SOLVED QUESTIONS 1 / 2. in a closed container at equilibrium. What would be the effect of addition of CaCO 3 on the equilibrium concentration of CO 2? SOLVED QUESTIONS Multile Choie Questions. and are the veloity onstants of forward and bakward reations. The equilibrium onstant k of the reation is (A) (B) (C) (D). Whih of the following reations will

More information

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E')

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') 22.54 Neutron Interations and Appliations (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') Referenes -- J. R. Lamarsh, Introdution to Nulear Reator Theory (Addison-Wesley, Reading, 1966),

More information

The Procedure of Finding the Stress-Energy. Tensor and Equations of Vector Field of Any Form

The Procedure of Finding the Stress-Energy. Tensor and Equations of Vector Field of Any Form Advaned Studies in Theoretial Phsis Vol. 8, 14, no. 18, 771-779 HIKARI Ltd, www.m-hikari.om htt://d.doi.org/1.1988/ast.14.4711 The Proedure of Finding the Stress-Energ Tensor and Equations of Vetor Field

More information

3. THE SOLUTION OF TRANSFORMATION PARAMETERS

3. THE SOLUTION OF TRANSFORMATION PARAMETERS Deartment of Geosatial Siene. HE SOLUION OF RANSFORMAION PARAMEERS Coordinate transformations, as used in ratie, are models desribing the assumed mathematial relationshis between oints in two retangular

More information

Non-Markovian study of the relativistic magnetic-dipole spontaneous emission process of hydrogen-like atoms

Non-Markovian study of the relativistic magnetic-dipole spontaneous emission process of hydrogen-like atoms NSTTUTE OF PHYSCS PUBLSHNG JOURNAL OF PHYSCS B: ATOMC, MOLECULAR AND OPTCAL PHYSCS J. Phys. B: At. Mol. Opt. Phys. 39 ) 7 85 doi:.88/953-75/39/8/ Non-Markovian study of the relativisti magneti-dipole spontaneous

More information

EE451/551: Digital Control. Relationship Between s and z Planes. The Relationship Between s and z Planes 11/10/2011

EE451/551: Digital Control. Relationship Between s and z Planes. The Relationship Between s and z Planes 11/10/2011 /0/0 EE45/55: Digital Control Chater 6: Digital Control System Design he Relationshi Between s and Planes As noted reviously: s j e e e e r s j where r e and If an analog system has oles at: s n jn a jd

More information

EconS 503 Homework #8. Answer Key

EconS 503 Homework #8. Answer Key EonS 503 Homework #8 Answer Key Exerise #1 Damaged good strategy (Menu riing) 1. It is immediate that otimal rie is = 3 whih yields rofits of ππ = 3/ (the alternative being a rie of = 1, yielding ππ =

More information

Chapter 9. The excitation process

Chapter 9. The excitation process Chapter 9 The exitation proess qualitative explanation of the formation of negative ion states Ne and He in He-Ne ollisions an be given by using a state orrelation diagram. state orrelation diagram is

More information

Fast, Approximately Optimal Solutions for Single and Dynamic MRFs

Fast, Approximately Optimal Solutions for Single and Dynamic MRFs Fast, Aroximately Otimal Solutions for Single and Dynami MRFs Nikos Komodakis, Georgios Tziritas University of Crete, Comuter Siene Deartment {komod,tziritas}@sd.uo.gr Nikos Paragios MAS, Eole Centrale

More information

MOLECULAR ORBITAL THEORY- PART I

MOLECULAR ORBITAL THEORY- PART I 5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.

More information

2.2 BUDGET-CONSTRAINED CHOICE WITH TWO COMMODITIES

2.2 BUDGET-CONSTRAINED CHOICE WITH TWO COMMODITIES Essential Miroeonomis -- 22 BUDGET-CONSTRAINED CHOICE WITH TWO COMMODITIES Continuity of demand 2 Inome effets 6 Quasi-linear, Cobb-Douglas and CES referenes 9 Eenditure funtion 4 Substitution effets and

More information

Fundamental Theorem of Calculus

Fundamental Theorem of Calculus Chater 6 Fundamental Theorem of Calulus 6. Definition (Nie funtions.) I will say that a real valued funtion f defined on an interval [a, b] is a nie funtion on [a, b], if f is ontinuous on [a, b] and integrable

More information

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge The Conept of Mass as Interfering Photons, and the Originating Mehanism of Gravitation D.T. Froedge V04 Formerly Auburn University Phys-dtfroedge@glasgow-ky.om Abstrat For most purposes in physis the onept

More information

A General Approach for Analysis of Actuator Delay Compensation Methods for Real-Time Testing

A General Approach for Analysis of Actuator Delay Compensation Methods for Real-Time Testing The th World Conferene on Earthquake Engineering Otober -7, 8, Beijing, China A General Aroah for Analysis of Atuator Delay Comensation Methods for Real-Time Testing Cheng Chen and James M. Riles Post-dotoral

More information

Physics of Relaxation. Outline

Physics of Relaxation. Outline Physis of Relaxation Weiguo Li Outline Fundamental relaxation Mehanisms Magneti dipole-dipole oupling» Stati oupling» Dynami oupling Frequeny dependene of relaxation Rate Temperature dependene of relaxation

More information

Towards Applications of Bubble Acceleration: Bubble-driven Free-Electron-Lasers F. Grüner (LMU+MPQ), July 15, 2005, Photonics Lecture, Prof.

Towards Applications of Bubble Acceleration: Bubble-driven Free-Electron-Lasers F. Grüner (LMU+MPQ), July 15, 2005, Photonics Lecture, Prof. Towards Aliations of Bubble Aeleration: Bubble-driven Free-Eletron-Lasers F. Grüner (LMU+MPQ), July 5, 5, Photonis Leture, Prof. Krausz Motivation for FELs What is an FEL? Undulators and Wigglers SASE

More information

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

23.1 Tuning controllers, in the large view Quoting from Section 16.7: Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output

More information

Relativistic Dynamics

Relativistic Dynamics Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable

More information

Green s function for the wave equation

Green s function for the wave equation Green s funtion for the wave equation Non-relativisti ase January 2019 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 43 and 44): 1 2 A 2 2 2 A = µ 0

More information

Simple theory of Plasmons, surface plasmons, Polaritons and all that

Simple theory of Plasmons, surface plasmons, Polaritons and all that Simle theory of Plasmons, surfae lasmons, Polaritons and all that First, onsier a transverse eletromagneti wave in onduting (J=sE) homogeneous non magneti medium Consider the lane wave going uwards i.e.

More information

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six

More information

Summary of EUV Optics Contamination Modeling Meeting

Summary of EUV Optics Contamination Modeling Meeting Summary of EUV Otis Contamination Modeling Meeting IEUVI Otis Contamination/Lifetime TWG San Jose, CA, USA February 28, 2008 Andrea Wüest, SEMATECH Shannon Hill, NIST Lennie Klebanoff, Sandia Nat. Lab.

More information

USING GENETIC ALGORITHMS FOR OPTIMIZATION OF TURNING MACHINING PROCESS

USING GENETIC ALGORITHMS FOR OPTIMIZATION OF TURNING MACHINING PROCESS Journal of Engineering Studies and Researh Volume 19 (2013) No. 1 47 USING GENETIC ALGORITHMS FOR OPTIMIZATION OF TURNING MACHINING PROCESS DUSAN PETKOVIC 1, MIROSLAV RADOVANOVIC 1 1 University of Nis,

More information

Real-time Hand Tracking Using a Sum of Anisotropic Gaussians Model

Real-time Hand Tracking Using a Sum of Anisotropic Gaussians Model Real-time Hand Traking Using a Sum of Anisotroi Gaussians Model Srinath Sridhar 1, Helge Rhodin 1, Hans-Peter Seidel 1, Antti Oulasvirta 2, Christian Theobalt 1 1 Max Plank Institute for Informatis Saarbrüken,

More information

General Equilibrium. What happens to cause a reaction to come to equilibrium?

General Equilibrium. What happens to cause a reaction to come to equilibrium? General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember

More information

Simple Cyclic loading model based on Modified Cam Clay

Simple Cyclic loading model based on Modified Cam Clay Simle li loading model based on Modified am la Imlemented in RISP main rogram version 00. and higher B Amir Rahim, The RISP onsortium Ltd Introdution This reort resents a simle soil model whih rovides

More information

INFORMATION TRANSFER THROUGH CLASSIFIERS AND ITS RELATION TO PROBABILITY OF ERROR

INFORMATION TRANSFER THROUGH CLASSIFIERS AND ITS RELATION TO PROBABILITY OF ERROR IFORMATIO TRAFER TROUG CLAIFIER AD IT RELATIO TO PROBABILITY OF ERROR Deniz Erdogmus, Jose C. Prini Comutational euroengineering Lab (CEL, University of Florida, Gainesville, FL 6 [deniz,rini]@nel.ufl.edu

More information

A special reference frame is the center of mass or zero momentum system frame. It is very useful when discussing high energy particle reactions.

A special reference frame is the center of mass or zero momentum system frame. It is very useful when discussing high energy particle reactions. High nergy Partile Physis A seial referene frame is the enter of mass or zero momentum system frame. It is very useful when disussing high energy artile reations. We onsider a ollision between two artiles

More information

Determination of the reaction order

Determination of the reaction order 5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order

More information

On a Networked Automata Spacetime and its Application to the Dynamics of States in Uniform Translatory Motion

On a Networked Automata Spacetime and its Application to the Dynamics of States in Uniform Translatory Motion On a Networked Automata Saetime and its Aliation to the Dynamis of States in Uniform Translatory Motion Othon Mihail Deartment of Comuter Siene, University of Liverool, Liverool, UK Othon.Mihail@liverool.a.uk

More information

arxiv: v2 [econ.gn] 12 Oct 2018

arxiv: v2 [econ.gn] 12 Oct 2018 A new and stable estimation method of ountry eonomi fitness and rodut omlexity Vito D. P. Servedio 1, Paolo Buttà 2, Dario Mazzilli 3, Andrea Tahella 4, Luiano Pietronero 3 arxiv:1807.276v2 [eon.gn] 12

More information

NONLINEAR MODEL: CELL FORMATION

NONLINEAR MODEL: CELL FORMATION ational Institute of Tehnology Caliut Deartment of Mehanial Engineering OLIEAR MODEL: CELL FORMATIO System Requirements and Symbols 0-Integer Programming Formulation α i - umber of interell transfers due

More information

Temperature-Gradient-Driven Tearing Modes

Temperature-Gradient-Driven Tearing Modes 1 TH/S Temperature-Gradient-Driven Tearing Modes A. Botrugno 1), P. Buratti 1), B. Coppi ) 1) EURATOM-ENEA Fusion Assoiation, Frasati (RM), Italy ) Massahussets Institute of Tehnology, Cambridge (MA),

More information

INCOME AND SUBSTITUTION EFFECTS

INCOME AND SUBSTITUTION EFFECTS 3076-C5.df 3/12/04 10:56 AM Page 121 121 Chater 5 INCOME AND SUBSTITUTION EFFECTS In this hater we will use the utility-maimization model to study how the quantity of a good that an individual hooses is

More information

1 sin 2 r = 1 n 2 sin 2 i

1 sin 2 r = 1 n 2 sin 2 i Physis 505 Fall 005 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.5, 7.8, 7.16 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with

More information

Application of Drucker-Prager Plasticity Model for Stress- Strain Modeling of FRP Confined Concrete Columns

Application of Drucker-Prager Plasticity Model for Stress- Strain Modeling of FRP Confined Concrete Columns Available online at www.sienediret.om Proedia Engineering 14 (211) 687 694 The Twelfth East Asia-Paifi Conferene on Strutural Engineering and Constrution Aliation of Druker-Prager Plastiity Model for Stress-

More information

Line Radiative Transfer

Line Radiative Transfer http://www.v.nrao.edu/ourse/astr534/ineradxfer.html ine Radiative Transfer Einstein Coeffiients We used armor's equation to estimate the spontaneous emission oeffiients A U for À reombination lines. A

More information

Revista Mexicana de Ingeniería Química

Revista Mexicana de Ingeniería Química Revista exiana de Ingeniería Químia Vol. 8 o. 3 (009) 3-43 ERIVTIO PPLITIO OF THE STEF-XWELL EQUTIOS ESRROLLO Y PLIIÓ E LS EUIOES E STEF-XWELL bstrat Stehen Whitaker * eartment of hemial Engineering &

More information

Velocity Addition in Space/Time David Barwacz 4/23/

Velocity Addition in Space/Time David Barwacz 4/23/ Veloity Addition in Spae/Time 003 David arwaz 4/3/003 daveb@triton.net http://members.triton.net/daveb Abstrat Using the spae/time geometry developed in the previous paper ( Non-orthogonal Spae- Time geometry,

More information

Supplementary information for: All-optical signal processing using dynamic Brillouin gratings

Supplementary information for: All-optical signal processing using dynamic Brillouin gratings Supplementary information for: All-optial signal proessing using dynami Brillouin gratings Maro Santagiustina, Sanghoon Chin 2, Niolay Primerov 2, Leonora Ursini, Lu Thévena 2 Department of Information

More information

Hankel Optimal Model Order Reduction 1

Hankel Optimal Model Order Reduction 1 Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both

More information

Theory of the Integer and Fractional Quantum Hall Effects

Theory of the Integer and Fractional Quantum Hall Effects Theory of the Integer and Frational Quantum all Effets Shosuke SASAKI Center for Advaned igh Magneti Field Siene, Graduate Shool of Siene, Osaka University, - Mahikaneyama, Toyonaka, Osaka 56-, Jaan Prefae

More information

Aharonov-Bohm effect. Dan Solomon.

Aharonov-Bohm effect. Dan Solomon. Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that

More information

Berry s phase for coherent states of Landau levels

Berry s phase for coherent states of Landau levels Berry s phase for oherent states of Landau levels Wen-Long Yang 1 and Jing-Ling Chen 1, 1 Theoretial Physis Division, Chern Institute of Mathematis, Nankai University, Tianjin 300071, P.R.China Adiabati

More information

Controller Design Based on Transient Response Criteria. Chapter 12 1

Controller Design Based on Transient Response Criteria. Chapter 12 1 Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors

More information

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install

More information

A model for measurement of the states in a coupled-dot qubit

A model for measurement of the states in a coupled-dot qubit A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:

More information

Lecture 15 (Nov. 1, 2017)

Lecture 15 (Nov. 1, 2017) Leture 5 8.3 Quantum Theor I, Fall 07 74 Leture 5 (Nov., 07 5. Charged Partile in a Uniform Magneti Field Last time, we disussed the quantum mehanis of a harged partile moving in a uniform magneti field

More information

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013 Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it

More information

Research on Components Characteristics of Single Container Friction Hoist

Research on Components Characteristics of Single Container Friction Hoist 6 3 rd International Conferene on Mehanial, Industrial, and Manufaturing Engineering (MIME 6) ISB: 978--6595-33-7 Researh on Comonents Charateristis of Single Container Frition oist irong Xie, Jing Cheng,

More information

Takeshi Kurata Jun Fujiki Katsuhiko Sakaue. 1{1{4 Umezono, Tsukuba-shi, Ibaraki , JAPAN. fkurata, fujiki,

Takeshi Kurata Jun Fujiki Katsuhiko Sakaue. 1{1{4 Umezono, Tsukuba-shi, Ibaraki , JAPAN. fkurata, fujiki, Ane Eiolar Geometry via Fatorization Method akeshi Kurata Jun Fujiki Katsuhiko Sakaue Eletrotehnial Laboratory {{4 Umezono, sukuba-shi, Ibaraki 305-8568, JAPAN fkurata, fujiki, sakaueg@etl.go.j Abstrat

More information

Anelastic MHD Equations : Theory

Anelastic MHD Equations : Theory Anelasti MHD Equations : Theory Abhishek Kr. Srivastava Armagh Observatory Northern Ireland Different zones of the Sun Comlex magneti field is the key behind all tyes of solar ativity and formation of

More information

The synergy between the insect-inspired claws and adhesive pads increases the attachment ability on various rough surfaces

The synergy between the insect-inspired claws and adhesive pads increases the attachment ability on various rough surfaces The synergy between the inset-insired laws and adhesive ads inreases the attahment ability on various rough surfaes Yi Song 1, 2, Zhendong Dai 1, Zhouyi Wang 1, Aihong Ji 1 1, 3,*, and Stanislav Gorb 1

More information

arxiv: v2 [astro-ph] 1 Oct 2007

arxiv: v2 [astro-ph] 1 Oct 2007 Mon. Not. R. Astron. So. 000, 1 25 (2007) Printed 20 June 2018 (MN LATEX style file v2.2) Long-term tidal evolution of short-eriod lanets with omanions arxiv:0706.0224v2 [astro-h] 1 Ot 2007 Rosemary A.

More information

ELECTROMAGNETIC WAVES

ELECTROMAGNETIC WAVES ELECTROMAGNETIC WAVES Now we will study eletromagneti waves in vauum or inside a medium, a dieletri. (A metalli system an also be represented as a dieletri but is more ompliated due to damping or attenuation

More information

The Effectiveness of the Linear Hull Effect

The Effectiveness of the Linear Hull Effect The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports

More information

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry 9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived

More information

System Modeling Concepts

System Modeling Concepts 1. E+2 1. E+1 1. E+4 1. E+3. 1 1.. 1 1. 2. 3. 4. 5. L i n e a r F r e q u e n y ( H z ) 2. 3. 4. 5. L i n e a r F r e q u e n y ( H z ) 2 1 1 2 2 1 2 1 Strutural Dynami odeling ehniques & odal nalysis

More information

Problem 1(a): At equal times or in the Schrödinger picture the quantum scalar fields ˆΦ a (x) and ˆΠ a (x) satisfy commutation relations

Problem 1(a): At equal times or in the Schrödinger picture the quantum scalar fields ˆΦ a (x) and ˆΠ a (x) satisfy commutation relations PHY 396 K. Solutions for homework set #7. Problem 1a: At equal times or in the Shrödinger iture the quantum salar fields ˆΦ a x and ˆΠ a x satisfy ommutation relations ˆΦa x, ˆΦ b y 0, ˆΠa x, ˆΠ b y 0,

More information

3 Tidal systems modelling: ASMITA model

3 Tidal systems modelling: ASMITA model 3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level

More information

Introduction to Relaxation Theory James Keeler

Introduction to Relaxation Theory James Keeler EUROMAR Zürich, 24 Introduction to Relaxation Theory James Keeler University of Cambridge Department of Chemistry What is relaxation? Why might it be interesting? relaxation is the process which drives

More information

REVIEW OF PROPELLER-WING AERODYNAMIC INTERFERENCE

REVIEW OF PROPELLER-WING AERODYNAMIC INTERFERENCE 24 TH INTERNATIONA CONGRESS OF THE AERONAUTICA SCIENCES REVIEW OF PROPEER-WING AERODYNAMIC INTERFERENCE..M. Veldhuis Delft University of Tehnology Keywords: roeller-wing interation, numerial modeling,

More information

max min z i i=1 x j k s.t. j=1 x j j:i T j

max min z i i=1 x j k s.t. j=1 x j j:i T j AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be

More information

arxiv:gr-qc/ v2 6 Feb 2004

arxiv:gr-qc/ v2 6 Feb 2004 Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this

More information

WEIGHTED MAXIMUM OVER MINIMUM MODULUS OF POLYNOMIALS, APPLIED TO RAY SEQUENCES OF PADÉ APPROXIMANTS

WEIGHTED MAXIMUM OVER MINIMUM MODULUS OF POLYNOMIALS, APPLIED TO RAY SEQUENCES OF PADÉ APPROXIMANTS WEIGHTED MAXIMUM OVER MINIMUM MODULUS OF POLYNOMIALS, APPLIED TO RAY SEQUENCES OF PADÉ APPROXIMANTS D.S. LUBINSKY Abstrat. Let a 0; " > 0. We use otential theory to obtain a shar lower bound for the linear

More information

arxiv:hep-ph/ v1 10 Feb 2003

arxiv:hep-ph/ v1 10 Feb 2003 Two-loo HTL Thermodynamis with uarks Jens O. Andersen Institute for Theoretial hysis, University of Utreht, Leuvenlaan 4, 584 CE Utreht, The Netherlands Emmanuel etitgirard hysis Deartment, Ohio State

More information

11.4 Molecular Orbital Description of the Hydrogen Molecule Electron Configurations of Homonuclear Diatomic Molecules

11.4 Molecular Orbital Description of the Hydrogen Molecule Electron Configurations of Homonuclear Diatomic Molecules Chap Moleular Eletroni Struture Table of Contents. The orn-oppenheimer pproximation -. The Hydrogen Moleule Ion.3 Calulation of the Energy of the Hydrogen Moleule Ion.4 Moleular Orbital Desription of the

More information

INTUITIONISTIC FUZZY SOFT MATRIX THEORY IN MEDICAL DIAGNOSIS USING MAX-MIN AVERAGE COMPOSITION METHOD

INTUITIONISTIC FUZZY SOFT MATRIX THEORY IN MEDICAL DIAGNOSIS USING MAX-MIN AVERAGE COMPOSITION METHOD Journal of Theoretial and lied Information Tehnology 10 th Setember 014. Vol. 67 No.1 005-014 JTIT & LLS. ll rights reserved. ISSN: 199-8645 www.atit.org E-ISSN: 1817-3195 INTUITIONISTI FUZZY SOFT MTRIX

More information

Announcements. Lecture 5 Chapter. 2 Special Relativity. The Doppler Effect

Announcements. Lecture 5 Chapter. 2 Special Relativity. The Doppler Effect Announements HW1: Ch.-0, 6, 36, 41, 46, 50, 51, 55, 58, 63, 65 *** Lab start-u meeting with TA yesterday; useful? *** Lab manual is osted on the ourse web *** Physis Colloquium (Today 3:40m anelled ***

More information

Critical Reflections on the Hafele and Keating Experiment

Critical Reflections on the Hafele and Keating Experiment Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As

More information

top coils plasma bottom coils

top coils plasma bottom coils 1 Stabiliation of the ertial Mode in Tokamaks by Loalied Nonaxisymmetri Fields A. Reiman Prineton Plasma Physis Laboratory, Prineton, NJ 8543 reiman@l.gov Abstrat. ertial instability of a tokamak lasma

More information

Brazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle

Brazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle Brazilian Journal of Physis, vol. 9, no. 3, September, 1999 51 Classial and Quantum Mehanis of a Charged Partile in Osillating Eletri and Magneti Fields V.L.B. de Jesus, A.P. Guimar~aes, and I.S. Oliveira

More information

Detection and Estimation Theory

Detection and Estimation Theory ESE 524 Detetion and Estimation heory Joseh A. O Sullivan Samuel C. Sahs Professor Eletroni Systems and Signals Researh Laboratory Eletrial and Systems Engineering Washington University 2 Urbauer Hall

More information

A simple expression for radial distribution functions of pure fluids and mixtures

A simple expression for radial distribution functions of pure fluids and mixtures A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.

More information

Wave Propagation through Random Media

Wave Propagation through Random Media Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene

More information

Applications of Complex Analysis to Celestial Navigation

Applications of Complex Analysis to Celestial Navigation Aliations of Comlex Analysis to Celestial Navigation Robin G. Stuart Valhalla, New York, USA Abstrat Maing oints on the shere by stereograhi rojetion to oints on the lane of omlex numbers transforms the

More information

A Characterization of Wavelet Convergence in Sobolev Spaces

A Characterization of Wavelet Convergence in Sobolev Spaces A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics IC-W08D2-11 Jumping Off as Flatcar Solution

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics IC-W08D2-11 Jumping Off as Flatcar Solution N eole, eah o mass MASSACHUSETTS INSTITUTE OF TECHNOLOGY Deartment o Physis Physis 8.01 IC-W08D2-11 Juming O as Flatar Solution m, stand on a railway latar o mass m. They jum o one end o the latar with

More information

Advanced Computational Fluid Dynamics AA215A Lecture 4

Advanced Computational Fluid Dynamics AA215A Lecture 4 Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas

More information

arxiv: v2 [cs.ai] 16 Feb 2016

arxiv: v2 [cs.ai] 16 Feb 2016 Cells in Multidimensional Reurrent Neural Networks arxiv:42.2620v2 [s.ai] 6 Feb 206 Gundram Leifert Tobias Sauß Tobias Grüning Welf Wustlih Roger Labahn University of Rostok Institute of Mathematis 805

More information

Nuclear Shell Structure Evolution Theory

Nuclear Shell Structure Evolution Theory Nulear Shell Struture Evolution Theory Zhengda Wang (1) Xiaobin Wang () Xiaodong Zhang () Xiaohun Wang () (1) Institute of Modern physis Chinese Aademy of SienesLan Zhou P. R. China 70000 () Seagate Tehnology

More information

On the Quantum Theory of Radiation.

On the Quantum Theory of Radiation. Physikalishe Zeitshrift, Band 18, Seite 121-128 1917) On the Quantum Theory of Radiation. Albert Einstein The formal similarity between the hromati distribution urve for thermal radiation and the Maxwell

More information

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution. arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat

More information

Quantum Mechanics: Wheeler: Physics 6210

Quantum Mechanics: Wheeler: Physics 6210 Quantum Mehanis: Wheeler: Physis 60 Problems some modified from Sakurai, hapter. W. S..: The Pauli matries, σ i, are a triple of matries, σ, σ i = σ, σ, σ 3 given by σ = σ = σ 3 = i i Let stand for the

More information

Stress and Displacement Estimates for Arches

Stress and Displacement Estimates for Arches Stress and Dislaement Estimates for Arhes Clive L. Dym, F.ASCE 1 ; and Harry E. Williams Abstrat: This aer resents analytial estimates of the behavior exhibited by urved, arhlike strutures under radially

More information

A MULTIRATE APPROACH TO DDS ALGORITHM VERSUS TAYLOR-SERIES EXPANSION TECHNIQUE

A MULTIRATE APPROACH TO DDS ALGORITHM VERSUS TAYLOR-SERIES EXPANSION TECHNIQUE A ULTIRATE APPROACH TO DDS ALGORITH VERSUS TAYLOR-SERIES EXPANSION TECHNIQUE Ewa Hermanowiz 1 and iroslaw Rojewski Faulty of Eletronis, Teleommuniations and Informatis, Gdansk University of Tehnology,

More information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple

More information

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. All systems with interation of some type have normal modes. One may desribe them as solutions in absene of soures; they are exitations of the system

More information

Properties of Quarks

Properties of Quarks PHY04 Partile Physis 9 Dr C N Booth Properties of Quarks In the earlier part of this ourse, we have disussed three families of leptons but prinipally onentrated on one doublet of quarks, the u and d. We

More information

Analysis of discretization in the direct simulation Monte Carlo

Analysis of discretization in the direct simulation Monte Carlo PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,

More information