DOIT: a program to calculate thermal rate constants and mode‐specific tunneling splittings directly from quantum‐chemical calculations

In this contribution we discuss computational aspects of a recently introduced method for the calculation of proton tunneling rate constants, and tunneling splittings, which has been applied to molecules and complexes, and should apply equally well to bulk materials. The method is based on instanton theory, adapted so as to permit a direct link to the output of quantum‐chemical codes. It is implemented in the DOIT (dynamics of instanton tunneling) code, which calculates temperature‐dependent tunneling rate constants and mode‐specific tunneling splittings. As input, it uses the structure, energy, and vibrational force field of the stationary configurations along the reaction coordinate, computed by conventional quantum‐chemical programs. The method avoids the difficult problem of calculating the exact least‐action trajectory, known as the instanton path, and instead focusses on the corresponding instanton action, because it governs the dynamic properties. To approximate this action for a multidimensional system, the program starts from the one‐dimensional instanton action along the reaction coordinate, which can be obtained without difficulty. It then applies correction terms for the coupling to the other vibrational degrees of freedom, which are treated as harmonic oscillators (transverse normal modes). The couplings are assumed linear in these modes. Depending on the frequency and the character of the transverse modes, they may either decrease or increase the action, i.e., help or hinder the transfer. A number of tests have shown that the program is at least as accurate as alternative programs based on transition‐state theory with tunneling corrections, and is also much less demanding in computer time, thus allowing application to much larger systems. An outline of the instanton formalism is presented, some new developments are introduced, and special attention is paid to the connection with quantum‐chemical codes. Possible sources of error are investigated. To show the program in action, calculations are presented of tunneling rates and splittings associated with triple proton transfer in the chiral water trimer. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 787–801, 2001

[1]  B. C. Garrett,et al.  Current status of transition-state theory , 1983 .

[2]  Gert-Ludwig Ingold,et al.  Quantum Brownian motion: The functional integral approach , 1988 .

[3]  T. Sewell,et al.  SEMICLASSICAL CALCULATIONS OF TUNNELING SPLITTING IN MALONALDEHYDE , 1995 .

[4]  M. Zgierski,et al.  Photochromism of salicylideneaniline (SA). How the photochromic transient is created: A theoretical approach , 2000 .

[5]  W. Miller,et al.  A semiclassical tunneling model for use in classical trajectory simulations , 1989 .

[6]  B. C. Garrett,et al.  A general small-curvature approximation for transition-state-theory transmission coefficients , 1981 .

[7]  John E. Adams,et al.  Reaction path Hamiltonian for polyatomic molecules , 1980 .

[8]  P. Bunker,et al.  Molecular symmetry and spectroscopy , 1979 .

[9]  B. C. Garrett,et al.  A least‐action variational method for calculating multidimensional tunneling probabilities for chemical reactions , 1983 .

[10]  R. Saykally,et al.  Measurement of quantum tunneling between chiral isomers of the cyclic water trimer. , 1992, Science.

[11]  A. Fernández-Ramos,et al.  A comparison of two methods for direct tunneling dynamics: Hydrogen exchange in the glycolate anion as a test case , 1997 .

[12]  D. Clary,et al.  Calculations of the tunneling splittings in water dimer and trimer using diffusion Monte Carlo , 1995 .

[13]  W. Siebrand,et al.  Mode‐specific hydrogen tunneling in tropolone: An instanton approach , 1996 .

[14]  Curtis G. Callan,et al.  Fate of the false vacuum. II. First quantum corrections , 1977 .

[15]  D. Thompson,et al.  Initial conditions and paths in semiclassical tunneling , 1997 .

[16]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[17]  A. Fernández-Ramos,et al.  Double proton transfer in the complex of acetic acid with methanol: Theory versus experiment , 2001 .

[18]  A. Fernández-Ramos,et al.  A direct-dynamics study of the zwitterion-to-neutral interconversion of glycine in aqueous solution , 2000 .

[19]  W. Tung.  A comment on Melosh transformation and the photoproduction of resonances , 1977 .

[20]  D. Thompson,et al.  A multidimensional semiclassical method for treating tunneling in molecular collisions , 1996 .

[21]  B. C. Garrett,et al.  Application of the large-curvature tunneling approximation to polyatomic molecules: Abstraction of H or D by methyl radical , 1989 .

[22]  D. Makarov,et al.  Quantum chemical dynamics in two dimensions , 1993 .

[23]  J. Leszczynski,et al.  A direct-dynamics study of proton transfer through water bridges in guanine and 7-azaindole , 2000 .

[24]  W. Miller.  Semiclassical limit of quantum mechanical transition state theory for nonseparable systems , 1975 .

[25]  Hänggi,et al.  Numerical study of tunneling in a dissipative system. , 1990, Physical Review B (Condensed Matter).

[26]  J. Langer.  Statistical theory of the decay of metastable states , 1969 .

[27]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[28]  A. Fernández-Ramos,et al.  Mode-specific tunneling splittings in 9-hydroxyphenalenone: Comparison of two methods for direct tunneling dynamics , 1998 .

[29]  D. C. Clary,et al.  The Theory of Chemical Reaction Dynamics , 1986 .

[30]  A. Fernández-Ramos,et al.  Direct-dynamics approach to catalytic effects: The tautomerization of 3-hydroxyisoquinoline as a test case , 2000 .

[31]  F. Zerbetto,et al.  Dynamics of molecular inversion: An instanton approach , 1995 .

[32]  F. Zerbetto,et al.  An Instanton approach to hindered torsions: methyl glycolate — a case study , 1997 .

[33]  W. Miller,et al.  Model studies of mode specificity in unimolecular reaction dynamics , 1980 .

[34]  F. Zerbetto,et al.  Tunneling splittings from ab initio data: indoline, a test case , 1995 .

[35]  A. Fernández-Ramos,et al.  Proton tunnelling in polyatomic molecules: A direct-dynamics instanton approach , 1999 .

[36]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[37]  Carl Eckart,et al.  Some Studies Concerning Rotating Axes and Polyatomic Molecules , 1935 .

[38]  B. Rode,et al.  Predictions of rate constants and estimates for tunneling splittings of concerted proton transfer in small cyclic water clusters , 1998 .

[39]  P. Kozlowski,et al.  Dynamics of tautomerism in porphine: An instanton approach , 1998 .

[40]  W. Siebrand,et al.  An instanton approach to intramolecular hydrogen exchange: Tunneling splittings in malonaldehyde and the hydrogenoxalate anion , 1995 .

[41]  D. Makarov,et al.  Low-temperature chemical reactions. Effect of symmetrically coupled vibrations in collinear exchange reactions , 1991 .

[42]  D. Makarov,et al.  Quantum dynamics in low-temperature chemistry , 1993 .