Concerted Proton Tunneling Research Articles

Significance Of Nuclear Quantum Effects In Hydrogen Bonded Molecular Chains.

In hydrogen-bonded systems, nuclear quantum effects such as zero-point motion and tunneling can significantly affect their material properties through underlying physical and chemical processes. Presently, direct observation of the influence of nuclear quantum effects on the strength of hydrogen bonds with resulting structural and electronic implications remains elusive, leaving opportunities for deeper understanding to harness their fascinating properties. We studied hydrogen-bonded one-dimensional quinonediimine molecular networks which may adopt two isomeric electronic configurations via proton transfer. Herein, we demonstrate that concerted proton transfer promotes a delocalization of π-electrons along the molecular chain, which enhances the cohesive energy between molecular units, increasing the mechanical stability of the chain and giving rise to distinctive electronic in-gap states localized at the ends. These findings demonstrate the identification of a class of isomeric hydrogen-bonded molecular systems where nuclear quantum effects play a dominant role in establishing their chemical and physical properties. This identification is a step toward the control of mechanical and electronic properties of low-dimensional molecular materials via concerted proton tunneling.

Ice Ih revisited: No proton tunneling observed in a quasielastic neutron scattering experiment

A large broadening ($\ensuremath{\sim}0.4\phantom{\rule{0.16em}{0ex}}\mathrm{meV}$) of quasielastic neutron scattering (QENS) signal in ${\mathrm{H}}_{2}\mathrm{O}$ ice ${I}_{h}$ at $T=5\phantom{\rule{0.16em}{0ex}}\mathrm{K}$ was observed by Bove et al. [Phys. Rev. Lett. 103, 165901 (2009)] and explained by a model of concerted proton tunneling. This result was rather unexpected, as prior studies never showed significant mobility in water at low temperatures and ambient pressure. There were a few attempts of theoretical understanding of the effect. Recent path-integral simulations as well as quantum lattice-gauge theory supported the possibility of the collective tunneling of protons in ice ${I}_{h}$, however other studies stated that concerted tunneling in ice ${I}_{h}$ should have very low frequency. Here, we report on QENS measurements of ${\mathrm{H}}_{2}\mathrm{O}$ ice ${I}_{h}$ at 1.8 and 5 K by using neutron scattering spectrometers with the energy resolution similar to and four times better than the energy resolution in the original experiment of Bove and co-workers. We did not observe any QENS broadening, and the measured spectra for the ice ${I}_{h}$ and the reference vanadium sample were almost identical. Therefore, we conclude that there is no proton tunneling in ice ${I}_{h}$ at temperatures down to 1.8 K measurable on an energy scale of 3.5 $\ensuremath{\mu}\mathrm{eV}$ and above. The literature data on low-temperature heat capacity of ice ${I}_{h}$ support this conclusion.

Microwave measurements of the tropolone-formic acid doubly hydrogen bonded dimer.

The microwave spectrum was measured for the doubly hydrogen bonded dimer formed between tropolone and formic acid. The predicted symmetry of this dimer was C2v(M), and it was expected that the concerted proton tunneling motion would be observed. After measuring 25 a- and b-type rotational transitions, no splittings which could be associated with a concerted double proton tunneling motion were observed. The calculated barrier to the proton tunneling motion is near 15,000 cm(-1), which would likely make the tunneling frequencies too small to observe in the microwave spectra. The rotational and centrifugal distortion constants determined from the measured transitions were A = 2180.7186(98) MHz, B = 470.873 90(25) MHz, C = 387.689 84(22) MHz, DJ = 0.0100(14) kHz, DJK = 0.102(28) kHz, and DK = 13.2(81) kHz. The B3LYP/aug-cc-pVTZ calculated rotational constants were within 1% of the experimentally determined values.

Classification of hydrogen bond flips in small water polyhedra applied to concerted proton tunneling.

Recently a new mechanism of proton tunneling in a prism-like water hexamer was revealed [Richardson et al., Science, 2016, 351, 1310]. The tunneling motion involves the concerted breaking of two hydrogen bonds and rotations of two nearest water molecules. Eventually, this structural transformation means flipping one of the hydrogen bonds without the creation of defects in the hydrogen bond network. On the surface of polyhedral water clusters, there are five essentially different types of hydrogen bonds, and only two of them can be changed in this manner. In this article, the topological classification of such transformations for five small water polyhedra: triangular, pentagonal, and hexagonal prisms as well as cube and polyhedron 4454, consisting of four square and four pentagonal faces, is presented. Our classification includes the enumeration of all possible one-bond-flips with consideration of the types of hydrogen bonds on the polyhedral surface. Attention is paid to the most stable proton configurations which can be studied in experiments. It was established that a number of one-bond-flip transitions between the low energy configurations are possible in clusters in the shape of triangular and pentagonal prisms.

Pressure-induced localisation of the hydrogen-bond network in KOH-VI.

Using a combination of ab initio crystal structure prediction and neutron diffraction techniques, we have solved the full structure of KOH-VI at 7 GPa. Rather than being orthorhombic and proton-ordered as had previously be proposed, we find that this high-pressure phase of potassium hydroxide is tetragonal (space group I4/mmm) and proton disordered. It has an unusual hydrogen bond topology, where the hydroxyl groups form isolated hydrogen-bonded square planar (OH)4 units. This structure is stable above 6.5 GPa and, despite being macroscopically proton-disordered, local ice rules enforce microscopic order of the hydrogen bonds. We suggest the use of this novel type of structure to study concerted proton tunneling in the solid state, while the topology of the hydrogen bond network could conceivably be exploited in data storage applications based solely on the manipulations of hydrogen bonds. The unusual localisation of the hydrogen bond network under applied pressure is found to be favored by a more compact packing of the constituents in a distorted cesium chloride structure.

Calculations and measurements of the deuterium tunneling frequency in the propiolic acid-formic acid dimer and description of a newly constructed Fourier transform microwave spectrometer

The concerted proton tunneling frequency for the propiolic acid-formic acid dimer was calculated using a relaxed ab initio double-well potential in the imaginary-frequency mode of the saddle point, and new measurements were made for the deuterated propiolic acid-formic acid (ProOD-FAOD) isotopologue. It is important to have consistent calculated tunneling frequency values between normal and deuterated isotopologues since parameters can be readily adjusted to get good agreement with one isotopologue. High-resolution rotational spectra of deuterated (ProOD-FAOD) dimer were measured using a newly constructed Fourier Transform microwave spectrometer. The new spectrometer has mirror size: 30 cm in diameter with a radius of curvature of 59 cm and is equipped with multiple-FID data collection (5-10 FID's for each gas pulse). For the deuterated (ProOD-FAOD) isotopologue, 45 rotational lines (a type: 34; b type: 11) were measured in the lowest tunneling states range between 6.5 GHz and 15.5 GHz. With the new high-resolution measurements of the tunneling doublets (b-dipole transitions), the double potential well responsible for the deuterium tunneling was depicted much more precisely. The two tunneling states are separated by 3.48 MHz. The rotational constants obtained in this work are quite helpful for further structure analysis as well.

Microwave measurements of proton tunneling and structural parameters for the propiolic acid–formic acid dimer

Microwave spectra of the propiolic acid-formic acid doubly hydrogen bonded complex were measured in the 1 GHz to 21 GHz range using four different Fourier transform spectrometers. Rotational spectra for seven isotopologues were obtained. For the parent isotopologue, a total of 138 a-dipole transitions and 28 b-dipole transitions were measured for which the a-dipole transitions exhibited splittings of a few MHz into pairs of lines and the b-type dipole transitions were split by ~580 MHz. The transitions assigned to this complex were fit to obtain rotational and distortion constants for both tunneling levels: A(0+) = 6005.289(8), B(0+) = 930.553(8), C(0+) = 803.9948(6) MHz, Δ(0+)(J) = 0.075(1), Δ(0+)(JK) = 0.71(1), and δ(0+)(j) = -0.010(1) kHz and A(0-) = 6005.275(8), B(0-) = 930.546(8), C(0-) = 803.9907(5) MHz, Δ(0-)(J) = 0.076(1), Δ(0-)(JK) = 0.70(2), and δ(0-)(j) = -0.008(1) kHz. Double resonance experiments were used on some transitions to verify assignments and to obtain splittings for cases when the b-dipole transitions were difficult to measure. The experimental difference in energy between the two tunneling states is 291.428(5) MHz for proton-proton exchange and 3.35(2) MHz for the deuterium-deuterium exchange. The vibration-rotation coupling constant between the two levels, F(ab), is 120.7(2) MHz for the proton-proton exchange. With one deuterium atom substituted in either of the hydrogen-bonding protons, the tunneling splittings were not observed for a-dipole transitions, supporting the assignment of the splitting to the concerted proton tunneling motion. The spectra were obtained using three Flygare-Balle type spectrometers and one chirped-pulse machine at the University of Virginia. Rotational constants and centrifugal distortion constants were obtained for HCOOH···HOOCCCH, H(13)COOH···HOOCCCH, HCOOD···HOOCCCH, HCOOH···DOOCCCH, HCOOD···DOOCCCH, DCOOH···HOOCCCH, and DCOOD···HOOCCCH. High-level ab initio calculations provided initial rotational constants for the complex, structural parameters, and some details of the proton tunneling potential energy surface. A least squares fit to the isotopic data reveals a planar structure that is slightly asymmetric in the OH distances. The formic OH···O propiolic hydrogen bond length is 1.8 Å and the propiolic OH···O formic hydrogen bond length is 1.6 Å, for the equilibrium configuration. The magnitude of the dipole moment was experimentally determined to be 1.95(3) × 10(-30) C m (0.584(8) D) for the 0(+) states and 1.92(5) × 10(-30) C m (0.576(14) D) for the 0(-) states.

Anomalous Proton Dynamics in Ice at Low Temperatures

We present incoherent quasielastic neutron scattering measurements on ice Ih (ordinary ice) and Ic (cubic ice) which show the existence of nonharmonic motion of hydrogen at low temperatures, down to 5 K. We show that this dynamics is localized, nonvibrational, and related to the hydrogen disorder since it is absent in ordered ice VIII. A main jump distance of 0.75 A is identified, hence close to the distance between the two possible proton sites along the oxygen-oxygen bond. The dynamics is non-Arrhenius, has a large time rate of 2.7x10(11) s-1, and affects only a few percent of the total number of hydrogen atoms in the crystal. These results give evidence for the existence of concerted proton tunneling in these ice phases.

Vibrational spectrum of the formic acid dimer in the OH stretch region. A model 3D study

The vibrational spectrum of the formic acid dimer in the OH stretch region is calculated by solving a 3D vibrational Schrödinger equation and employing an electronic potential and dipole moment calculated at the MP2 level of theory. Two Raman-active bands of the symmetric OH stretch fundamental are located around 2600 cm −1 ; the gap between them (∼ 300 cm −1 ) is caused by concerted proton tunneling. IR-active bands around 3100 cm −1 are assigned to the antisymmetric OH stretch fundamental, its first overtone and two complex transitions. The calculated H/ D isotope effect on the intensity of the antisymmetric OH fundamental is 1.92.

Vibrational spectra of the methanol tetramer in the OH stretch region. Two cyclic isomers and concerted proton tunneling

Second order Mo/ller–Plesset perturbation theory and density functional theory are employed to localize several stationary points on the potential energy surface of the cyclic methanol tetramer. Two cyclic isomers are identified: one of S4 symmetry, with methyl groups in up-down-up-down configuration, and a second one of Ci symmetry, with the methyl groups in up-up-down-down configuration. The latter minimum is 360 cm−1 above the S4 minimum, with a barrier of 475 cm−1 separating them. These isomers give rise to four asymmetric OH modes around 3300 cm−1. A model of the concerted proton transfer, S4→D2d→S4 (D2d is the transition structure), yields an estimate of 0.7 cm−1 for the tunneling splitting of the totally symmetric OH stretch vibrational fundamental. Raman spectra would show evidence of this fundamental and help to identify admixtures of the Ci isomer.

New Insights into the Dynamics of Concerted Proton Tunneling in Cyclic Water and Hydrogen Fluoride Clusters

The concerted proton transfer hypersurface of cyclic water and hydrogen fluoride clusters has been described by high-level molecular quantum mechanical calculations. For the cyclic water clusters the concerted proton transfer transition states have been investigated for the first time with methods including treatment of dynamic electron correlation. The crucial importance of dynamic electron correlation for the barrier heights is demonstrated. A detailed analysis of the minimum energy path has been performed. The reaction swath of the concerted proton transfer was examined indicating a reasonable description by harmonic approximation of the energy hypersurface. Transfer rates have been calculated by means of variational transition state theory with interpolated corrections (VTST-IC) and dual-level direct dynamics (DLDD), both with semiclassical tunneling corrections. Tunneling is very efficient in the concerted proton exchange reaction of the cyclic hydrogen-bonded clusters under investigation. Rate const...