Example: Vibrational dephasing and energy relaxation next up previous
Next: About this document Up: No Title Previous: Mori-Zwanzig Theory: A more

Example: Vibrational dephasing and energy relaxation

Recall that the Fourier transform of a time correlation function can be related to some kind of frequency spectrum. For example, the Fourier transform of the velocity autocorrelation function of a particular degree of freedom q of interest

displaymath439

where tex2html_wrap_inline1267 , gives the relevant frequencies contributing to the dynamics of q, but does not give amplitudes. This ``frequency'' spectrum tex2html_wrap_inline1271 is simply given by

displaymath446

That is, we take the Laplace transform of tex2html_wrap_inline1273 using tex2html_wrap_inline1275 . Since tex2html_wrap_inline1273 carries information about the relevant frequencies of the system, the decay of tex2html_wrap_inline1273 in time is a measure of how strongly coupled the motion of q is to the rest of the bath, i.e., how much of an overlap there is between the relevant frequencies of the bath and those of q. The more of an overlap there is, the more mixing there will be between the system and the bath, and hence, the more rapidly the motion of the system will become vibrationally ``out of phase'' or decorrelated with itself. Thus, the decay time of tex2html_wrap_inline1273 , which is denoted tex2html_wrap_inline1287 is called the vibrational dephasing time.

Another measure of the strength of the coupling between the system and the bath is the time required for the system to dissipate energy into the bath when it is excited away from equilibrium. This time can be obtained by studying the decay of the energy autocorrelation function:

displaymath456

where tex2html_wrap_inline1289 is defined to be

displaymath459

The decay time of this correlation function is denoted tex2html_wrap_inline1291 .

The question then becomes: what are these characteristic decay times and how are they related? To answer this, we will take a phenomenological approach. We will assume the validity of the GLE for q:

displaymath463

and use it to calculate tex2html_wrap_inline1291 and tex2html_wrap_inline1287 .

Suppose the potential tex2html_wrap_inline1077 is harmonic and takes the form

displaymath468

Substituting into the GLE and dividing through by m gives

displaymath471

where

displaymath475

An equation of motion for tex2html_wrap_inline1273 can be obtained directly by multiplying both sides of the GLE by tex2html_wrap_inline1305 and averaging over a canonical ensemble:

displaymath481

Recall that

displaymath489

and note that

displaymath494

also

displaymath502

Thus,

displaymath509

Combining these results gives an equation for tex2html_wrap_inline1273

eqnarray514

which is known as the memory function equation and the kernel K(t) is known as the memory function or memory kernel. This type of integro-differential equation is called a Volterra equation and it can be solved by Laplace transforms.

Taking the Laplace transform of both sides gives

displaymath524

However, it is clear that tex2html_wrap_inline1311 and also

displaymath533

Thus, it follows that

eqnarray538

In order to perform the inverse Laplace transform, we need the poles of the integrand, which will be determined by the solutions of

displaymath549

which we could solve directly if we knew the explicit form of tex2html_wrap_inline1313 .

However, if tex2html_wrap_inline1315 is sufficiently larger than tex2html_wrap_inline1317 , then it is possible to develop a perturbation solution to this equation. Let us assume the solutions for s can be written as

displaymath554

Substituting in this ansatz gives

displaymath556

Since we are assuming tex2html_wrap_inline1321 is small, then to lowest order, we have

displaymath560

so that tex2html_wrap_inline1323 . The first order equation then becomes

displaymath562

or

displaymath565

Note, however, that

eqnarray569

Thus, stopping the first order result, the poles of the integrand occur at

displaymath575

Define

eqnarray579

Then

displaymath583

and tex2html_wrap_inline1273 is then given by the contour integral

displaymath589

Taking the residue at each pole, we find

displaymath594

which can be simplified to give

displaymath599

Thus, we see that the GLE predicts tex2html_wrap_inline1273 oscillates with a frequency tex2html_wrap_inline1329 and decays exponentially. From the exponential decay, we can directly read off the time tex2html_wrap_inline1287 :

displaymath605

That is, the value of the real part of the Fourier (Laplace) transform of the friction kernel evaluated at the renormalized frequency divided by 2m gives the vibrational dephasing time! By a similar scheme, one can easily show that the position autocorrelation function tex2html_wrap_inline1335 decays with the same dephasing time. It's explicit form is

displaymath611

The energy autocorrelation function tex2html_wrap_inline1337 can be expressed in terms of the more primitive correlation functions tex2html_wrap_inline1339 and tex2html_wrap_inline1273 . It is a straightforward, although extremely tedious, matter to show that the relation, valid for the harmonic potential of mean force, is

displaymath618

Substituting in the expressions for tex2html_wrap_inline1339 and tex2html_wrap_inline1273 gives

displaymath629

so that the decay time tex2html_wrap_inline1291 can be seen to be

displaymath633

and therefore, the relation between tex2html_wrap_inline1291 and tex2html_wrap_inline1287 can be seen immediately to be

displaymath637

The incredible fact is that this result is also true quantum mechanically. That is, by doing a simple, purely classical treatment of the problem, we obtained a result that turns out to be the correct quantum mechanical result!

Just how big are these times? If tex2html_wrap_inline1315 is very large compared to any typical frequency relevant to the bath, then the friction kernel evaluated at this frequency will be extremely small, giving rise to a long decay time. This result is expect, since, if tex2html_wrap_inline1315 is large compared to the bath, there are very few ways in which the system can dissipate energy into the bath. The situation changes dramatically, however, if a small amount of anharmonicity is added to the potential of mean force. The figure below illustrates the point for a harmonic diatomic molecule interacting with a Lennard-Jones bath. The top figure shows the velocity autocorrelation function for an oscillator whose frequency is approximately 3 times the characteristic frequency of the bath, while the bottom one shows the velocity autocorrelation function for the case that the frequency disparity is a factor of 6.

   figure643
Figure 1:


next up previous
Next: About this document Up: No Title Previous: Mori-Zwanzig Theory: A more

Mark Tuckerman
Thu Apr 27 14:44:56 EDT 2000