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, 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
|
3
|
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
|
4
|
Manzhos S, Ihara M. Computational vibrational spectroscopy of molecule-surface interactions: what is still difficult and what can be done about it. Phys Chem Chem Phys 2022; 24:15158-15172. [DOI: 10.1039/d2cp01389d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Interactions of molecules with solid surfaces are responsible for key functionalities for a range of currently actively pursued technologies, including heterogeneous catalysis for synthesis or decomposition of molecules, sensitization, surface...
Collapse
|
5
|
Meng Q, Schröder M, Meyer HD. High-Dimensional Quantum Dynamics Study on Excitation-Specific Surface Scattering Including Lattice Effects of a Five-Atom Surface Cell. J Chem Theory Comput 2021; 17:2702-2713. [PMID: 33904716 DOI: 10.1021/acs.jctc.1c00241] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, high-dimensional (21D) quantum dynamics calculations on the mode-specific surface scattering of a carbon monoxide molecule on a copper(100) surface with lattice effects of a five-atom surface cell are performed through the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. We employ a surface model in which five surface atoms near the impact site are treated as fully flexible quantum particles, while all other more distant atoms are kept at fixed locations. To efficiently perform the 21D ML-MCTDH wave packet propagation, the potential energy surface is transferred to a canonical polyadic decomposition form with the aid of a Monte Carlo-based method. Excitation-specific sticking probabilities of CO on Cu(100) are computed, and lattice effects caused by the flexible surface atoms are demonstrated by comparison with sticking probabilities computed for a rigid surface. The dependence of the sticking probability of the initial state of the system is studied, and it is found that the sticking probability is reduced when the surface atom on the impact site is initially vibrationally excited.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Markus Schröder
- Theoretische Chemie, Physikalisch-Chemisches Institut, Ruprecht-Karls Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Ruprecht-Karls Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
6
|
Song Q, Zhang Q, Meng Q. Neural-network potential energy surface with small database and high precision: A benchmark of the H + H2 system. J Chem Phys 2019; 151:114302. [PMID: 31542037 DOI: 10.1063/1.5118692] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingfei Song
- Department of Applied Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Ministry-of-Education Key Laboratory of Materials Physics and Chemistry under Extraordinary Conditions, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| | - Qiuyu Zhang
- Department of Applied Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Ministry-of-Education Key Laboratory of Materials Physics and Chemistry under Extraordinary Conditions, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| | - Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Ministry-of-Education Key Laboratory of Materials Physics and Chemistry under Extraordinary Conditions, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| |
Collapse
|
7
|
Meng Q. Ring-Polymer Molecular Dynamics with Coarse-Grained Treatment of the Rate Coefficients of Chlorine Atom Reactions with Methane, Ethane, and Propane. J Phys Chem A 2018; 122:8320-8325. [PMID: 30281311 DOI: 10.1021/acs.jpca.8b08052] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, 710072 Xi’an, China
- State Key Laboratory
of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| |
Collapse
|
8
|
Meng Q, Meyer HD. Lattice effects of surface cell: Multilayer multiconfiguration time-dependent Hartree study on surface scattering of CO/Cu(100). J Chem Phys 2017. [DOI: 10.1063/1.4982962] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, 710072 Xi’an, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Ruprecht-Karls Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|