TY - JOUR

T1 - Evaluation of thermal rates for reactions with intermediate wells

T2 - Removal of bound state contributions to quantum flux correlation functions

AU - Thachuk, Mark

AU - Mayne, Howard R.

AU - Schatz, George C

PY - 1993

Y1 - 1993

N2 - We consider the calculation of rates of chemical reactions that have bound intermediate states (i.e., wells along the reaction path) using flux correlation function methods. When time-dependent wave packets are used to evaluate the propagator matrix elements, and the dividing surface is located at a point where bound states have nonzero probability density, the standard expression for the flux correlation function shows infinitely long lived oscillations due to these bound states, making the evaluation of rate constants numerically ill-behaved. However, if the bound part of the initial wave packet is projected out, the resulting continuum-only propagators produce rapidly decaying correlation functions, and numerically well-behaved rate constants. We illustrate this projection operator approach by considering a one-dimensional reaction path model in which the potential is taken to be an Eckart well, and the dividing surface is located at the minimum. In another application, we consider a two degree of freedom model of H2 dissociative chemisorption on a rigid metal surface. This application is sufficiently complex that it is impractical to calculate the chemisorption rate using conventional flux correlation function methods, but with the projected wavepacket approach, the problem is made relatively, easy. We also consider a second approach to the treatment of bound states in which the flux correlation function is altered to remove implicitly the bound state contributions to the propagator at long times. This second expression can be used with the full propagator, eliminating the need to construct and project explicitly the bound states. This should be advantageous when many bound states are present.

AB - We consider the calculation of rates of chemical reactions that have bound intermediate states (i.e., wells along the reaction path) using flux correlation function methods. When time-dependent wave packets are used to evaluate the propagator matrix elements, and the dividing surface is located at a point where bound states have nonzero probability density, the standard expression for the flux correlation function shows infinitely long lived oscillations due to these bound states, making the evaluation of rate constants numerically ill-behaved. However, if the bound part of the initial wave packet is projected out, the resulting continuum-only propagators produce rapidly decaying correlation functions, and numerically well-behaved rate constants. We illustrate this projection operator approach by considering a one-dimensional reaction path model in which the potential is taken to be an Eckart well, and the dividing surface is located at the minimum. In another application, we consider a two degree of freedom model of H2 dissociative chemisorption on a rigid metal surface. This application is sufficiently complex that it is impractical to calculate the chemisorption rate using conventional flux correlation function methods, but with the projected wavepacket approach, the problem is made relatively, easy. We also consider a second approach to the treatment of bound states in which the flux correlation function is altered to remove implicitly the bound state contributions to the propagator at long times. This second expression can be used with the full propagator, eliminating the need to construct and project explicitly the bound states. This should be advantageous when many bound states are present.

UR - http://www.scopus.com/inward/record.url?scp=0040140956&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0040140956&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0040140956

VL - 99

SP - 3516

EP - 3525

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 5

ER -