Scattering of CO with H2O: Statistical and classical alternatives to close-coupling calculations

Adiabatic Trajectory Approximation: A New General Method in the Toolbox of Mixed Quantum/Classical Theory for Collisional Energy Transfer

A new version of the MQCT program is presented, which includes an important addition, adiabatic trajectory approximation (AT-MQCT), in which the equations of motion for the classical and quantum parts of the system are decoupled. This method is much faster, which permits calculations for larger molecular systems and at higher collision energies than was possible before. AT-MQCT is general and can be applied to any molecule + molecule inelastic scattering problem. A benchmark study is presented for H₂O + H₂O rotational energy transfer, an important asymmetric-top rotor + asymmetric-top rotor collision process, a very difficult problem unamenable to the treatment by other codes that exist in the community. Our results indicate that AT-MQCT represents a reliable computational tool for prediction of collisional energy transfer between the individual rotational states of two molecules, and this is valid for all combinations of state symmetries (such as para and ortho states of each collision partner).

Collisonal energy transfer in the CO-CO system

An accurate determination of the physical conditions in astrophysical environments relies on the modeling of molecular spectra. In such environments, densities can be so low ($n << 10^{10}$ cm$^{-3}$) that...

HD-H+ collisions: statistical and quantum state-to-state studies

In the early Universe, the cooling mechanisms of the gas significantly rely on the HD abundance and excitation conditions. A proper modeling of its formation and destruction paths as well as its excitation by both radiative and collisional processes is then required to accurately describe the cooling mechanisms of the pristine gas. In such media, ion-molecule reactions are dominant. Their theoretical study is challenging and state-of-the-art quantum time-independent methods are computationally limited to collisions involving light molecules. Here, we report a state-to-state scattering study of the HD-H⁺ collisional system using two different methods: an exact quantum time-independent approach and a recently developed fast and efficient statistical method. Reactive and inelastic rate coefficients were obtained for temperatures up to 300 K. The statistical method is able to reproduce exact calculations with an accuracy reaching the astrophysical needs while drastically reducing the computational resources requirements. Such results suggest that this new statistical method should be considered to provide the astrophysical community collisional data for which quantum calculations are impossible.

Benchmarking an improved statistical adiabatic channel model for competing inelastic and reactive processes

Inelastic collisions and elementary chemical reactions proceeding through the formation and subsequent decay of an intermediate collision complex, with an associated deep well on the potential energy surface, pose a challenge for accurate fully quantum mechanical approaches, such as the close-coupling method. In this study, we report on the theoretical prediction of temperature-dependent state-to-state rate coefficients for these complex-mode processes, using a statistical quantum method. This statistical adiabatic channel model is benchmarked by a direct comparison using accurate rate coefficients from the literature for a number of systems (H₂ + H⁺, HD + H⁺, SH⁺ + H, and CH⁺ + H) of interest in astrochemistry and astrophysics. For all of the systems considered, an error of less than factor 2 was found, at least for the dominant transitions and at low temperatures, which is sufficiently accurate for applications in the above mentioned disciplines.

Potential energy surface and bound states of the H2O-HF complex

We present the first global five-dimensional potential energy surface for the H₂O-HF dimer, a prototypical hydrogen bonded complex. Large scale ab initio calculations were carried out using the explicitly correlated coupled cluster approach with single- and double-excitations together with non-iterative perturbative treatment of triple excitations with the augmented correlation-consistent triple zeta basis sets, in which the water and hydrogen fluoride monomers were frozen at their vibrationally averaged geometries. The ab initio data points were fitted to obtain a global potential energy surface for the complex. The equilibrium geometry of the complex corresponds to the formation of a hydrogen bond with water acting as a proton acceptor and a binding energy of D_e = 3059 cm^-1 (8.75 kcal/mol). The energies and wavefunctions of the lowest bound states of the complex were computed using a variational approach, and the dissociation energies of both ortho-H₂O-HF (D₀ = 2089.4 cm^-1 or 5.97 kcal/mol) and para-H₂O-HF (D₀ = 2079.6 cm^-1 or 5.95 kcal/mol) were obtained. The rotational constant of the complex was found to be in good agreement with the available experimental data.

When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111)

The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.

State-to-state inelastic rotational cross sections in five-atom systems with the multiconfiguration time dependent Hartree method

We present a MultiConfiguration Time Dependent Hartree (MCTDH) method as an attractive alternative approach to the usual quantum close-coupling method that approaches some computational limits in the calculation of rotational excitation (and de-excitation) between polyatomic molecules (here collisions between triatomic and diatomic rigid molecules). We have performed a computational investigation of the rotational (de-)excitation of the benchmark rigid rotor H₂O-H₂ system on a recently developed Potential Energy Surface of the complex using the MCTDH method. We focus here on excitations and de-excitations from the 0₀₀, 1₁₁, and 1₁₀ states of H₂O with H₂ in its ground rotational state, looking at all the potential transitions in the energy range 1-200 cm^-1. This work follows a recently completed study on the H₂O-H₂ cluster where we characterized its spectroscopy and more generally serves a broader goal to describe inelastic collision processes of high dimensional systems using the MCTDH method. We find that the cross sections obtained from the MCTDH calculations are in excellent agreement with time independent calculations from previous studies but does become challenging for the lower kinetic energy range of the de-excitation process: that is, below approximately 20 cm^-1 of collision energy, calculations with a relative modest basis become unreliable. The MCTDH method therefore appears to be a useful complement to standard approaches to study inelastic collision for various collision partners, even at low energy, though performing better for rotational excitation than for de-excitation.

Theoretical Study of Barrierless Chemical Reactions Involving Nearly Elastic Rebound: The Case of S(1D) + X2, X = H, D

For some values of the total angular momentum consistent with reaction, the title processes involve nonreactive trajectories proceeding through a single rebound mechanism during which the internal motion of the reagent diatom is nearly unperturbed. When such paths are in a significant amount, the classical reaction probability is found to be markedly lower than the quantum mechanical one. This finding was recently attributed to an unusual quantum effect called diffraction-mediated trapping, and a semiclassical correction was proposed in order to take into account this effect in the classical trajectory method. In the present work, we apply the resulting approach to the calculation of opacity functions as well as total and state-resolved integral cross sections (ICSs) and compare the values obtained with exact quantum ones, most of which are new. As the title reactions proceed through a deep insertion well, mean potential statistical calculations are also presented. Seven values of the collision energy, ranging from 30 to 1127 K, are considered. Two remarkable facts stand out: (i) The corrected classical treatment strongly improves the accuracy of the opacity function as compared to the usual classical treatment. When the entrance transition state is tight, however, those trajectories crossing it with a bending vibrational energy below the zero point energy must be discarded. (ii) The quantum opacity function, particularly its cutoff, is finely reproduced by the statistical approach. Consequently, the total ICS is also very well described by the two previous approximate methods. These, however, do not predict state-resolved ICSs with the same accuracy, proving thereby that (i) one or several genuine quantum effects involved in the dynamics are missed by the corrected classical treatment and (ii) the dynamics are not fully statistical.

An accurate full-dimensional permutationally invariant potential energy surface for the interaction between H2O and CO

The first full-dimensional accurate potential energy surface was developed for the CO + H₂O system based onca.102 000 points calculated at the CCSD(T)-F12a/AVTZ level using a permutation invariant polynomial-neural network (PIP-NN) method.

The water–carbon monoxide dimer: new infrared spectra, ab initio rovibrational energy level calculations, and an interesting in-termolecular mode

Bound state rovibrational energy level calculations using a high-level intermolecular potential surface are reported for H₂O–CO and D₂O–CO.