Calculations of quantum tunnelling rates for muonium reactions with methane, ethane and propane

A single resonance Regge pole dominates the forward-angle scattering of the state-to-state F + H2 → FH + H reaction at Etrans = 62.09 meV

The aim of the present paper is to bring clarity, through simplicity, to the important and long-standing problem: does a resonance contribute to the forward-angle scattering of the F + H2 reaction? We reduce the problem to its essentials and present a well-defined, yet rigorous and unambiguous, investigation of structure in the differential cross sections (DCSs) of the following three state-to-state reactions at a translational energy of 62.09 meV: F + H2(vi = 0, ji = 0, mi = 0) → FH(vf = 3, jf = 0, 1, 2, mf = 0) + H, where vi, ji, mi and vf, jf, mf are the initial and final vibrational, rotational and helicity quantum numbers respectively. Firstly, we carry out quantum-scattering calculations for the Fu-Xu-Zhang potential energy surface, obtaining accurate numerical scattering matrix elements for indistinguishable H2. The calculations use a time-independent method, with hyperspherical coordinates and an enhanced Numerov method. Secondly, the following theoretical techniques are employed to analyse structures in the DCSs: (a) full and Nearside-Farside (NF) partial wave series (PWS) and local angular momentum theory, including resummations of the full PWS up to second order. (b) The recently introduced "CoroGlo" test, which lets us distinguish between glory and corona scattering at forward angles for a Legendre PWS. (c) Six asymptotic (semiclassical) forward-angle glory theories and three asymptotic farside rainbow theories, valid for rainbows at sideward-scattering angles. (d) Complex angular momentum (CAM) theories of forward and backward scattering, with the Regge pole positions and residues computed by Thiele rational interpolation. Thirdly, our conclusions for the three PWS DCSs are: (a) the forward-angle peaks arise from glory scattering. (b) A broad (hidden) farside rainbow is present at sideward angles. (c) A single Regge pole contributes to the DCS across the whole angular range, being most prominent at forward angles. This proves that a resonance contributes to the DCSs for the three transitions. (d) The diffraction oscillations in the DCSs arise from NF interference, in particular, interference between the Regge pole and direct subamplitudes.

Time evolution of ML-MCTDH wavefunctions. I. Gauge conditions, basis functions, and singularities

We derive a family of equations-of-motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than those used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOM approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.

Large Anharmonic Effects on Tunneling and Kinetics: Reaction of Propane with Muonium

Calculations of kinetic isotope effects (KIEs) provide challenging tests of quantal mass effects on reaction rates, and muonium KIEs are the most challenging. Here, we show that it can be very important to include reaction-coordinate-dependent vibrational anharmonicity along the whole reaction path to calculate tunneling probabilities and KIEs. For the reaction of propane with Mu, this decreases both the height and width of the vibrationally adiabatic potential barrier, with both effects increasing the rate constants. Our results agree well with the experimental observations.

Microcanonical Tunneling Rates from Density-of-States Instanton Theory

Semiclassical instanton theory is a form of quantum transition-state theory which can be applied to the computation of thermal reaction rates in complex molecular systems including quantum tunneling effects. There have been a number of attempts to extend the theory to treat microcanonical rates. However, the previous formulations are either computationally unfeasible for large systems due to an explicit sum over states or they involve extra approximations, which make them less reliable. We propose a robust and practical microcanonical formulation called density-of-states instanton theory, which avoids the sum over states altogether. In line with the semiclassical approximations inherent to the instanton approach, we employ the stationary-phase approximation to the inverse Laplace transform to obtain the densities of states. This can be evaluated using only post-processing of the data available from a small set of instanton calculations, such that our approach remains computationally efficient. We show that the new formulation predicts results that agree well with quantum scattering theory for an atom-diatom reaction and with experiments for a photoexcited unimolecular hydrogen transfer in a Criegee intermediate. When the thermal rate is evaluated from a Boltzmann average over our new microcanonical formalism, it can overcome some problems of conventional instanton theory. In particular, it predicts a smooth transition at the crossover temperature and is able to describe bimolecular reactions with pre-reactive complexes such as CH₃OH + OH.

Instanton theory of tunneling in molecules with asymmetric isotopic substitutions

We consider quantum tunneling in asymmetric double-well systems for which the local minima in the two wells have the same energy, but the frequencies differ slightly. In a molecular context, this situation can arise if the symmetry is broken by isotopic substitutions. We derive a generalization of instanton theory for these asymmetric systems, leading to a semiclassical expression for the tunneling matrix element and hence the energy-level splitting. We benchmark the method using a set of one- and two-dimensional models, for which the results compare favorably with numerically exact quantum calculations. Using the ring-polymer instanton approach, we apply the method to compute the level splittings in various isotopomers of malonaldehyde in full dimensionality and analyze the relative contributions from the zero-point energy difference and tunneling effects.