1
|
Guo CM, Zhang H, Cheng XL. Exhaustive State-to-State Cross Sections and Rate Coefficients for Inelastic N 2-N 2 Collisions using QCT Combined with Neural Network Models. J Phys Chem A 2024; 128:5435-5444. [PMID: 38953499 DOI: 10.1021/acs.jpca.4c00590] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Using the quasi-classical trajectory method, we systematically studied the state-to-state vibrational relaxation process of N2(v1) + N2(v2) collisions over a wide temperature range (5000-30,000 K). Different temperature dependencies of the single- and multiquantum VV and VT events in various (v1,v2) collisions are captured, with the dominant channel being related to the initial vibrational energy levels (vmax = 50). At a specified relative translational energy, there is a monotonic relationship of the VT cross sections with the vibrational energy level, particularly in high-energy collisions. Additionally, we constructed well-trained neural network models (R-values reaching 0.99) using limited quasi-classical trajectory (QCT) data sets, which can be used to predict the state-to-state cross sections and rate coefficients of the VV processes N2(v1) + N2(v2) → N2(v1 - Δv) + N2(v2 + Δv) and VT processes N2(v1) + N2(v2) → N2(v1 - Δv) + N2(v2) (Δv = ±1, ±2, ±3) for collisions with arbitrary initial vibrational states. This work not only significantly reduces computational resources but also serves as a reference for the study of the state-to-state dynamics of all four-atom collision systems in hypersonic flows.
Collapse
Affiliation(s)
- Chang-Min Guo
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Hong Zhang
- College of Physics, Sichuan University, Chengdu 610065, China
| | - Xin-Lu Cheng
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| |
Collapse
|
2
|
Yang J, Li J, Li J, Li J. Gaussian Process Regression for State-to-State Integral Cross Sections: The Case of the O + O 2 Collision Dissociation Reactions. J Phys Chem A 2024; 128:4966-4975. [PMID: 38869143 DOI: 10.1021/acs.jpca.4c01445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2024]
Abstract
Research on hypersonic vehicles has become increasingly important worldwide in recent years. However, accurately simulating the dynamics of the nonequilibrium high-temperature reactions that are in the hypersonic flow around the vehicles presents a significant challenge as a large number of states and transitions are accessible even for the smallest atom-diatom reaction systems. It is quite difficult, sometimes even impossible, to exhaustively investigate all relevant combinations or determine high-dimensional analytical representations for the state-to-state reaction probabilities. In this study, we used Gaussian process regression (GPR) to fit a model based on only 807 QCT data for training. The confidence interval of the GPR prediction and the Kullback-Leibler (KL) divergence were used to help minimize the sampling amount of data for fitting the converged GPR model. The model aims to predict the state-to-state integral cross section (ICS) of the O + O2 → 3O dissociation reaction under random initial conditions (Et, v, j). In total, it took almost a month to obtain this converged GPR model, but it took only a few seconds to predict the ICS value for any initial condition. For 330 initial conditions not included in the training set, the mean-square error (MSE) between the QCT-calculated ICSs and the GPR-predicted ones is only 0.08 Å2 and the R2 is 0.9986, indicating that the GPR model can replace the direct expensive QCT calculation with high accuracy. Finally, we calculated the equilibrium dissociation rate coefficients based on the StS ICS values predicted by the GPR model, and the results were in good agreement with available experimental and theoretical results. Thus, this study provides an effective and accurate approach to the extensive direct state-to-state reaction dynamic calculations.
Collapse
Affiliation(s)
- Jiawei Yang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Jia Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Junhong Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| |
Collapse
|
3
|
Hong Q, Storchi L, Coletti C, Li J, Sun Q, Li J. Quantum-classical rate coefficient datasets of vibrational energy transfer in carbon monoxide based on highly accurate potential energy surface. J Chem Phys 2024; 160:084305. [PMID: 38411231 DOI: 10.1063/5.0189772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Accepted: 02/02/2024] [Indexed: 02/28/2024] Open
Abstract
A merged potential energy surface (PES) is introduced for CO + CO collisions by combining a recent full-dimensional ab initio PES [Chen et al. J. Chem. Phys. 153, 054310 (2020)] and analytical long-range multipolar interactions. This merged PES offers a double advantage: it retains the precision of the ab initio PES in describing the van der Waals well and repulsive short range while providing an accurate physical description of long-range interaction; it significantly reduces the computational time required for trajectory integration since the long-range portion of the ab initio PES (involving numerous neural network fitting parameters) is now replaced by the analytical model potential. Based on the present merged PES, mixed Quantum-Classical (MQC) calculations, which capture quantum effects related to vibrational motion, align with a range of experimental data, including transport properties, vibrational energy transfer between CO and its isotoplogues, as well as rate coefficients for V-V and V-T/R processes. Notably, the original ab initio PES yields V-T/R rate coefficients at low temperatures that are significantly higher than the experimental data due to the artificial contribution of its unphysical long-range potential. In addition to conducting extensive MQC calculations to obtain raw data for V-V and V-T/R rate coefficients, we employ Gaussian process regression to predict processes lacking computed MQC data, thereby completing the considered V-V and V-T/R datasets. These extensive rate coefficient datasets, particularly for V-T/R processes, are unprecedented and reveal the significant role played by V-T/R processes at high temperatures, emphasizing the necessity of incorporating both V-V and V-T/R processes in the applications.
Collapse
Affiliation(s)
- Qizhen Hong
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190 Beijing, China
| | - Loriano Storchi
- Dipartimento di Farmacia, Università G. d'Annunzio Chieti-Pescara, via dei Vestini, 66100 Chieti, Italy
| | - Cecilia Coletti
- Dipartimento di Farmacia, Università G. d'Annunzio Chieti-Pescara, via dei Vestini, 66100 Chieti, Italy
| | - Jia Li
- School of Chemistry and Chemical Engineering and Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Quanhua Sun
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190 Beijing, China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jun Li
- School of Chemistry and Chemical Engineering and Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| |
Collapse
|
4
|
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] [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
|
5
|
He D, Hong Q, Li F, Sun Q, Si T, Luo X. Experimental and numerical studies on the thermal nonequilibrium behaviors of CO with Ar, He, and H2. J Chem Phys 2023; 159:234302. [PMID: 38108486 DOI: 10.1063/5.0176176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Accepted: 11/24/2023] [Indexed: 12/19/2023] Open
Abstract
The time-dependent rotational and vibrational temperatures were measured to study the shock-heated thermal nonequilibrium behaviors of CO with Ar, He, and H2 as collision partners. Three interference-free transition lines in the fundamental vibrational band of CO were applied to the fast, in situ, and state-specific measurements. Vibrational relaxation times of CO were summarized over a temperature range of 1110-2820 K behind reflected shocks. The measured rotational temperature instantaneously reached an equilibrium state behind shock waves. The measured vibrational temperature experienced a relaxation process before reaching the equilibrium state. The measured vibrational temperature time histories were compared with predictions based on the Landau-Teller model and the state-to-state approach. The state-to-state approach treats the vibrational energy levels of CO as pseudo-species and accurately describes the detailed thermal nonequilibrium processes behind shock waves. The datasets of state-specific inelastic rate coefficients of CO-Ar, CO-He, CO-CO, and CO-H2 collisions were calculated in this study using the mixed quantum-classical method and the semiclassical forced harmonic oscillator model. The predictions based on the state-to-state approach agreed well with the measured data and nonequilibrium (non-Boltzmann) vibrational distributions were found in the post-shock regions, while the Landau-Teller model predicted slower vibrational temperature time histories than the measured data. Modifications were applied to the Millikan-White vibrational relaxation data of the CO-Ar and CO-H2 systems to improve the performance of the Landau-Teller model. In addition, the thermal nonequilibrium processes behind incident shocks, the acceleration effects of H2O on the relaxation process of CO, and the characterization of vibrational temperature were highlighted.
Collapse
Affiliation(s)
- Dong He
- Deep Space Exploration Laboratory, Department of Modern Mechanics, University of Science and Technology of China, Hefei 230026, People's Republic of China
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | - Qizhen Hong
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | - Fei Li
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | - Quanhua Sun
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Ting Si
- Deep Space Exploration Laboratory, Department of Modern Mechanics, University of Science and Technology of China, Hefei 230026, People's Republic of China
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | - Xisheng Luo
- Deep Space Exploration Laboratory, Department of Modern Mechanics, University of Science and Technology of China, Hefei 230026, People's Republic of China
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| |
Collapse
|
6
|
Hong Q, Storchi L, Sun Q, Bartolomei M, Pirani F, Coletti C. Improved Quantum-Classical Treatment of N 2-N 2 Inelastic Collisions: Effect of the Potentials and Complete Rate Coefficient Data Sets. J Chem Theory Comput 2023; 19:8557-8571. [PMID: 38007713 PMCID: PMC10720385 DOI: 10.1021/acs.jctc.3c01103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 11/11/2023] [Accepted: 11/13/2023] [Indexed: 11/28/2023]
Abstract
In this study, complete (i.e., including all vibrational quantum numbers in an N2 vibrational ladder) data sets of vibration-to-vibration and vibration-to-translation rate coefficients for N2-N2 collisions are explicitly computed along with transport properties (shear and bulk viscosity, thermal conductivity, and self-diffusion) in the temperature range 100-9000 K. To reach this goal, we improved a mixed quantum-classical (MQC) dynamics approach by lifting the constraint of a Morse treatment of the vibrational wave function and intramolecular potential and permitting the use of more realistic and flexible representations. The new formulation has also allowed us to separately analyze the role of intra- and intermolecular potentials on the calculated rates and properties. Ab initio intramolecular potentials are indispensable for highly excited vibrational states, though the Morse potential still gives reasonable values up to v = 20. An accurate description of the long-range interaction and the van der Waals well is a requisite for the correct reproduction of qualitative and quantitative rate coefficients, particularly at low temperatures, making physically meaningful analytical representations still the best choice compared to currently available ab initio potential energy surfaces. These settings were used to directly compute the MQC rates corresponding to a large number of initial vibrational quantum numbers, and the missing intermediate values were predicted using a machine learning technique (i.e., the Gaussian process regression approach). The obtained values are reliable in the wide temperature range employed and are therefore valuable data for many communities dealing with nonlocal thermal equilibrium conditions in different environments.
Collapse
Affiliation(s)
- Qizhen Hong
- State
Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190 Beijing, China
| | - Loriano Storchi
- Dipartimento
di Farmacia, Università Gabriele
d’Annunzio Chieti-Pescara, via dei Vestini, 66100 Chieti, Italy
| | - Quanhua Sun
- State
Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190 Beijing, China
- School
of Engineering Science, University of Chinese
Academy of Sciences, 100049 Beijing, China
| | | | - Fernando Pirani
- Dipartimento
di Chimica, Biologia e Biotecnologie, Università
di Perugia, via Elce di Sotto 8, 06123 Perugia, Italy
| | - Cecilia Coletti
- Dipartimento
di Farmacia, Università Gabriele
d’Annunzio Chieti-Pescara, via dei Vestini, 66100 Chieti, Italy
| |
Collapse
|
7
|
Del Fré S, Santamaría AR, Duflot D, Basalgète R, Féraud G, Bertin M, Fillion JH, Monnerville M. Mechanism of Ultraviolet-Induced CO Desorption from CO Ice: Role of Vibrational Relaxation Highlighted. PHYSICAL REVIEW LETTERS 2023; 131:238001. [PMID: 38134796 DOI: 10.1103/physrevlett.131.238001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 10/13/2023] [Indexed: 12/24/2023]
Abstract
Although UV photon-induced CO ice desorption is clearly observed in many cold regions of the Universe as well as in the laboratory, the fundamental question of the mechanisms involved at the molecular scale remains debated. In particular, the exact nature of the involved energy transfers in the indirect desorption pathway highlighted in previous experiments is not explained. Using ab initio molecular dynamics simulations, we explore a new indirect desorption mechanism in which a highly vibrationally excited CO (v=40) within an aggregate of 50 CO molecules triggers the desorption of molecules at the surface. The desorption originates first from a mutual attraction between the excited molecule and the surrounding molecule(s), followed by a cascade of energy transfers, ultimately resulting in the desorption of vibrationally cold CO (∼95% in v=0). The theoretical vibrational distribution, along with the kinetic energy one, which peaks around 25 meV for CO with low rotational levels (v=0, J<7), is in excellent agreement with the results obtained from VUV laser induced desorption (157 nm) of CO (v=0, 1) probed using REMPI.
Collapse
Affiliation(s)
- Samuel Del Fré
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | | | - Denis Duflot
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Romain Basalgète
- Sorbonne Université, Observatoire de Paris, Université PSL, CNRS, LERMA, F-75005 Paris, France
| | - Géraldine Féraud
- Sorbonne Université, Observatoire de Paris, Université PSL, CNRS, LERMA, F-75005 Paris, France
| | - Mathieu Bertin
- Sorbonne Université, Observatoire de Paris, Université PSL, CNRS, LERMA, F-75005 Paris, France
| | - Jean-Hugues Fillion
- Sorbonne Université, Observatoire de Paris, Université PSL, CNRS, LERMA, F-75005 Paris, France
| | - Maurice Monnerville
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| |
Collapse
|
8
|
He N, Huang M, Evangelista FA. CO Inversion on a NaCl(100) Surface: A Multireference Quantum Embedding Study. J Phys Chem A 2023; 127:1975-1987. [PMID: 36799901 PMCID: PMC9986868 DOI: 10.1021/acs.jpca.2c05844] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
We develop a multireference quantum embedding model to investigate a recent experimental observation of the isomerization of vibrationally excited CO molecules on a NaCl(100) surface [Science 2020, 367, 175-178]. To explore this mechanism, we built a reduced potential energy surface of CO interacting with NaCl(100) using a second-order multireference perturbation theory, modeling the adsorbate-surface interaction with our previously developed active space embedding theory (ASET). We considered an isolated CO molecule on NaCl(100) and a high-coverage CO monolayer (1/1), and for both we generated potential energy surfaces parametrized by the CO stretching, adsorption, and inversion coordinates. These surfaces are used to determine stationary points and adsorption energies and to perform a vibrational analysis of the states relevant to the inversion mechanism. We found that for near-equilibrium bond lengths, CO adsorbed in the C-down configuration is lower in energy than in the O-down configuration. Stretching of the C-O bond reverses the energetic order of these configurations, supporting the accepted isomerization mechanism. The vibrational constants obtained from these potential energy surfaces show a small (< 10 cm-1) blue- and red-shift for the C-down and O-down configurations, respectively, in agreement with experimental assignments and previous theoretical studies. Our vibrational analysis of the monolayer case suggests that the O-down configuration is energetically more stable than the C-down one beyond the 16th vibrational excited state of CO, a value slightly smaller than the one from quasi-classical trajectory simulations (22nd) and consistent with the experiment. Our analysis suggests that CO-CO interactions in the monolayer play an important role in stabilizing highly vibrationally excited states in the O-down configuration and reducing the barrier between the C-down and O-down geometries, therefore playing a crucial role in the inversion mechanism.
Collapse
Affiliation(s)
- Nan He
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
9
|
Töpfer K, Upadhyay M, Meuwly M. Quantitative molecular simulations. Phys Chem Chem Phys 2022; 24:12767-12786. [PMID: 35593769 PMCID: PMC9158373 DOI: 10.1039/d2cp01211a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 04/30/2022] [Indexed: 11/21/2022]
Abstract
All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.
Collapse
Affiliation(s)
- Kai Töpfer
- 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.
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| |
Collapse
|
10
|
Żółtowski M, Loreau J, Lique F. Collisonal energy transfer in the CO-CO system. Phys Chem Chem Phys 2022; 24:11910-11918. [PMID: 35510882 PMCID: PMC9116445 DOI: 10.1039/d2cp01065h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An accurate determination of the physical conditions in astrophysical environments relies on the modeling of molecular spectra. In such environments, densities can be so low ($n << 10^{10}$ cm$^{-3}$) that...
Collapse
Affiliation(s)
- Michał Żółtowski
- LOMC - UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, F-76063 Le Havre, France.
- 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
| | - François Lique
- LOMC - UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, F-76063 Le Havre, France.
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| |
Collapse
|
11
|
Lu D, Chen J, Guo H, Li J. Vibrational energy pooling via collisions between asymmetric stretching excited CO 2: a quasi-classical trajectory study on an accurate full-dimensional potential energy surface. Phys Chem Chem Phys 2021; 23:24165-24174. [PMID: 34671798 DOI: 10.1039/d1cp03687d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In low temperature plasmas, energy transfer between asymmetric stretching excited CO2 molecules can be highly efficient, which leads to further excitation (and de-excitation) of the CO2 molecules: CO2(vas) + CO2(vas) → CO2(vas + 1) + CO2(vas - 1). Through such a vibrational ladder climbing mechanism, CO2 can be activated and eventually dissociates. To gain mechanistic insight of such processes, a full-dimensional accurate potential energy surface (PES) for the CO2 + CO2 system is developed using the permutational invariant polynomial-neural network method based on CCSD(T)-F12a/AVTZ energies at about 39 000 geometries. This PES is used in quasi-classical trajectory (QCT) studies of the vibrational energy transfer between CO2 molecules excited in the asymmetric stretching mode. A machine learning algorithm is used to determine state-specific rate coefficients for the vibrational transfer processes from a limited data set. In addition to the CO2(vas + 1) + CO2(vas - 1) channel, the QCT simulations revealed significant contributions from the CO2(vas + 2,3) + CO2(vas - 2,3) channels, particularly at low collision energies/temperatures. These multi-vibrational-quantum processes are attributed to enhanced energy flow in the collisional complex formed by enhanced dipole-dipole interaction between asymmetric stretching excited CO2 molecules.
Collapse
Affiliation(s)
- Dandan Lu
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China. .,Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131, USA
| | - Jun Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, Fujian 350002, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131, USA
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| |
Collapse
|
12
|
Quasi-classical trajectory study of inelastic collision energy transfer between H2CO and H2 on a full-dimensional potential energy surface. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.139014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
13
|
Liu Q, Liu L, An F, Huang J, Zhou Y, Xie D. A full-dimensional ab initio intermolecular potential energy surface and rovibrational spectra for OC-HF and OC-DF. J Chem Phys 2021; 155:084302. [PMID: 34470366 DOI: 10.1063/5.0061291] [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/14/2022] Open
Abstract
We present a full-dimensional ab initio intermolecular potential energy surface (IPES) for the OC-HF van der Waals complex. 3167 ab initio points were computed at the frozen-core (FC) explicitly correlated coupled cluster [FC-CCSD(T)-F12b] level, with the augmented correlation-consistent polarized valence quadruple-zeta basis set plus bond functions. Basis set superposition error correction was also considered by the full counterpoise procedure. Gaussian process regression (GPR) was used to map out the potential energy surface, while a multipole expansion method was employed to smooth the ab initio noise of intermolecular potential in the long range. The global minimum of -1248.364 cm-1 was located at the linear configuration with the C atom pointing toward the H atom of the HF molecule. In addition, a local minimum of -602.026 cm-1 was found at another linear configuration with the O atom pointing toward the H atom of the HF molecule. The eigenstates were calculated on the vibrational averaged four-dimensional IPESs with the mixed radial discrete variable representation/angular finite basis representation method and Lanczos propagation algorithm. The dissociation energy D0 was calculated to be 701.827 cm-1, well reproducing the experimental value of 732 ± 2 cm-1. The dipole moment surfaces were also fitted by GPR from 3132 ab initio points calculated using the coupled cluster method [CCSD(T)] with AVTZ basis set plus bond functions. The frequencies and relative line intensities of rovibrational transitions in the HF (DF) and CO stretching bands were further calculated and compared well with the experimental results. These results indicate the high fidelity of the new IPES.
Collapse
Affiliation(s)
- Qiong Liu
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Lu Liu
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Feng An
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Jing Huang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yanzi Zhou
- 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
| |
Collapse
|
14
|
Auerbach DJ, Tully JC, Wodtke AM. Chemical dynamics from the gas‐phase to surfaces. ACTA ACUST UNITED AC 2021. [DOI: 10.1002/ntls.10005] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Daniel J. Auerbach
- Institut für physikalische Chemie Georg‐August Universität Göttingen Göttingen Germany
- Abteilung für Dynamik an Oberflächen Max‐Planck‐Institut für biophysikalische Chemie Göttingen Germany
| | - John C. Tully
- Department of Chemistry Yale University New Haven Connecticut USA
| | - Alec M. Wodtke
- Institut für physikalische Chemie Georg‐August Universität Göttingen Göttingen Germany
- Abteilung für Dynamik an Oberflächen Max‐Planck‐Institut für biophysikalische Chemie Göttingen Germany
| |
Collapse
|
15
|
Ceriotti M, Clementi C, Anatole von Lilienfeld O. Machine learning meets chemical physics. J Chem Phys 2021; 154:160401. [PMID: 33940847 DOI: 10.1063/5.0051418] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Over recent years, the use of statistical learning techniques applied to chemical problems has gained substantial momentum. This is particularly apparent in the realm of physical chemistry, where the balance between empiricism and physics-based theory has traditionally been rather in favor of the latter. In this guest Editorial for the special topic issue on "Machine Learning Meets Chemical Physics," a brief rationale is provided, followed by an overview of the topics covered. We conclude by making some general remarks.
Collapse
Affiliation(s)
- Michele Ceriotti
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Cecilia Clementi
- Department of Physics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | | |
Collapse
|
16
|
Morita M, Zuo J, Guo H, Balakrishnan N. Rainbow scattering in rotationally inelastic collisions of HCl and H 2. J Chem Phys 2021; 154:104304. [PMID: 33722024 DOI: 10.1063/5.0043658] [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
We examine rotational transitions of HCl in collisions with H2 by carrying out quantum mechanical close-coupling and quasi-classical trajectory (QCT) calculations on a recently developed globally accurate full-dimensional ab initio potential energy surface for the H3Cl system. Signatures of rainbow scattering in rotationally inelastic collisions are found in the state resolved integral and differential cross sections as functions of the impact parameter (initial orbital angular momentum) and final rotational quantum number. We show the coexistence of distinct dynamical regimes for the HCl rotational transition driven by the short-range repulsive and long-range attractive forces whose relative importance depends on the collision energy and final rotational state, suggesting that the classification of rainbow scattering into rotational and l-type rainbows is effective for H2 + HCl collisions. While the QCT method satisfactorily predicts the overall behavior of the rotationally inelastic cross sections, its capability to accurately describe signatures of rainbow scattering appears to be limited for the present system.
Collapse
Affiliation(s)
- Masato Morita
- Department of Chemistry and Biochemistry, University of Nevada, Las Vegas, Las Vegas, Nevada 89154, USA
| | - Junxiang Zuo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Naduvalath Balakrishnan
- Department of Chemistry and Biochemistry, University of Nevada, Las Vegas, Las Vegas, Nevada 89154, USA
| |
Collapse
|
17
|
Quasiclassical simulations based on cluster models reveal vibration-facilitated roaming in the isomerization of CO adsorbed on NaCl. Nat Chem 2021; 13:249-254. [PMID: 33462381 DOI: 10.1038/s41557-020-00612-y] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2020] [Accepted: 11/18/2020] [Indexed: 11/08/2022]
Abstract
The desire to better understand the quantum nature of isomerization led to recent experimental observations of the vibrationally induced isomerization of OC-NaCl(100) to CO-NaCl(100). To investigate the mechanism of this isomerization, we performed dynamics calculations using finite (CO-NaCl)n cluster models. We constructed new potential energy surfaces for CO-NaCl and CO-CO interactions using high-level ab initio data and report key properties of the bare CO-NaCl potential energy surface, which show much in common with the experiment. We investigated the isomerization dynamics using several cluster models and, in all cases, isomerization was seen for highly excited CO vibrational states, in agreement with experiments. A detailed examination of the reaction trajectories indicates that isomerization occurs when the distance between CO and NaCl is larger than the distance at the conventional isomerization saddle point, which is a strong indicator of 'roaming'.
Collapse
|
18
|
Mandal B, Semenov A, Babikov D. Adiabatic Trajectory Approximation within the Framework of Mixed Quantum/Classical Theory. J Phys Chem A 2020; 124:9877-9888. [DOI: 10.1021/acs.jpca.0c07547] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bikramaditya Mandal
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, United States
| | - Alexander Semenov
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, United States
| | - Dmitri Babikov
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, United States
| |
Collapse
|