Rotational excitation cross sections of para-H2+para-H2 collisions. A full-dimensional wave-packet propagation study using an exact form of the kinetic energy

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.

New Full-Dimensional Reactive Potential Energy Surface for the H4 System

As the most abundant molecule in the universe, collisions involving H2 have important implications in astrochemistry. Collisions between hydrogen molecules also represent a prototype for assessing various dynamic methods for understanding fundamental few-body processes. In this work, we develop a new and highly accurate full-dimensional potential energy surface (PES) covering all reactive channels of the H2 + H2 system, which extends our previously reported H2 + H2 nonreactive PES [J. Chem. Theory Comput., 2021, 17, 6747] by adding 39,538 additional ab initio points calculated at the MRCI/AV5Z level in the reactive channels. The global PES is represented with high fidelity (RMSE = 0.6 meV for a total of 79,000 points) by a permutation invariant polynomial neural network (PIP-NN) and is suitable for studying collision-induced dissociation, single-exchange, as well as four-center exchange reactions. Preliminary quasi-classical trajectory studies on the new PIP-NN PES reveal strong vibrational enhancement of all reaction channels.

Construction of a Mode-Combination Hamiltonian under the Grid-Based Representation for the Quantum Dynamics of OH + HO2 → O2 + H2O

In this work, a systematic construction framework on a mode-combination Hamiltonian operator of a typical polyatomic reaction, OH + HO2 → O2 + H2O, is developed. First, a set of Jacobi coordinates are employed to construct the kinetic energy operator (KEO) through the polyspherical approach ( Phys. Rep. 2009, 484, 169). Second, due to the multiconfigurational electronic structure of this system, a non-adiabatic potential energy surface (PES) is constructed where the first singlet and triplet states are involved with spin-orbital coupling. To improve the training database, the training set of random energy data was optimized through a popular iterative optimization approach with extensive trajectories. Here, we propose an automatic trajectory method, instead of the classical trajectory on a crude PES, where the gradients are directly computed by the present ab initio calculations. Third, on the basis of the training set, the potential function is directly constructed in the canonical polyadic decomposition (CPD) form ( J. Chem. Theory Comput. 2021, 17, 2702-2713) which is helpful in propagating the nuclear wave function under the grid-based representation. To do this, the Gaussian process regression (GPR) approach for building the CPD form, called the CPD-GPR method ( J. Phys. Chem. Lett. 2022, 13, 11128-11135) is adopted where we further revise CPD-GPR by introducing the mode-combination (mc) scheme leading to the present CPD-mc-GPR approach. Constructing the full-dimension non-adiabatic Hamiltonian operator with mode combination, as test calculations, the nuclear wave function is propagated to preliminarily compute the reactive probability of OH + HO2 → O2 + H2O where the reactants are prepared in vibrational ground states and in the first triplet electronic state.

Cold Collisions of Ro-Vibrationally Excited D2 Molecules

The H₂ + H₂ system has long been considered a benchmark system for ro-vibrational energy transfer in bimolecular collisions. However, most studies thus far have focused on collisions involving H₂ molecules in the ground vibrational level or in the first excited vibrational state. While H₂ + H₂/HD collisions have received wide attention due to the important role they play in astrophysics, D₂ + D₂ collisions have received much less attention. Recently, Zhou et al. [ Nat. Chem. 2022, 14, 658-663, DOI: 10.1038/s41557-022-00926-z] examined stereodynamic aspects of rotational energy transfer in collisions of two aligned D₂ molecules prepared in the v = 2 vibrational level and j = 2 rotational level. Here, we report quantum calculations of rotational and vibrational energy transfer in collisions of two D₂ molecules prepared in vibrational levels up to v = 2 and identify key resonance features that contribute to the angular distribution in the experimental results of Zhou et al. The quantum scattering calculations were performed in full dimensionality and using the rigid-rotor approximation using a recently developed highly accurate six-dimensional potential energy surface for the H₄ system that allows descriptions of collisions involving highly vibrationally excited H₂ and its isotopologues.

Recent advances in quantum theory on ro-vibrationally inelastic scattering

Molecular collisions are of fundamental importance in understanding intermolecular interaction and dynamics. Its importance is accentuated in cold and ultra-cold collisions because of the dominant quantum mechanical nature of the scattering. We review recent advances in the time-independent approach to quantum mechanical characterization of non-reactive scattering in tetratomic systems, which is ideally suited for large collisional de Broglie wavelengths characteristic in cold and ultracold conditions. We discuss quantum scattering algorithms between two diatoms and between a triatom and an atom and their implementation, as well as various approximate schemes. They not only enable the characterization of collision dynamics in realistic systems but also serve as benchmarks for developing more approximate methods.

Wave Packet Approach to Adiabatic and Nonadiabatic Dynamics of Cold Inelastic Scatterings

Due to the extremely large de Broglie wavelength of cold molecules, cold inelastic scattering is always characterized by the time-independent close-coupling (TICC) method. However, the TICC method is difficult to apply to collisions of large molecular systems. Here, we present a new strategy for characterizing cold inelastic scattering using wave packet (WP) method. In order to deal with the long de Broglie wavelength of cold molecules, the total wave function is divided into interaction, asymptotic and long-range regions (IALR). The three regions use different numbers of ro-vibrational basis functions, especially the long-range region, which uses only one function corresponding to the initial ro-vibrational state. Thus, a very large grid range can be used to characterize long de Broglie wavelengths in scattering coordinates. Due to its better numerical scaling law, the IALR-WP method has great potential in studying the inelastic scatterings of larger collision systems at cold and ultracold regimes.

Direct product-type grid representations for angular coordinates in extended space and their application in the MCTDH approach

Multi-configurational time-dependent Hartree (MCTDH) calculations using time-dependent grid representations can be used to accurately simulate high-dimensional quantum dynamics on general ab initio potential energy surfaces. Employing the correlation discrete variable representation, sets of direct product type grids are employed in the calculation of the required potential energy matrix elements. This direct product structure can be a problem if the coordinate system includes polar and azimuthal angles that result in singularities in the kinetic energy operator. In the present work, a new direct product-type discrete variable representation (DVR) for arbitrary sets of polar and azimuthal angles is introduced. It employs an extended coordinate space where the range of the polar angles is taken to be [-π, π]. The resulting extended space DVR resolves problems caused by the singularities in the kinetic energy operator without generating a very large spectral width. MCTDH calculations studying the F·CH₄ complex are used to investigate important properties of the new scheme. The scheme is found to allow for more efficient integration of the equations of motion compared to the previously employed cot-DVR approach [G. Schiffel and U. Manthe, Chem. Phys. 374, 118 (2010)] and decreases the required central processing unit times by about an order of magnitude.

Rotationally inelastic scattering of O3–Ar: state-to-state rates with the multiconfigurational time dependent Hartree method

The rates of state-changing collisions are compared for different isotopologues of ozone from quantum scattering calculations with the MCTDH method.

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.

Controlling rotational quenching rates in cold molecular collisions

The relative orientation and alignment of colliding molecules plays a key role in determining the rates of chemical processes. Here, we examine in detail a prototypical example: rotational quenching of HD in cold collisions with H₂. We show that the rotational quenching rate from j = 2 → 0, in the v = 1 vibrational level, can be maximized by aligning the HD along the collision axis and can be minimized by aligning the HD at the so called magic angle. This follows from quite general helicity considerations and suggests that quenching rates for other similar systems can also be controlled in this manner.

Unraveling the Stereodynamics of Cold Controlled HD-H_{2} Collisions

Measuring inelastic rates with partial-wave resolution requires temperatures close to a Kelvin or below, even for the lightest molecule. In a recent experiment, Perreault, Mukherjee, and Zare [Nat. Chem. 10, 561 (2018).NCAHBB1755-433010.1038/s41557-018-0028-5] studied collisional relaxation of excited HD molecules in the v=1, j=2 state by para- and ortho-H_{2} at a temperature of about 1 K, extracting the angular distribution of scattered HD in the v=1, j=0 state. By state preparation of the HD molecules, control of the angular distribution of scattered HD was demonstrated. Here, we report a first-principles simulation of that experiment which enables us to attribute the main features of the observed angular distribution to a single L=2 partial-wave shape resonance. Our results demonstrate important stereodynamical insights that can be gained when numerically exact quantum scattering calculations are combined with experimental results in the few-partial-wave regime.

Inelastic Scattering of Identical Molecules within Framework of the Mixed Quantum/Classical Theory: Application to Rotational Excitations in H2 + H2

Theoretical foundation is laid out for description of permutation symmetry in the inelastic scattering processes that involve collisions of two identical molecules, within the framework of the mixed quantum/classical theory (MQCT). In this approach, the rotational (and vibrational) states of two molecules are treated quantum-mechanically, whereas their translational motion (responsible for scattering) is treated classically. This theory is applied to H2 + H2 system, and the state-to-state transition cross sections are compared versus those obtained from the full-quantum calculations and experimental results from the literature. Good agreement is found in all cases. It is also found that results of MQCT, where the Coriolis coupling is included classically, are somewhat closer to exact full-quantum results than results of the other approximate quantum methods, where those coupling terms are neglected. These new developments allow applications of MQCT to a broad variety of molecular systems and processes.

Rovibrational energy transfer in ortho-H2+para-H2 collisions

We present the results of a full-dimensional quantum mechanical study of the rovibrational energy transfer in the collision between ortho-H2 and para-H2 in the energy range of 0.1-1.0 eV. The multiconfiguration time-dependent Hartree algorithm has been used to propagate the wave packets on the global potential energy surface by Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] and on a modification of this surface where the short range anisotropy is reduced. State-to-state attributes such as probabilities or integral cross sections are obtained using the formalism of Tannor and Weeks [J. Chem. Phys. 98, 3884 (1993)] by Fourier transforming the correlation functions. The effect of initial rotation of the diatoms on the inelastic and de-excitation processes is investigated.

Using neural networks to represent potential surfaces as sums of products

By using exponential activation functions with a neural network (NN) method we show that it is possible to fit potentials to a sum-of-products form. The sum-of-products form is desirable because it reduces the cost of doing the quadratures required for quantum dynamics calculations. It also greatly facilitates the use of the multiconfiguration time dependent Hartree method. Unlike potfit product representation algorithm, the new NN approach does not require using a grid of points. It also produces sum-of-products potentials with fewer terms. As the number of dimensions is increased, we expect the advantages of the exponential NN idea to become more significant.

Theoretical investigation of highly excited vibrational states in DFCO: Calculation of the out-of-plane bending states and simulation of the intramolecular vibrational energy redistribution

A previously developed modified Davidson scheme [C. Iung and F. Ribeiro, J. Chem. Phys. 121, 174105 (2005)] is applied to compute and analyze highly excited (nu2,nu6) eigenstates in DFCO. The present paper is also devoted to the simulations of the intramolecular vibrational energy redistribution (IVR) initiated by an excitation of the out-of-plane bending vibration (nnu6, n=2,4,6, . . . ,18, and 20). The multiconfiguration time-dependent Hartree method is exploited to propagate the corresponding six-dimensional wave packets. A comprehensive comparison with experimental data as well as with previous simulations of IVR in HFCO [G. Pasin et al. J. Chem. Phys. 124, 194304 (2006)] is presented.

Direct quantum dynamics using variational multi-configuration Gaussian wavepackets. Implementation details and test case

We present here a direct quantum dynamics method using variational multi-configuration Gaussian wavepackets. Based on the efficient multi-configuration time-dependent Hartree wavepacket propagation algorithm, it uses on-the-fly quantum chemical calculation of the potential energy and its derivatives rather than fitted surfaces. Intermediate results are stored in a database so that expensive quantum chemical computations can be recycled. This method is intended to treat quantum effects in the photochemistry of large molecules and the use of Cartesian coordinates to perform direct dynamics is discussed with a comparison between Cartesian coordinates of Jacobi vectors and Cartesian coordinates of the nuclei, using various free and constrained approaches depending on the way the rotation is treated. As a test calculation to be compared to full quantum dynamics it is applied here to the computation of the photodissociation spectrum of nitrosyl chloride (NOCl).

Rotranslational state-to-state rates and spectral representation of inelastic collisions in low-temperature molecular hydrogen

Inelastic collisions in natural H2 are studied from the experimental and theoretical points of view between 10 and 140 K. Rotational populations and number densities measured by Raman spectroscopy along supersonic expansions of H2 provide the link between experimental and theoretical rotranslational state-to-state rate coefficients of H2 in the vibrational ground state. These rates are calculated in the close-scattering approach with the MOLSCAT code employing a recent ab initio H2-H2 potential. The calculated rates are assessed by means of a master equation describing the time evolution of the experimental rotational populations. The feasibility for obtaining the rates on the sole basis of the experiment is discussed. The dominant processes j(1)j(2)-->j'(1)j'(2) in the investigated thermal range are found to be 21-->01 >30-->12 >31-->11, proving the importance of double processes such as 30-->12. Good agreement is found between theory and experiment, as well as with earlier ultrasonic measurements of relaxation times. A spectral representation is proposed in order to visualize quantitatively the collisional contributions in any nonequilibrium time evolving process.

State-to-state rotational transitions in H2+H2 collisions at low temperatures

We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the superthermal region. The two most recently published potential energy surfaces for the H(2)-H(2) complex, the so-called Diep-Johnson (DJ) [J. Chem. Phys. 112, 4465 (2000); 113, 3480 (2000)] and Boothroyd-Martin-Keogh-Peterson (BMKP) [J. Chem. Phys. 116, 666 (2002)] surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H(2)+H(2) collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H(2)-H(2) interaction. We found significant differences in rotational excitation/deexcitation cross sections computed on the two surfaces in collisions between two para-H(2) molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment B. Mate et al., [J. Chem. Phys. 122, 064313 (2005)] for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.