Generalized correlations in terms of polarizability for van der Waals interaction potential parameter calculations

Four-Dimensional Scaling of Dipole Polarizability: From Single-Particle Models to Atoms and Molecules

Scaling laws enable the determination of physicochemical properties of molecules and materials as a function of their size, density, number of electrons or other easily accessible descriptors. Such relations can be counterintuitive and nonlinear, and ultimately yield much needed insight into quantum mechanics of many-particle systems. In this work, we show on the basis of single-particle models, multielectron atoms and molecules that the dipole polarizability of quantum systems is generally proportional to the fourth power of a characteristic length, computed from the ground-state wave function. This four-dimensional (4D) scaling is independent of the ratio of bound-to-bound and bound-to-continuum electronic transitions and applies to many-electron atoms when a correlated length metric is used. Finally, this scaling law is applied to predict the polarizability of molecules by electrostatically coupled atoms-in-molecules approach, obtaining approximately 8% absolute and relative accuracy with respect to hybrid density functional theory (DFT) on the QM7-X data set of organic molecules, providing an efficient and scalable model for the anisotropic polarizability tensors of extended (bio)molecules.

Formation of Products of Incomplete Destruction (PID) from the Thermal Oxidative Decomposition of Perfluorooctanoic Acid (PFOA): Measurement, Modeling, and Reaction Pathways

The thermal decomposition of perfluorooctanoic acid (PFOA) under oxidative conditions was investigated using air (O2) and N2O as oxidants over temperatures ranging from 400 to 1000 °C in an α-alumina reactor. In the presence of air, PFOA was found to decompose into perfluorohept-1-ene (C7F14) and perfluoroheptanoyl fluoride (C7F14O) in addition to HF, CO, and CO2. At temperatures above 800 °C, both C7F14 and C7F14O were no longer detected. A comprehensive analysis of the reaction mechanisms through quantum chemical analysis and kinetic modeling in combination with experimental observations was utilized to identify key reaction pathways. Quantum chemical analysis led to the conclusion that oxygen atoms are crucial in decomposing perfluoroalk-1-enes, especially the stable perfluorohept-1-ene (C7F14). Under oxidative conditions, it was found that significant quantities of C2F6 and CF4 were formed. Further quantum chemical analysis suggests that the O atoms facilitate the formation of volatile fluorinated compounds (VFCs) such as tetrafluoromethane (CF4) and hexafluoroethane (C2F6), particularly at higher temperatures. By elucidating these key reactions, an improved understanding of the potential formation products of incomplete combustion (PICs) or products of incomplete destruction (PIDs) is made.

The dawn of hydrogen and halogen bonds and their crucial role in collisional processes probing long-range intermolecular interactions

This perspective review focuses on the results of an internally consistent study developed in the Perugia laboratory, centered on the fundamental interaction components that, at large intermolecular distances, determine the formation of weak intermolecular hydrogen (HB) and halogen (XB) bonds. This investigation exploits old and novel molecular beam scattering experiments involving several gaseous prototypical systems. In particular, we focus on the kinetic energy dependence of the total (elastic + inelastic) integral cross-sections. Of particular interest is the measure of quantum interference patterns in the energy dependence of cross-sections of targeted systems and their shift compared to that of known reference systems. We interpreted these findings as interaction energy stabilization components, such as charge transfer, σ-hole, and polar flattening, that emerge at intermediate separation distance ranges and selectively manifest for specific geometries of collision complexes. Another significant observable we discuss is the absolute value of the cross-section and its dependence on permanent multipole moments of the collisional partners. Specifically, we show how the spontaneous orientation of rotationally cold and polar molecules, due to the electric field gradient associated with the interaction between permanent multipole moments, can significantly modify the magnitude of the total cross-section, even at high values of the impact parameter. We are confident that the present results can help extend the force field formulation to various interacting systems and carry out molecular dynamics simulations under conditions of application interest.

The Sticking of N2 on W(100) Surface: An Improvement in the Description of the Adsorption Dynamics Further Reconciling Theory and Experiment

The adsorption of nitrogen molecules on a (100) tungsten surface has been studied using a new potential energy surface in which long-range interactions are suitably characterized and represented by the Improved Lennard-Jones function. The new potential energy surface is used to carry out molecular dynamics simulations by adopting a semiclassical collisional method that explicitly includes the interaction with the surface phonons. The results of the sticking probability, evaluated as a function of the collision energy, are in good agreement with those obtained in the experiments and improve the already good comparison recently obtained with calculations performed using interactions from the Density Functional Theory method and corrected for long-range van der Waals contributions. The dependence of trapping probability on the surface temperature for a well-defined collision energy has also been investigated.

Electronic and dynamical properties of non-covalent diatomic aggregates formed by He with neutral and ionic Li and Be

CONTEXT AND RESULTS

This work aims to study the influence of the absence and presence of permanent charges on the electronic and dynamical properties of the non-covalent bound diatomic systems involving He and Li, Be as neutral and ionic partners. The charge displacement results suggest that in the formation of HeLi[Formula: see text], HeBe[Formula: see text], and HeBe[Formula: see text], the neutral He atom undergoes, in the electric field of the ion, a pronounced electronic polarization, and the natural bond order theoretical approach indicates that in the formation of the molecular orbital He acts as a weak electron donor. The energy decomposition analysis provides the dispersion and induction components as the attractive leading terms controlling the stability of all systems, confirming that the formed bond substantially maintains a non-covalent nature which is also supported by the Quantum Theory of Atoms in Molecules (QTAIM) analysis. Finally, it was found that the HeLi and HeBe neutral systems are unstable under any condition, HeLi[Formula: see text] and HeBe[Formula: see text] ionic systems are stable below 317K and 138K, respectively, while the HeBe[Formula: see text] system becomes unstable only after 3045K.

COMPUTATIONAL AND THEORETICAL TECHNIQUES

The potential energy curves and interactions in all systems were studied theoretically based on coupled-cluster singles and doubles method with perturbative inclusion of triples CCSD(T) method with an aug-cc-pV5Z basis set. More precisely, it was determined the potential energy curves describing the stability of the HeLi, HeLi[Formula: see text], HeBe, HeBe[Formula: see text], and HeBe[Formula: see text] systems, the charge displacement within the formed adducts, the decomposition of their total interaction energy, the topological analysis of their bonds, their rovibrational energies, their spectroscopic constants and lifetimes.

A Detailed Study of Electronic and Dynamic Properties of Noble Gas-Oxygen Molecule Adducts

In this work, the binding features of adducts formed by a noble gas (Ng = He, Ne, Ar, Kr, Xe, and Rn) atom and the oxygen molecule (O2) in its ground Σg-3, in the past target of several experimental studies, have been characterized under different theoretical points of view to clarify fundamental aspects of the intermolecular bond. For the most stable configuration of all Ng-O2 systems, binding energy has been calculated at the theory's CCSD(T)/aug-cc-pVTZ level and compared with the experimental findings. Rovibrational energies, spectroscopic constants, and lifetime as a function of temperature were also evaluated by adopting properly formulated potential energy curves. The nature of the interaction involved was deeply investigated using charge displacement analysis, symmetry-adapted perturbation theory (SAPT), and natural bond orbital (NBO) methods. In all adducts, it was found that the charge transfer plays a minor role, although O2 is an open shell species exhibiting a positive electron affinity. Obtained results also indicate that the dispersion attraction contribution is the main responsible for the complex stability.

Molecular hydrogen isotope separation by a graphdiyne membrane: a quantum-mechanical study

Graphdiyne (GDY) has emerged as a very promising two-dimensional (2D) membrane for gas separation technologies. One of the most challenging goals is the separation of deuterium (D₂) and tritium (T₂) from a mixture with the most abundant hydrogen isotope, H₂, an achievement that would be of great value for a number of industrial and scientific applications. In this work we study the separation of hydrogen isotopes in their transport through a GDY membrane due to mass-dependent quantum effects that are enhanced by the confinement provided by its intrinsic sub-nanometric pores. A reliable improved Lennard-Jones force field, optimized on accurate ab initio calculations, has been built to describe the molecule-membrane interaction, where the molecule is treated as a pseudoatom. The quantum dynamics of the molecules impacting on the membrane along a complete set of incidence directions have been rigorously addressed by means of wave packet calculations in the 3D space, which have allowed us to obtain transmission probabilities and, in turn, permeances, as the thermal average of the molecular flux per unit pressure. The effect of the different incidence directions on the probabilities is analyzed in detail and it is concluded that restricting the simulations to a perpendicular incidence leads to reasonable results. Moreover, it is found that a simple 1D model-using a zero-point energy-corrected interaction potential-provides an excellent agreement with the 3D probailities for perpendicular incidence conditions. Finally, D₂/H₂ and T₂/H₂ selectivities are found to reach maximum values of about 6 and 21 at ≈50 and 45 K, respectively, a feature due to a balance between zero-point energy and tunneling effects in the transport dynamics. Permeances at these temperatures are below recommended values for practical applications, however, at slightly higher temperatures (77 K) they become acceptable while the selectivities preserve promising values, particularly for the separation of tritium.

Unsaturated Dinitriles Formation Routes in Extraterrestrial Environments: A Combined Experimental and Theoretical Investigation of the Reaction between Cyano Radicals and Cyanoethene (C2H3CN)

The reaction between cyano radicals (CN, X²Σ⁺) and cyanoethene (C₂H₃CN) has been investigated by a combined approach coupling crossed molecular beam (CMB) experiments with mass spectrometric detection and time-of-flight analysis at a collision energy of 44.6 kJ mol^–1 and electronic structure calculations to determine the relevant potential energy surface. The experimental results can be interpreted by assuming the occurrence of a dominant reaction pathway leading to the two but-2-enedinitrile (1,2-dicyanothene) isomers (E- and Z-NC–CH=CH–CN) in a H-displacement channel and, to a much minor extent, to 1,1-dicyanoethene, CH₂C(CN)₂. In order to derive the product branching ratios under the conditions of the CMB experiments and at colder temperatures, including those relevant to Titan and to cold interstellar clouds, we have carried out RRKM statistical calculations using the relevant potential energy surface of the investigated reaction. We have also estimated the rate coefficient at very low temperatures by employing a semiempirical method for the treatment of long-range interactions. The reaction has been found to be barrierless and fast also under the low temperature conditions of cold interstellar clouds and the atmosphere of Titan. Astrophysical implications and comparison with literature data are also presented. On the basis of the present work, 1,2-dicyanothene and 1,1-dicyanothene are excellent candidates for the search of dinitriles in the interstellar medium.

Potential Energy Surfaces for Noble Gas (Ar, Kr, Xe, Rn)–Propylene Oxide Systems: Analytical Formulation and Binding

Multidimensional potential energy surfaces for heavy noble gas–propylene oxide systems are obtained by applying the phenomenological method successfully used to describe homologous systems involving He and Ne atoms. Such potential energy surfaces, where the interaction exclusively arises from the anisotropic van der Waals interaction components, are given in an analytical form. Therefore, they can be easily used as force fields to carry out molecular simulations to evaluate spectroscopic features and the dynamical selectivity of weakly bound complexes formed by propylene oxide (a prototype chiral species) with a noble gas atom (a prototype isotropic partner) by two-body collisions under a variety of conditions. Several potential energy minima are identified on the surfaces, which are confirmed and characterized by high level ab initio calculations. The next step to further generalize this methodology is its extension to systems involving propylene oxide-diatomic molecules (as H2, O2 and N2), as well as to propylene oxide dimers.

Helium structures around SF5+ and SF6+: novel intermolecular potential and mass spectrometry experiments

Helium clusters around the recently experimentally observed sulphur hexafluoride SF₆⁺ and sulphur pentafluoride SF₅⁺ ions are investigated using a combined experimental and theoretical effort. Mass spectrometry ion yields are obtained and the energetics and structure of the corresponding He_N–SF₆⁺ and He_N–SF₅⁺ clusters are analyzed using path integral molecular dynamics calculations as a function of N, the number of He atoms, employing a new intermolecular potential describing the interaction between the dopant and the surrounding helium. The new force field is optimized on benchmark potential energy ab initio calculations and represented by improved Lennard-Jonnes analytical expressions. This procedure improves the previous potentials employed in similar simulations for neutral SF₆ attached to helium nanodroplets. The theoretical analysis explains the characteristic features observed in the experimental ion yields which suggest the existence of stable configurations at specific sizes.

The structure of the He atoms around SF₅⁺ and SF₆⁺ is investigated both experimentally and theoretically.

Master Equation Study of Hydrogen Abstraction from HCHO by OH Via a Chemically Activated Intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction...

A computational study of the HO2 + SO3 → HOSO2 + 3O2 reaction catalyzed by water monomer, water dimer and small clusters of sulfuric acid: kinetics and atmospheric implications

Herein, the reaction mechanisms and kinetics for the HO2 + SO3 → HOSO2 + 3O2 reaction catalyzed by water monomer, water dimer and small clusters of sulfuric acid have been...

Molecular beam scattering experiments on noble gas-propylene oxide: Total integral cross sections and potential energy surfaces of He- and Ne-C3H6O

The interactions of He and Ne with propylene oxide have been investigated with the molecular beam technique by measuring the total (elastic + inelastic) integral cross section as a function of collision velocity. Starting from the analysis of these experimental data, potential energy surfaces, formulated as a function of the separation distance and orientation of propylene oxide with respect to the interacting partners, have been built: The average depth of potential wells (located at intermediate separation distances) has been characterized by analyzing the observed "glory" quantum effects, and the strength of long-range attractions has been obtained from the magnitude and the velocity dependence of the smooth component of measured cross sections. The surfaces, tested and improved against new ab initio calculations of minima interaction energies at the complete basis set level of theory, are defined in the full space of relative configurations. This represents a crucial condition to provide force fields useful to carry out, in general, important molecular property simulations and to evaluate, in the present case, the spectroscopic features and the dynamical selectivity of weakly bound complexes formed by propylene oxide, a prototype chiral species, during collisions in interstellar clouds and winds, in the space and planetary atmospheres. The adopted formulation of the interaction can be readily extended to similar systems, involving heavier noble gases or diatomic molecules (H₂, O₂, and N₂) as well as to propylene oxide dimers.

Influence of the solvent on the Lewis acidity of antimony pentahalides and group 13 Lewis acids toward acetonitrile and pyridine

Energetic effects of solvation of SbF₅ , SbCl₅ , and 21 group 13 Lewis acids (LA) and their molecular complexes with acetonitrile and pyridine are evaluated using SMD approach. Compared to the gas phase, solvation increases the stability of boron- and aluminum-containing complexes but decreases the stability of gallium and indium-containing homologs due to larger solation energies of free LA. New Lewis acidity scales, based on the Gibbs energy of dissociation of the molecular complexes LA·pyridine and LA·acetonitrile in the gas phase, in benzene and acetonitrile solutions, are proposed.

Molecular Collision Diameters and Electronic Polarizabilities: Inherent Relationship and Fast Evaluation

The collision diameter σ for a large set of molecular species is related to the static electronic polarizability α_el. A remarkable correlation between these quantities conceptually similar to the analogous one previously identified for atoms is revealed. Our recommended model is the function σ(α_el) = p₁ + p₂α_el^1/3, where p₁ = 0.768 Å and p₂ = 2.168 are the fitting parameters providing the best overall match to the reference data for collision diameter (181 data points). The obtained correlation allows one to easily find the collision diameter of molecules from the known polarizability and vice versa. These findings can be useful for many applications, where there is a need for inexpensive assessments of the collision diameters or electronic polarizabilities, for example, when developing the transport property databases for modeling of chemically reacting flows.

Spectroscopy, lifetime, and charge-displacement of the methanol-noble gas complexes: An integrated experimental-theoretical investigation

An integrated experimental-theoretical investigation was employed to determine rovibrational energies, spectroscopic constants, lifetime as a function of temperature in gas phase complexes of methanol with noble gas (NgHe, Ne, Ar, Kr, Xe, and Rn). Beside that, a parallel effort has been addressed to theoretically characterize the nature of intermolecular interactions determining the dissociation energy and equilibrium distance of the formed adducts. Dynamics and lifetime results reveal that, except for the CH₃OH-He aggregate, all other methanol-Ng compounds are sufficiently stable under thermal conditions. Their lifetimes are larger than 1 ps for the temperature of the bulk in the range between 200 and 500 K. In addition, the current lifetime results suggest that the aggregates formed by methanol and Ng are globally more stable than corresponding complexes formed by water with Ng. From the point of view of the CCSD(T)/aug-cc-pVTZ level calculation, in all compounds, the electron densities of Ng partners are weakly polarized in the presence of CH₃OH molecule. The charge-displacement curves and NBO analysis indicate that the charge transfer from Ng to methanol molecule, in general, plays a minor role, being appreciable only in the aggregate involving Ar. Finally, it was verified from the SAPT2 + (CCD)-δMP2/aug-cc-pVTZ calculations and NCI analysis that the dispersion is the essential long-range attractive contribution to the interaction energy for all studied complexes. This feature strongly suggests that these compounds are held bonded substantially by van der Waals forces. Then non-covalent intermolecular bonds are effectively formed in the gas phase, which is disturbed by small stabilizing charge-transfer contributions.

On the formation of van der Waals complexes through three-body recombination

In this work, we show that van der Waals molecules X-RG (where RG is the rare gas atom) may be created through direct three-body recombination collisions, i.e., X + RG + RG → X-RG + RG. In particular, the three-body recombination rate at temperatures relevant for buffer gas cell experiments is calculated via a classical trajectory method in hyperspherical coordinates [Pérez-Ríos et al., J. Chem. Phys. 140, 044307 (2014)]. As a result, it is found that the formation of van der Waals molecules in buffer gas cells (1 K ≲ T ≲ 10 K) is dominated by the long-range tail (distances larger than the LeRoy radius) of the X-RG interaction. For higher temperatures, the short-range region of the potential becomes more significant. Moreover, we notice that the rate of formation of van der Walls molecules is of the same order of the magnitude independent of the chemical properties of X. As a consequence, almost any X-RG molecule may be created and observed in a buffer gas cell under proper conditions.

Nature and role of the weak intermolecular bond in enantiomeric conformations of H2O2–noble gas adducts: a chiral prototypical model

The role and nature of the weak intermolecular bond in the H₂O₂–noble gas enantiomeric conformations are presented. Charge transfer associated with the formation of a weak intermolecular hydrogen bond tends to stabilize the cis-barrier conformation.

Leading Interaction Components in the Structure and Reactivity of Noble Gases Compounds

The nature, strength, range and role of the bonds in adducts of noble gas atoms with both neutral and ionic partners have been investigated by exploiting a fine-tuned integrated phenomenological–theoretical approach. The identification of the leading interaction components in the noble gases adducts and their modeling allows the encompassing of the transitions from pure noncovalent to covalent bound aggregates and to rationalize the anomalous behavior (deviations from noncovalent type interaction) pointed out in peculiar cases. Selected adducts affected by a weak chemical bond, as those promoting the formation of the intermolecular halogen bond, are also properly rationalized. The behavior of noble gas atoms excited in their long-life metastable states, showing a strongly enhanced reactivity, has been also enclosed in the present investigation.

Pyrolysis of Glyphosate and Its Toxic Products

With health concerns developing about the use of glyphosate (phosphonomethylglycine, PMG), the world's most used herbicide, the possibility of destruction of stockpiles via incineration arises. Little is known, however, about the possible toxic products of decomposition. We have performed a quantum chemical computation of the mechanism of thermal decomposition of PMG. Two initiation channels, one producing sarcosine and the other producing N-methylaminomethylphosphonic acid, have been located. Both the initial products decompose to dimethylamine (DMA), and the mechanism of further decomposition and toxic products is explored. Global potential energy surfaces for the initial decomposition of DMA are presented together with chemical kinetic modeling wherein the rate constants employed have been calculated from the quantum chemical data. Time and temperature evolution of the expected toxic products are presented and discussed.

Structure Distribution of Gaseous Ions in Strong Electrostatic Fields

In drift tube experiments, where ions move in gases under the action of an electrostatic field, collision excitation is implemented for the study of the energy partitioning in the molecular degrees of freedom and the corresponding relaxation rates when the excitation is turned off. In the case of flexible ions, the vibration modes related to metastable molecular structures have been activated in ion mobility spectrometry and their population has been probed with respect to the field strength and the gas temperature. Here, we study the angular vibrational excitation and relaxation of such systems by examining the motion of molecular ions with one bending mode at strong fields using a nonequilibrium molecular dynamics simulation method. The relatively stable structures are introduced through the use of an intramolecular angular potential with minima at the position of the most stable conformations. We calculate the first few moments of the velocity and angular velocity distribution functions as well as the distribution of the conformers, and find that they follow unified curves when plotted with respect to the relative ion-atom collision energy. At high field strengths, the angular vibration is excited and a portion of the ions interchanges conformations continuously in time with the populations of the molecular structures to attain limiting values. In addition, orientational alignment, with the perpendicular angular momentum being greater than the one parallel to the field, is observed. Our observations, although based on a specific system, must be rather general for the case of large flexible molecular ions.

Molecular Beam Scattering Experiments as a Sensitive Probe of the Interaction in Bromine-Noble Gas Complexes

This paper reports for the first time molecular beam experiments for the scattering of He, Ne, and Ar by the Br₂ molecule, with the aim of probing in detail the intermolecular interaction. Measurements have been performed under the experimental condition to resolve the glory pattern, a quantum interference effect observable in the collision velocity dependence of the integral cross section. We analyzed the experimental data with a reliable potential model defined as a combination of an anisotropic van der Waals component with the additional contribution due to charge transfer and polar flattening effects related to the formation of an intermolecular halogen bond. The model involves few parameters, whose values are related to fundamental physical properties of the interacting partners, and it allows an internally consistent comparison of the stability of the gas-phase adducts formed by Br₂ moiety with different noble gases as well as homologous complexes with the Cl₂ molecule. The same model appears to be also easily generalized to describe the interaction of diatomic halogen molecules in the excited B(³Π) electronic state where the halogen bond contribution tends to vanish and more anisotropic van der Waals components dominate the structure of the complexes with noble gases.

Dynamics and spectroscopy of van der Waals complexes composed of ammonia and noble gases

In this work, we calculate the rovibrational energies and spectroscopic constants for the systems formed by ammonia (NH₃) and noble gases (Ng=He, Ne, Ar, Kr and Xe). For the spectroscopic constant calculations, we used two different methods: Dunham and another one that use rovibrational energies (here calculated by discrete variable method). In both cases, we used the improved Lennard-Jones potential energy curves (PECs). These PECs, which describe very well van der Waals systems, were built using the dissociation and equilibrium distance obtained from experiments of crossed molecular beams. The spectroscopic constant results, obtained by both methods were in excellent agreement with each other for all NH₃-Ng studied systems. Also in relation to NH₃-He system, we realize that although this system has a relatively small dissociation energy, it has one vibrational level. Finally, the spectroscopic constants and fundamental rovibrational energy results were used to verify the stability of each system through the lifetime decomposition.

Modelling Charge Transfer in Weak Chemical Bonds: Insights from the Chemistry of Helium

We studied the nature of the interaction of the weakly bound Be-He adduct by means of an integrated theoretical approach based on high-level quantum chemical calculations for the characterization of the potential energy surfaces and charge displaced upon adduct formation, together with the development of a semi-empirical analytical formulation of the interaction potential. Our results show that Be is able to form a stable adduct with He when the Be(¹ D) (1s² 2s² →1s² 2s⁰ 2p² ) excited state is involved, with a binding energy of as much as 10.2 kcal/mol, an astonishingly large value for He in neutral systems. The analysis of the leading interaction components in the Be*-He adduct proves the relevance of the charge transfer to the overall stability, which contributes to decreasing the intermolecular distance, thus strengthening the induction-energy component.

Ab initio kinetics of the HOSO2 + 3O2 → SO3 + HO2 reaction

The detailed kinetic mechanism of the HOSO₂ + ³O₂ reaction, which plays a pivotal role in the atmospheric oxidation of SO₂, was investigated using accurate electronic structure calculations and novel statistical thermodynamic/kinetic models. Explored using the accurate composite method W1U, the detailed potential energy surface (PES) revealed that the addition of O₂ to a HOSO₂ radical to form the adduct (HOSO₄) proceeds via a transition state with a slightly positive barrier (i.e., 0.7 kcal mol^-1 at 0 K). Such a finding compromises a long-term hypothesis about this channel of being a barrierless process. Moreover, the overall reaction was found to be slightly exothermic by 1.7 kcal mol^-1 at 0 K, which is in good agreement with recent studies. On the newly-constructed PES, the temperature- and pressure-dependent behaviors of the title reaction were characterized in a wide range of conditions (T = 200-1000 K & P = 10-760 Torr) using the integrated deterministic and stochastic master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) rate model in which corrections for hindered internal rotation (HIR) and tunneling treatments were included. The calculated numbers were found to be in excellent agreement with literature data. The sensitivity analyses on the derived rate coefficients with respect to the ab initio input parameters (i.e., barrier height and energy transfer) were also performed to further understand the kinetic behaviors of the title reaction. The detailed kinetic mechanism, consisting of thermodynamic and kinetic data (in NASA polynomial and modified Arrhenius formats, respectively), was also provided at different T & P for further use in the modeling/simulation of any related systems.

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.

Modeling of Manganese Atom and Dimer Isolated in Solid Rare Gases: Structure, Stability, and Effect on Spin Coupling

Structures and energies of the trapping sites of manganese atom and dimer in solid Ar, Kr, and Xe are investigated within the classical model, which balances local distortion and long-range crystal order of the host and provides a means to estimate the relative site stabilities. The model is implemented with the additive pairwise potential field based on the ab initio and best empirical interatomic potential functions. In agreement with experiment, Mn single substitution (SS) and tetrahedral vacancy (TV) occupation are identified as stable for Ar and Kr, whereas the SS site is only found for Xe. Stable trapping sites of the weakly bound Mn₂ dimer are shown to be the mergers of SS and/or TV atomic sites. For Ar, (SS + SS) and (TV + TV) sites are close in energy, whereas (SS + TV) site lies higher. The (SS + SS) accommodation is identified as the only stable site in Kr and Xe at low energies. The results are compared with the resonance Raman, electron spin resonance, and absorption spectroscopy data. Reproducing the numbers of stable sites, the calculations tend to underestimate the matrix effect on the dimer vibrational frequency and spin-spin coupling constant. Nonetheless, the level of agreement is found to be informative for tentative assignments of the complex features seen in Mn₂ matrix isolation spectroscopy.

Methanol–methanol and methanol–water systems: the intermolecular interactions controlling the transition from small clusters to the liquid phase

The present paper focuses on the characterization of the properties of methanol and water molecules in gas and liquid enviroments.

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.

Interaction of O2 with CH4, CF4, and CCl4 by Molecular Beam Scattering Experiments and Theoretical Calculations

Gas phase collisions of O2 by CH4, CF4, and CCl4 have been investigated with the molecular beam technique by measuring both the integral cross section value, Q, and its dependence on the collision velocity, v. The adopted experimental conditions have been appropriate to resolve the oscillating "glory" pattern, a quantum interference effect controlled by the features of the intermolecular interaction, for all the three case studies. The analysis of the Q(v) data, performed by adopting a suitable representation of the intermolecular potential function, provided the basic features of the anisotropic potential energy surfaces at intermediate and large separation distances and information on the relative role of the physically relevant types of contributions to the global interaction. The present work demonstrates that while O2-CH4 and O2-CF4 are basically bound through the balance between size (Pauli) repulsion and dispersion attraction, an appreaciable intermolecular bond stabilization by charge transfer is operative in O2-CCl4. Ab initio calculations of the strength of the interaction, coupled with detailed analysis of electronic charge displacement promoted by the formation of the dimer, fully rationalizes the experimental findings. This investigation indicates that the interactions of O2, when averaged over its relative orientations, are similar to that of a noble gas (Ng), specifically Ar. We also show that the binding energy in the basic configurations of the prototypical Ng-CF4,CCl4 systems [ Cappelletti , D. ; Chem. Eur. J. 2015 , 21 , 6234 - 6240 ] can be reconstructed by using the interactions in Ng-F and Ng-Cl systems, previously characterized by molecular beam scattering experiments of state-selected halogen atom beams. This information is fundamental to approach the modeling of the weak intermolecular halogen bond. On the basis of the electronic polarizability, this also confirms [ Aquilanti , V. ; Angew. Chem., Int. Ed. 2005 , 44 , 2356 - 2360 ] that O2 can be taken as a proper reference partner for simulating the behavior of some basic noncovalent components of the interactions involving water. Present results are of fundamental relevance to build up the force field controlling the hydrophobic behavior of prototypical apolar CX4 (X = H, F, Cl) molecules.

Chirality of weakly bound complexes: the potential energy surfaces for the hydrogen-peroxide-noble-gas interactions

We consider the analytical representation of the potential energy surfaces of relevance for the intermolecular dynamics of weakly bound complexes of chiral molecules. In this paper we study the H2O2-Ng (Ng=He, Ne, Ar, Kr, and Xe) systems providing the radial and the angular dependence of the potential energy surface on the relative position of the Ng atom. We accomplish this by introducing an analytical representation which is able to fit the ab initio energies of these complexes in a wide range of geometries. Our analysis sheds light on the role that the enantiomeric forms and the symmetry of the H2O2 molecule play on the resulting barriers and equilibrium geometries. The proposed theoretical framework is useful to study the dynamics of the H2O2 molecule, or other systems involving O-O and S-S bonds, interacting by non-covalent forces with atoms or molecules and to understand how the relative orientation of the O-H bonds changes along collisional events that may lead to a hydrogen bond formation or even to selectivity in chemical reactions.

H2O–CH4and H2S–CH4complexes: a direct comparison through molecular beam experiments and ab initio calculations

Electron density redistribution upon the formation of the water–methane complex arises from polarisation and charge transfer effects.

The stereodynamics of the Penning ionization of water by metastable neon atoms

The stereodynamics of the Penning ionization of water molecules by collision with metastable neon atoms, occurring in the thermal energy range, is of great relevance for the understanding of fundamental aspects of the physical chemistry of water. This process has been studied by analyzing the energy spectrum of the emitted electrons previously obtained in our laboratory in a crossed beam experiment [B. G. Brunetti, P. Candori, D. Cappelletti, S. Falcinelli, F. Pirani, D. Stranges, and F. Vecchiocattivi, Chem. Phys. Lett. 539-540, 19 (2012)]. For the spectrum analysis, a novel semiclassical method is proposed, that assumes ionization events as mostly occurring in the vicinities of the collision turning points. The potential energy driving the system in the relevant configurations of the entrance and exit channels, used in the spectrum simulation, has been formulated by the use of a semiempirical method. The analysis puts clearly in evidence how different approaches of the metastable atom to the water molecule lead to ions in different electronic states. In particular, it provides the angular acceptance cones where the selectivity of the process leading to the specific formation of each one of the two energetically possible ionic product states of H2O(+) emerges. It is shown how the ground state ion is formed when neon metastable atoms approach water mainly perpendicularly to the molecular plane, while the first excited electronic state is formed when the approach occurs preferentially along the C2v axis, on the oxygen side. An explanation is proposed for the observed vibrational excitation of the product ions.