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.

Theoretical reaction kinetics astride the transition between moderate and deep tunneling regimes: the F + HD case.

For the reaction between F and HD, giving HF + D and DF + H, the rate constants, obtained from rigorous quantum scattering calculations at temperatures ranging from 350 K down to 100 K, show deviations from the Arrhenius behavior that have been interpreted in terms of tunneling of either H or D atoms through a potential energy barrier. The interval of temperature investigated extends from above to below a crossover value Tc, a transition temperature separating the moderate and deep quantum tunneling regimes. Below Tc, the rate of the H or D exchange reaction is controlled by the prevalence of tunneling over the thermal activation mechanism. In this temperature range, Bell's early treatment of quantum tunneling, based on a semiclassical approximation for the barrier permeability, provides a reliable tool to quantitatively account for the contribution of the tunneling effect. This treatment is here applied for extracting from rate constants properties of the effective tunneling path, such as the activation barrier height and width. We show that this is a way of parametrizing the dependence of the apparent activation energy on temperature useful for both calculated and experimental rate constants in an ample interval of temperature, from above to below Tc, relevant for modelization of astrophysical and in general very low-temperature environments.

Quantum stereodynamics for the two product channels of the F + HD reaction from the complete scattering matrix in the stereodirected representation.

For the two exit arrangements of the F + HD reaction, the full scattering matrix is obtained by exact quantum dynamics on an accurate potential energy surface. The S matrix is expressed in the stereodirected representation, for the first time, for all channels of a triatomic reaction. We analyze a collision energy where the dominant reaction mechanism is direct and a total angular momentum J = 0. It is found that the introduction of steric quantum numbers (correlated in the vector model to the angles measuring directions of approaching reactants and of separating products) provides a sharp description of stereodynamical features for both exit channels.

On the role of scattering resonances in the F+HD reaction dynamics.

We study scattering resonances in the F+HD-->HF+D reaction using a new method for direct evaluation of the lifetime Q-matrix [Aquilanti et al., J. Chem. Phys. 2005, 123, 054314]. We show that most of the resonances are due to van der Waals states in the entrance and exit reaction channels. The metastable states observed in the product reaction channel are assigned by calculating the energy levels and wave functions of the HF...D van der Waals complex. The behavior of resonance energies, widths, and decay branching ratios as functions of total angular momentum is analyzed. The effect of isotopic substitution on resonance energies and lifetimes is elucidated by comparison with previous results for the F+H2 reaction. It is demonstrated that HF(v'=3) products near threshold are formed by decay of the narrow resonances supported by van der Waals wells in the exit channel. State-to-state differential cross sections in the HF(v'=3) channel exhibit characteristic forward-backward peaks due to the formation of a long-lived metastable complex. The role of the exit-channel resonances in the interpretation of molecular beam experiments is discussed.

On the origin of the forward peak and backward oscillations in the F + H(2)(v = 0) --> HF(v' = 2) + H reaction.

We present a semiclassical complex angular momentum (CAM) analysis of the forward scattering peak which occurs at a translational collision energy around 32 meV in the quantum mechanical calculations for the F + H(2)(v = 0, j = 0) --> HF(v' = 2, j' = 0) + H reaction on the Stark-Werner potential energy surface. The semiclassical CAM theory is modified to cover the forward and backward scattering angles. The peak is shown to result from constructive/destructive interference of the two Regge states associated with two resonances, one in the transition state region and the other in the exit channel van der Waals well. In addition, we demonstrate that the oscillations in the energy dependence of the backward differential cross section are caused by the interference between the direct backward scattering and the decay of the two resonance complexes returning to the backward direction after one full rotation.

Normal and hyperspherical mode analysis of NO-doped Kr crystals upon Rydberg excitation of the impurity.

Molecular dynamics simulations and both normal mode and hyperspherical mode analyses of NO-doped Kr solid are carried out in order to get insights into the structural relaxation of the medium upon electronic excitation of the NO molecule. A combined study is reported on the time evolution of the cage radius and on the density of vibrational states, according to the hyperspherical and normal mode analyses. For the hyperspherical modes, hyper-radial and grand angular contributions are considered. For the normal modes, radial and tangential contributions are examined. Results show that the first shell radius dynamics is driven by modes with frequencies at approximately 47 and approximately 15 cm-1. The first one is related to the ultrafast regime where a large part of the energy is transmitted to the lattice and the second one to relaxation and slow redistribution of the energy. The density of vibrational states gamma(omega) is characterized by a broad distribution of bands peaking around the frequencies of approximately 13, approximately 19, approximately 25, approximately 31, approximately 37, approximately 47, and approximately 103 cm-1 (very small band). The dominant modes in the relaxation process were at 14.89, 23.49, and 53.78 cm-1; they present the largest amplitudes and the greatest energy contributions. The mode at 14.89 cm-1 is present in both the fit of the first shell radius and in the hyper-radial kinetic energy spectrum and resulted the one with the largest amplitude, although could not be revealed by the total kinetic energy power spectrum.

Overlapping resonances and Regge oscillations in the state-to-state integral cross sections of the F+H2 reaction.

A Regge pole analysis is employed to explain the oscillatory patterns observed in numerical simulations of integral cross section for the F+H(2)(v=0,j=0)-->HF(v(')=2,j(')=0)+H reaction in the translational collision energy range 25-50 meV. In this range the integral cross section for the transition, affected by two overlapping resonances, shows nearly sinusoidal oscillations below 38 meV and a more structured oscillatory pattern at larger energies. The two types of oscillations are related to the two Regge trajectories which (pseudo) cross near the energy where the resonances are aligned. Simple estimates are given for the periods of the oscillations.

Interacting resonances in the F+H2 reaction revisited: complex terms, Riemann surfaces, and angular distributions.

We study the effect of overlapping resonances on the angular distributions of the reaction F+H2(v=0,j=0)-->HF(v=2,j=0)+H in the collision energy range from 5 to 65 meV, i.e., under the reaction barrier. Reactive scattering calculations were performed using the hyperquantization algorithm on the potential energy surface of Stark and Werner [J. Chem. Phys. 104, 6515 (1996)]. The positions of the Regge and complex energy poles are obtained by Pade reconstruction of the scattering matrix element. The Sturmian theory is invoked to relate the Regge and complex energy terms. For two interacting resonances, a two-sheet Riemann surface is contracted and inverted. The semiclassical complex angular momentum analysis is used to decompose the scattering amplitude into the direct and resonance contributions.

Observed molecular alignment in gaseous streams and possible chiral effects in vortices and in surface scattering.

Extensive work in this laboratory has been devoted to the study of intermolecular interactions from scattering experiments, in order to provide ingredients for modelling forces acting in systems involving hydrocarbons, the components of atmospheres, and water. Our detection of aligned oxygen in gaseous streams and further evidence on simple molecules has been extended to benzene and various hydrocarbons. Chiral effects can be seen in the differential scattering of oriented molecules, in particular from surfaces. It is pointed out that it may be of pre-biotical interest that we focus on possible mechanisms for chiral bio-stereochemistry of oriented reactants, for example when flowing in atmospheres of rotating bodies, specifically the planet earth, as well as in vortex motions of celestial objects. Molecular dynamics simulations and experimental verifications are in progress.

Isomerization dynamics and thermodynamics of ionic argon clusters.

The dynamics and thermodynamics of small Ar(n) (+) clusters, n=3, 6, and 9, are investigated using molecular dynamics (MD) and exchange Monte Carlo (MC) simulations. A diatomic-in-molecule Hamiltonian provides an accurate model for the electronic ground state potential energy surface. The microcanonical caloric curves calculated from MD and MC methods are shown to agree with each other, provided that the rigorous conservation of angular momentum is accounted for in the phase space density of the MC simulations. The previously proposed projective partition of the kinetic energy is used to assist MD simulations in interpreting the cluster dynamics in terms of inertial, internal, and external modes. The thermal behavior is correlated with the nature of the charged core in the cluster by computing a dedicated charge localization order parameter. We also perform systematic quenches to establish a connection with the various isomers. We find that the Ar(3) (+) cluster is very stable in its linear ground state geometry up to about 300 K, and then isomerizes to a T-shaped isomer in which a quasineutral atom lies around a charged dimer. In Ar(6) (+) and Ar(9) (+), the covalent trimer core is solvated by neutral atoms, and the weakly bound solvent shell melts at much lower energies, occasionally leading to a tetramer or pentamer core with weakly charged extremities. At high energies the core itself becomes metastable and the cluster transforms into Ar(2) (+) solvated by a fluid of neutral argon atoms.

Orientation of benzene in supersonic expansions, probed by IR-laser absorption and by molecular beam scattering.

This work represents the first experimental demonstration that planar molecules tend to travel as a "frisbee" when a gaseous mixture with lighter carriers expands into a vacuum, the orientation being due to collisions. The molecule is benzene, the prototype of aromatic chemistry. The demonstration is via two complementary experiments: interrogating benzene by IR-laser light and controlling its orientation by selective scattering on rare gas targets. The results cast new light on the microscopic mechanisms of collisional alignment and suggest a useful way to produce intense beams of aligned molecules, permitting studies of steric effects in gas-phase processes and in surface catalysis.