1
|
Zhang X, Huang K, Fu Y, Zhang N, Kong X, Cheng Y, Zheng M, Cheng Y, Zhu T, Fu B, Feng S, Chen H. Demethylation C-C coupling reaction facilitated by the repulsive Coulomb force between two cations. Nat Commun 2024; 15:5881. [PMID: 38997250 DOI: 10.1038/s41467-024-49946-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2023] [Accepted: 06/26/2024] [Indexed: 07/14/2024] Open
Abstract
Carbon chain elongation (CCE) is normally carried out using either chemical catalysts or bioenzymes. Herein we demonstrate a catalyst-free approach to promote demethylation C-C coupling reactions for advanced CCE constructed with functional groups under ambient conditions. Accelerated by the electric field, two organic cations containing a methyl group (e.g., ketones, acids, and aldehydes) approach each other with such proximity that the energy of the repulsive Coulomb interaction between these two cations exceeds the bond energy of the methyl group. This results in the elimination of a methyl cation and the coupling of the residual carbonyl carbon groups. As confirmed by high-resolution mass spectrometry and isotope-labeling experiments, the C-C coupling reactions (yields up to 76.5%) were commonly observed in the gas phase or liquid phase, for which the mechanism was further studied using molecular dynamics simulations and stationary-point calculations, revealing deep insights and perspectives of chemistry.
Collapse
Affiliation(s)
- Xiaoping Zhang
- Jiangxi Key Laboratory for Mass Spectrometry and Instrumentation, East China University of Technology, Nanchang, 330013, P. R. China
| | - Keke Huang
- State Key Laboratory of Inorganic Synthesis and Preparative Chemistry, Jilin University, Changchun, 130012, P. R. China
| | - Yanlin Fu
- State Key Laboratory of Molecular Reaction Dynamics, Center for Theoretical and Computational Chemistry, Dalian Institution of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, P. R. China
| | - Ni Zhang
- School of Pharmacy, Jiangxi University of Chinese Medicine, Nanchang, 330004, P. R. China
| | - Xianglei Kong
- State Key Laboratory of Elemento-organic Chemistry, College of Chemistry, Nankai University, Tianjin, 300071, P. R. China
| | - Yuanyuan Cheng
- Jiangxi Key Laboratory for Mass Spectrometry and Instrumentation, East China University of Technology, Nanchang, 330013, P. R. China
| | - Mingyu Zheng
- Jiangxi Key Laboratory for Mass Spectrometry and Instrumentation, East China University of Technology, Nanchang, 330013, P. R. China
| | - Yihao Cheng
- Jiangxi Key Laboratory for Mass Spectrometry and Instrumentation, East China University of Technology, Nanchang, 330013, P. R. China
| | - Tenggao Zhu
- School of Pharmacy, Jiangxi University of Chinese Medicine, Nanchang, 330004, P. R. China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics, Center for Theoretical and Computational Chemistry, Dalian Institution of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, P. R. China.
- Hefei National Laboratory, Hefei, 230088, China.
| | - Shouhua Feng
- State Key Laboratory of Inorganic Synthesis and Preparative Chemistry, Jilin University, Changchun, 130012, P. R. China.
| | - Huanwen Chen
- Jiangxi Key Laboratory for Mass Spectrometry and Instrumentation, East China University of Technology, Nanchang, 330013, P. R. China.
- School of Pharmacy, Jiangxi University of Chinese Medicine, Nanchang, 330004, P. R. China.
| |
Collapse
|
2
|
Sharma M, Banik S, Roy TK. Performance of Effective Harmonic Oscillator Approach for the Calculations of Vibrational Transition Energies of Large Molecules. J Phys Chem A 2024. [PMID: 38979981 DOI: 10.1021/acs.jpca.4c01583] [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
The accuracy and performance of the effective harmonic oscillator approximation for the description of anharmonic vibrational structure calculations are tested for large molecular systems and compared with experimental values along with vibrational self-consistent field and second-order perturbation theories. The effective harmonic oscillator approach is an effective single-particle approximation where the variational parameters are the centroids and widths of the multidimensional Gaussian product functions posited as the vibrational wave functions. A comprehensive calculation for 849 transitions that include the fundamentals, two and three quanta overtone transitions, and several combination bands of three polyaromatic hydrocarbons and one DNA nucleobase with a total of 231 normal modes are assessed. A comparison of EHO results with the experimental values is done for the polyaromatic hydrocarbons, and a close agreement is found between the two results. It also offers anharmonic eigenstates and eigenfunctions that are nearly identical with vibrational self-consistent field theory. An extensive analysis on the resultant wave functions of the excited states is performed. The overall root-mean-square deviation (RMSD) between these two methods for 849 transitions understudy is only about 8.3 cm-1, suggesting the effective harmonic oscillator as a viable alternative for the reliable calculations of transition energies of large molecular systems.
Collapse
Affiliation(s)
- Mokshi Sharma
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, Jammu and Kashmir 181143, India
| | - Subrata Banik
- Department of Chemistry, School of Chemical and Biotechnology, SASTRA Deemed University, Thanjavur, Tamil Nadu 613401, India
| | - Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, Jammu and Kashmir 181143, India
| |
Collapse
|
3
|
Xiong F, Hou S, Li J, Wang Z, Xie C. Computational determination of the S1(Ã1A″) absorption spectra of HONO and DONO using full-dimensional neural network potential energy surfaces. J Chem Phys 2024; 161:014305. [PMID: 38953448 DOI: 10.1063/5.0216840] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2024] [Accepted: 06/14/2024] [Indexed: 07/04/2024] Open
Abstract
The Ã1A″ ← X̃1A' absorption spectra of HONO and DONO were simulated by a full six-dimensional quantum mechanical method based on the newly constructed potential energy surfaces for the ground and excited electronic states, which were represented by the neural network method utilizing over 36 000 ab initio energy points calculated at the multireference configuration interaction level with Davidson correction. The absorption spectrum of HONO/DONO comprises a superposition of the spectra from two isomers, namely, trans- and cis-HONO/DONO, due to their coexistence in the ground X̃1A' state. Our calculated spectra of both HONO and DONO were found to be in fairly good agreement with the experiment, including the energy positions and widths of the peaks. The dominant progression was assigned to the N=O stretch mode (20n) associated with trans-HONO/DONO, which can be attributed to the promotion of an electron to the π* orbital of N=O. Specifically, the resonances with higher vibrational quanta were found to be in the domain of the Feshbach-type resonances. The assignments of the spectra and mode specificity therein are discussed.
Collapse
Affiliation(s)
- Feng Xiong
- Institute of Modern Physics, Northwest University, Xi'an 710127, China
| | - Siting Hou
- Institute of Modern Physics, Northwest University, Xi'an 710127, China
| | - Jiayuan Li
- Institute of Modern Physics, Northwest University, Xi'an 710127, China
| | - Zhimo Wang
- Institute of Modern Physics, Northwest University, Xi'an 710127, China
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an 710127, China
- Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an 710127, China
| |
Collapse
|
4
|
Yang J, Li J, Li J, Li J. Gaussian Process Regression for State-to-State Integral Cross Sections: The Case of the O + O 2 Collision Dissociation Reactions. J Phys Chem A 2024; 128:4966-4975. [PMID: 38869143 DOI: 10.1021/acs.jpca.4c01445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2024]
Abstract
Research on hypersonic vehicles has become increasingly important worldwide in recent years. However, accurately simulating the dynamics of the nonequilibrium high-temperature reactions that are in the hypersonic flow around the vehicles presents a significant challenge as a large number of states and transitions are accessible even for the smallest atom-diatom reaction systems. It is quite difficult, sometimes even impossible, to exhaustively investigate all relevant combinations or determine high-dimensional analytical representations for the state-to-state reaction probabilities. In this study, we used Gaussian process regression (GPR) to fit a model based on only 807 QCT data for training. The confidence interval of the GPR prediction and the Kullback-Leibler (KL) divergence were used to help minimize the sampling amount of data for fitting the converged GPR model. The model aims to predict the state-to-state integral cross section (ICS) of the O + O2 → 3O dissociation reaction under random initial conditions (Et, v, j). In total, it took almost a month to obtain this converged GPR model, but it took only a few seconds to predict the ICS value for any initial condition. For 330 initial conditions not included in the training set, the mean-square error (MSE) between the QCT-calculated ICSs and the GPR-predicted ones is only 0.08 Å2 and the R2 is 0.9986, indicating that the GPR model can replace the direct expensive QCT calculation with high accuracy. Finally, we calculated the equilibrium dissociation rate coefficients based on the StS ICS values predicted by the GPR model, and the results were in good agreement with available experimental and theoretical results. Thus, this study provides an effective and accurate approach to the extensive direct state-to-state reaction dynamic calculations.
Collapse
Affiliation(s)
- Jiawei Yang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Jia Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Junhong Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| |
Collapse
|
5
|
Jiang J, Yang J, Hong Q, Sun Q, Li J. Global Potential Energy Surfaces by Compressed-State Multistate Pair-Density Functional Theory for Hyperthermal Collisions in the O 2+O 2 System. Chemphyschem 2024; 25:e202400078. [PMID: 38526528 DOI: 10.1002/cphc.202400078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Revised: 03/02/2024] [Accepted: 03/25/2024] [Indexed: 03/26/2024]
Abstract
Interactions between oxygen molecules play an important role in atmospheric chemistry and hypersonic flow chemistry in atmospheric entries. Recently, high-quality ab initio potential energy surface (PES) of the quintet O4 was reported by Paukku et al. [J. Chem. Phys. 147, 034301 (2017)]. 10543 configurations were sampled and calculated at the level of MS-CASPT2/maug-cc-pVTZ with scaled external correlation. The PES was fitted to a many-body (MB) form with the many-body part described by the permutationally invariant polynomial approach (MB-PIP). In this work, the PIP-Neural Network (PIP-NN) and MB-PIP-NN methods were used to refit the PES based on the same data by Paukku et al. Three PESs were compared. It was found that the performances differ significantly in the O+O3 region as well as in the long-range region. Therefore, additional 1300 points were sampled, and the efficient compressed-state multistate pair-density functional theory (CMS-PDFT) was used to calculate the electronic structure of these 1300 points and 10543 points by Paukku et al. Then, a completely new quintet PES was fitted using the MB-PIP-NN method. Based on this PES, the quasi-classical trajectory (QCT) approach was used to reveal all possible reaction channels for hyperthermal O2-O2 collisions.
Collapse
Affiliation(s)
- Jie Jiang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing, 401331, China
| | - Jiawei Yang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing, 401331, China
| | - Qizhen Hong
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190, Beijing, China
| | - Quanhua Sun
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190, Beijing, China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing, 401331, China
| |
Collapse
|
6
|
Shu Y, Truhlar DG. Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients. J Chem Theory Comput 2024; 20:4396-4426. [PMID: 38819014 DOI: 10.1021/acs.jctc.4c00424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
7
|
Czakó G, Gruber B, Papp D, Tajti V, Tasi DA, Yin C. First-principles mode-specific reaction dynamics. Phys Chem Chem Phys 2024; 26:15818-15830. [PMID: 38639072 DOI: 10.1039/d4cp00417e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/20/2024]
Abstract
Controlling the outcome of chemical reactions by exciting specific vibrational and/or rotational modes of the reactants is one of the major goals of modern reaction dynamics studies. In the present Perspective, we focus on first-principles vibrational and rotational mode-specific dynamics computations on reactions of neutral and anionic systems beyond six atoms such as X + C2H6 [X = F, Cl, OH], HX + C2H5 [X = Br, I], OH- + CH3I, and F- + CH3CH2Cl. The dynamics simulations utilize high-level ab initio analytical potential energy surfaces and the quasi-classical trajectory method. Besides initial state specificity and the validity of the Polanyi rules, mode-specific vibrational-state assignment for polyatomic product species using normal-mode analysis and Gaussian binning is also discussed and compared with experiment.
Collapse
Affiliation(s)
- Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Balázs Gruber
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Viktor Tajti
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Domonkos A Tasi
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Cangtao Yin
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
8
|
Spencer RJ, Zhanserkeev AA, Yang EL, Steele RP. The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings. J Am Chem Soc 2024; 146:15376-15392. [PMID: 38771156 DOI: 10.1021/jacs.4c03198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.
Collapse
Affiliation(s)
- Ryan J Spencer
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Asylbek A Zhanserkeev
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Emily L Yang
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| |
Collapse
|
9
|
Rodriguez L, Natalizio M, Sode O. Theoretical Insights into the Vibrational Structure of Carbon Dioxide Rare-Gas Complexes. J Phys Chem A 2024; 128:4199-4205. [PMID: 38770817 DOI: 10.1021/acs.jpca.4c00639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar-CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.
Collapse
Affiliation(s)
- Larry Rodriguez
- Department of Chemistry and Biochemistry, California State University, Los Angeles, 5151 State University Drive, Los Angeles, California 90032, United States
| | - Michael Natalizio
- Department of Chemistry and Biochemistry, California State University, Los Angeles, 5151 State University Drive, Los Angeles, California 90032, United States
| | - Olaseni Sode
- Department of Chemistry and Biochemistry, California State University, Los Angeles, 5151 State University Drive, Los Angeles, California 90032, United States
| |
Collapse
|
10
|
Weike N, Fritsch F, Eisfeld W. Compensation States Approach in the Hybrid Diabatization Scheme: Extension to Multidimensional Data and Properties. J Phys Chem A 2024; 128:4353-4368. [PMID: 38748493 DOI: 10.1021/acs.jpca.4c01134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2024]
Abstract
The diabatization of reactive systems for more than just a couple of states is a very demanding problem and generally requires advanced diabatization techniques. Especially for dissociative processes, the drastic changes in the adiabatic wave functions often would require large diabatic state bases, which quickly become impractical. Recently, we addressed this problem by the compensation states approach developed in the context of our hybrid diabatization scheme. This scheme utilizes wave function as well as energy data in combination with a diabatic potential model. In regions where the initial diabatic state basis becomes insufficient for an appropriate representation of the adiabatic states, new model states are generated. The new model states compensate for the state space not spanned by the initial diabatic basis. Such a compensation state is obtained by projecting the initial diabatic state space out of the adiabatic wave function. This yields a very efficient basis representation of the electronic Hamiltonian. The present work presents two new aspects. First, it is shown how other operators like the spin-orbit operator in the framework of the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) can be evaluated in this compact model state space without losing the correct wave function information and accuracy. Second, the extension of the approach to multidimensional potential energy surface models is presented for methyl iodide including the C-I dissociation coordinate and the angular H3C-I bending coordinates.
Collapse
Affiliation(s)
- Nicole Weike
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Fabian Fritsch
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
11
|
Liu Y, Ončák M, Meyer J, Ard SG, Shuman NS, Viggiano AA, Guo H. Multistate Dynamics and Kinetics of CO 2 Activation by Ta + in the Gas Phase: Insights into Single-Atom Catalysis. J Am Chem Soc 2024; 146:14182-14193. [PMID: 38741473 DOI: 10.1021/jacs.4c03192] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
The activation of carbon dioxide (CO2) by a transition-metal cation in the gas phase is a unique model system for understanding single-atom catalysis. The mechanism of such reactions is often attributed to a "two-state reactivity" model in which the high-energy barrier of a spin state correlating with ground-state reactants is avoided by intersystem crossing (ISC) to a different spin state with a lower barrier. However, such a "spin-forbidden" mechanism, along with the corresponding dynamics, has seldom been rigorously examined theoretically, due to the lack of global potential energy surfaces (PESs). In this work, we report full-dimensional PESs of the lowest-lying quintet, triplet, and singlet states of the TaCO2+ system, machine-learned from first-principles data. These PESs and the corresponding spin-orbit couplings enable us to provide an extensive theoretical characterization of the dynamics and kinetics of the reaction between the tantalum cation (Ta+) and CO2, which have recently been investigated experimentally at high collision energies using crossed beams and velocity map imaging, as well as at thermal energies using a selected-ion flow tube apparatus. The multistate quasi-classical trajectory simulations with surface hopping reproduce most of the measured product translational and angular distributions, shedding valuable light on the nonadiabatic reaction dynamics. The calculated rate coefficients from 200 to 600 K are also in good agreement with the latest experimental measurements. More importantly, these calculations revealed that the reaction is controlled by intersystem crossing, rather than potential barriers.
Collapse
Affiliation(s)
- Yang Liu
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Milan Ončák
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstra. 25, 6020 Innsbruck, Austria
| | - Jennifer Meyer
- Fachbereich Chemie und Forschungszentrum OPTIMAS, RPTU Kaiserslautern-Landau, Erwin-Schrödinger Str. 52, 67663 Kaiserslautern, Germany
| | - Shaun G Ard
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Nicholas S Shuman
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Albert A Viggiano
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
12
|
Gibbas B, Kaledin M, Kaledin AL. Quantum Monte Carlo Simulations of the Vibrational Wavefunction of the Aromatic Cyclo[10]carbon Using a Full Dimensional Permutationally Invariant Potential Energy Surface. J Phys Chem Lett 2024; 15:5070-5075. [PMID: 38701515 PMCID: PMC11103689 DOI: 10.1021/acs.jpclett.4c00893] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2024] [Revised: 04/24/2024] [Accepted: 05/01/2024] [Indexed: 05/05/2024]
Abstract
New experimental measurements [Sun et al., Nature 2023, 623, 972] of the cyclic C10 reveal a cumulenic pentagon-like D5h structure at ∼5 K. However, the long-standing presumption that a large zero-point vibrational energy combined with an extremely flat D5h ↔ D10h ↔ D5h isomerization pathway washes out the pentagonal D5h structure and yields a symmetric D10h decagon remains at odds with the experiment. We resolve this issue with our fitting approach based on a bond-order charge-density matrix expressed in permutationally invariant polynomials. We train the model on τHCTH/cc-pVQZ data morphed to reproduce a relativistic all-electron CCSDT(Q)/CBS D5h-D10h potential energy barrier (benchmarked previously by others). Large scale diffusion Monte Carlo simulations in full dimensionality show that the vibrational ground state of C10 has compositional character of more than 96% D5h, fully reflecting the experimental imaging data. Quantum mechanical variational calculations in 1-D further suggest persistence of the D5h symmetry structure at higher temperatures.
Collapse
Affiliation(s)
- Benjamin
D. Gibbas
- Department
of Chemistry & Biochemistry, Kennesaw
State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department
of Chemistry & Biochemistry, Kennesaw
State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Alexey L. Kaledin
- Cherry
L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| |
Collapse
|
13
|
Shi Z, Lele AD, Jasper AW, Klippenstein SJ, Ju Y. Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data. J Phys Chem A 2024; 128:3449-3457. [PMID: 38642065 DOI: 10.1021/acs.jpca.4c00750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/22/2024]
Abstract
Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. The accuracy of the force field is demonstrated by calculating rate constants for the H + HO2 → H2 + 3O2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.
Collapse
Affiliation(s)
- Zhiyu Shi
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Aditya Dilip Lele
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ahren W Jasper
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Stephen J Klippenstein
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Yiguang Ju
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
14
|
Pandey P, Arandhara M, Houston PL, Qu C, Conte R, Bowman JM, Ramesh SG. Assessing Permutationally Invariant Polynomial and Symmetric Gradient Domain Machine Learning Potential Energy Surfaces for H 3O 2. J Phys Chem A 2024; 128:3212-3219. [PMID: 38624168 PMCID: PMC11056970 DOI: 10.1021/acs.jpca.4c01044] [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/16/2024] [Revised: 03/15/2024] [Accepted: 03/20/2024] [Indexed: 04/17/2024]
Abstract
The singly hydrated hydroxide anion OH-(H2O) is of central importance to a detailed molecular understanding of water; therefore, there is strong motivation to develop a highly accurate potential to describe this anion. While this is a small molecule, it is necessary to have an extensive data set of energies and, if possible, forces to span several important stationary points. Here, we assess two machine-learned potentials, one using the symmetric gradient domain machine learning (sGDML) method and one based on permutationally invariant polynomials (PIPs). These are successors to a PIP potential energy surface (PES) reported in 2004. We describe the details of both fitting methods and then compare the two PESs with respect to precision, properties, and speed of evaluation. While the precision of the potentials is similar, the PIP PES is much faster to evaluate for energies and energies plus gradient than the sGDML one. Diffusion Monte Carlo calculations of the ground vibrational state, using both potentials, produce similar large anharmonic downshift of the zero-point energy compared to the harmonic approximation of the PIP and sGDML potentials. The computational time for these calculations using the sGDML PES is roughly 300 times greater than using the PIP one.
Collapse
Affiliation(s)
- Priyanka Pandey
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Mrinal Arandhara
- Department
of Inorganic and Physical Chemistry, Indian
Institute of Science, Bangalore 560012, India
| | - Paul L. Houston
- Department
of Chemistry and Chemical Biology, Cornell
University, Ithaca, New York 14853, United States
- Department
of Chemistry and Biochemistry, Georgia Institute
of Technology, Atlanta, Georgia 30332, United States
| | - Chen Qu
- Independent
Researcher, Toronto, Ontario M9B0E3, Canada
| | - Riccardo Conte
- Dipartimento
di Chimica, Università degli Studi
di Milano, Milano 20133, Italy
| | - Joel M. Bowman
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Sai G. Ramesh
- Department
of Inorganic and Physical Chemistry, Indian
Institute of Science, Bangalore 560012, India
| |
Collapse
|
15
|
Ge F, Wang R, Qu C, Zheng P, Nandi A, Conte R, Houston PL, Bowman JM, Dral PO. Tell Machine Learning Potentials What They Are Needed For: Simulation-Oriented Training Exemplified for Glycine. J Phys Chem Lett 2024; 15:4451-4460. [PMID: 38626460 DOI: 10.1021/acs.jpclett.4c00746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2024]
Abstract
Machine learning potentials (MLPs) are widely applied as an efficient alternative way to represent potential energy surfaces (PESs) in many chemical simulations. The MLPs are often evaluated with the root-mean-square errors on the test set drawn from the same distribution as the training data. Here, we systematically investigate the relationship between such test errors and the simulation accuracy with MLPs on an example of a full-dimensional, global PES for the glycine amino acid. Our results show that the errors in the test set do not unambiguously reflect the MLP performance in different simulation tasks, such as relative conformer energies, barriers, vibrational levels, and zero-point vibrational energies. We also offer an easily accessible solution for improving the MLP quality in a simulation-oriented manner, yielding the most precise relative conformer energies and barriers. This solution also passed the stringent test by diffusion Monte Carlo simulations.
Collapse
Affiliation(s)
- Fuchun Ge
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Ran Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Chen Qu
- Independent Researcher, Toronto, Ontario M9B0E3, Canada
| | - Peikun Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
- Department of Physics and Materials Science, University of Luxembourg, Luxembourg City L-1511, Luxembourg
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
16
|
Houston PL, Qu C, Yu Q, Pandey P, Conte R, Nandi A, Bowman JM. No Headache for PIPs: A PIP Potential for Aspirin Runs Much Faster and with Similar Precision Than Other Machine-Learned Potentials. J Chem Theory Comput 2024; 20:3008-3018. [PMID: 38593438 DOI: 10.1021/acs.jctc.4c00054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Assessments of machine-learning (ML) potentials are an important aspect of the rapid development of this field. We recently reported an assessment of the linear-regression permutationally invariant polynomial (PIP) method for ethanol, using the widely used (revised) rMD17 data set. We demonstrated that the PIP approach outperformed numerous other methods, e.g., ANI, PhysNet, sGDML, and p-KRR, with respect to precision and notably with respect to speed [Houston et al., J. Chem. Phys. 2022, 156, 044120]. Here, we extend this assessment to the 21-atom aspirin molecule, using the rMD17 data set, with a focus on the speed of evaluation. Both energies and forces are used for training, and the precision of several PIPs is examined for both. Normal mode frequencies, the methyl torsional potential, and 1d vibrational energies for an OH stretch are presented. We show that the PIP approach achieves the level of precision obtained from other ML methods, e.g., atom-centered neural network methods, linear regression ACE, and kernel methods, as reported by Kovács et al. in J. Chem. Theory Comput. 2021, 17, 7696-7711. More significantly, we show that the PIP PESs run much faster than all other ML methods, whose timings were evaluated in that paper. We also show that the PIP PES extrapolates well enough to describe several internal motions of aspirin, including an OH stretch.
Collapse
Affiliation(s)
- Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Chen Qu
- Independent Researcher, Toronto, Ontario M9B0E3, Canada
| | - Qi Yu
- Department of Chemistry, Fudan University, Shanghai 200438, P. R. China
| | - Priyanka Pandey
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
- Department of Physics and Materials Science, University of Luxembourg, Luxembourg City L-1511, Luxembourg
| | - Joel M Bowman
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
17
|
Guan Y, Chen Q, Varandas AJC. Accurate diabatization based on combined-hyperbolic-inverse-power-representation: 1,2 2A' states of BeH2. J Chem Phys 2024; 160:154105. [PMID: 38624109 DOI: 10.1063/5.0200732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2024] [Accepted: 03/28/2024] [Indexed: 04/17/2024] Open
Abstract
A diabatic potential energy matrix (DPEM) for the two lowest states of BeH2+ has been constructed using the combined-hyperbolic-inverse-power-representation (CHIPR) method. By imposing symmetry constraints on the coefficients of polynomials, the complete nuclear permutation inversion symmetry is correctly preserved in the CHIPR functional form. The symmetrized CHIPR functional form is then used in the diabatization by ansatz procedure. The ab initio energies are reproduced with satisfactory accuracy. In addition, the CHIPR-based DPEM also reproduces the local topology of a conical intersection. Future work will focus on a complete four-state diabatic representation with emphasis on the long-range interactions and spin-orbit couplings, which will enable accurate quantum scattering calculations for the Be+(2P) + H2 → BeH+(X1Σ+) + H(2S) reaction.
Collapse
Affiliation(s)
- Yafu Guan
- 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
| | - Qun Chen
- 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
| | - António J C Varandas
- School of Physics and Physical Engineering, Qufu Normal University, 273165 Qufu, People's Republic of China
- Department of Physics, Universidade Federal do Espírito Santo, 29075-910 Vitória, Brazil
- Department of Chemistry and Chemistry Centre, University of Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
18
|
Allen AEA, Csányi G. Toward transferable empirical valence bonds: Making classical force fields reactive. J Chem Phys 2024; 160:124108. [PMID: 38526105 DOI: 10.1063/5.0196952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 03/07/2024] [Indexed: 03/26/2024] Open
Abstract
The empirical valence bond technique allows classical force fields to model reactive processes. However, parametrization from experimental data or quantum mechanical calculations is required for each reaction present in the simulation. We show that the parameters present in the empirical valence bond method can be predicted using a neural network model and the SMILES strings describing a reaction. This removes the need for quantum calculations in the parametrization of the empirical valence bond technique. In doing so, we have taken the first steps toward defining a new procedure for enabling reactive atomistic simulations. This procedure would allow researchers to use existing classical force fields for reactive simulations, without performing additional quantum mechanical calculations.
Collapse
Affiliation(s)
- Alice E A Allen
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Gábor Csányi
- Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| |
Collapse
|
19
|
Dékány AÁ, Czakó G. Detailed quasiclassical dynamics of the F - + SiH 3Cl multi-channel reaction. Phys Chem Chem Phys 2024; 26:10008-10020. [PMID: 38482549 DOI: 10.1039/d4cp00048j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/28/2024]
Abstract
We report a detailed quasiclassical trajectory study on the F- + SiH3Cl multi-channel reaction using a full-dimensional ab initio analytical potential energy surface. Reaction probabilities, cross sections, initial attack and scattering angle distributions as well as product relative translational, internal, vibrational, and rotational energy distributions are obtained in the collision energy range of 1-40 kcal mol-1 for the following channels: SiH3F + Cl-, SiH2Cl- + HF, SiH2F- + HCl, SiH2FCl + H-, SiH2 + FHCl-, and SiHFCl- + H2. All the channels are translationally cold indicating indirect mechanisms, except proton transfer (SiH2Cl- + HF), which shows mixed direct-indirect character. The angular distributions vary depending on collision energy and inversion/retention for SiH3F + Cl-. In the case of SiH2Cl- + HF front-side/back-side attack backward-forward/forward scattering preference is found at low/high collision energy. SiH2F- + HCl is formed with isotropic scattering and the preferred angle of attack is similar to the SiH3F + Cl- channel. SiH2FCl + H-/SiH2 + FHCl- favors back-side attack and isotropic/backward scattering, whereas SiHFCl- + H2 does not show significant angular preference.
Collapse
Affiliation(s)
- Attila Á Dékány
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
20
|
Ellerbrock R, Johnson KG, Seritan S, Hoppe H, Zhang JH, Lenzen T, Weike T, Manthe U, Martínez TJ. QuTree: A tree tensor network package. J Chem Phys 2024; 160:112501. [PMID: 38497471 DOI: 10.1063/5.0180233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2023] [Accepted: 02/20/2024] [Indexed: 03/19/2024] Open
Abstract
We present QuTree, a C++ library for tree tensor network approaches. QuTree provides class structures for tensors, tensor trees, and related linear algebra functions that facilitate the fast development of tree tensor network approaches such as the multilayer multiconfigurational time-dependent Hartree approach or the density matrix renormalization group approach and its various extensions. We investigate the efficiency of relevant tensor and tensor network operations and show that the overhead for managing the network structure is negligible, even in cases with a million leaves and small tensors. QuTree focuses on providing simple, high-level routines while retaining easy access to the backend to facilitate novel developments. We demonstrate the capabilities of the package by computing the eigenstates of coupled harmonic oscillator Hamiltonians and performing random circuit simulations on a virtual quantum computer.
Collapse
Affiliation(s)
- Roman Ellerbrock
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - K Grace Johnson
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Stefan Seritan
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Hannes Hoppe
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - J H Zhang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Tim Lenzen
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Thomas Weike
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Uwe Manthe
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| |
Collapse
|
21
|
Montero de Hijes P, Dellago C, Jinnouchi R, Schmiedmayer B, Kresse G. Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks. J Chem Phys 2024; 160:114107. [PMID: 38506284 DOI: 10.1063/5.0197105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 03/03/2024] [Indexed: 03/21/2024] Open
Abstract
In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests.
Collapse
Affiliation(s)
- Pablo Montero de Hijes
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- University of Vienna, Faculty of Earth Sciences, Geography and Astronomy, Josef-Holaubuek-Platz 2, 1090 Vienna, Austria
| | - Christoph Dellago
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
| | - Ryosuke Jinnouchi
- Toyota Central R&D Labs., Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | | | - Georg Kresse
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- VASP Software GmbH, Berggasse 21, A-1090 Vienna, Austria
| |
Collapse
|
22
|
Kim SS, Rhee YM. Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water. Phys Chem Chem Phys 2024; 26:9021-9036. [PMID: 38440829 DOI: 10.1039/d3cp05643k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.
Collapse
Affiliation(s)
- Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
23
|
Houston PL, Qu C, Yu Q, Pandey P, Conte R, Nandi A, Bowman JM, Kukolich SG. Formic Acid-Ammonia Heterodimer: A New Δ-Machine Learning CCSD(T)-Level Potential Energy Surface Allows Investigation of the Double Proton Transfer. J Chem Theory Comput 2024; 20:1821-1828. [PMID: 38382541 DOI: 10.1021/acs.jctc.3c01273] [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
The formic acid-ammonia dimer is an important example of a hydrogen-bonded complex in which a double proton transfer can occur. Its microwave spectrum has recently been reported and rotational constants and quadrupole coupling constants were determined. Calculated estimates of the double-well barrier and the internal barriers to rotation were also reported. Here, we report a full-dimensional potential energy surface (PES) for this complex, using two closely related Δ-machine learning methods to bring it to the CCSD(T) level of accuracy. The PES dissociates smoothly and accurately. Using a 2d quantum model the ground vibrational-state tunneling splitting is estimated to be less than 10-4 cm-1. The dipole moment along the intrinsic reaction coordinate is calculated along with a Mullikan charge analysis and supports the mildly ionic character of the minimum and strongly ionic character at the double-well barrier.
Collapse
Affiliation(s)
- Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, U.S.A. and Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Chen Qu
- Independent Researcher, Toronto, Ontario M9B0E3, Canada
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Priyanka Pandey
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, Via Golgi 19, Milano 20133, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
- Department of Physics and Materials Science, University of Luxembourg, Luxembourg City L-1511, Luxembourg
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Stephen G Kukolich
- Department of Chemistry and Biochemistry, University of Arizona, 1306 E. University Avenue, Tucson, Arizona 85721, United States
| |
Collapse
|
24
|
Liu Y, Jambrina PG, Croft JFE, Balakrishnan N, Aoiz FJ, Guo H. New Full-Dimensional Reactive Potential Energy Surface for the H 4 System. J Chem Theory Comput 2024; 20:1829-1837. [PMID: 38354106 DOI: 10.1021/acs.jctc.3c01379] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2024]
Abstract
As the most abundant molecule in the universe, collisions involving H2 have important implications in astrochemistry. Collisions between hydrogen molecules also represent a prototype for assessing various dynamic methods for understanding fundamental few-body processes. In this work, we develop a new and highly accurate full-dimensional potential energy surface (PES) covering all reactive channels of the H2 + H2 system, which extends our previously reported H2 + H2 nonreactive PES [J. Chem. Theory Comput., 2021, 17, 6747] by adding 39,538 additional ab initio points calculated at the MRCI/AV5Z level in the reactive channels. The global PES is represented with high fidelity (RMSE = 0.6 meV for a total of 79,000 points) by a permutation invariant polynomial neural network (PIP-NN) and is suitable for studying collision-induced dissociation, single-exchange, as well as four-center exchange reactions. Preliminary quasi-classical trajectory studies on the new PIP-NN PES reveal strong vibrational enhancement of all reaction channels.
Collapse
Affiliation(s)
- Yang Liu
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Pablo G Jambrina
- Departamento de Quimica Fisica, Universidad de Salamanca, Salamanca 37008, Spain
| | - James F E Croft
- The Dodd-Walls Centre for Photonic and Quantum Technologies, New Zealand and Department of Physics, University of Otago, Dunedin 9054, New Zealand
| | - Naduvalath Balakrishnan
- Department of Chemistry and Biochemistry, University of Nevada, Las Vegas, Nevada 89154, United States
| | - F Javier Aoiz
- Departamento de Quimica Fisica, Universidad Complutense, Madrid 28040, Spain
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
25
|
Stachowiak M, Grabowska E, Wang XG, Carrington T, Szalewicz K, Jankowski P. Theory cracks old data: Rovibrational energy levels of orthoH 2-CO derived from experiment. SCIENCE ADVANCES 2024; 10:eadj8632. [PMID: 38394212 PMCID: PMC10889352 DOI: 10.1126/sciadv.adj8632] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Accepted: 01/22/2024] [Indexed: 02/25/2024]
Abstract
Measurements of rovibrational spectra of clusters provide physical insight only if spectral lines can be assigned to pairs of quantum states, and further insight is obtained if one can deduce the quantitative energy-level pattern. Both steps can be so difficult that some measured spectra remain unassigned, one example is orthoH2-CO. To extend the scope of spectroscopic insights, we propose to use theoretical information in interpretation of spectra. We first performed high accuracy, full-dimensional calculations of the orthoH2-CO spectrum, at the highest practically achievable levels of electronic structure theory and quantum nuclear dynamics. Then, an iterative, theory-guided method developed here allowed us to fully interpret the spectrum of orthoH2-CO, extending the range of van der Waals clusters for which spectroscopy can provide physical insights.
Collapse
Affiliation(s)
- Marcin Stachowiak
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
| | - Ewelina Grabowska
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
| | - Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
| | - Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
| |
Collapse
|
26
|
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
|
27
|
Weike N, Eisfeld W. The effective relativistic coupling by asymptotic representation approach for molecules with multiple relativistic atoms. J Chem Phys 2024; 160:064104. [PMID: 38341788 DOI: 10.1063/5.0191529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2023] [Accepted: 01/18/2024] [Indexed: 02/13/2024] Open
Abstract
The Effective Relativistic Coupling by Asymptotic Representation (ERCAR) approach is a method to generate fully coupled diabatic potential energy surfaces (PESs) including relativistic effects, especially spin-orbit coupling. The spin-orbit coupling of a full molecule is determined only by the atomic states of selected relativistically treated atoms. The full molecular coupling effect is obtained by a diabatization with respect to asymptotic states, resulting in the correct geometry dependence of the spin-orbit effect. The ERCAR approach has been developed over the last decade and initially only for molecules with a single relativistic atom. This work presents its extension to molecules with more than a single relativistic atom using the iodine molecule as a proof-of-principle example. The theory for the general multiple atomic ERCAR approach is given. In this case, the diabatic basis is defined at the asymptote where all relativistic atoms are separated from the remaining molecular fragment. The effective spin-orbit operator is then a sum of spin-orbit operators acting on isolated relativistic atoms. PESs for the iodine molecule are developed within the new approach and it is shown that the resulting fine structure states are in good agreement with spin-orbit ab initio calculations.
Collapse
Affiliation(s)
- Nicole Weike
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
28
|
Briling K, Calvino Alonso Y, Fabrizio A, Corminboeuf C. SPA HM(a,b): Encoding the Density Information from Guess Hamiltonian in Quantum Machine Learning Representations. J Chem Theory Comput 2024; 20:1108-1117. [PMID: 38227222 PMCID: PMC10867806 DOI: 10.1021/acs.jctc.3c01040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Revised: 12/20/2023] [Accepted: 12/26/2023] [Indexed: 01/17/2024]
Abstract
Recently, we introduced a class of molecular representations for kernel-based regression methods─the spectrum of approximated Hamiltonian matrices (SPAHM)─that takes advantage of lightweight one-electron Hamiltonians traditionally used as a self-consistent field initial guess. The original SPAHM variant is built from occupied-orbital energies (i.e., eigenvalues) and naturally contains all of the information about nuclear charges, atomic positions, and symmetry requirements. Its advantages were demonstrated on data sets featuring a wide variation of charge and spin, for which traditional structure-based representations commonly fail. SPAHM(a,b), as introduced here, expand the eigenvalue SPAHM into local and transferable representations. They rely upon one-electron density matrices to build fingerprints from atomic and bond density overlap contributions inspired from preceding state-of-the-art representations. The performance and efficiency of SPAHM(a,b) is assessed on the predictions for data sets of prototypical organic molecules (QM7) of different charges and azoheteroarene dyes in an excited state. Overall, both SPAHM(a) and SPAHM(b) outperform state-of-the-art representations on difficult prediction tasks such as the atomic properties of charged open-shell species and of π-conjugated systems.
Collapse
Affiliation(s)
- Ksenia
R. Briling
- Laboratory
for Computational Molecular Design, Institute of Chemical Sciences
and Engineering, École Polytechnique
Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Yannick Calvino Alonso
- Laboratory
for Computational Molecular Design, Institute of Chemical Sciences
and Engineering, École Polytechnique
Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Alberto Fabrizio
- Laboratory
for Computational Molecular Design, Institute of Chemical Sciences
and Engineering, École Polytechnique
Fédérale de Lausanne, 1015 Lausanne, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - Clemence Corminboeuf
- Laboratory
for Computational Molecular Design, Institute of Chemical Sciences
and Engineering, École Polytechnique
Fédérale de Lausanne, 1015 Lausanne, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
29
|
Moghaddasi Fereidani R, Vaníček JJL. High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics. J Chem Phys 2024; 160:044113. [PMID: 38284658 DOI: 10.1063/5.0180070] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Accepted: 12/22/2023] [Indexed: 01/30/2024] Open
Abstract
Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but they still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multidimensional, nonseparable coupled Morse potential.
Collapse
Affiliation(s)
- Roya Moghaddasi Fereidani
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří J L Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
30
|
Iyengar SS, Ricard TC, Zhu X. Reformulation of All ONIOM-Type Molecular Fragmentation Approaches and Many-Body Theories Using Graph-Theory-Based Projection Operators: Applications to Dynamics, Molecular Potential Surfaces, Machine Learning, and Quantum Computing. J Phys Chem A 2024; 128:466-478. [PMID: 38180503 DOI: 10.1021/acs.jpca.3c05630] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2024]
Abstract
We present a graph-theory-based reformulation of all ONIOM-based molecular fragmentation methods. We discuss applications to (a) accurate post-Hartree-Fock AIMD that can be conducted at DFT cost for medium-sized systems, (b) hybrid DFT condensed-phase studies at the cost of pure density functionals, (c) reduced cost on-the-fly large basis gas-phase AIMD and condensed-phase studies, (d) post-Hartree-Fock-level potential surfaces at DFT cost to obtain quantum nuclear effects, and (e) novel transfer machine learning protocols derived from these measures. Additionally, in previous work, the unifying strategy discussed here has been used to construct new quantum computing algorithms. Thus, we conclude that this reformulation is robust and accurate.
Collapse
Affiliation(s)
- Srinivasan S Iyengar
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Timothy C Ricard
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Xiao Zhu
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
31
|
Manzhos S, Ihara M. Degeneration of kernel regression with Matern kernels into low-order polynomial regression in high dimension. J Chem Phys 2024; 160:021101. [PMID: 38189605 DOI: 10.1063/5.0187867] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Accepted: 12/17/2023] [Indexed: 01/09/2024] Open
Abstract
Kernel methods such as kernel ridge regression and Gaussian process regression with Matern-type kernels have been increasingly used, in particular, to fit potential energy surfaces (PES) and density functionals, and for materials informatics. When the dimensionality of the feature space is high, these methods are used with necessarily sparse data. In this regime, the optimal length parameter of a Matern-type kernel may become so large that the method effectively degenerates into a low-order polynomial regression and, therefore, loses any advantage over such regression. This is demonstrated theoretically as well as numerically in the examples of six- and fifteen-dimensional molecular PES using squared exponential and simple exponential kernels. The results shed additional light on the success of polynomial approximations such as PIP for medium-size molecules and on the importance of orders-of-coupling-based models for preserving the advantages of kernel methods with Matern-type kernels of on the use of physically motivated (reproducing) kernels.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|
32
|
Yu Y, Yang D, Zhou Y, Xie D. A New Full-Dimensional Ab Initio Intermolecular Potential Energy Surface and Rovibrational Energies of the H 2O-H 2 Complex. J Phys Chem A 2024; 128:170-181. [PMID: 38109882 DOI: 10.1021/acs.jpca.3c06805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
H2O-H2 is a prototypical five-atom van der Waals system, and the interaction between H2O and H2 plays an important role in many physical and chemical environments. However, previous full-dimensional intermolecular potential energy surfaces (IPESs) cannot accurately describe the H2O-H2 interaction in the repulsive or van der Waals minimum region. In this work, we constructed a full-dimensional IPES for the title system with a small root-mean-square error of 0.252 cm-1 by using the permutation invariant polynomial neural network method. The ab initio calculations were performed by employing the explicitly corrected coupled cluster [CCSD(T)-F12a] method with the augmented correlation-consistent polarized valence quintuple-ζ basis set. Based on the newly developed IPES, the bound states of the H2O-H2 complex were calculated within the rigid-rotor approximation. The transition frequencies and band origins agreed well with the experimental values [Weida, M. J.; Nesbitt, D. J. J. Chem. Phys. 1999, 110, 156-167] with errors less than 0.1 cm-1 for most transitions. Those results demonstrate the high accuracy of our new IPES, which would build a solid foundation for the collisional dynamics of H2O-H2 at low temperatures.
Collapse
Affiliation(s)
- Yipeng Yu
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Dongzheng Yang
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yanzi Zhou
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
33
|
Herringer NSM, Dasetty S, Gandhi D, Lee J, Ferguson AL. Permutationally Invariant Networks for Enhanced Sampling (PINES): Discovery of Multimolecular and Solvent-Inclusive Collective Variables. J Chem Theory Comput 2024; 20:178-198. [PMID: 38150421 DOI: 10.1021/acs.jctc.3c00923] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2023]
Abstract
The typically rugged nature of molecular free-energy landscapes can frustrate efficient sampling of the thermodynamically relevant phase space due to the presence of high free-energy barriers. Enhanced sampling techniques can improve phase space exploration by accelerating sampling along particular collective variables (CVs). A number of techniques exist for the data-driven discovery of CVs parametrizing the important large-scale motions of the system. A challenge to CV discovery is learning CVs invariant to the symmetries of the molecular system, frequently rigid translation, rigid rotation, and permutational relabeling of identical particles. Of these, permutational invariance has proved a persistent challenge in frustrating the data-driven discovery of multimolecular CVs in systems of self-assembling particles and solvent-inclusive CVs for solvated systems. In this work, we integrate permutation invariant vector (PIV) featurizations with autoencoding neural networks to learn nonlinear CVs invariant to translation, rotation, and permutation and perform interleaved rounds of CV discovery and enhanced sampling to iteratively expand the sampling of configurational phase space and obtain converged CVs and free-energy landscapes. We demonstrate the permutationally invariant network for enhanced sampling (PINES) approach in applications to the self-assembly of a 13-atom argon cluster, association/dissociation of a NaCl ion pair in water, and hydrophobic collapse of a C45H92 n-pentatetracontane polymer chain. We make the approach freely available as a new module within the PLUMED2 enhanced sampling libraries.
Collapse
Affiliation(s)
| | - Siva Dasetty
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Diya Gandhi
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Junhee Lee
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
34
|
Wang W, Liu X, Pérez-Ríos J. AlF-AlF Reaction Dynamics between 200 K and 1000 K: Reaction Mechanisms and Intermediate Complex Characterization. Molecules 2023; 29:222. [PMID: 38202805 PMCID: PMC10780286 DOI: 10.3390/molecules29010222] [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: 11/21/2023] [Revised: 12/20/2023] [Accepted: 12/22/2023] [Indexed: 01/12/2024] Open
Abstract
AlF is a relevant molecule in astrochemistry as a tracer of F-bearing molecules. Additionally, AlF presents diagonal Franck-Condon factors and can be created very efficiently in the lab, which makes it a prototypical molecular for laser cooling. However, very little is known about the reaction dynamics of AlF. In this work, we report on the reaction dynamics of AlF-AlF between 200 and 1000 K using ab initio molecular dynamics and a highly efficient active learning approach for the potential energy surface, including all degrees of freedom. As a result, we identify the main reaction mechanisms and the lifetime of the intermediate complex AlF-AlF relevant to astrochemistry environments and regions in buffer gas cells.
Collapse
Affiliation(s)
- Weiqi Wang
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany; (W.W.); (X.L.)
| | - Xiangyue Liu
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany; (W.W.); (X.L.)
| | - Jesús Pérez-Ríos
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, NY 11794-3800, USA
| |
Collapse
|
35
|
Stark W, Westermayr J, Douglas-Gallardo OA, Gardner J, Habershon S, Maurer RJ. Machine Learning Interatomic Potentials for Reactive Hydrogen Dynamics at Metal Surfaces Based on Iterative Refinement of Reaction Probabilities. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:24168-24182. [PMID: 38148847 PMCID: PMC10749455 DOI: 10.1021/acs.jpcc.3c06648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 11/12/2023] [Accepted: 11/15/2023] [Indexed: 12/28/2023]
Abstract
The reactive chemistry of molecular hydrogen at surfaces, notably dissociative sticking and hydrogen evolution, plays a crucial role in energy storage and fuel cells. Theoretical studies can help to decipher underlying mechanisms and reaction design, but studying dynamics at surfaces is computationally challenging due to the complex electronic structure at interfaces and the high sensitivity of dynamics to reaction barriers. In addition, ab initio molecular dynamics, based on density functional theory, is too computationally demanding to accurately predict reactive sticking or desorption probabilities, as it requires averaging over tens of thousands of initial conditions. High-dimensional machine learning-based interatomic potentials are starting to be more commonly used in gas-surface dynamics, yet robust approaches to generate reliable training data and assess how model uncertainty affects the prediction of dynamic observables are not well established. Here, we employ ensemble learning to adaptively generate training data while assessing model performance with full uncertainty quantification (UQ) for reaction probabilities of hydrogen scattering on different copper facets. We use this approach to investigate the performance of two message-passing neural networks, SchNet and PaiNN. Ensemble-based UQ and iterative refinement allow us to expose the shortcomings of the invariant pairwise-distance-based feature representation in the SchNet model for gas-surface dynamics.
Collapse
Affiliation(s)
- Wojciech
G. Stark
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Julia Westermayr
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | | | - James Gardner
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Scott Habershon
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Reinhard J. Maurer
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
- Department
of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| |
Collapse
|
36
|
Manzhos S, Ihara M. A controlled study of the effect of deviations from symmetry of the potential energy surface (PES) on the accuracy of the vibrational spectrum computed with collocation. J Chem Phys 2023; 159:211103. [PMID: 38038200 DOI: 10.1063/5.0182373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Accepted: 11/12/2023] [Indexed: 12/02/2023] Open
Abstract
Symmetry, in particular permutational symmetry, of a potential energy surface (PES) is a useful property in quantum chemical calculations. It facilitates, in particular, state labelling and identification of degenerate states. In many practically important applications, however, these issues are unimportant. The imposition of exact symmetry and the perception that it is necessary create additional methodological requirements narrowing or complicating algorithmic choices that are thereby biased against methods and codes that by default do not incorporate symmetry, including most off-the-shelf machine learning methods that cannot be directly used if exact symmetry is demanded. By introducing symmetric and unsymmetric errors into the PES of H2CO in a controlled way and computing the vibrational spectrum with collocation using symmetric and nonsymmetric collocation point sets, we show that when the deviations from an ideal PES are random, imposition of exact symmetry does not bring any practical advantages. Moreover, a calculation ignoring symmetry may be more accurate. We also compare machine-learned PESs with and without symmetrization and demonstrate that there is no advantage of imposing exact symmetry for the accuracy of the vibrational spectrum.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|
37
|
Fu B, Zhang DH. Accurate fundamental invariant-neural network representation of ab initio potential energy surfaces. Natl Sci Rev 2023; 10:nwad321. [PMID: 38274241 PMCID: PMC10808953 DOI: 10.1093/nsr/nwad321] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2023] [Revised: 11/01/2023] [Accepted: 11/02/2023] [Indexed: 01/27/2024] Open
Abstract
Highly accurate potential energy surfaces are critically important for chemical reaction dynamics. The large number of degrees of freedom and the intricate symmetry adaption pose a big challenge to accurately representing potential energy surfaces (PESs) for polyatomic reactions. Recently, our group has made substantial progress in this direction by developing the fundamental invariant-neural network (FI-NN) approach. Here, we review these advances, demonstrating that the FI-NN approach can represent highly accurate, global, full-dimensional PESs for reactive systems with even more than 10 atoms. These multi-channel reactions typically involve many intermediates, transition states, and products. The complexity and ruggedness of this potential energy landscape present even greater challenges for full-dimensional PES representation. These PESs exhibit a high level of complexity, molecular size, and accuracy of fit. Dynamics simulations based on these PESs have unveiled intriguing and novel reaction mechanisms, providing deep insights into the intricate dynamics involved in combustion, atmospheric, and organic chemistry.
Collapse
Affiliation(s)
- Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- Hefei National Laboratory, Hefei 230088, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- Hefei National Laboratory, Hefei 230088, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
38
|
Xia J, Zhang Y, Jiang B. Accuracy Assessment of Atomistic Neural Network Potentials: The Impact of Cutoff Radius and Message Passing. J Phys Chem A 2023; 127:9874-9883. [PMID: 37943102 DOI: 10.1021/acs.jpca.3c06024] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2023]
Abstract
Atomistic neural network potentials have achieved great success in accelerating atomistic simulations in complicated systems in recent years. They are typically based on the atomic decomposition of total properties, truncating the interatomic correlations to a local environment within a given cutoff radius. A more recently developed message passing (MP) neural network framework can, in principle, incorporate nonlocal effects through iteratively correlating some atoms outside the cutoff sphere with atoms inside, a process referred to as MP. However, how the model accuracy depends on the cutoff radius and the MP process has rarely been discussed. In this work, we investigate this dependence using a recursively embedded atom neural network method that possesses both local and MP features, in two representative systems: liquid H2O and solid Al2O3. We focus on how these settings influence predictions for structural and vibrational properties, namely, radial distribution functions (RDFs) and vibrational density of states (VDOSs). We find that while MP lowers test errors of energy and forces in general, it may not improve the prediction for RDFs and/or VDOSs if direct interatomic correlations in the local environment are insufficiently described. A cutoff radius exceeding the first neighbor shell is necessary, beyond which involving MP quickly enhances the model accuracy until convergence. This is a potentially more efficient way to increase the model accuracy than directly increasing the cutoff radius, especially with more memory savings in the GPU implementation. Our findings also suggest that using the mean test error as the measure of the model accuracy alone is inadequate.
Collapse
Affiliation(s)
- Junfan Xia
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yaolong Zhang
- École Polytechnique FFlytech de Lausanne, 1015 Lausanne, Switzerland
| | - Bin Jiang
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
39
|
Shu Y, Varga Z, Zhang D, Truhlar DG. ChemPotPy: A Python Library for Analytic Representations of Potential Energy Surfaces and Diabatic Potential Energy Matrices. J Phys Chem A 2023; 127:9635-9640. [PMID: 37916790 DOI: 10.1021/acs.jpca.3c05899] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2023]
Abstract
Constructing analytic representations of global and semiglobal potential energy surfaces is difficult and can be laborious, and it is even harder when one needs coupled potential energy surfaces and their electronically nonadiabatic couplings. When accomplished, however, the resulting potential functions are a valuable resource. To facilitate the convenient use of potentials that have been developed, we provide a collection of existing surfaces in a library with consistent units and formats. A potential energy surface library of this type, namely PotLib, was built more than 20 years ago. However, that library only provided pristine Fortran subroutines for each potential energy surface, and therefore, it is not as user-friendly as would be desirable. Here, we report the creation of ChemPotPy, a CHEMical library of POTential energy surfaces in PYthon. ChemPotPy is a user-friendly library for analytic representation of single-state and multistate potential energy surfaces and couplings. A given entry in the library contains an analytic potential energy function or analytic functions for a set of coupled potential energy surfaces, and depending on the case, it may also include analytic or numerical gradients, nonadiabatic coupling vectors, and/or diabatic potential energy matrices and their gradients. Only three inputs, namely, the chemical formula of the system, the name of the potential energy surface or surface set, and the Cartesian geometry, are required. ChemPotPy uses the same units for input and output quantities of all surfaces and surface sets to facilitate general interfaces with the dynamics programs. The initial version of the library contains 338 entries, and we anticipate that more will be added in the future.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Dayou Zhang
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
40
|
Caracciolo A, San Vicente Veliz JC, Lu D, Guo H, Meuwly M, Minton TK. Experimental and Theoretical Studies of Hyperthermal N + O 2 Collisions. J Phys Chem A 2023; 127:8834-8848. [PMID: 37843300 DOI: 10.1021/acs.jpca.3c04516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2023]
Abstract
The dynamics of hyperthermal N(4S) + O2 collisions were investigated both experimentally and theoretically. Crossed molecular beams experiments were performed at an average center-of-mass (c.m.) collision energy of ⟨Ecoll⟩ = 77.5 kcal mol-1, with velocity- and angle-resolved product detection by a rotatable mass spectrometer detector. Nonreactive (N + O2) and reactive (NO + O) product channels were identified. In the c.m. reference frame, the nonreactively scattered N atoms and reactively scattered NO molecules were both directed into the forward direction with respect to the initial direction of the reagent N atoms. On average, more than 90% of the available energy (⟨Eavl⟩ = 77.5 kcal mol-1) was retained in translation of the nonreactive products (N + O2), whereas a much smaller fraction of the available energy for the reactive pathway (⟨Eavl⟩ = 109.5 kcal mol-1) went into translation of the NO + O products, and the distribution of translational energies for this channel was broad, indicating extensive internal excitation in the nascent NO molecules. The experimentally derived c.m. translational energy and angular distributions of the reactive products suggested at least two dynamical pathways to the formation of NO + O. Quasiclassical trajectory (QCT) calculations were performed with a collision energy of Ecoll = 77 kcal mol-1 using two sets of potential energy surfaces, denoted as PES-I and PES-II, and these theoretical results were compared to each other and to the experimental results. PES-I is a reproducing kernel Hilbert space (RKHS) representation of multireference configurational interaction (MRCI) energies, while PES-II is a many-body permutation invariant polynomial (MB-PIP) fit of complete active space second order perturbation (CASPT2) points. The theoretical investigations were both consistent with the experimental suggestion of two dynamical pathways to produce NO + O, where reactive collisions may proceed on the doublet (12A') and quartet (14A') surfaces. When analyzed with this theoretical insight, the experimental c.m. translational energy and angular distributions were in reasonably good agreement with those predicted by the QCT calculations, although minor differences were observed which are discussed. Theoretical translational energy and angular distributions for the nonreactive N + O2 products matched the experimental translational energy and angular distributions almost quantitatively. Finally, relative yields for the nonreactive and reactive scattering channels were determined from the experiment and from both theoretical methods, and all results are in reasonable agreement.
Collapse
Affiliation(s)
- Adriana Caracciolo
- Ann and H.J. Smead Department of Aerospace Engineering Sciences, University of Colorado, Boulder, Colorado 80303, United States
| | | | - Dandan Lu
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Markus Meuwly
- Department of Chemistry, University of Basel, CH-4056 Basel, Switzerland
| | - Timothy K Minton
- Ann and H.J. Smead Department of Aerospace Engineering Sciences, University of Colorado, Boulder, Colorado 80303, United States
| |
Collapse
|
41
|
Feldmann R, Baiardi A, Reiher M. Symmetry-Projected Nuclear-Electronic Hartree-Fock: Eliminating Rotational Energy Contamination. J Phys Chem A 2023; 127:8943-8954. [PMID: 37831620 PMCID: PMC10614303 DOI: 10.1021/acs.jpca.3c04822] [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] [Revised: 08/31/2023] [Indexed: 10/15/2023]
Abstract
We present a symmetry projection technique for enforcing rotational and parity symmetries in nuclear-electronic Hartree-Fock wave functions, which treat electrons and nuclei on equal footing. The molecular Hamiltonian obeys rotational and parity inversion symmetries, which are, however, broken by expanding in Gaussian basis sets that are fixed in space. We generate a trial wave function with the correct symmetry properties by projecting the wave function onto representations of the three-dimensional rotation group, i.e., the special orthogonal group in three dimensions SO(3). As a consequence, the wave function becomes an eigenfunction of the angular momentum operator which (i) eliminates the contamination of the ground-state wave function by highly excited rotational states arising from the broken rotational symmetry and (ii) enables the targeting of specific rotational states of the molecule. We demonstrate the efficiency of the symmetry projection technique by calculating the energies of the low-lying rotational states of the H2 and H3+ molecules.
Collapse
Affiliation(s)
- Robin Feldmann
- ETH Zürich, Department of Chemistry
and Applied Biosciences, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| | - Alberto Baiardi
- ETH Zürich, Department of Chemistry
and Applied Biosciences, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| | - Markus Reiher
- ETH Zürich, Department of Chemistry
and Applied Biosciences, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| |
Collapse
|
42
|
Liu Y, Guo H. A Gaussian Process Based Δ-Machine Learning Approach to Reactive Potential Energy Surfaces. J Phys Chem A 2023; 127:8765-8772. [PMID: 37815868 DOI: 10.1021/acs.jpca.3c05318] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/12/2023]
Abstract
The Gaussian process (GP) is an efficient nonparametric machine learning (ML) method. A distinct advantage of the GP is its ability to provide an estimate of statistical uncertainties. This is particularly useful in constructing high-dimensional potential energy surfaces (PESs) from ab initio data as it offers an optimal way to add new geometries to reduce the overall error. In this work, GP is employed in the context of Δ-machine learning (Δ-ML), in which a correction PES to an inaccurate low-level PES is constructed using a small number of high-level ab initio calculations. This new method is tested in three prototypical reactive systems, namely, the H + H2 → H + H2, OH + H2 → H2O + H, and H + CH4 → H2 + CH3 reactions. The results show that the GP-based Δ-ML approach is more efficient than its direct application in constructing high-level PESs. We also compare the new method to a previously proposed neural-network-based Δ-ML approach [Liu and Li J. Phys. Chem. Lett. 2022, 13, 4729-4738]. The results indicate that the two Δ-ML methods have comparable efficiencies in constructing accurate PESs.
Collapse
Affiliation(s)
- Yang Liu
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
43
|
Giricz A, Czakó G, Papp D. Alternating Stereospecificity upon Central-Atom Change: Dynamics of the F - +PH 2 Cl S N 2 Reaction Compared to its C- and N-Centered Analogues. Chemistry 2023; 29:e202302113. [PMID: 37698297 DOI: 10.1002/chem.202302113] [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: 07/03/2023] [Indexed: 09/13/2023]
Abstract
Central-atom effects on bimolecular nucleophilic substitution (SN 2) reactions are well-known in chemistry, however, the atomic-level SN 2 dynamics at phosphorous (P) centers has never been studied. We investigate the dynamics of the F- +PH2 Cl reaction with the quasi-classical trajectory method on a novel full-dimensional analytical potential energy surface fitted on high-level ab initio data. Our computations reveal intermediate dynamics compared to the F- +CH3 Cl and the F- +NH2 Cl SN 2 reactions: phosphorus as central atom leads to a more indirect SN 2 reaction with extensive complex-formation with respect to the carbon-centered one, however, the title reaction is more direct than its N-centered pair. Stereospecificity, characteristic at C-center, does not appear here either, due to the submerged front-side-attack retention path and the repeated entrance-channel inversional motion, whereas the multi-inversion mechanism discovered at nitrogen center is also undermined by the deep Walden-well. At low collision energies, 6 % of the PH2 F products form with retained configuration, mostly through complex-mediated mechanisms, while this ratio reaches 24 % at the highest energy due to the increasing dominance of the direct front-side mechanism and the smaller chance for hitting the deep Walden-inversion minimum. Our results suggest pronounced central-atom effects in SN 2 reactions, which can fundamentally change their (stereo)dynamics.
Collapse
Affiliation(s)
- Anett Giricz
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary
| | - Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary
| |
Collapse
|
44
|
Liu X, Wang W, Pérez-Ríos J. Molecular dynamics-driven global potential energy surfaces: Application to the AlF dimer. J Chem Phys 2023; 159:144103. [PMID: 37811831 DOI: 10.1063/5.0169080] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2023] [Accepted: 09/20/2023] [Indexed: 10/10/2023] Open
Abstract
In this work, we present a full-dimensional potential energy surface for AlF-AlF. We apply a general machine learning approach for full-dimensional potential energy surfaces, employing an active learning scheme trained on ab initio points, whose size grows based on the accuracy required. The training points are selected based on molecular dynamics simulations, choosing the most suitable configurations for different collision energy and mapping the most relevant part of the potential energy landscape of the system. The present approach does not require long-range information and is entirely general. As a result, it is possible to provide the full-dimensional AlF-AlF potential energy surface, requiring ≲0.01% of the configurations to be calculated ab initio. Furthermore, we analyze the general properties of the AlF-AlF system, finding critical differences with other reported results on CaF or bi-alkali dimers.
Collapse
Affiliation(s)
- Xiangyue Liu
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Weiqi Wang
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Jesús Pérez-Ríos
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, USA
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794-3800, USA
| |
Collapse
|
45
|
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
|
46
|
Herman KM, Stone AJ, Xantheas SS. Accurate Calculation of Many-Body Energies in Water Clusters Using a Classical Geometry-Dependent Induction Model. J Chem Theory Comput 2023; 19:6805-6815. [PMID: 37703063 DOI: 10.1021/acs.jctc.3c00575] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/14/2023]
Abstract
We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The moment expansion is carried out up to the hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA classical interaction potentials for water. The classical induction model with no adjustable parameters reproduces the ab initio 3-/4-body terms contained in these two databases with a root-mean-square error (RMSE) of 0.104/0.058 and a mean-absolute error (MAE) of 0.054/0.026 kcal/mol, respectively. These results are on par with the ones obtained by fitting the same data using over 14,000 (for the 3-body) and 200 (for the 4-body) parameters via Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy, included in the calculation of the 3-body energy, was found to be small but not quite negligible. The model represents a practical, efficient, and transferable approach for obtaining accurate nonadditive interactions for multicomponent systems without the need to perform tens of thousands of high level electronic structure calculations and fitting them with PIPs.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98185, United States
| | - Anthony J Stone
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, U.K
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98185, United States
- Advanced Computing Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MSIN J7-10, Richland, Washington 99352, United States
| |
Collapse
|
47
|
Szűcs T, Czakó G. ManyHF-based full-dimensional potential energy surface development and quasi-classical dynamics for the Cl + CH3NH2 reaction. J Chem Phys 2023; 159:134306. [PMID: 37791624 DOI: 10.1063/5.0166680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2023] [Accepted: 09/06/2023] [Indexed: 10/05/2023] Open
Abstract
A full-dimensional spin-orbit (SO)-corrected potential energy surface (PES) is developed for the Cl + CH3NH2 multi-channel system. Using the new PES, a comprehensive reaction dynamics investigation is performed for the most reactive hydrogen-abstraction reactions forming HCl + CH2NH2/CH3NH. Hartree-Fock (HF) convergence problems in the reactant region are handled by the ManyHF method, which finds the lowest-energy HF solution considering several different initial guess orbitals. The PES development is carried out with the Robosurfer program package, which iteratively improves the surface. Energy points are computed at the ManyHF-UCCSD(T)-F12a/cc-pVDZ-F12 level of theory combined with basis set (ManyHF-RMP2-F12/cc-pVTZ-F12 - ManyHF-RMP2-F12/cc-pVDZ-F12) and SO (MRCI+Q/aug-cc-pwCVDZ) corrections. Quasi-classical trajectory simulations show that the CH3-side hydrogen abstraction occurs more frequently in contrast to the NH2-side reaction. In both cases, the integral cross sections decrease with increasing collision energy (Ecoll). A reaction mechanism shifting from indirect to direct stripping can be observed from the opacity functions, scattering angle, and translation energy distributions as Ecoll increases. Initial attack angle distributions reveal that chlorine prefers to abstract hydrogen from the approached functional group. The collision-energy dependence of the product energy distributions shows that the initial translational energy mainly transfers to product recoil. The HCl vibrational and rotational energy values are comparable and nearly independent of collision energy, while the CH2NH2 and CH3NH co-products' vibrational energy values are higher than the rotational energy values with more significant Ecoll dependence. The HCl(v = 0) rotational distributions are compared with experiment, setting the direction for future investigations.
Collapse
Affiliation(s)
- Tímea Szűcs
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
48
|
Chen J, Zhang H, Zhou L, Hu X, Xie D. New accurate diabatic potential energy surfaces for the two lowest 1A'' states of H 2S and photodissociation dynamics in its first absorption band. Phys Chem Chem Phys 2023; 25:26032-26042. [PMID: 37750311 DOI: 10.1039/d3cp03026a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
In this work, state-to-state photodissociation dynamics of H2S in its first absorption band has been studied quantum mechanically with a new set of coupled potential energy surfaces (PESs) for the first two 1A'' excited states, which were developed at the explicitly correlated internally contracted multi-reference configuration interaction level with the cc-pVQZ-F12 basis set and a large active space. The calculated absorption spectrum, product state distributions, and angular distributions are in excellent agreement with available experimental data, validating the accuracy of the PESs and the non-adiabatic couplings. Detailed analysis of the dynamics reveals that there are strong non-adiabatic couplings between the bound 11B1 and dissociative 11A2 states around the Franck-Condon region, leading to very fast predissociation to ro-vibrationally cold SH(X̃) fragments, during which marginal angular anisotropy of the PESs is involved. This study provides quantitatively accurate characterization of the electronic structure and detailed fragmentation dynamics of this prototypical photodissociation system, which is desirable for improving astrochemical modelling.
Collapse
Affiliation(s)
- Junjie Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Hanzi Zhang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Linsen Zhou
- Institute of Materials, China Academy of Engineering Physics, Mianyang 621907, China.
| | - Xixi Hu
- Kuang Yaming Honors School, Institute for Brain Sciences, Nanjing University, Nanjing 210023, China.
- Hefei National Laboratory, Hefei 230088, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
49
|
Song H, Guo H. Theoretical Insights into the Dynamics of Gas-Phase Bimolecular Reactions with Submerged Barriers. ACS PHYSICAL CHEMISTRY AU 2023; 3:406-418. [PMID: 37780541 PMCID: PMC10540288 DOI: 10.1021/acsphyschemau.3c00009] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 06/02/2023] [Accepted: 06/06/2023] [Indexed: 10/03/2023]
Abstract
Much attention has been paid to the dynamics of both activated gas-phase bimolecular reactions, which feature monotonically increasing integral cross sections and Arrhenius kinetics, and their barrierless capture counterparts, which manifest monotonically decreasing integral cross sections and negative temperature dependence of the rate coefficients. In this Perspective, we focus on the dynamics of gas-phase bimolecular reactions with submerged barriers, which often involve radicals or ions and are prevalent in combustion, atmospheric chemistry, astrochemistry, and plasma chemistry. The temperature dependence of the rate coefficients for such reactions is often non-Arrhenius and complex, and the corresponding dynamics may also be quite different from those with significant barriers or those completely dominated by capture. Recent experimental and theoretical studies of such reactions, particularly at relatively low temperatures or collision energies, have revealed interesting dynamical behaviors, which are discussed here. The new knowledge enriches our understanding of the dynamics of these unusual reactions.
Collapse
Affiliation(s)
- Hongwei Song
- State
Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science
and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University
of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
50
|
Manzhos S, Ihara M. Neural Network with Optimal Neuron Activation Functions Based on Additive Gaussian Process Regression. J Phys Chem A 2023; 127:7823-7835. [PMID: 37698519 DOI: 10.1021/acs.jpca.3c02949] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/13/2023]
Abstract
Feed-forward neural networks (NNs) are a staple machine learning method widely used in many areas of science and technology, including physical chemistry, computational chemistry, and materials informatics. While even a single-hidden-layer NN is a universal approximator, its expressive power is limited by the use of simple neuron activation functions (such as sigmoid functions) that are typically the same for all neurons. More flexible neuron activation functions would allow the use of fewer neurons and layers and thereby save computational cost and improve expressive power. We show that additive Gaussian process regression (GPR) can be used to construct optimal neuron activation functions that are individual to each neuron. An approach is also introduced that avoids nonlinear fitting of neural network parameters by defining them with rules. The resulting method combines the advantage of robustness of a linear regression with the higher expressive power of an NN. We demonstrate the approach by fitting the potential energy surfaces of the water molecule and formaldehyde. Without requiring any nonlinear optimization, the additive-GPR-based approach outperforms a conventional NN in the high-accuracy regime, where a conventional NN suffers more from overfitting.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|