1
|
Welsh BA, Urbina AS, Ho TA, Rempe SL, Slipchenko LV, Zwier TS. Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2) n, n = 1,2. J Phys Chem A 2024; 128:1124-1133. [PMID: 38306293 DOI: 10.1021/acs.jpca.3c07930] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2024]
Abstract
Pyrimidine has two in-plane CH(δ+)/N̈(δ-)/CH(δ+) binding sites that are complementary to the (δ-/2δ+/δ-) quadrupole moment of CO2. We recorded broadband microwave spectra over the 7.5-17.5 GHz range for pyrimidine-(CO2)n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO2 complex is planar, with CO2 occupying one of the quadrupolar binding sites, forming a structure in which the CO2 is stabilized in the plane by interactions with the C-H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191-195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C2v symmetry in which the second CO2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO2 is -13.7 kJ mol-1, including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO2)2 is almost exactly double that size, indicating little interaction between the two CO2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.
Collapse
Affiliation(s)
- Blair A Welsh
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| | - Andres S Urbina
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907-1393, United States
| | - Tuan A Ho
- Geochemistry Department, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Susan L Rempe
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907-1393, United States
| | - Timothy S Zwier
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| |
Collapse
|
2
|
Mendonça BHS, de Moraes EE, Kirch A, Batista RJC, de Oliveira AB, Barbosa MC, Chacham H. Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models. J Phys Chem B 2023; 127:8634-8643. [PMID: 37754781 DOI: 10.1021/acs.jpcb.3c02889] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/28/2023]
Abstract
In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.
Collapse
Affiliation(s)
- Bruno H S Mendonça
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, Belo Horizonte 30123-970, MG, Brazil
| | - Elizane E de Moraes
- Instituto de Física, Universidade Federal da Bahia, Campus Universitário de Ondina, Salvador 40210-340, BA, Brazil
| | - Alexsandro Kirch
- Instituto de Física, Universidade de São Paulo, CP 66318, São Paulo 05315-970, SP, Brazil
| | - Ronaldo J C Batista
- Departamento de Física, Universidade Federal de Ouro Preto, Campus Morro do Cruzeiro, Ouro Preto 35400-000, MG, Brazil
| | - Alan B de Oliveira
- Departamento de Física, Universidade Federal de Ouro Preto, Campus Morro do Cruzeiro, Ouro Preto 35400-000, MG, Brazil
| | - Marcia C Barbosa
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil
| | - Hélio Chacham
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, Belo Horizonte 30123-970, MG, Brazil
| |
Collapse
|
3
|
Riera M, Knight C, Bull-Vulpe EF, Zhu X, Agnew H, Smith DGA, Simmonett AC, Paesani F. MBX: A many-body energy and force calculator for data-driven many-body simulations. J Chem Phys 2023; 159:054802. [PMID: 37526156 PMCID: PMC10550339 DOI: 10.1063/5.0156036] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 07/11/2023] [Indexed: 08/02/2023] Open
Abstract
Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.
Collapse
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Argonne National Laboratory, Computational Science Division, Lemont, Illinois 60439, USA
| | - Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Xuanyu Zhu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | | | - Andrew C. Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | | |
Collapse
|
4
|
Méndez E, Videla PE, Laria D. Collective Proton Transfers in Cyclic Water-Ammonia Tetramers: A Path Integral Machine-Learning Study. J Phys Chem A 2023; 127:1839-1848. [PMID: 36794937 DOI: 10.1021/acs.jpca.2c07994] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
We present results from machine-learning-based path integral molecular dynamics simulations that describe isomerization paths articulated via collective proton transfers along mixed, cyclic tetramers combining water and ammonia at cryogenic conditions. The net result of such isomerizations is a reverse of the chirality of the global hydrogen-bonding architecture along the different cyclic moieties. In monocomponent tetramers, the classical free energy profiles associated with these isomerizations present the usual symmetric double-well characteristics whereas the reactive paths exhibit full concertedness among the different intermolecular transfer processes. Contrastingly, in mixed water/ammonia tetramers, the incorporation of a second component introduces imbalances in the strengths of the different hydrogen bonds leading to a partial loss of concertedness, most notably at the vicinity of the transition state. As such, the highest and lowest degrees of progression are registered along OH···N and O···HN coordinations, respectively. These characteristics lead to polarized transition state scenarios akin to solvent-separated ion-pair configurations. The explicit incorporation of nuclear quantum effects promotes drastic depletions in the activation free energies and modifications in the overall shape of the profiles which include central plateau-like stages, indicating the prevalence of deep tunneling regimes. On the other hand, the quantum treatment of the nuclei partially restores the degree of concertedness among the evolutions of the individual transfers.
Collapse
Affiliation(s)
- Emilio Méndez
- Departamento de Química Inorgánica, Analítica y Química-Física and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E Videla
- Department of Chemistry and Energy Sciences Institute, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Daniel Laria
- Departamento de Química Inorgánica, Analítica y Química-Física and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina.,Departamento de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
5
|
Lamaire A, Cools-Ceuppens M, Bocus M, Verstraelen T, Van Speybroeck V. Quantum Free Energy Profiles for Molecular Proton Transfers. J Chem Theory Comput 2022; 19:18-24. [PMID: 36563337 PMCID: PMC9835831 DOI: 10.1021/acs.jctc.2c00874] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Although many molecular dynamics simulations treat the atomic nuclei as classical particles, an adequate description of nuclear quantum effects (NQEs) is indispensable when studying proton transfer reactions. Herein, quantum free energy profiles are constructed for three typical proton transfers, which properly take NQEs into account using the path integral formalism. The computational cost of the simulations is kept tractable by deriving machine learning potentials. It is shown that the classical and quasi-classical centroid free energy profiles of the proton transfers deviate substantially from the exact quantum free energy profile.
Collapse
|
6
|
Zhuang D, Riera M, Zhou R, Deary A, Paesani F. Hydration Structure of Na + and K + Ions in Solution Predicted by Data-Driven Many-Body Potentials. J Phys Chem B 2022; 126:9349-9360. [PMID: 36326071 DOI: 10.1021/acs.jpcb.2c05674] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
The hydration structure of Na+ and K+ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na+ and K+ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.
Collapse
Affiliation(s)
- Debbie Zhuang
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Ruihan Zhou
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Alexander Deary
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California92093, United States
| |
Collapse
|
7
|
Bull-Vulpe EF, Riera M, Bore SL, Paesani F. Data-Driven Many-Body Potential Energy Functions for Generic Molecules: Linear Alkanes as a Proof-of-Concept Application. J Chem Theory Comput 2022. [PMID: 36113028 DOI: 10.1021/acs.jctc.2c00645] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a generalization of the many-body energy (MB-nrg) theoretical/computational framework that enables the development of data-driven potential energy functions (PEFs) for generic covalently bonded molecules, with arbitrary quantum mechanical accuracy. The "nearsightedness of electronic matter" is exploited to define monomers as "natural building blocks" on the basis of their distinct chemical identity. The energy of generic molecules is then expressed as a sum of individual many-body energies of incrementally larger subsystems. The MB-nrg PEFs represent the low-order n-body energies, with n = 1-4, using permutationally invariant polynomials derived from electronic structure data carried out at an arbitrary quantum mechanical level of theory, while all higher-order n-body terms (n > 4) are represented by a classical many-body polarization term. As a proof-of-concept application of the general MB-nrg framework, we present MB-nrg PEFs for linear alkanes. The MB-nrg PEFs are shown to accurately reproduce reference energies, harmonic frequencies, and potential energy scans of alkanes, independently of their length. Since, by construction, the MB-nrg framework introduced here can be applied to generic covalently bonded molecules, we envision future computer simulations of complex molecular systems using data-driven MB-nrg PEFs, with arbitrary quantum mechanical accuracy.
Collapse
Affiliation(s)
- Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Sigbjørn L. Bore
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
8
|
Méndez E, Videla PE, Laria D. Equilibrium and Dynamical Characteristics of Hydrogen Bond Bifurcations in Water-Water and Water-Ammonia Dimers: A Path Integral Molecular Dynamics Study. J Phys Chem A 2022; 126:4721-4733. [PMID: 35834556 DOI: 10.1021/acs.jpca.2c02525] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present path integral molecular dynamics results that describe the effects of nuclear quantum fluctuations on equilibrium and dynamical characteristics pertaining to bifurcation pathways in hydrogen bonded dimers combining water and ammonia, at cryogenic temperatures of the order of 20 K. Along these isomerizations, the hydrogen atoms in the molecules acting as hydrogen-bond donors interchange their original dangling/connective characters. Our results reveal that the resulting quantum transition paths comprise three stages: the initial and final ones involve overall rotations during which the two protons retain their classical-like characteristics. Effects from quantum fluctuation are clearly manifested in the changes operated at the intermediate passages over transition states, as the spatial extents of the protons stretch over typical lengths comparable to the distances between connective and dangling basins of attractions. Consequently, the classical over-the-hill path is replaced by a tunneling controlled mechanism which, within the path integral perspective, can be cast in terms of concerted inter-basin migrations of polymer beads from dangling-to-connective and from connective-to-dangling, at practically no energy costs. We also estimated the characteristic timescales describing such interconversions within the approximate ring polymer rate theory. Effects derived from full and partial deuteration are also discussed.
Collapse
Affiliation(s)
- Emilio Méndez
- Departamento de Química Inorgánica, Analítica y Química-Física and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E Videla
- Department of Chemistry and Energy Sciences Institute, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Daniel Laria
- Departamento de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
9
|
Bull-Vulpe EF, Riera M, Götz AW, Paesani F. MB-Fit: Software infrastructure for data-driven many-body potential energy functions. J Chem Phys 2021; 155:124801. [PMID: 34598567 DOI: 10.1063/5.0063198] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.
Collapse
Affiliation(s)
- Ethan F Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
10
|
Caruso A, Paesani F. Data-driven many-body models enable a quantitative description of chloride hydration from clusters to bulk. J Chem Phys 2021; 155:064502. [PMID: 34391363 DOI: 10.1063/5.0059445] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
We present a new data-driven potential energy function (PEF) describing chloride-water interactions, which is developed within the many-body-energy (MB-nrg) theoretical framework. Besides quantitatively reproducing low-order many-body energy contributions, the new MB-nrg PEF is able to correctly predict the interaction energies of small chloride-water clusters calculated at the coupled cluster level of theory. Importantly, classical and quantum molecular dynamics simulations of a single chloride ion in water demonstrate that the new MB-nrg PEF predicts x-ray spectra in close agreement with the experimental results. Comparisons with an popular empirical model and a polarizable PEF emphasize the importance of an accurate representation of short-range many-body effect while demonstrating that pairwise additive representations of chloride-water and water-water interactions are inadequate for correctly representing the hydration structure of chloride in both gas-phase clusters and solution. We believe that the analyses presented in this study provide additional evidence for the accuracy and predictive ability of the MB-nrg PEFs, which can then enable more realistic simulations of ionic aqueous systems in different environments.
Collapse
Affiliation(s)
- Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
11
|
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
|
12
|
Lambros E, Paesani F. How good are polarizable and flexible models for water: Insights from a many-body perspective. J Chem Phys 2020; 153:060901. [DOI: 10.1063/5.0017590] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
13
|
Méndez E, Laria D. Nuclear quantum effects on the hydrogen bond donor-acceptor exchange in water-water and water-methanol dimers. J Chem Phys 2020; 153:054302. [PMID: 32770908 DOI: 10.1063/5.0016122] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present results from path integral molecular dynamics simulations that describe effects from the explicit incorporation of nuclear quantum fluctuations on the topology of the free energy associated with the geared exchange of hydrogen bonds in the water-water dimer. Compared to the classical treatment, our results reveal important reductions in the free energy barriers and changes at a qualitative level in the overall profile. Most notable are those manifested by a plateau behavior, ascribed to nuclear tunneling, which bridges reactant and product states, contrasting with the usual symmetric double-well profile. The characteristics of the proton localizations along the pathway are examined. An imaginary time analysis of the rotational degrees of freedom of the partners in the dimer at the vicinities of transition states shows a clear "anticorrelation" between intermolecular interactions coupling beads localized in connective and dangling basins of attractions. As such, the transfer is operated by gradual concerted inter-basin migrations in opposite directions, at practically no energy costs. Modifications operated by partial deuteration and by the asymmetries in the hydrogen bonding characteristics prevailing in water-methanol heterodimers are also examined.
Collapse
Affiliation(s)
- Emilio Méndez
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Daniel Laria
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| |
Collapse
|
14
|
Zhai Y, Caruso A, Gao S, Paesani F. Active learning of many-body configuration space: Application to the Cs+–water MB-nrg potential energy function as a case study. J Chem Phys 2020; 152:144103. [DOI: 10.1063/5.0002162] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Yaoguang Zhai
- Department of Computer Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
| | - Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Sicun Gao
- Department of Computer Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
15
|
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
|
16
|
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
|