Publications by authors named "Marko T Cvitas"

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

The recently developed ring-polymer instanton (RPI) method [J. Chem. Phys.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

We extend a recently developed wave packet method for computing the state-to-state quantum dynamics of AB + CD → ABC + D reactions [M. T. Cvitaš and S.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF