1
|
Montes de Oca-Estévez MJ, Darna B, García-Ruiz B, Prosmiti R, González-Lezana T, Koner D. Ar+ ArH + Reactive Collisions of Astrophysical Interest: The Case of 36 Ar. Chemphyschem 2023; 24:e202300450. [PMID: 37477047 DOI: 10.1002/cphc.202300450] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Revised: 07/19/2023] [Accepted: 07/20/2023] [Indexed: 07/22/2023]
Abstract
The reactive collision between 36 Ar and the 36 ArH+ species has been investigated by means of quantum mechanical (QM), quasiclassical trajectories (QCT) and statistical quantum mechanical (SQM) approaches. Reaction probabilities, cross sections as a function of the energy and rate constants in terms of the temperature have been obtained. Cumulative distributions as a function of the collision time and the inspection of selected QCT corresponding to specific dynamical mechanisms have been analysed. Predictions by means of the SQM method are in good agreement with the QM results, thus supporting the complex-forming nature of the process.
Collapse
Affiliation(s)
- María Judit Montes de Oca-Estévez
- Instituto de Física Fundamental, IFF-CSIC, Serrano 123, 28006, Madrid, Spain
- Atelgraphics S.L., Mota de Cuervo 42, 28043, Madrid, Spain
- Doctoral Programme in Theoretical Chemistry and Computational Modelling, Doctoral School, UAM, Madrid, Spain
| | - Beatriz Darna
- Faculty of Science, KU Leuven, Celestijnenlaan 200H, box 2100, 3001 Heverlee, Leuven, Belgium
| | - Borja García-Ruiz
- Instituto de Física Fundamental, IFF-CSIC, Serrano 123, 28006, Madrid, Spain
| | - Rita Prosmiti
- Instituto de Física Fundamental, IFF-CSIC, Serrano 123, 28006, Madrid, Spain
| | | | - Debasish Koner
- Department of Chemistry, Indian Institute of Technology Hyderabad Kandi, Sangareddy, 502284, Telangana, India
| |
Collapse
|
2
|
Goswami S, San Vicente Veliz JC, Upadhyay M, Bemish RJ, Meuwly M. Quantum and quasi-classical dynamics of the C( 3P) + O 2( 3Σ-g) → CO( 1Σ +) + O( 1D) reaction on its electronic ground state. Phys Chem Chem Phys 2022; 24:23309-23322. [PMID: 36165004 PMCID: PMC9533374 DOI: 10.1039/d2cp02840a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dynamics of the C(3P) + O2(3Σ−g) → CO(1Σ+) + O(1D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations. For the moderate collision energies considered (Ec = 0.001 to 0.4 eV, corresponding to a range from 10 K to 4600 K) the total reaction probabilities from the two different treatments of the nuclear dynamics agree very favourably. The undulations present in P(E) from the quantum mechanical treatment can be related to stabilization of the intermediate CO2 complex with lifetimes on the 0.05 ps time scale. This is also confirmed from direct analysis of the TDWP simulations and QCT trajectories. Product diatom vibrational and rotational level resolved state-to-state reaction probabilities from TDWP and QCT simulations agree well except for the highest product vibrational states (v′ ≥ 15) and for the lowest product rotational states (j′ ≤ 10). Opening of the product vibrational level CO(v′ = 17) requires ∼0.2 eV from QCT and TDWP simulations with O2(j = 0) and decreases to 0.04 eV if all initial rotational states are included in the QCT analysis, compared with Ec > 0.04 eV obtained from experiments. It is thus concluded that QCT simulations are suitable for investigating and realistically describe the C(3P) + O2(3Σ−g) → CO(1Σ+) + O(1D) reaction down to low collision energies when compared with results from a quantum mechanical treatment using TDWPs. The dynamics of the C(3P) + O2(3Σ−g) → CO(1Σ+) + O(1D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations.![]()
Collapse
Affiliation(s)
- Sugata Goswami
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | | | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Raymond J Bemish
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland AFB, New Mexico 87117, USA
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland. .,Department of Chemistry, Brown University, RI, USA
| |
Collapse
|
3
|
Patra S, San Vicente Veliz JC, Koner D, Bieske EJ, Meuwly M. Photodissociation dynamics of N3+. J Chem Phys 2022; 156:124307. [DOI: 10.1063/5.0085081] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The photodissociation dynamics of [Formula: see text] excited from its (linear) 3[Formula: see text]/(bent) 3A″ ground to the first excited singlet and triplet states is investigated. Three-dimensional potential energy surfaces for the 1A′, 1A″, and 3A′ electronic states, correlating with the 1Δg and 3Πu states in linear geometry, for [Formula: see text] are constructed using high-level electronic structure calculations and represented as reproducing kernels. The reference ab initio energies are calculated at the MRCI+Q/aug-cc-pVTZ level of theory. For following the photodissociation dynamics in the excited states, rotational and vibrational distributions P( v′) and P( j′) for the N2 product are determined from vertically excited ground state distributions. Due to the different shapes of the ground state 3A″ potential energy surface and the excited states, appreciable angular momentum j′ ∼ 60 is generated in diatomic fragments. The lifetimes in the excited states extend to at least 50 ps. Notably, results from sampling initial conditions from a thermal ensemble and from the Wigner distribution of the ground state wavefunction are comparable.
Collapse
Affiliation(s)
- Sarbani Patra
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | | | - Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Mangalam, Tirupati 517507, Andhra Pradesh, India
| | - Evan J. Bieske
- Department of Chemistry, University of Melbourne, Parkville 3010, Australia
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
4
|
Arnold J, San Vicente Veliz JC, Koner D, Singh N, Bemish RJ, Meuwly M. Machine learning product state distributions from initial reactant states for a reactive atom–diatom collision system. J Chem Phys 2022; 156:034301. [DOI: 10.1063/5.0078008] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Julian Arnold
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland
| | | | - Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Mangalam, Tirupati 517507, Andhra Pradesh, India
| | - Narendra Singh
- Department of Mechanical Engineering, Stanford University, Stanford, California 94305, USA
| | - Raymond J. Bemish
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland AFB, New Mexico 87117, USA
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
5
|
Koner D, Barrios L, González-Lezana T, Panda AN. Atom-Diatom Reactive Scattering Collisions in Protonated Rare Gas Systems. Molecules 2021; 26:4206. [PMID: 34299481 PMCID: PMC8304066 DOI: 10.3390/molecules26144206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2021] [Revised: 07/07/2021] [Accepted: 07/08/2021] [Indexed: 11/17/2022] Open
Abstract
The study of the dynamics of atom-diatom reactions involving two rare gas (Rg) atoms and protons is of crucial importance given the astrophysical relevance of these processes. In a series of previous studies, we have been investigating a number of such Rg(1)+ Rg(2)H+→ Rg(2)+ Rg(1)H+ reactions by means of different numerical approaches. These investigations comprised the construction of accurate potential energy surfaces by means of ab initio calculations. In this work, we review the state-of-art of the study of these protonated Rg systems making special emphasis on the most relevant features regarding the dynamical mechanisms which govern these reactive collisions. The aim of this work therefore is to provide an as complete as possible description of the existing information regarding these processes.
Collapse
Affiliation(s)
- Debasish Koner
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Tirupati 517507, Andhra Pradesh, India;
| | - Lizandra Barrios
- Department of Chemistry, CMS—Centre for Molecular Simulation, IQST—Institute for Quantum Science and Technology and Quantum Alberta, University of Calgary, 2500 University Drive N.W., Calgary, AB T2N 1N4, Canada;
| | | | - Aditya N. Panda
- Department of Chemistry, Indian Institute of Technology Guwahati, Assam 781039, India;
| |
Collapse
|
6
|
San Vicente Veliz JC, Koner D, Schwilk M, Bemish RJ, Meuwly M. The C( 3P) + O 2( 3Σ g-) → CO 2 ↔ CO( 1Σ +) + O( 1D)/O( 3P) reaction: thermal and vibrational relaxation rates from 15 K to 20 000 K. Phys Chem Chem Phys 2021; 23:11251-11263. [PMID: 33949507 PMCID: PMC8133592 DOI: 10.1039/d1cp01101d] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Thermal rates for the C(3P) + O2(3Σg−) ↔ CO(1Σ+)+ O(1D)/O(3P) reaction are investigated over a wide temperature range based on quasi classical trajectory (QCT) simulations on 3-dimensional, reactive potential energy surfaces (PESs) for the 1A′, (2)1A′, 1A′′, 3A′ and 3A′′ states. These five states are the energetically low-lying states of CO2 and their PESs are computed at the MRCISD+Q/aug-cc-pVTZ level of theory using a state-average CASSCF reference wave function. Analysis of the different electronic states for the CO2 → CO + O dissociation channel rationalizes the topography of this region of the PESs. The forward rates from QCT simulations match measurements between 15 K and 295 K whereas the equilibrium constant determined from the forward and reverse rates is consistent with that derived from statistical mechanics at high temperature. Vibrational relaxation, O + CO(ν = 1,2) → O + CO(ν = 0), is found to involve both, non-reactive and reactive processes. The contact time required for vibrational relaxation to take place is τ ≥ 150 fs for non-reacting and τ ≥ 330 fs for reacting (oxygen atom exchange) trajectories and the two processes are shown to probe different parts of the global potential energy surface. In agreement with experiments, low collision energy reactions for the C(3P) + O2(3Σg−, ν = 0) → CO(1Σ+) + O(1D) lead to CO(1Σ+, ν′ = 17) with an onset at Ec ∼ 0.15 eV, dominated by the 1A′ surface with contributions from the 3A′ surface. Finally, the barrier for the COA(1Σ+) + OB(3P) → COB(1Σ+) + OA(3P) atom exchange reaction on the 3A′ PES yields a barrier of ∼7 kcal mol−1 (0.300 eV), consistent with an experimentally reported value of 6.9 kcal mol−1 (0.299 eV). Reaction and vibrational relaxation rate computed for C(3P) + O2(3Σg−) ↔ CO(1Σ+) + O(1D)/O(3P) for a wide range of temperatures using quasiclassical trajectory calculations on five new potential energy surfaces for different electronic states.![]()
Collapse
Affiliation(s)
| | - Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Max Schwilk
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland. and University of Vienna, Faculty of Physics, 1090 Vienna, Austria
| | - Raymond J Bemish
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland AFB, New Mexico 87117, USA
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland. and Brown University, Providence, RI 02912, USA
| |
Collapse
|
7
|
On the Proton-Bound Noble Gas Dimers (Ng-H-Ng) + and (Ng-H-Ng') + (Ng, Ng'= He-Xe): Relationships betweenStructure, Stability, and Bonding Character. Molecules 2021; 26:molecules26051305. [PMID: 33671081 PMCID: PMC7957648 DOI: 10.3390/molecules26051305] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2021] [Revised: 02/23/2021] [Accepted: 02/24/2021] [Indexed: 11/17/2022] Open
Abstract
The structure, stability, and bonding character of fifteen (Ng-H-Ng)+ and (Ng-H-Ng’)+ (Ng, Ng’ = He-Xe) compounds were explored by theoretical calculations performed at the coupled cluster level of theory. The nature of the stabilizing interactions was, in particular, assayed using a method recently proposed by the authors to classify the chemical bonds involving the noble-gas atoms. The bond distances and dissociation energies of the investigated ions fall in rather large intervals, and follow regular periodic trends, clearly referable to the difference between the proton affinity (PA) of the various Ng and Ng’. These variations are nicely correlated with the bonding situation of the (Ng-H-Ng)+ and (Ng-H-Ng’)+. The Ng-H and Ng’-H contacts range, in fact, between strong covalent bonds to weak, non-covalent interactions, and their regular variability clearly illustrates the peculiar capability of the noble gases to undergo interactions covering the entire spectrum of the chemical bond.
Collapse
|
8
|
Koner D. Quantum and quasiclassical dynamical simulations for the Ar 2H + on a new global analytical potential energy surface. J Chem Phys 2021; 154:054303. [PMID: 33557552 DOI: 10.1063/5.0039252] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
A new analytical potential energy surface (PES) has been constructed for the Ar2H+ system from a dataset consisting of a large number of ab initio energies computed using the coupled-cluster singles, doubles and perturbative triples method and aug-cc-pVQZ basis set. The long-range interaction is added to the diatomic potentials using a standard long range expansion form to better describe the asymptotic regions. The vibrational states for the most stable structures of the Ar2H+ system have been calculated, and few low lying states are assigned to quantum numbers. Reactive scattering studies have been performed for the Ar + Ar'H+ → Ar' + ArH+ proton exchange reaction on the newly generated PES. Reaction probability, cross sections, and rate constants are calculated for the Ar + Ar'H+(v = 0, j = 0) collisions within 0.01 eV-0.6 eV of relative translational energy using exact quantum dynamical simulations as well as quasiclassical trajectory (QCT) calculations. The effect of vibrational excitation of the reactants is also explored for the reaction. State averaged rate constants are calculated for the proton exchange reaction at different temperatures using the QCT method. The mechanistic pathways for the reaction are understood by analyzing the quasiclassical trajectories.
Collapse
Affiliation(s)
- Debasish Koner
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Mangalam, Tirupati 517507, Andhra Pradesh, India
| |
Collapse
|
9
|
Arnold J, Koner D, Käser S, Singh N, Bemish RJ, Meuwly M. Machine Learning for Observables: Reactant to Product State Distributions for Atom–Diatom Collisions. J Phys Chem A 2020; 124:7177-7190. [DOI: 10.1021/acs.jpca.0c05173] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Julian Arnold
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Silvan Käser
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Narendra Singh
- Department of Mechanical Engineering, Stanford University Stanford, California 94305, United States
| | - Raymond J. Bemish
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland AFB, New Mexico 87117, United States
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
10
|
Koner D, Bemish RJ, Meuwly M. Dynamics on Multiple Potential Energy Surfaces: Quantitative Studies of Elementary Processes Relevant to Hypersonics. J Phys Chem A 2020; 124:6255-6269. [DOI: 10.1021/acs.jpca.0c01870] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Raymond J. Bemish
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland AFB, New Mexico 87117, United States
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
11
|
San Vicente Veliz JC, Koner D, Schwilk M, Bemish RJ, Meuwly M. The N(4S) + O2(X3Σ−g) ↔ O(3P) + NO(X2Π) reaction: thermal and vibrational relaxation rates for the 2A′, 4A′ and 2A′′ states. Phys Chem Chem Phys 2020; 22:3927-3939. [DOI: 10.1039/c9cp06085e] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cross sections, rates, equilibrium constants and vibrational relaxation times for the N(4S) + O2(X3Σ−g) ↔ O(3P) + NO(X2Π) reaction from simulations on new, RKHS-based surfaces for the three lowest electronic states.
Collapse
Affiliation(s)
| | - Debasish Koner
- Department of Chemistry
- University of Basel
- CH-4056 Basel
- Switzerland
| | - Max Schwilk
- Department of Chemistry
- University of Basel
- CH-4056 Basel
- Switzerland
| | - Raymond J. Bemish
- Air Force Research Laboratory
- Space Vehicles Directorate
- Kirtland AFB
- USA
| | - Markus Meuwly
- Department of Chemistry
- University of Basel
- CH-4056 Basel
- Switzerland
| |
Collapse
|
12
|
González-Lezana T, Bossion D, Scribano Y, Bhowmick S, Suleimanov YV. Dynamics of H + HeH +( v = 0, j = 0) → H 2+ + He: Insight on the Possible Complex-Forming Behavior of the Reaction. J Phys Chem A 2019; 123:10480-10489. [PMID: 31725286 DOI: 10.1021/acs.jpca.9b06122] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The H + HeH+→ He + H2+ reaction has been studied by means of a combination of theoretical approaches: a statistical quantum method (SQM), ring polymer molecular dynamics (RPMD), and the quasiclassical trajectory (QCT) method. Cross sections and rate constants have been calculated in an attempt to investigate the dynamics of the process. The comparison with previous calculations and experimental results reveals that despite the fact that statistical predictions seem to reproduce some of the overall observed features, the analysis at a more detailed state-to-state level shows noticeable deviations from a complex-forming dynamics. We find some differences in cross sections and rate constants obtained in the QCT calculation with a Gaussian binning procedure with respect to previous works in which the standard histogram binning was employed.
Collapse
Affiliation(s)
| | - Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - Somnath Bhowmick
- Computation-based Science and Technology Research Center , The Cyprus Institute , 20 Konstantinou Kavafi Street , Nicosia 2121 , Cyprus
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center , The Cyprus Institute , 20 Konstantinou Kavafi Street , Nicosia 2121 , Cyprus
| |
Collapse
|
13
|
Simulation of experimental imaging results for four isotopic variants of the OH + CH4 reaction with a simple and relatively accurate theoretical approach. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.09.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
14
|
|