Optimized calculations of reaction paths and reaction-path functions for chemical reactions
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431 (United States)
- Molecular Science Research Center, Pacific Northwest Laboratory, Richland, Washington 99352 (United States)
In this paper we optimize several algorithms for the computation of reaction rates based on information calculated along minimum energy reaction paths and we evaluate the efficiencies of the optimized algorithms. The investigations are based on the calculation of chemical reaction rate constants using variational transition state theory and multidimensional semiclassical transmission coefficients including reaction path curvature. Several methods are evaluated and compared by a systematic set of applications to test cases involving the hydrogen-atom transfer reactions CH{sub 3}+H{sub 2}{r arrow}CH{sub 4}+H and OH+H{sub 2}{r arrow}H{sub 2}O+H. For each method we present general recommendations for all algorithmic choices other than gradient step size so that future calculations may be carried out reasonably efficiently by varying only one parameter. In the process of these optimizations we have found that the accuracy of the Euler stabilization method can be significantly increased by choosing the auxiliary parameters differently than in previous work; the optimized algorithm is called ES1*. Our final recommendations for future work are (i) when the Hessian/gradient computational cost ratio is low ({approx lt}3): the Page--McIver algorithm with the Hessian recalculated at every step, with a cubic starting step, and with curvature calculated from the derivative of the gradient, and (ii) when the Hessian/gradient computational cost ratio is moderate or large: the ES1* algorithm with a Hessian step size three times larger than the gradient step size, with a quadratic starting step, and with curvature calculated from the derivative of the gradient.
- DOE Contract Number:
- FG02-86ER13579; AC06-76RL01830
- OSTI ID:
- 7281581
- Journal Information:
- Journal of Chemical Physics; (United States), Vol. 96:8; ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Quantum mechanics/molecular mechanics minimum free-energy path for accurate reaction energetics in solution and enzymes: Sequential sampling and optimization on the potential of mean force surface
The definition of reaction coordinates for reaction-path dynamics
Related Subjects
ORGANIC
PHYSICAL AND ANALYTICAL CHEMISTRY
CHEMICAL REACTIONS
ALGORITHMS
HYDROGEN
HYDROXYL RADICALS
METHYL RADICALS
ATOMS
EFFICIENCY
FUNCTIONS
METHANE
OPTIMIZATION
VARIATIONAL METHODS
WATER
ALKANES
ALKYL RADICALS
ELEMENTS
HYDROCARBONS
HYDROGEN COMPOUNDS
MATHEMATICAL LOGIC
NONMETALS
ORGANIC COMPOUNDS
OXYGEN COMPOUNDS
RADICALS
400201* - Chemical & Physicochemical Properties