Li + HF and Li + HCl Reactions Revisited I: QCT Calculations and Simulation of Experimental Results

The Li + HF and Li + HCl reactions share some common features. They have the same kinematics, relatively small barrier heights, bent transition states, and are both exothermic when the zero point energy is considered. Nevertheless, the pioneering crossed beam experiments by Lee and co-workers in the 80s (Becker et al., J. Chem. Phys. 1980, 73, 2833) revealed that the dynamics of the two reactions differ significantly, especially at low collision energies. In this work, we present theoretical simulations of their results in the laboratory frame (LAB), based on quasiclassical trajectories and obtained using accurate potential energy surfaces. The calculated LAB angular distributions and time-of-flight spectra agree well with the raw experimental data, although our simulations do not reproduce the experimentally derived center-of-mass (CM) differential cross section and velocity distributions. The latter were derived by forward convolution fitting under the questionable assumption that the CM recoil velocity and scattering angle distribution were uncoupled, while our results show that the coupling between them is relevant. Some important insights into the reaction mechanism discussed in the article by Becker et al. had not been contrasted with those that can be extracted from the theoretical results. Among them, the correlation between the angular momenta involved in the reactions has also been examined. Given the kinematics of both systems, the reagent orbital angular momentum, l , is almost completely transformed into the rotation of the product diatom, j'. However, contrary to the coplanar mechanism proposed in the original paper, we find that the initial and final relative orbital angular momenta are not necessarily parallel. Both reactions are found to be essentially direct, although about 15% of the LiFH complexes live longer than 200 fs.

OH(2Π) + C2H4 Reaction: A Combined Crossed Molecular Beam and Theoretical Study

The reaction between the ground-state hydroxyl radical, OH(²Π), and ethylene, C₂H₄, has been investigated under single-collision conditions by the crossed molecular beam scattering technique with mass-spectrometric detection and time-of-flight analysis at the collision energy of 50.4 kJ/mol. Electronic structure calculations of the underlying potential energy surface (PES) and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of product branching fractions on the derived PES for the addition pathway have been performed. The theoretical results indicate a temperature-dependent competition between the anti-/syn-CH₂CHOH (vinyl alcohol) + H, CH₃CHO (acetaldehyde) + H, and H₂CO (formaldehyde) + CH₃ product channels. The yield of the H-abstraction channel could not be quantified with the employed methods. The RRKM results predict that under our experimental conditions, the anti- and syn-CH₂CHOH + H product channels account for 38% (in similar amounts) of the addition mechanism yield, the H₂CO + CH₃ channel for ∼58%, while the CH₃CHO + H channel is formed in negligible amount (<4%). The implications for combustion and astrochemical environments are discussed.

Universal crossed beam imaging studies of polyatomic reaction dynamics

Crossed-beam imaging studies of polyatomic reactions show surprising dynamics not anticipated by extrapolation from smaller model systems.

Combined Experimental-Theoretical Study of the OH + CO → H + CO2 Reaction Dynamics

A combined experimental-theoretical study is performed to advance our understanding of the dynamics of the prototypical tetra-atom, complex-forming reaction OH + CO → H + CO₂, which is also of great practical relevance in combustion, Earth's atmosphere, and, potentially, Mars's atmosphere and interstellar chemistry. New crossed molecular beam experiments with mass spectrometric detection are analyzed together with the results from previous experiments and compared with quasi-classical trajectory (QCT) calculations on a new, full-dimensional potential energy surface (PES). Comparisons between experiment and theory are carried out both in the center-of-mass and laboratory frames. Good agreement is found between experiment and theory, both for product angular and translational energy distributions, leading to the conclusion that the new PES is the most accurate at present in elucidating the dynamics of this fundamental reaction. Yet, small deviations between experiment and theory remain and are presumably attributable to the QCT treatment of the scattering dynamics.

The Rovibrational Spectra of trans- and cis-HOCO, Calculated by MULTIMODE with ab Initio Potential Energy and Dipole Moment Surfaces

The code MULTIMODE is used in its reaction path version, along with ab initio potential energy and dipole moment surfaces introduced earlier, to predict the infrared spectra of both trans and cis forms of HOCO at temperatures 296 and 15 K. All six fundamentals are isolated for each isomer and temperature, and their main features examined, paying particular attention to the OH stretch fundamental, whose spectrum has been reported experimentally for trans-HOCO. The current spectra for cis-HOCO, while not of "spectroscopic" accuracy, should be sufficient to aid in new experimental efforts to record the spectrum of this isomer.

Crossed beam polyatomic reaction dynamics: recent advances and new insights

This review summarizes the developments in polyatomic reaction dynamics, focusing on reactions of unsaturated hydrocarbons with O-atoms and methane with atoms/radicals.

A Combined Experimental and Theoretical Study on the Formation of the 2-Methyl-1-silacycloprop-2-enylidene Molecule via the Crossed Beam Reactions of the Silylidyne Radical (SiH; X(2)Π) with Methylacetylene (CH3CCH; X(1)A1) and D4-Methylacetylene (CD3CCD; X(1)A1)

The bimolecular gas-phase reactions of the ground-state silylidyne radical (SiH; X(2)Π) with methylacetylene (CH3CCH; X(1)A1) and D4-methylacetylene (CD3CCD; X(1)A1) were explored at collision energies of 30 kJ mol(-1) under single-collision conditions exploiting the crossed molecular beam technique and complemented by electronic structure calculations. These studies reveal that the reactions follow indirect scattering dynamics, have no entrance barriers, and are initiated by the addition of the silylidyne radical to the carbon-carbon triple bond of the methylacetylene molecule either to one carbon atom (C1; [i1]/[i2]) or to both carbon atoms concurrently (C1-C2; [i3]). The collision complexes [i1]/[i2] eventually isomerize via ring-closure to the c-SiC3H5 doublet radical intermediate [i3], which is identified as the decomposing reaction intermediate. The hydrogen atom is emitted almost perpendicularly to the rotational plane of the fragmenting complex resulting in a sideways scattering dynamics with the reaction being overall exoergic by -12 ± 11 kJ mol(-1) (experimental) and -1 ± 3 kJ mol(-1) (computational) to form the cyclic 2-methyl-1-silacycloprop-2-enylidene molecule (c-SiC3H4; p1). In line with computational data, experiments of silylidyne with D4-methylacetylene (CD3CCD; X(1)A1) depict that the hydrogen is emitted solely from the silylidyne moiety but not from methylacetylene. The dynamics are compared to those of the related D1-silylidyne (SiD; X(2)Π)-acetylene (HCCH; X(1)Σg(+)) reaction studied previously in our group, and from there, we discovered that the methyl group acts primarily as a spectator in the title reaction. The formation of 2-methyl-1-silacycloprop-2-enylidene under single-collision conditions via a bimolecular gas-phase reaction augments our knowledge of the hitherto poorly understood silylidyne (SiH; X(2)Π) radical reactions with small hydrocarbon molecules leading to the synthesis of organosilicon molecules in cold molecular clouds and in carbon-rich circumstellar envelopes.

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.

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.

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.

A high-level ab initio study of the N2 + N2 reaction channel

A new six-dimensional (6D) global potential energy surface (PES) is proposed for the full range description of the interaction of the N2(1Σg+)+N2(1Σg+) system governing collisional processes, including N atom exchange. The related potential energy values were determined using high-level ab initio methods. The calculations were performed at a coupled-cluster with single and double and perturbative triple excitations level of theory in order to have a first full range picture of the PES. Subsequently, in order to accurately describe the stretching of the bonds of the two interacting N2 molecules by releasing the constraints of being considered as rigid rotors, for the same molecular geometries higher level of theory multi reference calculations were performed. Out of the calculated values a 6D 4-atoms global PES was produced for use in dynamical calculations. The ab initio calculations were made possible by the combined use of High Throughput Computing and High Performance Computing techniques within the frame of a computing grid empowered molecular simulator.