1
|
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
|
2
|
Beckmann R, Topolnicki R, Marx D. Deciphering the Impact of Helium Tagging on Flexible Molecules: Probing Microsolvation Effects of Protonated Acetylene by Quantum Configurational Entropy. J Phys Chem A 2023; 127:2460-2471. [PMID: 36917575 DOI: 10.1021/acs.jpca.2c08967] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Helium, the lightest and most weakly interacting noble gas, is well-known for its unsurpassed chemical inertness. In many applications of helium in experimental techniques, such as tagging, messenger, or nanodroplet isolation action spectroscopy of molecules or complexes, it is assumed that the interaction of helium with the respective species, and thus the resulting interaction-induced perturbation, is small enough not to affect their structure and dynamics. Here, we probe the impact of one up to many attached helium atoms on protonated acetylene─an important nonclassical carbocation subject to three-center two-electron bonding in its ground state structure─using highly accurate interaction potentials in conjunction with entropy-based higher-order nonlinear correlation analysis. In particular, using neural network potentials at CCSD(T) accuracy, we disclose the specific structural perturbations due to the tagging of C2H3+ with up to 20 He atoms at a temperature of 1 K. Analysis reveals that microsolvation by helium influences the structure of C2H3+ noticeably, while our investigation of the quantum configurational information entropy additionally shows that correlations between individual orientational degrees of freedom are affected as a function of cluster size. In particular, it is found that the most probable bridge-like structure of the ro-vibrational quantum ground state of C2H3+, which is nonplanar and trans-bent in contrast to the perfectly planar equilibrium structure, becomes increasingly more localized upon adding helium atoms. The remarkably nonlinear behavior of the angular correlations as a function of cluster size is traced back to the buildup of the quantum microsolvation shell that enhances anisotropy up to NHe = 6 while more and more isotropic solvation takes over beyond six. Our approach is general and thus sets the stage to investigate the salient effects on the structure of flexible molecules due to tagging beyond the specific case.
Collapse
Affiliation(s)
- Richard Beckmann
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Rafal Topolnicki
- Dioscuri Center in Topological Data Analysis, Institute of Mathematics, Polish Academy of Sciences, Śniadeckich 8, Warsaw 00-656, Poland
- Institute of Experimental Physics, University of Wrocław, 50-204 Wrocław, Poland
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
3
|
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
|
4
|
Daru J, Forbert H, Behler J, Marx D. Coupled Cluster Molecular Dynamics of Condensed Phase Systems Enabled by Machine Learning Potentials: Liquid Water Benchmark. PHYSICAL REVIEW LETTERS 2022; 129:226001. [PMID: 36493459 DOI: 10.1103/physrevlett.129.226001] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/17/2022] [Revised: 09/05/2022] [Accepted: 10/05/2022] [Indexed: 06/17/2023]
Abstract
Coupled cluster theory is a general and systematic electronic structure method, but in particular the highly accurate "gold standard" coupled cluster singles, doubles and perturbative triples, CCSD(T), can only be applied to small systems. To overcome this limitation, we introduce a framework to transfer CCSD(T) accuracy of finite molecular clusters to extended condensed phase systems using a high-dimensional neural network potential. This approach, which is automated, allows one to perform high-quality coupled cluster molecular dynamics, CCMD, as we demonstrate for liquid water including nuclear quantum effects. The machine learning strategy is very efficient, generic, can be systematically improved, and is applicable to a variety of complex systems.
Collapse
Affiliation(s)
- János Daru
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
5
|
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
|
6
|
Impact of nuclear quantum effects on the $$^{13}$$C nuclear shielding of linear carboxylates in water. Theor Chem Acc 2022. [DOI: 10.1007/s00214-022-02894-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
7
|
Simulation of Nuclear Quantum Effects in Condensed Matter Systems via Quantum Baths. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12094756] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/10/2022]
Abstract
This paper reviews methods that aim at simulating nuclear quantum effects (NQEs) using generalized thermal baths. Generalized (or quantum) baths simulate statistical quantum features, and in particular zero-point energy effects, through non-Markovian stochastic dynamics. They make use of generalized Langevin Equations (GLEs), in which the quantum Bose–Einstein energy distribution is enforced by tuning the random and friction forces, while the system degrees of freedom remain classical. Although these baths have been formally justified only for harmonic oscillators, they perform well for several systems, while keeping the cost of the simulations comparable to the classical ones. We review the formal properties and main characteristics of classical and quantum GLEs, in relation with the fluctuation–dissipation theorems. Then, we describe the quantum thermostat and quantum thermal bath, the two generalized baths currently most used, providing several examples of applications for condensed matter systems, including the calculation of vibrational spectra. The most important drawback of these methods, zero-point energy leakage, is discussed in detail with the help of model systems, and a recently proposed scheme to monitor and mitigate or eliminate it—the adaptive quantum thermal bath—is summarised. This approach considerably extends the domain of application of generalized baths, leading, for instance, to the successful simulation of liquid water, where a subtle interplay of NQEs is at play. The paper concludes by overviewing further development opportunities and open challenges of generalized baths.
Collapse
|
8
|
Depondt P, Huppert S, Finocchi F. The quantum taste of hydrogen. EPJ WEB OF CONFERENCES 2022. [DOI: 10.1051/epjconf/202226301014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Electronic properties of materials are dominated by quantum effects, but nuclei, being much heavier, are usually treated as classical particles. This approximation, although tremendously convenient, is not always valid, even in close to ambient pressure and temperature conditions, especially when light nuclei such as hydrogen are involved. Zero point energy and proton tunneling can be relevant. Isotopic effects, obtained by replacing hydrogen with deuterium, are observed experimentally and are a clear indication of Nuclear Quantum Effects (NQE) since mean values obtained through classical statistical physics do not depend on mass. Introducing NQEs into simulations at an acceptable computational cost raises fundamental questions and yields subtle and unexpected results.
Collapse
|
9
|
Mauger N, Plé T, Lagardère L, Bonella S, Mangaud É, Piquemal JP, Huppert S. Nuclear Quantum Effects in Liquid Water at Near Classical Computational Cost Using the Adaptive Quantum Thermal Bath. J Phys Chem Lett 2021; 12:8285-8291. [PMID: 34427440 DOI: 10.1021/acs.jpclett.1c01722] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQEs) in molecular simulations: the adaptive quantum thermal bath (adQTB). In this method, zero-point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation-dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to that of classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Thomas Plé
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Sara Bonella
- CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne, Batochimie, Avenue Forel 2, 1015 Lausanne, Switzerland
| | - Étienne Mangaud
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, 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
| | - Simon Huppert
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| |
Collapse
|
10
|
Schran C, Brieuc F, Marx D. Transferability of machine learning potentials: Protonated water neural network potential applied to the protonated water hexamer. J Chem Phys 2021; 154:051101. [DOI: 10.1063/5.0035438] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - 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
|
11
|
Topolnicki R, Brieuc F, Schran C, Marx D. Deciphering High-Order Structural Correlations within Fluxional Molecules from Classical and Quantum Configurational Entropy. J Chem Theory Comput 2020; 16:6785-6794. [DOI: 10.1021/acs.jctc.0c00642] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Rafał Topolnicki
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum 44780, Germany
- Institute of Experimental Physics, University of Wrocław, Wrocław 50-204, Poland
| | - Fabien Brieuc
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum 44780, Germany
| | - Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum 44780, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum 44780, Germany
| |
Collapse
|
12
|
Brieuc F, Schran C, Uhl F, Forbert H, Marx D. Converged quantum simulations of reactive solutes in superfluid helium: The Bochum perspective. J Chem Phys 2020; 152:210901. [DOI: 10.1063/5.0008309] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Fabien Brieuc
- 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
| | - Felix Uhl
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
13
|
Kapil V, Wilkins DM, Lan J, Ceriotti M. Inexpensive modeling of quantum dynamics using path integral generalized Langevin equation thermostats. J Chem Phys 2020; 152:124104. [PMID: 32241150 DOI: 10.1063/1.5141950] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The properties of molecules and materials containing light nuclei are affected by their quantum mechanical nature. Accurate modeling of these quantum nuclear effects requires computationally demanding path integral techniques. Considerable success has been achieved in reducing the cost of such simulations by using generalized Langevin dynamics to induce frequency-dependent fluctuations. Path integral generalized Langevin equation methods, however, have this far been limited to the study of static, thermodynamic properties due to the large perturbation to the system's dynamics induced by the aggressive thermostatting. Here, we introduce a post-processing scheme, based on analytical estimates of the dynamical perturbation induced by the generalized Langevin dynamics, which makes it possible to recover meaningful time correlation properties from a thermostatted trajectory. We show that this approach yields spectroscopic observables for model and realistic systems that have an accuracy comparable to much more demanding approximate quantum dynamics techniques based on full path integral simulations.
Collapse
Affiliation(s)
- Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - David M Wilkins
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Jinggang Lan
- Department of Chemistry, University of Zürich, Zürich, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| |
Collapse
|
14
|
Drewsen M, Imparato A. Quantum duets working as autonomous thermal motors. Phys Rev E 2019; 100:042138. [PMID: 31770990 DOI: 10.1103/physreve.100.042138] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Indexed: 11/07/2022]
Abstract
We study the dynamic properties of a thermal autonomous machine made up of two quantum Brownian particles, each of which is in contact with an environment at different temperature and moves on a periodic sinusoidal track. When such tracks are shifted, the center of mass of the system exhibits a nonvanishing velocity, for which we provide an exact expression in the limit of small track undulations. We discuss the role of the broken spatial symmetry in the emergence of directed motion in thermal machines. We then consider the case in which external deterministic forces are applied to the system, and we characterize its steady-state velocity. If the applied external force opposes the system motion, work can be extracted from such a steady-state thermal machine, without any external cyclic protocol. When the two particles are not interacting, our results reduce to those of Fisher and Zwerger [Phys. Rev. B 32, 6190 (1985)PRBMDO0163-182910.1103/PhysRevB.32.6190] and Aslangul, Pottier, and Saint-James [J. Phys. France 48, 1093 (1987)JOPQAG0302-073810.1051/jphys:019870048070109300] for a single particle moving in a periodic tilted potential. We finally use our results for the motor velocity to check the validity of the quantum molecular dynamics algorithm in the nonlinear, nonequilibrium regime.
Collapse
Affiliation(s)
- Michael Drewsen
- Department of Physics and Astronomy, University of Aarhus Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark
| | - Alberto Imparato
- Department of Physics and Astronomy, University of Aarhus Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark
| |
Collapse
|
15
|
Mangaud E, Huppert S, Plé T, Depondt P, Bonella S, Finocchi F. The Fluctuation–Dissipation Theorem as a Diagnosis and Cure for Zero-Point Energy Leakage in Quantum Thermal Bath Simulations. J Chem Theory Comput 2019; 15:2863-2880. [DOI: 10.1021/acs.jctc.8b01164] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Etienne Mangaud
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Thomas Plé
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Philippe Depondt
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Sara Bonella
- CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne, Batochimie, Avenue Forel 2, 1015 Lausanne, Switzerland
| | - Fabio Finocchi
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| |
Collapse
|
16
|
Dammak H, Brieuc F, Geneste G, Torrent M, Hayoun M. Isotope effect on hydrogen bond symmetrization in hydrogen and deuterium fluoride crystals by molecular dynamics simulation. Phys Chem Chem Phys 2019; 21:3211-3217. [PMID: 30681084 DOI: 10.1039/c8cp06949b] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The isotope effect on the collective proton/deuteron transfer in hydrogen and deuterium fluoride crystals has been investigated at 100 K by ab initio quantum-thermal-bath path-integral molecular dynamics (QTB-PIMD) simulation. The deuterons within a planar zigzag chain of the orthorhombic structure simultaneously flip between covalent and hydrogen bonds due to the barrier crossing through tunnelling. The height of the corresponding static barrier normalized for one deuteron is 29.2 meV. In the HF crystal, all the protons are located at the center of the heavy-atom distance. This evidences the symmetrization of the H-bonds, and indicates that the proton zero-point energy is above the barrier top. The decrease of the heavy-atom distance due to quantum fluctuations in both HF and DF crystals corresponds to a large decrease and an increase of the hydrogen and covalent bond lengths, respectively. Upon deuteration, the increase of the heavy-atom distance (Ubbelohde effect) is in agreement with experimental data.
Collapse
Affiliation(s)
- Hichem Dammak
- Laboratoire Structures Propriétés et Modélisation des Solides, CentraleSupélec, CNRS, Université Paris-Saclay, F 91190 Gif-sur-Yvette, France.
| | | | | | | | | |
Collapse
|
17
|
Schran C, Brieuc F, Marx D. Converged Colored Noise Path Integral Molecular Dynamics Study of the Zundel Cation Down to Ultralow Temperatures at Coupled Cluster Accuracy. J Chem Theory Comput 2018; 14:5068-5078. [PMID: 30217111 DOI: 10.1021/acs.jctc.8b00705] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
For a long time, performing converged path integral simulations at ultralow but finite temperatures of a few Kelvin has been a nearly impossible task. However, recent developments in advanced colored noise thermostatting schemes for path integral simulations, namely, the Path Integral Generalized Langevin Equation Thermostat (PIGLET) and the Path Integral Quantum Thermal Bath (PIQTB), have been able to greatly reduce the computational cost of these simulations, thus making the ultralow temperature regime accessible in practice. In this work, we investigate the influence of these two thermostatting schemes on the description of hydrogen-bonded systems at temperatures down to a few Kelvin as encountered, for example, in helium nanodroplet isolation or tagging photodissociation spectroscopy experiments. For this purpose, we analyze the prototypical hydrogen bond in the Zundel cation (H5O2+) as a function of both oxygen-oxygen distance and temperature in order to elucidate how the anisotropic quantum delocalization and, thus, the shape of the shared proton adapts depending on the donor-acceptor distance. The underlying electronic structure of the Zundel cation is described in terms of Behler's Neural Network Potentials of essentially converged Coupled Cluster accuracy, CCSD(T*)-F12a/AVTZ. In addition, the performances of the PIQTB and PIGLET methods for energetic, structural, and quantum delocalization properties are assessed and directly compared. Overall, our results emphasize the validity and practical usefulness of these two modern thermostatting approaches for path integral simulations of hydrogen-bonded systems even at ultralow temperatures.
Collapse
Affiliation(s)
- 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
|
18
|
|
19
|
Affiliation(s)
- Xinzijian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| | - Jian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| |
Collapse
|
20
|
Mouhat F, Sorella S, Vuilleumier R, Saitta AM, Casula M. Fully Quantum Description of the Zundel Ion: Combining Variational Quantum Monte Carlo with Path Integral Langevin Dynamics. J Chem Theory Comput 2017; 13:2400-2417. [DOI: 10.1021/acs.jctc.7b00017] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Félix Mouhat
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Sandro Sorella
- International
School for Advanced Studies (SISSA), Via Bonomea 26, 34136 Trieste, Italy
- INFM Democritos
National Simulation Center, 34151 Trieste, Italy
| | - Rodolphe Vuilleumier
- PASTEUR,
Département de chimie, École normale supérieure,
UPMC Univ Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ Paris 06, École normale supérieure, CNRS, Processus d’activation
sélective par transfert d’énergie uni-électronique
ou radiatif (PASTEUR), 75005 Paris, France
| | - Antonino Marco Saitta
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Michele Casula
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| |
Collapse
|
21
|
Calvo F, Yurtsever E. Solvation of carbonaceous molecules by para-H2 and ortho-D2 clusters. II. Fullerenes. J Chem Phys 2016; 145:084304. [PMID: 27586919 DOI: 10.1063/1.4961159] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The coating of various fullerenes by para-hydrogen and ortho-deuterium molecules has been computationally studied as a function of the solvent amount. Rotationally averaged interaction potentials for structureless hydrogen molecules are employed to model their interaction with neutral or charged carbonaceous dopants containing between 20 and 240 atoms, occasionally comparing different fullerenes having the same size but different shapes. The solvation energy and the size of the first solvation shell obtained from path-integral molecular dynamics simulations at 2 K show only minor influence on the dopant charge and on the possible deuteration of the solvent, although the shell size is largest for ortho-D2 coating cationic fullerenes. Nontrivial finite size effects have been found with the shell size varying non-monotonically close to its completion limit. For fullerenes embedded in large hydrogen clusters, the shell size and solvation energy both follow linear scaling with the fullerene size. The shell sizes obtained for C60 (+) and C70 (+) are close to 49 and 51, respectively, and agree with mass spectrometry experiments.
Collapse
Affiliation(s)
- F Calvo
- Université Grenoble Alpes, LIPHY, F-38000 Grenoble, France and CNRS, LIPHy, F-38000 Grenoble, France
| | - E Yurtsever
- Koç University, Rumelifeneriyolu, Sariyer, Istanbul 34450, Turkey
| |
Collapse
|
22
|
Liu J, Li D, Liu X. A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat. J Chem Phys 2016; 145:024103. [DOI: 10.1063/1.4954990] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Dezhang Li
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Xinzijian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|