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
|
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
|
3
|
Song Q, Zhang X, Gatti F, Miao Z, Zhang Q, Meng Q. Multilayer Multiconfiguration Time-Dependent Hartree Study on the Mode-/Bond-Specific Quantum Dynamics of Water Dissociation on Cu(111). J Phys Chem A 2022; 126:6047-6058. [PMID: 36054932 DOI: 10.1021/acs.jpca.2c03092] [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/28/2022]
Abstract
In this work, full-dimensional (9D) quantum dynamics calculations on mode-/bond-specific surface scattering of a water molecule on a copper (111) rigid surface are performed through the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. To easily perform the ML-MCTDH calculations on such a triatomic molecule-surface system, we first choose specific Jacobi coordinates as a set of coordinates of water. Next, to efficiently perform the 9D ML-MCTDH wavepacket propagation, the potential energy surface is transferred to a canonical polyadic decomposition form with the aid of a Monte Carlo-based method. Excitation-specific dissociation probabilities of H2O on Cu(111) are computed, and mode-/bond-specific dynamics are demonstrated by comparison with a probability curve computed for a water molecule in the ground state. The dependence of the dissociation probability of the initial state of H2O is studied, and it is found that the excitation-specific dissociation probabilities can be divided into three groups. We find that the vibrationally excited states enhance the dissociation reactivity of H2O, while the rotationally excited states hardly influence it.
Collapse
Affiliation(s)
- Qingfei Song
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.,Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405 Orsay, France
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Fabien Gatti
- Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405 Orsay, France
| | - Zekai Miao
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Qiuyu Zhang
- Department of Chemistry, 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
|
4
|
Han S, Schröder M, Gatti F, Meyer HD, Lauvergnat D, Yarkony DR, Guo H. Representation of Diabatic Potential Energy Matrices for Multiconfiguration Time-Dependent Hartree Treatments of High-Dimensional Nonadiabatic Photodissociation Dynamics. J Chem Theory Comput 2022; 18:4627-4638. [PMID: 35839299 DOI: 10.1021/acs.jctc.2c00370] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Conventional quantum mechanical characterization of photodissociation dynamics is restricted by steep scaling laws with respect to the dimensionality of the system. In this work, we examine the applicability of the multi-configurational time-dependent Hartree (MCTDH) method in treating nonadiabatic photodissociation dynamics in two prototypical systems, taking advantage of its favorable scaling laws. To conform to the sum-of-product form, elements of the ab initio diabatic potential energy matrix (DPEM) are re-expressed using the recently proposed Monte Carlo canonical polyadic decomposition method, with enforcement of proper symmetry. The MCTDH absorption spectra and product branching ratios are shown to compare well with those calculated using conventional grid-based methods, demonstrating its promise for treating high-dimensional nonadiabatic photodissociation problems.
Collapse
Affiliation(s)
- Shanyu Han
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Markus Schröder
- Theoretische Chemie, Physikalisch Chemisches Institut, Ruprecht-Karls Universität Heidelberg, D-69120 Heidelberg, Germany
| | - Fabien Gatti
- ISMO, Institut des Sciences Moléculaires d'Orsay─UMR 8214 CNRS/Université Paris-Saclay, F-91405 Orsay, France
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch Chemisches Institut, Ruprecht-Karls Universität Heidelberg, D-69120 Heidelberg, Germany
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
5
|
Meng Q, Chen J, Ma J, Zhang X, Chen J. Adiabatic models for the quantum dynamics of surface scattering with lattice effects. Phys Chem Chem Phys 2022; 24:16415-16436. [PMID: 35766107 DOI: 10.1039/d2cp01560a] [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
In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Junbo Chen
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China. .,Xi'an Modern Chemistry Research Institute, China North Industries Group Corp., Ltd., East Zhangba Road 168, 710065 Xi'an, China
| | - Jianxing Ma
- 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.
| | - Jun Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Yangqiao Road West 155, 350002 Fuzhou, China.,Fujian Science & Technology Innovation Laboratory for Optoelectronic Information of China, Optoelectronic Industry Base at High-tech Zone, 350108 Fuzhou, China
| |
Collapse
|