1
|
Jasper AW, Harding LB, Knight C, Georgievskii Y. Anharmonic Rovibrational Partition Functions at High Temperatures: Tests of Reduced-Dimensional Models for Systems with up to Three Fluxional Modes. J Phys Chem A 2019; 123:6210-6228. [DOI: 10.1021/acs.jpca.9b03592] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
2
|
Jasper AW, Gruey ZB, Harding LB, Georgievskii Y, Klippenstein SJ, Wagner AF. Anharmonic Rovibrational Partition Functions for Fluxional Species at High Temperatures via Monte Carlo Phase Space Integrals. J Phys Chem A 2018; 122:1727-1740. [DOI: 10.1021/acs.jpca.7b11722] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ahren W. Jasper
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Zackery B. Gruey
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Lawrence B. Harding
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Yuri Georgievskii
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Stephen J. Klippenstein
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Albert F. Wagner
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
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
|
6
|
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]
|
7
|
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
|
8
|
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
|
9
|
Curotto E, Freeman DL, Doll JD. A stereographic projection path integral study of the coupling between the orientation and the bending degrees of freedom of water. J Chem Phys 2008; 128:204107. [DOI: 10.1063/1.2925681] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Avilés MW, McCandless ML, Curotto E. Stereographic projection path integral simulations of (HCl)n clusters (n=2–5): Evidence of quantum induced melting in small hydrogen bonded networks. J Chem Phys 2008; 128:124517. [DOI: 10.1063/1.2837802] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
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
|
12
|
Langley SF, Curotto E, Freeman DL, Doll JD. Rigid quantum Monte Carlo simulations of condensed molecular matter: Water clusters in the n=2→8 range. J Chem Phys 2007; 126:084506. [PMID: 17343457 DOI: 10.1063/1.2484229] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The numerical advantage of quantum Monte Carlo simulations of rigid bodies relative to the flexible simulations is investigated for some simple systems. The results show that if high frequency modes in molecular condensed matter are predominantly in the ground state, the convergence of path integral simulations becomes nonuniform. Rigid body quantum parallel tempering simulations are necessary to accurately capture thermodynamic phenomena in the temperature range where the dynamics are influenced by intermolecular degrees of freedom; the stereographic projection path integral adapted for quantum simulations of asymmetric tops is a significantly more efficient strategy compared with Cartesian coordinate simulations for molecular condensed matter under these conditions. The reweighted random series approach for stereographic path integral Monte Carlo is refined and implemented for the quantum simulation of water clusters treated as an assembly of rigid asymmetric tops.
Collapse
Affiliation(s)
- Stephen F Langley
- Department of Chemistry, Arcadia University, Glenside, Pennsylvania 19038-3295, 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
|
Patchkovskii S, Heine T. Evaluation of the adsorption free energy of light guest molecules in nanoporous host structures. Phys Chem Chem Phys 2007; 9:2697-705. [PMID: 17627313 DOI: 10.1039/b617657g] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present a simple method for estimation of the adsorption free energy of lightweight featureless molecules in nanoporous structures. The partition function of the system is computed quantum mechanically within the ideal gas approximation, and leads directly to the free adsorption energy. Storage capacities can be estimated by combining the ideal gas results with the real gas equation of state. Two possible implementations of the approach are discussed. The assumptions and potential sources of errors in the calculations are analyzed. The performance of the method is illustrated for graphitic slit pores and solid C60.
Collapse
Affiliation(s)
- Serguei Patchkovskii
- Steacie Institute for Molecular Sciences, NRC Canada, 100 Sussex Drive, Ottawa, Ontario, Canada K1A 0R6.
| | | |
Collapse
|
15
|
Avilés MW, Gray PT, Curotto E. Stereographic projection path-integral simulations of (HF)n clusters. J Chem Phys 2006; 124:174305. [PMID: 16689569 DOI: 10.1063/1.2192773] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We perform several quantum canonical ensemble simulations of (HF)(n) clusters. The HF stretches are rigid, and the stereographic projection path-integral method is employed for the simulation in the resulting curved configuration space. We make use of the reweighted random series techniques to accelerate the convergence of the path-integral simulation with respect to the number of path coefficients. We develop and test estimators for the total energy and heat capacity based on a finite difference approach for non-Euclidean spaces. The quantum effects at temperatures below 400 K are substantial for all sizes. We observe interesting thermodynamic behaviors in the quantum simulations of the octamer and the heptamer.
Collapse
Affiliation(s)
- Michael W Avilés
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA
| | | | | |
Collapse
|
16
|
Pav DM, Curotto E. Parallel tempering simulations of the 13-center Lennard-Jones dipole-dipole cluster (muD=0-->0.5 a.u.). J Chem Phys 2005; 123:144301. [PMID: 16238385 DOI: 10.1063/1.2049279] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the thermodynamic behavior of the thirteen center uniform Lennard-Jones dipole-dipole cluster [(LJDD)(13)] for a wide range of dipole moment strengths. We find a relatively wide range of potential parameters where solid-solid coexistence manifests itself. Using structural characterization methods we determine the shape of the few isomers that contribute to the solid-solid coexistence region. The thermal distributions of the size of the net dipole moment are broad even at the coldest temperatures of the simulation where the (LJDD)(13) cluster is solid.
Collapse
Affiliation(s)
- Dawn M Pav
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA
| | | |
Collapse
|
17
|
Curotto E. A reweighted random series method for stereographic projection path integrals. J Chem Phys 2005; 123:134102. [PMID: 16223270 DOI: 10.1063/1.2036970] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
A set of general reweighted random series methods for metric affine spaces is developed. The extension of the theorems to metric affine spaces demands the introduction of a configuration-independent reference metric tensor; this geometric object is used to treat the path expansion coefficients beyond the core path, in both the partial averaging and the reweighted random series approach. Numerical tests are conducted by simulating a particle in a ring. The reweighted random series results show better convergence properties and better statistical quality at a fraction of the cost compared with the related partial averaging simulation.
Collapse
Affiliation(s)
- E Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA.
| |
Collapse
|
18
|
Avilés MW, Curotto E. Partial averaging and the centroid virial estimator for stereographic projection path-integral simulations in curved spaces. J Chem Phys 2005; 122:164109. [PMID: 15945674 DOI: 10.1063/1.1884109] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop and test three different partial averaging theories for the stereographic projection path integral in curved spaces. Additionally, we adapt and test the centroid virial estimator for the kinetic energy in curved spaces. We tested both a confining as well as a nonconfining potential for the particle in a ring. All three partial averaging theories are demonstrated to converge linearly in the asymptotic region with k(-2)max, the number of Fourier coefficients. We use an error estimator to determine the optimal parameters for an extrapolation to infinite kmax. We verify that the centroid virial estimator (derived from the primitive DeWitt path-integral formula) converges to the kinetic energy for both confining and nonconfining potentials.
Collapse
Affiliation(s)
- Michael W Avilés
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA
| | | |
Collapse
|
19
|
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
|
20
|
Mielke SL, Truhlar DG. A ‘path-by-path’ monotone extrapolation sequence for Feynman path integral calculations of quantum mechanical free energies. Chem Phys Lett 2003. [DOI: 10.1016/j.cplett.2003.07.007] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
21
|
Miller TF, Clary DC. Torsional path integral Monte Carlo method for calculating the absolute quantum free energy of large molecules. J Chem Phys 2003. [DOI: 10.1063/1.1568727] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
|
22
|
Prudente FV, Varandas AJC. A Direct Evaluation of the Partition Function and Thermodynamic Data for Water at High Temperatures. J Phys Chem A 2002. [DOI: 10.1021/jp020797b] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
23
|
Miller TF, Clary DC. Torsional path integral Monte Carlo method for the quantum simulation of large molecules. J Chem Phys 2002. [DOI: 10.1063/1.1467342] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
24
|
Riganelli A, Prudente FV, Varandas AJC. On the Rovibrational Partition Function of Molecular Hydrogen at High Temperatures. J Phys Chem A 2001. [DOI: 10.1021/jp011330o] [Citation(s) in RCA: 20] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Antonio Riganelli
- Departamento de Química, Universidade de Coimbra, P-3049 Coimbra Codex, Portugal
| | | | | |
Collapse
|
25
|
Mielke SL, Truhlar DG. Displaced-points path integral method for including quantum effects in the Monte Carlo evaluation of free energies. J Chem Phys 2001. [DOI: 10.1063/1.1378043] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Mielke SL, Truhlar DG. A new Fourier path integral method, a more general scheme for extrapolation, and comparison of eight path integral methods for the quantum mechanical calculation of free energies. J Chem Phys 2001. [DOI: 10.1063/1.1290476] [Citation(s) in RCA: 43] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
27
|
Vidler M, Tennyson J. Accurate partition function and thermodynamic data for water. J Chem Phys 2000. [DOI: 10.1063/1.1321769] [Citation(s) in RCA: 103] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|