1
|
Montes de Oca-Estévez MJ, Valdés Á, Prosmiti R. A kernel-based machine learning potential and quantum vibrational state analysis of the cationic Ar hydride (Ar 2H +). Phys Chem Chem Phys 2024; 26:7060-7071. [PMID: 38345626 DOI: 10.1039/d3cp05865d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
One of the most fascinating discoveries in recent years, in the cold and low pressure regions of the universe, was the detection of ArH+ and HeH+ species. The identification of such noble gas-containing molecules in space is the key to understanding noble gas chemistry. In the present work, we discuss the possibility of [Ar2H]+ existence as a potentially detectable molecule in the interstellar medium, providing new data on possible astronomical pathways and energetics of this compound. As a first step, a data-driven approach is proposed to construct a full 3D machine-learning potential energy surface (ML-PES) via the reproducing kernel Hilbert space (RKHS) method. The training and testing data sets are generated from CCSD(T)/CBS[56] computations, while a validation protocol is introduced to ensure the quality of the potential. In turn, the resulting ML-PES is employed to compute vibrational levels and molecular spectroscopic constants for the cation. In this way, the most common isotopologue in ISM, [36Ar2H]+, was characterized for the first time, while simultaneously, comparisons with previously reported values available for [40Ar2H]+ are discussed. Our present data could serve as a benchmark for future studies on this system, as well as on higher-order cationic Ar-hydrides of astrophysical interest.
Collapse
Affiliation(s)
- María Judit Montes de Oca-Estévez
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
- Atelgraphics S.L., Mota de Cuervo 42, 28043, Madrid, Spain
| | - Álvaro Valdés
- Escuela de Física, Universidad Nacional de Colombia, Sede Medellín, A. A., 3840, Medellín, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
2
|
Riera M, Knight C, Bull-Vulpe EF, Zhu X, Agnew H, Smith DGA, Simmonett AC, Paesani F. MBX: A many-body energy and force calculator for data-driven many-body simulations. J Chem Phys 2023; 159:054802. [PMID: 37526156 PMCID: PMC10550339 DOI: 10.1063/5.0156036] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 07/11/2023] [Indexed: 08/02/2023] Open
Abstract
Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.
Collapse
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Argonne National Laboratory, Computational Science Division, Lemont, Illinois 60439, USA
| | - Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Xuanyu Zhu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | | | - Andrew C. Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | | |
Collapse
|
3
|
Bull-Vulpe EF, Riera M, Bore SL, Paesani F. Data-Driven Many-Body Potential Energy Functions for Generic Molecules: Linear Alkanes as a Proof-of-Concept Application. J Chem Theory Comput 2022. [PMID: 36113028 DOI: 10.1021/acs.jctc.2c00645] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a generalization of the many-body energy (MB-nrg) theoretical/computational framework that enables the development of data-driven potential energy functions (PEFs) for generic covalently bonded molecules, with arbitrary quantum mechanical accuracy. The "nearsightedness of electronic matter" is exploited to define monomers as "natural building blocks" on the basis of their distinct chemical identity. The energy of generic molecules is then expressed as a sum of individual many-body energies of incrementally larger subsystems. The MB-nrg PEFs represent the low-order n-body energies, with n = 1-4, using permutationally invariant polynomials derived from electronic structure data carried out at an arbitrary quantum mechanical level of theory, while all higher-order n-body terms (n > 4) are represented by a classical many-body polarization term. As a proof-of-concept application of the general MB-nrg framework, we present MB-nrg PEFs for linear alkanes. The MB-nrg PEFs are shown to accurately reproduce reference energies, harmonic frequencies, and potential energy scans of alkanes, independently of their length. Since, by construction, the MB-nrg framework introduced here can be applied to generic covalently bonded molecules, we envision future computer simulations of complex molecular systems using data-driven MB-nrg PEFs, with arbitrary quantum mechanical accuracy.
Collapse
Affiliation(s)
- Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Sigbjørn L. Bore
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
4
|
Prosmiti R, Valdés Á. The smallest proton-bound dimer H 5+: theoretical progress. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2019; 377:20180396. [PMID: 31378176 PMCID: PMC6710890 DOI: 10.1098/rsta.2018.0396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/19/2019] [Indexed: 06/10/2023]
Abstract
The protonated hydrogen dimer, H5+, is the smallest system including proton transfer, and has been of long-standing interest since its first laboratory observation in 1962. H5+ and its isotopologues are the intermediate complexes in deuterium fractionation reactions, and are of central importance in molecular astrophysics. The recently recorded infrared spectra of both H5+ and D5+ reveal a rich vibrational dynamics of the cations, which presents a challenge for standard theoretical approaches. Although H5+ is a four-electron ion, which makes highly accurate electronic structure calculations tractable, the construction of ab initio-based potential energy and dipole moment surfaces has proved a hard task. In the same vein, the difficulties in treating the nuclear motion could also become cumbersome due to their high dimensionality, floppiness and/or symmetry. These systems are prototypical examples for studying large-amplitude motions, as they are highly delocalized, interconverting between equivalent minima through internal rotation and proton transfer motions requiring state-of-the-art treatments. Recent advances in the computational vibrational spectroscopy of the H5+ cation and its isotopologues are reported from full quantum spectral simulations, providing important information in a rigorous manner, and open perspectives for further future investigations. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.
Collapse
Affiliation(s)
- Rita Prosmiti
- Departamento PAMS, Instituto de Física Fundamental (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain
| | - Álvaro Valdés
- Departamento de Física, Universidad Nacional de Colombia, Calle 26, Cra 39, Edificio 404, Bogotá, Colombia
| |
Collapse
|
5
|
Affiliation(s)
- Chen Qu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Qi Yu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M. Bowman
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
6
|
Ziegler B, Rauhut G. Efficient generation of sum-of-products representations of high-dimensional potential energy surfaces based on multimode expansions. J Chem Phys 2016; 144:114114. [DOI: 10.1063/1.4943985] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Benjamin Ziegler
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
7
|
Han H, Song H, Li J, Guo H. Near Spectroscopically Accurate Ab Initio Potential Energy Surface for NH4+ and Variational Calculations of Low-Lying Vibrational Levels. J Phys Chem A 2015; 119:3400-6. [DOI: 10.1021/acs.jpca.5b01835] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Huixian Han
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
- School
of Physics, Northwest University, Xi’an, Shaanxi 710069, China
| | - Hongwei Song
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Jun Li
- School
of Chemistry and Chemical Engineering, Chongqing University, Chongqing 400030, China
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
8
|
Qu C, Bowman JM. Diffusion Monte Carlo Calculations of Zero-Point Structures of Partially Deuterated Isotopologues of H7+. J Phys Chem B 2014; 118:8221-6. [DOI: 10.1021/jp501371z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Chen Qu
- Department of Chemistry and
Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Department of Chemistry and
Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
9
|
Glushkov V, Wilson S. On the Coulson–Fischer wave function for the X1A′ H +3molecular ion: parametrisation using distributed Gaussian basis sets. Mol Phys 2014. [DOI: 10.1080/00268976.2013.812256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
10
|
Barragán P, Pérez de Tudela R, Qu C, Prosmiti R, Bowman JM. Full-dimensional quantum calculations of the dissociation energy, zero-point, and 10 K properties of H7+/D7+ clusters using an ab initio potential energy surface. J Chem Phys 2014; 139:024308. [PMID: 23862944 DOI: 10.1063/1.4812557] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Diffusion Monte Carlo (DMC) and path-integral Monte Carlo computations of the vibrational ground state and 10 K equilibrium state properties of the H7 (+)/D7 (+) cations are presented, using an ab initio full-dimensional potential energy surface. The DMC zero-point energies of dissociated fragments H5 (+)(D5 (+))+H2(D2) are also calculated and from these results and the electronic dissociation energy, dissociation energies, D0, of 752 ± 15 and 980 ± 14 cm(-1) are reported for H7 (+) and D7 (+), respectively. Due to the known error in the electronic dissociation energy of the potential surface, these quantities are underestimated by roughly 65 cm(-1). These values are rigorously determined for first time, and compared with previous theoretical estimates from electronic structure calculations using standard harmonic analysis, and available experimental measurements. Probability density distributions are also computed for the ground vibrational and 10 K state of H7 (+) and D7 (+). These are qualitatively described as a central H3 (+)/D3 (+) core surrounded by "solvent" H2/D2 molecules that nearly freely rotate.
Collapse
Affiliation(s)
- Patricia Barragán
- Instituto de Física Fundamental, IFF-CSIC, Serrano 123, 28006 Madrid, Spain
| | | | | | | | | |
Collapse
|
11
|
Qu C, Prosmiti R, Bowman JM. MULTIMODE calculations of the infrared spectra of H 7 + and D 7 + using ab initio potential energy and dipole moment surfaces. Theor Chem Acc 2013. [DOI: 10.1007/s00214-013-1413-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
12
|
Valdés Á, Prosmiti R, Delgado-Barrio G. Vibrational dynamics of the H5(+) and its isotopologues from multiconfiguration time-dependent Hartree calculations. J Chem Phys 2013; 137:214308. [PMID: 23231232 DOI: 10.1063/1.4769081] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Full-dimensional multiconfiguration time-dependent Hartree (MCTDH) computations are reported for the vibrational states of the H(5)(+) and its H(4)D(+), H(3)D(2)(+), H(2)D(3)(+), HD(4)(+), D(5)(+) isotopologues employing two recent analytical potential energy surfaces of Xie et al. [J. Chem. Phys. 122, 224307 (2005)] and Aguado et al. [J. Chem. Phys. 133, 024306 (2010)]. The potential energy operators are constructed using the n-mode representation adapted to a four-combined mode cluster expansion, including up to seven-dimensional grids, chosen adequately to take advantage in representing the MCTDH wavefunction. An error analysis is performed to quantify the convergence of the potential expansion to reproduce the reference surfaces at the energies of interest. An extensive analysis of the vibrational ground state properties of these isotopes and comparison with the reference diffusion Monte Carlo results by Acioli et al. [J. Chem. Phys. 128, 104318 (2008)] are presented. It is found that these systems are highly delocalized, interconverting between equivalent minima through rotation and internal proton transfer motions even at their vibrational ground state. Isotopic substitution affects the zero-point energy and structure, showing preference in the arrangements of the H and D within the mixed clusters, and the most stable conformers of each isotopomer are the ones with the H in the central position. Vibrational excited states are also computed and by comparing the energies and structures predicted from the two surfaces, the effect of the potential topology on them is discussed.
Collapse
Affiliation(s)
- Álvaro Valdés
- Instituto de Física Fundamental, IFF-CSIC, Serrano 123, 28006 Madrid, Spain
| | | | | |
Collapse
|
13
|
Young JW, Cheng TC, Bandyopadhyay B, Duncan MA. IR photodissociation spectroscopy of H7(+), H9(+), and their deuterated analogues. J Phys Chem A 2013; 117:6984-90. [PMID: 23374094 DOI: 10.1021/jp312630x] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cluster ions of H7(+)/D7(+) and H9(+)/D9(+) produced in a supersonic molecular beam with a pulsed discharge source are mass selected and studied with infrared laser photodissociation spectroscopy. Photodissociation occurs by the loss of H2 (D2) from each cluster, producing resonances in the 2000-4500 cm(-1) region. Vibrational patterns indicate that these ions consist of an H3(+) (D3(+)) core ion solvated by H2 (D2) molecules. There is no evidence for the shared proton structure seen previously for H5(+). The H3(+) ion core vibrational bands are weakened and broadened significantly, presumably by enhanced rates of intramolecular vibrational relaxation. Computational studies at the DFT/B3LYP or MP2 levels of theory (including scaling) are adequate to reproduce qualitative details of the vibrational spectra, but neither provides quantitative agreement with vibrational frequencies.
Collapse
Affiliation(s)
- J W Young
- Department of Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | | | | | | |
Collapse
|