Transitions induced in a double minimum system by interaction with a quantum mechanical heat bath

1H/17O Chemical Shift Waves in Carboxyl-Bridged Hydrogen Bond Networks in Organic Solids

We report solid-state 1H and 17O NMR results for four 17O-labeled organic compounds each containing an extensive carboxyl-bridged hydrogen bond (CBHB) network in the crystal lattice: tetrabutylammonium hydrogen di-[17O2]salicylate (1), [17O4]quinolinic acid (2), [17O4]dinicotinic acid (3), and [17O2]Gly/[17O2]Gly·HCl cocrystal (4). The 1H isotropic chemical shifts found for protons involved in different CBHB networks are between 8.2 and 20.5 ppm, which reflect very different hydrogen-bonding environments. Similarly, the 17O isotropic chemical shifts found for the carboxylate oxygen atoms in CBHB networks, spanning a large range between 166 and 341 ppm, are also remarkably sensitive to the hydrogen-bonding environments. We introduced a simple graphical representation in which 1H and 17O chemical shifts are displayed along the H and O atomic chains that form the CBHB network. In such a depiction, because wavy patterns are often observed, we refer to these wavy patterns as 1H/17O chemical shift waves. Typical patterns of 1H/17O chemical shift waves in CBHB networks are discussed. The reported 1H and 17O NMR parameters for the CBHB network models examined in this study can serve as benchmarks to aid in spectral interpretation for CBHB networks in proteins.

Quantum and classical effects in DNA point mutations: Watson-Crick tautomerism in AT and GC base pairs

Proton transfer along the hydrogen bonds of DNA can lead to the creation of short-lived, but biologically relevant point mutations that can further lead to gene mutation and, potentially, cancer. In this work, the energy landscape of the canonical A-T and G-C base pairs (standard, amino-keto) to tautomeric A*-T* and G*-C* (non-standard, imino-enol) Watson-Crick DNA base pairs is modelled with density functional theory and machine-learning nudge-elastic band methods. We calculate the energy barriers and tunnelling rates of hydrogen transfer between and within each base monomer (A, T, G and C). We show that the role of tunnelling in A-T tautomerisation is statistically unlikely due to the presence of a small reverse reaction barrier. On the contrary, the thermal populations of the G*-C* point mutation could be non-trivial and propagate through the replisome. For the direct intramolecular transfer, the reaction is hindered by a substantial energy barrier. However, our calculations indicate that tautomeric bases in their monomeric form have remarkably long lifetimes.

A quantum mechanical alternative to the Arrhenius equation in the interpretation of proton spin-lattice relaxation data for the methyl groups in solids

The theory of nuclear spin-lattice relaxation in methyl groups in solids has been a recurring problem in nuclear magnetic resonance (NMR) spectroscopy. The current view is that, except for extreme cases of low torsional barriers where special quantum effects are at stake, the relaxation behaviour of the nuclear spins in methyl groups is controlled by thermally activated classical jumps of the methyl group between its three orientations. The temperature effects on the relaxation rates can be modelled by Arrhenius behaviour of the correlation time of the jump process. The entire variety of relaxation effects in protonated methyl groups have recently been given a consistent quantum mechanical explanation not invoking the jump model regardless of the temperature range. It exploits the damped quantum rotation (DQR) theory originally developed to describe NMR line shape effects for hindered methyl groups. In the DQR model, the incoherent dynamics of the methyl group include two quantum rate (i.e., coherence-damping) processes. For proton relaxation only one of these processes is relevant. In this paper, temperature-dependent proton spin-lattice relaxation data for the methyl groups in polycrystalline methyltriphenyl silane and methyltriphenyl germanium, both deuterated in aromatic positions, are reported and interpreted in terms of the DQR model. A comparison with the conventional approach exploiting the phenomenological Arrhenius equation is made. The present observations provide further indications that incoherent motions of molecular moieties in the condensed phase can retain quantum character over much broader temperature range than is commonly thought.

Modelling proton tunnelling in the adenine–thymine base pair

Environment-assisted quantum tunnelling is unlikely to play a significant role in adenine–thymine tautomerisation.

Molecular recognition and interaction between uracil and urea in solid-state studied by terahertz time-domain spectroscopy

Using terahertz time-domain spectroscopy characterization, we observe that urea is able to recognize and interact with uracil efficiently even in the solid phase without involving water or solvents. A cocrystal configuration linked by a pair of hydrogen bonds between uracil and urea was formed. The terahertz absorption spectrum of the cocrystal shows a distinct new absorption at 0.8 THz (26.7 cm(-1)), which originates from the intermolecular hydrogen bonding. Both mechanical milling and heating can accelerate the reaction efficiently. Density functional theory was adopted to simulate the vibrational modes of the cocrystal, and the results agree well with the experimental observation. Multiple techniques, including powder X-ray diffraction, scanning electron microscopy, and Fourier transform infrared spectroscopy, were performed to investigate the reaction process, and they presented supportive evidence. This work enables in-depth understanding of recognition and interaction of urea with nucleobases and comprehension of the denaturation related to RNA. We also demonstrate that terahertz spectroscopy is an effective and alternative tool for online measurement and quality control in pharmaceutical and chemical industries.

Proton tunnelling in the hydrogen bonds of the benzoic acid dimer: (18)O substitution and isotope effects of the heavy atom framework

Field-cycling (1)H NMR relaxometry has been used to measure the rate of concerted double proton transfer in the hydrogen bonds of (16)O and (18)O isotopologues of benzoic acid dimers. The experiments have been conducted in the solid state at low temperature 13.3 ≤ T ≤ 80 K where the dynamics are dominated by incoherent proton tunnelling. The low temperature tunnelling rate in the (16)O isotopologue is observed to be approximately 15% faster than in the (18)O isotopologue. The difference is attributed to an isotope effect of the heavy atom framework of the benzoic acid dimer resulting from displacements of the oxygen atoms that accompany the proton transfer. Sources of systematic uncertainty have been minimized in the design of the experimental protocols and the experiments are critically appraised in formally assigning the measured differences to an effect of mass on the tunnelling dynamics.

Scale-estimation of quantum coherent energy transport in multiple-minima systems

A generic and intuitive model for coherent energy transport in multiple minima systems coupled to a quantum mechanical bath is shown. Using a simple spin-boson system, we illustrate how a generic donor-acceptor system can be brought into resonance using a narrow band of vibrational modes, such that the transfer efficiency of an electron-hole pair (exciton) is made arbitrarily high. Coherent transport phenomena in nature are of renewed interest since the discovery that a photon captured by the light-harvesting complex (LHC) in photosynthetic organisms can be conveyed to a chemical reaction centre with near-perfect efficiency. Classical explanations of the transfer use stochastic diffusion to model the hopping motion of a photo-excited exciton. This accounts inadequately for the speed and efficiency of the energy transfer measured in a series of recent landmark experiments. Taking a quantum mechanical perspective can help capture the salient features of the efficient part of that transfer. To show the versatility of the model, we extend it to a multiple minima system comprising seven-sites, reminiscent of the widely studied Fenna-Matthews-Olson (FMO) light-harvesting complex. We show that an idealised transport model for multiple minima coupled to a narrow-band phonon can transport energy with arbitrarily high efficiency.

Kinetics of Coupled Double Proton and Deuteron Transfer in Hydrogen-Bonded Ribbons of Crystalline Pyrazole-4-carboxylic Acid

The proton tautomerism of pyrazole-4-carboxylic acid (PCA) has been studied by a combination of 15N CPMAS and 2H NMR spectroscopy and relaxometry. Down to 250 K, PCA forms a hydrogen bonded ribbon where adjacent carboxylic and pyrazole groups are linked by an OH···N and an O···HN hydrogen bond, forming either the tautomeric state A or B. Down to about 250 K, the tautomerism is fast on the NMR timescale and degenerate, corresponding to a phase exhibiting dynamic proton disorder. At lower temperatures, a transition to an ordered phase is observed with localized protons, assigned to an all-syn conformation adopting the sequence of tautomeric states ..ABABA.. The longitudinal 15N relaxation times T 1 of PCA-15N2 have been measured at 9.12 MHz (2.1 T). Because of the low field, a chemical shift anisotropy mechanism could be neglected, and the data were analyzed in terms of a dipolar 1H-15N relaxation mechanism, yielding the rate constants k HH. The rate constants k HD and k DD were obtained from the measurement and analysis of the 2H T 1 values of PCA-15N1-d0.9 and PCA-15N1-d0.1 measured at 46.03 MHz. Within the margin of error, no kinetic isotope effects could be detected, in contrast to previous results reported for the very fast tautomerism of solid benzoic acid dimers and the much slower tautomerism of solid 3,5-diphenyl-4-brompyrazole (DPBrP) dimers. The Arrhenius curves of all three systems were simulated using the Bell–Limbach tunneling model. Evidence for a major heavy atom motion for the tautomerism of PCA is obtained, associated with small angle reorientation of PCA molecules around the molecular axis. The observed proton order-disorder transition and the mechanism of the observed rate process are discussed in terms of a coupling of adjacent tautomeric states.

A model proton-transfer system in the condensed phase: NH4(+)OOH(-), a crystal with short intermolecular H-bonds

The crystal structure of NH(4)(+)OOH(-) is determined from single-crystal x-ray data obtained at 150 K. The crystal belongs to the space group P2(1)/c and has four molecules in a unit cell. The structure consists of discrete NH(4)(+) and OOH(-) ions. The OOH(-) ions are linked by short hydrogen bonds (2.533 Å) to form parallel infinite chains. The ammonium ions form links between these chains (the N⋯O distances vary from 2.714 to 2.855 Å) giving a three-dimensional network. The harmonic IR spectrum and H-bond energies are computed at the Perdew-Burke-Ernzerhof (PBE)/6-31G(∗∗) level with periodic boundary conditions. A detailed analysis of the shared (bridging) protons' dynamics is obtained from the CPMD simulations at different temperatures. PBE functional with plane-wave basis set (110 Ry) is used. At 10 K the shared proton sits near the oxygen atom, only a few proton jumps along the chain are detected at 70 K while at 270 K numerous proton jumps exist in the trajectory. The local-minimum structure of the space group Cc is localized. It appears as a result of proton transfer along a chain. This process is endothermic (∼2 kJ/mol) and is described as P2(1)/c↔2Cc. The computed IR spectrum at 10 K is close to the harmonic one, the numerous bands appear at 70 K while at 270 K it shows a very broad absorption band that covers frequencies from about 1000 to 3000 cm(-1). The advantages of the NH(4)(+)OOH(-) crystal as a promising model for the experimental and DFT based molecular dynamics simulation studies of proton transfer along the chain are discussed.

Complex methyl group and hydrogen-bonded proton motions in terms of the Arrhenius and Schrödinger equations

Equations for the temperature dependence of the spectral densities J(is)(m)(momega(I) +/-omega(T)), where m=1, 2, omega(I) and omega(T) are the resonance and tunnel splitting angular frequencies, in the presence of a complex motion, have been derived. The spin pairs of the protons or deuterons of the methyl group perform a complex motion consisting of three component motions. Two of them involve mass transportation over the barrier and through the barrier. They are characterized by k((H)) (Arrhenius) and k((T)) (Schrödinger) rate constants, respectively. The third motion causes fluctuations of the frequencies (nomega(I)+/-omega(T)) and it is related to the lifetime of the methyl spin at the energy level influenced by the rotor-bath interactions. These interactions induce rapid transitions, changing the symmetry of the torsional sublevels either from A to E or from E(a) to E(b). The correlation function for this third motion (k((omega)) rate constant) has been proposed by Müller-Warmuth et al. The spectral densities of the methyl group hindered rotation (k((H)), k((T)) and k((omega)) rate constants) differ from the spectral densities of the proton transfer (k((H)) and k((T)) rate constants) because three compound motions contribute to the complex motion of the methyl group. The recently derived equation [Formula: see text] , where [Formula: see text] and [Formula: see text] are the fraction and energy of particles with energies from zero to E(H), is taken into account in the calculations of the spectral densities. This equation follows from Maxwell's distribution of thermal energy. The spectral densities derived are applied to analyse the experimental temperature dependencies of proton and deuteron spin-lattice relaxation rate in solids containing the methyl group. A wide range of temperatures from zero Kelvin up to the melting point is considered. It has been established that the motion characterized by k((omega)) influences the spin-lattice relaxation up to the temperature T(tun) only. This temperature is directly determined by the equation C(p)T=E(H) (thermal energy=activation energy), where C(p) is the molar heat capacity. Probably the cessation of the third motion is a result of the de Broglie wavelength related to this motion becoming too short. As shown recently, the potential barrier can be an obstacle for the de Broglie wave. The theoretical equations derived in this paper are compared to those known in the literature.

Tunneling molecular dynamics in the light of the corpuscular-wave dualism theory

This paper presents the experimental demonstration of the corpuscular-wave dualism theory. The correlation between the de Broglie wavelength related to the thermal motion and the potential barrier width and height is reported. The stochastic jumps of light atoms (hydrogen, deuterium) between two equilibrium sites A and B (identical geometry) occur via different pathways; one pathway is over the barrier (classical dynamics), and the other one is through the barrier (tunneling). On the over-the-barrier pathway, there are no obstacles for the de Broglie waves, and this pathway exists from high to low temperatures up to 0 K because the thermal energy is subjected to the Maxwell distribution and a certain number of particles owns enough energy for the hopping over the barrier. On the tunneling pathway, the particles pass through the barrier, or they are reflected from the barrier. Only particles with the energy lower than barrier heights are able to perform a tunneling hopping. The de Broglie waves related to these energies are longer than the barrier width. The Schrödinger equation is applied to calculate the rate constant of tunneling dynamics. The Maxwell distribution of the thermal energy has been taken into account to calculate the tunneling rate constant. The equations for the total spectral density of complex motion derived earlier by us together with the expression for the tunneling rate constant, derived in the present paper, are used in analysis of the temperature dependence of deuteron spin-lattice relaxation of the ammonium ion in the deuterated analogue of ammonium hexachloroplumbate ((ND4)2PbCl6). It has been established that the equation CpTtun = EH (thermal energy equals activation energy), where Cp is the molar heat capacity (temperature-dependent, known from literature), determines directly the low temperature Ttun at which the de Broglie wavelength, lambdadeBroglie, related to the thermal energy, CpT, is equal to the potential barrier width, L. Above Ttun, the lambdadeBroglie wavelength related to the CpT energy is shorter than the potential barrier width and not able to overcome the barrier. The activation energy EH equals 7.5 kJ/mol, and therefore, the Ttun temperature for deuterons in ((ND4)2PbCl6 is 55.7 K. The agreement between the potential barrier width following from the simple geometrical calculations (L = 0.722 A) and de Broglie wavelength at Ttun (L = 0.752 A) is good. The temperature plots of the deuteron correlation times for (ND4)2PbCl6 reveal comparable values of the correlation times of the tunneling, (tau(T)), and over-the-barrier jumps (tau(H)) near 34.8 K. Matsuo, on the basis of the molar heat capacity study, found the first-order phase transition at this temperature.

NMR studies of double proton transfer in hydrogen bonded cyclic N,N'-diarylformamidine dimers: conformational control, kinetic HH/HD/DD isotope effects and tunneling

Using dynamic NMR spectroscopy, the kinetics of the degenerate double proton transfer in cyclic dimers of polycrystalline (15)N,(15)N'-di-(4-bromophenyl)-formamidine (DBrFA) have been studied including the kinetic HH/HD/DD isotope effects in a wide temperature range. This transfer is controlled by intermolecular interactions, which in turn are controlled by the molecular conformation and hence the molecular structure. At low temperatures, rate constants were determined by line shape analysis of (15)N NMR spectra obtained using cross-polarization (CP) and magic angle spinning (MAS). At higher temperatures, in the microsecond time scale, rate constants and kinetic isotope effects were obtained by a combination of longitudinal (15)N and (2)H relaxation measurements. (15)N CPMAS line shape analysis was also employed to study the non-degenerate double proton transfer of polycrystalline (15)N,(15)N'-diphenyl-formamidine (DPFA). The kinetic results are in excellent agreement with the kinetics of DPFA and (15)N,(15)N'-di-(4-fluorophenyl)-formamidine (DFFA) studied previously for solutions in tetrahydrofuran. Two large HH/HD and HD/DD isotope effects are observed in the whole temperature range which indicates a concerted double proton transfer mechanism in the domain of the reaction energy surface. The Arrhenius curves are non-linear indicating a tunneling mechanism. Arrhenius curve simulations were performed using the Bell-Limbach tunneling model. The role of the phenyl group conformation and hydrogen bond compression on the barrier of the proton transfer is discussed.

Application of Schrödinger Equation to Study the Tunnelling Dynamics of Proton Transfer in the Hydrogen Bond of 2,5-Dinitrobenzoic Acid: Proton T1, T1ρ, and Deuteron T1 Relaxation Methods

Temperature measurements of proton T1 (24.7 MHz), deuteron (deuterated hydroxyl group) T1 (55.2 MHz), and proton T1(rho) (B1 = 9 G) spin-lattice relaxation times of 2,5-dinitrobenzoic acid have been performed. An analysis of present experimental data together with previously published proton T1 (55.2 MHz) data has revealed the following molecular motions: proton/deuteron transfer in the hydrogen bond and two-site hopping of the whole dimer. It is shown that the proton-transfer dynamics are characterized by two correlation times tau(ov) and tau(tu), describing two fundamentally different motional processes, namely, thermally activated jumps over the barrier and tunneling through the barrier. The temperature dependence of 1/tau(tu) is the solution of Schrödinger's equation, which also yields the temperature T(tun), where begins the tunnel pathway for proton transfer. A new equation for the spectral density function of complex motion consisting of the three motions is derived. The third motion (two-site hopping of the whole dimer characterized by tau(lib) correlation time) is responsible for a proton T1(rho) minimum in high temperatures, just below the melting point. Such a minimum is not reached by T1 temperature dependencies. The minimum of T1(rho) assigned to the classical hopping of a hydrogen-bonded proton occurs in the same low-temperature regime in which the flattening of the temperature dependencies of T1 points to the dominance of incoherent tunneling. This experimental fact denies the known theories predicting the intermediate temperature regime where a smooth transition between classical and quantum tunneling dynamics is expected. The fit of the derived theoretical equations to the experimental data T1(rho) and T1 is satisfactory. The correlation times obtained for deuterons indicate deuteron-transfer dynamics much slower than proton-transfer dynamics. It is concluded that the classical proton transfer takes place over the whole temperature regime, while the incoherent tunneling occurs below 46.5 (hydrogen) or 87.2 K (deuterium) only.

Ab initio quantum chemical investigation of the ground and excited states of salicylic acid dimer

The structure and stability of different forms of salicylic acid dimer have been examined by Hartree-Fock and density functional theoretic calculations using 6-31G(d,p) and 6-311++g(d,p) basis sets. Vertical excitation energies for the monomer as well as the dimer have been computed using the time-dependent density functional theory using 6-311++G(d,p) basis set. The predicted absorption maxima for the first excited singlet state of salicylic acid monomer and the dimer of the primary form are in reasonable agreement with the experimental result. There is a slight red shift (approximately 6 nm) in the absorption maximum in going from the monomer to the dimer, in accord with the experimental observation. Configuration-interaction calculations including single excitation have been carried out to map the potential-energy profile for the intra- as well as the intermolecular proton transfer in different forms of the dimer. The barrier for proton transfer in the ground state as well as the excited states makes it clear that most of the processes take place in the primary form and largely by intramolecular proton transfer.

QTAIM study of strong H-bonds with the O-H...a fragment (A=O, N) in three-dimensional periodical crystals

The relationship between the d(H...A) distance (A=O, N) and the topological properties at the H...A bond critical point of 37 strong (short) hydrogen bonds occurring in 26 molecular crystals are analyzed using the quantum theory of atoms in molecules (QTAIM). Ground-state wave functions of the three-dimensional periodical structures representing the accurate experimental geometries calculated at the B3LYP/6-31G** level of approximation were used to obtain the QTAIM electron density characteristics. The use of an electron-correlated method allowed us to reach the quantitatively correct values of electron density rhob at the H...A bond critical point. However, quite significant differences can appear for small absolute values of the Laplacian (<0.5 au). The difference between the H...O and H...N interactions is described using the rhob versus d(H...A) dependence. It is demonstrated that the values of parameters in this dependence are defined by the nature of the heavy atom forming the H...A bond. An intermediate (or transit) region separating the shared and closed-shell interactions is observed for the H-bonded crystals in which the bridging proton can move from one heavy atom to another. The crystalline environment changes the location of the bridging proton in strong H-bonded systems; however, the d(O-H)/d(H...O) ratio is approximately the same for both the gas-phase complexes and molecular crystals with a linear or near-linear O-H...O bond.

1H/2H NMR Studies of Geometric H/D Isotope Effects on the Coupled Hydrogen Bonds in Porphycene Derivatives

The 1H and 2H NMR spectra of porphycene (1), 2,3,6,7,12,13,16,17-octaethylporphycene (2), 2,7,12,17-tetra-n-propylporphycene (3), and 2,7,12,17-tetra-(tert-butyl)-3,6-13,16-dibenzo[cde;mno]porphycene (4) partially deuterated in the mobile proton sites are reported. These compounds exhibit two intramolecular NHN hydrogen bonds of increasing strength representing models of the concerted HH transfer in the parent compound, porphycene. The 1H chemical shifts of the mobile protons are correlated with the difference of the energies of the amino- and imino-N1s orbitals reported by Ghosh A.; Moulder J.; Bröring M.; Vogel E. Angew. Chem., Int. Ed. 2001, 113, 445-448. The chemical shifts of 4 indicate a reduced contribution of the aromatic ring current as compared to the other compounds which is associated to the nonplanarity of this molecule. The primary H/D isotope effects on the chemical shifts give information about the primary, secondary, and vicinal geometric isotope effects of the two inner hydrogen bonds of porphycenes. The vicinal effects indicate a cooperative coupling of the two hydrogen bonds which may favor a concerted double proton-transfer mechanism.

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.

Evidence of conformational exchange averaging in the thermal rotational spectrum of ethyl cyanoformate

The Stark modulated low resolution microwave spectrum of ethyl cyanoformate between 21.5 and 24.0 GHz at 210, 300, and 358 K, which shows the J + 1 <-- J = 8 <-- 7 bands of three species, is compared to simulations based on electronic structure calculations at the MP2/6-311++G theory level. Calculations at this theory level reproduce the relative energies of the syn-anti and syn-gauche conformers, obtained in a previous study, and indicate that the barrier to conformer exchange is approximately 360 cm(-1) higher in energy than the syn-anti minimum. Simulated spectra of the eigenstates of the calculated O-ethyl torsional potential function reproduce the relative intensities and shapes of the lower and higher frequency bands which correspond to transitions of the syn-anti and syn-gauche conformers, respectively, but fail to reproduce the intense center band in the experimental spectra. A model incorporating exchange averaging reproduces the intensity of the center band and its temperature dependence. These simulations indicate that a large fraction of the thermal population at all three temperatures undergoes conformational exchange with an average energy specific rate constant, <k(E)>, of approximately 25 GHz. This model can explain anomalies present in rotational spectra of many other compounds composed of mixtures of conformers.

Fingerprints of Damped Quantum Rotation Observed in Solid-State Proton NMR Spectra

(1)H NMR spectra of the methyl group in an oriented crystal sample of methylmalonic acid with all three non-methyl protons replaced by deuterons are interpreted in terms of the damped quantum rotation (DQR) theory of NMR line shapes. The DQR approach offers a perfect theoretical reproduction of the observed spectra while the conventional Alexander-Binsch line-shape model shows evident defects in the present case. The temperature trends of the quantities characterizing the coherent and incoherent dynamics of the methyl group in the DQR approach (the effective tunnelling frequency and two coherence-damping rates) derived from the spectra are fairly reproduced using a model reported previously. The present findings provide further evidence of limitations to the validity of the common belief that molecular rate processes in condensed phases are necessarily classical.

1,4- vs 1,3-prototropic mechanism for intramolecular double proton transfer reaction in monothiooxalic acid. Theoretical investigation of potential energy surface

Proton transfer is a very common and important chemical step in many systems. Despite its apparent simplicity, a correct description of this chemical process is difficult from a theoretical point of view. It requires a correct and simultaneous description of a bond breaking and a bond formation. The situation is even much more complicated when two protons are implied. This is the case for monothiooxalic acid, for which two different types (1,3- and 1,4-prototropy) of proton transfers can be invoked. A further problem is the type of the reaction (concerted or not). This paper reports a complete investigation of the potential energy surfaces: characterization of equilibrium points and transitions states. The main conclusion is: the 1,4-prototropic mechanism, mainly considered as a one step concerted exchange of protons, is the most favored from an energetic point of view.

Probing Photophysical and Photochemical Processes of Benzoic Acid from ab Initio Calculations

The CASSCF and DFT methods have been used to determine geometric and electronic structures of the benzoic acid monomer in the S(0), S(1), S(2), T(1), and T(2) electronic states. The S(1)/T(2)/T(1) three-surface intersection was found by the state-averaged CASSCF calculations, which, in combination with features of the five lowest electronic states, provides new insights into photophysical processes of the benzoic acid monomer. The potential energy profiles of the alpha C-C and C-O bond fissions as well as decarboxylation reaction in different electronic states have been determined for the benzoic acid monomer. The alpha C-O bond cleavage starts from the T(2) state and leads to the fragments of C(6)H(5)CO(X(2)A') and OH(X(2)Pi) in the ground state, which is predicted to be the most possible channel upon photoexcitation of the benzoic acid monomer at 270 nm or shorter wavelengths.

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.

NMR Relaxation Study of the Protonated Form of 1,8-Bis(dimethylamino)naphthalene in Isotropic Solution: Anisotropic Motion outside of Extreme Narrowing and Ultrafast Proton Transfer

The protonated form of 1,8-bis(dimethylamino)naphthalene (DMANH(+)) consists of a rigid, aromatic framework, substituted by two amino groups that are connected by a strong, symmetric (on the NMR time-scale) hydrogen bond bridge. The reorientational motion of the molecule in dimethylformamide-d7 solution was characterized by T(1) and NOE measurements for aromatic (13)C nuclei. Treating the reorientation of DMANH(+) as anisotropic rotational diffusion of a rigid body, the diffusion tensor was determined with good accuracy. Measurements and interpretation of (15)N T(1) and NOE indicate that the proton transfer between potential minima in the hydrogen bond bridge is faster than the molecular reorientation.

Dynamic NMR Study of the Mechanisms of Double, Triple, and Quadruple Proton and Deuteron Transfer in Cyclic Hydrogen Bonded Solids of Pyrazole Derivatives

Using dynamic solid state (15)N CPMAS NMR spectroscopy (CP = cross polarization, MAS = magic angle spinning), the kinetics of the degenerate intermolecular double and quadruple proton and deuteron transfers in the cyclic dimer of (15)N labeled polycrystalline 3,5-diphenyl-4-bromopyrazole (DPBrP) and in the cyclic tetramer of (15)N labeled polycrystalline 3,5-diphenylpyrazole (DPP) have been studied in a wide temperature range at different deuterium fractions in the mobile proton sites. Rate constants were measured on a millisecond time scale by line shape analysis of the doubly (15)N labeled compounds, and by magnetization transfer experiments on a second timescale of the singly (15)N labeled compounds in order to minimize the effects of proton-driven (15)N spin diffusion. For DPBrP the multiple kinetic HH/HD/DD isotope effects could be directly obtained. By contrast, four rate constants k(1) to k(4) were obtained for DPP at different deuterium fractions. Whereas k(1) corresponds to the rate constant k(HHHH) of the HHHH isotopolog, an appropriate kinetic reaction model was needed for the kinetic assignment of the other rate constants. Using the model described by Limbach, H. H.; Klein, O.; Lopez Del Amo, J. M.; Elguero, J. Z. Phys. Chem. 2004,218, 17, a concerted quadruple proton-transfer mechanism as well as a stepwise consecutive single transfer mechanism could be excluded. By contrast, using the kinetic assignment k(2) approximately k(3) approximately k(HHHD) approximately k(HDHD) and k(3) approximately k(HDDD) approximately k(DDDD), the results could be explained in terms of a two-step process involving a zwitterionic intermediate. In this mechanism, each reaction step involves the concerted transfer of two hydrons, giving rise to primary kinetic HH/HD/DD isotope effects, whereas the nontransferred hydrons only contribute small secondary effects, which are not resolved experimentally. By contrast, the multiple kinetic isotope effects of the double proton transfer in DPBrP and of the triple proton proton transfer in cyclic pyrazole trimers studied previously indicate concerted transfer processes. Thus, between n = 3 and 4 a switch of the reaction mechanism takes place. This switch is rationalized in terms of hydrogen bond compression effects associated with the multiple proton transfers. The Arrhenius curves of all processes are nonlinear and indicate tunneling processes at low temperatures. In a preliminary analysis, they are modeled in terms of the Bell-Limbach tunneling model.