1
|
Song Q, Zhang X, Miao Z, Meng Q. Construction of a Mode-Combination Hamiltonian under the Grid-Based Representation for the Quantum Dynamics of OH + HO 2 → O 2 + H 2O. J Chem Theory Comput 2024; 20:597-613. [PMID: 38199964 DOI: 10.1021/acs.jctc.3c01090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
In this work, a systematic construction framework on a mode-combination Hamiltonian operator of a typical polyatomic reaction, OH + HO2 → O2 + H2O, is developed. First, a set of Jacobi coordinates are employed to construct the kinetic energy operator (KEO) through the polyspherical approach ( Phys. Rep. 2009, 484, 169). Second, due to the multiconfigurational electronic structure of this system, a non-adiabatic potential energy surface (PES) is constructed where the first singlet and triplet states are involved with spin-orbital coupling. To improve the training database, the training set of random energy data was optimized through a popular iterative optimization approach with extensive trajectories. Here, we propose an automatic trajectory method, instead of the classical trajectory on a crude PES, where the gradients are directly computed by the present ab initio calculations. Third, on the basis of the training set, the potential function is directly constructed in the canonical polyadic decomposition (CPD) form ( J. Chem. Theory Comput. 2021, 17, 2702-2713) which is helpful in propagating the nuclear wave function under the grid-based representation. To do this, the Gaussian process regression (GPR) approach for building the CPD form, called the CPD-GPR method ( J. Phys. Chem. Lett. 2022, 13, 11128-11135) is adopted where we further revise CPD-GPR by introducing the mode-combination (mc) scheme leading to the present CPD-mc-GPR approach. Constructing the full-dimension non-adiabatic Hamiltonian operator with mode combination, as test calculations, the nuclear wave function is propagated to preliminarily compute the reactive probability of OH + HO2 → O2 + H2O where the reactants are prepared in vibrational ground states and in the first triplet electronic state.
Collapse
Affiliation(s)
- Qingfei Song
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| | - Zekai Miao
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| | - Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| |
Collapse
|
2
|
Miao Z, Zhang X, Zhang Y, Wang L, Meng Q. Chemistry-Informed Generative Model for Classical Dynamics Simulations. J Phys Chem Lett 2024; 15:532-539. [PMID: 38194494 DOI: 10.1021/acs.jpclett.3c03114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2024]
Abstract
In this work, a chemistry-informed generative model was proposed, leading to the chemistry-informed generative adversarial network (CI-GAN) approach. To easily build the input database for complex molecular systems, an image-input algorithm is also implemented, leading to the capability to directly recognize the molecular image. Extensive test calculations and analysis on typical examples, H + H2, OH + HO2, and H2O/TiO2(110), find that the present CI-GAN approach generates distributions of geometry and energy. Calculations on the above examples show that the present CI-GAN approach is able to generate 50%-80% meaningful results among all of the generated data with chemistry constraints. Thus, it has the potential capability to predict classical dynamics simulations as well as ab initio calculations avoiding expensive calculations. These results and the power of CI-GANs in generating ab initio energies and MD trajectories are deeply discussed.
Collapse
Affiliation(s)
- Zekai Miao
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Yuyuan Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Lemei Wang
- Ministry-of-Education Engineering Center for Embedded System Integration, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| |
Collapse
|
3
|
Konings M, Harvey JN, Loreau J. Machine Learning Representations of the Three Lowest Adiabatic Electronic Potential Energy Surfaces for the ArH 2+ Reactive System. J Phys Chem A 2023; 127:8083-8094. [PMID: 37748085 DOI: 10.1021/acs.jpca.3c04015] [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, we present Gaussian process regression machine learning representations of the three lowest coupled 2A' adiabatic electronic potential energy surfaces of the ArH2+ reactive system in full dimensionality. Additionally, the nonadiabatic coupling matrix elements were calculated. These adiabatic potentials and their nonadiabatic couplings are necessary ingredients in the theoretical investigation of the nonadiabatic reaction dynamics of the Ar + H2+ → ArH+ + H and Ar+ + H2 → ArH+ + H reactions, as well as the competing charge transfer process, Ar + H2+↔ Ar+ + H2. Accurate ab initio electronic structure calculations (ic-MRCI+Q/aug-cc-pVQZ), whereby the effect of spin-orbit coupling in Ar+ has been accounted for through the state interaction method, serve as input for the machine learning training process. The potential energy surfaces are fitted with high accuracies, with root-mean-square errors on the order of 10-7 eV for the three surfaces, which meet the requirements for chemical dynamics at low temperature. It was found that quite a large number of training points (of the order of 5000 ab initio points) are needed in order to achieve these accuracies due to the complex topography of these electronic surfaces.
Collapse
Affiliation(s)
- Maarten Konings
- Division of Quantum Chemistry and Physical Chemistry, Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Jeremy N Harvey
- Division of Quantum Chemistry and Physical Chemistry, Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Jérôme Loreau
- Division of Quantum Chemistry and Physical Chemistry, Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| |
Collapse
|
4
|
Speak T, Blitz MA, Medeiros DJ, Seakins PW. New Measurements and Calculations on the Kinetics of an Old Reaction: OH + HO 2 → H 2O + O 2. JACS AU 2023; 3:1684-1694. [PMID: 37388696 PMCID: PMC10301680 DOI: 10.1021/jacsau.3c00110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 05/15/2023] [Accepted: 05/15/2023] [Indexed: 07/01/2023]
Abstract
Literature rate coefficients for the prototypical radical-radical reaction at 298 K vary by close to an order of magnitude; such variations challenge our understanding of fundamental reaction kinetics. We have studied the title reaction at room temperature via the use of laser flash photolysis to generate OH and HO2 radicals, monitoring OH by laser-induced fluorescence using two different approaches, looking at the direct reaction and also the perturbation of the slow OH + H2O2 reaction with radical concentration, and over a wide range of pressures. Both approaches give a consistent measurement of k1,298K ∼1 × 10-11 cm3 molecule-1 s-1, at the lowest limit of previous determinations. We observe, experimentally, for the first time, a significant enhancement in the rate coefficient in the presence of water, k1,H2O, 298K = (2.17 ± 0.09) × 10-28 cm6 molecule-2 s-1, where the error is statistical at the 1σ level. This result is consistent with previous theoretical calculations, and the effect goes some way to explaining some, but not all, of the variation in previous determinations of k1,298K. Supporting master equation calculations, using calculated potential energy surfaces at the RCCSD(T)-F12b/CBS//RCCSD/aug-cc-pVTZ and UCCSD(T)/CBS//UCCSD/aug-cc-pVTZ levels, are in agreement with our experimental observations. However, realistic variations in barrier heights and transition state frequencies give a wide range of calculated rate coefficients showing that the current precision and accuracy of calculations are insufficient to resolve the experimental discrepancies. The lower value of k1,298K is consistent with experimental observations of the rate coefficient of the related reaction, Cl + HO2 → HCl + O2. The implications of these results in atmospheric models are discussed.
Collapse
Affiliation(s)
- Thomas
H. Speak
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, U.K.
| | - Mark A. Blitz
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, U.K.
- National
Centre for Atmospheric Science, University
of Leeds, Leeds LS2 9JT, U.K.
| | | | - Paul W. Seakins
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, U.K.
| |
Collapse
|
5
|
Yang Z, Chen H, Buren B, Chen M. Globally Accurate Gaussian Process Potential Energy Surface and Quantum Dynamics Studies on the Li(2S) + Na2 → LiNa + Na Reaction at Low Collision Energies. Molecules 2023; 28:molecules28072938. [PMID: 37049701 PMCID: PMC10096016 DOI: 10.3390/molecules28072938] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Revised: 03/22/2023] [Accepted: 03/23/2023] [Indexed: 03/29/2023] Open
Abstract
The LiNa2 reactive system has recently received great attention in the experimental study of ultracold chemical reactions, but the corresponding theoretical calculations have not been carried out. Here, we report the first globally accurate ground-state LiNa2 potential energy surface (PES) using a Gaussian process model based on only 1776 actively selected high-level ab initio training points. The constructed PES had high precision and strong generalization capability. On the new PES, the quantum dynamics calculations on the Li(2S) + Na2(v = 0, j = 0) → LiNa + Na reaction were carried out in the 0.001–0.01 eV collision energy range using an improved time-dependent wave packet method. The calculated results indicate that this reaction is dominated by a complex-forming mechanism at low collision energies. The presented dynamics data provide guidance for experimental research, and the newly constructed PES could be further used for ultracold reaction dynamics calculations on this reactive system.
Collapse
|
6
|
Song Q, Zhang X, Peláez D, Meng Q. Direct Canonical-Polyadic-Decomposition of the Potential Energy Surface from Discrete Data by Decoupled Gaussian Process Regression. J Phys Chem Lett 2022; 13:11128-11135. [PMID: 36442084 DOI: 10.1021/acs.jpclett.2c03080] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
A Gaussian process regression (GPR) approach for directly constructing the canonical polyadic decomposition (CPD) of a multidimensional potential energy surface (PES) by discrete training energies is proposed and denoted by CPD-GPR. The present CPD-GPR method requires the kernel function in a product of a series of one-dimensional functions. To test CPD-GPR, the reactive probabilities of H + H2 as a function of kinetics energy are performed. Comparing the dynamics results computed by the CPD-GPR PES with those by the original PES, a good agreement between these results can be clearly found. Discussions on the previous algorithms for building the decomposed form are also given. We further show that the CPD-GPR method might be the general algorithm for building the decomposed form. However, further development is needed to reduce the CPD rank. Therefore, the present CPD-GPR method might be helpful to inspire ideas for developing new tools in building decomposed potential functions.
Collapse
Affiliation(s)
- Qingfei Song
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072Xi'an, China
- Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405Orsay, France
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072Xi'an, China
| | - Daniel Peláez
- Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405Orsay, France
| | - Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072Xi'an, China
| |
Collapse
|
7
|
Dai J, Krems RV. Quantum Gaussian process model of potential energy surface for a polyatomic molecule. J Chem Phys 2022; 156:184802. [PMID: 35568545 DOI: 10.1063/5.0088821] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
With gates of a quantum computer designed to encode multi-dimensional vectors, projections of quantum computer states onto specific qubit states can produce kernels of reproducing kernel Hilbert spaces. We show that quantum kernels obtained with a fixed ansatz implementable on current quantum computers can be used for accurate regression models of global potential energy surfaces (PESs) for polyatomic molecules. To obtain accurate regression models, we apply Bayesian optimization to maximize marginal likelihood by varying the parameters of the quantum gates. This yields Gaussian process models with quantum kernels. We illustrate the effect of qubit entanglement in the quantum kernels and explore the generalization performance of quantum Gaussian processes by extrapolating global six-dimensional PESs in the energy domain.
Collapse
Affiliation(s)
- J Dai
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, CanadaStewart Blusson Quantum Matter Institute, Vancouver, British Columbia V6T 1Z4, Canada
| | - R V Krems
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, CanadaStewart Blusson Quantum Matter Institute, Vancouver, British Columbia V6T 1Z4, Canada
| |
Collapse
|
8
|
Yang Z, Chen H, Chen M. Representing Globally Accurate Reactive Potential Energy Surfaces with Complex Topography by Combining Gaussian Process Regression and Neural Network. Phys Chem Chem Phys 2022; 24:12827-12836. [DOI: 10.1039/d2cp00719c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
There has been increasing attention in using machine learning technologies, such as neural network (NN) and Gaussian process regression (GPR), to model multidimensional potential energy surfaces (PESs). NN PES features...
Collapse
|
9
|
Ceriotti M, Clementi C, Anatole von Lilienfeld O. Machine learning meets chemical physics. J Chem Phys 2021; 154:160401. [PMID: 33940847 DOI: 10.1063/5.0051418] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Over recent years, the use of statistical learning techniques applied to chemical problems has gained substantial momentum. This is particularly apparent in the realm of physical chemistry, where the balance between empiricism and physics-based theory has traditionally been rather in favor of the latter. In this guest Editorial for the special topic issue on "Machine Learning Meets Chemical Physics," a brief rationale is provided, followed by an overview of the topics covered. We conclude by making some general remarks.
Collapse
Affiliation(s)
- Michele Ceriotti
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Cecilia Clementi
- Department of Physics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | | |
Collapse
|
10
|
Liu Y, Song H, Li J. Kinetic study of the OH + HO 2 → H 2O + O 2 reaction using ring polymer molecular dynamics and quantum dynamics. Phys Chem Chem Phys 2020; 22:23657-23664. [PMID: 33112305 DOI: 10.1039/d0cp04120c] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reaction OH + HO2 → H2O + O2 is a prototype of radical-radical reactions. It plays an important role in interstellar/atmospheric chemistry and combustion, and considerable attention has thus been dedicated to its kinetics. In our previous work, we reported an accurate full-dimensional potential energy surface for the title reaction on the ground triplet electronic state. The quasi-classical trajectory (QCT) approach was employed to investigate its kinetics. Although the QCT rate coefficients were in good agreement with some experimental and theoretical results, QCT cannot account for the quantum mechanical effects, such as zero-point vibrational energy, recrossing, and tunneling, which may significantly affect the rate coefficients, particularly at low temperatures. In this work, the reduced-dimensional quantum dynamics and ring polymer molecular dynamics calculations were carried out to examine these effects and their impact on rate coefficients over the temperature range of 300-1300 K.
Collapse
Affiliation(s)
- Yang Liu
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| | - Hongwei Song
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| |
Collapse
|
11
|
Bahlke MP, Mogos N, Proppe J, Herrmann C. Exchange Spin Coupling from Gaussian Process Regression. J Phys Chem A 2020; 124:8708-8723. [DOI: 10.1021/acs.jpca.0c05983] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Marc Philipp Bahlke
- Department of Chemistry, University of Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Natnael Mogos
- Department of Chemistry, University of Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Jonny Proppe
- Institute of Physical Chemistry, Georg-August University, Tammannstr. 6, 37077 Göttingen, Germany
| | - Carmen Herrmann
- Department of Chemistry, University of Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| |
Collapse
|
12
|
Yang W, Peng L, Zhu Y, Hong L. When machine learning meets multiscale modeling in chemical reactions. J Chem Phys 2020; 153:094117. [DOI: 10.1063/5.0015779] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Wuyue Yang
- Yau Mathematical Sciences Center, Tsinghua University, Beijing 100084, People’s Republic of China
| | - Liangrong Peng
- College of Mathematics and Data Science, Minjiang University, Fuzhou 350108, People’s Republic of China
| | - Yi Zhu
- Yau Mathematical Sciences Center, Tsinghua University, Beijing 100084, People’s Republic of China
| | - Liu Hong
- School of Mathematics, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| |
Collapse
|