Search

In this work, we present a global potential energy surface for the ground electronic state of the H4+ based on ab initio calculations. The final fit is based on triatomics-in-molecules (TRIM) approximation and it includes extra four-body terms for the better description of some discrepancies found on the TRIM model. The TRIM method itself allows a very accurate description of the asymptotic regions. The global fit uses more than 19 000 multireference configuration interaction ab initio points. The global potential energy surface has an overall root mean square error of 0.013 eV for energies up to 2 eV above the global minimum. This work presents an analysis of the stationary points, reactant and product channels, and crossing between the two lowest TRIM adiabatic states. It is as well included a brief description of the two first excited states of the TRIM matrix, concluding that TRIM method is a very good approximation not only for the ground state but also for at least two of the excited states of H4+ system. [ABSTRACT FROM AUTHOR]

Abstract: We present the first comprehensive ab initio study of the Ne–Br2 potential energy surfaces and the non-adiabatic couplings between the valence excited electronic states. These ab initio results are used to obtain 3-D approximate potentials for each electronic state, and these potentials are used in a wave packet calculation of the competing electronic predissociation and vibrational predissociation dynamics. The results of this calculation are in excellent agreement with both experimental results and a previous empirical fit to the experiments. The calculations allow us to observe not only the competition between vibrational and electronic dynamics for the dimer, but also the competition between two different electronic channels. Coupling to the 2 g state dominates for the levels studied here, but coupling to the C state is progressively more important for low vibrational levels, and may dominate at levels below which the current results pertain. The ability of ab initio surfaces and couplings to so accurately reproduce experimental data raises the hope of a complete understanding of the VP and EP dynamics for other Rg-halogen dimers. Success in the case presented here is largely due to the fact that the VP dynamics for the vibrational levels in this study are in the simple, direct regime. Understanding the simple case so thoroughly provides new hope that the more complicated examples, such as ArI2 and NeCl2, for which experiment and theory are not currently in accord, may yet yield to analysis. [Copyright &y& Elsevier]

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 H5+ 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 H5+ 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. [ABSTRACT FROM AUTHOR]