The N2Ar potential energy surface

Inelastic scattering in isotopologues of O2-Ar: the effects of mass, symmetry, and density of states

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.

Energy exchange rate coefficients from vibrational inelastic O2(Σg-3) + O2(Σg-3) collisions on a new spin-averaged potential energy surface

A new spin-averaged potential energy surface (PES) for non-reactive O₂(Σg-3) + O₂(Σg-3) collisions is presented. The potential is formulated analytically according to the nature of the principal interaction components, with the main van der Waals contribution described through the improved Lennard-Jones model. All the parameters involved in the formulation, having a physical meaning, have been modulated in restricted variation ranges, exploiting a combined analysis of experimental and ab initio reference data. The new PES is shown to be able to reproduce a wealth of different physical properties, ranging from the second virial coefficients to transport properties (shear viscosity and thermal conductivity) and rate coefficients for inelastic scattering collisions. Rate coefficients for the vibrational inelastic processes of O₂, including both vibration-to-vibration (V-V) and vibration-to-translation/rotation (V-T/R) energy exchanges, were then calculated on this PES using a mixed quantum-classical method. The effective formulation of the potential and its combination with an efficient, yet accurate, nuclear dynamics treatment allowed for the determination of a large database of V-V and V-T/R energy transfer rate coefficients in a wide temperature range.

Inelastic rate coefficients based on an improved potential energy surface for N2 + N2 collisions in a wide temperature range

Vibration-to-translation and vibration-to-vibration rate coefficients for N₂–N₂ inelastic scattering are calculated on an improved potential including high temperature regimes.

Insight into the halogen-bond nature of noble gas-chlorine systems by molecular beam scattering experiments, ab initio calculations and charge displacement analysis

We have carried out molecular-beam scattering experiments and high-level ab initio investigations on the potential energy surfaces of a series of noble-gas-Cl2 adducts. This effort has permitted the construction of a simple, reliable and easily generalizable analytical model potential formulation, which is based on a few physically meaningful parameters of the interacting partners and transparently shows the origin, strength, and stereospecificity of the various interaction components. The results demonstrate quantitatively beyond doubt that the interaction between a noble-gas (Ng) atom - even He - and Cl2 in a collinear configuration is characterized by weak halogen bond (XB) formation, accompanied by charge transfer (CT) from the Ng to chlorine. This characteristic, which stabilizes the adduct, rapidly disappears on going towards the T-shaped configuration, dominated by pure van der Waals (vdW) forces. Similarly, a pure vdW interaction takes place - with no CT component in any configuration - if Cl2 is present in the lowest πg* → σu* excited state, because the change in electron density that accompanies the excitation eliminates the Cl2 polar flattening and σ hole, making the XB interaction inaccessible.

Selective Emergence of the Halogen Bond in Ground and Excited States of Noble-Gas-Chlorine Systems

Molecular-beam scattering experiments and theoretical calculations prove the nature, strength, and selectivity of the halogen bonds (XB) in the interaction of halogen molecules with the series of noble gas (Ng) atoms. The XB, accompanied by charge transfer from the Ng to the halogen, is shown to take place in, and measurably stabilize, the collinear conformation of the adducts, which thus becomes (in contrast to what happens for other Ng-molecule systems) approximately as bound as the T-shaped form. It is also shown how and why XB is inhibited when the halogen molecule is in the ³ Π_u excited state. A general potential formulation fitting the experimental observables, based on few physically essential parameters, is proposed to describe the interaction accurately and is validated by ab initio computations.

Closer versus Long Range Interaction Effects on the Non-Arrhenius Behavior of Quasi-Resonant O2 + N2 Collisions

We report in this paper an investigation on energy transfer processes from vibration to vibration and/or translation in thermal and subthermal regimes for the O₂ + N₂ system performed using quantum-classical calculations on different empirical, semiempirical, and ab initio potential energy surfaces. In particular, the paper focuses on the rationalization of the non-Arrhenius behavior (inversion of the temperature dependence) of the quasi-resonant vibration-to-vibration energy transfer transition rate coefficients at threshold. To better understand the microscopic nature of the involved processes, we pushed the calculations to the detail of the related cross sections and analyzed the impact of the medium and long-range components of the interaction on them. Furthermore, the variation with temperature of the dependence of the quasi-resonant rate coefficient on the vibrational energy gap between initial and final vibrational states and the effectiveness of quantum-classical calculations to overcome the limitations of the purely classical treatments were also investigated. These treatments, handled in an open molecular science fashion by chaining data and competencies of the various laboratories using a grid empowered molecular simulator, have allowed a rationalization of the dependence of the computed rate coefficients in terms of the distortion of the O₂-N₂ configuration during the diatom-diatom collisions. A way of relating such distortions to a smooth and continuous progress variable, allowing a proper evolution from both long to closer range formulation of the interaction and from its entrance to exit channel (through the strong interaction region) relaxed graphical representations, is also discussed in the paper.

Enhanced Flexibility of the O2 + N2 Interaction and Its Effect on Collisional Vibrational Energy Exchange

Prompted by a comparison of measured and computed rate coefficients of Vibration-to-Vibration and Vibration-to-Translation energy transfer in O2 + N2 non-reactive collisions, extended semiclassical calculations of the related cross sections were performed to rationalize the role played by attractive and repulsive components of the interaction on two different potential energy surfaces. By exploiting the distributed concurrent scheme of the Grid Empowered Molecular Simulator we extended the computational work to quasiclassical techniques, investigated in this way more in detail the underlying microscopic mechanisms, singled out the interaction components facilitating the energy transfer, improved the formulation of the potential, and performed additional calculations that confirmed the effectiveness of the improvement introduced.

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.

Efficiency of Collisional O2 + N2 Vibrational Energy Exchange

By following the scheme of the Grid Empowered Molecular Simulator (GEMS), a new O2 + N2 intermolecular potential, built on ab initio calculations and experimental (scattering and second virial coefficient) data, has been coupled with an appropriate intramolecular one. On the resulting potential energy surface detailed rate coefficients for collision induced vibrational energy exchanges have been computed using a semiclassical method. A cross comparison of the computed rate coefficients with the outcomes of previous semiclassical calculations and kinetic experiments has provided a foundation for characterizing the main features of the vibrational energy transfer processes of the title system as well as a critical reading of the trajectory outcomes and kinetic data. On the implemented procedures massive trajectory runs for the proper interval of initial conditions have singled out structures of the vibrational distributions useful to formulate scaling relationships for complex molecular simulations.

New exchange-Coulomb N2-Ar potential-energy surface and its comparison with other recent N2-Ar potential-energy surfaces

The reliability of five N2-Ar potential-energy surfaces in representing the N2-Ar interaction has been investigated by comparing their abilities to reproduce a variety of experimental results, including interaction second viral coefficients, bulk transport properties, relaxation phenomena, differential scattering cross sections, and the microwave and infrared spectra of the van der Waals complexes. Four of the surfaces are the result of high-level ab initio quantal calculations; one of them utilized fine tuning by fitting to microwave data. To date, these four potential-energy surfaces have only been tested against experimental microwave data. The fifth potential-energy surface, based upon the exchange-Coulomb potential-energy model for the interaction of closed-shell species, is developed herein: it is a combination of a damped dispersion energy series and ab initio calculations of the Heitler-London interaction energy, and has adjustable parameters determined by requiring essentially simultaneous agreement with selected quality interaction second viral coefficient and microwave data. Comparisons are also made with the predictions of three other very good literature potential-energy surfaces, including the precursor of the new exchange-Coulomb potential-energy surface developed here. Based upon an analysis of a large body of information, the new exchange-Coulomb and microwave-tuned ab initio potential-energy surfaces provide the best representations of the N2-Ar interaction; nevertheless, the other potential-energy surfaces examined still have considerable merit with respect to the prediction of specific properties of the N2-Ar van der Waals complex. Of the two recommended surfaces, the new exchange-Coulomb surface is preferred on balance due to its superior predictions of the effective cross sections related to various relaxation phenomena, and to its reliable, and relatively simple, representation of the long-range part of the potential-energy surface. Moreover, the flexibility still inherent in the exchange-Coulomb potential form can be further exploited, if required, in future studies of the N2-Ar interaction.

Highly excited rovibrational states of small molecules

The use of first principles variational calculations for the calculation of high-lying energy levels, wavefunctions and transition intensities for triatomic molecules is considered. Theoretical developments are considered, including the use of generalized internal coordinates, the use of a two-step procedure for rotationally excited systems and a finite element method known as the discrete variable representation. Illustrative calculations are presented including ones for H 2, LiCN and the Ar-N2 Van der Waals molecule. A first principles ‘rotational’ spectrum of H 2D+ is computed using states up to J = 30. The transition intensities in this spectrum are reproduced accurately in a frozen dipole approximation but are poorly represented by models that involve approximating the wavefunction.