Mixed quantum/classical theory (MQCT) approach to the dynamics of molecule-molecule collisions in complex systems.

We developed a general theoretical approach and a user-ready computer code that permit study of the dynamics of collisional energy transfer and ro-vibrational energy exchange in complex molecule-molecule collisions. The method is a mixture of classical and quantum mechanics. The internal ro-vibrational motion of collision partners is treated quantum mechanically using a time-dependent Schrödinger equation that captures many quantum phenomena including state quantization and zero-point energy, propensity and selection rules for state-to-state transitions, quantum symmetry and interference phenomena. A significant numerical speed up is obtained by describing the translational motion of collision partners classically, using the Ehrenfest mean-field trajectory approach. Within this framework a family of approximate methods for collision dynamics is developed. Several benchmark studies for diatomic and triatomic molecules, such as H2O and ND3 collided with He, H2 and D2, show that the results of MQCT are in good agreement with full-quantum calculations in a broad range of energies, especially at high collision energies where they become nearly identical to the full quantum results. Numerical efficiency of the method and massive parallelism of the MQCT code permit us to embrace some of the most complicated collisional systems ever studied, such as C6H6 + He, CH3COOH + He and H2O + H2O. Application of MQCT to the collisions of chiral molecules such as CH3CHCH2O + He, and to molecule-surface collisions is also possible and will be pursued in the future.

Mixed quantum/classical calculations of rotationally inelastic scattering in the CO + CO system: a comparison with fully quantum results.

An updated version of the CO + CO potential energy surface from [R. Dawes, X. G. Wang and T. Carrington, J. Phys. Chem. A 2013, 117, 7612] is presented, that incorporates an improved treatment of the asymptotic behavior. It is found that this new surface is only slightly different from the other popular PES available for this system in the literature [G. W. M. Vissers, P. E. S. Wormer and A. Van Der Avoird, Phys. Chem. Chem. Phys. 2003, 5, 4767]. The differences are quantified by expanding both surfaces over a set of analytic functions and comparing the behavior of expansion coefficients along the molecule-molecule distance R. It is shown that all expansion coefficients behave similarly, except in the very high energy range at small R where the PES is repulsive. That difference has no effect on low collision-energy dynamics, which is explored via inelastic scattering calculations carried out using the MQCT program which implements the mixed quantum/classical theory for molecular energy exchange processes. The validity of MQCT predictions of state-to-state transition cross sections for CO + CO is also tested by comparison against full-quantum coupled-states calculations. In all cases MQCT gives reliable results, except at very low collision energy where the full-quantum calculations predict strong oscillations of state-to-state transition cross sections due to resonances. For strong transitions with large cross sections, the results of MQCT are reliable, especially at higher collision energy. For weaker transitions, and lower collision energies, the cross sections predicted by MQCT may be up to a factor of 2-3 different from those obtained by full-quantum calculations.

Symmetry Breaking: A Classic Example of Quantum Interference Captured by Mixed Quantum/Classical Theory.

The phenomena of propensity and inverse propensity are explored using time-dependent mixed quantum classical theory, MQCT, in which the rotational motion of the molecule is treated quantum mechanically, whereas the scattering process is described classically. Good agreement with the results of accurate full-quantum calculations is reported for a closed shell approximation to the NO + Ar system. It is shown that MQCT reproduces both phenomena in a broad range of the final states of the molecule and for various initial rotational states, offering a unique time-dependent insight. It permits seeing that both propensity and inverse propensity occur due to efficient depopulation of some states at the early postcollisional stage of the scattering process, when the molecule exists in a coherent superposition of many excited states that span a very broad range of angular momentum quantum numbers, populated by an efficient stepladder process of many consecutive transitions with small Δj.

Mixed quantum/classical theory for rotational energy exchange in symmetric-top-rotor + linear-rotor collisions and a case study of the ND3 + D2 system.

The extension of mixed quantum/classical theory (MQCT) to describe collisional energy transfer is developed for a symmetric-top-rotor + linear-rotor system and is applied to ND3 + D2. State-to-state transition cross sections are computed in a broad energy range for all possible processes: when both ND3 and D2 molecules are excited or both are quenched, when one is excited while the other is quenched and vice versa, when the ND3 state changes its parity while D2 is excited or quenched, and when ND3 is excited or quenched while D2 remains in the same state, ground or excited. In all these processes the results of MQCT are found to approximately satisfy the principle of microscopic reversibility. For a set of sixteen state-to-state transitions available from the literature for a collision energy of 800 cm-1 the values of cross sections predicted by MQCT are within 8% of accurate full-quantum results. A useful time-dependent insight is obtained by monitoring the evolution of state populations along MQCT trajectories. It is shown that, if before the collision, D2 is in its ground state, the excitation of ND3 rotational states proceeds through a two-step mechanism in which the kinetic energy of molecule-molecule collision is first used to excite D2 and only then is transferred to the excited rotational states of ND3. It is found that both potential coupling and Coriolis coupling play important roles in ND3 + D2 collisions.

Description of quantum interference using mixed quantum/classical theory of inelastic scattering.

Manifestation of the quantum interference effect in the oscillation of scattering cross section is explored using the N2 + O system as a case study. Calculations are carried out for two electronic PESs of the system, for various initial rotational states of N2, in a broad range of N2 + O collision energies and using three theoretical methods: two versions of the approximate mixed quantum/classical theory (MQCT and AT-MQCT) and the accurate full-quantum coupled-channel method (implemented in MOLSCAT). A good agreement between different methods is observed, especially at high energies. Elastic scattering cross-sections oscillate as a function of collision energy, which is the result of quantum interference. The effects of initial rotational excitation and of the PES properties are studied in detail. For the final (thermally averaged) cross sections, both MOLSCAT and MQCT calculations predict a rather regular pattern of quantum oscillations that persist through a broad range of collision energies and expand into the low-energy regime where quantum scattering resonances are common. The difference between cross sections predicted by MQCT and MOLSCAT decreases from ∼8% at low energies to ∼2% at high energies. Experimental data available at high collision energies are well reproduced.

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 H2O + H2O 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).