1
|
Bostan D, Mandal B, Joy C, Żółtowski M, Lique F, Loreau J, Quintas-Sánchez E, Batista-Planas A, Dawes R, Babikov D. Mixed quantum/classical calculations of rotationally inelastic scattering in the CO + CO system: a comparison with fully quantum results. Phys Chem Chem Phys 2024; 26:6627-6637. [PMID: 38115799 DOI: 10.1039/d3cp05369e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
An updated version of the CO + CO potential energy surface from [R. Dawes, X. G. Wang and T. Carrington, J. Phys. Chem. A 2013, 117, 7612] is presented, that incorporates an improved treatment of the asymptotic behavior. It is found that this new surface is only slightly different from the other popular PES available for this system in the literature [G. W. M. Vissers, P. E. S. Wormer and A. Van Der Avoird, Phys. Chem. Chem. Phys. 2003, 5, 4767]. The differences are quantified by expanding both surfaces over a set of analytic functions and comparing the behavior of expansion coefficients along the molecule-molecule distance R. It is shown that all expansion coefficients behave similarly, except in the very high energy range at small R where the PES is repulsive. That difference has no effect on low collision-energy dynamics, which is explored via inelastic scattering calculations carried out using the MQCT program which implements the mixed quantum/classical theory for molecular energy exchange processes. The validity of MQCT predictions of state-to-state transition cross sections for CO + CO is also tested by comparison against full-quantum coupled-states calculations. In all cases MQCT gives reliable results, except at very low collision energy where the full-quantum calculations predict strong oscillations of state-to-state transition cross sections due to resonances. For strong transitions with large cross sections, the results of MQCT are reliable, especially at higher collision energy. For weaker transitions, and lower collision energies, the cross sections predicted by MQCT may be up to a factor of 2-3 different from those obtained by full-quantum calculations.
Collapse
Affiliation(s)
- Dulat Bostan
- Chemistry Department, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Bikramaditya Mandal
- Chemistry Department, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Carolin Joy
- Chemistry Department, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Michał Żółtowski
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes)-UMR 6251, F-35000 Rennes, France
| | - François Lique
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes)-UMR 6251, F-35000 Rennes, France
| | - Jérôme Loreau
- KU Leuven, Department of Chemistry, B-3001 Leuven, Belgium
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Adrian Batista-Planas
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Dmitri Babikov
- Chemistry Department, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| |
Collapse
|
2
|
Tajouo Tela H, Quintas-Sánchez E, Dubernet ML, Scribano Y, Dawes R, Gatti F, Ndengué S. Rovibrational states calculations of the H 2O-HCN heterodimer with the multiconfiguration time dependent Hartree method. Phys Chem Chem Phys 2023; 25:31813-31824. [PMID: 37966067 DOI: 10.1039/d3cp03225f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
Water and hydrogen cyanide are two of the most common species in space and the atmosphere with the ability of binding to form dimers such as H2O-HCN. In the literature, while calculations characterizing various properties of the H2O-HCN cluster (equilibrium distance, vibrational frequencies and rotational constants) have been done in the past, extensive calculations of the rovibrational states of this system using a reliable quantum dynamical approach have yet to be reported. In this work, we intend to mend that by performing the first calculation of the rovibrational states of the H2O-HCN van der Waals complex on a recently developed potential energy surface. We use the block improved relaxation procedure implemented in the Heidelberg MultiConfiguration Time-Dependent Hartree (MCTDH) package to compute the states of the H2O-HCN isomer, from which we extract the transition frequencies and rotational constants of the complex. We further adapt an approach first suggested by Wang and Carrington-and supported here by analysis routines of the Heidelberg MCTDH package-to properly characterize the computed rovibrational states. The subsequent assignment of rovibrational states was done by theoretical analysis and visual inspection of the wavefunctions. Our simulations provide a Zero Point Energy (ZPE) and intermolecular vibrational frequencies in good agreement with past ab initio calculations. The transition frequencies and rotational constants obtained from our simulations match well with the available experimental data. This work has the broad aim to propose the MCTDH approach as a reliable option to compute and characterize rovibrational states of van der Waals complexes such as the current one.
Collapse
Affiliation(s)
- Hervé Tajouo Tela
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali, Rwanda.
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, 65409 Rolla, Missouri, USA
| | - Marie-Lise Dubernet
- LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne University, UPMC Univ Paris 06, 75014 Paris, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, 65409 Rolla, Missouri, USA
| | - Fabien Gatti
- Institut de Sciences Moleculaires d'Orsay, UMR 8214, Université Paris-Sud - Université Paris-Saclay, 91405 Orsay, France
| | - Steve Ndengué
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali, Rwanda.
- Department of Physics, Trinity College, 06106 Hartford, Connecticut, USA
| |
Collapse
|
3
|
Olejnik A, Jóźwiak H, Gancewski M, Quintas-Sánchez E, Dawes R, Wcisło P. Ab initio quantum scattering calculations and a new potential energy surface for the HCl(X1Σ+)-O2(X3Σg-) system: Collision-induced line shape parameters for O2-perturbed R(0) 0-0 line in H35Cl. J Chem Phys 2023; 159:134301. [PMID: 37782252 DOI: 10.1063/5.0169968] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2023] [Accepted: 09/01/2023] [Indexed: 10/03/2023] Open
Abstract
The remote sensing of abundance and properties of HCl-the main atmospheric reservoir of Cl atoms that directly participate in ozone depletion-is important for monitoring the partitioning of chlorine between "ozone-depleting" and "reservoir" species. Such remote studies require knowledge of the shapes of molecular resonances of HCl, which are perturbed by collisions with the molecules of the surrounding air. In this work, we report the first fully quantum calculations of collisional perturbations of the shape of a pure rotational line in H35Cl perturbed by an air-relevant molecule [as the first model system we choose the R(0) line in HCl perturbed by O2]. The calculations are performed on our new highly accurate HCl(X1Σ+)-O2(X3Σg-) potential energy surface. In addition to pressure broadening and shift, we also determine their speed dependencies and the complex Dicke parameter. This gives important input to the community discussion on the physical meaning of the complex Dicke parameter and its relevance for atmospheric spectra (previously, the complex Dicke parameter for such systems was mainly determined from phenomenological fits to experimental spectra and the physical meaning of its value in that context is questionable). We also calculate the temperature dependence of the line shape parameters and obtain agreement with the available experimental data. We estimate the total combined uncertainties of our calculations at 2% relative root-mean-square error in the simulated line shape at 296 K. This result constitutes an important step toward computational population of spectroscopic databases with accurate ab initio line shape parameters for molecular systems of terrestrial atmospheric importance.
Collapse
Affiliation(s)
- Artur Olejnik
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| | - Hubert Jóźwiak
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| | - Maciej Gancewski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409-0010, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409-0010, USA
| | - Piotr Wcisło
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| |
Collapse
|
4
|
Ndengué S, Quintas-Sánchez E, Dawes R, Blackstone CC, Osborn DL. Temperature Dependence of the Electronic Absorption Spectrum of NO 2. J Phys Chem A 2023. [PMID: 37384555 DOI: 10.1021/acs.jpca.3c02832] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/01/2023]
Abstract
The nitrogen dioxide (NO2) radical is composed of the two most abundant elements in the atmosphere, where it can be formed in a variety of ways including combustion, detonation of energetic materials, and lightning. Relevant also to smog and ozone cycles, together these processes span a wide range of temperatures. Remarkably, high-resolution NO2 electronic absorption spectra have only been reported in a narrow range below about 300 K. Previously, we reported [ J. Phys. Chem. A 2021, 125, 5519-5533] the construction of quasi-diabatic potential energy surfaces (PESs) for the lowest four electronic states (X̃, Ã, B̃, and C̃) of NO2. In addition to three-dimensional PESs based on explicitly correlated MRCI(Q)-F12/VTZ-F12 ab initio data, the geometry dependence of each component of the dipoles and transition dipoles was also mapped into fitted surfaces. The multiconfigurational time-dependent Hartree (MCTDH) method was then used to compute the 0 K electronic absorption spectrum (from the ground rovibrational initial state) employing those energy and transition dipole surfaces. Here, in an extension of that work, we report an investigation into the effects of elevated temperature on the spectrum, considering the effects of the population of rotationally and vibrationally excited initial states. The calculations are complemented by new experimental measurements. Spectral contributions from hundreds of rotational states up to N = 20 and from 200 individually-characterized vibrational states were computed. A spectral simulation tool was developed that enables modeling the spectrum at various temperatures─by weighting individual spectral contributions via the partition function, or for pure excited initial states, which can be probed via transient absorption spectroscopy. We validate these results against experimental absorption spectroscopy data at high temperatures, as well as via a new measurement from the (1,0,1) initial vibrational state.
Collapse
Affiliation(s)
- Steve Ndengué
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali 4285, Rwanda
| | | | - Richard Dawes
- Missouri University of Science and Technology, Rolla, Missouri 65409-0010, United States
| | - Christopher C Blackstone
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
| | - David L Osborn
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
- Department of Chemical Engineering, University of California, Davis, California 95616, United States
| |
Collapse
|
5
|
Dumouchel F, Quintas-Sánchez E, Balança C, Dawes R, Lique F, Feautrier N. Collisional excitation of C2H- by H2: New interaction potential and scattering calculations. J Chem Phys 2023; 158:2885323. [PMID: 37096848 DOI: 10.1063/5.0148119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 04/12/2023] [Indexed: 04/26/2023] Open
Abstract
Interstellar anions play an important role in astrochemistry as being tracers of the physical and chemical conditions in cold molecular clouds and circumstellar gas. The local thermodynamic equilibrium is generally not fulfilled in media where anions are detected and radiative and collisional data are required to model the observed lines. The C2H- anion has not yet been detected in the interstellar medium; however, collisional data could be used for non-LTE models that would help in identifying the most intense lines. For this purpose, we have computed the first 4D potential energy surface (PES) of the C2H--H2 complex using an explicitly correlated coupled-cluster approach. The PES is characterized by a single deep minimum with a well-depth of 924.96 cm-1. From this interaction potential, we derived excitation cross sections and rate coefficients of C2H- induced by collisions with para- and ortho-H2. The results obtained for collisions with para-H2 are compared to previous calculations performed using a 2D-PES obtained from an average over H2 rotations.
Collapse
Affiliation(s)
- Fabien Dumouchel
- LOMC - UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, 76 063 Le Havre cedex, France
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Christian Balança
- Observatoire de Paris, Université PSL, CNRS-UMR 8112, LERMA, F-92195 Meudon, France
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - François Lique
- University Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | - Nicole Feautrier
- Observatoire de Paris, Université PSL, CNRS-UMR 8112, LERMA, F-92195 Meudon, France
| |
Collapse
|
6
|
Zadrożny A, Jóźwiak H, Quintas-Sánchez E, Dawes R, Wcisło P. Publisher’s Note: “ Ab initio quantum scattering calculations for the CO–O 2 system and a new CO–O 2 potential energy surface: O 2 and air broadening of the R(0) line in CO” [J. Chem. Phys. 157, 174310 (2022)]. J Chem Phys 2022; 157:229901. [DOI: 10.1063/5.0135156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Adam Zadrożny
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| | - Hubert Jóźwiak
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409-0010, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409-0010, USA
| | - Piotr Wcisło
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| |
Collapse
|
7
|
Zadrożny A, Jóźwiak H, Quintas-Sánchez E, Dawes R, Wcisło PJ. Ab initio quantum scattering calculations for the CO-O 2 system and a new CO-O 2 potential energy surface: O 2 and air broadening of the R(0) line in CO. J Chem Phys 2022; 157:174310. [DOI: 10.1063/5.0115654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present ab initio calculations of the collisional broadening of the R(0) pure rotational line in CO (at 115 GHz) perturbed by O2. Our calculations are done in a fully quantum way by solving close-coupling quantum-scattering equations without any approximations. We also report a new highly accurate CO-O2 potential energy surface on which we did the quantum-scattering calculations. The calculated collisional broadening agrees with the available experimental data in a wide temperature range. The calculated collisional shift is negligible compared to the broadening, which is also consistent with the experimental data. We combine this result with our previous calculations for the same line in CO perturbed by N2 [J. Chem. Phys. 154, 054314 (2021)]; the obtained air-perturbed broadening of the R(0) pure rotational line in CO and its temperature dependence perfectly agree with the HITRAN database. This result constitutes an important step towards developing a methodology for providing accurate ab initio reference data on spectroscopic collisional line-shape parameters for molecular systems relevant for Earth atmosphere and for populating spectroscopic line-by-line databases.
Collapse
Affiliation(s)
- Adam Zadrożny
- Department of Atomic, Molecular and Optical Physics, Nicolaus Copernicus University in Torun Institute of Physics, Poland
| | - Hubert Jóźwiak
- Nicolaus Copernicus University Institute of Physics, Poland
| | | | - Richard Dawes
- Missouri University of Science and Technology, United States of America
| | - Piotr Jan Wcisło
- Nicolaus Copernicus University in Torun Institute of Physics, Poland
| |
Collapse
|
8
|
Ajili Y, Quintas-Sánchez E, Mehnen B, Żuchowski PS, Brzęk F, Kork NE, Gacesa M, Dawes R, Hochlaf M. Theoretical study of the CO2–O2 van der Waals complex: potential energy surface and applications. Phys Chem Chem Phys 2022; 24:28984-28993. [DOI: 10.1039/d2cp04101d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO2-O2) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single,...
Collapse
|
9
|
Quintas-Sánchez E, Dawes R, Denis-Alpizar O. Theoretical study of the HCS+–H2 van der Waals complex: potential energy surface, rovibrational bound states, and rotationally inelastic collisional cross sections. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1980234] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, USA
| | - Otoniel Denis-Alpizar
- Instituto de Ciencias Químicas Aplicadas, Facultad de Ingeniería, Universidad Autónoma de Chile, Santiago, Chile
| |
Collapse
|
10
|
Gancewski M, Jóźwiak H, Quintas-Sánchez E, Dawes R, Thibault F, Wcisło P. Fully quantum calculations of O 2-N 2 scattering using a new potential energy surface: Collisional perturbations of the oxygen 118 GHz fine structure line. J Chem Phys 2021; 155:124307. [PMID: 34598560 DOI: 10.1063/5.0063006] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A proper description of the collisional perturbation of the shapes of molecular resonances is important for remote spectroscopic studies of the terrestrial atmosphere. Of particular relevance are the collisions between the O2 and N2 molecules-the two most abundant atmospheric species. In this work, we report a new highly accurate O2(X3Σg -)-N2(X1Σg +) potential energy surface and use it for performing the first quantum scattering calculations addressing line shapes for this system. We use it to model the shape of the 118 GHz fine structure line in O2 perturbed by collisions with N2 molecules, a benchmark system for testing our methodology in the case of an active molecule in a spin triplet state. The calculated collisional broadening of the line agrees well with the available experimental data over a wide temperature range relevant for the terrestrial atmosphere. This work constitutes a step toward populating the spectroscopic databases with ab initio line shape parameters for atmospherically relevant systems.
Collapse
Affiliation(s)
- Maciej Gancewski
- Institute of Physics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| | - Hubert Jóźwiak
- Institute of Physics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409-0010, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409-0010, USA
| | - Franck Thibault
- Univ. Rennes, CNRS, IPR (Institut de Physique de Rennes)-UMR 6251, Rennes F-35000, France
| | - Piotr Wcisło
- Institute of Physics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
| |
Collapse
|
11
|
Ndengué S, Quintas-Sánchez E, Dawes R, Osborn D. The Low-Lying Electronic States of NO 2: Potential Energy and Dipole Surfaces, Bound States, and Electronic Absorption Spectrum. J Phys Chem A 2021; 125:5519-5533. [PMID: 34114826 DOI: 10.1021/acs.jpca.1c03482] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nitrogen dioxide, NO2, is a free radical composed of the two most abundant elements in Earth's atmosphere, nitrogen and oxygen, and is relevant to atmospheric and combustion chemistry. The electronic structure of even its lowest-lying states is remarkably complex, with various conical intersections and Renner-Teller pairings, giving rise to complex and perturbed vibronic states. Here we report some analysis of the 18 molecular states of doublet spin-multiplicity formed by combining ground-state N(4Su) and O(3Pg) atoms. Three-dimensional potential energy surfaces were fit at the MRCI(Q)-F12/VTZ-F12 level, describing the lowest four (X̃, Ã, B̃, and C̃) electronic states. A properties-based diabatization procedure was applied to accommodate the intersections, producing energies in a quasidiabatic representation and yielding couplings that were also fit into surfaces. The low-lying vibrational levels on the ground X̃ state were computed and compared with experimental measurements. Compared to experiment, the lowest 125 calculated vibrational levels (up to 8500 cm-1 above the zero-point energy) have a root-mean-squared error of 16.5 cm-1. In addition, dipole moments for each of the lowest four electronic states-and the transition dipoles between them-were also computed and fit. With the coupled energy and dipole surfaces, the electronic spectrum was calculated in absolute intensity and compared with experimental measurements. Detailed structure in the experimental spectrum was successfully reproduced, and the total integrated intensity matches experiment to an accuracy of ∼1.5% with no empirical adjustments.
Collapse
Affiliation(s)
- Steve Ndengué
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali, Rwanda
| | | | - Richard Dawes
- Missouri University of Science and Technology, Rolla, Missouri 65409-0010, United States
| | - David Osborn
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States.,Department of Chemical Engineering, University of California, Davis, California 95616, United States
| |
Collapse
|
12
|
Bop CT, Quintas-Sánchez E, Sur S, Robin M, Lique F, Dawes R. Inelastic scattering in isotopologues of O 2-Ar: the effects of mass, symmetry, and density of states. Phys Chem Chem Phys 2021; 23:5945-5955. [PMID: 33666616 DOI: 10.1039/d1cp00326g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The two species considered here, O2 (oxygen molecule) and Ar (argon-atom), are both abundant components of Earth's atmosphere and hence familiar collision partners in this medium. O2 is quite reactive and extensively involved in atmospheric chemistry, including Chapman's cycle of the formation and destruction of ozone; while Ar, like N2, typically plays the nevertheless crucial role of inert collider. Inert species can provide stabilization to metastable encounter-complexes through the energy transfer associated with inelastic collisions. The interplay of collision frequency and energy transfer efficiency, with state lifetimes and species concentrations, contributes to the rich and varied chemistry and dynamics found in diverse environments ranging from planetary atmospheres to the interstellar and circumstellar media. The nature and density of bound and resonance states, coupled electronic states, symmetry, and nuclear spin-statistics can all play a role. Here, we systematically investigate some of those factors by looking at the O2-Ar system, comparing rigorous quantum-scattering calculations for the 16O16O-40Ar, 18O16O-40Ar, and 18O18O-40Ar isotope combinations. A new accurate potential energy surface was constructed for this purpose holding the O2 bond distance at its vibrationally averaged distance.
Collapse
Affiliation(s)
- Cheikh T Bop
- Laboratoire Ondes et Milieux Complexes, UMR 6294, Centre National de la Recherche Scientifique-Universite du Havre, F-76063 Le Havre, France. and Universite du Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | | | - Sangeeta Sur
- Missouri University of Science and Technology, Rolla, MO 65409-0010, USA.
| | - Mathurin Robin
- Laboratoire Ondes et Milieux Complexes, UMR 6294, Centre National de la Recherche Scientifique-Universite du Havre, F-76063 Le Havre, France.
| | - François Lique
- Laboratoire Ondes et Milieux Complexes, UMR 6294, Centre National de la Recherche Scientifique-Universite du Havre, F-76063 Le Havre, France. and Universite du Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | - Richard Dawes
- Missouri University of Science and Technology, Rolla, MO 65409-0010, USA.
| |
Collapse
|
13
|
Abstract
The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H2O)2, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.
Collapse
Affiliation(s)
- Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA;
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA;
| |
Collapse
|
14
|
Desrousseaux B, Quintas-Sánchez E, Dawes R, Marinakis S, Lique F. Collisional excitation of interstellar PN by H 2: New interaction potential and scattering calculations. J Chem Phys 2021; 154:034304. [PMID: 33499633 DOI: 10.1063/5.0039145] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Rotational excitation of interstellar PN molecules induced by collisions with H2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN-H2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm-1. The dissociation energies were found to be 40.19 cm-1 and 75.05 cm-1 for complexes of PN with ortho-H2 and para-H2, respectively. Integral cross sections for rotational excitation in PN-H2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.
Collapse
Affiliation(s)
- Benjamin Desrousseaux
- LOMC, UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, 76063 Le Havre Cedex, France
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Sarantos Marinakis
- School of Health, Sport and Bioscience, University of East London, Stratford Campus, Water Lane, London E15 4LZ, United Kingdom
| | - François Lique
- LOMC, UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, 76063 Le Havre Cedex, France
| |
Collapse
|
15
|
Abstract
The CO2-N2 complex is formed from two key components of Earth's atmosphere, and as such, has received some attention from both experimental and theoretical studies. On the theory side, a potential energy surface (PES) based on high level ab initio data was reported [Nasri et al., J. Chem. Phys., 2015, 142, 174301] and then used in more recently reported rovibrational calculations [Lara-Moreno et al., Phys. Chem. Chem. Phys., 2019, 21, 3550]. Accuracy of about 1 percent was achieved for calculated rotational transitions of the ground vibrational state of the complex, compared with previously reported microwave spectra. However, a very recent measurement of the geared bending mode frequency [Barclay et al., J. Chem. Phys., 2020, 153, 014303] recorded a value of 21.4 cm-1, which is wildly different from the corresponding calculated value of 45.9 cm-1. To provide some insight into this discrepancy, we have constructed a new more accurate PES, and used it to perform highly converged variational rovibrational calculations. Our new results yield a value of 21.1 cm-1 for that bending frequency, in close agreement with the experiment. We also obtain significantly improved predicted rotational transitions. Finally, we note that a very shallow well, previously reported as a distinct second isomer, is not found on our new PES, but rather a transition structure is seen in that location.
Collapse
|
16
|
Quintas-Sánchez E, Dawes R, Lee K, McCarthy MC. Automated Construction of Potential Energy Surfaces Suitable to Describe van der Waals Complexes with Highly Excited Nascent Molecules: The Rotational Spectra of Ar-CS( v) and Ar-SiS( v). J Phys Chem A 2020; 124:4445-4454. [PMID: 32368913 DOI: 10.1021/acs.jpca.0c02685] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Some reactions produce extremely hot nascent products which nevertheless can form sufficiently long-lived van der Waals (vdW) complexes-with atoms or molecules from a bath gas-as to be observed via microwave spectroscopy. Theoretical calculations of such unbound resonance states can be much more challenging than ordinary bound-state calculations depending on the approach employed. One encounters not just the floppy, and perhaps multiwelled potential energy surface (PES) characteristic of vdWs complexes, but in addition, one must contend with excitation of the intramolecular modes and its corresponding influence on the PES. Straightforward computation of the (resonance) rovibrational levels of interest, involves the added complication of the unbound nature of the wave function, often treated with techniques such as introducing a complex absorbing potential. Here, we have demonstrated that a simplified approach of making a series of vibrationally effective PESs for the intermolecular coordinates-one for each reaction product vibrational quantum number of interest-can produce vdW levels for the complex with spectroscopic accuracy. This requires constructing a series of appropriately weighted lower-dimensional PESs for which we use our freely available PES-fitting code AUTOSURF. The applications of this study are the Ar-CS and Ar-SiS complexes, which are isovalent to Ar-CO and Ar-SiO, the latter of which we considered in a previously reported study. Using a series of vibrationally effective PESs, rovibrational levels and predicted microwave transition frequencies for both complexes were computed variationally. A series of shifting rotational transition frequencies were also computed as a function of the diatom vibrational quantum number. The predicted transitions were used to guide and inform an experimental effort to make complementary observations. Comparisons are given for the transitions that are within the range of the spectrometer and were successfully recorded. Calculations of the rovibrational level pattern agree to within 0.2% with experimental measurements.
Collapse
Affiliation(s)
- Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - Kelvin Lee
- Harvard-Smithsonian Center for Astrophysics, Cambridge, Massachusetts 02138, United States
| | - Michael C McCarthy
- Harvard-Smithsonian Center for Astrophysics, Cambridge, Massachusetts 02138, United States.,School of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
17
|
Ben Khalifa M, Quintas-Sánchez E, Dawes R, Hammami K, Wiesenfeld L. Rotational quenching of an interstellar gas thermometer: CH 3CNHe collisions. Phys Chem Chem Phys 2020; 22:17494-17502. [PMID: 32716451 DOI: 10.1039/d0cp02985h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Among all the molecular species found in the interstellar medium, molecules with threefold symmetry axes play a special role, as their rotational spectroscopy allows them to act as practical gas thermometers. Methyl-cyanide (CH3CN) is the second most abundant of those (after ammonia). We compute in this paper the collisional dynamics of methyl-cyanide in collision with helium, for both the A- and the E-symmetries of CH3CN. The potential energy surface is determined using the CCSD(T)-F12b formalism and fit with convenient analytic functions. We compute the rotationally inelastic cross sections for all levels up to 510 cm-1 of collision energy, employing at low energy exact Coupled Channels methods, and at higher energies, approximate Coupled States methods. For temperatures from 7 K up to 300 K, rates of quenching are computed and most are found to differ from those reported earlier (up to a factor of a thousand), calling for a possible reexamination of the temperatures assigned to low density gasses.
Collapse
Affiliation(s)
- M Ben Khalifa
- Laboratoire Aimé-Cotton, CNRS and Université Paris-Saclay, Orsay, France. and Laboratoire de Spectroscopie Atomique Moléculaire et Applications, Université Tunis El-Manar, Tunis, Tunisie.
| | - E Quintas-Sánchez
- Missouri University of Science and Technology, Rolla, MO 65409-0010, USA
| | - R Dawes
- Missouri University of Science and Technology, Rolla, MO 65409-0010, USA
| | - K Hammami
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications, Université Tunis El-Manar, Tunis, Tunisie.
| | - L Wiesenfeld
- Laboratoire Aimé-Cotton, CNRS and Université Paris-Saclay, Orsay, France.
| |
Collapse
|
18
|
Sur S, Ndengué SA, Quintas-Sánchez E, Bop C, Lique F, Dawes R. Rotationally inelastic scattering of O3–Ar: state-to-state rates with the multiconfigurational time dependent Hartree method. Phys Chem Chem Phys 2020; 22:1869-1880. [DOI: 10.1039/c9cp06501f] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The rates of state-changing collisions are compared for different isotopologues of ozone from quantum scattering calculations with the MCTDH method.
Collapse
Affiliation(s)
- Sangeeta Sur
- Department of Chemistry
- Missouri University of Science & Technology
- Rolla
- USA
| | - Steve A. Ndengué
- Department of Chemistry
- Missouri University of Science & Technology
- Rolla
- USA
- ICTP-East African Institute for Fundamental Research
| | | | - Cheikh Bop
- LOMC – UMR 6294
- CNRS-Université du Havre
- F-76063 Le Havre
- France
| | - François Lique
- LOMC – UMR 6294
- CNRS-Université du Havre
- F-76063 Le Havre
- France
| | - Richard Dawes
- Department of Chemistry
- Missouri University of Science & Technology
- Rolla
- USA
| |
Collapse
|
19
|
Castro-Juárez E, Wang XG, Carrington T, Quintas-Sánchez E, Dawes R. Computational study of the ro-vibrational spectrum of CO-CO 2. J Chem Phys 2019; 151:084307. [PMID: 31470713 DOI: 10.1063/1.5119762] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An accurate ab initio ground-state intermolecular potential energy surface (PES) was determined for the CO-CO2 van der Waals dimer. The Lanczos algorithm was used to compute rovibrational energies on this PES. For both the C-in and O-in T-shaped isomers, the fundamental transition frequencies agree well with previous experimental results. We confirm that the in-plane states previously observed are geared states. In addition, we have computed and assigned many other vibrational states. The rotational constants we determine from J = 1 energy levels agree well with their experimental counterparts. Planar and out-of-plane cuts of some of the wavefunctions we compute are quite different, indicating strong coupling between the bend and torsional modes. Because the stable isomers are T-shaped, vibration along the out-of-plane coordinates is very floppy. In CO-CO2, when the molecule is out-of-plane, interconversion of the isomers is possible, but the barrier height is higher than the in-plane geared barrier height.
Collapse
Affiliation(s)
| | - Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| |
Collapse
|
20
|
Desrousseaux B, Quintas-Sánchez E, Dawes R, Lique F. Collisional Excitation of CF+ by H2: Potential Energy Surface and Rotational Cross Sections. J Phys Chem A 2019; 123:9637-9643. [DOI: 10.1021/acs.jpca.9b05538] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Benjamin Desrousseaux
- LOMC - UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, F-76063 Le Havre, France
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - François Lique
- LOMC - UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, F-76063 Le Havre, France
| |
Collapse
|
21
|
Sur S, Quintas-Sánchez E, Ndengué SA, Dawes R. Development of a potential energy surface for the O3–Ar system: rovibrational states of the complex. Phys Chem Chem Phys 2019; 21:9168-9180. [DOI: 10.1039/c9cp01044k] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Collisional stabilization is an important step in the process of atmospheric formation of ozone.
Collapse
Affiliation(s)
- Sangeeta Sur
- Department of Chemistry
- Missouri University of Science & Technology
- Rolla
- USA
| | | | - Steve A. Ndengué
- Department of Chemistry
- Missouri University of Science & Technology
- Rolla
- USA
| | - Richard Dawes
- Department of Chemistry
- Missouri University of Science & Technology
- Rolla
- USA
| |
Collapse
|
22
|
Affiliation(s)
- Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65401, United States
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65401, United States
| |
Collapse
|
23
|
Dubernet ML, Quintas-Sánchez E, Tuckey P. New potential energy surface for the HCS(+)-He system and inelastic rate coefficients. J Chem Phys 2015; 143:044315. [PMID: 26233137 DOI: 10.1063/1.4926839] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
A new high quality potential energy surface is calculated at a coupled-cluster single double triple level with an aug-cc-pV5Z basis set for the HCS(+)-He system. This potential energy surface is used in low energy quantum scattering calculations to provide a set of (de)-excitation cross sections and rate coefficients among the first 20 rotational levels of HCS(+) by He in the range of temperature from 5 K to 100 K. The paper discusses the impact of the new ab initio potential energy surface on the cross sections at low energy and provides a comparison with the HCO(+)-He system. The HCS(+)-He rate coefficients for the strongest transitions differ by factors of up to 2.5 from previous rate coefficients; thus, analysis of astrophysical spectra should be reconsidered with the new rate coefficients.
Collapse
Affiliation(s)
- Marie-Lise Dubernet
- LERMA, UMR8112, PSL Research University, Observatoire de Paris, Sorbonne Universités, UPMC Univ Paris 06, CNRS, 5 Place Janssen, 92195 Meudon, France
| | - Ernesto Quintas-Sánchez
- LERMA, UMR8112, PSL Research University, Observatoire de Paris, Sorbonne Universités, UPMC Univ Paris 06, CNRS, 5 Place Janssen, 92195 Meudon, France
| | - Philip Tuckey
- LNE-SYRTE, PSL Research University, Observatoire de Paris, Sorbonne Universités, UPMC Univ Paris 06, CNRS, LNE, 61 Av. de l'Observatoire, 75014 Paris, France
| |
Collapse
|
24
|
Quintas-Sánchez E, Crespos C, Larrégaray P, Rayez JC, Martin-Gondre L, Rubayo-Soneira J. Surface temperature effects on the dynamics of N2 Eley-Rideal recombination on W(100). J Chem Phys 2013; 138:024706. [PMID: 23320712 DOI: 10.1063/1.4774024] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quasiclassical trajectories simulations are performed to study the influence of surface temperature on the dynamics of a N atom colliding a N-preadsorbed W(100) surface under normal incidence. A generalized Langevin surface oscillator scheme is used to allow energy transfer between the nitrogen atoms and the surface. The influence of the surface temperature on the N(2) formed molecules via Eley-Rideal recombination is analyzed at T = 300, 800, and 1500 K. Ro-vibrational distributions of the N(2) molecules are only slightly affected by the presence of the thermal bath whereas kinetic energy is rather strongly decreased when going from a static surface model to a moving surface one. In terms of reactivity, the moving surface model leads to an increase of atomic trapping cross section yielding to an increase of the so-called hot atoms population and a decrease of the direct Eley-Rideal cross section. The energy exchange between the surface and the nitrogen atoms is semi-quantitatively interpreted by a simple binary collision model.
Collapse
|
25
|
Quintas-Sánchez E, Larrégaray P, Crespos C, Martin-Gondre L, Rubayo-Soneira J, Rayez JC. Dynamical reaction pathways in Eley-Rideal recombination of nitrogen from W(100). J Chem Phys 2012; 137:064709. [DOI: 10.1063/1.4742815] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|