Cálculos mecano-cuánticos y cuasiclásicos de la dinámica y etereodinámica de reacciones bimoleculares y procesos inelásticos.

 

Vicente Sáez Rábanos

 

Javier Aoiz (U.C.M.)

Enrique Verdasco (U.C.M.)

Victor Herrero (C.S.I.C.)

Marcelo P. de Miranda (U.Leeds, UK)

 
                                                                                                     

 

 

En los últimos años se han aplicado métodos mecano-cuánticos (QM) y de trayectorias cuasi-clásicas (QCT) para el estudio de la dinámica de una serie de reacciones de abstracción, de inserción, así como otras cuyos mecanismos son, en general, más complejos, tales como Li+HF y O(1D)+HCl . El comportamiento del método QCT ha sido contrastado en numerosas ocasiones con resultados QM exactos, parte de los cuales han sido realizados en nuestro grupo. En cualquier caso, el método QCT permite la identificación y caracterización de efectos cuánticos relevantes en la dinámica de reacciones.

 

 Se han desarrollado métodos para el cálculo de probabilidades de reacción para cualquier estado rotacional inicial y probabilidades cumulativas de reacción (CRP), comparando resultados QM y QCT. El acuerdo entre los resultados QM exactos y QCT es notablemente bueno para las reacciones examinadas. Como una extensión de este trabajo, se ha desarrollado un nuevo modelo estadístico QCT (SQCT) en el último año . El modelo garantiza que se cumpla el principio de balance detallado y conservación de la paridad, lo que no ocurre en los cálculos QCT convencionales. El acuerdo con los resultados obtenidos por el modelo SQM es extraordinariamente bueno, hasta el punto de que es capaz de mostrar las inexactitudes relativamente pequeñas resultantes del uso de la aproximación CS (centrifugal sudden) en los cálculos QM

 

La estereodinámica de reacciones elementales es de especial interés dados el progreso experimental y los avances metodológicos en la interpretación y análisis de los datos. Este particular subcampo trata del papel de los factores direccionales en la reactividad. Específicamente de cómo las distribuciones de ejes moleculares o de los momentos angulares de rotación (correlaciones vectoriales k-r-, k-j-, donde k, y j son las velocidades y momentos angulares iniciales y finales, respectivamente, y r es la dirección del eje molecular) influencian la reactividad y cuáles son las distribuciones de los momentos angulares de los productos de la reacción y su dependencia con el ángulo de dispersión cuando la reacción tiene lugar (correlación vectorial kk´j´). En colaboración con el Dr. Miranda (Universidad de Leeds,UK) se  ha conseguido un gran progreso en el desarrollo teórico y aplicación de métodos QM y en la unificación de los enfoques QM y QCT para la descripción exacta de la estereodinámica de reacciones átomo-diátomo, así como en el desarrollo de una teoría rigurosa en términos de distribuciones de poblaciones de estados de mínima incertidumbre para visualizar las distribuciones de momentos angulares (‘retratos’) que dan lugar a reacción . En particular, por primera vez, se ha estudiado la estereodinámica de reacciones a baja y ultra-baja energía de colisión e investigado la capacidad potencial de controlar esas reacciones.

 

Quantum mechanical and Quasi-classical Trajectory dynamical and stereodynamical calculations of bimolecular chemical reactions and inelastic processes

 

Over the last years we have applied quantum mechanical (QM) and Quasi-Classical Trajectory (QCT) methods to study the dynamics of a series of abstraction reactions, insertion reactions and other reactions whose mechanisms are, in general, more complex, such Li+HF and O(1D)+HCl. The performance of the QCT method has been contrasted in many occasions with exact QM calculations, part of them carried out also in our group. In any case, the QCT method permits the identification and characterization of relevant quantum effects in reaction dynamics.

 

We have developed methods to calculate and analyze the reaction probabilities (RP) for any initial rotational angular momentum and cumulative reactive probabilities (CRP) comparing the QM and QCT results. The agreement between QM exact results and QCT data is remarkably good for the reactions examined. As an extension of this work, a new statistical QCT (SQCT) model has been developed in the last year. The model guarantees the detailed balance and the conservation of parity. The agreement with the SQM is extraordinarily good, to such extent that the SQCT model is sensitive enough to show the minor inaccuracies when the the centrifugal sudden approximation is used in the QM calculations.

 

The stereodynamics of elementary chemical reactions is of special interest given the experimental progress and methodological advances in the interpretation and analysis of the data. This particular subfield deals with the role of directional factors on the reactivity. Specifically, how the spatial distributions of molecular axes or rotational angular momenta (k-r-k´, k-j-k´ vector correlations, where k, k´ and j j´ are initial and final relative velocities rotational angular momenta, respectively, and r is the direction of the molecular axis) influence the reactivity , and which are the distributions of product angular momentum and their dependence with the scattering angle when the reaction has taken place (kk'j' vector correlation). In collaboration with Dr. Marcelo P. de Miranda (University of Leeds, UK)  we have achieved a great progress in the theoretical development and application of QM methods and the unification of the QM and QCT approaches for the exact description of the stereodynamics of atom-diatom reactions, as well as in the development of a rigorous theory in terms of distributions of minimum-uncertainty angular momentum populations to visualize the spatial angular momentum distributions (‘portraits’) leading to reaction. In particular, for the first time, we have study the stereodynamics of reactions at low and ultra-low collision energies and investigate the potential capabilities of controlling these reactions the F(2P3/2 , 2P1/2)+H2 and Cl(2P3/2 , 2P1/2)+H2 for which calculations on new potential PESs have become recently available [39].

 

 

Publicaciones recientes /Recent publications

 

-F. J. Aoiz, J.E. Verdasco, V.J. Herrero,  V. Sáez Rábanos, M.A. Alexander. Attractive and repulsive interactions in the inelastic scattering of NO by Ar : A comparison between classical trajectory and close-coupling quantum mechanical results.  J. Chem. Phys 119,  5860- 58666, 2003

-Teresa Martínez, María Luz Hernandéz, José María Alvariño, F.J. Aoiz, V. Sáez Rábanos. A detailed study of the dynamics of the  O( 1D )+HCl  + Cl , ClO + H reactions. J. Chem. Phys. 119, 7871- 7885, 2003

-F. J. Aoiz, V.J. Herrero,  V. Sáez Rábanos, J.E. Verdasco. Classical stereodynamics in Ar + NO inelastic collisions.

PCCP.  6, 4407-4415, 2004

- Andrew E. Pomerantz, Florian Ausfelder, Richard N. Zare, Juan Carlos Juanes-Marcos, Stuart C. Althorpe,  V. Sáez Rábanos, F. J. Aoiz, L. Bañares, J. F. Castillo. Rovibrational product state distribution for inelastic H + D2 collisions.

J. Chem. Phys. 121,6587- 6590, 2004.

-Marcelo P. De Miranda, F. J. Aoiz, V. Sáez Rábanos, Mark Brouard. Spatal distributions of angular momenta in quantum and quasiclassical stereodynamics..J. Chem. Phys. 121, 9830- 9843, 2004.

-Jesús Aldegunde, Marcelo P. de Miranda, James M. Haigh, Brian K. Kendrick,  V. Sáez Rábanos,  F. J. Aoiz. How reactants polarization can be used to change and unravel chemical reactivity. J. Phys. Chem. A. 109,  6200-6217, 2005

- F. J. Aoiz,  V. Sáez Rábanos, Bruno Martinez-Haya,  Tomas Gonzalez-Lezana.Quasiclassical determination of reaction probabilities as a function of the total angular momentum. J. Chem. Phys. 123, 094101-1- 094101-14, 2005

-Jesús Aldegunde, J.M. Alvariño, Marcelo P. de Miranda, V. Sáez Rábanos,  F. J. Aoiz.

 Mechamisn and Control of the F + H2 reaction at low and ultralow collision energies.

 J. Chem. Phys. 125, 133104-1-133104-12, 2006

- F. J. Aoiz , M. Brouard, C.J. Eyles, J.F. Castillo, V. Sáez Rábanos.  Cumulative reaction probabilities. A comparison between quasiclassical and quantum mechanical results. J. Chem. Phys.125, 144105-1- 44105-14, 2006

-J. Aldegunde, J.M. Alvariño, B.K. Kendrick,  V. Sáez Rábanos, M. P. de Miranda, F. J. Aoiz. Analysis of the H+D2 reaction mechanism through consideration of the intrinsic reactant polarisation.  PCCP. 8,  4881-4896,  2006

-F. J. Aoiz,  V. Sáez Rábanos,  Tomas Gonzalez-Lezana, D. E. Manolopoulos. A statistical quasiclassical trajectory model for atom-diatom isertion reactions. J. Chem. Phys. 126, 161101-1- 161101-5,  2007

- F. J. Aoiz, V. J. Herrero, M. P. de Miranda, V. Sáez Rábanos  Constraints at the transition state of the D + H2 reaction: quantum bottlenecks vs. stereodynamics. PCCP.9, 5367- 5373, 2007

-F. J.  Aoiz,  T. Gonzalez-Lezama,  V. Sáez Rábanos . A stringent test of the statistical quasiclassi cal trajectory model for the H+3 exchange reaction:  A comparison with rigorous statistical quantum mechanical results. J. Chem. Phys. 127, 174109-1- 174109-15, 2007

 -Jesús Aldegunde,  F. Javier Aoiz,  Vicente Sáez-Rábanos,   Brian  K. Kendrick, Marcelo P. de Miranda. The canonical and  other  mechanisms of elementary chemical reactions. PCCP.9,  5794- 5808, 2007