Tunneling splittings using modified WKB method in Cartesian coordinates: The test case of vinyl radical.

Modified WKB theory for calculating tunneling splittings in symmetric multi-well systems in full dimensionality is re-derived using Cartesian coordinates. It is explicitly shown that the theory rests on the wavefunction that is exact for harmonic potentials. The theory was applied to calculate tunneling splittings in vinyl radical and some of its deuterated isotopologues in their vibrational ground states and the low-lying vibrationally excited states and compared to exact variational results. The exact results are reproduced within a factor of 2 in most states. Remarkably, all large enhancements of tunneling splittings relative to the ground state, up to three orders in magnitude in some excited mode combinations, are well reproduced. It is also shown that in the asymmetrically deuterated vinyl radical, the theory correctly predicts the states that are localized in a single well and the delocalized tunneling states. Modified WKB theory on the minimum action path is computationally inexpensive and can also be applied without modification to much larger systems in full dimensionality; the results of this test case serve to give insight into the expected accuracy of the method.

Tunneling splittings in the vibrationally excited states of water trimer.

Tunneling splitting (TS) patterns in vibrationally excited states of the water trimer are calculated, taking into account six tunneling pathways that describe the flips of free OH bonds and five bifurcation mechanisms that break and reform hydrogen bonds in the trimer ring. A tunneling matrix (TM) model is used to derive the energy shifts due to tunneling in terms of the six distinct TM elements in symbolic form. TM elements are calculated using the recently-developed modified WKB (Wentzel-Kramers-Brillouin) method in full dimensionality. Convergence was achieved for the lowest six excited vibrational modes. Bifurcation widths of the pseudorotational quartets are shown to be of comparable size to the ground-state widths, obtained using instanton theory, or increased for some particular modes of vibration. The largest increase is obtained for the excited out-of-phase flip of two adjacent water monomers with free OH bonds pointing in opposite directions relative to the ring plane. Bifurcation widths in (D2O)3 are found to be two orders of magnitude smaller than in (H2O)3. Geometrical arguments were used to explain the order of states in some TS multiplets in vibrationally excited water trimers.

Vibrational Tunneling Spectra of Molecules with Asymmetric Wells: A Combined Vibrational Configuration Interaction and Instanton Approach.

A combined approach that uses the vibrational configuration interaction (VCI) and semiclassical instanton theory was developed to study vibrational tunneling spectra of molecules with multiple wells in full dimensionality. The method can be applied to calculate low-lying vibrational states in the systems with an arbitrary number of minima, which are not necessarily equal in energy or shape. It was tested on a two-dimensional double-well model system and on malonaldehyde, and the calculations reproduced the exact quantum mechanical (QM) results with high accuracy. The method was subsequently applied to calculate the vibrational spectrum of the asymmetrically deuterated malonaldehyde with nondegenerate vibrational frequencies in the two wells. The spectrum is obtained at a cost of single-well VCI calculations used to calculate the local energies. The interactions between states of different wells are computed semiclassically using the instanton theory at a comparatively negligible computational cost. The method is particularly suited to systems in which the wells are separated by large potential barriers and tunneling splittings are small, for example, in some water clusters, when the exact QM methods come at a prohibitive computational cost.

Tunnelling splitting patterns in some partially deuterated water trimers.

We apply our recently developed semiclassical method for calculating tunnelling splittings (TS) in asymmetric systems to make the first characterization of the ground-state TS pattern of some partially deuterated water trimers. Similarly to homoisotopic water trimers, the ground-state TS patterns are explained in terms of six distinct rearrangement mechanisms. TS patterns in (D2O)(H2O)2 and (H2O)(D2O)2 are composed of sextets induced by the dynamics of flips, and each of its levels is further finely split into a quartet of doublets and a doublet of quartets, respectively, due to various bifurcation dynamics. The TS pattern is obtained using 18 distinct tunnelling matrix elements. TS patterns of (HOD)(H2O)2 and (HOD)(D2O)2 each consists of two sextets, belonging to in-bond and out-of-bond substituted isomers. These sextet levels are further split into quartets by bifurcations. The TS pattern is computed in terms of 13 matrix elements. We also derive analytic expressions for bifurcation tunnelling splittings in terms of tunnelling matrix elements using symmetry. The present approach can be applied to other water clusters and also to the low-lying vibrationally excited states and should help in the interpretation and assignment of experimental spectra in the future.

Vibrationally resolved valence and core photoionization and photoexcitation spectra of an electron-deficient trivalent boron compound: the case of catecholborane.

Compounds containing trivalent boron (TB) as the electron-deficient site(s) find numerous practical uses ranging from Lewis bases in organic synthesis to high-tech industry, with a number of novel applications anticipated. We present an experimental and theoretical study of the gas-phase valence photoionization (VUV-PES), core photoionization (XPS) and photoexcitation (NEXAFS) spectra of a representative TB compound catecholborane (CB). For modelling and assigning the spectra we used the ΔDFT and restricted single excitation space TD-DFT methods for the XPS and NEXAFS, and OVGF and EOM-CCSD for the VUV-PES. The vibrationally resolved structure was computed in the Franck-Condon (FC) and Herzberg-Teller (FCHT) approximations generally resulting in a good agreement with the observed spectral features. For the prediction of core-electron binding energies (CEBEs) several density functionals were tested. The best performance overall was furnished by ωB97X-D suggesting that including the dispersion correction is beneficial. The FCHT vibronic intensities are in clear discrepancy with the B 1s NEXAFS spectrum if the harmonic approximation is used for the B-H wagging mode both in the ground and in the first core-excited state. Instead, a much better agreement is obtained if the excited state potential is approximated to a symmetric double-well. The observed vibronic pattern could be a general fingerprint of the presence of TB centre(s), specifically, the transfer of the (core) density to the vacant boron p-orbital in the excited state.

Tunneling splittings of vibrationally excited states using general instanton paths.

A multidimensional semiclassical method for calculating tunneling splittings in vibrationally excited states of molecules using Cartesian coordinates is developed. It is an extension of the theory by Mil'nikov and Nakamura [J. Chem. Phys. 122, 124311 (2005)] to asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. Additionally, new terms are introduced in the description of the semiclassical wavefunction that drastically improves the splitting estimates for certain systems. The method is based on the instanton theory and builds the semiclassical wavefunction of the vibrationally excited states from the ground-state instanton wavefunction along the minimum action path and its harmonic neighborhood. The splittings of excited states are thus obtained at a negligible added numerical effort. The cost is concentrated, as for the ground-state splittings, in the instanton path optimization and the hessian evaluation along the path. The method can thus be applied without modification to many mid-sized molecules in full dimensionality and in combination with on-the-fly evaluation of electronic potentials. The tests were performed on several model potentials and on the water dimer.

Instanton theory of ground-state tunneling splittings with general paths.

We derive a multidimensional instanton theory for calculating ground-state tunneling splittings in Cartesian coordinates for general paths. It is an extension of the method by Mil'nikov and Nakamura [J. Chem. Phys. 115, 6881 (2001)] to include asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. The approach avoids multiple expensive matrix diagonalizations to converge the fluctuation prefactor in the ring-polymer instanton (RPI) method, and instead replaces them by an integration of a Riccati differential equation. When combined with the string method for locating instantons, we avoid the need to converge the calculation with respect to the imaginary time period of the semiclassical orbit, thereby reducing the number of convergence parameters of the optimized object to just one: the number of equally spaced system replicas used to represent the instanton path. The entirety of the numerical effort is thus concentrated in optimizing the shape of the path and evaluating hessians along the path, which is a dramatic improvement over RPI. In addition to the standard instanton approximations, we neglect the coupling of vibrational modes to external rotations. The method is tested on the model potential of malonaldehyde and on the water dimer and trimer, giving close agreement with RPI at a much-reduced cost.

Quantum tunnelling pathways of the water pentamer.

We apply the semiclassical instanton method to calculate all feasible tunnelling pathways in the water pentamer. Similarly to the water trimer, there are labile flip dynamics as well as a number of different bifurcation pathways coupled to flips. In contrast to the trimer, the puckering motion of the oxygen ring makes the ring-polymer instanton approach difficult to converge, a problem which is resolved by using a recently developed time-independent formalism of the method. We use the results to predict the complete ground-state tunnelling splitting pattern of 320 states, which should help in the continuing effort to assign the experimental spectrum. A comparison between the rearrangement pathways in the water trimer and pentamer sheds light on the many-body cooperative effects of hydrogen bonding which are important for a full understanding of the liquid state.

Rotation-tunneling spectrum of the water dimer from instanton theory.

We present an extension to the ring polymer instanton (RPI) method that includes overall rotations in the tunneling pathways, allowing a calculation of the full rotation-tunneling spectrum of the water dimer. Our method gives a qualitative description of the rotational and tunneling processes underpinning the spectrum, and shows the drastic reduction of the largest splitting with increasing rotational excitation. We show that this reduction is due to the strong coupling between the rotational motion about the principal axis and the acceptor tunneling motion (where the acceptor monomer rotates about the hydrogen bond), which results in a large increase in the path length and hence the quantum action. Our method gives a clear physical understanding of the behaviour of tunneling splittings in the presence of rotations, and scales linearly with system size.

Quadratic String Method for Locating Instantons in Tunneling Splitting Calculations.

The ring-polymer instanton (RPI) method is an efficient technique for calculating approximate tunneling splittings in high-dimensional molecular systems. In the RPI method, tunneling splitting is evaluated from the properties of the minimum action path (MAP) connecting the symmetric wells, whereby the extensive sampling of the full potential energy surface of the exact quantum-dynamics methods is avoided. Nevertheless, the search for the MAP is usually the most time-consuming step in the standard numerical procedures. Recently, nudged elastic band (NEB) and string methods, originaly developed for locating minimum energy paths (MEPs), were adapted for the purpose of MAP finding with great efficiency gains [ J. Chem. Theory Comput. 2016 , 12 , 787 ]. In this work, we develop a new quadratic string method for locating instantons. The Euclidean action is minimized by propagating the initial guess (a path connecting two wells) over the quadratic potential energy surface approximated by means of updated Hessians. This allows the algorithm to take many minimization steps between the potential/gradient calls with further reductions in the computational effort, exploiting the smoothness of potential energy surface. The approach is general, as it uses Cartesian coordinates, and widely applicable, with computational effort of finding the instanton usually lower than that of determining the MEP. It can be combined with expensive potential energy surfaces or on-the-fly electronic-structure methods to explore a wide variety of molecular systems.

Locating Instantons in Calculations of Tunneling Splittings: The Test Case of Malonaldehyde.

The recently developed ring-polymer instanton (RPI) method [J. Chem. Phys. 2011, 134, 054109] is an efficient technique for calculating approximate tunneling splittings in high-dimensional molecular systems. The key step is locating the instanton tunneling-path at zero temperature. Here, we show that techniques previously designed for locating instantons in finite-temperature rate calculations can be adapted to the RPI method, where they become extremely efficient, reducing the number of potential energy calls by 2 orders of magnitude. We investigate one technique that employs variable time steps to minimize the action integral, and two that employ equally spaced position steps to minimize the abbreviated (i.e., Jacobi) action integral, using respectively the nudged elastic band (NEB) and string methods. We recommend use of the latter because it is parameter-free, but all three methods give comparable efficiency savings. Having located the instanton pathway, we then interpolate the instanton path onto a fine grid of imaginary time points, allowing us to compute the fluctuation prefactor. The crucial modification needed to the original finite-temperature algorithms is to allow the end points of the zero-temperature instanton path to describe overall rotations, which is done using a standard quaternion algorithm. These approaches will allow the RPI method to be combined effectively with expensive potential energy surfaces or on-the-fly electronic structure methods.

Shallow-tunnelling correction factor for use with Wigner-Eyring transition-state theory.

We obtain a shallow-tunnelling correction factor for use with Wigner-Eyring transition-state theory (TST). Our starting point is quantum transition state theory (QTST), which approximates the accurate quantum rate as the instantaneous flux through a delocalised transition-state ensemble of ring-polymers. Expanding the ring-polymer potential to second order gives the well-known Wigner tunnelling-factor which diverges at the cross-over temperature between deep and shallow tunnelling. Here, we show how to remove this divergence by integrating numerically over the two softest ring-polymer normal modes. This results in a modified Wigner correction factor involving a one-dimensional integral evaluated along a straight line on the potential energy surface. Comparisons with accurate quantum calculations indicate that the newly derived correction factor gives realistic estimates of quantum rate coefficients in the shallow-tunnelling regime.

Which Is Better at Predicting Quantum-Tunneling Rates: Quantum Transition-State Theory or Free-Energy Instanton Theory?

Quantum transition-state theory (QTST) and free-energy instanton theory (FEIT) are two closely related methods for estimating the quantum rate coefficient from the free-energy at the reaction barrier. In calculations on one-dimensional models, FEIT typically gives closer agreement than QTST with the exact quantum results at all temperatures below the crossover to deep tunneling, suggesting that FEIT is a better approximation than QTST in this regime. Here we show that this simple trend does not hold for systems of greater dimensionality. We report tests on several collinear and three-dimensional reactions, in which QTST outperforms FEIT over a range of temperatures below crossover, which can extend down to half the crossover temperature (below which FEIT outperforms QTST). This suggests that QTST-based methods such as ring-polymer molecular dynamics (RPMD) may often give closer agreement with the exact quantum results than FEIT.

A Chebyshev method for state-to-state reactive scattering using reactant-product decoupling: OH + H2 → H2O + H.

We extend a recently developed wave packet method for computing the state-to-state quantum dynamics of AB + CD → ABC + D reactions [M. T. Cvitas and S. C. Althorpe, J. Phys. Chem. A 113, 4557 (2009)] to include the Chebyshev propagator. The method uses the further partitioned approach to reactant-product decoupling, which uses artificial decoupling potentials to partition the coordinate space of the reaction into separate reactant, product, and transition-state regions. Separate coordinates and basis sets can then be used that are best adapted to each region. We derive improved Chebyshev partitioning formulas which include Mandelshtam-and-Taylor-type decoupling potentials, and which are essential for the non-unitary discrete variable representations that must be used in 4-atom reactive scattering calculations. Numerical tests on the fully dimensional OH + H2 → H2O + H reaction for J = 0 show that the new version of the method is as efficient as the previously developed split-operator version. The advantages of the Chebyshev propagator (most notably the ease of parallelization for J > 0) can now be fully exploited in state-to-state reactive scattering calculations on 4-atom reactions.

State-to-state reactive scattering in six dimensions using reactant-product decoupling: OH + H2 → H2O + H (J = 0).

We extend to full dimensionality a recently developed wave packet method [M. T. Cvitas and S. C. Althorpe, J. Phys. Chem. A 113, 4557 (2009)] for computing the state-to-state quantum dynamics of AB + CD → ABC + D reactions and also increase the computational efficiency of the method. This is done by introducing a new set of product coordinates, by applying the Crank-Nicholson approximation to the angular kinetic energy part of the split-operator propagator and by using a symmetry-adapted basis-to-grid transformation to evaluate integrals over the potential energy surface. The newly extended method is tested on the benchmark OH + H(2) → H(2)O + H reaction, where it allows us to obtain accurately converged state-to-state reaction probabilities (on the Wu-Schatz-Fang-Lendvay-Harding potential energy surface) with modest computational effort. These methodological advances will make possible efficient calculations of state-to-state differential cross sections on this system in the near future.

Quantum wave packet method for state-to-state reactive scattering calculations on AB + CD --> ABC + D reactions.

We describe a quantum wave packet method for computing the state-to-state quantum dynamics of 4-atom AB + CD --> ABC + D reactions. The approach is an extension to 4-atom reactions of a version of the reactant-product decoupling (RPD) approach, applied previously to 3-atom reactions ( J. Chem. Phys. 2001, 114 , 1601 ). The approach partitions the coordinate space of the reaction into separate reagent, strong-interaction, and product regions, using a system of artificial absorbing and reflecting potentials. It employs a partitioned version of the split-operator propagator, which is more efficient than partitioning the (exact) time-dependent Schrodinger equation. The wave packet bounces off a reflecting potential in the entrance channel, which generates a source term; this is transformed efficiently from reagent to product Jacobi coordinates by exploiting some simple angular momentum properties. The efficiency and accuracy of the method is demonstrated by numerical tests on the benchmark OH + H(2) --> H(2)O + H reaction.

Interactions and dynamics in Li+Li2 ultracold collisions.

A potential energy surface for the lowest quartet electronic state ((4)A(')) of lithium trimer is developed and used to study spin-polarized Li+Li(2) collisions at ultralow kinetic energies. The potential energy surface allows barrierless atom exchange reactions. Elastic and inelastic cross sections are calculated for collisions involving a variety of rovibrational states of Li(2). Inelastic collisions are responsible for trap loss in molecule production experiments. Isotope effects and the sensitivity of the results to details of the potential energy surface are investigated. It is found that for vibrationally excited states, the cross sections are only quite weakly dependent on details of the potential energy surface.

Ultracold collisions involving heteronuclear alkali metal dimers.

We carry out the first quantum dynamics calculations on ultracold atom-diatom collisions in isotopic mixtures. The systems studied are spin-polarized 7Li + 6Li7Li, 7Li + 6Li2, 6Li + 6Li7Li, and 6Li + 7Li2. Reactive scattering can occur for the first two systems even when the molecules are in their ground rovibrational states, but is slower than vibrational relaxation in homonuclear systems. Implications for sympathetic cooling of heteronuclear molecules are discussed.

Ultracold Li + Li2 collisions: bosonic and fermionic cases.

We have carried out quantum dynamical calculations of vibrational quenching in Li + Li(2) collisions for both bosonic (7)Li and fermionic (6)Li. These are the first ever such calculations involving fermionic atoms. We find that for the low initial vibrational states considered here (v < or = 3), the quenching rates are not suppressed for fermionic atoms. This contrasts with the situation found experimentally for molecules formed via Feshbach resonances in very high vibrational states.

Quantum dynamics of ultracold Na+ Na2 collisions.

Ultracold collisions between spin-polarized Na atoms and vibrationally excited Na2 molecules are investigated theoretically, using a reactive scattering formalism (including atom exchange). Calculations are carried out on both pairwise additive and nonadditive potential energy surfaces for the quartet electronic state. The Wigner threshold laws are followed for energies below 10(-5) K. Vibrational relaxation processes dominate elastic processes for temperatures below 10(-3)-10(-4) K. For temperatures below 10(-5) K, the rate coefficients for vibrational relaxation (v=1-->0) are 4.8x10(-11) and 5.2x10(-10) cm(3) s(-1) for the additive and nonadditive potentials, respectively. The large difference emphasizes the importance of using accurate potential energy surfaces for such calculations.