1
|
Beckmann R, Schran C, Brieuc F, Marx D. Theoretical infrared spectroscopy of protonated methane isotopologues. Phys Chem Chem Phys 2024; 26:22846-22852. [PMID: 39171731 DOI: 10.1039/d4cp02295e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
The vibrational spectroscopy of protonated methane and its mixed hydrogen/deuterium isotopologues remains a challenge to both experimental and computational spectroscopy due to the iconic floppiness of CH5+. Here, we compute the finite-temperature broadband infrared spectra of CH5+ and all its isotopologues, i.e. CHnD5-n+ up to CD5+, from path integral molecular dynamics in conjunction with interactions and dipoles computed consistently at CCSD(T) coupled cluster accuracy. The potential energy and dipole moment surfaces have been accurately represented in full dimensionality in terms of high-dimensional neural networks. The resulting computational efficiency allows us to establish CCSD(T) accuracy at the level of converged path integral simulations. For all six isotopologues, the computed broadband spectra compare very favorably to the available experimental broadband spectra obtained from laser induced reactions action vibrational spectroscopy. The current approach is found to consistently and significantly improve on previous calculations of these broadband vibrational spectra and defines the new cutting-edge for what has been dubbed the "enfant terrible" of molecular spectroscopy in view of its pronounced large-amplitude motion that involves all intramolecular degrees of freedom.
Collapse
Affiliation(s)
- Richard Beckmann
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Fabien Brieuc
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| |
Collapse
|
2
|
Gomez A, Thompson WH, Laage D. Neural-network-based molecular dynamics simulations reveal that proton transport in water is doubly gated by sequential hydrogen-bond exchange. Nat Chem 2024:10.1038/s41557-024-01593-y. [PMID: 39164581 DOI: 10.1038/s41557-024-01593-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2024] [Accepted: 07/03/2024] [Indexed: 08/22/2024]
Abstract
The transport of excess protons in water is central to acid-base chemistry, biochemistry and energy production. However, elucidating its mechanism has been challenging. Recent nonlinear vibrational spectroscopy experiments could not be explained by existing models. Here we use both vibrational spectroscopy calculations and neural-network-based molecular dynamics simulations that account for nuclear quantum effects for all atoms to determine the proton transport mechanism. Our simulations reveal an equilibrium between two stable proton-localized structures with distinct Eigen-like and Zundel-like hydrogen-bond motifs. Proton transport follows a three-step mechanism gated by two successive hydrogen-bond exchanges: the first reduces the proton-acceptor water coordination, leading to proton transfer, and the second, the rate-limiting step, prevents rapid back-transfer by increasing the proton-donor coordination. This sequential mechanism is consistent with experimental characterizations of proton diffusion, explaining the low activation energy and the prolonged intermediate lifetimes in vibrational spectroscopy. These results are crucial for understanding proton dynamics in biochemical and technological systems.
Collapse
Affiliation(s)
- Axel Gomez
- PASTEUR, Department of Chemistry, École normale supérieure, PSL University, Sorbonne Université, CNRS, Paris, France
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, KS, USA
| | - Damien Laage
- PASTEUR, Department of Chemistry, École normale supérieure, PSL University, Sorbonne Université, CNRS, Paris, France.
| |
Collapse
|
3
|
Litman Y, Kapil V, Feldman YMY, Tisi D, Begušić T, Fidanyan K, Fraux G, Higer J, Kellner M, Li TE, Pós ES, Stocco E, Trenins G, Hirshberg B, Rossi M, Ceriotti M. i-PI 3.0: A flexible and efficient framework for advanced atomistic simulations. J Chem Phys 2024; 161:062504. [PMID: 39140447 DOI: 10.1063/5.0215869] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2024] [Accepted: 07/11/2024] [Indexed: 08/15/2024] Open
Abstract
Atomic-scale simulations have progressed tremendously over the past decade, largely thanks to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD, and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows for deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities.
Collapse
Affiliation(s)
- Yair Litman
- Y. Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Venkat Kapil
- Y. Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Department of Physics and Astronomy, University College London, 17-19 Gordon St, London WC1H 0AH, United Kingdom
- Thomas Young Centre and London Centre for Nanotechnology, 19 Gordon St, London WC1H 0AH, United Kingdom
| | | | - Davide Tisi
- Laboratory of Computational Science and Modeling, Institut des Matériaux, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Tomislav Begušić
- Div. of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Karen Fidanyan
- MPI for the Structure and Dynamics of Matter, Hamburg, Germany
| | - Guillaume Fraux
- Laboratory of Computational Science and Modeling, Institut des Matériaux, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jacob Higer
- School of Physics, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Matthias Kellner
- Laboratory of Computational Science and Modeling, Institut des Matériaux, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Tao E Li
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Eszter S Pós
- MPI for the Structure and Dynamics of Matter, Hamburg, Germany
| | - Elia Stocco
- MPI for the Structure and Dynamics of Matter, Hamburg, Germany
| | - George Trenins
- MPI for the Structure and Dynamics of Matter, Hamburg, Germany
| | - Barak Hirshberg
- School of Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Mariana Rossi
- MPI for the Structure and Dynamics of Matter, Hamburg, Germany
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institut des Matériaux, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
4
|
Silarski M, Dziedzic-Kocurek K, Drużbicki K, Reterski R, Grabowski P, Krzystyniak M. Non-invasive detection of hazardous materials with a thermal-to-epithermal neutron station: a feasibility study towards practical application. Sci Rep 2024; 14:18584. [PMID: 39127754 DOI: 10.1038/s41598-024-69290-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2024] [Accepted: 08/01/2024] [Indexed: 08/12/2024] Open
Abstract
The growing scale of the devastation that even a single terrorist attack can cause requires more effective methods for the detection of hazardous materials. In particular, there are no solutions for effectively monitoring threats at sea, both for the off-shore infrastructure and ports. Currently, state-of-the-art detection methods determine the density distribution and the shapes of tested subjects but only allow for a limited degree of substance identification. This work aims to present a feasibility study of the possible usage of several methods available on the thermal-to-epithermal neutron station, VESUVIO, at the ISIS neutron and muon spallation source, UK, for the detection of hazardous materials. To this end, we present the results of a series of experiments performed concurrently employing neutron transmission and Compton scattering using melamine, a commonly used explosive surrogate, in order to determine its signal characteristics and limits of detection and quantitation. The experiments are supported by first-principles modelling, providing detailed scrutiny of the material structure and the nuclear dynamics behind the neutron scattering observables.
Collapse
Affiliation(s)
- Michał Silarski
- M. Smoluchowski Institute of Physics of the Jagiellonian University, Łojasiewicza 11, 30-348, Cracow, Poland
| | - Katarzyna Dziedzic-Kocurek
- M. Smoluchowski Institute of Physics of the Jagiellonian University, Łojasiewicza 11, 30-348, Cracow, Poland
| | - Kacper Drużbicki
- Centre of Molecular and Macromolecular Studies, Polish Academy of Sciences, Sienkiewicza 112, 90-363, Lodz, Poland
| | - Radosław Reterski
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387, Cracow, Poland
| | - Patryk Grabowski
- M. Smoluchowski Institute of Physics of the Jagiellonian University, Łojasiewicza 11, 30-348, Cracow, Poland
| | - Matthew Krzystyniak
- ISIS Neutron and Muon Source, STFC Rutherford Appleton Laboratory, Oxford, OX11 0QX, UK.
| |
Collapse
|
5
|
Arandhara M, Ramesh SG. Nuclear quantum effects in gas-phase ethylene glycol. Phys Chem Chem Phys 2024; 26:19529-19542. [PMID: 38979630 DOI: 10.1039/d4cp00700j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Path integral molecular simulations are used to explore the nuclear quantum effects (NQEs) on the structure, dihedral landscape and infrared spectrum of ethylene glycol. The simulations are carried out on a new reaction surface Hamiltonian-based model potential energy surface, with special focus on the role of the OCCO and HOCC dihedrals. In contrast with classical simulations, we analyse how the intramolecular interactions between the OH groups change due to zero-point effects as well as temperature. These are found to be weak. The NQEs on the free energy profile along the OCCO dihedral are analysed, where notable effects are seen at low temperatures and found to be correlated with the radii of gyration of the atoms. Finally, the power spectrum of the molecule from path integral simulations is compared with the experimental infrared spectrum, yielding a good agreement of band positions.
Collapse
Affiliation(s)
- Mrinal Arandhara
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India.
| | - Sai G Ramesh
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India.
| |
Collapse
|
6
|
Murakami T, Takayanagi T. Computational study of the post-transition state dynamics for the OH + CH 3OH reaction probed by photodetachment of the CH 3O -(H 2O) anion. Phys Chem Chem Phys 2024; 26:19195-19206. [PMID: 38956990 DOI: 10.1039/d4cp01466a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Dissociative photodetachment dynamics simulations were conducted to study the CH3O-(H2O) → CH3O + H2O + e- reaction using classical molecular dynamics (MD) and ring-polymer molecular dynamics (RPMD) techniques on two newly formulated neutral potential energy surfaces (PES1 and PES2) by different research groups. While the dissociation dynamics exhibited similarities between classical MD and RPMD, there were noticeable differences in the fluctuation of probability densities for the internal modes due to nuclear quantum effects. Upon comparison of our findings with experimental data concerning the electron binding energy distribution and photofragment relative energy, it suggests that the potential energy landscapes of PES2 are reasonably precise. The time evolution of occupied vibrational states of the H2O photofragment is presented in this study.
Collapse
Affiliation(s)
- Tatsuhiro Murakami
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
- Department of Materials & Life Sciences, Faculty of Science & Technology, Sophia University, 7-1 Kioicho, Chiyoda-ku, Tokyo, 102-8554, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| |
Collapse
|
7
|
Zhang X, Li C, Ye HZ, Berkelbach TC, Chan GKL. Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics. J Chem Phys 2024; 161:014109. [PMID: 38949583 DOI: 10.1063/5.0212274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.
Collapse
Affiliation(s)
- Xing Zhang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Chenghan Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
8
|
Angiolari F, Mandelli G, Huppert S, Aieta C, Spezia R. Eludicating Heavy-Atom-Tunneling Kinetics in the Cope Rearrangement of Semibullvalene. Chemistry 2024:e202401000. [PMID: 38924666 DOI: 10.1002/chem.202401000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2024] [Revised: 06/07/2024] [Accepted: 06/26/2024] [Indexed: 06/28/2024]
Abstract
In this work, we characterize the temperature dependence of kinetic properties in heavy atom tunneling reactions by means of molecular dynamics simulations, including nuclear quantum effects (NQEs) via Path Integral theory. To this end, we consider the prototypical Cope rearrangement of semibullvalene. The reaction was studied in the 25-300 K temperature range observing that the inclusion of NQEs modifies the temperature behavior of both free energy barriers and dynamical recrossing factors with respect to classical dynamics. Notably, while in classical simulations the activation free energy shows a very weak temperature dependence, it becomes strongly dependent on temperature when NQEs are included. This temperature behavior shows a transition from a regime where the quantum effects are limited and can mainly be traced back to zero point energy, to a low temperature regime where tunneling plays a dominant role. In this regime, the free energy curve tunnels below the potential energy barrier along the reaction coordinate, allowing much faster reaction rates. Finally, the temperature dependence of the rate constants obtained from molecular dynamics simulations was compared with available experimental data and with semi-classical transition state theory calculations, showing comparable behaviors and similar transition temperatures from thermal to (deep) tunneling regime.
Collapse
Affiliation(s)
- Federica Angiolari
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 4 Place Jussieu, 75005, Paris, France
| | - Giacomo Mandelli
- Dipartimento di Chimica, Universitá degli Studi di Milano, 20133, Milano, Italy
| | - Simon Huppert
- Sorbonne Université, Institut de Nanosciences de Paris, UMR 7588 CNRS, 4 Place Jussieu, 75005, Paris, France
| | - Chiara Aieta
- Dipartimento di Chimica, Universitá degli Studi di Milano, 20133, Milano, Italy
| | - Riccardo Spezia
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 4 Place Jussieu, 75005, Paris, France
| |
Collapse
|
9
|
Li TE. Mesoscale Molecular Simulations of Fabry-Pérot Vibrational Strong Coupling. J Chem Theory Comput 2024. [PMID: 38912683 DOI: 10.1021/acs.jctc.4c00349] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
Developing theoretical frameworks for vibrational strong coupling (VSC) beyond the single-mode approximation is crucial for a comprehensive understanding of experiments with planar Fabry-Pérot cavities. Herein, a generalized cavity molecular dynamics (CavMD) scheme is developed to simulate VSC of a large ensemble of realistic molecules coupled to an arbitrary 1D or 2D photonic environment. This approach is built upon the Power-Zienau-Woolley Hamiltonian in the normal mode basis and uses a grid representation of the molecular ensembles to reduce the computational cost. When simulating the polariton dispersion relation for a homogeneous distribution of molecules in planar Fabry-Pérot cavities, our data highlight the importance of preserving the in-plane translational symmetry of the molecular distribution. In this homogeneous limit, CavMD yields the consistent polariton dispersion relation as an analytic theory, i.e., incorporating many cavity modes with varying in-plane wave vectors (k∥) produces the same spectrum as the system with a single cavity mode. Furthermore, CavMD reveals that the validity of the single-mode approximation is challenged when nonequilibrium polariton dynamics are considered, as polariton-polariton scattering occurs between modes with the nearest neighbor k∥. The procedure for numerically approaching the macroscopic limit is also demonstrated with CavMD by increasing the system size. Looking forward, our generalized CavMD approach may facilitate understanding vibrational polariton transport and condensation.
Collapse
Affiliation(s)
- Tao E Li
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
10
|
Althorpe SC. Path Integral Simulations of Condensed-Phase Vibrational Spectroscopy. Annu Rev Phys Chem 2024; 75:397-420. [PMID: 38941531 DOI: 10.1146/annurev-physchem-090722-124705] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2024]
Abstract
Recent theoretical and algorithmic developments have improved the accuracy with which path integral dynamics methods can include nuclear quantum effects in simulations of condensed-phase vibrational spectra. Such methods are now understood to be approximations to the delocalized classical Matsubara dynamics of smooth Feynman paths, which dominate the dynamics of systems such as liquid water at room temperature. Focusing mainly on simulations of liquid water and hexagonal ice, we explain how the recently developed quasicentroid molecular dynamics (QCMD), fast-QCMD, and temperature-elevated path integral coarse-graining simulations (Te PIGS) methods generate classical dynamics on potentials of mean force obtained by averaging over quantum thermal fluctuations. These new methods give very close agreement with one another, and the Te PIGS method has recently yielded excellent agreement with experimentally measured vibrational spectra for liquid water, ice, and the liquid-air interface. We also discuss the limitations of such methods.
Collapse
Affiliation(s)
- Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom;
| |
Collapse
|
11
|
Linker TM, Krishnamoorthy A, Daemen LL, Ramirez-Cuesta AJ, Nomura K, Nakano A, Cheng YQ, Hicks WR, Kolesnikov AI, Vashishta PD. Neutron scattering and neural-network quantum molecular dynamics investigation of the vibrations of ammonia along the solid-to-liquid transition. Nat Commun 2024; 15:3911. [PMID: 38724541 PMCID: PMC11082248 DOI: 10.1038/s41467-024-48246-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2023] [Accepted: 04/25/2024] [Indexed: 05/12/2024] Open
Abstract
Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs.
Collapse
Affiliation(s)
- T M Linker
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, CA, 90089-0242, USA
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California, 94025, USA
| | - A Krishnamoorthy
- Department of Mechanical Engineering Texas A&M, 400 Bizzell St, College Station, TX, 77843, USA
| | - L L Daemen
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - A J Ramirez-Cuesta
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - K Nomura
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, CA, 90089-0242, USA
| | - A Nakano
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, CA, 90089-0242, USA
| | - Y Q Cheng
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA.
| | - W R Hicks
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - A I Kolesnikov
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA.
| | - P D Vashishta
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, CA, 90089-0242, USA.
| |
Collapse
|
12
|
Zheng J, Frisch MJ. Multiple-time scale integration method based on an interpolated potential energy surface for ab initio path integral molecular dynamics. J Chem Phys 2024; 160:144111. [PMID: 38597307 DOI: 10.1063/5.0196634] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 03/25/2024] [Indexed: 04/11/2024] Open
Abstract
A new multiple-time scale integration method is presented that propagates ab initio path integral molecular dynamics (PIMD). This method uses a large time step to generate an approximate geometrical configuration whose energy and gradient are evaluated at the level of an ab initio method, and then, a more precise integration scheme, e.g., the Bulirsch-Stoer method or velocity Verlet integration with a smaller time step, is used to integrate from the previous step using the computationally efficient interpolated potential energy surface constructed from two consecutive points. This method makes the integration of PIMD more efficient and accurate compared with the velocity Verlet integration. A Nosé-Hoover chain thermostat combined with this new multiple-time scale method has good energy conservation even with a large time step, which is usually challenging in velocity Verlet integration for PIMD due to the very small chain mass when a large number of beads are used. The new method is used to calculate infrared spectra and free energy profiles to demonstrate its accuracy and capabilities.
Collapse
Affiliation(s)
- Jingjing Zheng
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Michael J Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| |
Collapse
|
13
|
London N, Limbu DK, Momeni MR, Shakib FA. DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations. J Chem Phys 2024; 160:132501. [PMID: 38557837 DOI: 10.1063/5.0197822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2024] [Accepted: 03/12/2024] [Indexed: 04/04/2024] Open
Abstract
DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman's path integral formalism of quantum mechanics. The flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.
Collapse
Affiliation(s)
- Nathan London
- Division of Energy, Matter and Systems, School of Science and Engineering, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA
| | - Dil K Limbu
- Division of Energy, Matter and Systems, School of Science and Engineering, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - Mohammad R Momeni
- Division of Energy, Matter and Systems, School of Science and Engineering, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA
| | - Farnaz A Shakib
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| |
Collapse
|
14
|
Lan J, Chergui M, Pasquarello A. Dynamics of the charge transfer to solvent process in aqueous iodide. Nat Commun 2024; 15:2544. [PMID: 38514610 PMCID: PMC11258362 DOI: 10.1038/s41467-024-46772-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Accepted: 03/05/2024] [Indexed: 03/23/2024] Open
Abstract
Charge-transfer-to-solvent states in aqueous halides are ideal systems for studying the electron-transfer dynamics to the solvent involving a complex interplay between electronic excitation and solvent polarization. Despite extensive experimental investigations, a full picture of the charge-transfer-to-solvent dynamics has remained elusive. Here, we visualise the intricate interplay between the dynamics of the electron and the solvent polarization occurring in this process. Through the combined use of ab initio molecular dynamics and machine learning methods, we investigate the structure, dynamics and free energy as the excited electron evolves through the charge-transfer-to-solvent process, which we characterize as a sequence of states denoted charge-transfer-to-solvent, contact-pair, solvent-separated, and hydrated electron states, depending on the distance between the iodine and the excited electron. Our assignment of the charge-transfer-to-solvent states is supported by the good agreement between calculated and measured vertical binding energies. Our results reveal the charge transfer process in terms of the underlying atomic processes and mechanisms.
Collapse
Affiliation(s)
- Jinggang Lan
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland.
- Department of Chemistry, New York University, New York, NY, 10003, USA.
- Simons Center for Computational Physical Chemistry at New York University, New York, NY, 10003, USA.
| | - Majed Chergui
- Lausanne Centre for Ultrafast Science (LACUS), ISIC, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
- Elettra - Sincrotrone Trieste, Area Science Park I - 34149, Trieste, Italy
| | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| |
Collapse
|
15
|
de la Puente M, Gomez A, Laage D. Neural Network-Based Sum-Frequency Generation Spectra of Pure and Acidified Water Interfaces with Air. J Phys Chem Lett 2024; 15:3096-3102. [PMID: 38470065 DOI: 10.1021/acs.jpclett.4c00113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2024]
Abstract
The affinity of hydronium ions (H3O+) for the air-water interface is a crucial question in environmental chemistry. While sum-frequency generation (SFG) spectroscopy has been instrumental in indicating the preference of H3O+ for the interface, key questions persist regarding the molecular origin of the SFG spectral changes in acidified water. Here we combine nanosecond long neural network (NN) reactive simulations of pure and acidified water slabs with NN predictions of molecular dipoles and polarizabilities to calculate SFG spectra of long reactive trajectories including proton transfer events. Our simulations show that H3O+ ions cause two distinct changes in phase-resolved SFG spectra: first, a low-frequency tail due to the vibrations of H3O+ and its first hydration shell, analogous to the bulk proton continuum, and second, an enhanced hydrogen-bonded band due to the ion-induced static field polarizing molecules in deeper layers. Our calculations confirm that changes in the SFG spectra of acidic solutions are caused by hydronium ions preferentially residing at the interface.
Collapse
Affiliation(s)
- Miguel de la Puente
- PASTEUR, Department of Chemistry, École Normale Supérieur, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Axel Gomez
- PASTEUR, Department of Chemistry, École Normale Supérieur, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Damien Laage
- PASTEUR, Department of Chemistry, École Normale Supérieur, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
16
|
Kim DS, Xu M, LeBeau JM. Modeling Temperature-Dependent Electron Thermal Diffuse Scattering via Machine-Learned Interatomic Potentials and Path-Integral Molecular Dynamics. PHYSICAL REVIEW LETTERS 2024; 132:086301. [PMID: 38457736 DOI: 10.1103/physrevlett.132.086301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 01/12/2024] [Accepted: 01/25/2024] [Indexed: 03/10/2024]
Abstract
Electron thermal diffuse scattering is shown to be sensitive to subtle changes in atomic vibrations and shows promise in assessing lattice dynamics at nanometer resolution. Here, we demonstrate that machine-learned interatomic potentials (MLIPs) and path-integral molecular dynamics can accurately capture the potential energy landscape and lattice dynamics needed to describe electron thermal diffuse scattering. Using SrTiO_{3} as a test bed at cryogenic and room temperatures, we compare electron thermal diffuse scattering simulations using different approximations to incorporate thermal motion. Only when the simulations are based on quantum mechanically accurate MLIPs in combination with path-integral molecular dynamics that include nuclear quantum effects is there excellent agreement with experiments.
Collapse
Affiliation(s)
- Dennis S Kim
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Michael Xu
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - James M LeBeau
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
17
|
Arandhara M, Ramesh SG. Nuclear quantum effects in gas-phase 2-fluoroethanol. Phys Chem Chem Phys 2024; 26:6885-6902. [PMID: 38333949 DOI: 10.1039/d3cp05657k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Torsional motions along the FCCO and HOCC dihedrals lead to the five unique conformations of 2-fluoroethanol, of which the conformer that is gauche along both dihedrals has the lowest energy. In this work, we explore how nuclear quantum effects (NQEs) manifest in the structural parameters of the lowest energy conformer, in the intramolecular free energy landscape along the FCCO and HOCC dihedrals, and also in the infrared spectrum of the title molecule, through the use of path integral simulations. We have first developed a full dimensional potential energy surface using the reaction surface Hamiltonian framework. On this potential, we have carried out path integral molecular dynamics simulations at several temperatures starting from the minimum energy well to explore structural influences of NQEs including geometrical markers of the interaction between the OH and F groups. From the computed free energy landscapes, significant reduction of the torsional barrier is found at low temperature near the cis region of the dihedrals, which can be understood through the trends in the radii of gyration of the atomic ring polymers. We find that the inclusion of NQEs in the computation of infrared spectrum is important to obtain good agreement with the experimental band positions.
Collapse
Affiliation(s)
- Mrinal Arandhara
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India.
| | - Sai G Ramesh
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India.
| |
Collapse
|
18
|
Kapil V, Kovács DP, Csányi G, Michaelides A. First-principles spectroscopy of aqueous interfaces using machine-learned electronic and quantum nuclear effects. Faraday Discuss 2024; 249:50-68. [PMID: 37799072 PMCID: PMC10845015 DOI: 10.1039/d3fd00113j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Accepted: 07/18/2023] [Indexed: 10/07/2023]
Abstract
Vibrational spectroscopy is a powerful approach to visualising interfacial phenomena. However, extracting structural and dynamical information from vibrational spectra is a challenge that requires first-principles simulations, including non-Condon and quantum nuclear effects. We address this challenge by developing a machine-learning enhanced first-principles framework to speed up predictive modelling of infrared, Raman, and sum-frequency generation spectra. Our approach uses machine learning potentials that encode quantum nuclear effects to generate quantum trajectories using simple molecular dynamics efficiently. In addition, we reformulate bulk and interfacial selection rules to express them unambiguously in terms of the derivatives of polarisation and polarisabilities of the whole system and predict these derivatives efficiently using fully-differentiable machine learning models of dielectric response tensors. We demonstrate our framework's performance by predicting the IR, Raman, and sum-frequency generation spectra of liquid water, ice and the water-air interface by achieving near quantitative agreement with experiments at nearly the same computational efficiency as pure classical methods. Finally, to aid the experimental discovery of new phases of nanoconfined water, we predict the temperature-dependent vibrational spectra of monolayer water across the solid-hexatic-liquid phases transition.
Collapse
Affiliation(s)
- Venkat Kapil
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.
| | | | - Gábor Csányi
- Engineering Laboratory, University of Cambridge, Cambridge, CB2 1PZ, UK
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.
| |
Collapse
|
19
|
Thakur AC, Remsing RC. Nuclear quantum effects in the acetylene:ammonia plastic co-crystal. J Chem Phys 2024; 160:024502. [PMID: 38189604 DOI: 10.1063/5.0179161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Accepted: 12/20/2023] [Indexed: 01/09/2024] Open
Abstract
Organic molecular solids can exhibit rich phase diagrams. In addition to structurally unique phases, translational and rotational degrees of freedom can melt at different state points, giving rise to partially disordered solid phases. The structural and dynamic disorder in these materials can have a significant impact on the physical properties of the organic solid, necessitating a thorough understanding of disorder at the atomic scale. When these disordered phases form at low temperatures, especially in crystals with light nuclei, the prediction of material properties can be complicated by the importance of nuclear quantum effects. As an example, we investigate nuclear quantum effects on the structure and dynamics of the orientationally disordered, translationally ordered plastic phase of the acetylene:ammonia (1:1) co-crystal that is expected to exist on the surface of Saturn's moon Titan. Titan's low surface temperature (∼90 K) suggests that the quantum mechanical behavior of nuclei may be important in this and other molecular solids in these environments. By using neural network potentials combined with ring polymer molecular dynamics simulations, we show that nuclear quantum effects increase orientational disorder and rotational dynamics within the acetylene:ammonia (1:1) co-crystal by weakening hydrogen bonds. Our results suggest that nuclear quantum effects are important to accurately model molecular solids and their physical properties in low-temperature environments.
Collapse
Affiliation(s)
- Atul C Thakur
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Richard C Remsing
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| |
Collapse
|
20
|
Dufils T, Schran C, Chen J, Geim AK, Fumagalli L, Michaelides A. Origin of dielectric polarization suppression in confined water from first principles. Chem Sci 2024; 15:516-527. [PMID: 38179530 PMCID: PMC10763014 DOI: 10.1039/d3sc04740g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Accepted: 11/23/2023] [Indexed: 01/06/2024] Open
Abstract
It has long been known that the dielectric constant of confined water should be different from that in bulk. Recent experiments have shown that it is vanishingly small, however the origin of the phenomenon remains unclear. Here we used ab initio molecular dynamics simulations (AIMD) and AIMD-trained machine-learning potentials to understand water's structure and electronic properties underpinning this effect. For the graphene and hexagonal boron-nitride substrates considered, we find that it originates in the spontaneous anti-parallel alignment of the water dipoles in the first two water layers near the solid interface. The interfacial layers exhibit net ferroelectric ordering, resulting in an overall anti-ferroelectric arrangement of confined water. Together with constrained hydrogen-bonding orientations, this leads to much reduced out-of-plane polarization. Furthermore, we directly contrast AIMD and simple classical force-field simulations, revealing important differences. This work offers insight into a property of water that is critical in modulating surface forces, the electric-double-layer formation and molecular solvation, and shows a way to compute it.
Collapse
Affiliation(s)
- T Dufils
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - C Schran
- Cavendish Laboratory, Department of Physics, University of Cambridge Cambridge CB3 0HE UK
- Lennard-Jones Centre, University of Cambridge Trinity Ln Cambridge CB2 1TN UK
| | - J Chen
- School of Physics, Peking University Beijing 100871 China
| | - A K Geim
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - L Fumagalli
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - A Michaelides
- Lennard-Jones Centre, University of Cambridge Trinity Ln Cambridge CB2 1TN UK
- Yusuf Hamied Department of Chemistry, University of Cambridge Lensfield Road Cambridge CB2 1EW UK
| |
Collapse
|
21
|
Zhang Y, Wang Y, Xu X, Chen Z, Yang Y. Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods. J Chem Theory Comput 2023; 19:9358-9368. [PMID: 38096546 DOI: 10.1021/acs.jctc.3c01037] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological, and materials systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.
Collapse
Affiliation(s)
- Yuzhe Zhang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
22
|
Lawrence JE, Lieberherr AZ, Fletcher T, Manolopoulos DE. Fast Quasi-Centroid Molecular Dynamics for Water and Ice. J Phys Chem B 2023; 127:9172-9180. [PMID: 37830934 PMCID: PMC10614180 DOI: 10.1021/acs.jpcb.3c05028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 09/21/2023] [Indexed: 10/14/2023]
Abstract
We describe how the fast quasi-centroid molecular dynamics (f-QCMD) method can be applied to condensed-phase systems by approximating the quasi-centroid potential of mean force as a sum of inter- and intramolecular corrections to the classical interaction potential. The corrections are found by using a regularized iterative Boltzmann inversion procedure to recover the inter- and intramolecular quasi-centroid distribution functions obtained from a path integral molecular dynamics simulation. The resulting methodology is found to give good agreement with a previously published QCMD dipole absorption spectrum for liquid water and satisfactory agreement for ice. It also gives good agreement with spectra from a recent implementation of CMD that uses a precomputed elevated temperature potential of mean force. Modern centroid molecular dynamics methods, therefore, appear to be reaching a consensus regarding the impact of nuclear quantum effects on the vibrational spectra of water and ice.
Collapse
Affiliation(s)
| | - Annina Z. Lieberherr
- Physical
and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United
Kingdom
| | - Theo Fletcher
- Physical
and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United
Kingdom
| | - David E. Manolopoulos
- Physical
and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United
Kingdom
| |
Collapse
|
23
|
Zhang Y, Jiang B. Universal machine learning for the response of atomistic systems to external fields. Nat Commun 2023; 14:6424. [PMID: 37827998 PMCID: PMC10570356 DOI: 10.1038/s41467-023-42148-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2023] [Accepted: 10/01/2023] [Indexed: 10/14/2023] Open
Abstract
Machine learned interatomic interaction potentials have enabled efficient and accurate molecular simulations of closed systems. However, external fields, which can greatly change the chemical structure and/or reactivity, have been seldom included in current machine learning models. This work proposes a universal field-induced recursively embedded atom neural network (FIREANN) model, which integrates a pseudo field vector-dependent feature into atomic descriptors to represent system-field interactions with rigorous rotational equivariance. This "all-in-one" approach correlates various response properties like dipole moment and polarizability with the field-dependent potential energy in a single model, very suitable for spectroscopic and dynamics simulations in molecular and periodic systems in the presence of electric fields. Especially for periodic systems, we find that FIREANN can overcome the intrinsic multiple-value issue of the polarization by training atomic forces only. These results validate the universality and capability of the FIREANN method for efficient first-principles modeling of complicated systems in strong external fields.
Collapse
Affiliation(s)
- Yaolong Zhang
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui, 230026, China
- École Polytechnique Fédérale de Lausanne, 1015, Lausanne, Switzerland
| | - Bin Jiang
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui, 230026, China.
- Hefei National Laboratory, University of Science and Technology of China, Hefei, 230088, China.
| |
Collapse
|
24
|
Limbu DK, Shakib FA. Real-Time Dynamics and Detailed Balance in Ring Polymer Surface Hopping: The Impact of Frustrated Hops. J Phys Chem Lett 2023; 14:8658-8666. [PMID: 37732811 DOI: 10.1021/acs.jpclett.3c02085] [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/22/2023]
Abstract
Ring polymer surface hopping (RPSH) has been recently introduced as a well-tailored method for incorporating nuclear quantum effects, such as zero-point energy and tunneling, into nonadiabatic molecular dynamics simulations. The practical widespread usage of RPSH demands a comprehensive benchmarking of different reaction regimes and conditions with equal emphasis on demonstrating both the cons and the pros of the method. Here, we investigate the fundamental questions related to the conservation of energy and detailed balance in the context of RPSH. Using Tully's avoided crossing model as well as a 2-state quantum system coupled to a classical bath undergoing Langevin dynamics, we probe the critical problem of the proper treatment of the classically forbidden transitions stemming from the surface hopping algorithm. We show that proper treatment of these frustrated hops is key to the accurate description of real-time dynamics as well as reproducing the correct quantum Boltzmann populations.
Collapse
Affiliation(s)
- Dil K Limbu
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, New Jersey 07102, United States
| | - Farnaz A Shakib
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, New Jersey 07102, United States
| |
Collapse
|
25
|
Cook L, Runeson JE, Richardson JO, Hele TJH. Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics? J Chem Theory Comput 2023; 19:6109-6125. [PMID: 37704193 PMCID: PMC10536990 DOI: 10.1021/acs.jctc.3c00709] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Indexed: 09/15/2023]
Abstract
A common strategy to simulate mixed quantum-classical dynamics is by propagating classical trajectories with mapping variables, often using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian or the related spin-mapping approach. When mapping the quantum subsystem, the coupled dynamics reduce to a set of equations of motion to integrate. Several numerical algorithms have been proposed, but a thorough performance comparison appears to be lacking. Here, we compare three time-propagation algorithms for the MMST Hamiltonian: the Momentum Integral (MInt) (J. Chem. Phys., 2018, 148, 102326), the Split-Liouvillian (SL) (Chem. Phys., 2017, 482, 124-134), and the algorithm in J. Chem. Phys., 2012, 136, 084101 that we refer to as the Degenerate Eigenvalue (DE) algorithm due to the approximation required during derivation. We analyze the accuracy of individual trajectories, correlation functions, energy conservation, symplecticity, Liouville's theorem, and the computational cost. We find that the MInt algorithm is the only rigorously symplectic algorithm. However, comparable accuracy at a lower computational cost can be obtained with the SL algorithm. The approximation implicitly made within the DE algorithm conserves energy poorly, even for small timesteps, and thus leads to slightly different results. These results should guide future mapping-variable simulations.
Collapse
Affiliation(s)
- Lauren
E. Cook
- Department
of Chemistry, University College London, Christopher Ingold Building, London WC1H 0AJ, U.K.
| | - Johan E. Runeson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, Zürich 8093, Switzerland
| | - Jeremy O. Richardson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, Zürich 8093, Switzerland
| | - Timothy J. H. Hele
- Department
of Chemistry, University College London, Christopher Ingold Building, London WC1H 0AJ, U.K.
| |
Collapse
|
26
|
Hasegawa T. Nuclear Quantum Dynamics of Three-Dimensional Condensed-Phase Systems by Constant Uncertainty Molecular Dynamics. J Phys Chem Lett 2023; 14:8043-8049. [PMID: 37651720 DOI: 10.1021/acs.jpclett.3c01254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Constant uncertainty molecular dynamics (CUMD), which was developed to include nuclear quantum effects (NQEs) in molecular dynamics (MD) simulations, is extended to three-dimensional condensed-phase systems. Its applicability was verified via CUMD simulations of bulk water and ice Ih with the q-TIP4P/F potential model. The simulated radial distribution functions, proton momentum distributions, and infrared spectra were compared with those of previous studies in which the simulations were carried out on the basis of the imaginary-time path-integral formalism. We found that CUMD can function as an alternative trajectory-based simulation method for including static and dynamic NQEs in MD simulations of three-dimensional condensed-phase systems.
Collapse
Affiliation(s)
- Taisuke Hasegawa
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| |
Collapse
|
27
|
Khuu T, Schleif T, Mohamed A, Mitra S, Johnson MA, Valdiviezo J, Heindel JP, Head-Gordon T. Intra-cluster Charge Migration upon Hydration of Protonated Formic Acid Revealed by Anharmonic Analysis of Cold Ion Vibrational Spectra. J Phys Chem A 2023; 127:7501-7509. [PMID: 37669457 DOI: 10.1021/acs.jpca.3c03971] [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/07/2023]
Abstract
The rates of many chemical reactions are accelerated when carried out in micron-sized droplets, but the molecular origin of the rate acceleration remains unclear. One example is the condensation reaction of 1,2-diaminobenzene with formic acid to yield benzimidazole. The observed rate enhancements have been rationalized by invoking enhanced acidity at the surface of methanol solvent droplets with low water content to enable protonation of formic acid to generate a cationic species (protonated formic acid or PFA) formed by attachment of a proton to the neutral acid. Because PFA is the key feature in this reaction mechanism, vibrational spectra of cryogenically cooled, microhydrated PFA·(H2O)n=1-6 were acquired to determine how the extent of charge localization depends on the degree of hydration. Analysis of these highly anharmonic spectra with path integral ab initio molecular dynamics simulations reveals the gradual displacement of the excess proton onto the water network in the microhydration regime at low temperatures with n = 3 as the tipping point for intra-cluster proton transfer.
Collapse
Affiliation(s)
- Thien Khuu
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Tim Schleif
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Ahmed Mohamed
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Sayoni Mitra
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Mark A Johnson
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Jesús Valdiviezo
- Pitzer Theory Center, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Joseph P Heindel
- Pitzer Theory Center, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Pitzer Theory Center, Department of Chemistry, University of California, Berkeley, California 94720, United States
| |
Collapse
|
28
|
Litman Y, Lan J, Nagata Y, Wilkins DM. Fully First-Principles Surface Spectroscopy with Machine Learning. J Phys Chem Lett 2023; 14:8175-8182. [PMID: 37671886 PMCID: PMC10510433 DOI: 10.1021/acs.jpclett.3c01989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Accepted: 08/29/2023] [Indexed: 09/07/2023]
Abstract
Our current understanding of the structure and dynamics of aqueous interfaces at the molecular level has grown substantially due to the continuous development of surface-specific spectroscopies, such as vibrational sum-frequency generation (VSFG). As in other vibrational spectroscopies, we must turn to atomistic simulations to extract all of the information encoded in the VSFG spectra. The high computational cost associated with existing methods means that they have limitations in representing systems with complex electronic structure or in achieving statistical convergence. In this work, we combine high-dimensional neural network interatomic potentials and symmetry-adapted Gaussian process regression to overcome these constraints. We show that it is possible to model VSFG signals with fully ab initio accuracy using machine learning and illustrate the versatility of our approach on the water/air interface. Our strategy allows us to identify the main sources of theoretical inaccuracy and establish a clear pathway toward the modeling of surface-sensitive spectroscopy of complex interfaces.
Collapse
Affiliation(s)
- Yair Litman
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K.
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Jinggang Lan
- Department
of Chemistry, New York University, New York, New York 10003, United States
- Simons
Center for Computational Physical Chemistry at New York University, New York, New York 10003, United States
| | - Yuki Nagata
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - David M. Wilkins
- Centre
for Quantum Materials and Technologies School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| |
Collapse
|
29
|
Atsango AO, Morawietz T, Marsalek O, Markland TE. Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations. J Chem Phys 2023; 159:074101. [PMID: 37581418 DOI: 10.1063/5.0162066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Accepted: 07/31/2023] [Indexed: 08/16/2023] Open
Abstract
The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions.
Collapse
Affiliation(s)
- Austin O Atsango
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Tobias Morawietz
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Ondrej Marsalek
- Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
30
|
Chen MS, Mao Y, Snider A, Gupta P, Montoya-Castillo A, Zuehlsdorff TJ, Isborn CM, Markland TE. Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure. J Phys Chem Lett 2023; 14:6610-6619. [PMID: 37459252 DOI: 10.1021/acs.jpclett.3c01444] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. We show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Andrew Snider
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Prachi Gupta
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Andrés Montoya-Castillo
- Department of Chemistry, University of Colorado, Boulder, Boulder, Colorado 80309, United States
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Christine M Isborn
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
31
|
Lieberherr AZ, Furniss STE, Lawrence JE, Manolopoulos DE. Vibrational strong coupling in liquid water from cavity molecular dynamics. J Chem Phys 2023; 158:234106. [PMID: 37326163 DOI: 10.1063/5.0156808] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Accepted: 05/30/2023] [Indexed: 06/17/2023] Open
Abstract
We assess the cavity molecular dynamics method for the calculation of vibrational polariton spectra using liquid water as a specific example. We begin by disputing a recent suggestion that nuclear quantum effects may lead to a broadening of polariton bands, finding instead that they merely result in anharmonic red shifts in the polariton frequencies. We go on to show that our simulated cavity spectra can be reproduced to graphical accuracy with a harmonic model that uses just the cavity-free spectrum and the geometry of the cavity as input. We end by showing that this harmonic model can be combined with the experimental cavity-free spectrum to give results in good agreement with optical cavity measurements. Since the input to our harmonic model is equivalent to the input to the transfer matrix method of applied optics, we conclude that cavity molecular dynamics cannot provide any more insight into the effect of vibrational strong coupling on the absorption spectrum than this transfer matrix method, which is already widely used by experimentalists to corroborate their cavity results.
Collapse
Affiliation(s)
- Annina Z Lieberherr
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Seth T E Furniss
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Joseph E Lawrence
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - David E Manolopoulos
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
32
|
Wang Y, Chen Z, Yang Y. Calculating Vibrational Excited State Absorptions with Excited State Constrained Minimized Energy Surfaces. J Phys Chem A 2023. [PMID: 37335973 DOI: 10.1021/acs.jpca.3c01420] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
The modeling and interpretation of vibrational spectra are crucial for studying reaction dynamics using vibrational spectroscopy. Most prior theoretical developments focused on describing fundamental vibrational transitions while fewer developments focused on vibrational excited state absorptions. In this study, we present a new method that uses excited state constrained minimized energy surfaces (CMESs) to describe vibrational excited state absorptions. The excited state CMESs are obtained similarly to the previous ground state CMES development in our group but with additional wave function orthogonality constraints. Using a series of model systems, including the harmonic oscillator, Morse potential, double-well potential, quartic potential, and two-dimensional anharmonic potential, we demonstrate that this new procedure provides good estimations of the transition frequencies for vibrational excited state absorptions. These results are significantly better than those obtained from harmonic approximations using conventional potential energy surfaces, demonstrating the promise of excited state CMES-based methods for calculating vibrational excited state absorptions in real systems.
Collapse
Affiliation(s)
- Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
33
|
Shepherd S, Tribello GA, Wilkins DM. A fully quantum-mechanical treatment for kaolinite. J Chem Phys 2023; 158:2892274. [PMID: 37220200 DOI: 10.1063/5.0152361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Accepted: 05/03/2023] [Indexed: 05/25/2023] Open
Abstract
Neural network potentials for kaolinite minerals have been fitted to data extracted from density functional theory calculations that were performed using the revPBE + D3 and revPBE + vdW functionals. These potentials have then been used to calculate the static and dynamic properties of the mineral. We show that revPBE + vdW is better at reproducing the static properties. However, revPBE + D3 does a better job of reproducing the experimental IR spectrum. We also consider what happens to these properties when a fully quantum treatment of the nuclei is employed. We find that nuclear quantum effects (NQEs) do not make a substantial difference to the static properties. However, when NQEs are included, the dynamic properties of the material change substantially.
Collapse
Affiliation(s)
- Sam Shepherd
- Centre for Quantum Materials and Technologies, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| | - Gareth A Tribello
- Centre for Quantum Materials and Technologies, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| | - David M Wilkins
- Centre for Quantum Materials and Technologies, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| |
Collapse
|
34
|
Videla PE, Batista VS. Matsubara dynamics approximation for generalized multi-time correlation functions. J Chem Phys 2023; 158:2889027. [PMID: 37154285 DOI: 10.1063/5.0146654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/21/2023] [Indexed: 05/10/2023] Open
Abstract
We introduce a semi-classical approximation for calculating generalized multi-time correlation functions based on Matsubara dynamics, a classical dynamics approach that conserves the quantum Boltzmann distribution. This method is exact for the zero time and harmonic limits and reduces to classical dynamics when only one Matsubara mode is considered (i.e., the centroid). Generalized multi-time correlation functions can be expressed as canonical phase-space integrals, involving classically evolved observables coupled through Poisson brackets in a smooth Matsubara space. Numerical tests on a simple potential show that the Matsubara approximation exhibits better agreement with exact results than classical dynamics, providing a bridge between the purely quantum and classical descriptions of multi-time correlation functions. Despite the phase problem that prevents practical applications of Matsubara dynamics, the reported work provides a benchmark theory for the future development of quantum-Boltzmann-preserving semi-classical approximations for studies of chemical dynamics in condensed phase systems.
Collapse
Affiliation(s)
- Pablo E Videla
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
35
|
Schütt KT, Hessmann SSP, Gebauer NWA, Lederer J, Gastegger M. SchNetPack 2.0: A neural network toolbox for atomistic machine learning. J Chem Phys 2023; 158:144801. [PMID: 37061495 DOI: 10.1063/5.0138367] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
SchNetPack is a versatile neural network toolbox that addresses both the requirements of method development and the application of atomistic machine learning. Version 2.0 comes with an improved data pipeline, modules for equivariant neural networks, and a PyTorch implementation of molecular dynamics. An optional integration with PyTorch Lightning and the Hydra configuration framework powers a flexible command-line interface. This makes SchNetPack 2.0 easily extendable with a custom code and ready for complex training tasks, such as the generation of 3D molecular structures.
Collapse
Affiliation(s)
- Kristof T Schütt
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| | | | - Niklas W A Gebauer
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| | - Jonas Lederer
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| | - Michael Gastegger
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| |
Collapse
|
36
|
Begušić T, Blake GA. Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality. Nat Commun 2023; 14:1950. [PMID: 37029146 PMCID: PMC10082090 DOI: 10.1038/s41467-023-37667-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 03/22/2023] [Indexed: 04/09/2023] Open
Abstract
Two-dimensional spectroscopic techniques combining terahertz (THz), infrared (IR), and visible pulses offer a wealth of information about coupling among vibrational modes in molecular liquids, thus providing a promising probe of their local structure. However, the capabilities of these spectroscopies are still largely unexplored due to experimental limitations and inherently weak nonlinear signals. Here, through a combination of equilibrium-nonequilibrium molecular dynamics (MD) and a tailored spectrum decomposition scheme, we identify a relationship between the tetrahedral order of liquid water and its two-dimensional IR-IR-Raman (IIR) spectrum. The structure-spectrum relationship can explain the temperature dependence of the spectral features corresponding to the anharmonic coupling between low-frequency intermolecular and high-frequency intramolecular vibrational modes of water. In light of these results, we propose new experiments and discuss the implications for the study of tetrahedrality of liquid water.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| | - Geoffrey A Blake
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, 91125, USA.
| |
Collapse
|
37
|
Prada A, Pós ES, Althorpe SC. Comparison of Matsubara dynamics with exact quantum dynamics for an oscillator coupled to a dissipative bath. J Chem Phys 2023; 158:114106. [PMID: 36948794 DOI: 10.1063/5.0138250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023] Open
Abstract
We report the first numerical calculations in which converged Matsubara dynamics is compared directly with exact quantum dynamics with no artificial damping of the time-correlation functions (TCFs). The system treated is a Morse oscillator coupled to a harmonic bath. We show that, when the system-bath coupling is sufficiently strong, the Matsubara calculations can be converged by explicitly including up to M = 200 Matsubara modes, with the remaining modes included as a harmonic "tail" correction. The resulting Matsubara TCFs are in near-perfect agreement with the exact quantum TCFs, for non-linear as well as linear operators, at a temperature at which the TCFs are dominated by quantum thermal fluctuations. These results provide compelling evidence that incoherent classical dynamics can arise in the condensed phase at temperatures at which the statistics are dominated by quantum (Boltzmann) effects, as a result of smoothing of imaginary-time Feynman paths. The techniques developed here may also lead to efficient methods for benchmarking system-bath dynamics in the overdamped regime.
Collapse
Affiliation(s)
- Adam Prada
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Eszter S Pós
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
38
|
Plé T, Mauger N, Adjoua O, Inizan TJ, Lagardère L, Huppert S, Piquemal JP. Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: From Force Fields to Machine Learning Potentials. J Chem Theory Comput 2023; 19:1432-1445. [PMID: 36856658 DOI: 10.1021/acs.jctc.2c01233] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2023]
Abstract
We report the implementation of a multi-CPU and multi-GPU massively parallel platform dedicated to the explicit inclusion of nuclear quantum effects (NQEs) in the Tinker-HP molecular dynamics (MD) package. The platform, denoted Quantum-HP, exploits two simulation strategies: the Ring-Polymer Molecular Dynamics (RPMD) that provides exact structural properties at the cost of a MD simulation in an extended space of multiple replicas and the adaptive Quantum Thermal Bath (adQTB) that imposes the quantum distribution of energy on a classical system via a generalized Langevin thermostat and provides computationally affordable and accurate (though approximate) NQEs. We discuss some implementation details, efficient numerical schemes, and parallelization strategies and quickly review the GPU acceleration of our code. Our implementation allows an efficient inclusion of NQEs in MD simulations for very large systems, as demonstrated by scaling tests on water boxes with more than 200,000 atoms (simulated using the AMOEBA polarizable force field). We test the compatibility of the approach with Tinker-HP's recently introduced Deep-HP machine learning potentials module by computing water properties using the DeePMD potential with adQTB thermostatting. Finally, we show that the platform is also compatible with the alchemical free energy estimation capabilities of Tinker-HP and fast enough to perform simulations. Therefore, we study how NQEs affect the hydration free energy of small molecules solvated with the recently developed Q-AMOEBA water force field. Overall, the Quantum-HP platform allows users to perform routine quantum MD simulations of large condensed-phase systems and will help to shed new light on the quantum nature of important interactions in biological matter.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | | | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Simon Huppert
- Institut des Nanosciences de Paris (INSP), CNRS UMR 7588, and Sorbonne Université, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France.,Institut Universitaire de France, 75005 Paris, France.,Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
39
|
Kurapothula PJ, Shepherd S, Wilkins DM. Competing Nuclear Quantum Effects and Hydrogen-Bond Jumps in Hydrated Kaolinite. J Phys Chem Lett 2023; 14:1542-1547. [PMID: 36745462 PMCID: PMC9940297 DOI: 10.1021/acs.jpclett.2c03896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 02/01/2023] [Indexed: 06/18/2023]
Abstract
Recent work has shown that the dynamics of hydrogen bonds in pure clays are affected by nuclear quantum fluctuations, with different effects for the hydrogen bonds holding different layers of the clay together and for those within the same layer. At the clay-water interface there is an even wider range of types of hydrogen bond, suggesting that the quantum effects may be yet more varied. We apply classical and thermostated ring polymer molecular dynamics simulations to show that nuclear quantum effects accelerate hydrogen-bond dynamics to varying degrees. By interpreting the results in terms of the extended jump model of hydrogen-bond switching, we can understand the origins of these effects in terms of changes in the quantum kinetic energy of hydrogen atoms during an exchange. We also show that the extended jump mechanism is applicable not only to the hydrogen bonds involving water, but also those internal to the clay.
Collapse
|
40
|
Chen Z, Yang Y. Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface. J Phys Chem Lett 2023; 14:279-286. [PMID: 36595586 DOI: 10.1021/acs.jpclett.2c02905] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The accurate incorporation of nuclear quantum effects in large-scale molecular dynamics (MD) simulations remains a significant challenge. Recently, we combined constrained nuclear-electronic orbital (CNEO) theory with classical MD and obtained a new approach (CNEO-MD) that can accurately and efficiently incorporate nuclear quantum effects into classical simulations. In this Letter, we provide the theoretical foundation for CNEO-MD by developing an alternative formulation of the equations of motion for MD. In this new formulation, the expectation values of quantum nuclear positions evolve classically on an effective energy surface that is obtained from a constrained energy minimization procedure when solving for the quantum nuclear wave function, thus enabling the incorporation of nuclear quantum effects in classical MD simulations. For comparison with other existing approaches, we examined a series of model systems and found that this new MD approach is significantly more accurate than the conventional way of performing classical MD and generally outperforms centroid MD and ring-polymer MD in describing vibrations in these model systems.
Collapse
Affiliation(s)
- Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| |
Collapse
|
41
|
Ikeda T. First principles molecular dynamics study of proton disorder in C1′ phase of H2 hydrate. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2022.140252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
42
|
Brünig FN, Hillmann P, Kim WK, Daldrop JO, Netz RR. Proton-transfer spectroscopy beyond the normal-mode scenario. J Chem Phys 2022; 157:174116. [DOI: 10.1063/5.0116686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single [Formula: see text] cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.
Collapse
Affiliation(s)
- Florian N. Brünig
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Paul Hillmann
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Won Kyu Kim
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, Republic of Korea
| | - Jan O. Daldrop
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Roland R. Netz
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
43
|
Mauger N, Plé T, Lagardère L, Huppert S, Piquemal JP. Improving Condensed-Phase Water Dynamics with Explicit Nuclear Quantum Effects: The Polarizable Q-AMOEBA Force Field. J Phys Chem B 2022; 126:8813-8826. [PMID: 36270033 DOI: 10.1021/acs.jpcb.2c04454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We introduce a new parametrization of the AMOEBA polarizable force field for water denoted Q-AMOEBA, for use in simulations that explicitly account for nuclear quantum effects (NQEs). This study is made possible thanks to the recently introduced adaptive Quantum Thermal Bath (adQTB) simulation technique which computational cost is comparable to classical molecular dynamics. The flexible Q-AMOEBA model conserves the initial AMOEBA functional form, with an intermolecular potential including an atomic multipole description of electrostatic interactions (up to quadrupole), a polarization contribution based on the Thole interaction model and a buffered 14-7 potential to model van der Waals interactions. It has been obtained by using a ForceBalance fitting strategy including high-level quantum chemistry reference energies and selected condensed-phase properties targets. The final Q-AMOEBA model is shown to accurately reproduce both gas-phase and condensed-phase properties, notably improving the original AMOEBA water model. This development allows the fine study of NQEs on water liquid phase properties such as the average H-O-H angle compared to its gas-phase equilibrium value, isotope effects, and so on. Q-AMOEBA also provides improved infrared spectroscopy prediction capabilities compared to AMOEBA03. Overall, we show that the impact of NQEs depends on the underlying model functional form and on the associated strength of hydrogen bonds. Since adQTB simulations can be performed at near classical computational cost using the Tinker-HP package, Q-AMOEBA can be extended to organic molecules, proteins, and nucleic acids opening the possibility for the large-scale study of the importance of NQEs in biophysics.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Thomas Plé
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| |
Collapse
|
44
|
Loose T, Sahrmann PG, Voth GA. Centroid Molecular Dynamics Can Be Greatly Accelerated through Neural Network Learned Centroid Forces Derived from Path Integral Molecular Dynamics. J Chem Theory Comput 2022; 18:5856-5863. [PMID: 36103576 PMCID: PMC9558744 DOI: 10.1021/acs.jctc.2c00706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Indexed: 11/29/2022]
Abstract
For nearly the past 30 years, centroid molecular dynamics (CMD) has proven to be a viable classical-like phase space formulation for the calculation of quantum dynamical properties. However, calculation of the centroid effective force remains a significant computational cost and limits the ability of CMD to be an efficient approach to study condensed phase quantum dynamics. In this paper, we introduce a neural network-based methodology for first learning the centroid effective force from path integral molecular dynamics data, which is subsequently used as an effective force field to evolve the centroids directly with the CMD algorithm. This method, called machine-learned centroid molecular dynamics (ML-CMD), is faster and far less costly than both standard "on the fly" CMD and ring polymer molecular dynamics (RPMD). The training aspect of ML-CMD is also straightforwardly implemented utilizing the DeepMD software kit. ML-CMD is then applied to two model systems to illustrate the approach: liquid para-hydrogen and water. The results show comparable accuracy to both CMD and RPMD in the estimation of quantum dynamical properties, including the self-diffusion constant and velocity time correlation function, but with significantly reduced overall computational cost.
Collapse
Affiliation(s)
- Timothy
D. Loose
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, James Franck Institute,
and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Patrick G. Sahrmann
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, James Franck Institute,
and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, James Franck Institute,
and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
45
|
Shiga M. Path integral Brownian chain molecular dynamics: A simple approximation of quantum vibrational dynamics. J Comput Chem 2022; 43:1864-1879. [PMID: 36094104 DOI: 10.1002/jcc.26989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2022] [Revised: 08/06/2022] [Accepted: 08/09/2022] [Indexed: 11/05/2022]
Abstract
An approximate approach to quantum vibrational dynamics, "Brownian chain molecular dynamics (BCMD)," is proposed to alleviate the chain resonance and curvature problems in the imaginary time-based path integral (PI) simulation. Here the non-centroid velocity is randomized at each step when solving the equation of motion of path integral molecular dynamics. This leads to a combination of the Newton equation and the overdamped Langevin equation for the centroid and non-centroid variables, respectively. BCMD shares the basic properties of other PI approaches such as centroid and ring polymer molecular dynamics: It gives the correct Kubo-transformed correlation function at short times, conserves the time symmetry, has the correct high-temperature/classical limits, gives exactly the position and velocity autocorrelations of harmonic oscillator systems, and does not have the zero-point leakage problem. Numerical tests were done on simple molecular models and liquid water. On-the-fly ab initio BCMD simulations were performed for the protonated water cluster, H 5 O 2 + $$ {\mathrm{H}}_5{\mathrm{O}}_2^{+} $$ , and its isotopologue, D 5 O 2 + $$ {\mathrm{D}}_5{\mathrm{O}}_2^{+} $$ .
Collapse
Affiliation(s)
- Motoyuki Shiga
- Center for Computational Science and e-Systems, Japan Atomic Energy Agency, Chiba, Japan
| |
Collapse
|
46
|
Ryu WH, Voth GA. Coarse-Graining of Imaginary Time Feynman Path Integrals: Inclusion of Intramolecular Interactions and Bottom-up Force-Matching. J Phys Chem A 2022; 126:6004-6019. [PMID: 36007243 PMCID: PMC9466601 DOI: 10.1021/acs.jpca.2c04349] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Revised: 08/16/2022] [Indexed: 11/29/2022]
Abstract
Feynman's imaginary time path integral formalism of quantum statistical mechanics and the corresponding quantum-classical isomorphism provide a tangible way of incorporating nuclear quantum effect (NQE) in the simulation of condensed matter systems using well-developed classical simulation techniques. Our previous work has presented the many-body coarse-graining of path integral (CG-PI) theory that builds an isomorphism between the quantum partition function of N distinguishable particles and the classical partition function of 2N pseudoparticles. In this present work, we develop a generalized version of the many-body CG-PI theory that incorporates many-body interactions in the force field. Based on the new derivation, we provide a numerical CG-PI (n-CG-PI) modeling strategy parametrized from the underlying path integral molecular dynamics (PIMD) trajectories using force matching and Boltzmann inversion. The n-CG-PI models for two liquid systems are shown to capture well both the intramolecular and intermolecular structural correlations of the reference PIMD simulations. The generalized derivation of the many-body CG-PI theory and the n-CG-PI model presented in this work extend the scope of the CG-PI formalism by generalizing the previously limited theory to incorporate force fields of realistic molecular systems.
Collapse
Affiliation(s)
- Won Hee Ryu
- Department of Chemistry,
James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry,
James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
47
|
Jung KA, Markland TE. 2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions. J Chem Phys 2022; 157:094111. [DOI: 10.1063/5.0107087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The third-order response lies at the heart of simulating and interpreting nonlinear spectroscopies ranging from two dimensional infrared (2D-IR) to 2D electronic (2D-ES), and 2D sum frequency generation (2D-SFG). The extra time and frequency dimensions in these spectroscopies provides access to rich information on the electronic and vibrational states present, the coupling between them, and the resulting rates at which they exchange energy that are obscured in linear spectroscopy, particularly for condensed phase systems that usually contain many overlapping features. While the exact quantum expression for the third-order response is well established it is incompatible with the methods that are practical for calculating the atomistic dynamics of large condensed phase systems. These methods, which include both classical mechanics and quantum dynamics methods that retain quantum statistical properties while obeying the symmetries of classical dynamics, such as LSC-IVR, Centroid Molecular Dynamics (CMD) and Ring Polymer Molecular Dynamics (RPMD) naturally provide short-time approximations to the multi-time symmetrized Kubo transformed correlation function. Here, we show how the third-order response can be formulated in terms of equilibrium symmetrized Kubo transformed correlation functions. We demonstrate the utility and accuracy of our approach by showing how it can be used to obtain the third-order response of a series of model systems using both classical dynamics and RPMD. In particular, we show that this approach captures features such as anharmonically induced vertical splittings and peak shifts while providing a physically transparent framework for understanding multidimensional spectroscopies.
Collapse
|
48
|
Chu W, Tan S, Zheng Q, Fang W, Feng Y, Prezhdo OV, Wang B, Li XZ, Zhao J. Ultrafast charge transfer coupled to quantum proton motion at molecule/metal oxide interface. SCIENCE ADVANCES 2022; 8:eabo2675. [PMID: 35714193 DOI: 10.1126/sciadv.abo2675] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Understanding how the nuclear quantum effects (NQEs) in the hydrogen bond (H-bond) network influence the photoexcited charge transfer at semiconductor/molecule interface is a challenging problem. By combining two kinds of emerging molecular dynamics methods at the ab initio level, the path integral-based molecular dynamics and time-dependent nonadiabatic molecular dynamics, and choosing CH3OH/TiO2 as a prototypical system to study, we find that the quantum proton motion in the H-bond network is strongly coupled with the ultrafast photoexcited charge dynamics at the interface. The hole trapping ability of the adsorbed methanol molecule is notably enhanced by the NQEs, and thus, it behaves as a hole scavenger on titanium dioxide. The critical role of the H-bond network is confirmed by in situ scanning tunneling microscope measurements with ultraviolet light illumination. It is concluded the quantum proton motion in the H-bond network plays a critical role in influencing the energy conversion efficiency based on photoexcitation.
Collapse
Affiliation(s)
- Weibin Chu
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, CA 90089, USA
- Key Laboratory of Computational Physical Sciences (Ministry of Education), Institute of Computational Physical Sciences, Fudan University, Shanghai 200433, People's Republic of China
| | - Shijing Tan
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Qijing Zheng
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Wei Fang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
- Department of Chemistry, Fudan University, Shanghai 200438, People's Republic of China
- Laboratory of Physical Chemistry, ETH Zurich, CH-8093 Zürich, Switzerland
| | - Yexin Feng
- School of Physics and Electronics, Hunan University, Changsha 410082, People's Republic of China
| | - Oleg V Prezhdo
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, CA 90089, USA
| | - Bing Wang
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Xin-Zheng Li
- Interdisciplinary Institute of Light-Element Quantum Materials, Research Center for Light-Element Advanced Materials, State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Frontier Science Center for Nano-optoelectronics and School of Physics, Peking University, Beijing 100871, People's Republic of China
- Peking University Yangtze Delta Institute of Optoelectronics, Nantong, Jiangsu 226010, People's Republic of China
| | - Jin Zhao
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
- Hefei National Laboratory, Hefei 230088, People's Republic of China
| |
Collapse
|
49
|
Töpfer K, Upadhyay M, Meuwly M. Quantitative molecular simulations. Phys Chem Chem Phys 2022; 24:12767-12786. [PMID: 35593769 PMCID: PMC9158373 DOI: 10.1039/d2cp01211a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 04/30/2022] [Indexed: 11/21/2022]
Abstract
All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.
Collapse
Affiliation(s)
- Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| |
Collapse
|
50
|
Eltareb A, Lopez GE, Giovambattista N. Nuclear quantum effects on the dynamics and glass behavior of a monatomic liquid with two liquid states. J Chem Phys 2022; 156:204502. [PMID: 35649856 PMCID: PMC9132595 DOI: 10.1063/5.0087680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Accepted: 05/02/2022] [Indexed: 11/14/2022] Open
Abstract
We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid-liquid phase transition and liquid-liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., Tg LDL(P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., Tg HDL(P) increases upon compression. (iii) NQE shift both Tg LDL(P) and Tg HDL(P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA-HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA-HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA-HDA transformation.
Collapse
|