Description
Abstract
Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or of the VTST. In the current version, Pilgrim can evaluate these thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reactions by multi-path canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through the MSTor and Q2DTor programs. Dual-level calculations are also available: automatic high-level single point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate the whole process using KMC. This algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.
Funding information
Sponsorship: U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award DE-SC0015997; Consellería de Cultura, Educación e Ordenación Universitaria (Axuda para Consolidación e Estructuración de unidades de investigación competitivas do Sistema Universitario de Galicia, Xunta de Galicia ED431C 2017/17 & Centro singular de investigación de Galicia acreditación 2019-2022, ED431G 2019/03) and the European Regional Development Fund (ERDF)
Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or of the VTST. In the current version, Pilgrim can evaluate these thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reactions by multi-path canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through the MSTor and Q2DTor programs. Dual-level calculations are also available: automatic high-level single point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate the whole process using KMC. This algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.
Funding information
Sponsorship: U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award DE-SC0015997; Consellería de Cultura, Educación e Ordenación Universitaria (Axuda para Consolidación e Estructuración de unidades de investigación competitivas do Sistema Universitario de Galicia, Xunta de Galicia ED431C 2017/17 & Centro singular de investigación de Galicia acreditación 2019-2022, ED431G 2019/03) and the European Regional Development Fund (ERDF)
| Date made available | 2020 |
|---|---|
| Publisher | Data Repository for the University of Minnesota |