An instanton approach to intramolecular hydrogen exchange: Tunneling splittings in malonaldehyde and the hydrogenoxalate anion

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.

A Structure-Based Gaussian Expansion for Quantum Reaction Dynamics in Molecules: Application to Hydrogen Tunneling in Malonaldehyde

We developed an approximate method for quantum reaction dynamics simulations, namely, a structure-based Gaussian (SBG) expansion approach, where SBG bases for the expansion of the wave function Ψ, expressed by a product of single-atom Cartesian Gaussians centered at the positions of respective nuclei, are mainly placed around critical structures on reaction pathways such as on the intrinsic reaction coordinate (IRC) through a transition state. In the present approach, the "pseudo-lattice points" at which SBGs are deployed are selected in a perturbative manner so as to make moderate the expansion length. We first applied the SBG idea to a two-dimensional quadruple-well model and obtained accurate tunneling splitting values between the lowest four states. We then applied it to hydrogen tunneling in malonaldehyde and achieved a tunneling splitting of 27.1 cm^-1 with only 875 SBGs at the MP2/6-31G(d,p) level of theory, in good agreement with 25 cm^-1 by the more elaborate multiconfiguration time-dependent Hartree method. Reasonable results were also obtained for singly and doubly deuterated malonaldehyde. We analyzed the tunneling states by utilizing expansion coefficients of individual SBGs and found that 40-45% of the SBGs in Ψ are nonplanar structures and SBGs away from the IRC contribute a little to hydrogen transfer.

MC-QTAIM analysis reveals an exotic bond in coherently quantum superposed malonaldehyde

The proton between the two oxygen atoms of the malonaldehyde molecule experiences an effective double-well potential in which the proton's wavefunction is delocalized between the two wells. Herein we employ a state-of-the-art multi-component quantum theory of atoms in molecules partitioning scheme to obtain the molecular structure, i.e. atoms in molecules and bonding network, from the superposed ab initio wavefunctions of malonaldehyde. In contrast to the familiar clamped-proton portrayal of malonaldehyde, in which the proton forms a hydrogen basin, for the superposed states the hydrogen basin disappears and two novel hybrid oxygen-hydrogen basins appear instead, with an even distribution of the proton population between the two basins. The interaction between the hybrid basins is stabilizing thanks to an unprecedented mechanism. This involves the stabilizing classical Coulomb interaction of the one-proton density in one of the basins with one-electron density in the other basin. This stabilizing mechanism yields a bond foreign to the known bonding modes in chemistry.

Wavepacket dynamical study of H-atom tunneling in catecholate monoanion: the role of intermode couplings and energy flow

We present a study of H-atom tunneling in catecholate monoanion through wavepacket dynamical simulations. In our earlier study of this symmetrical double-well system [Phys. Chem. Chem. Phys., 2022, 24, 10887], a limited number of transition state modes were identified as being important for the tunneling process. These include the imaginary frequency mode Q₁, the CO scissor mode Q₁₀, and the OHO bending mode Q₂₉. In this work, starting from non-stationary initial states prepared with excitations in these modes, we have carried out wavepacket dynamics in two and three dimensional spaces. We analyse the dynamical effects of the intermode couplings, in particular the role of energy flow between the studied modes on H-atom tunneling. We find that while Q₁₀ strongly modulates the donor-acceptor distance, it does not exchange energy with Q₁. However, excitation in Q₂₉ or Q₁ does lead to rapid energy exchange between these modes, which modifies the tunneling rate at early times.

Multidimensional H-atom tunneling in the catecholate monoanion

We present the catecholate monoanion as a new model system for the study of multidimensional tunneling. It has a symmetrical O-H double-well structure, and the H atom motion between the two wells is coupled to both low and high frequency modes with different strengths. With a view to studying mode-specific tunneling in the catecholate monoanion, we have developed a full (33) dimensional potential energy surface in transition state (TS) normal modes using a Distributed Gaussian Empirical Valence Bond (DGEVB) based approach. We have computed eigenstates in different subspaces using both unrelaxed and relaxed potentials based on the DGEVB model. With unrelaxed potentials, we present results up to 7D subspaces that include the imaginary frequency mode and six modes coupled to it. With relaxed potentials, we focus on the two most important coupling modes. The structures of the ground and vibrationally excited eigenstates are discussed for both approaches and mode-specific tunneling splitting and their trends are presented.

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 Mechanical Investigation of Mode-Specific Tunneling upon Fundamental Excitation in Malonaldehyde

We present a quantum mechanical study of mode-specific tunneling upon fundamental excitation in malonaldehyde with a multidimensional theory that utilizes the saddle-point normal coordinates. We find that a ring-deformation normal mode is as essential as the well-known imaginary-frequency normal mode in the multidimensional investigation. The changes in tunneling splittings upon fundamental excitation are calculated. The results are competitive with those from a recently developed mixed classical-quantum method. Moreover, the results are qualitatively consistent with experiment for about half of all the modes.

Multidimensional Hamiltonian for tunneling with position-dependent mass

A multidimensional Hamiltonian for tunneling is formulated, based on the mode with imaginary frequency of the transition state as a reaction coordinate. To prepare it for diagonalization, it is transformed into a lower-dimension Hamiltonian by incorporating modes that move faster than the tunneling into a coordinate-dependent kinetic energy operator, for which a Hermitian form is chosen and tested for stability of the eigenvalues. After transformation to a three-dimensional form, which includes two normal modes strongly coupled to the tunneling mode, this Hamiltonian is diagonalized in terms of a basis set of harmonic oscillator functions centered at the transition state. This involves a sparse matrix which is easily partially diagonalized to yield tunneling splittings for the zero-point level and the two fundamental levels of the coupled modes. The method is tested on the well-known benchmark molecule malonaldehyde and a deuterium isotopomer, for which these splittings have been measured. Satisfactory agreement with experiment results is obtained.

The ground state tunneling splitting and the zero point energy of malonaldehyde: A quantum Monte Carlo determination

Quantum dynamics calculations of the ground state tunneling splitting and of the zero point energy of malonaldehyde on the full dimensional potential energy surface proposed by Yagi et al. [J. Chem. Phys. 1154, 10647 (2001)] are reported. The exact diffusion Monte Carlo and the projection operator imaginary time spectral evolution methods are used to compute accurate benchmark results for this 21-dimensional ab initio potential energy surface. A tunneling splitting of 25.7+/-0.3 cm-1 is obtained, and the vibrational ground state energy is found to be 15 122+/-4 cm-1. Isotopic substitution of the tunneling hydrogen modifies the tunneling splitting down to 3.21+/-0.09 cm-1 and the vibrational ground state energy to 14 385+/-2 cm-1. The computed tunneling splittings are slightly higher than the experimental values as expected from the potential energy surface which slightly underestimates the barrier height, and they are slightly lower than the results from the instanton theory obtained using the same potential energy surface.

Quantum Tunneling in the Midrange Vibrational Fundamentals of Tropolone

The Fourier transform infrared spectrum of tropolone(OH) vapor in the 1175-1700 cm(-1) region is reported at 0.0025 and 0.10 cm(-1) spectral resolutions. The 12 vibrational fundamentals in this region of rapidly rising vibrational state density are dominated by mixtures of the CC, CO, CCH, and COH internal coordinates. Estimates based on the measurement of sharp Q branch peaks are reported for 11 of the spectral doublet component separations DS(v) = |Delta(v) +/- Delta(0)|. Delta(0) = 0.974 cm(-1) is the known zero-point splitting, and three a(1) modes show tunneling splittings Delta(v) approximately Delta(0), four b(2) modes show splittings Delta(v) approximately 0.90Delta(0), and the remaining four modes show splittings Delta(v) falling 5-14% from Delta(0.) Significantly, the splitting for the nominal COH bending mode nu(8) (a(1)) is small, that is, 10% from Delta(0). Many of the vibrational excited states demonstrate strong anharmonic behavior, but there are only mild perturbations on the tautomerization mechanism driving Delta(0). The data suggest, especially for the higher frequency a(1) fundamentals, the onset of selective intramolecular vibrational energy redistribution processes that are fast on the time scale of the tautomerization process. These appear to delocalize and smooth out the topographical modifications of the zero-point potential energy surface that are anticipated to follow absorption of the nu(v) photon. Further, the spectra show the propensity for the Delta(v) splittings of b(2) and other complex vibrations to be damped relative to Delta(0).

Simple and accurate method to evaluate tunneling splitting in polyatomic molecules

A practical and accurate semiclassical method for calculating the tunneling splitting of the ground state in polyatomic molecules is presented based on a recent version of the instanton theory [J. Chem. Phys. 115, 6881 (2001)]. The method uses ab initio quantum chemical data for the potential energy surface without any concomitant extrapolation and requires only a small number of ab initio data points to get convergence even for large molecules. This enables one to use an advanced level of electronic structure theory and achieve a high accuracy of the result. The method is applied to the 9-atomic malonaldehyde molecule by making use of the potential energy surface at the level of CCSD(T) with the hybrid basis set of aug-cc-pVTZ (for oxygen atoms and the transferred hydrogen atom) and cc-pVTZ (for other atoms).

Ground-state and vibrationally assisted tunneling in the formic acid dimer

The previously developed instanton theory [G. V. Mil'nikov and H. Nakamura, J. Chem. Phys. 122, 124311 (2005)] is applied to the calculation of vibrationally assisted tunneling splitting of the deuterated formic acid dimer (DCOOH)2 with all the degrees of freedom taken into account. The ground-state tunnel splitting is determined by the density-functional theory combined with coupled cluster level of quantum chemistry to be 0.0038 cm(-1) which is comparable to the experimental value of 0.0029 cm(-1). Further, the tunnel splittings of fundamental excitations are estimated for frequencies below 300 cm(-1). In this energy range it is found that the excitation modes may either enhance or suppress tunneling as compared to the ground state. For the higher-frequency modes a rapid growth of the tunnel splitting is observed. At frequencies above 1000 cm(-1) the semiclassical solution becomes unstable and no reliable tunneling splittings can be obtained. This is in vast contrast to the adiabatic approximation to the instanton theory in which the tunnel splittings can be retrieved up to 3000 cm(-1). We discuss this disparity from the viewpoint of the multidimensional character of tunneling in hydrogen bonds and the adiabatic approximation is concluded to be inaccurate.

Tunneling dynamics of double proton transfer in formic acid and benzoic acid dimers

Direct dynamics calculations based on instanton techniques are reported of tunneling splittings due to double proton transfer in formic and benzoic acid dimers. The results are used to assign the observed splittings to levels for which the authors of the high-resolution spectra could not provide a definitive assignment. In both cases the splitting is shown to be due mainly to the zero-point level rather than to the vibrationally or electronically excited level whose spectrum was investigated. This leads to zero-point splittings of 375 MHz for (DCOOH)(2) and 1107 MHz for the benzoic acid dimer. Thus, contrary to earlier calculations, it is found that the splitting is considerably larger in the benzoic than in the formic acid dimer. The calculations are extended to solid benzoic acid where the asymmetry of the proton-transfer potential induced by the crystal can be overcome by suitable doping. This has allowed direct measurement of the interactions responsible for double proton transfer, which were found to be much larger than those in the isolated dimer. To account for this observation both static and dynamic effects of the crystal forces on the intradimer hydrogen bonds are included in the calculations. The same methodology, extended to higher temperatures, is used to calculate rate constants for HH, HD, and DD transfers in neat benzoic acid crystals. The results are in good agreement with reported experimental rate constants measured by NMR relaxometry and, if allowance is made for small structural changes induced by doping, with the transfer matrix elements observed in doped crystals. Hence the method used allows a unified description of tunneling splittings in the gas phase and in doped crystals as well as of transfer rates in neat crystals.

Instanton theory for the tunneling splitting of low vibrationally excited states

We develop the instanton theory for calculating the tunneling splitting of excited states. For the case of low vibrational quantum states we derive a canonically invariant formula which is applicable to a multidimensional system of arbitrary Riemannian metric. The effect of multidimensionality in relation to the vibrational excitation is explained in terms of the effective frequencies along the instanton trajectory. The theory is demonstrated to work well by taking HO2 molecule as an example.

The ground state tunneling splitting of malonaldehyde: Accurate full dimensional quantum dynamics calculations

Benchmark calculations of the tunneling splitting in malonaldehyde using the full dimensional potential proposed by Yagi et al. are reported. Two exact quantum dynamics methods are used: the multiconfigurational time-dependent Hartree (MCTDH) approach and the diffusion Monte Carlo based projection operator imaginary time spectral evolution (POITSE) method. A ground state tunneling splitting of 25.7+/-0.3 cm(-1) is calculated using POITSE. The MCTDH computation yields 25 cm(-1) converged to about 10% accuracy. These rigorous results are used to evaluate the accuracy of approximate dynamical approaches, e.g., the instanton theory.

The ground-state tunneling splitting of various carboxylic acid dimers

Carboxylic acid dimers in gas phase reveal ground-state tunneling splittings due to a double proton transfer between the two subunits. In this study we apply a recently developed accurate semiclassical method to determine the ground-state tunneling splittings of eight different carboxylic acid derivative dimers (formic acid, benzoic acid, carbamic acid, fluoro formic acid, carbonic acid, glyoxylic acid, acrylic acid, and N,N-dimethyl carbamic acid) and their fully deuterated analogs. The calculated splittings range from 5.3e-4 to 0.13 cm(-1) (for the deuterated species from 2.8e-7 to 3.3e-4 cm(-1)), thus indicating a strong substituent dependence of the splitting, which varies by more than two orders of magnitude. One reason for differences in the splittings could be addressed to different barriers heights, which vary from 6.3 to 8.8 kcal/mol, due to different mesomeric stabilization of the various transition states. The calculated splittings were compared to available experimental data and good agreement was found. A correlation could be found between the tunneling splitting and the energy barrier of the double proton transfer, as the splitting increases with increased strength of the hydrogen bonds. From this correlation an empirical formula was derived, which allows the prediction of the ground-state tunneling splitting of carboxylic acid dimers at a very low cost and the tunneling splittings for parahalogen substituted benzoic acid dimers is predicted.