1
|
Naskar K, Mukherjee S, Ghosh S, Adhikari S. Coupled 3D ( J ≥ 0) Time-Dependent Wave Packet Calculation for the F + H 2 Reaction on Accurate Ab Initio Multi-State Diabatic Potential Energy Surfaces. J Phys Chem A 2024; 128:1438-1456. [PMID: 38359800 DOI: 10.1021/acs.jpca.3c05590] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
We had calculated adiabatic potential energy surfaces (PESs), nonadiabatic, and spin-orbit (SO) coupling terms among the lowest three electronic states (12A', 22A', and 12A″) of the F + H2 system using the multireference configuration interaction (MRCI) level of theory, and the adiabatic-to-diabatic transformation equations were solved to formulate the diabatic Hamiltonian matrix [J. Chem. Phys. 2020, 153, 174301] for the entire region of the nuclear configuration space. The accuracy of such diabatic PESs is explored by performing scattering calculations to evaluate integral cross sections (ICSs) and rate constants. The nonadiabatic and SO effects are studied by utilizing coupled 3D time-dependent wave packet formalism with zero and nonzero total angular momentum on multiple adiabatic/diabatic surfaces calculation. We depict the convergence profiles of reaction probabilities for the reactive as well as nonreactive processes on various electronic states at different collision energies with respect to total angular momentum including all helicity quantum numbers. Finally, total ICSs are calculated as functions of collision energies for the initial rovibrational state (v = 0, j = 0) of the H2 molecule along with the temperature-dependent rate coefficient, where those quantities are compared with previous theoretical and experimental results.
Collapse
Affiliation(s)
- Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India
| | - Soumya Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India
| | - Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India
- Department of Chemical Sciences, Indian Institute of Science Education and Research, Kolkata 741246, West Bengal, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India
| |
Collapse
|
2
|
Naskar K, Ghosh S, Adhikari S. Accurate Calculation of Rate Constant and Isotope Effect for the F + H 2 Reaction by the Coupled 3D Time-Dependent Wave Packet Method on the Newly Constructed Ab Initio Ground Potential Energy Surface. J Phys Chem A 2022; 126:3311-3328. [PMID: 35594416 DOI: 10.1021/acs.jpca.2c01209] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We employ coupled three-dimensional (3D) time dependent wave packet formalism in hyperspherical coordinates for reactive scattering problem on the newly constructed ab initio calculated ground adiabatic potential energy surface for the F + H2/D2 reaction. The convergence profiles for various reactive channels are depicted at low collision energy regimes with respect to the total angular momentum (J) quantum numbers. For two different reactant diatomic molecules (H2 and D2) initially at their respective ground roto-vibrational state (v = 0, j = 0), calculated state-to-state as well as total integral cross sections as a function of collision energy, temperature dependent rate constants, and the kinetic isotope effect for various reactivity profiles of F + H2 and F + D2 reactions are presented along with previous theoretical and experimental results.
Collapse
Affiliation(s)
- Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India
| | - Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India.,Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata West Bengal-741246, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India
| |
Collapse
|
3
|
Ghosh S, Sharma R, Adhikari S, Varandas AJC. Dynamical calculations of O( 3P) + OH( 2Π) reaction on the CHIPR potential energy surface using the fully coupled time-dependent wave-packet approach in hyperspherical coordinates. Phys Chem Chem Phys 2021; 23:21784-21796. [PMID: 34550126 DOI: 10.1039/d1cp02488d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have carried out quantum dynamics calculations for the O + OH → H + O2 reaction on the CHIPR [A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117] potential energy surface (PES) for ground state HO2 using the fully coupled 3D time-dependent wavepacket formalism [S. Adhikari and A. J. C. Varandas, Comput. Phys. Commun., 2013, 184, 270] in hyperspherical coordinates. Reaction probabilities for J > 0 are calculated for different initial rotational states of the OH radical (v = 0; j = 0, 1). State-to-state as well as total integral cross sections and rate-coefficients are evaluated and compared with previous theoretical calculations and available experimental studies. Using the rate constant for the forward (hereinafter considered to be H + O2 → O + OH) and backward (O + OH → H + O2) reactions of this reactive system, the equilibrium constant of the reversible process [H + O2 ⇌ O + OH] is calculated as a function of temperature and compared with previous experimental measurements.
Collapse
Affiliation(s)
- Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India.
| | - Rahul Sharma
- Department of Chemistry, St. Xaviers' College, Kolkata-700016, West Bengal, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India.
| | - António J C Varandas
- School of Physics and Physical Engineering, Qufu Normal University, 273165 Qufu, China.,Department of Physics, Universidade Federal do Espírito Santo, 29075-910 Victória, Brazil.,Departamento de Qumica, and Centro de Qumica, Universidade de Coimbra, 3004-535 Coimbra, Portugal.
| |
Collapse
|
4
|
Pham TV, Lin M. Ab initio quantum-chemical and kinetic studies of the O(1D) + N2(X1Σg+) spin-forbidden quenching process. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138955] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
5
|
Ghosh S, Sahoo T, Baer M, Adhikari S. Charge Transfer Processes for H + H 2+ Reaction Employing Coupled 3D Wavepacket Approach on Beyond Born-Oppenheimer Based Ab Initio Constructed Diabatic Potential Energy Surfaces. J Phys Chem A 2021; 125:731-745. [PMID: 33461293 DOI: 10.1021/acs.jpca.0c08975] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The dynamics of the H + H2+ reaction has been analyzed from the electronically first excited state of diabatic potential energy surfaces constructed by employing the Beyond Born-Oppenheimer theory [J. Chem. Phys. 2014, 141, 204306]. We have employed the coupled 3D time-dependent wavepacket formalism in hyperspherical coordinates for multisurface reactive scattering problems. To be specific, the charge transfer processes have been investigated extensively by calculating state-to-state as well as total reaction probabilities and integral cross sections, when the reaction process is initiated from the first excited electronic state (21A'). We have depicted the convergence profiles of reaction probabilities for the competing charge transfer processes, namely, reactive charge transfer (RCT) and nonreactive charge transfer (NRCT) processes for different total energies with respect to total angular momentum, J. Total and state-to-state integral cross sections are calculated as a function of total energy for the initial rovibrational state, namely, v = 0, j = 0 level of H2+ (2Σg+) molecule and are compared with previous theoretical calculations. Finally, we have calculated temperature-dependent rate constants using our presently evaluated cross sections and compared their average with the experimentally measured one.
Collapse
Affiliation(s)
- Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata, 700 032, India
| | - Tapas Sahoo
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Michael Baer
- The Fritz Haber Center for Molecular Dynamics, The Hebrew University of Jerusalem, Jerusalem, 91904, Israel
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata, 700 032, India
| |
Collapse
|
6
|
Sathyamurthy N, Mahapatra S. Time-dependent quantum mechanical wave packet dynamics. Phys Chem Chem Phys 2020; 23:7586-7614. [PMID: 33306771 DOI: 10.1039/d0cp03929b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Starting from a model study of the collinear (H, H2) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH4. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H2) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.
Collapse
|
7
|
Ghosh S, Sharma R, Adhikari S, Varandas AJC. Fully coupled (J > 0) time-dependent wave-packet calculations using hyperspherical coordinates for the H + O2 reaction on the CHIPR potential energy surface. Phys Chem Chem Phys 2019; 21:20166-20176. [DOI: 10.1039/c9cp03171e] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
ICS calculation by time dependent wavepacket approach for H + O2 reaction using non-zero J values.
Collapse
Affiliation(s)
- Sandip Ghosh
- School of Chemical Sciences
- Indian Association for the Cultivation of Science
- Kolkata-700032
- India
| | - Rahul Sharma
- Department of Chemistry
- St. Xaviers' College
- Kolkata-700016
- India
| | - Satrajit Adhikari
- School of Chemical Sciences
- Indian Association for the Cultivation of Science
- Kolkata-700032
- India
| | - António J. C. Varandas
- School of Physics and Physical Engineering
- Qufu Normal University
- 273165 Qufu
- China
- Departamento de Química, and Centro de Química
| |
Collapse
|
8
|
An F, Han S, Hu X, Xie D, Guo H. First-principles dynamics of collisional intersystem crossing: resonance enhanced quenching of C(1D) by N2. Phys Chem Chem Phys 2019; 21:8645-8653. [DOI: 10.1039/c8cp07171c] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Intersystem crossing is a common and important nonadiabatic process in molecular systems. Revealed with a first-principles investigation, the quenching of C(1D) by N2 is efficient due to multiple passages via long-lived collisional resonances, despite relatively small spin–orbit couplings.
Collapse
Affiliation(s)
- Feng An
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University
- Nanjing 210023
- China
| | - Shanyu Han
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University
- Nanjing 210023
- China
- Department of Chemistry and Chemical Biology, University of New Mexico
- Albuquerque
| | - Xixi Hu
- 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
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico
- Albuquerque
- USA
| |
Collapse
|
9
|
Ghosh S, Sharma R, Adhikari S, Varandas AJC. 3D time-dependent wave-packet approach in hyperspherical coordinates for the H + O2 reaction on the CHIPR and DMBE IV potential energy surfaces. Phys Chem Chem Phys 2018; 20:478-488. [DOI: 10.1039/c7cp06254k] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
3D wavepacket quantum dynamics methodology ICS calculation of H + O2 reaction on the CHIPR and DMBE IV PESs by J-shifting scheme.
Collapse
Affiliation(s)
- Sandip Ghosh
- Department of Physical Chemistry
- Indian Association for the Cultivation of Science
- Kolkata-700032
- India
| | - Rahul Sharma
- Department of Chemistry
- St. Xaviers’ College
- Kolkata-700016
- India
| | - Satrajit Adhikari
- Department of Physical Chemistry
- Indian Association for the Cultivation of Science
- Kolkata-700032
- India
| | - António J. C. Varandas
- Departamento de Química
- and Centro de Química
- Universidade de Coimbra
- 3004-535 Coimbra
- Portugal
| |
Collapse
|
10
|
Grondin R, Loison JC, Hickson KM. Low Temperature Rate Constants for the Reactions of O(1D) with N2, O2, and Ar. J Phys Chem A 2016; 120:4838-44. [DOI: 10.1021/acs.jpca.5b12358] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Romain Grondin
- Université de Bordeaux and ‡CNRS, Institut des Sciences Moléculaires, UMR 5255, F-33400 Talence, France
| | - Jean-Christophe Loison
- Université de Bordeaux and ‡CNRS, Institut des Sciences Moléculaires, UMR 5255, F-33400 Talence, France
| | - Kevin M. Hickson
- Université de Bordeaux and ‡CNRS, Institut des Sciences Moléculaires, UMR 5255, F-33400 Talence, France
| |
Collapse
|
11
|
Defazio P, Gamallo P, Petrongolo C. Nonadiabatic dynamics of O(1D) + N2(XΣg+1)→O(3P) + N2(XΣg+1) on three coupled potential surfaces: Symmetry, Coriolis, spin-orbit, and Renner-Teller effects. J Chem Phys 2012; 136:054308. [DOI: 10.1063/1.3682467] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Brouard M, Quadrini F, Vallance C. The photodissociation dynamics of OCS at 248nm: The S(PJ3) atomic angular momentum polarization. J Chem Phys 2007; 127:084305. [PMID: 17764246 DOI: 10.1063/1.2757619] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dissociation of OCS has been investigated subsequent to excitation at 248 nm using velocity map ion imaging. Speed distributions, speed dependent translational anisotropy parameters, and the atomic angular momentum orientation and alignment are reported for the channel leading to S((3)P(J)). The speed distributions and beta parameters are in broad agreement with previous work and show behavior that is highly sensitive to the S-atom spin-orbit state. The data are shown to be consistent with the operation of at least two triplet production mechanisms. Interpretation of the angular momentum polarization data in terms of an adiabatic picture has been used to help identify a likely dissociation pathway for the majority of the S((3)P(J)) products, which strongly favors production of J=2 fragment atoms, correlated, it is proposed, with rotationally hot and vibrationally cold CO cofragments. For these fragments, optical excitation to the 2 (1)A(') surface is thought to constitute the first step, as for the singlet dissociation channel. This is followed by crossing, via a conical intersection, to the ground 1 (1)A(') state, from where intersystem crossing occurs, populating the 1 (3)A(')1 (3)A(")((3)Pi) states. The proposed mechanism provides a qualitative rationale for the observed spin-orbit populations, as well as the S((3)P(J)) quantum yield and angular momentum polarization. At least one other production mechanism, leading to a more statistical S-atom spin-orbit state distribution and rotationally cold, vibrationally hot CO cofragments, is thought to involve direct excitation to either the (3)Sigma(-) or (3)Pi states.
Collapse
Affiliation(s)
- M Brouard
- The Physical and Theoretical Chemistry Laboratory, The Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom.
| | | | | |
Collapse
|
13
|
Quan WL, Song Q, Tang BY. The effects of the rotational excitation on the Br + H2 reaction and its dependence on the potential energy surfaces. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.05.105] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
14
|
Lu RF, Zhang PY, Chu TS, Xie TX, Han KL. Spin-orbit effect in the energy pooling reaction O2(aΔ1)+O2(aΔ1)→O2(bΣ1)+O2(XΣ3). J Chem Phys 2007; 126:124304. [PMID: 17411120 DOI: 10.1063/1.2713399] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Five-dimensional nonadiabatic quantum dynamics studies have been carried out on two new potential energy surfaces of S(2)((1)A(')) and T(7)((3)A(")) states for the title oxygen molecules collision with coplanar configurations, along with the spin-orbit coupling between them. The ab initio calculations are based on complete active state second-order perturbation theory with the 6-31+G(d) basis set. The calculated spin-orbit induced transition probability as a function of collision energy is found to be very small for this energy pooling reaction. The rate constant obtained from a uniform J-shifting approach is compared with the existing theoretical and experimental data, and the spin-orbit effect is also discussed in this electronic energy-transfer process.
Collapse
Affiliation(s)
- Rui-Feng Lu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | | | | | | | | |
Collapse
|
15
|
Vissers GWM, McCoy AB. Time-Dependent Wave Packet Studies on the Cl + HCl Hydrogen Exchange Reaction. J Phys Chem A 2006; 110:5978-81. [PMID: 16671664 DOI: 10.1021/jp061196d] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The initiation of the hydrogen exchange reaction Cl((2)P)+HCl --> ClH+Cl((2)P) by excitation of the HCl molecular stretch to v=2 is studied for total angular momentum quantum number J=(1)/(2) and both even and odd parity. The calculations were performed using a time-dependent propagation from an initial quasi-bound state and employed all three relevant potential energy surfaces and the nonadiabatic couplings between them. Coriolis and spin-orbit coupling were also taken into account. The electronic and HCl rotational distributions of the products in both dissociation channels are analyzed, and the results are interpreted using features of the potential energy surfaces.
Collapse
Affiliation(s)
- Gé W M Vissers
- Department of Chemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | |
Collapse
|
16
|
Chu TS, Zhang Y, Han KL. The time-dependent quantum wave packet approach to the electronically nonadiabatic processes in chemical reactions. INT REV PHYS CHEM 2006. [DOI: 10.1080/01442350600677929] [Citation(s) in RCA: 416] [Impact Index Per Article: 23.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|