1
|
Szekeres B, Pártay LB, Mátyus E. Direct Computation of the Quantum Partition Function by Path-Integral Nested Sampling. J Chem Theory Comput 2018; 14:4353-4359. [PMID: 29944376 DOI: 10.1021/acs.jctc.8b00368] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the present work we introduce a computational approach to the absolute rovibrational quantum partition function using the path-integral formalism of quantum mechanics in combination with the nested sampling technique. The numerical applicability of path-integral nested sampling is demonstrated for small molecules of spectroscopic interest. The computational cost of the method is determined by the evaluation time of a point on the potential energy surface (PES). For efficient PES implementations, the path-integral nested sampling method can be a viable alternative to the direct-Boltzmann-summation technique of variationally computed rovibrational energies, especially for medium-sized molecules and at elevated temperatures.
Collapse
Affiliation(s)
- Béla Szekeres
- Institute of Chemistry , Eötvös Loránd University , Pázmány Péter sétány 1/A , Budapest H-1117 , Hungary.,Department of Numerical Analysis, Faculty of Informatics , Eötvös Loránd University , Pázmány Péter sétány 1/C , Budapest H-1117 , Hungary
| | - Lívia B Pártay
- Department of Chemistry , University of Reading , Whiteknights, Reading RG6 6AD , United Kingdom
| | - Edit Mátyus
- Institute of Chemistry , Eötvös Loránd University , Pázmány Péter sétány 1/A , Budapest H-1117 , Hungary
| |
Collapse
|
2
|
Mielke SL, Truhlar DG. A whole-path importance-sampling scheme for Feynman path integral calculations of absolute partition functions and free energies. J Chem Phys 2016; 144:034110. [PMID: 26801023 DOI: 10.1063/1.4939869] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Using Feynman path integrals, a molecular partition function can be written as a double integral with the inner integral involving all closed paths centered at a given molecular configuration, and the outer integral involving all possible molecular configurations. In previous work employing Monte Carlo methods to evaluate such partition functions, we presented schemes for importance sampling and stratification in the molecular configurations that constitute the path centroids, but we relied on free-particle paths for sampling the path integrals. At low temperatures, the path sampling is expensive because the paths can travel far from the centroid configuration. We now present a scheme for importance sampling of whole Feynman paths based on harmonic information from an instantaneous normal mode calculation at the centroid configuration, which we refer to as harmonically guided whole-path importance sampling (WPIS). We obtain paths conforming to our chosen importance function by rejection sampling from a distribution of free-particle paths. Sample calculations on CH4 demonstrate that at a temperature of 200 K, about 99.9% of the free-particle paths can be rejected without integration, and at 300 K, about 98% can be rejected. We also show that it is typically possible to reduce the overhead associated with the WPIS scheme by sampling the paths using a significantly lower-order path discretization than that which is needed to converge the partition function.
Collapse
Affiliation(s)
- Steven L Mielke
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis, Minnesota 55455-0431, USA
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
3
|
Major DT, Gao J. An Integrated Path Integral and Free-Energy Perturbation-Umbrella Sampling Method for Computing Kinetic Isotope Effects of Chemical Reactions in Solution and in Enzymes. J Chem Theory Comput 2015; 3:949-60. [PMID: 26627415 DOI: 10.1021/ct600371k] [Citation(s) in RCA: 79] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
An integrated centroid path integral and free-energy perturbation-umbrella sampling (PI-FEP/UM) method for computing kinetic isotope effects (KIEs) for chemical reactions in solution and in enzymes is presented. The method is based on the bisection sampling in centroid path integral simulations to include nuclear quantum effects to the classical potential of mean force. The required accuracy for computing kinetic isotope effects is achieved by coupled free-energy perturbation and umbrella sampling for reactions involving different isotopes. The use of FEP with respect to different masses results in relatively small statistical uncertainties, whereas if KIEs are computed directly by the difference in free energies obtained from the quantum mechanical potentials of mean force for different isotopes, the statistical errors are significantly greater. The PI-FEP/UM method is illustrated in two applications. The first reaction is the decarboxylation of N-methyl picolinate in water, for which the primary (13)C and secondary (15)N KIEs have been determined. The second reaction is the proton-transfer reaction between nitroethane and an acetate ion in water. In both cases, the computational results are in accord with experimental data, and the findings provide further insight into the mechanism of these reactions in water.
Collapse
Affiliation(s)
- Dan Thomas Major
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455, and Centro Nacional de Supercomputación, Programa Biología Computacional, C/ Jordi Girona 29, 08034 Barcelona, Spain
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455, and Centro Nacional de Supercomputación, Programa Biología Computacional, C/ Jordi Girona 29, 08034 Barcelona, Spain
| |
Collapse
|
4
|
Major DT, Garcia-Viloca M, Gao J. Path Integral Simulations of Proton Transfer Reactions in Aqueous Solution Using Combined QM/MM Potentials. J Chem Theory Comput 2015; 2:236-45. [PMID: 26626510 DOI: 10.1021/ct050257t] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
A bisection sampling method was implemented in path integral simulations of chemical reactions in solution in the framework of the quantized classical path approach. In the present study, we employ a combined quantum mechanical and molecular mechanical (QM/MM) potential to describe the potential energy surface and the path integral method to incorporate nuclear quantum effects. We examine the convergence of the bisection method for two proton-transfer reactions in aqueous solution at room temperature. The first reaction involves the symmetrical proton transfer between an ammonium ion and an ammonia molecule. The second reaction is the ionization of nitroethane by an acetate ion. To account for nuclear quantum mechanical corrections, it is sufficient to quantize the transferring light atom in the ammonium ion-ammonia reaction, while it is necessary to also quantize the donor and acceptor atoms in the nitroethane-acetate ion reaction. Kinetic isotope effects have been computed for isotopic substitution of the transferring proton by a deuteron in the nitroethane-acetate reaction. In all computations, it is important to employ a sufficient number of polymer beads along with a large number of configurations to achieve convergence in these simulations.
Collapse
Affiliation(s)
- Dan Thomas Major
- Department of Chemistry and Supercomputing Insititute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| | - Mireia Garcia-Viloca
- Department of Chemistry and Supercomputing Insititute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| | - Jiali Gao
- Department of Chemistry and Supercomputing Insititute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| |
Collapse
|
5
|
Mielke SL, Truhlar DG. Improved methods for Feynman path integral calculations and their application to calculate converged vibrational–rotational partition functions, free energies, enthalpies, entropies, and heat capacities for methane. J Chem Phys 2015; 142:044105. [PMID: 25637967 DOI: 10.1063/1.4905526] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Steven L. Mielke
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
6
|
Buchowiecki M, Vaníček J. Monte Carlo evaluation of the equilibrium isotope effects using the Takahashi–Imada factorization of the Feynman path integral. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.09.070] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
7
|
Mielke SL, Dinpajooh M, Siepmann JI, Truhlar DG. Efficient methods for including quantum effects in Monte Carlo calculations of large systems: Extension of the displaced points path integral method and other effective potential methods to calculate properties and distributions. J Chem Phys 2013; 138:014110. [PMID: 23298031 DOI: 10.1063/1.4772667] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Steven L Mielke
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, USA.
| | | | | | | |
Collapse
|
8
|
Mielke SL, Truhlar DG. Accelerating the Convergence and Reducing the Variance of Path Integral Calculations of Quantum Mechanical Free Energies by Using Local Reference Potentials. J Chem Theory Comput 2012; 8:1589-96. [DOI: 10.1021/ct300098p] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Steven L. Mielke
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis,
Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis,
Minnesota 55455-0431, United States
| |
Collapse
|
9
|
Xu X, Papajak E, Zheng J, Truhlar DG. Multi-structural variational transition state theory: kinetics of the 1,5-hydrogen shift isomerization of the 1-butoxyl radical including all structures and torsional anharmonicity. Phys Chem Chem Phys 2012; 14:4204-16. [DOI: 10.1039/c2cp23692c] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
10
|
Mielke SL, Truhlar DG. Improved Methods for Feynman Path Integral Calculations of Vibrational−Rotational Free Energies and Application to Isotopic Fractionation of Hydrated Chloride Ions. J Phys Chem A 2009; 113:4817-27. [DOI: 10.1021/jp900834u] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Steven L. Mielke
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
11
|
Anderson KE, Mielke SL, Siepmann JI, Truhlar DG. Bond Angle Distributions of Carbon Dioxide in the Gas, Supercritical, and Solid Phases. J Phys Chem A 2009; 113:2053-9. [DOI: 10.1021/jp808711y] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Kelly E. Anderson
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431
| | - Steven L. Mielke
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431
| | - J. Ilja Siepmann
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
12
|
Lynch VA, Mielke SL, Truhlar DG. High-precision quantum thermochemistry on nonquasiharmonic potentials: converged path-integral free energies and a systematically convergent family of generalized Pitzer-Gwinn approximations. J Phys Chem A 2007; 109:10092-9. [PMID: 16838929 DOI: 10.1021/jp051742n] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accurate quantum mechanical (QM) vibrational-rotational partition functions for HOOD, D(2)O(2), H(18)OOH, H(2)(18)O(2), D(18)OOH, and H(18)OOD are determined using a realistic potential energy surface for temperatures ranging from 300 to 2400 K by using the TT-FPI-ESPE path-integral Monte Carlo method. These data, together with our prior results for H(2)O(2), provide benchmarks for testing approximate methods of estimating isotope effects for systems with torsional motions. Harmonic approximations yield poor accuracy for these systems, and although the well-known Pitzer-Gwinn (PG) approximation provides better results for absolute partition functions, it yields the same results as the harmonic approximation for isotope effects because these are intrinsically quantal phenomena. We present QM generalizations of the PG approximation that can provide high accuracy for both isotope effects and absolute partition functions. These approximations can be systematically improved until they approach the accurate result and converge rapidly. These methods can also be used to obtain affordable estimates of zero-point energies from accurate partition functions-even those at relatively high temperatures.
Collapse
Affiliation(s)
- Vanessa Audette Lynch
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, USA
| | | | | |
Collapse
|
13
|
Fernandez-Ramos A, Miller JA, Klippenstein SJ, Truhlar DG. Modeling the kinetics of bimolecular reactions. Chem Rev 2007; 106:4518-84. [PMID: 17091928 DOI: 10.1021/cr050205w] [Citation(s) in RCA: 393] [Impact Index Per Article: 23.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Antonio Fernandez-Ramos
- Departamento de Quimica Fisica, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | | | | | | |
Collapse
|
14
|
Pu J, Gao J, Truhlar DG. Multidimensional tunneling, recrossing, and the transmission coefficient for enzymatic reactions. Chem Rev 2006; 106:3140-69. [PMID: 16895322 PMCID: PMC4478620 DOI: 10.1021/cr050308e] [Citation(s) in RCA: 288] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jingzhi Pu
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431
| | - Jiali Gao
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431
| |
Collapse
|
15
|
Chempath S, Predescu C, Bell AT. Quantum mechanical single molecule partition function from path integral Monte Carlo simulations. J Chem Phys 2006; 124:234101. [PMID: 16821901 DOI: 10.1063/1.2196885] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
An algorithm for calculating the partition function of a molecule with the path integral Monte Carlo method is presented. Staged thermodynamic perturbation with respect to a reference harmonic potential is utilized to evaluate the ratio of partition functions. Parallel tempering and a new Monte Carlo estimator for the ratio of partition functions are implemented here to achieve well converged simulations that give an accuracy of 0.04 kcal/mol in the reported free energies. The method is applied to various test systems, including a catalytic system composed of 18 atoms. Absolute free energies calculated by this method lead to corrections as large as 2.6 kcal/mol at 300 K for some of the examples presented.
Collapse
Affiliation(s)
- Shaji Chempath
- Department of Chemical Engineering, University of California-Berkeley, Berkeley, CA 94720, USA
| | | | | |
Collapse
|
16
|
|
17
|
Abstract
A bisection sampling method is implemented in the framework of a quantized classical path algorithm to include nuclear quantum effects in path integral simulations. The present study examines the convergence of these calculations on two model systems with respect to the number of beads used in the polymer chain and the number of configurations both in free-particle sampling and in classical configuration sampling. The results will be useful for future studies of kinetic isotope effects in enzymatic reactions.
Collapse
Affiliation(s)
- Dan Thomas Major
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA.
| | | |
Collapse
|
18
|
Lynch VA, Mielke SL, Truhlar DG. Accurate vibrational-rotational partition functions and standard-state free energy values for H2O2 from Monte Carlo path-integral calculations. J Chem Phys 2004; 121:5148-62. [PMID: 15352807 DOI: 10.1063/1.1782511] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate quantum mechanical partition functions and absolute free energies of H(2)O(2) are determined using a realistic potential energy surface [J. Koput, S. Carter, and N. C. Handy, J. Phys. Chem. A 102, 6325 (1998)] for temperatures ranging from 300 to 2,400 K by using Monte Carlo path integral calculations with new, efficient polyatomic importance sampling methods. The path centroids are sampled in Jacobi coordinates via a set of independent ziggurat schemes. The calculations employed enhanced-same-path extrapolation of trapezoidal Trotter Fourier path integrals, and the paths were constructed using fast Fourier sine transforms. Importance sampling was also used in Fourier coefficient space, and adaptively optimized stratified sampling was used in configuration space. The free energy values obtained from the path-integral calculations are compared to separable-mode approximations, to the Pitzer-Gwinn approximation, and to values in thermodynamic tables. Our calculations support the recently proposed revisions to the JANAF tables.
Collapse
Affiliation(s)
- Vanessa Audette Lynch
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | | | | |
Collapse
|