arxiv:cond-mat/ v4 [cond-mat.stat-mech] 14 Jan 2007

Transcription

1 Exact analytical evaluation of time dependent transmission coefficient from the method of reactive flux for an inverted parabolic barrier arxiv:cond-mat/79v4 [cond-mat.stat-mech] 4 Jan 7 Raarshi Chakrabarti Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore, 56, India In this paper we derive a general expression for the transmission coefficient using the method of reactive flux for a particle coupled to a harmonic bath surmounting a one dimensional inverted parabolic barrier. Unlike Kohen and Tannor [J. Chem. Phys. 3, ] we use a normal mode analysis where the unstable and the other modes have a complete physical meaning. Importantly our approach results a very general expression for the time dependent transmission coefficient not restricted to overdamped limit. Once the spectral density for the problem is known one can use our formula to evaluate the time dependent transmission coefficient. We have done the calculations with time dependent friction used by Xie [Phys. Rev. Lett 93, 863 4] and also the one used by Kohen and Tannor [J. Chem. Phys. 3, ]. Like the formula of Kohen and Tannor our formula also reproduces the results of transition state theory as well as the Kramers theory in the limits t and t respectively. I. INTRODUCTION The dynamics of a particle surmounting a barrier is an important and interesting problem in chemistry and physics. In chemistry, chemical reactions are most common examples where one studies the dynamics of barrier crossing along a suitably defined reaction coordinate. While in physics nucleation phenomena, electrical transport etc. involve barrier crossing. A review by Hänggi, Talkner and Borkovec [] gives a comprehensive account on the subect of barrier crossing. The simplest possible approach to calculate the rate constant for such processes is given by transition state theory [, ]. As the simplest version of the problem one can describe the rate process along a single reaction co-ordinate. Then the one dimensional version of the transition state rate constant k TST reads as, k TST = ω π exp[ E b k ], where ω BT is the angular frequency in the reactant well and E b is the barrier height or the activation energy. The derivation of transition state rate constant is based on two assumptions. Firstly an overall thermal equilibrium is assumed and secondly frictional effect is not taken into account by incorporating what is known as no recrossing assumption. Subsequently, Kramers [,, 3, 4] theory takes into account of the friction in a phenomenological way. His formulation was based on Markovian dynamics of a Brownian particle escaping from a metastable state. Further extensions of Kramers s theory were carried out by Grote and Hynes [5] followed by Hänggi and Motabi [6]. They took into account of the non-markovian nature of the dynamics. Another approach to calculate the rate constant is due to Chandler [7, 8] known as the method of reactive flux in chemical physics literature. In this approach the thermal rate constant is written as a correlation function which is nothing but the ensemble average over infinite number of traectories starting at the barrier top and ending on the product side at time t. The rate constant at time t is written as a product of the time dependent transmission coefficient κt and the transition state rate constant k TST. Thus one writes, kt = κtk TST. In the limit t, κt and one obtains k = k TST. Here wefollowthemethodofreactivefluxtocalculatethetime dependent rate constant for a particle bilinearly coupled to a harmonic bath [9,, ]. The bilinear coupling allows us to get an analytical expression for the time dependent transmission coefficient. Although the method ofreactivefluxandthe hamiltonianwehaveusedarewell known in the literature, as far as our knowledge nobody has used such hamiltonian to calculate the time dependent transmission coefficient analytically. The calculation is of interest as the transmission coefficient for the overdamped limit of the model was obtained recently by simulation []. Our approach directly gives expression for the time dependent transmission coefficient in terms of the spectral density Jω and not restricted to overdamped limit. A normal mode analysis of the coupled Hamiltonian was made by Pollak [3, 4, 5] and others to analyze the Kramers problem. Here we point out that the approximation is very powerful and can be used to get the exact expression for κt. Further we believe that it can also be used to get the quantum expression for κt. An exact calculation of κt using the phase space distribution function formulation was carried out in classical regime by Kohen and Tannor [6, 7], Bao [8], and more recently in the context of single enzyme kinetics by Chowdhury and Cherayil[9]. Also the calculation of the transmission coefficient in the quantum domain is carried out by Ray s group []. Our analysis in terms of normal modes has the advantage that it makes the physical ideas very clear. The paper is arranged as follows. In section II we introduce our system plus harmonic bath hamiltonian and the Generalized Langevin Equation GLE. In section III we mention the method of reactive flux and using our hamiltonian derive an analytical expression for the time dependent transmission coefficient, κt. Section IV deals with the calculation of κt for Markovian and Non Markovian dynamics. Section V is conclusion.

2 II. GENERALIZED LANGEVIN EQUATION AND NORMAL MODE ANALYSIS We consider a particle coupled to a harmonic bath [9, ] with the total hamiltonian H = P M +VQ+ P +m ω Q c Q m m ω Here P,Q are the momenta and coordinates of the th bath oscillator whose mass and frequency are m,ω respectively. P, Q is the momentum and coordinate of the system. c couples the bath oscillator to the system. Now in terms of the mass weighted coordinates our Hamiltonian becomes H = p +Vq+ q = MQ,q = m Q p c +ω q q Mm 3 and then with this Hamiltonian one can derive what is know as the Generalized Langevin Equation GLE which reads as t qt+v q+ dt γt t qt = ζt 4 One can also write the time dependent friction in terms of the spectral density of the bath Jω as follows where γt = Θt M Jω = π π One can easily show that dω Jω ω cosωt 5 c m ω δω ω 6 ζt ζt = Mk B Tγt t 7 We take the barrier top to be at q = and introduce ωb = V q. For small amplitude motion around the q= barriertop, one can use the normal modes, η = N U i q i i= where U is an orthogonal matrix such that U D U T is diagonal. D is a dynamical matrix given by ωb + N c Mm ω.. c Mm D = c ω Mm The reaction co-ordinate q may then be written as q = U η i = 8 Of the modes η, we take η to be the unstable mode. Wedenotethefrequencyofη asλ. η hastheimaginary frequency and we write it as λ = iλ, where Λ is real. As η is a normal mode, its time development is given by η t = η cosλ t+p η sinλ t λ with =,,,...N. Now the Hamiltonian when written in terms of the normal modes the coupling is no more there and is given by H nor = = p + λ η 9 In the next section we analytically derive the rate constant in the normal mode description using the above Hamiltonian. III. METHOD OF REACTIVE FLUX AND THE RATE CONSTANT In this section we briefly describe the method of reactiveflux [7, 8] tocalculatethe rateconstant. The method of reactive flux expresses the thermal rate constantwritten in the mass-weighted coordinates as kt = qtδqθ Pqt θ R q = k nt k d t where the angular bracket indicates the thermal average over the system as well as the bath degrees of freedom. The top of the barrier is chosen as q =. θ R q is if q < and otherwise. θ p is ust -θ R, q is the velocity of the particle. The numerator of Eq. could be written as k n t = i r expir qδqθ P qt r =

3 3 Expressing the step function as an integral over delta function and using the Fourier integral representation of delta function the expression for k n t becomes k n t = π dh i r ds ds exp is h expir q+is qt+is q r = The angular bracket indicates the thermal average over all degrees of freedom of the problem and it is done by multiplying the quantities within the angular bracket by exp βh nor, and carrying out the integration over all coordinates and momenta degrees of freedom. Writing k n t in terms of η leads to k n t = exp βe b π dh i r ds expir U η expis dp p ds exp is hexp β η Λ η expir U p η U η coshλt+p η sinhλt Λ expir U η expir U p η p dη exp β + λ η expis U η cosλ t+p η sinλt λ r = where E b is the barrier height. After one carries out the integrations the expression for k n t can be written as k n t = π βλ ċt exp βe βλ ct c b 3 Similarly the denominator of the Eq. can be writ- ten as k d t = = dp p dη exp β + λ η = π = βλ 4 Then the time dependent rate constant becomes kt = k nt k d t = πλ λ = λ ċt ct c exp βe b 5

4 Here we use the following identity [], and this simplifies our kt NQ λ = = ω Q Λ N ω b, λ In the next section we calculate ct and from that κt for Markovian as well as for Non Markovian noise. 4 kt = ω ċt π ω b ct c where, ct = defined by = U cosλt λ exp βe b 6. As usual kt [6] is kt = κtk 7 We show in the Appendix that k = ω π exp βe b, which is nothing but the transition state rate constant k TST. In the Appendix we also show how one gets Kramers rate constant in the t limit. Using the above expression for k our time dependent transmission coefficient, κt becomes κt = ċt 8 ω b ct c IV. CALCULATION OF ct AND κt Calculation of κt involves the calculation of ct. First we describe a general derivation of ct then we evaluate κt for different cases. A. Calculation of ct ct is defined by ct = N = by writing N U δω λ = π Im N evaluate N U ω +iη λ U cosλt. One proceeds λ U ω +iη λ. To here we adopt notations similar to quantum mechanics and also use the fact that D is a matrix with eigen values λ, =,,...N. U ω +iη λ = ω +iη λ = ω 9 +iη D Then by introducing in between the sum and using the resolution of identity N = I. U δω λ = ω +iη D = π Im Gω +iη Using partitioning technique one can evaluate G to get Gω +iη = ω +ω b M c ω m ω ω +iη ω Then using the definition of Jω Eq.6 it becomes Gω +iη = ω +ωb ω Mπ J ωd ω ωω +iη ω As π Im Gω = N U δω λ = ρω. We can write ct = = dω ω cosωtρω. B. Calculation of κt with γt = γh H t H First we perform the calculation with the time dependent friction γt used by Xie [], which is

5 5 Transmission coefficient t t. FIG. : Time dependent transmission coefficient, κt against time, for H = dashed line and H = 3 solid line, parameters used are ω b =,γ = 4. a inverse Fourier transform of Eq. 5 and write Jω in terms of γt. Jω = Mω Using the definition of γt we get dtγt cosωt 4 ω J ωd ω Mπ ωω +iη ω = γγh +Λ H First we consider the case ω = Λ which means we consider the unstable mode of the problem. Then Transmission Coefficient t t FIG. : Time dependent transmission coefficient, κt against time for Non-Markovian friction, γt = γh exp ht. Three different regimes, caging regime solid line with parameters, γ = 5,h =.,ω b =, intermediate regime dotted line with parameters γ = 5,h =,ω b =, and nonadiabatic regime dashed line with parameters γ = 5,h =.,ω b =. γt = γhh t H 3 G Λ = Λ +ω b γγh +Λ H 5 Our next task is to find out the poles, i.e. the unstable modes for which G Λ blows up, i.e. the solutions of the equation Λ +ω b γγh +Λ H = 6 On the other hand, to find ρω we have to consider a situation where ω is positive. In other words the stable modes. Now the Evaluation of ρω involves the integral ω J ωd ω Mπ ωω +iη ω. Using the identity ω +iη ω = P ω ω iπδ ω ω, wherep ω ω istheprincipal value of ω ω we find ω J ωd ω Mπ ωω +iη ω = expihπγh+ω H γ 7 Thus Gω = ω +ω b +expihπγh +ω H γ with, < H <, thus the thermal fluctuations are described by fractional Gaussian noise []. Now we take and 8 ρω = π Im Gω = γγh+ω H sinhπ π ω +ωb coshπ +γγh+ω H +γ sin HπΓH+ ω 4 4H 9

6 6. The Case H = First we consider H =, for which the γt = γδt and with this the noise-noise correlation function defined in Eq. 7 becomes delta function correlated, reading ζt ζt = Mγk B Tδt t 3 Then the dynamics is described by usual Langevin equation rather than a GLE. Thus the dynamics is then Markovian i.e. with no memory. For the Markovian case Eq. 6 becomes, Λ + ωb γλ = and has two roots. But only one of them is physically acceptable. The root is, Λ H= = γ + γ +ω b. Similarly from Eq. 9 one gets ρω H= ω γ +ω +ω b. = ωγ π Then Λ H= and ρω H= are used to calculate ct and taking t limit one gets c. In this case we could do these analytically. The time dependent transmission coefficient in this case is given by κt = +e tf γ 4e ft +e tf γ +e tf γ γ ω b f fe tf γ +γ +e ft where f = γ +4ω b +γ /. κt is then calculated using the above expression and plotted against time Fig... The Case H = 3 4 When H then the noise is no longer delta function correlated, i.e. not a white noise with Gaussian distribution. This regime is know as Non Markovian regime which has memory. We have similarly evaluated κt for H = 3 4 and plotted against time in Fig.. But in this case the integrations are done numerically since analytical evaluation of integrals were not possible. C. Calculation of κt with γt = γhexp ht In this section we carryout the calculation of κt with a time dependent friction γt = γhexp ht. The same friction is used by Kohen and Tannor [6, 7]. First we calculate Jω with this friction as follows Jω = Mγhω exp htcosωtdt = γmh ω h +ω Then the following integral is evaluated as follows ω J ωd ω Mπ α + ω ω +iη ω = ω h π d ω α + ω ω +iη ω = γhωω ih h +ω We have used the identity, ω +iη ω = P ω ω iπδ ω ω, where P ω ω is the principal value of ω ω. Thus and Gω = ω +ω b γhωω ih h +ω ρω = π Im Gω = ωγh h +ω π h +ω h +ω b γhω +ω γ h 4

7 7 One should notice that lim h ωγh h +ω = π h +ω h +ω b γhω +ω γ h 4 γω πω +ω b +γ ω = ρω H= Thus in this limit the result obtained is identical with the friction γt = γhh t H with H =. In other words the Markovian dynamics is recovered. To calculate κt we proceed as in the previous case. First considering the case ω = Λ, one gets Thus ω J ωd ω Mπ α + ω ω +iη ω = γhλ h+λ Gω Λ = +ω b γhλ h+λ Hencetofindoutthepolesonehastosolvethesolution of the equation method of reactive flux [7, 8] to calculate the time dependent transmission coefficient, κt. Analytically we derive a general formula for κt Eq. 8. For the Markovian case with no memory H = for the friction γt = γhh t H with H = which is same as with h = for the friction γt = γhexp ht, κt is calculated analytically. Whereas in the case with memory H = 3 4 and with any finite value of h the calculations were done numerically. In all the cases a plot of κt vs t starts from at t = and reaches a plateau in the long time limit. As expected our formula for the time dependent rate constant, kt becomes equal to the transition state rate constant when one takes the limit t [6, 7]. Similarly Kramers rate constant is obtained by taking t [6, 7]. In future we would like to extend our formulation to quantum domain. Λ +ω b γhλ h+λ = 3 To explore the non-markovian dynamics we choose three set of parameters corresponding to three different regimes of non-markovian dynamics following Kohen and Tannor [6, 7]. These regimes are called non-adiabatic, caging and intermediate according to Kohen and Tannor [6, 7]. The time dependent coefficients in these regimes are calculated numerically and plotted against time. We observe similar kind of oscillatory behavior of κt in the caging regime as obtained by them. Similarly in the intermediate regime the plot looks similar with a nonmonotonic decay Fig.. VI. ACKNOWLEDGEMENTS The author is grateful to Prof. K. L. Sebastian for his useful comments and encouragements and thanks Prof. B. J. Cherayil for his comments on the manuscript. The author also acknowledges Council of Scientific and Industrial Research CSIR, India for financial support. V. CONCLUSIONS VII. APPENDIX : TRANSITION STATE RATE AND KRAMERS RATE FROM kt The paper shows an elegant way of combining the traditional system plus reservoir model [9, ] and the First we take the limit t. Then from Eq. 6 one writes

8 8 k = lim t k+t = ω π ω b lim t ċ+t c+t c exp βe b = ω π ω b lim t ċ ct c c+tċ+ t c «exp βe b = ω π ω b lim t ct c c +c ct exp βe b = ω c π ω b exp βe b = ω c c π ω b exp βe b c Where we have used the facts, ċ = and c =. This is because ċt = dωρω sinωt making ċ =. Similarly ct = dωρω sinωtω = U = N as follows =. Now one can also show, c = ω b c = N = U λ = N = U λ U = N = λ = lim ω ω D = lim ω = = ω D = lim ω Gω From Eq. one can see that lim ω Gω = = c. Then k is nothing but the transition ωb state rate constant, k TST. Similarly taking t one gets the Kramers rate constant as follows k = ω π exp βe b = k TST 3 k = ω π ω b U Λ sinhλt s U Λ coshλt «t exp βe b = ω π ω b ΛtanhΛt t exp βe b Then k = Λ ω exp βe b π ω b γ Now if one replaces Λ by Λ H= = γ + +ω b H = means white noise one gets k = γ γ + +ω ω b π ω b exp βe b Which is nothing but the usual Kramers rate. Also in the overdamped limit i.e. when γ >> ωb one gets

9 9 k = γ + γ + ω 4 b ω γ exp βe b = ω ω b πω b πγ exp βe b which is Kramers rate expression in overdamped regime. [] P. Hänggi, P. Talkner and M. Borkovec, Rev. Mod. Phys. 6, [] A. Nitzan, Chemical Dynamics in condensed phases Oxford University Press, New York, 6. [3] H. A. Kramers, Physica Amsterdam 7, [4] W. T. Coffey, Y. P. Kalmykov, and J. T. Waldron, The Langevin Equation, Second Edition World Scientific, 4. [5] R. F. Grote and J. T. Hynes, J. Chem. Phys. 73, [6] P. Hänggi and F. Motabai, Phys. Rev. A 6, [7] D. Chandler, J. Stat. Phys. 4, [8] D. Chandler, Introduction to Modern Statistical Mechanics Oxford University Press, New York, 987. [9] U. Weiss, Quantum Dissipative Systems World Scientific, 999. [] R. Zwanzig, Nonequilibrium Statistical Mechanics Oxford University Press, New York,. [] H. Zhou and R. Zwanzig, J. Phys. Chem. A 6, 756. [] W. Min and X. S. Xie, Phys. Rev. E 73, 9 6. [3] A. M. Levin, M. Shapiro and E. Pollak, J. Chem. Phys. 88, [4] E. Pollak, J. Chem. Phys. 85, [5] E. Pollak, H. Grabert, and P. Hänggi, J. Chem. Phys. 9, [6] D. Kohen and D. J. Tannor, J. Chem. Phys. 3, [7] D. Kohen and D. J. Tannor, Adv. Chem. Phys., 9. [8] J. Bao, J. Chem. Phys. 4, [9] S. Chowdhury and B. J. Cherayil, J. Chem. Phys. 5, [] D. Barik, S K. Banik and D. S. Ray, J. Chem. Phys. 9, [] S. C. Kou and X. S. Xie, Phys. Rev. Lett. 93,