POLYRATE 4

A new version of a computer program for the calculation of chemical reaction rates for polyatomics

Da hong Lu, Thanh N. Truong, Vasilios S. Melissas, Gillian C. Lynch, Yi Ping Liu, Bruce C. Garrett, Rozeanne Steckler, Alan D. Isaacson, Sachchida N. Rai, Gene C. Hancock, J. G. Lauderdale, Tomi Joseph, Donald G. Truhlar

Research output: Contribution to journalArticle

491 Citations (Scopus)

Abstract

POLYRATE is a computer program for the calculation of chemical reaction rates of polyatomic species (and also atoms and diatoms as special cases). Version 1.1 was submitted to the CPC Program Library in 1987, and since that time we have considerably improved the program and made it more portable, and we have added several new capabilities, resulting in the present improved version 4. The methods used are variational or conventional transition state theory and multidimensional semiclassical adiabatic and large-curvature approximations for tunneling and nonclassical reflection. Rate constants may be calculated for canonical or microcanonical ensembles or for specific vibrational states of selected modes with translational, rotational, and other vibrational modes treated thermally. Bimolecular and unimolecular reactions and gas-phase, solid-state, and gas-solid interface reactions are all included. Potential energy surfaces may be global analytic functions or implicit functions defined by interpolation from input energies, gradients, and force constants (Hessian matrices) at selected points on a reaction path. The program calculates reaction paths by the Euler, Euler stabilization, or Page-McIver methods. Variational transition states are optimized from among a one-parameter sequence of generalized transition states orthogonal to the reaction path. Tunneling probabilities are calculated by numerical quadrature, using either the centrifugal-dominat-small-curvature approximation, the large-curvature-version-3 approximation, and / or methods that were available earlier. In the large curvature case the tunneling probabilities may be summed over final vibrational states for exoergic reactions or initial vibrational states for endoergic reactions.

Original languageEnglish (US)
Pages (from-to)235-262
Number of pages28
JournalComputer Physics Communications
Volume71
Issue number3
DOIs
StatePublished - 1992
Externally publishedYes

Fingerprint

Reaction rates
Computer program listings
Chemical reactions
chemical reactions
reaction kinetics
computer programs
Potential energy surfaces
Gases
Rate constants
curvature
Interpolation
vibrational states
Stabilization
Atoms
approximation
gas-solid interfaces
Hessian matrices
analytic functions
algae
quadratures

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

Lu, D. H., Truong, T. N., Melissas, V. S., Lynch, G. C., Liu, Y. P., Garrett, B. C., ... Truhlar, D. G. (1992). POLYRATE 4: A new version of a computer program for the calculation of chemical reaction rates for polyatomics. Computer Physics Communications, 71(3), 235-262. https://doi.org/10.1016/0010-4655(92)90012-N

POLYRATE 4 : A new version of a computer program for the calculation of chemical reaction rates for polyatomics. / Lu, Da hong; Truong, Thanh N.; Melissas, Vasilios S.; Lynch, Gillian C.; Liu, Yi Ping; Garrett, Bruce C.; Steckler, Rozeanne; Isaacson, Alan D.; Rai, Sachchida N.; Hancock, Gene C.; Lauderdale, J. G.; Joseph, Tomi; Truhlar, Donald G.

In: Computer Physics Communications, Vol. 71, No. 3, 1992, p. 235-262.

Research output: Contribution to journalArticle

Lu, DH, Truong, TN, Melissas, VS, Lynch, GC, Liu, YP, Garrett, BC, Steckler, R, Isaacson, AD, Rai, SN, Hancock, GC, Lauderdale, JG, Joseph, T & Truhlar, DG 1992, 'POLYRATE 4: A new version of a computer program for the calculation of chemical reaction rates for polyatomics', Computer Physics Communications, vol. 71, no. 3, pp. 235-262. https://doi.org/10.1016/0010-4655(92)90012-N
Lu, Da hong ; Truong, Thanh N. ; Melissas, Vasilios S. ; Lynch, Gillian C. ; Liu, Yi Ping ; Garrett, Bruce C. ; Steckler, Rozeanne ; Isaacson, Alan D. ; Rai, Sachchida N. ; Hancock, Gene C. ; Lauderdale, J. G. ; Joseph, Tomi ; Truhlar, Donald G. / POLYRATE 4 : A new version of a computer program for the calculation of chemical reaction rates for polyatomics. In: Computer Physics Communications. 1992 ; Vol. 71, No. 3. pp. 235-262.
@article{2dadbe79762b433c8fde18325f04a3ba,
title = "POLYRATE 4: A new version of a computer program for the calculation of chemical reaction rates for polyatomics",
abstract = "POLYRATE is a computer program for the calculation of chemical reaction rates of polyatomic species (and also atoms and diatoms as special cases). Version 1.1 was submitted to the CPC Program Library in 1987, and since that time we have considerably improved the program and made it more portable, and we have added several new capabilities, resulting in the present improved version 4. The methods used are variational or conventional transition state theory and multidimensional semiclassical adiabatic and large-curvature approximations for tunneling and nonclassical reflection. Rate constants may be calculated for canonical or microcanonical ensembles or for specific vibrational states of selected modes with translational, rotational, and other vibrational modes treated thermally. Bimolecular and unimolecular reactions and gas-phase, solid-state, and gas-solid interface reactions are all included. Potential energy surfaces may be global analytic functions or implicit functions defined by interpolation from input energies, gradients, and force constants (Hessian matrices) at selected points on a reaction path. The program calculates reaction paths by the Euler, Euler stabilization, or Page-McIver methods. Variational transition states are optimized from among a one-parameter sequence of generalized transition states orthogonal to the reaction path. Tunneling probabilities are calculated by numerical quadrature, using either the centrifugal-dominat-small-curvature approximation, the large-curvature-version-3 approximation, and / or methods that were available earlier. In the large curvature case the tunneling probabilities may be summed over final vibrational states for exoergic reactions or initial vibrational states for endoergic reactions.",
author = "Lu, {Da hong} and Truong, {Thanh N.} and Melissas, {Vasilios S.} and Lynch, {Gillian C.} and Liu, {Yi Ping} and Garrett, {Bruce C.} and Rozeanne Steckler and Isaacson, {Alan D.} and Rai, {Sachchida N.} and Hancock, {Gene C.} and Lauderdale, {J. G.} and Tomi Joseph and Truhlar, {Donald G.}",
year = "1992",
doi = "10.1016/0010-4655(92)90012-N",
language = "English (US)",
volume = "71",
pages = "235--262",
journal = "Computer Physics Communications",
issn = "0010-4655",
publisher = "Elsevier",
number = "3",

}

TY - JOUR

T1 - POLYRATE 4

T2 - A new version of a computer program for the calculation of chemical reaction rates for polyatomics

AU - Lu, Da hong

AU - Truong, Thanh N.

AU - Melissas, Vasilios S.

AU - Lynch, Gillian C.

AU - Liu, Yi Ping

AU - Garrett, Bruce C.

AU - Steckler, Rozeanne

AU - Isaacson, Alan D.

AU - Rai, Sachchida N.

AU - Hancock, Gene C.

AU - Lauderdale, J. G.

AU - Joseph, Tomi

AU - Truhlar, Donald G.

PY - 1992

Y1 - 1992

N2 - POLYRATE is a computer program for the calculation of chemical reaction rates of polyatomic species (and also atoms and diatoms as special cases). Version 1.1 was submitted to the CPC Program Library in 1987, and since that time we have considerably improved the program and made it more portable, and we have added several new capabilities, resulting in the present improved version 4. The methods used are variational or conventional transition state theory and multidimensional semiclassical adiabatic and large-curvature approximations for tunneling and nonclassical reflection. Rate constants may be calculated for canonical or microcanonical ensembles or for specific vibrational states of selected modes with translational, rotational, and other vibrational modes treated thermally. Bimolecular and unimolecular reactions and gas-phase, solid-state, and gas-solid interface reactions are all included. Potential energy surfaces may be global analytic functions or implicit functions defined by interpolation from input energies, gradients, and force constants (Hessian matrices) at selected points on a reaction path. The program calculates reaction paths by the Euler, Euler stabilization, or Page-McIver methods. Variational transition states are optimized from among a one-parameter sequence of generalized transition states orthogonal to the reaction path. Tunneling probabilities are calculated by numerical quadrature, using either the centrifugal-dominat-small-curvature approximation, the large-curvature-version-3 approximation, and / or methods that were available earlier. In the large curvature case the tunneling probabilities may be summed over final vibrational states for exoergic reactions or initial vibrational states for endoergic reactions.

AB - POLYRATE is a computer program for the calculation of chemical reaction rates of polyatomic species (and also atoms and diatoms as special cases). Version 1.1 was submitted to the CPC Program Library in 1987, and since that time we have considerably improved the program and made it more portable, and we have added several new capabilities, resulting in the present improved version 4. The methods used are variational or conventional transition state theory and multidimensional semiclassical adiabatic and large-curvature approximations for tunneling and nonclassical reflection. Rate constants may be calculated for canonical or microcanonical ensembles or for specific vibrational states of selected modes with translational, rotational, and other vibrational modes treated thermally. Bimolecular and unimolecular reactions and gas-phase, solid-state, and gas-solid interface reactions are all included. Potential energy surfaces may be global analytic functions or implicit functions defined by interpolation from input energies, gradients, and force constants (Hessian matrices) at selected points on a reaction path. The program calculates reaction paths by the Euler, Euler stabilization, or Page-McIver methods. Variational transition states are optimized from among a one-parameter sequence of generalized transition states orthogonal to the reaction path. Tunneling probabilities are calculated by numerical quadrature, using either the centrifugal-dominat-small-curvature approximation, the large-curvature-version-3 approximation, and / or methods that were available earlier. In the large curvature case the tunneling probabilities may be summed over final vibrational states for exoergic reactions or initial vibrational states for endoergic reactions.

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

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

U2 - 10.1016/0010-4655(92)90012-N

DO - 10.1016/0010-4655(92)90012-N

M3 - Article

VL - 71

SP - 235

EP - 262

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

IS - 3

ER -