Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene

Accelerating equilibrium isotope effect calculations. II. Stochastic implementation of direct estimators

Path integral calculations of equilibrium isotope effects and isotopic fractionation are expensive due to the presence of path integral discretization errors, statistical errors, and thermodynamic integration errors. Whereas the discretization errors can be reduced by high-order factorization of the path integral and statistical errors by using centroid virial estimators, two recent papers proposed alternative ways to completely remove the thermodynamic integration errors: Cheng and Ceriotti [J. Chem. Phys. 141, 244112 (2015)] employed a variant of free-energy perturbation called "direct estimators," while Karandashev and Vaníček [J. Chem. Phys. 143, 194104 (2017)] combined the thermodynamic integration with a stochastic change of mass and piecewise-linear umbrella biasing potential. Here, we combine the former approach with the stochastic change in mass in order to decrease its statistical errors when applied to larger isotope effects and perform a thorough comparison of different methods by computing isotope effects first on a harmonic model and then on methane and methanium, where we evaluate all isotope effects of the form CH_4-xD_x/CH₄ and CH_5-xD_x ⁺/CH₅ ⁺, respectively. We discuss the reasons for a surprising behavior of the original method of direct estimators, which performed well for a much larger range of isotope effects than what had been expected previously, as well as some implications of our work for the more general problem of free energy difference calculations.

Can Path Integral Molecular Dynamics Make a Good Approximation for Vapor Pressure Isotope Effects Prediction for Organic Solvents? A Comparison to ONIOM QM/MM and QM Cluster Calculation

Isotopic fractionation of volatile organic compounds (VOCs), which are under strict measures of control because of their potential harm to the environment and humans, has an important ecological aspect, as the isotopic composition of compounds may depend on the conditions in which such compounds are distributed in Nature. Therefore, detailed knowledge on isotopic fractionation, not only experimental but also based on theoretical models, is crucial to follow conditions and pathways within which these contaminants are spread throughout the ecosystems. In this work, we present carbon and, for the first time, bromine vapor pressure isotope effect (VPIE) on the evaporation process from pure-phase systems-dibromomethane and bromobenzene, the representatives of aliphatic and aromatic brominated VOCs. We combine isotope effects measurements with their theoretical prediction using three computational techniques, namely path integral molecular dynamics, QM cluster, and hybrid ONIOM models. While evaporation of both compounds resulted in normal bromine VPIEs, the difference in the direction of carbon isotopic fractionation is observed for the aliphatic and aromatic compounds, where VPIEs are inverse and normal, respectively. Even though theoretical models tested here turned out to be insufficient for quantitative agreement with the experimental values, cluster electronic structure calculations, as well as two-layer ONIOM computations, provided better reproduction of experimental trends.

The role of tachysterol in vitamin D photosynthesis – a non-adiabatic molecular dynamics study

To investigate the role of tachysterol in the regulation of vitamin D photosynthesis, we studied its absorption properties and photodynamics by ab initio methods and non-adiabatic molecular dynamics.

Quantum Instanton Evaluations of the Thermal Rate Constants for Complex Systems

Quantum instanton (QI) approximation is recently proposed for the evaluations of the chemical reaction rate constants with use of full dimensional potential energy surfaces. Its strategy is to use the instanton mechanism and to approximate time-dependent quantum dynamics to the imaginary time propagation of the quantities of partition function. It thus incorporates the properties of the instanton idea and the quantum effect of partition function and can be applied to chemical reactions of complex systems. In this paper, we present the QI approach and its applications to several complex systems mainly done by us. The concrete systems include, (1) the reaction of H+CH4→H2+CH3, (2) the reaction of H+SiH4→H2+SiH3, (3) H diffusion on Ni(100) surface; and (4) surface-subsurface transport and interior migration for H/Ni. Available experimental and other theoretical data are also presented for the purpose of comparison.

THE RATE CONSTANT OF THE O(3P) + HCl → OH + Cl REACTION WITHIN QUANTUM INSTANTON APPROXIMATION

Temperature dependence of the rate constant of the O (3 P ) + HCl → OH + Cl reaction was calculated for 200–700 K range within the quantum instanton approximation using the 3A″ potential energy surface. Curvature of the Arrhenius plot arising from quantum effects in the hydrogen exchange reaction was observed. Results of the present method were found to significantly improve upon those obtained within the ICVT/μOMT method. The errors pertaining to the present method were estimated.