1
|
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
|
2
|
Pacifici L, Verdicchio M, Lago NF, Lombardi A, Costantini A. A high-level ab initio study of the N2 + N2 reaction channel. J Comput Chem 2013; 34:2668-76. [PMID: 24037708 DOI: 10.1002/jcc.23415] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2013] [Revised: 07/26/2013] [Accepted: 07/27/2013] [Indexed: 11/08/2022]
Abstract
A new six-dimensional (6D) global potential energy surface (PES) is proposed for the full range description of the interaction of the N2(1Σg+)+N2(1Σg+) system governing collisional processes, including N atom exchange. The related potential energy values were determined using high-level ab initio methods. The calculations were performed at a coupled-cluster with single and double and perturbative triple excitations level of theory in order to have a first full range picture of the PES. Subsequently, in order to accurately describe the stretching of the bonds of the two interacting N2 molecules by releasing the constraints of being considered as rigid rotors, for the same molecular geometries higher level of theory multi reference calculations were performed. Out of the calculated values a 6D 4-atoms global PES was produced for use in dynamical calculations. The ab initio calculations were made possible by the combined use of High Throughput Computing and High Performance Computing techniques within the frame of a computing grid empowered molecular simulator.
Collapse
Affiliation(s)
- Leonardo Pacifici
- Department of Chemistry, University of Perugia, via Elce di Sotto, 8, Perugia, 06123, Italy
| | | | | | | | | |
Collapse
|
3
|
Song H, Lu Y, Lee SY. Fully converged integral cross sections of collision induced dissociation, four-center, and single exchange reactions, and accuracy of the centrifugal sudden approximation in H2 + D2 reaction. J Chem Phys 2012; 136:114307. [DOI: 10.1063/1.3693762] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
4
|
Song H, Lu Y, Lee SY. Full-dimensional time-dependent wave packet dynamics of H2 + D2 reaction. J Chem Phys 2011; 135:014305. [PMID: 21744900 DOI: 10.1063/1.3607965] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Collision induced dissociation (CID), four center reaction (4C), and single exchange reaction (SE) in H(2) (v(1) = high) + D(2) (v(2) = low) were studied by means of time-dependent wave packet approach within a full-dimensional model. Initial state-selected total reaction probabilities for the three competitive processes have been computed on two realistic global potential energy surfaces of Aguado-Suárez-Paniagua and Boothroyd-Martin-Keogh-Peterson (BMKP) with the total angular momentum J = 0. The role of both vibrationally excited and rotationally excited reagents was examined by varying the initial vibrational and rotational states. The vibrational excitation of the hot diatom gives an enhancement effect on the CID process, while the vibrational excitation of the cold diatom gives an inhibition effect. The rotational excitation of both reagents has a significant effect on the reaction process. The 4C and SE probabilities are at least one order of magnitude smaller than the CID probabilities over the energy range considered. Isotope substitution effects were also studied by substituting the collider D(2) by H(2) and HD on the BMKP potential energy surfaces. The CID process is most efficient for the H(2) + D(2) combination and least efficient for the H(2) + H(2) combination and is different for the 4C and SE processes.
Collapse
Affiliation(s)
- Hongwei Song
- Division of Physics and Applied Physics, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371, Singapore
| | | | | |
Collapse
|
5
|
|
6
|
Garcia E, Saracibar A, Sánchez C, Laganà A. Effect of the Total Angular Momentum on the Dynamics of the H2 + H2 System. J Phys Chem A 2009; 113:14312-20. [DOI: 10.1021/jp903072j] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Ernesto Garcia
- Departamento de Quimica Fisica, Universidad del Pais Vasco, 01006 Vitoria, Spain, and Dipartimento di Chimica, Università di Perugia, 06123 Perugia, Italy
| | - Amaia Saracibar
- Departamento de Quimica Fisica, Universidad del Pais Vasco, 01006 Vitoria, Spain, and Dipartimento di Chimica, Università di Perugia, 06123 Perugia, Italy
| | - Carlos Sánchez
- Departamento de Quimica Fisica, Universidad del Pais Vasco, 01006 Vitoria, Spain, and Dipartimento di Chimica, Università di Perugia, 06123 Perugia, Italy
| | - Antonio Laganà
- Departamento de Quimica Fisica, Universidad del Pais Vasco, 01006 Vitoria, Spain, and Dipartimento di Chimica, Università di Perugia, 06123 Perugia, Italy
| |
Collapse
|
7
|
Carmona-Novillo E, Bartolomei M, Hernández MI, Campos-Martínez J. Quasiclassical trajectory study of reactive and dissociative processes in H2+H2: Comparison with quantum-mechanical calculations. J Chem Phys 2007; 126:124315. [PMID: 17411131 DOI: 10.1063/1.2712842] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quasiclassical trajectory calculations have been carried out for H(2)(v(1)=high)+H(2)(v(2)=low) collisions within a three degrees of freedom model where five different geometries of the colliding complex were considered. Within this approach, probabilities for different competitive processes are studied: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. The purpose is to compare in detail with equivalent quantum-mechanical wave packet calculations [Bartolomei et al., J. Chem. Phys 122, 064305 (2005)], especially the behavior of the probabilities near reaction thresholds. Quasiclassical calculations compare quite well with the quantum-mechanical ones for collision induced dissociation as well as for the four center reaction, although quantum effects become very important near thresholds, particularly for lower v(1)'s and for the four center process. Less quantitative agreement is found for reactive dissociation and three-body complex formation. It is found that most quantum effects are due to differences between quantum and classical vibrational distributions of H(2)(v(1)=high). Zero point energy violation has been found in the classical reactive-dissociative probabilities. Extension of these findings to full-dimensional treatments is examined.
Collapse
Affiliation(s)
- Estela Carmona-Novillo
- Instituto de Matemáticas y Física Fundamental CSIC, C/ Serrano 123, E-28006 Madrid, Spain
| | | | | | | |
Collapse
|
8
|
Bartolomei M, Hernández MI, Campos-Martínez J. Wave packet dynamics of H2(v1=8-14)+H2(v2=0-2): the role of the potential energy surface on different reactive and dissociative processes. J Chem Phys 2005; 122:064305. [PMID: 15740370 DOI: 10.1063/1.1846691] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A time-dependent wave packet method has been used to study different competing products of H(2)+H(2) collisions: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. A three-degree-of-freedom reduced dimensionality model has been used for five different geometries of the colliding complex (parallel H, crossed X, collinear L, and two T-shaped geometries T(I) and T(II)), with reactants in selected vibrational states with one diatom vibrationally "hot" and the other one vibrationally "cold." Product probabilities have been calculated using two potential energy surfaces [J. Chem. Phys. 101, 4004 (1994); J. Chem. Phys. 116, 666 (2002)] in order to compare their performance in the dynamics. The regions of the potential energy surfaces responsible of the threshold behavior of the probabilities have been identified. Overall, we have found that the most recent potential energy surface is less anisotropic, provides a smaller propensity for insertion-type processes, and gives lower energy thresholds.
Collapse
|
9
|
Hernández * MI, Campos-Martínez J, Van Caillie C, Di Domenico † D. Wave packet calculations for H2+H2collisions: isotopic substitution effects. Mol Phys 2004. [DOI: 10.1080/00268970412331294793] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
10
|
Mandy ME, Pogrebnya SK. Inelastic collisions of molecular hydrogen: A comparison of results from quantum and classical mechanics. J Chem Phys 2004; 120:5585-91. [PMID: 15267434 DOI: 10.1063/1.1649722] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Full-dimensional quantum and classical calculations have been carried out for inelastic (nonreactive) energy transfer in H2+H2 on the ab initio potential energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)]. State-to-state cross sections are determined and compared for transitions from H2(0,j(ab))+H2(1,j(cd)). While there is excellent agreement for transitions involving small Deltaj, for larger Deltaj and for vibrational relaxation, significant differences are observed which exhibit no systematic trends. Reasons for this disagreement are discussed.
Collapse
Affiliation(s)
- M E Mandy
- Program in Chemistry, University of Northern British Columbia, Prince George, British Columbia, Canada V2N 4Z9
| | | |
Collapse
|
11
|
Albertí M, Lucas JM, Aguilar A. A Quasiclassical Trajectory Study of the Multichannel H(1) + H(2)BO → BO + H2, H(1)BO + H(2), HOB + H Reaction. J Phys Chem A 2003. [DOI: 10.1021/jp026173n] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- M. Albertí
- Departament de Química Física i Centre de Recerca en Química Teòrica, Universitat de Barcelona, C/ Martí i Franquès, 1. 08028 Barcelona, Spain
| | - J. M. Lucas
- Departament de Química Física i Centre de Recerca en Química Teòrica, Universitat de Barcelona, C/ Martí i Franquès, 1. 08028 Barcelona, Spain
| | - A. Aguilar
- Departament de Química Física i Centre de Recerca en Química Teòrica, Universitat de Barcelona, C/ Martí i Franquès, 1. 08028 Barcelona, Spain
| |
Collapse
|
12
|
Ying Lin S, Guo H. Full-dimensional quantum wave packet study of rotationally inelastic transitions in H2+H2 collision. J Chem Phys 2002. [DOI: 10.1063/1.1500731] [Citation(s) in RCA: 90] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Ceballos A, Garcia E, Laganà A. Reaction and dissociation mechanism control: the H2 + H2system. Phys Chem Chem Phys 2002. [DOI: 10.1039/b205377b] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
14
|
di Domenico D, Hernández MI, Campos Martı́nez J. Wave packet calculations for H2(v1=10–14)+H2(v2=0–2): Reaction and dissociation mechanisms. J Chem Phys 2001. [DOI: 10.1063/1.1409357] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
15
|
Progress in validating the potential energy surface of the OH+H2 reaction: product vibrational distributions. Chem Phys Lett 2001. [DOI: 10.1016/s0009-2614(01)00917-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
16
|
|