Modeling of Energy Transfer From Vibrationally Excited CO2 Molecules: Cross Sections and Probabilities for Kinetic Modeling of Atmospheres, Flows, and Plasmas

Vibrational energy pooling via collisions between asymmetric stretching excited CO2: a quasi-classical trajectory study on an accurate full-dimensional potential energy surface

In low temperature plasmas, energy transfer between asymmetric stretching excited CO₂ molecules can be highly efficient, which leads to further excitation (and de-excitation) of the CO₂ molecules: CO₂(v_as) + CO₂(v_as) → CO₂(v_as + 1) + CO₂(v_as - 1). Through such a vibrational ladder climbing mechanism, CO₂ can be activated and eventually dissociates. To gain mechanistic insight of such processes, a full-dimensional accurate potential energy surface (PES) for the CO₂ + CO₂ 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 CO₂ 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 CO₂(v_as + 1) + CO₂(v_as - 1) channel, the QCT simulations revealed significant contributions from the CO₂(v_as + 2,3) + CO₂(v_as - 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 CO₂ molecules.

Collision-Induced Rotational Excitation of CO2 by N(4S) Atoms: A New Ab Initio Potential Energy Surface and Scattering Calculations

Collisional excitations of CO₂ molecules are significant to fully understand the physical and chemical processes of astrophysical and atmospheric environments. Rotational excitations of CO₂ molecules induced by N(⁴S) atoms have been studied for the first time. First, we have computed a new highly accurate ab initio potential energy surface (PES) of a CO₂-N(⁴S) van der Waals complex. The PES has been obtained by employing the partially spin-restricted coupled cluster with open-shell single, double, and perturbative triple excitation method with aug-cc-pVQZ basis sets. The full close-coupling calculations have been performed to compute cross sections for kinetic energies up to 800 cm^-1. For all of the excitations, rotational cross sections exhibit an overall decrease with the increase of the energy gaps. Rate coefficients are calculated by averaging the cross sections over a Maxwell-Boltzmann distribution for temperatures ranging from 1 to 150 K. The trends in rate coefficients are in good agreement with those of similar collision systems. The decrease in energy gaps and the increase in temperature are the key factors to enhance the rate coefficients of CO₂ excitation. Our study will be useful for accurately establishing the atmospheric model of terrestrial planets and determining the abundance of CO₂ and N(⁴S) in space.

Chirality in molecular collision dynamics

Chirality is a phenomenon that permeates the natural world, with implications for atomic and molecular physics, for fundamental forces and for the mechanisms at the origin of the early evolution of life and biomolecular homochirality. The manifestations of chirality in chemistry and biochemistry are numerous, the striking ones being chiral recognition and asymmetric synthesis with important applications in molecular sciences and in industrial and pharmaceutical chemistry. Chiral discrimination phenomena, due to the existence of two enantiomeric forms, very well known in the case of interaction with light, but still nearly disregarded in molecular collision studies. Here we review some ideas and recent advances about the role of chirality in molecular collisions, designing and illustrating molecular beam experiments for the demonstration of chiral effects and suggesting a scenario for a stereo-directional origin of chiral selection.

Construction of a coarse-grain quasi-classical trajectory method. I. Theory and application to N2-N2 system

This work aims to construct a reduced order model for energy transfer and dissociation in non-equilibrium nitrogen mixtures. The objective is twofold: to present the Coarse-Grain Quasi-Classical Trajectory (CG-QCT) method, a novel framework for constructing a reduced order model for diatom-diatom systems; and to analyze the physics of non-equilibrium relaxation of the nitrogen molecules undergoing dissociation in an ideal chemical reactor. The CG-QCT method couples the construction of the reduced order model under the coarse-grain model framework with the quasi-classical trajectory calculations to directly construct the reduced model without the need for computing the individual rovibrational specific kinetic data. In the coarse-grain model, the energy states are lumped together into groups containing states with similar properties, and the distribution of states within each of these groups is prescribed by a Boltzmann distribution at the local translational temperature. The required grouped kinetic properties are obtained directly by the QCT calculations. Two grouping strategies are considered: energy-based grouping, in which states of similar internal energy are lumped together, and vibrational grouping, in which states with the same vibrational quantum number are grouped together. A zero-dimensional chemical reactor simulation, in which the molecules are instantaneously heated, forcing the system into strong non-equilibrium, is used to study the differences between the two grouping strategies. The comparison of the numerical results against available experimental data demonstrates that the energy-based grouping is more suitable to capture dissociation, while the energy transfer process is better described with a vibrational grouping scheme. The dissociation process is found to be strongly dependent on the behavior of the high energy states, which contribute up to 50% of the dissociating molecules. Furthermore, up to 40% of the energy required to dissociate the molecules comes from the rotational mode, underscoring the importance of accounting for this mode when constructing non-equilibrium kinetic models. In contrast, the relaxation process is governed primarily by low energy states, which exhibit significantly slower transitions in the vibrational binning model due to the prevalence of mode separation in these states.

Energy transfer dynamics and kinetics of elementary processes (promoted) by gas-phase CO2 -N2 collisions: Selectivity control by the anisotropy of the interaction

In this work, we exploit a new formulation of the potential energy and of the related computational procedures, which embodies the coupling between the intra and intermolecular components, to characterize possible propensities of the collision dynamics in energy transfer processes of interest for simulation and control of phenomena occurring in a variety of equilibrium and nonequilibrium environments. The investigation reported in the paper focuses on the prototype CO2 -N2 system, whose intramolecular component of the interaction is modeled in terms of a many body expansion while the intermolecular component is modeled in terms of a recently developed bonds-as-interacting-molecular-centers' approach. The main advantage of this formulation of the potential energy surface is that of being (a) truly full dimensional (i.e., all the variations of the coordinates associated with the molecular vibrations and rotations on the geometrical and electronic structure of the monomers, are explicitly taken into account without freezing any bonds or angles), (b) more flexible than other usual formulations of the interaction and (c) well suited for fitting procedures better adhering to accurate ab initio data and sensitive to experimental arrangement dependent information. Specific attention has been given to the fact that a variation of vibrational and rotational energy has a higher (both qualitative and quantitative) impact on the energy transfer when a more accurate formulation of the intermolecular interaction (with respect to that obtained when using rigid monomers) is adopted. This makes the potential energy surface better suited for the kinetic modeling of gaseous mixtures in plasma, combustion and atmospheric chemistry computational applications. © 2016 Wiley Periodicals, Inc.

Energy transfer upon collision of selectively excited CO2 molecules: State-to-state cross sections and probabilities for modeling of atmospheres and gaseous flows

Carbon dioxide molecules can store and release tens of kcal/mol upon collisions, and such an energy transfer strongly influences the energy disposal and the chemical processes in gases under the extreme conditions typical of plasmas and hypersonic flows. Moreover, the energy transfer involving CO2 characterizes the global dynamics of the Earth-atmosphere system and the energy balance of other planetary atmospheres. Contemporary developments in kinetic modeling of gaseous mixtures are connected to progress in the description of the energy transfer, and, in particular, the attempts to include non-equilibrium effects require to consider state-specific energy exchanges. A systematic study of the state-to-state vibrational energy transfer in CO2 + CO2 collisions is the focus of the present work, aided by a theoretical and computational tool based on quasiclassical trajectory simulations and an accurate full-dimension model of the intermolecular interactions. In this model, the accuracy of the description of the intermolecular forces (that determine the probability of energy transfer in molecular collisions) is enhanced by explicit account of the specific effects of the distortion of the CO2 structure due to vibrations. Results show that these effects are important for the energy transfer probabilities. Moreover, the role of rotational and vibrational degrees of freedom is found to be dominant in the energy exchange, while the average contribution of translations, under the temperature and energy conditions considered, is negligible. Remarkable is the fact that the intramolecular energy transfer only involves stretching and bending, unless one of the colliding molecules has an initial symmetric stretching quantum number greater than a threshold value estimated to be equal to 7.

An innovative synergistic grid approach to the computational study of protein aggregation mechanisms

Thanks to the advances in grid technologies, we are able to propose here an evolution of our molecular simulator that, when moving to larger systems, instead of reducing the granularity of the dynamical treatment (as is often done in molecular dynamics studies of such systems) exploits the extra power of the grid approach to the end of preserving the detailed nature of theatomistic formulation of the interaction. Key steps of such evolution are: (1) the assemblage of the interaction based on a composition of the ab initio intramolecular data and a portable parameterization of the intermolecular potential linking ab initio evaluation of intramolecular potentials and the partitioning of molecular polarizability; (2) the exploitation of an efficient coordinated porting and running of molecular dynamics codes on the European grid distributed computing infrastructure. As a prototype case study, the N-methylacetamide dimer in vacuo has been considered and the formation of possible conformers is analyzed.