Lola González-Sánchez

Dynamically biased statistical model for the ortho/para conversion in the H2+H3+ → H3++ H2 reaction

In this work we present a dynamically biased statistical model to describe the evolution of the title reaction from statistical to a more direct mechanism, using quasi-classical trajectories (QCT). The method is based on the one previously proposed by Park and Light [J. Chem. Phys. 126, 044305 (2007)]. A recent global potential energy surface is used here to calculate the capture probabilities, instead of the long-range ion-induced dipole interactions. The dynamical constraints are introduced by considering a scrambling matrix which depends on energy and determine the probability of the identity/hop/exchange mechanisms. These probabilities are calculated using QCT. It is found that the high zero-point energy of the fragments is transferred to the rest of the degrees of freedom, what shortens the lifetime of H(5)(+) complexes and, as a consequence, the exchange mechanism is produced with lower proportion. The zero-point energy (ZPE) is not properly described in quasi-classical trajectory calculations and an approximation is done in which the initial ZPE of the reactants is reduced in QCT calculations to obtain a new ZPE-biased scrambling matrix. This reduction of the ZPE is explained by the need of correcting the pure classical level number of the H(5)(+) complex, as done in classical simulations of unimolecular processes and to get equivalent quantum and classical rate constants using Rice-Ramsperger-Kassel-Marcus theory. This matrix allows to obtain a ratio of hop/exchange mechanisms, α(T), in rather good agreement with recent experimental results by Crabtree et al. [J. Chem. Phys. 134, 194311 (2011)] at room temperature. At lower temperatures, however, the present simulations predict too high ratios because the biased scrambling matrix is not statistical enough. This demonstrates the importance of applying quantum methods to simulate this reaction at the low temperatures of astrophysical interest.

Quantum stereodynamics of the F+OH(v, j) reactive collision on the 13A″ state

The reactivity dependence on the k, j vector correlation for the title reaction is studied for j=1, 2, 3 by a quantum wavepacket method in a large range of translational energy on a recently obtained 3A″ adiabatic potential energy surface. It is found that the reaction is extraordinarily enhanced when initially k and j are perpendicular, or equivalently, when the OH internuclear vector is parallel to the initial relative velocity in the rotating body-fixed plane. Such a result is interpreted in terms of the kinematic constraints imposed by the small skew angle in heavy–light + heavy reactions, being a general situation in this kind of HL + H → H + LH reactions. The total reaction cross-sections obtained are in rather good agreement with quasiclassical results, especially as j increases. Finally, the rate constants obtained are compared with available experimental data at nearly room temperature.

Dynamical regimes on the Cl + H2 collisions: inelastic rainbow scattering.

While Cl + H(2) reactive collisions have been a subject of numerous experimental and theoretical studies, inelastic collisions leading to rotational energy transfer and/or vibrational excitation have been largely ignored. In this work, extensive quantum mechanical calculations covering the 0.5-1.5 eV total energy range and various initial rovibrational states have been carried out and used to perform a joint study of inelastic and reactive Cl + H(2) collisions. Quasiclassical trajectories calculations complement the quantum mechanical results. The analysis of the inelastic transition probabilities has revealed the existence of two distinct dynamical regimes that correlate with low and high impact parameters, b, and are neatly separated by glory scattering. It has been found that while high-b collisions are mainly responsible for |Δj| = 2 transitions which dominate the inelastic scattering, they are very inefficient in promoting higher |Δj| transitions. The effectiveness of this type of collision also drops with rotational excitation of H(2). In contrast, reactive scattering, that competes with |Δj| > 2 inelastic transitions, is exclusively caused by low-b collisions, and it is greatly favored when the reactants get rotationally excited. Previous studies focusing on the reactivity of the Cl + H(2) system established that the van der Waals well located in the entrance channel play a key role in determining the mechanism of the collisions. Our results prove this to be also a case for inelastic processes, where the origin of the double dynamical regime can be traced back to the influence exerted by this well that shapes the topology of the entrance channel of the Cl-H(2) system.

Rotational ‘cooling’ and ‘heating’ of OH+(3Σ−) by collisions with He: quantum dynamics revealing propensity rules under ion trap conditions

ABSTRACT Multichannel scattering calculations are presented for the low-energy collisions of the OH+ cation and He atoms, using an ab initio evaluation of the interaction potential, which had been obtained in earlier work, and a time-independent, multichannel treatment of the quantum dynamics carried out in this study using our in-house scattering code ASPIN. Given the presence of spin-rotation coupling effects, within an essentially electrostatic formulation of the interaction forces with He atoms in the trap, the ensuing propensity rules which control the relative size of the state-changing cross sections and of the corresponding inelastic rates, also computed at the most likely temperatures in an ion trap, are presented and analysed in detail. GRAPHICAL ABSTRACT

Modeling Quantum Kinetics in Ion Traps: State-changing Collisions for OH+(3Σ-) Ions with He as a Buffer Gas

Abstract We present quantum scattering calculations for rotational state‐changing cross sections and rates, up to about 50 K of ion translational temperatures, for the OH+ molecular ion in collision with He atoms as the buffer gas in the trap. The results are obtained both by using the correct spin‐rotation coupling of angular momenta and also within a recoupling scheme that treats the molecular target as a pseudo‐(1Σ+ ) state, and then compares our findings with similar data for the OH−(1Σ+ ) molecular partner under the same conditions. This comparison intends to link the cation behaviour to the one found earlier for the molecular anion. The full calculations including the spin‐rotation angular momenta coupling effects have been recently reported (L. González‐Sánchez and R. Wester and F.A. Gianturco, Mol.Phys.2018, DOI 10.1080/00268976.2018.14425971) with the aim of extracting specific propensity rules controlling the size of the cross sections. The present study is instead directed to modelling trap cooling dynamics by further obtaining the solutions of the corresponding kinetics equations under different trap schemes so that, using the presently computed rates can allow us to indicate specific optimal conditions for the experimental setup of the collisional rotational cooling in an ion trap for the system under study.

Influence of the Reactants Rotational Excitation on the H + D2(v = 0, j) Reactivity

We have analyzed the influence of the rotational excitation on the H + D2(v = 0, j) reaction through quantum mechanical (QM) and quasiclassical trajectories (QCT) calculations at a wide range of total energies. The agreement between both types of calculations is excellent. We have found that the rotational excitation largely increases the reactivity at large values of the total energy. Such an increase cannot be attributed to a stereodynamical effect but to the existence of recrossing trajectories that become reactive as the target molecule gets rotationally excited. At low total energies, however, recrossing is not significant and the reactivity evolution is dominated by changes in the collision energy; the reactivity decreases with the collision energy as it shrinks the acceptance cone. When state-to-state results are considered, rotational excitation leads to cold products rovibrational distributions, so that most of the energy is released as recoil energy.

Reaction dynamics and mechanism of the Cl + HD(v = 1) reaction: a quantum mechanical study.

Time-independent quantum mechanical calculations have been performed in order to characterize the dynamics and stereodynamics of Cl + HD reactive collisions. Calculations have been carried out at two different total energy values and for various initial states using the adiabatic potential energy surface by Bian and Werner [J. Chem. Phys. 2000, 112, 220]. Special attention has been paid to the reaction with HD(v = 1) for which integral and differential cross-sections have been calculated and the effect of vibrational vs translational energy on the reactivity has been examined. In addition, the reactant polarization parameters and polarization-dependent differential cross-sections have been determined. From these results, the spatial preferences of the reaction and the extent of the control of the cross sections achievable through a suitable preparation of the reactants have been also studied. The directional requirements are tighter for the HCl channel than for the DCl one. Formation of the products takes place preferentially when the rotational angular momentum of the HD molecule is perpendicular to the reactants approach direction. Cross-sections and polarization moments computed from the scattering calculations have been compared with experimental results by Kandel et al. [J. Chem. Phys. 2000, 112, 670] for the reaction with HD(v = 1) produced by stimulated Raman pumping. The agreement so obtained is good, and it improves the accordance found in previous calculations with other methodologies and potential energy surfaces.