1
|
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
|
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
|
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: 15] [Impact Index Per Article: 7.5] [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
|
4
|
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
|
5
|
Thomas DA, Taccone M, Ober K, Mucha E, Meijer G, von Helden G. Helium Nanodroplet Infrared Action Spectroscopy of the Proton-Bound Dimer of Hydrogen Sulfate and Formate: Examining Nuclear Quantum Effects. J Phys Chem A 2021; 125:9279-9287. [PMID: 34652165 PMCID: PMC8558860 DOI: 10.1021/acs.jpca.1c05705] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The proton-bound dimer of hydrogen sulfate and formate is an archetypal structure for ionic hydrogen-bonding complexes that contribute to biogenic aerosol nucleation. Of central importance for the structure and properties of this complex is the location of the bridging proton connecting the two conjugate base moieties. The potential energy surface for bridging proton translocation features two local minima, with the proton localized at either the formate or hydrogen sulfate moiety. However, electronic structure methods reveal a shallow potential energy surface governing proton translocation, with a barrier on the order of the zero-point energy. This shallow potential complicates structural assignment and necessitates a consideration of nuclear quantum effects. In this work, we probe the structure of this complex and its isotopologues, utilizing infrared (IR) action spectroscopy of ions captured in helium nanodroplets. The IR spectra indicate a structure in which a proton is shared between the hydrogen sulfate and formate moieties, HSO4-···H+···-OOCH. However, because of the nuclear quantum effects and vibrational anharmonicities associated with the shallow potential for proton translocation, the extent of proton displacement from the formate moiety remains unclear, requiring further experiments or more advanced theoretical treatments for additional insight.
Collapse
Affiliation(s)
- Daniel A Thomas
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Martín Taccone
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Katja Ober
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Eike Mucha
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Gerard Meijer
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Gert von Helden
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
6
|
Unke O, Chmiela S, Sauceda HE, Gastegger M, Poltavsky I, Schütt KT, Tkatchenko A, Müller KR. Machine Learning Force Fields. Chem Rev 2021; 121:10142-10186. [PMID: 33705118 PMCID: PMC8391964 DOI: 10.1021/acs.chemrev.0c01111] [Citation(s) in RCA: 347] [Impact Index Per Article: 115.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Indexed: 12/27/2022]
Abstract
In recent years, the use of machine learning (ML) in computational chemistry has enabled numerous advances previously out of reach due to the computational complexity of traditional electronic-structure methods. One of the most promising applications is the construction of ML-based force fields (FFs), with the aim to narrow the gap between the accuracy of ab initio methods and the efficiency of classical FFs. The key idea is to learn the statistical relation between chemical structure and potential energy without relying on a preconceived notion of fixed chemical bonds or knowledge about the relevant interactions. Such universal ML approximations are in principle only limited by the quality and quantity of the reference data used to train them. This review gives an overview of applications of ML-FFs and the chemical insights that can be obtained from them. The core concepts underlying ML-FFs are described in detail, and a step-by-step guide for constructing and testing them from scratch is given. The text concludes with a discussion of the challenges that remain to be overcome by the next generation of ML-FFs.
Collapse
Affiliation(s)
- Oliver
T. Unke
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- DFG
Cluster of Excellence “Unifying Systems in Catalysis”
(UniSysCat), Technische Universität Berlin, 10623 Berlin, Germany
| | - Stefan Chmiela
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Huziel E. Sauceda
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- BASLEARN,
BASF-TU Joint Lab, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Michael Gastegger
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- DFG
Cluster of Excellence “Unifying Systems in Catalysis”
(UniSysCat), Technische Universität Berlin, 10623 Berlin, Germany
- BASLEARN,
BASF-TU Joint Lab, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Igor Poltavsky
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Kristof T. Schütt
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Alexandre Tkatchenko
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- BIFOLD−Berlin
Institute for the Foundations of Learning and Data, Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul 02841, Korea
- Max Planck
Institute for Informatics, Stuhlsatzenhausweg, 66123 Saarbrücken, Germany
- Google
Research, Brain Team, Berlin, Germany
| |
Collapse
|
7
|
Affiliation(s)
- Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstraße 6, 37077 Göttingen, Germany
| |
Collapse
|
8
|
Hu Y, Ounkham P, Marsalek O, Markland TE, Krishmoorthy B, Clark AE. Persistent Homology Metrics Reveal Quantum Fluctuations and Reactive Atoms in Path Integral Dynamics. Front Chem 2021; 9:624937. [PMID: 33748074 PMCID: PMC7973227 DOI: 10.3389/fchem.2021.624937] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2020] [Accepted: 01/22/2021] [Indexed: 11/13/2022] Open
Abstract
Nuclear quantum effects (NQEs) are known to impact a number of features associated with chemical reactivity and physicochemical properties, particularly for light atoms and at low temperatures. In the imaginary time path integral formalism, each atom is mapped onto a "ring polymer" whose spread is related to the quantum mechanical uncertainty in the particle's position, i.e., its thermal wavelength. A number of metrics have previously been used to investigate and characterize this spread and explain effects arising from quantum delocalization, zero-point energy, and tunneling. Many of these shape metrics consider just the instantaneous structure of the ring polymers. However, given the significant interest in methods such as centroid molecular dynamics and ring polymer molecular dynamics that link the molecular dynamics of these ring polymers to real time properties, there exists significant opportunity to exploit metrics that also allow for the study of the fluctuations of the atom delocalization in time. Here we consider the ring polymer delocalization from the perspective of computational topology, specifically persistent homology, which describes the 3-dimensional arrangement of point cloud data, (i.e. atomic positions). We employ the Betti sequence probability distribution to define the ensemble of shapes adopted by the ring polymer. The Wasserstein distances of Betti sequences adjacent in time are used to characterize fluctuations in shape, where the Fourier transform and associated principal components provides added information differentiating atoms with different NQEs based on their dynamic properties. We demonstrate this methodology on two representative systems, a glassy system consisting of two atom types with dramatically different de Broglie thermal wavelengths, and ab initio molecular dynamics simulation of an aqueous 4 M HCl solution where the H-atoms are differentiated based on their participation in proton transfer reactions.
Collapse
Affiliation(s)
- Yunfeng Hu
- Department of Mathematics and Statistics, Washington State University, Pullman, WA, United States
| | - Phonemany Ounkham
- Department of Chemistry, Washington State University, Pullman, WA, United States
| | - Ondrej Marsalek
- Faculty of Mathematics and Physics, Charles University, Prague, Czech
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, CA, United States
| | - Bala Krishmoorthy
- Department of Mathematics and Statistics, Washington State University, Vancouver, WA, United States
| | - Aurora E Clark
- Department of Chemistry, Washington State University, Pullman, WA, United States
| |
Collapse
|
9
|
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
|
10
|
Sauceda HE, Vassilev-Galindo V, Chmiela S, Müller KR, Tkatchenko A. Dynamical strengthening of covalent and non-covalent molecular interactions by nuclear quantum effects at finite temperature. Nat Commun 2021; 12:442. [PMID: 33469007 PMCID: PMC7815839 DOI: 10.1038/s41467-020-20212-1] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Accepted: 11/12/2020] [Indexed: 11/08/2022] Open
Abstract
Nuclear quantum effects (NQE) tend to generate delocalized molecular dynamics due to the inclusion of the zero point energy and its coupling with the anharmonicities in interatomic interactions. Here, we present evidence that NQE often enhance electronic interactions and, in turn, can result in dynamical molecular stabilization at finite temperature. The underlying physical mechanism promoted by NQE depends on the particular interaction under consideration. First, the effective reduction of interatomic distances between functional groups within a molecule can enhance the n → π* interaction by increasing the overlap between molecular orbitals or by strengthening electrostatic interactions between neighboring charge densities. Second, NQE can localize methyl rotors by temporarily changing molecular bond orders and leading to the emergence of localized transient rotor states. Third, for noncovalent van der Waals interactions the strengthening comes from the increase of the polarizability given the expanded average interatomic distances induced by NQE. The implications of these boosted interactions include counterintuitive hydroxyl-hydroxyl bonding, hindered methyl rotor dynamics, and molecular stiffening which generates smoother free-energy surfaces. Our findings yield new insights into the versatile role of nuclear quantum fluctuations in molecules and materials.
Collapse
Affiliation(s)
- Huziel E Sauceda
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg.
- Machine Learning Group, Technische Universität Berlin, 10587, Berlin, Germany.
- BASLEARN, BASF-TU joint Lab, Technische Universität Berlin, 10587, Berlin, Germany.
| | - Valentin Vassilev-Galindo
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg
| | - Stefan Chmiela
- Machine Learning Group, Technische Universität Berlin, 10587, Berlin, Germany
| | - Klaus-Robert Müller
- Machine Learning Group, Technische Universität Berlin, 10587, Berlin, Germany.
- Department of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea.
- Max Planck Institute for Informatics, Stuhlsatzenhausweg, 66123, Saarbrücken, Germany.
- Google Research, Brain team, Berlin, Germany.
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg.
| |
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
|
Schran C, Brezina K, Marsalek O. Committee neural network potentials control generalization errors and enable active learning. J Chem Phys 2020; 153:104105. [DOI: 10.1063/5.0016004] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Christoph Schran
- Charles University, Faculty of Mathematics and Physics, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | - Krystof Brezina
- Charles University, Faculty of Mathematics and Physics, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | - Ondrej Marsalek
- Charles University, Faculty of Mathematics and Physics, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| |
Collapse
|
13
|
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
|
14
|
Schran C, Behler J, Marx D. Automated Fitting of Neural Network Potentials at Coupled Cluster Accuracy: Protonated Water Clusters as Testing Ground. J Chem Theory Comput 2019; 16:88-99. [DOI: 10.1021/acs.jctc.9b00805] [Citation(s) in RCA: 56] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Christoph Schran
- Lehrstuhl für Theoretische Chemie, 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
|
15
|
Bertaina G, Di Liberto G, Ceotto M. Reduced rovibrational coupling Cartesian dynamics for semiclassical calculations: Application to the spectrum of the Zundel cation. J Chem Phys 2019; 151:114307. [PMID: 31542046 DOI: 10.1063/1.5114616] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from the work of Huang et al. [J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows for the physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid initially exciting lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of ∼29 cm-1 with respect to available multiconfiguration time-dependent Hartree calculations and MAE ∼ 14 cm-1 when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster systems.
Collapse
Affiliation(s)
- G Bertaina
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - G Di Liberto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - M Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
16
|
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
|
17
|
Schran C, Marx D. Quantum nature of the hydrogen bond from ambient conditions down to ultra-low temperatures. Phys Chem Chem Phys 2019; 21:24967-24975. [DOI: 10.1039/c9cp04795f] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Quantum simulations reveal strong temperature effects for weak hydrogen bonds and differences in quantum delocalization between various hydrogen-bonded systems.
Collapse
Affiliation(s)
- Christoph Schran
- 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
|