1
|
Yu Q, Qu C, Houston PL, Conte R, Nandi A, Bowman JM. q-AQUA: A Many-Body CCSD(T) Water Potential, Including Four-Body Interactions, Demonstrates the Quantum Nature of Water from Clusters to the Liquid Phase. J Phys Chem Lett 2022; 13:5068-5074. [PMID: 35652912 DOI: 10.1021/acs.jpclett.2c00966] [Citation(s) in RCA: 37] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Many model potential energy surfaces (PESs) have been reported for water; however, none are strictly from "first-principles". Here we report such a potential, based on a many-body representation at the CCSD(T) level of theory up to the four-body interaction. The new PES is benchmarked for the isomers of the water hexamer for dissociation energies, harmonic frequencies, and unrestricted diffusion Monte Carlo (DMC) calculations of the zero-point energies of the Prism, Book, and Cage isomers. Dissociation energies of several isomers of the 20-mer agree well with recent benchmark energies. Exploratory DMC calculations on this cluster verify the robustness of the new PES for quantum simulations. The accuracy and speed of the new PES are demonstrated for standard condensed phase properties, i.e., the radial distribution function and the self-diffusion constant. Quantum effects are shown to be substantial for these observables and also needed to bring theory into excellent agreement with experiment.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Chen Qu
- Independent Researcher, Toronto, Ontario, Canada
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
2
|
Schmidt M, Millar J, Roy PN. Path integral simulations of confined parahydrogen molecules within clathrate hydrates: Merging low temperature dynamics with the zero-temperature limit. J Chem Phys 2022; 156:014303. [PMID: 34998330 DOI: 10.1063/5.0076386] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Clathrate hydrates, or cages comprised solely of water molecules, have long been investigated as a clean storage facility for hydrogen molecules. A breakthrough occurred when hydrogen molecules were experimentally placed within a structure-II clathrate hydrate, which sparked much interest to determine their feasibility for energy storage [Mao et al., Science 297, 2247-2249 (2002)]. We use Path Integral Molecular Dynamics (PIMD) and Langevin equation Path Integral Ground State (LePIGS) for finite temperature and zero-temperature studies, respectively, to determine parahydrogen occupancy properties in the small dodecahedral (512) and large hexakaidecahedral (51264) sized cages that comprise the structure-II unit cell. We look at energetic and structural properties of small clusters of hydrogen, treated as point-like particles, confined within each of the different sized clathrates, and treated as rigid, to determine energetic and structural properties in the zero-temperature limit. Our predicted hydrogen occupancy within these two cage sizes is consistent with previous literature values. We then calculate the energies as a function of temperature and merge the low temperature results calculated using finite temperature PIMD with the zero-temperature results using LePIGS, demonstrating that the two methods are compatible.
Collapse
Affiliation(s)
- Matthew Schmidt
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Jayme Millar
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
3
|
DiRisio RJ, Lu F, McCoy AB. GPU-Accelerated Neural Network Potential Energy Surfaces for Diffusion Monte Carlo. J Phys Chem A 2021; 125:5849-5859. [PMID: 34165989 DOI: 10.1021/acs.jpca.1c03709] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Diffusion Monte Carlo (DMC) provides a powerful method for understanding the vibrational landscape of molecules that are not well-described by conventional methods. The most computationally demanding step of these calculations is the evaluation of the potential energy. In this work, a general approach is developed in which a neural network potential energy surface is trained by using data generated from a small-scale DMC calculation. Once trained, the neural network can be evaluated by using highly parallelizable calls to a graphics processing unit (GPU). The power of this approach is demonstrated for DMC simulations on H2O, CH5+, and (H2O)2. The need to include permutation symmetry in the neural network potentials is explored and incorporated into the molecular descriptors of CH5+ and (H2O)2. It is shown that the zero-point energies and wave functions obtained by using the neural network potentials are nearly identical to the results obtained when using the potential energy surfaces that were used to train the neural networks at a substantial savings in the computational requirements of the simulations.
Collapse
Affiliation(s)
- Ryan J DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Fenris Lu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
4
|
Sahoo T, Serwatka T, Roy PN. A path integral ground state approach for asymmetric top rotors with nuclear spin symmetry: Application to water chains. J Chem Phys 2021; 154:244305. [PMID: 34241367 DOI: 10.1063/5.0053051] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A path integral ground state (PIGS) approach for the simulation of asymmetric top rotors is presented. The method is based on Monte Carlo sampling of angular degrees of freedom. A symmetry-adapted rotational density matrix is used to account for nuclear spin statistics. To illustrate the method, ground-state properties of collections of para-water molecules confined to a one-dimensional lattice are computed. Those include energetic and structural observables. An advantage of the PIGS method is that expectation values can be obtained directly since the square of the wavefunction is sampled during a simulation. To benchmark the method, ground state energies and orientational distributions are computed using exact diagonalization for a single para-water molecule in an external field using a finite basis of symmetric top eigenfunctions. Benchmark results are also provided for N = 2 para-water molecules pinned to lattice sites at various distances to sample the crossover from hydrogen bonding to the dipole-dipole interaction regime. Excellent agreement between the PIGS results and the finite basis set calculations is observed. A thorough analysis of the convergence in terms of the imaginary time propagation length and systematic Trotter error is performed. The PIGS approach is then applied to a chain of N = 11 water molecules, and an equation of state is constructed in terms of the intermolecular separation. Ordering effects are also studied, and a transition between hydrogen bonding to dipole-dipole alignment is observed. The method is scalable and can also be applied in higher dimensions.
Collapse
Affiliation(s)
- Tapas Sahoo
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Tobias Serwatka
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
5
|
Finney JM, DiRisio RJ, McCoy AB. Guided Diffusion Monte Carlo: A Method for Studying Molecules and Ions That Display Large Amplitude Vibrational Motions. J Phys Chem A 2020; 124:9567-9577. [DOI: 10.1021/acs.jpca.0c07181] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jacob M. Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Ryan J. DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
6
|
Lee VGM, Vetterli NJ, Boyer MA, McCoy AB. Diffusion Monte Carlo Studies on the Detection of Structural Changes in the Water Hexamer upon Isotopic Substitution. J Phys Chem A 2020; 124:6903-6912. [DOI: 10.1021/acs.jpca.0c05686] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Victor G. M. Lee
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Nicholas J. Vetterli
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mark A. Boyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
7
|
Lee VGM, McCoy AB. An Efficient Approach for Studies of Water Clusters Using Diffusion Monte Carlo. J Phys Chem A 2019; 123:8063-8070. [DOI: 10.1021/acs.jpca.9b06444] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Victor G. M. Lee
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
8
|
Mallory JD, Mandelshtam VA. Quantum-induced solid-solid transitions and melting in the Lennard-Jones LJ 38 cluster. J Chem Phys 2018; 149:104305. [DOI: 10.1063/1.5050410] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Joel D. Mallory
- Department of Chemistry, University of California, Irvine, California 92617, USA
| | | |
Collapse
|
9
|
Mallory JD, Mandelshtam VA. Nuclear Quantum Effects and Thermodynamic Properties for Small (H2O)1–21X– Clusters (X– = F–, Cl–, Br–, I–). J Phys Chem A 2018; 122:4167-4180. [DOI: 10.1021/acs.jpca.8b00917] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Joel D. Mallory
- Department of Chemistry, University of California, Irvine, California 92697, United States
| | | |
Collapse
|
10
|
Schmidt M, Roy PN. Path integral molecular dynamic simulation of flexible molecular systems in their ground state: Application to the water dimer. J Chem Phys 2018; 148:124116. [DOI: 10.1063/1.5017532] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Matthew Schmidt
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
11
|
Curotto E, Mella M. Diffusion Monte Carlo simulations of gas phase and adsorbed D 2-(H 2) n clusters. J Chem Phys 2018; 148:102315. [PMID: 29544319 DOI: 10.1063/1.5000372] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We have computed ground state energies and analyzed radial distributions for several gas phase and adsorbed D2(H2)n and HD(H2)n clusters. An external model potential designed to mimic ionic adsorption sites inside porous materials is used [M. Mella and E. Curotto, J. Phys. Chem. A 121, 5005 (2017)]. The isotopic substitution lowers the ground state energies by the expected amount based on the mass differences when these are compared with the energies of the pure clusters in the gas phase. A similar impact is found for adsorbed aggregates. The dissociation energy of D2 from the adsorbed clusters is always much higher than that of H2 from both pure and doped aggregates. Radial distributions of D2 and H2 are compared for both the gas phase and adsorbed species. For the gas phase clusters, two types of hydrogen-hydrogen interactions are considered: one based on the assumption that rotations and translations are adiabatically decoupled and the other based on nonisotropic four-dimensional potential. In the gas phase clusters of sufficiently large size, we find the heavier isotopomer more likely to be near the center of mass. However, there is a considerable overlap among the radial distributions of the two species. For the adsorbed clusters, we invariably find the heavy isotope located closer to the attractive interaction source than H2, and at the periphery of the aggregate, H2 molecules being substantially excluded from the interaction with the source. This finding rationalizes the dissociation energy results. For D2-(H2)n clusters with n≥12, such preference leads to the desorption of D2 from the aggregate, a phenomenon driven by the minimization of the total energy that can be obtained by reducing the confinement of (H2)12. The same happens for (H2)13, indicating that such an effect may be quite general and impact on the absorption of quantum species inside porous materials.
Collapse
Affiliation(s)
- E Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | - M Mella
- Dipartimento di Scienza ed Alta Tecnologia, Università degli Studi dell'Insubria, Via Valleggio 11, 22100 Como, Italy
| |
Collapse
|
12
|
Rugeles EJ, Ujevic S, Vitiello SA. Solid ^{4}He and the diffusion Monte Carlo method: A study of their properties. Phys Rev E 2017; 96:043306. [PMID: 29347541 DOI: 10.1103/physreve.96.043306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2017] [Indexed: 06/07/2023]
Abstract
Properties of helium atoms in the solid phase are investigated using the multiweight diffusion Monte Carlo method. Two different importance function transformations are used in two series of independent calculations. The kinetic energy is estimated for both the solid and liquid phases of ^{4}He. We estimate the melting and freezing densities, among other properties of interest. Our estimates are compared with experimental values. We discuss why walkers biased by two distinctly different guiding functions do not lead to noticeable changes in the reported results. Criticisms concerning the bias introduced into our estimates by population control and system size effects are considered.
Collapse
Affiliation(s)
- E J Rugeles
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas-UNICAMP, 13083-970 Campinas, SP, Brazil
| | - Sebastian Ujevic
- Departamento de Ciências Exatas-EEIMVR, Universidade Federal Fluminense, 27255-125 Volta Redonda, RJ, Brazil
| | - S A Vitiello
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas-UNICAMP, 13083-970 Campinas, SP, Brazil
| |
Collapse
|
13
|
Mallory JD, Mandelshtam VA. Quantum Melting and Isotope Effects from Diffusion Monte Carlo Studies of p-H2 Clusters. J Phys Chem A 2017; 121:6341-6348. [DOI: 10.1021/acs.jpca.7b06649] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Joel D. Mallory
- Department of Chemistry, University of California, Irvine, California 92697, United States
| | | |
Collapse
|
14
|
Mallory JD, Mandelshtam VA. Diffusion Monte Carlo studies of MB-pol (H2O)2−6 and (D2O)2−6 clusters: Structures and binding energies. J Chem Phys 2016. [DOI: 10.1063/1.4960610] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joel D. Mallory
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | | |
Collapse
|
15
|
Christensen HM, Jake LC, Curotto E. Smart darting diffusion Monte Carlo: Applications to lithium ion-Stockmayer clusters. J Chem Phys 2016; 144:174115. [PMID: 27155633 DOI: 10.1063/1.4948562] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In a recent investigation [K. Roberts et al., J. Chem. Phys. 136, 074104 (2012)], we have shown that, for a sufficiently complex potential, the Diffusion Monte Carlo (DMC) random walk can become quasiergodic, and we have introduced smart darting-like moves to improve the sampling. In this article, we systematically characterize the bias that smart darting moves introduce in the estimate of the ground state energy of a bosonic system. We then test a simple approach to eliminate completely such bias from the results. The approach is applied for the determination of the ground state of lithium ion-n-dipoles clusters in the n = 8-20 range. For these, the smart darting diffusion Monte Carlo simulations find the same ground state energy and mixed-distribution as the traditional approach for n < 14. In larger systems we find that while the ground state energies agree quantitatively with or without smart darting moves, the mixed-distributions can be significantly different. Some evidence is offered to conclude that introducing smart darting-like moves in traditional DMC simulations may produce a more reliable ground state mixed-distribution.
Collapse
Affiliation(s)
- H M Christensen
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | - L C Jake
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | - E Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| |
Collapse
|
16
|
Affiliation(s)
- Matúš Dubecký
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
17
|
Mallory JD, Mandelshtam VA. Binding energies from diffusion Monte Carlo for the MB-pol H2O and D2O dimer: A comparison to experimental values. J Chem Phys 2015; 143:144303. [DOI: 10.1063/1.4932596] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joel D. Mallory
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697, USA
| | - Vladimir A. Mandelshtam
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697, USA
| |
Collapse
|