Evaluation of thermal rates for reactions with intermediate wells

Removal of bound state contributions to quantum flux correlation functions

Mark Thachuk, Howard R. Mayne, George C Schatz

Research output: Contribution to journalArticle

13 Citations (Scopus)

Abstract

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.

Original languageEnglish
Pages (from-to)3516-3525
Number of pages10
JournalJournal of Chemical Physics
Volume99
Issue number5
Publication statusPublished - 1993

Fingerprint

Fluxes
evaluation
Wave packets
propagation
Chemisorption
wave packets
chemisorption
Rate constants
metal surfaces
Chemical reactions
chemical reactions
degrees of freedom
projection
Metals
Hot Temperature
continuums
operators
oscillations
matrices

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Evaluation of thermal rates for reactions with intermediate wells : Removal of bound state contributions to quantum flux correlation functions. / Thachuk, Mark; Mayne, Howard R.; Schatz, George C.

In: Journal of Chemical Physics, Vol. 99, No. 5, 1993, p. 3516-3525.

Research output: Contribution to journalArticle

@article{c6578f82f70a406a90f28cfbb1e05f96,
title = "Evaluation of thermal rates for reactions with intermediate wells: Removal of bound state contributions to quantum flux correlation functions",
abstract = "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.",
author = "Mark Thachuk and Mayne, {Howard R.} and Schatz, {George C}",
year = "1993",
language = "English",
volume = "99",
pages = "3516--3525",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "5",

}

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

VL - 99

SP - 3516

EP - 3525

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 5

ER -