POLYRATE-6.5 - A NEW VERSION OF A COMPUTER-PROGRAM FOR THE CALCULATION OF CHEMICAL-REACTION RATES FOR POLYATOMICS

Citation
R. Steckler et al., POLYRATE-6.5 - A NEW VERSION OF A COMPUTER-PROGRAM FOR THE CALCULATION OF CHEMICAL-REACTION RATES FOR POLYATOMICS, Computer physics communications, 88(2-3), 1995, pp. 341-343
Citations number
7
Categorie Soggetti
Mathematical Method, Physical Science","Physycs, Mathematical","Computer Science Interdisciplinary Applications
ISSN journal
00104655
Volume
88
Issue
2-3
Year of publication
1995
Pages
341 - 343
Database
ISI
SICI code
0010-4655(1995)88:2-3<341:P-ANVO>2.0.ZU;2-G
Abstract
POLYRATE is a computer program for the calculation of chemical reactio n rates of polyatomic species (and also atoms and diatoms as special c ases), Version 1.1 was submitted to the CPC Program Library in 1987, a nd version 4.0.1 was submitted in 1992, Since that time many new capab ilities have been added, old ones have been improved, and the code has been made more portable and user-friendly, resulting in the present i mproved version 6.5. The methods used are variational or conventional transition state theory and multidimensional semiclassical adiabatic a nd large-curvature approximations for tunneling and nonclassical refle ction, Rate constants may be calculated for canonical or microcanonica l ensembles or for specific vibrational states of selected modes with translational, rotational, and other vibrational modes treated thermal ly. 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 define d by interpolation from input energies, gradients, and force constants (Hessian matrices) at selected points on a reaction path. The data ne eded for the dynamics calculations may also be calculated from a globa l potential energy surface with more accurate calculations at stationa ry points. The program calculates reaction paths by the Euler, Euler s tabilization, or Page-McIver methods. Variational transition states ar e optimized from among a one-parameter sequence of generalized transit ion states orthogonal to the reaction path. Tunneling probabilities ar e calculated by numerical quadrature, using either the centrifugal-dom inant-small-curvature approximation, the large-curvature-version-3 app roximation, and/or optimized multidimensional tunneling approximations . In the large-curvature case the tunneling probabilities may be summe d over final vibrational states for exoergic reactions or initial vibr ational states for endoergic reactions.