Scattering of aligned molecules. The potential energy surfaces for the Kr-O2 and Xe-O2 systems

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.

New Aspects of the Airglow Problem and Reactivity of the Dioxygen Quintet O2(5Πg) State in the MLT Region as Predicted by DFT Calculations

Dioxygen in the quintet O₂(⁵Π_g) state is a weakly bound species near the entrance of the O(³P) + O(³P) recombination channel. It was predicted by ab initio calculations in 1977 and detected experimentally in 1999. Meantime, the O₂(⁵Π_g) species was tentatively assumed as intermediate in transport properties calculations for the rarefied gases of the Earth's upper atmosphere, though its potential energy curve is still debated. Besides six other strongly bound low-lying states of dioxygen, the O₂(⁵Π_g) state is an important potential candidate for modeling energy transfer and airglow of the upper atmosphere. A number of photochemical kinetic schemes designed to simulate energy flow upon atomic and molecular oxygen collisions in the rarefied mesosphere take into account a participation of the O₂(⁵Π_g) state in energy relaxation processes responsible for terrestrial nightglow. All mechanisms of energy redistribution are based on the hard-sphere collision models. The possibility of chemical interactions between the quintet excited state of dioxygen and other atmospheric components has not been considered so far in photochemistry of the upper atmosphere. In the present paper, the chemical reactivity of the quintet O₂(⁵Π_g) species is calculated for the first time in the framework of the density functional theory. Definitely, O₂(⁵Π_g) is the most reactive species among all other metastable dioxygen states below 5.1 eV. Quintet products of the O₂(⁵Π_g) state association with heavy inert gases, H₂O, N₂, and CO₂ are predicted to be chemically significant, while the complexes with abundant H₂ and He species are rather weak and not important even in the mesopause low-temperature region. The complex with N₂ molecule is unexpectedly stable with dissociation energy 4 kJ/mol, which can strongly influence the abundant termolecular association O + O + N₂ → O₂ + N₂ process. Reaction with meteoritic ablated Mg atom produces metastable ⁵A₁ excited state of MgO₂ being more strongly bound than the ground ³A₂ state of magnesium peroxide.

Stereodynamical Effects by Anisotropic Intermolecular Forces

Electric and magnetic field gradients, arising from sufficiently strong anisotropic intermolecular forces, tend to induce molecular polarization which can often modify substantially the results of molecular collisions, especially at low rotational temperatures and low collision energies. The knowledge of these phenomena, today still not fully understood, is of general relevance for the control of the stereo-dynamics of elementary chemical-physical processes, involving neutral and ionic species under a variety of conditions. This paper reports on results obtained by combining information from scattering, spectroscopic and reactivity experiments, within a collaboration between the research groups in Perugia and Trento. We addressed particular attention to the reactions of small atomic ions with polar neutrals for their relevance in several environments, including interstellar medium, planetary atmospheres, and laboratory plasmas. In the case of ion-molecule reactions, alignment/orientation is a general phenomenon due to the electric field generated by the charged particle. Such phenomenon originates critical stereo-dynamic effects that can either suppress (when the orientation drives the collision complex into non-reactive or less reactive configurations), or enhance the reactivity (when orientation confines reagents in the most appropriate configuration for reaction). The associated rate coefficients show the propensity to follow an Arrhenius and a non-Arrhenius behavior, respectively.

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.

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.

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.

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.

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.

Cluster-enhanced X–O2 photochemistry (X=CH3I, C3H6, C6H12, and Xe)

The effect of a local environment on the photodissociation of molecular oxygen is investigated in the van der Waals complex X-O(2) (X=CH(3)I, C(3)H(6), C(6)H(12), and Xe). A single laser operating at wavelengths around 226 nm is used for both photodissociation of the van der Waals complex and simultaneous detection of the O((3)P(J),J=2,1,0) atom photoproduct via (2+1) resonance enhanced multiphoton ionization. The kinetic energy distribution (KED) and angular anisotropy of the product O atom recoil in this dissociation are measured using the velocity map imaging technique configured for either full ("crush") or partial ("slice") detection of the three-dimensional O((3)P(J)) atom product Newton sphere. The measured KED and angular anisotropy reveal a distinct difference in the mechanism of O atom generation from an X-O(2) complex compared to a free O(2) molecule. The authors identify two one-photon excitation pathways, the relative importance of which depends on IPx, the ionization potential of the X partner. One pathway, observed for all complexes independent of IPx, involves a direct transition to the perturbed covalent state X-O(2)(A'(3)Delta(u)) with excitation localized on the O(2) subunit. The predominantly perpendicular character of this channel relative to the laser polarization detection, together with data on the structure of the complex, allows us to confirm that X partner induced admixing of an X(+)-O(2) (-) charge transfer (CT) state is the perturbing factor resulting in the well-known enhancement of photoabsorption within the Herzberg continuum of molecular oxygen. The second excitation pathway, observed for X-O(2) complexes with X=CH(3)I and C(3)H(6), involves direct excitation into the (3)(X(+)-O(2) (-)) CT state of the complex. The subsequent photodissociation of this CT state by the same laser pulse gives rise to the superoxide anion O(2) (-), which then photodissociates, providing fast (0.69 eV) O atoms with a parallel image pattern. Products from the photodissociation of singlet oxygen O(2)(b (1)Sigma(g) (+)) are also observed when the CH(3)I-O(2) complex was irradiated. Potential energy surfaces (PES) for the ground and relevant excited states of the X-O(2) complex have been constructed for CH(3)I-O(2) using the results of CASSCF calculations for the ground and CT states of the complex as well as literature data on PES of the subunits. These model potential energy surfaces allowed us to interpret all of the observed O((3)P(J)) atom production channels.

Intermolecular interaction potentials for the Ar–C2H2, Kr–C2H2, and Xe–C2H2 weakly bound complexes: Information from molecular beam scattering, pressure broadening coefficients, and rovibrational spectroscopy

Integral cross sections and pressure broadening coefficients have been measured for the acetylene-krypton complex, by molecular beam scattering and by high resolution IR spectroscopy, respectively. A new potential energy surface (PES) is proposed to describe structure and dynamical properties of this prototypical weakly bound complex. The PES has been parametrized exploiting a novel atom-bond pairwise additive scheme and has been fitted to the experimental data. A similar PES has been obtained for the acetylene-xenon system by a proper scaling of the interaction parameters of the krypton case, based on empirical considerations. These PESs together with that recently proposed by the same authors [J. Phys. Chem. 109, 8471 (2005)] for the acetylene-argon case have been employed for close coupling calculations of the pressure broadening cross sections and for a characterization of the rovibrational structure of the complexes.

Intermolecular interactions of H2S with rare gases from molecular beam scattering in the glory regime and from ab initio calculations

Integral cross sections for collisions of rotationally hot H2S molecules with rare gas atoms (Ne, Ar, and Kr) have been measured, in the collision energy range of 10-60 kJ mol(-1), using a molecular beam apparatus operating under high resolution both in angle and in velocity. A well resolved glory pattern has been measured which permitted the accurate characterization of the intermolecular potentials both at long range (in the attractive region) and at intermediate distances (in the well region). Considering the conditions used in the experiments, the obtained potentials must be considered very close to the spherical averages of the full intermolecular potential energy surfaces. Extensive ab initio calculations have also been carried out in parallel in order to characterize energy minima in the potential energy surfaces and energy barriers associated to the motion of the rare gas atoms around H2S. An assessment of the relative role of the various interaction components has been also attempted: the combined analysis of experimental and theoretical results suggests that H2S-rare gas aggregates are mainly bound by nearly isotropic noncovalent interactions of the van der Waals type.

Xe nuclear magnetic resonance line shapes in channels decorated with paramagnetic centers

To make predictions of the Xe NMR line shapes for Xe in channels decorated with paramagnetic centers, we consider a model system using the O(2) molecule as the paramagnetic center. The previously calculated quantum mechanical Xe@O(2) hyperfine tensor for various configurations of Xe in the presence of O(2) provides a model for the hyperfine response of Xe atom to the presence of a paramagnetic center. The averaging is carried out using the same grand canonical Monte Carlo methodology as for calculating NMR line shapes for Xe in diamagnetic channels, modified to include the effects of the hyperfine tensor response. We explore the temperature dependence of the Xe line shapes, the dependence on the concentration, and the symmetry of distribution of embedded paramagnetic centers, on the orientation of the paramagnetic center axis with respect to the channel axis, and on the radial distance of the paramagnetic center from the axis of the channel. We predict Xe line shape signatures of the presence and orientation of paramagnetic centers and deduce which tensor elements provide measures of concentration and radial distance of paramagnetic centers from the channel axis.

Collision-induced alignment of H2O+ drifting in helium

The collision-induced alignment of H(2)O(+) drifting in helium is studied with a molecular dynamics method that has been extended to treat nonlinear rigid ions. Rotational distribution functions and averaged quantities are presented in terms of the rho formalism [M. Thachuk, Phys. Rev. A 72, 032722 (2005)], and it is shown that this description gives a very good agreement with simulation results. In addition to velocity and angular velocity autocorrelation functions, a velocity-angular velocity cross correlation function is introduced. This cross correlation function provides insight into the dynamical nature of the alignment mechanism.

Collision-induced dissociation of transition metal-oxide ions: dynamics of VO+ collision with Xe

The collision-induced dissociation of VO(+) by Xe has been studied by the use of classical dynamics procedures on London-Eyring-Polanyi-Sato potential-energy surfaces in the collision energy range of 5.0-30 eV. The dissociation threshold behavior and the dependence of reaction cross sections on the collision energy closely follow the observed data with the threshold energy of 6.00 eV. The principal reaction pathway is VO(+) + Xe --> V(+)+ O + Xe and the minor pathway is VO(+) + Xe--> VXe(+) + O. At higher collision energies (E > 8.0 eV), the former reaction preferentially occurs near the O-V(+)...Xe collinear and perpendicular alignments, but the latter only occurs near the perpendicular alignment. At lower energies close to the threshold, the reactions are found to occur near the collinear configuration. No reaction occurs in the collinear alignment V(+)-O...Xe. The high and low energy-transfer efficiencies of the collinear alignments O-V(+)...Xe and V(+)-O...Xe are attributed to the effects of mass distribution. The activation of the VO(+) bond toward the dissociation threshold occurs through a translation-to-vibration energy transfer in a strong collision on a time scale of about 50 fs.

New insights on the stereodynamics of ethylene adsorption on an oxygen-precovered silver surface

The control of spatial orientation of molecules has a great influence on the stereodynamics of elementary processes occurring both in homogeneous and heterogeneous phases. Nonpolar molecules have so far escaped direct experimental investigations because of their poor sensitivity to several external constraints. Recently, it has been shown that the collisional alignment produced in supersonic expansions coupled with molecular-beam velocity selection can help solve such problems. Here we show that the sticking probability of ethylene, a nonpolar molecule prototypical of unsaturated hydrocarbons, on an O(2)-precovered Ag(001) surface is larger for molecules approaching in a helicopter-like motion than for those cartwheeling. A mechanism involving a weakly bound precursor state is suggested, with helicopter molecules having a lower chance of being scattered back into the gas phase than cartwheels when colliding with preadsorbed ethylene.

Collision Cross Sections, Pressure-Broadening Coefficients and Second Virial Coefficients for the Acetylene-Argon Complex: Experiments and Calculations on a New Potential Energy Surface

Integral cross sections and pressure-broadening coefficients have been measured by molecular beam scattering and by high-resolution infrared spectroscopy, respectively, for the acetylene-argon system. A new potential energy surface (PES) is proposed to describe structure and dynamical properties of this prototypical weakly bound complex. The PES has been parametrized exploiting a novel atom-bond pairwise additive scheme and has been fitted to the experimental data. Calculations of the scattering cross sections (both differential and integral), pressure-broadening, and second virial coefficients have been performed using both the present and also the most recent ab initio PES available in the literature. Analysis of the new experimental data indicates that the anisotropy of the interaction in the well region should be larger than that obtained in ab initio calculations. This is also in line with previous spectroscopic results.

Threshold collision-induced dissociation of diatomic molecules: A case study of the energetics and dynamics of O2− collisions with Ar and Xe

The energetics and dynamics of collision-induced dissociation of O2- with Ar and Xe targets are studied experimentally using guided ion-beam tandem mass spectrometry. The cross sections and the collision dynamics are modeled theoretically by classical trajectory calculations. Experimental apparent threshold energies are 2.1 and 1.1 eV in excess of the thermochemical O2- bond dissociation energy for argon and xenon, respectively. Classical trajectory calculations confirm the observed threshold behavior and the dependence of cross sections on the relative kinetic energy. Representative trajectories reveal that the bond dissociation takes place on a short time scale of about 50 fs in strong direct collisions. Collision-induced dissociation is found to be remarkably restricted to the perpendicular approach of ArXe to the molecular axis of O2-, while collinear collisions do not result in dissociation. The higher collisional energy-transfer efficiency of xenon compared with argon is attributed to both mass and polarizability effects.

Molecular-beam study of the water-helium system: Features of the isotropic component of the intermolecular interaction and a critical test for the available potential-energy surfaces

We report molecular-beam measurements of the total integral cross sections for the scattering of water molecules by helium atoms. A combined analysis of the new experimental data together with available differential cross section results has allowed an accurate determination of the isotropic component of the interaction potential for this prototypical system. The potential well shows a depth of 0.265 +/- 0.010 kJ/mol at a distance between He and the center of mass of the water molecule of 0.345 +/- 0.02 nm. An effective isotropic long-range attraction constant C(LR) = (6.3+/-0.3) x 10(-4) kJ mol(-1) nm(-6), including both dispersion and induction contributions, has also been determined. The most recent and accurate ab initio potential-energy surfaces have been tested against these new experimental results.