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.     

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcalmol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively.     

The all-Cartesian reaction plane Hamiltonian: formulation and application to the H-atom transfer in tropolone

In this work we present an all-Cartesian reaction surface approach, where the large amplitude coordinates span the so-called reaction plane, that is, the unique plane defined by the two minima and the saddle-point structure of an isomerization reaction. Orthogonal modes are treated within harmonic approximation which gives the total Hamiltonian an almost separable form that is suitable for multidimensional quantum dynamics calculations. The reaction plane Hamiltonian is constructed for the H-atom transfer in tropolone as an example for a system with an intramolecular O...H-O hydrogen bond. We find ground-state tunneling splittings of 3.5 and 0.16 cm(-1) for the normal and deuterated species, respectively. We calculated infrared-absorption spectra for a four-dimensional model focusing on the low-frequency region. Here, we identify a reaction mode which is closely connected to the tautomerization that is reflected in the increase of tunneling splitting to 18 cm(-1) upon excitation.     

The Barrier to Proton Transfer in the Dimer of Formic Acid: A Pure Rotational Study

The rotational spectra of three C‐deuterated isotopologues of the dimer of formic acid have been measured, thanks to the small dipole moment induced by asymmetric H→D substitution(s). For the DCOOH–HCOOH species, the concerted double proton transfer of the two hydroxy hydrogen atoms takes place between two equivalent minima and generates a tunneling splitting of 331.2(6) MHz. This splitting can be reproduced by a 3D model with a barrier of 2559 cm^?1 (30.6 kJ mol^?1) as obtained from theoretical calculations.     

Symmetric double proton tunneling in formic acid dimer: a diabatic basis approach

A model of double proton tunneling in formic acid dimer is developed using a reaction surface Hamiltonian. The surface includes the symmetric OH stretch plus the in-plane stretch and bend interdimer vibrations. The surface Hamiltonian is coupled to a bath of five A1g and B3g normal modes obtained at the D2h transition state structure. Eigenstates are calculated using Davidson and block-Davidson iterative methods. Strong mode specific effects are found in the tunneling splittings for the reaction surface, where splittings are enhanced upon excitation of the interdimer bend motion. The results are interpreted within the framework of a diabatic representation of reaction surface modes. The splitting patterns observed for the reaction surface eigenstates are only slightly modified upon coupling to the bath states. Splitting patterns for the bath states are also determined. It is found that predicting these splittings is greatly complicated by subtle mixings with the inter-dimer bend states.     

Ground and excited states of the monomer and dimer of certain carboxylic acids

The ground-state properties of the monomer and the dimer of formic acid, acetic acid, and benzoic acid have been investigated using Hartree-Fock (HF) and density functional theory (DFT) methods using the 6-311++G(d,p) basis set. Some of the low-lying excited states have been studied using the time-dependent density functional theory (TDDFT) with LDA and B3LYP functionals and also employing complete-active-space-self-consistent-field (CASSCF) and multireference configuration interaction (MRCI) methodologies. DFT calculations predict the ground-state geometries in quantitative agreement with the available experimental results. The computed binding energies for the three carboxylic acid dimers are also in accord with the known thermodynamic data. The TDDFT predicted wavelengths corresponding to the lowest energy n-pi* transition in formic acid (214 nm) and acetic acid (214 nm) and the pi-pi* transition in benzoic acid (255 nm) are comparable to the experimentally observed absorption maxima. In addition, TDDFT calculations predict qualitatively correctly the blue shift (4-5 nm) in the excitation energy for the pi-pi* transition in going from the monomer to the dimer of formic acid and acetic acid and the red shift (approximately 19 nm) in pi-pi* transition in going from benzoic acid monomer to dimer. This also indicates that the electronic interaction arising from the hydrogen bonds between the monomers is marginal in all three carboxylic acids investigated.     

NMR studies of solid state—solvent and H/D isotope effects on hydrogen bond geometries of 1:1 complexes of collidine with carboxylic acids

¹H and ¹⁵N NMR spectra of 10 complexes exhibiting strong OHN hydrogen bonds formed by ¹⁵N-labeled collidine and different proton donors, partially deuterated in mobile proton sites, have been observed by low-temperature NMR spectroscopy using a low-freezing CDF₃/CDF₂Cl mixture as polar aprotic solvent. The following proton donors have been used: HCl, formic acid, acetic acid, various substituted benzoic acids and HBF₄. The slow hydrogen bond exchange regime could be reached below 140 K, which allowed us to resolve ¹⁵N signal splittings due to H/D isotopic substitution. The valence bond order model is used to link the observed NMR parameters to hydrogen bond geometries. The results are compared to those obtained previously [Magn. Reson. Chem. 39 (2001) S18] for the same complexes in the organic solids. The increase of the dielectric constant from the organic solids to the solution (30 at 130 K) leads to a change of the hydrogen bond geometries along the geometric correlation line towards the zwitterionic structures, where the proton is partially transferred from oxygen to nitrogen. Whereas the changes of spectroscopic and, hence, geometric parameters are small for the systems which are already zwitterionic in the solid state, large changes are observed for molecular complexes which exhibit almost a full proton transfer from oxygen to nitrogen in the polar liquid solvent.     

Vibrational quenching of excitonic splittings in H-bonded molecular dimers: the electronic Davydov splittings cannot match experiment

The S(1)/S(2) state exciton splittings of symmetric doubly hydrogen-bonded gas-phase dimers provide spectroscopic benchmarks for the excited-state electronic couplings between UV chromophores. These have important implications for electronic energy transfer in multichromophoric systems ranging from photosynthetic light-harvesting antennae to photosynthetic reaction centers, conjugated polymers, molecular crystals, and nucleic acids. We provide laser spectroscopic data on the S(1)/S(2) excitonic splitting Δ(exp) of the doubly H-bonded o-cyanophenol (oCP) dimer and compare to the splittings of the dimers of (2-aminopyridine)(2), [(2AP)(2)], (2-pyridone)(2), [(2PY)(2)], (benzoic acid)(2), [(BZA)(2)], and (benzonitrile)(2), [(BN)(2)]. The experimental S(1)/S(2) excitonic splittings are Δ(exp) = 16.4 cm(-1) for (oCP)(2), 11.5 cm(-1) for (2AP)(2), 43.5 cm(-1) for (2PY)(2), and <1 cm(-1) for (BZA)(2). In contrast, the vertical S(1)/S(2) energy gaps Δ(calc) calculated by the approximate second-order coupled cluster (CC2) method for the same dimers are 10-40 times larger than the Δ(exp) values. The qualitative failure of this and other ab initio methods to reproduce the exciton splitting Δ(exp) arises from the Born-Oppenheimer (BO) approximation, which implicitly assumes the strong-coupling case and cannot be employed to evaluate excitonic splittings of systems that are in the weak-coupling limit. Given typical H-bond distances and oscillator strengths, the majority of H-bonded dimers lie in the weak-coupling limit. In this case, the monomer electronic-vibrational coupling upon electronic excitation must be accounted for; the excitonic splittings arise between the vibronic (and not the electronic) transitions. The discrepancy between the BO-based splittings Δ(calc) and the much smaller experimental Δ(exp) values is resolved by taking into account the quenching of the BO splitting by the intramolecular vibronic coupling in the monomer S(1) ← S(0) excitation. The vibrational quenching factors Γ for the five dimers (oCP)(2), (2AP)(2), (2AP)(2), (BN)(2), and (BZA)(2) lie in the range Γ = 0.03-0.2. The quenched excitonic splittings Γ[middle dot]Δ(calc) are found to be in very good agreement with the observed splittings Δ(exp). The vibrational quenching approach predicts reliable Δ(exp) values for the investigated dimers, confirms the importance of vibrational quenching of the electronic Davydov splittings, and provides a sound basis for predicting realistic exciton splittings in multichromophoric systems.     

Vibrational coupling in carboxylic acid dimers

The vibrational level splitting in the ground electronic state of carboxylic acid dimers mediated by the doubly hydrogen-bonded networks are investigated using pure and mixed dimers of benzoic acid with formic acid as molecular prototypes. Within the 0-2000-cm(-1) range, the frequencies for the fundamental and combination vibrations of the two dimers are experimentally measured by using dispersed fluorescence spectroscopy in a supersonic jet expansion. Density-functional-theory calculations predict that most of the dimer vibrations are essentially in-phase and out-of-phase combinations of the monomer modes, and many of such combinations show significantly large splitting in vibrational frequencies. The infrared spectrum of the jet-cooled benzoic acid dimer, reported recently by Bakker et al. [J. Chem. Phys. 119, 11180 (2003)], has been used along with the dispersed fluorescence spectra to analyze the coupled g-u vibrational levels. Assignments of the dispersed fluorescence spectra of the mixed dimer are suggested by comparing the vibronic features with those in the homodimer spectrum and the predictions of density-functional-theory calculation. The fluorescence spectra measured by excitations of the low-lying single vibronic levels of the mixed dimer reveal that the hydrogen-bond vibrations are extensively mixed with the ring modes in the S1 surface.     

Infrared spectroscopy of the intramolecular hydrogen bond in acethylacetone: a computational approach

The intramolecular hydrogen bond in the enol-acethylacetone (ACAC) is investigated by performing reduced-dimensional quantum calculations. To analyze the shared proton vibrations, two sets of coordinates were employed: normal mode coordinates describing the motion in the vicinity of the most stable configuration, and internal coordinates accounting for the double minimum proton motion. It is proved that the extreme broadness of the OH-stretch band in ACAC is a consequence of the coexistence of two enol-ACAC structures: the global minimum and the transition state for rotation of the distal methyl group. Further, a ground-state tunneling splitting of 116 cm(-1) is found, and it is shown that the inclusion of the kinematic coupling is mandatory when treating large-amplitude proton motion. In the OH-stretch direction a splitting of 853 cm(-1) was predicted.     

An infrared linear dichroism study of carboxylic acids oriented in stretched polyethylene

Polarized IR spectra of four carboxylic acids have been recorded using the stretched polymer method. Benzoic, propanoic and propynoic acids form cyclic dimers in polyethylene, while formic acid forms an α type polymer. The transition moment directions of 23 vibrations of benzoic acid dimer are given and the carboxylic ring in-plane vibrations of the studied acids are discussed in the light of the obtained results.     

Excited-state double proton transfer dynamics of model DNA base pairs: 7-hydroxyquinoline dimers

The excited-state double proton transfer of model DNA base pairs, 7-hydroxyquinoline dimers, in benzene has been investigated using picosecond time-resolved fluorescence spectroscopy. Upon excitation, whereas singly hydrogen-bonded noncyclic dimers do not go through tautomerization within the relaxation time of 1400 ps, doubly hydrogen-bonded cyclic dimers undergo excited-state double proton transfer on the time scale of 25 ps to form tautomeric dimers, which subsequently undergo a conformational change in 180 ps to produce singly hydrogen-bonded tautomers. The rate constant of the double proton transfer reaction is temperature-independent, showing a large kinetic isotope effect of 5.2, suggesting that the rate is governed mostly by tunneling.     

Isotope effects associated with tunneling and double proton transfer in the hydrogen bonds of benzoic acid

The isotope effects associated with double proton transfer in the hydrogen bonds of benzoic acid (BA) dimers have been measured using field-cycling (1)H NMR relaxometry and quasielastic neutron scattering. By studying mixed isotope (hydrogen and deuterium) samples, the dynamics of three isotopologues, BA-HH, BA-HD, and BA-DD, have been investigated. Low temperature measurements provide accurate measurements of the incoherent tunneling rate, k(0). This parameter scales accurately with the mass number, m, according to the formula k(0)=(E/m)e(-Fm) providing conclusive evidence that the proton transfer process is a strongly correlated motion of two hydrons. Furthermore, we conclude that the tunneling pathway is the same for the three isotopologue species. Measurements at higher temperatures illuminate the through barrier processes that are mediated via intermediate or excited vibrational states. In parallel with the investigation of proton transfer dynamics, the theoretical and experimental aspects of studying spin-lattice relaxation in single crystals of mixed isotope samples are investigated in depth. Heteronuclear dipolar interactions between (1)H and (2)H isotopes contribute significantly to the overall proton spin-lattice relaxation and it is shown that these must be modeled correctly to obtain accurate values for the proton transfer rates. Since the sample used in the NMR measurements was a single crystal, full account of the orientation dependence of the spin-lattice relaxation with respect to the applied B field was incorporated into the data analysis.     

Splitting of carbonyl stretching frequencies from transition dipole coupling

The observed splitting of the carbonyl stretching frequencies of hydrogen bonded carboxylic acid dimers has been explained by other authors in terms of transition dipole coupling. In this paper it is shown that the transition dipole interaction for carbonyl groups gives very small splittings when reasonable dipole moment changes are used.     

Concerted hydrogen exchange tunneling in formic acid dimer

The effect of conformational relaxation on the quantum dynamics of the hydrogen exchange tunneling is studied in the D2h subspace of formic acid dimer. The fully coupled quantum dynamics in up to six dimensions are derived for potential energy hypersurfaces interpolated directly from hybrid density functional calculations with and without geometry relaxation. For a calculated electronic barrier height of 35.0 kJ/mol the vibrational ground state shows a tunneling splitting of 0.0013 cm(-1). The results support the vibrational assignment of Madeja and Havenith [J. Chem. Phys. 2002, 117, 7162-7168]. Fully coupled ro-vibrational calculations demonstrate the compatibility of experimentally observed inertia defects with in-plane hydrogen exchange tunneling dynamics in formic acid dimer.     

Excited-state double proton transfer of 7-azaindole dimers in a low-temperature organic glass

The excited-state double proton transfer of model DNA base pairs, 7-azaindole (7AI) dimers, is explored in a low-temperature organic glass of n-dodecane using picosecond time-resolved fluorescence spectroscopy. Reaction mechanisms are found to depend on the conformations of 7AI dimers at the moment of excitation; whereas planar conformers tautomerize rapidly (<10 ps), twisted conformers undergo double proton transfer to form tautomeric dimers on the time scale of 250 ps at 8 K. The proton transfer is found to consist of two orthogonal steps: precursor-configurational optimization and intrinsic proton transfer via tunneling. The rate is almost isotope independent at cryogenic temperatures because configurational optimization is the rate-determining step of the overall proton transfer. This optimization is assisted by lattice vibrations below 150 K or by librational motions above 150 K.     

An ab initio study of tunneling splittings in the water dimer

Tunneling splittings in the water dimer have been determined by the semiclassical WKB method, based on pathways characterized at the computational level of second-order M?ller-Plesset (MP2) theory with basis sets of aug-cc-pVTZ quality. This calculation takes into account all three acceptor tunneling, donor-acceptor interchange, and bifurcation tunneling rearrangements of the water dimer. The tunneling splittings were evaluated as 7.73 cm(-1) (large splitting) and 0.42 cm(-1) (small splitting), which are in good agreement with the corresponding experimental values of 11.18 cm(-1) and 0.70 cm(-1), respectively.     

Generalized approximation to the reaction path: the formic acid dimer case

A set of mass-weighted internal coordinates was derived and applied to the double proton transfer reaction in the formic acid dimer (FAD). The coordinate set was obtained starting from the Hirschfelder "mobile" by an optimization procedure consisting of a sequence of kinematic rotations. In FAD, the optimization procedure leads to three coordinates that do change significantly along the reaction path. These coordinates span the reaction space, whereas the remaining modes are treated in a harmonic approximation. The effect that the dimer dissociative motion has on the ground and excited vibrational states dynamics was explored. In the frequency region corresponding to the symmetric OH-stretch vibration four doublets have been identified with splittings of 2.76, 0.07, 0.60, and 4.03 cm(-1).     

Hydrogen bonding in substituted formic acid dimers

The hydrogen-bonded dimers of formic acid derivatives XCOOH (X = H, F, Cl, and CH3) have been investigated using density functional theory (B3LYP) and second-order M?ller-Plesset perturbation (MP2) methods, with the geometry optimization carried out using 6-311++G(2d,2p) basis set. The dimerization energies calculated using aug-cc-pVXZ (with X = D and T) basis have been extrapolated to infinite basis set limit using the standard methodology. The results indicate that the fluorine-substituted formic acid dimer is the most stable one in comparison to the others. Topological analysis carried out using Bader's atoms in molecules (AIM) theory shows good correlation of the values of electron density and its Laplacian at the bond critical points (BCP) with the hydrogen bond length in the dimers. Natural bond orbital (NBO) analysis carried out to study the charge transfer from the proton acceptor to the antibonding orbital of the X-H bond in the complexes reveals that most of the dimers are associated with conventional H-bonding except a few, where improper blue-shifting hydrogen bonds are found to be present.