Quantum dynamics of non-rigid systems comprising two polyatomic fragments

Full-dimensional coupled-channel statistical approach to atom-triatom systems and applications to H/D + O3 reaction

The statistical quantum model (SQM), which assumes that the reactivity is controlled by entrance/exit channel quantum capture probabilities, is well suited for chemical reactions with a long-lived intermediate complex. In this work, a time-independent coupled-channel implementation of the SQM approach is developed for atom-triatom systems in full dimensionality. As SQM treats the capture dynamics quantum mechanically, it is capable of handling quantum effects such as tunneling. A detailed study of the H/D + O3 capture dynamics was performed by applying the newly developed SQM method on an accurate global potential energy surface. Agreement with previous ring polymer molecular dynamics (RPMD) results on the same potential energy surface is excellent except for very low temperatures. The SQM results are also in reasonably good agreement with available experimental rate coefficients. The strong H/D kinetic isotope effect underscores the dominant role of quantum tunneling under an entrance channel barrier at low temperatures.

Full quantum calculations of the line shape for H2O perturbed by Ar at temperatures from 20 to 300 K

This work theoretically studied the spectral line shape of H2O perturbed by Ar in the temperature range of 20-300 K for the pure rotational lines below 360 cm-1, as well as three lines (31, 2 ← 44, 1, 54, 2 ← 41, 3, and 73, 5 ← 60, 6) in the v2 band. In order to perform precise dynamical calculations at low collision energies, a full-dimensional long-range potential energy surface was constructed for the H2O-Ar system for the first time to correct the long range of our newly developed intermolecular potential energy surface. Subsequently, the six line-shape parameters (pressure-broadening and -shifting parameters, their speed dependencies, and the complex Dicke parameters) were determined from the generalized spectroscopic cross section by the full quantum time-independent close-coupling approach on this new potential energy surface. Our theoretical results are in good agreement with the available experimental observations. Furthermore, the influence of the speed-dependence and Dicke narrowing effects on the line contour was revealed by comparing the differences among the Hartmann-Tran, quadratic-speed-dependent Voigt, and Voigt profiles. The temperature dependence of each line-shape parameter was further parameterized using the triplet-power-law for three pure rotational 61, 6 ← 52, 3, 41, 4 ← 32, 1, and 31, 3 ← 22, 0 lines. These line-shape parameters will provide a comprehensive set of theoretical references for subsequent experimental measurements.

Methane dimer rovibrational states and Raman transition moments

Benchmark-quality rovibrational data are reported for the methane dimer from variational nuclear motion computations using an ab initio intermolecular potential energy surface reported by [M. P. Metz et al., Phys. Chem. Chem. Phys., 2019, 21, 13504-13525]. A simple polarizability model is used to compute Raman transition moments that may be relevant for future direct observation of the intermolecular dynamics. Non-negligible ΔK ≠ 0 transition moments arise in this symmetric top system due to strong rovibrational couplings.

A New Full-Dimensional Ab Initio Intermolecular Potential Energy Surface and Rovibrational Energies of the H2O-H2 Complex

H2O-H2 is a prototypical five-atom van der Waals system, and the interaction between H2O and H2 plays an important role in many physical and chemical environments. However, previous full-dimensional intermolecular potential energy surfaces (IPESs) cannot accurately describe the H2O-H2 interaction in the repulsive or van der Waals minimum region. In this work, we constructed a full-dimensional IPES for the title system with a small root-mean-square error of 0.252 cm-1 by using the permutation invariant polynomial neural network method. The ab initio calculations were performed by employing the explicitly corrected coupled cluster [CCSD(T)-F12a] method with the augmented correlation-consistent polarized valence quintuple-ζ basis set. Based on the newly developed IPES, the bound states of the H2O-H2 complex were calculated within the rigid-rotor approximation. The transition frequencies and band origins agreed well with the experimental values [Weida, M. J.; Nesbitt, D. J. J. Chem. Phys. 1999, 110, 156-167] with errors less than 0.1 cm-1 for most transitions. Those results demonstrate the high accuracy of our new IPES, which would build a solid foundation for the collisional dynamics of H2O-H2 at low temperatures.

Rovibrational states calculations of the H2O-HCN heterodimer with the multiconfiguration time dependent Hartree method

Water and hydrogen cyanide are two of the most common species in space and the atmosphere with the ability of binding to form dimers such as H2O-HCN. In the literature, while calculations characterizing various properties of the H2O-HCN cluster (equilibrium distance, vibrational frequencies and rotational constants) have been done in the past, extensive calculations of the rovibrational states of this system using a reliable quantum dynamical approach have yet to be reported. In this work, we intend to mend that by performing the first calculation of the rovibrational states of the H2O-HCN van der Waals complex on a recently developed potential energy surface. We use the block improved relaxation procedure implemented in the Heidelberg MultiConfiguration Time-Dependent Hartree (MCTDH) package to compute the states of the H2O-HCN isomer, from which we extract the transition frequencies and rotational constants of the complex. We further adapt an approach first suggested by Wang and Carrington-and supported here by analysis routines of the Heidelberg MCTDH package-to properly characterize the computed rovibrational states. The subsequent assignment of rovibrational states was done by theoretical analysis and visual inspection of the wavefunctions. Our simulations provide a Zero Point Energy (ZPE) and intermolecular vibrational frequencies in good agreement with past ab initio calculations. The transition frequencies and rotational constants obtained from our simulations match well with the available experimental data. This work has the broad aim to propose the MCTDH approach as a reliable option to compute and characterize rovibrational states of van der Waals complexes such as the current one.

H2O-HCN complex: A new potential energy surface and intermolecular rovibrational states from rigorous quantum calculations

In this work the H2O-HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O-HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O-HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O-HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former.

State-to-state rovibrational transition rates for CO2 in the bend mode in collisions with He atoms

Modeling environments that are not in local thermal equilibrium, such as protoplanetary disks or planetary atmospheres, with molecular spectroscopic data from space telescopes requires knowledge of the rate coefficients of rovibrationally inelastic molecular collisions. Here, we present such rate coefficients in a temperature range from 10 to 500 K for collisions of CO2 with He atoms in which CO2 is (de)excited in the bend mode. They are obtained from numerically exact coupled-channel (CC) calculations as well as from calculations with the less demanding coupled-states approximation (CSA) and the vibrational close-coupling rotational infinite-order sudden (VCC-IOS) method. All of the calculations are based on a newly calculated accurate ab initio four-dimensional CO2-He potential surface including the CO2 bend (ν2) mode. We find that the rovibrationally inelastic collision cross sections and rate coefficients from the CSA and VCC-IOS calculations agree to within 50% with the CC results at the rotational state-to-state level, except for the smaller ones and in the low energy resonance region, and to within 20% for the overall vibrational quenching rates except for temperatures below 50 K where resonances provide a substantial contribution. Our CC quenching rates agree with the most recent experimental data within the error bars. We also compared our results with data from Clary et al. calculated in the 1980s with the CSA [A. J. Banks and D. C. Clary, J. Chem. Phys. 86, 802 (1987)] and VCC-IOS [D. C. Clary, J. Chem. Phys. 78, 4915 (1983)] methods and a simple atom-atom model potential based on ab initio Hartree-Fock calculations and found that their cross sections agree fairly well with ours for collision energies above 500 cm-1, but that the inclusion of long range attractive dispersion interactions is crucial to obtain reliable cross sections at lower energies and rate coefficients at lower temperatures.

Computing excited OH stretch states of water dimer in 12D using contracted intermolecular and intramolecular basis functions

Due to the ubiquity and importance of water, water dimer has been intensively studied. Computing the (ro-)vibrational spectrum of water dimer is challenging. The potential has eight wells separated by low barriers, which makes harmonic approximations of limited utility. A variational approach is imperative, but difficult because there are 12 coupled vibrational coordinates. In this paper, we use a product contracted basis whose functions are products of intramolecular and intermolecular functions computed using an iterative eigensolver. An intermediate matrix F facilitates calculating matrix elements. Using F, it is possible to do calculations on a general potential without storing the potential on the full quadrature grid. We find that surprisingly many intermolecular functions are required. This is due to the importance of coupling between inter- and intra-molecular coordinates. The full G₁₆ symmetry of water dimer is exploited. We calculate, for the first time, monomer excited stretch states and compare P(1) transition frequencies with their experimental counterparts. We also compare with experimental vibrational shifts and tunneling splittings. Surprisingly, we find that the largest tunneling splitting, which does not involve the interchange of the two monomers, is smaller in the asymmetric stretch excited state than in the ground state. Differences between levels we compute and those obtained with a [6+6]D adiabatic approximation [Leforestier et al. J. Chem. Phys. 137 014305 (2012)] are ∼0.6 cm^-1 for states without monomer excitation, ∼4 cm^-1 for monomer excited bend states, and as large as ∼10 cm^-1 for monomer excited stretch states.

ABC+D: A time-independent coupled-channel quantum dynamics program for elastic and ro-vibrational inelastic scattering between atoms and triatomic molecules in full dimensionality

We discuss the details of a time-independent quantum mechanical method and its implementation for full-dimensional non-reactive scattering between a closed-shell triatomic molecule and a closed-shell atom. By solving the time-independent Schrödinger equation within the coupled-channel framework using a log-derivative method, the state-to-state scattering matrix (S-matrix) can be determined for inelastic scattering involving both the rotational and vibrational modes of the molecule. Various approximations are also implemented. The ABC+D code provides an important platform for understanding an array of physical phenomena involving collisions between atoms and molecules.

Adiabatic models for the quantum dynamics of surface scattering with lattice effects

In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.

Full-dimensional quantum studies of vibrational energy transfer dynamics between H2O and Ar: theory assessing experiment

We report the first full-dimensional quantum mechanical calculations of the ro-vibrational inelastic scattering dynamics between water molecules and argon atoms on an accurate potential energy surface, using a recently developed time-independent quantum method based on the close-coupling approach. The state-to-state integral cross-sections and rate coefficients show strong observance of gap laws. The calculated thermal rate coefficients for the relaxation of the stretching fundamental states of H₂O are in good agreement with experimental values, while those for the bending overtone state are approximately five times smaller than the values extracted through a previous kinetic modeling of fluorescence decay data. Our state-specific quantum scattering results suggest the need to reassess the kinetic modeling of the experimental data. This work advanced our understanding of the quantum dynamics of vibrationally inelastic energy transfer processes involving polyatomic molecules.

Reduced-dimensional vibrational models of the water dimer

A model based on the finite-basis representation of a vibrational Hamiltonian expressed in internal coordinates is developed. The model relies on a many-mode, low-order expansion of both the kinetic energy operator and the potential energy surface (PES). Polyad truncations and energy ceilings are used to control the size of the vibrational basis to facilitate accurate computations of the OH stretch and HOH bend intramolecular transitions of the water dimer (H₂ ¹⁶O)₂. Advantages and potential pitfalls of the applied approximations are highlighted. The importance of choices related to the treatment of the kinetic energy operator in reduced-dimensional calculations and the accuracy of different water dimer PESs are discussed. A range of different reduced-dimensional computations are performed to investigate the wavenumber shifts in the intramolecular transitions caused by the coupling between the intra- and intermolecular modes. With the use of symmetry, full 12-dimensional vibrational energy levels of the water dimer are calculated, predicting accurately the experimentally observed intramolecular fundamentals. It is found that one can also predict accurate intramolecular transition wavenumbers for the water dimer by combining a set of computationally inexpensive reduced-dimensional calculations, thereby guiding future effective-Hamiltonian treatments.

Vibrational Spectroscopy of the Water Dimer at Jet-Cooled and Atmospheric Temperatures

The vibrational spectroscopy of the water dimer provides an understanding of basic hydrogen bonding in water clusters, and with about one water dimer for every 1,000 water molecules, it plays a critical role in atmospheric science. Here, we review how the experimental and theoretical progress of the past decades has improved our understanding of water dimer vibrational spectroscopy under both cold and warm conditions. We focus on the intramolecular OH-stretching transitions of the donor unit, because these are the ones mostly affected by dimer formation and because their assignment has proven a challenge. We review cold experimental results from early matrix isolation to recent mass-selected jet expansion techniques and, in parallel, the improvements in the theoretical anharmonic models. We discuss and illustrate changes in the vibrational spectra of complexes upon increasing temperature, and the difficulties in recording and calculating these spectra. In the atmosphere, water dimer spectra at ambient temperature are crucial. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Noncovalently bound molecular complexes beyond diatom–diatom systems: full-dimensional, fully coupled quantum calculations of rovibrational states

The methodological advances made in recent years have significantly extended the range and dimensionality of noncovalently bound molecular complexes for which full-dimensional quantum calculations of their rovibrational states are feasible.

Femtochemistry of bimolecular reactions from weakly bound complexes: computational study of the H + H'OD → H'OH + D or HOD + H' exchange reactions

A full-dimensional wavepacket propagation describing the bimolecular exchange reactions H + H'OD → H'OH + D or HOD + H' initiated by photolysis of HCl in the hydrogen-bound complex (HCl)⋯(HOD) is reported. The dynamics of this reaction is carried out with the MCTDH method on an ab initio potential energy surface (PES) of H₃O and the initial state is derived from the ground state wavefunction of the complex obtained by relaxation on its own electronic ground state ab initio PES. The description of the system makes use of polyspherical coordinates parametrizing a set of Radau and Jacobi vectors. The calculated energy- and time-resolved reaction probabilities show, owing to the large collision energies at play stemming from the (almost full) photolysis of HCl, that the repulsion between oxygen in the H'OD molecule and the incoming hydrogen atom is the main feature of the collision and leads to non-reactive scattering. No abstraction reaction products are observed. However, both exchange processes are still observable, with a preference in O-H' bond dissociation over that of O-D. The selectivity is reversed upon vibrational pre-excitation of the O-D stretching mode in the H'OD molecule. It is shown that, after the collision, the hydrogen atom of HCl does most likely not encounter the almost stationary chlorine atom again but we also consider the limit case where the H atom is forced to collide multiple times against H'OD as a result of being pushed back by the Cl atom.

A Time-Independent Quantum Approach to Ro-vibrationally Inelastic Scattering between Atoms and Triatomic Molecules

A full-dimensional time-independent quantum mechanical theory for ro-vibrationally inelastic scattering of triatomic molecules with atoms is formulated. The Jacobi-Radau coordinate system used in the calculation allows not only a near perfect description of the vibrational problem but also the adaptation of the exchange symmetry for A₂B type triatoms. The S-matrix elements are obtained by solving the close-coupling equations with contracted basis using the log-derivative method. This method is applied to the inelastic scattering of the water molecule by a chlorine atom, which sheds light on the energy gap law in energy transfer in atom-triatom collisions.

DCl-H2O, HCl-D2O, and DCl-D2O Dimers: Inter- and Intramolecular Vibrational States and Frequency Shifts from Fully Coupled Quantum Calculations on a Full-Dimensional Neural Network Potential Energy Surface

We report full-dimensional and fully coupled quantum calculations of the inter- and intramolecular vibrational states of three isotopologues of the hydrogen chloride-water dimer: DCl-H₂O (DH), HCl-D₂O (HD), and DCl-D₂O (DD). The present study extends our recent theoretical investigation of the nine-dimensional (9D) vibrational level structure of the HCl-H₂O (HH) dimer [Liu, Y.; Li, J.; Felker, P. M.; Bačić, Z. Phys. Chem. Chem. Phys. 2021, 23, 7101-7114]. It employs the same accurate 9D permutation invariant polynomial-neural network potential energy surface and the highly efficient bound-state methodology. The objective of this work is to elucidate the isotopologue variations of a range of bound-state properties of the hydrogen chloride-water dimer and compare them to those of the HH dimer. In order to achieve this, for the isotopologues considered, the rigorous 9D quantum calculations performed encompass all intramolecular vibrational fundamentals, and their frequency shifts relative to the isolated monomer values, together with the low-lying intermolecular vibrational states in each of the intramolecular vibrational manifolds of interest. Moreover, for the ground state of each isotopologue, several informative vibrationally averaged intermolecular geometric properties of the dimer are computed, as well as the three rotational constants. The energies of the intermolecular inversion and rock modes, which mainly involve the motions of the water moiety, differ greatly for H₂O and D₂O, but are much less sensitive to whether the hydrogen chloride isotopologue is HCl or DCl. On the other hand, the excitation of the HCl/DCl stretch changes significantly the energies of the water inversion and rock modes. The DCl stretch frequency shift computed in 9D for the DD dimer, -114.91 cm^-1, agrees extremely well with the corresponding experimental value of -115.20 cm^-1 measured by Saykally and co-workers.

Performance of a black-box-type rovibrational method in comparison with a tailor-made approach: Case study for the methane-water dimer

The present work intends to join and respond to the excellent and thoroughly documented rovibrational study of X. G. Wang and T. Carrington, Jr. [J. Chem. Phys. 154, 124112 (2021)] that used an approach tailored for floppy dimers with an analytic dimer Hamiltonian and a non-product basis set including Wigner D functions. It is shown in the present work that the GENIUSH black-box-type rovibrational method can approach the performance of the tailor-made computation for the example of the floppy methane-water dimer. Rovibrational transition energies and intensities are obtained in the black-box-type computation with a twice as large basis set and in excellent numerical agreement in comparison with the more efficient tailor-made approach.

Using nondirect product Wigner D basis functions and the symmetry-adapted Lanczos algorithm to compute the ro-vibrational spectrum of CH4-H2O

By doing calculations on the methane-water van der Waals complex, we demonstrate that highly converged energy levels and wavefunctions can be obtained using Wigner D basis functions and the Symmetry-Adapted Lanczos (SAL) method. The Wigner D basis is a nondirect product basis and, therefore, efficient when the kinetic energy operator has accessible singularities. The SAL method makes it possible to exploit symmetry to label energy levels and reduce the cost of the calculation, without explicitly using symmetry-adapted basis functions. Line strengths are computed, and new bands are identified. In particular, we find unusually strong transitions between states associated with the isomers of the global minimum and the secondary minimum.

Full-dimensional quantum dynamics study of isotope effects for the H2 + NH2/ND2/NHD and H2/D2/HD + NH2 reactions

A full-dimensional quantum dynamical study for the bimolecular reactions of hydrogen molecules with amino radicals for different isotopologues is reported. The nonreactive amino radical is described by two Radau vectors that are very close to the valence bond coordinates. Potential-optimized discrete variable representation basis is used for the vibrational coordinates of the amino radical. Starting from the reaction H₂ + NH₂, we study the isotope effects for the two reagents separately, i.e., H₂ + NH₂/ND₂/NHD and H₂/D₂/HD + NH₂. The effects of different vibrational mode excitations of the reagents on the reactivities are studied. Physical explanations about the isotope effects are also provided thoroughly including the influence of vibrational energy differences between the different isotopologues and the impact of the tunneling effect.

HCl–H2O dimer: an accurate full-dimensional potential energy surface and fully coupled quantum calculations of intra- and intermolecular vibrational states and frequency shifts

The present work reports a new full-dimensional potential energy surface (PES) of the HCl–H₂O dimer, and the first fully coupled 9D quantum calculations of the intra- and intermolecular vibrational states of the complex, utilizing this PES.

Full-dimensional quantum mechanical calculations of the reaction probability of the H + CH4 reaction based on a mixed Jacobi and Radau description

A full-dimensional time-dependent wave packet study using mixed polyspherical Jacobi and Radau coordinates for the title reaction has been reported. The non-reactive moiety CH₃ has been described using three Radau vectors, whereas two Jacobi vectors have been used for the bond breaking/formation process. A potential-optimized discrete variable representation basis has been employed to describe the vibrational coordinates of the reagent CH₄. About one hundred billion basis functions have been necessary to achieve converged results. The reaction probabilities for some initial vibrational states are given. A comparison between the present approach and other methods, including reduced and full-dimensional ones, is also presented.

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.

General mathematical formulation of scattering processes in atom-diatomic collisions in the RmatReact methodology

Accurately modelling cold and ultracold reactive collisions occurring over deep potential wells, such as [Formula: see text], requires the development of new theoretical and computational methodologies. One potentially useful framework is the R-matrix method adopted widely for electron-molecule collisions which has more recently been applied to non-reactive heavy-particle collisions such as Ar-Ar. The existing treatment of non-reactive elastic and inelastic scattering needs to be substantially extended to enable modelling of reactive collisions: this is the subject of this paper. Herein, we develop the general mathematical formulation for non-reactive elastic and inelastic scattering, photoassociation, photodissociation, charge exchange and reactive scattering using the R-matrix method. Of particular note is that the inner region, of central importance to calculable R-matrix methodologies, must be finite in all scattering coordinates rather than a single scattering coordinate as for non-reactive scattering. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H₃⁺, H₅⁺ and beyond'.

Computational study of the ro-vibrational spectrum of CO-CO2

An accurate ab initio ground-state intermolecular potential energy surface (PES) was determined for the CO-CO₂ van der Waals dimer. The Lanczos algorithm was used to compute rovibrational energies on this PES. For both the C-in and O-in T-shaped isomers, the fundamental transition frequencies agree well with previous experimental results. We confirm that the in-plane states previously observed are geared states. In addition, we have computed and assigned many other vibrational states. The rotational constants we determine from J = 1 energy levels agree well with their experimental counterparts. Planar and out-of-plane cuts of some of the wavefunctions we compute are quite different, indicating strong coupling between the bend and torsional modes. Because the stable isomers are T-shaped, vibration along the out-of-plane coordinates is very floppy. In CO-CO₂, when the molecule is out-of-plane, interconversion of the isomers is possible, but the barrier height is higher than the in-plane geared barrier height.

Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements

The rovibrational energy levels and intensities of the CH₃F-H₂ dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H₂O-H₂ dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H₂ are more strongly bound than dimers containing para-H₂. Energies and wavefunctions of the discrete rovibrational levels of CH₃F-paraH₂ complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH₃F-paraH₂ complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH₃F-orthoH₂ dimer has been predicted for the first time.

Computing vibration-rotation-tunnelling levels of HOD dimer

Using an accurate 6D water dimer potential energy surface, we compute vibration-rotation-tunnelling levels of HOD dimer, by assuming that the two monomers are rigid. HOD dimer has two isomers, a D-bonded isomer and an H-bonded isomer, and the wavefunctions of both isomers have amplitude in four wells. HOD dimer is important because, unlike the case of H2O dimer or D2O dimer, it is possible to measure the largest tunnelling splitting. Results for HOD dimer, therefore facilitate the testing of H2O dimer potentials. In J. Chem. Phys., 1995, 102, 1114, experimental results were interpreted in terms of 1D models. Experimental splittings of both isomers, obtained by fitting an energy level equation to spectra, are in good agreement with those we compute.

Analytic Morse/long-range potential energy surfaces and "adiabatic-hindered-rotor" treatment for a symmetric top-linear molecule dimer: A case study of CH3F-H2

A first effective six-dimensional ab initio potential energy surface (PES) for CH₃F-H₂ which explicitly includes the intramolecular Q₃ stretching normal mode of the CH₃F monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster level of theory [CCSD(T)-F12a] with an augmented correlation-consistent triple zeta basis set. Five-dimensional analytical intermolecular PESs for ν₃(CH₃F) = 0 and 1 are then obtained by fitting the vibrationally averaged potentials to the Morse/Long-Range (MLR) potential function form. The MLR function form is applied to the nonlinear molecule-linear molecule case for the first time. These fits to 25 015 points have root-mean-square deviations of 0.74 cm^-1 and 0.082 cm^-1 for interaction energies less than 0.0 cm^-1. Using the adiabatic hindered-rotor approximation, three-dimensional PESs for CH₃F-paraH₂ are generated from the 5D PESs over all possible orientations of the hydrogen monomer. The infrared and microwave spectra for CH₃F-paraH₂ dimer are predicted for the first time. These analytic PESs can be used for modeling the dynamical behavior in CH₃F-(H₂)_N clusters, including the possible appearance of microscopic superfluidity.

Rovibrational quantum dynamical computations for deuterated isotopologues of the methane–water dimer

Rovibrational states of methane–water isotopologues are computed in a variational procedure and the wave functions are analyzed in terms of the rigid-rotor and coupled-rotors models.

Classical calculation of the equilibrium constants for true bound dimers using complete potential energy surface

The present paper aims at deriving classical expressions which permit calculation of the equilibrium constant for weakly interacting molecular pairs using a complete multidimensional potential energy surface. The latter is often available nowadays as a result of the more and more sophisticated and accurate ab initio calculations. The water dimer formation is considered as an example. It is shown that even in case of a rather strongly bound dimer the suggested expression permits obtaining quite reliable estimate for the equilibrium constant. The reliability of our obtained water dimer equilibrium constant is briefly discussed by comparison with the available data based on experimental observations, quantum calculations, and the use of RRHO approximation, provided the latter is restricted to formation of true bound states only.

The vibration-rotation-tunneling levels of N2-H2O and N2-D2O

In this paper, we report vibration-rotation-tunneling levels of the van der Waals clusters N2-H2O and N2-D2O computed from an ab initio potential energy surface. The only dynamical approximation is that the monomers are rigid. We use a symmetry adapted Lanczos algorithm and an uncoupled product basis set. The pattern of the cluster's levels is complicated by splittings caused by H-H exchange tunneling (larger splitting) and N-N exchange tunneling (smaller splitting). An interesting result that emerges from our calculation is that whereas in N2-H2O, the symmetric H-H tunnelling state is below the anti-symmetric H-H tunnelling state for both K = 0 and K = 1, the order is reversed in N2-D2O for K = 1. The only experimental splitting measurements are the D-D exchange tunneling splittings reported by Zhu et al. [J. Chem. Phys. 139, 214309 (2013)] for N2-D2O in the v2 = 1 region of D2O. Due to the inverted order of the split levels, they measure the sum of the K = 0 and K = 1 tunneling splittings, which is in excellent agreement with our calculated result. Other splittings we predict, in particular those of N2-H2O, may guide future experiments.