Derivation of a true (t->0+) quantum transition-state theory. II. Recovery of the exact quantum rate in the absence of recrossing

In Part I [J. Chem. Phys. 138, 084108 (2013)] we derived a quantum transition-state theory by taking the t->0+ (short-time) limit of a new form of quantum flux-side time-correlation function containing a ring-polymer dividing surface. This t->0+ limit appears to be unique in giving positive-definite Boltzmann statistics, and is identical to ring-polymer molecular dynamics (RPMD) TST. Here, we show that quantum TST (i.e. RPMD-TST) is exact if there is no recrossing (by the real-time quantum dynamics) of the ring-polymer dividing surface, nor of any surface orthogonal to it in the space describing fluctuations in the polymer-bead positions along the reaction coordinate. In practice, this means that RPMD-TST gives a good approximation to the exact quantum rate for direct reactions, provided the temperature is not too far below the cross-over to deep tunnelling. We derive these results by comparing the long-time limit of the ring-polymer flux-side time-correlation function with that of a hybrid flux-side time-correlation function (containing a ring-polymer flux operator and a Miller-Schwarz-Tromp side function), and by representing the resulting ring-polymer momentum integrals as hypercubes. Together with Part I, the results of this article validate a large number of RPMD calculations of reaction rates.Comment: 14 pages, 4 figures. Argument and wording clarified, typographical errors corrected and references adde

On the uniqueness of t->0+ quantum transition-state theory

It was shown recently that there exists a true quantum transition-state theory (QTST) corresponding to the t->0+ limit of a (new form of) quantum flux-side time-correlation function. Remarkably, this QTST is identical to ring-polymer molecular dynamics (RPMD) TST. Here we provide evidence which suggests very strongly that this QTST (= RPMD-TST) is unique, in the sense that the t->0+ limit of any other flux-side time-correlation function gives either non-positive-definite quantum statistics or zero. We introduce a generalized flux-side time-correlation function which includes all other (known) flux-side time-correlation functions as special limiting cases. We find that the only non-zero t->0+ limit of this function that contains positive-definite quantum statistics is RPMD-TST.Comment: 10 pages, 1 figure. Typographical errors corrected, references updated and adde

Calculating splittings between energy levels of different symmetry using path-integral methods.

It is well known that path-integral methods can be used to calculate the energy splitting between the ground and the first excited state. Here we show that this approach can be generalized to give the splitting patterns between all the lowest energy levels from different symmetry blocks that lie below the first-excited totally symmetric state. We demonstrate this property numerically for some two-dimensional models. The approach is likely to be useful for computing rovibrational energy levels and tunnelling splittings in floppy molecules and gas-phase clusters.E.M. and S.C.A. acknowledge funding from the UK Engineering and Physical Sciences Research Council.This is the author accepted manuscript. The final version is available from the American Institute of Physics via http://dx.doi.org/10.1063/1.494398

The method of Gaussian weighted trajectories. V. On the 1GB procedure for polyatomic processes

In recent years, many chemical reactions have been studied by means of the quasi-classical trajectory (QCT) method within the Gaussian binning (GB) procedure. The latter consists in "quantizing" the final vibrational actions in Bohr spirit by putting strong emphasis on the trajectories reaching the products with vibrational actions close to integer values. A major drawback of this procedure is that if N is the number of product vibrational modes, the amount of trajectories necessary to converge the calculations is ~ 10^N larger than with the standard QCT method. Applying it to polyatomic processes is thus problematic. In a recent paper, however, Czako and Bowman propose to quantize the total vibrational energy instead of the vibrational actions [G. Czako and J. M. Bowman, J. Chem. Phys., 131, 244302 (2009)], a procedure called 1GB here. The calculations are then only ~ 10 times more time-consuming than with the standard QCT method, allowing thereby for considerable numerical saving. In this paper, we propose some theoretical arguments supporting the 1GB procedure and check its validity on model test cases as well as the prototype four-atom reaction OH+D_2 -> HOD+D

From angle-action to Cartesian coordinates: A key transformation for molecular dynamics

The transformation from angle-action variables to Cartesian coordinates is a crucial step of the (semi) classical description of bimolecular collisions and photo-fragmentations. The basic reason is that dynamical conditions corresponding to experiments are ideally generated in angle-action variables whereas the classical equations of motion are ideally solved in Cartesian coordinates by standard numerical approaches. To our knowledge, the previous transformation is available in the literature only for triatomic systems. The goal of the present work is to derive it for polyatomic ones.Comment: 10 pages, 11 figures, submitted to J. Chem. Phy

Quantum tunneling splittings from path-integral molecular dynamics.

We illustrate how path-integral molecular dynamics can be used to calculate ground-state tunnelling splittings in molecules or clusters. The method obtains the splittings from ratios of density matrix elements between the degenerate wells connected by the tunnelling. We propose a simple thermodynamic integration scheme for evaluating these elements. Numerical tests on fully dimensional malonaldehyde yield tunnelling splittings in good overall agreement with the results of diffusion Monte Carlo calculations.E.M., D.J.W., and S.C.A. acknowledge funding from the UK Engineering and Physical Sciences Research Council.This is the author accepted manuscript. The final version is available from the American Institute of Physics via http://dx.doi.org/10.1063/1.494386

Non-equilibrium dynamics from RPMD and CMD

We investigate the calculation of approximate non-equilibrium quantum time correlation functions (TCFs) using two popular path-integral-based molecular dynamics methods, ring-polymer molecular dynamics (RPMD) and centroid molecular dynamics (CMD). It is shown that for the cases of a sudden vertical excitation and an initial momentum impulse, both RPMD and CMD yield non-equilibrium TCFs for linear operators that are exact for high temperatures, in the t = 0 limit, and for harmonic potentials; the subset of these conditions that are preserved for non-equilibrium TCFs of non-linear operators is also discussed. Furthermore, it is shown that for these non-equilibrium initial conditions, both methods retain the connection to Matsubara dynamics that has previously been established for equilibrium initial conditions. Comparison of non-equilibrium TCFs from RPMD and CMD to Matsubara dynamics at short times reveals the orders in time to which the methods agree. Specifically, for the position-autocorrelation function associated with sudden vertical excitation, RPMD and CMD agree with Matsubara dynamics up to O(t^4) and O(t^1), respectively; for the position-autocorrelation function associated with an initial momentum impulse, RPMD and CMD agree with Matsubara dynamics up to O(t^5) and O(t^2), respectively. Numerical tests using model potentials for a wide range of non-equilibrium initial conditions show that RPMD and CMD yield non-equilibrium TCFs with an accuracy that is comparable to that for equilibrium TCFs. RPMD is also used to investigate excited-state proton transfer in a system-bath model, and it is compared to numerically exact calculations performed using a recently developed version of the Liouville space hierarchical equation of motion approach; again, similar accuracy is observed for non-equilibrium and equilibrium initial conditions

Boltzmann-conserving classical dynamics in quantum time-correlation functions: "Matsubara dynamics".

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or "classical Wigner approximation") results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e., a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads N → ∞, such that the lowest normal-mode frequencies take their "Matsubara" values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of ħ(2) at ħ(0) (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting "Matsubara" dynamics is inherently classical (since all terms O(ħ(2)) disappear from the Matsubara Liouvillian in the limit N → ∞) and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes and gives better agreement than LSC-IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods.T.J.H.H., M.J.W., and S.C.A. acknowledge funding from the U.K. Engineering and Physical Sciences Research Council. A.M. acknowledges the European Lifelong Learning Programme (LLP) for an Erasmus student placement scholarship. T.J.H.H. also acknowledges a Research Fellowship from Jesus College, Cambridge and helpful discussions with Dr. Adam Harper.This is the author accepted manuscript. The final version is available from AIP via http://dx.doi.org/10.1063/1.491631

Slow cross-symmetry phase relaxation in complex collisions

We discuss the effect of slow phase relaxation and the spin off-diagonal $S$-matrix correlations on the cross section energy oscillations and the time evolution of the highly excited intermediate systems formed in complex collisions. Such deformed intermediate complexes with strongly overlapping resonances can be formed in heavy ion collisions, bimolecular chemical reactions and atomic cluster collisions. The effects of quasiperiodic energy dependence of the cross sections, coherent rotation of the hyperdeformed $\simeq (3:1)$ intermediate complex, Schr\"odinger cat states and quantum-classical transition are studied for $^{24}$Mg+$^{28}$Si heavy ion scattering.Comment: 10 pages including 2 color ps figures. To be published in Physics of Atomic Nuclei (Yadernaya fizika