Hydrogen Bonds with Large Proton Polarizability and Proton Transfer Processes in Electrochemistry and Biology

Stabilization of hydrogen-bonded molecular chains by carbon nanotubes

We study numerically nonlinear dynamics of several types of molecular systems composed of hydrogen-bonded chains placed inside carbon nanotubes with open edges. We demonstrate that carbon nanotubes provide a stabilization mechanism for quasi-one-dimensional molecular chains via the formation of their secondary structures. In particular, a polypeptide chain (Gly)N placed inside a carbon nanotube can form a stable helical chain (310-, α-, π-, and β-helix) with parallel chains of hydrogen-bonded peptide groups. A chain of hydrogen fluoride molecules ⋯FH⋯FH⋯FH can form a hydrogen-bonded zigzag chain. Remarkably, we demonstrate that for molecular complexes (Gly)N∈CNT and (FH)N∈CNT, the hydrogen-bonded chains will remain stable even at T=500 K. Thus, our results suggest that the use of carbon nanotubes with encapsulated hydrogen fluoride molecules may be important for the realization of high proton conductivity at high temperatures.

Correlations between the Structures and Spectra of Protonated Water Clusters

Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.

Electronic Fingerprint of the Protonated Imidazole Dimer Probed by X-ray Absorption Spectroscopy

Protons in low-barrier superstrong hydrogen bonds are typically delocalized between two electronegative atoms. Conventional methods to characterize such superstrong hydrogen bonds are vibrational spectroscopy and diffraction techniques. We introduce soft X-ray spectroscopy to uncover the electronic fingerprints for proton sharing in the protonated imidazole dimer, a prototypical building block enabling effective proton transport in biology and high-temperature fuel cells. Using nitrogen core excitations as a sensitive probe for the protonation status, we identify the X-ray signature of a shared proton in the solvated imidazole dimer in a combined experimental and theoretical approach. The degree of proton sharing is examined as a function of structural variations that modify the shape of the low-barrier potential in the superstrong hydrogen bond. We conclude by showing how the sensitivity to the quantum distribution of proton motion in the double-well potential is reflected in the spectral signature of the shared proton.

ZundEig: The Structure of the Proton in Liquid Water from Unsupervised Learning

The structure of the excess proton in liquid water has been the subject of lively debate on both experimental and theoretical fronts for the last century. Fluctuations of the proton are typically interpreted in terms of limiting states referred to as the Eigen and Zundel species. Here, we put these ideas under the microscope, taking advantage of recent advances in unsupervised learning that use local atomic descriptors to characterize environments of acidic water combined with advanced clustering techniques. Our agnostic approach leads to the observation of only one charged cluster and two neutral ones. We demonstrate that the charged cluster involving the excess proton is best seen as an ionic topological defect in water's hydrogen bond network, forming a single local minimum on the global free-energy landscape. This charged defect is a highly fluxional moiety, where the idealized Eigen and Zundel species are neither limiting configurations nor distinct thermodynamic states. Instead, the ionic defect enhances the presence of neutral water defects through strong interactions with the network. We dub the combination of the charged and neutral defect clusters as ZundEig, demonstrating that the fluctuations between these local environments provide a general framework for rationalizing more descriptive notions of the proton in the existing literature.

Isotope Effects in the Zundel-Eigen Isomerization of H+(H2O)6

The isomerization pathway between the energetically low-lying Zundel and Eigen isomers of the protonated water hexamer was investigated using high-level ab initio calculations including a treatment of zero-point corrections. On the basis of these calculations, the Zundel-Eigen isomerization was found to proceed through a stable intermediate isomer, which consists of a four-membered ring with two single acceptor water molecules. The inclusion of vibrational zero-point energy is shown to be important for accurately establishing the relative energies of the three relevant isomers involved in the Zundel-Eigen isomerization. Diffusion Monte Carlo calculations including anharmonic vibrational effects show that all three isomers of H⁺(H₂O)₆ and D⁺(D₂O)₆ have well-defined structures. The energetic ordering of the three isomers changes upon deuteration. The implications of these results for the vibrational spectra of H⁺(H₂O)₆ and D⁺(D₂O)₆ are also discussed.

On the Properties of Nafion Membranes Recast from Dispersion in N-Methyl-2-Pyrrolidone

Perfluorosulfonic acid Nafion membranes are widely used as an electrolyte in electrolysis processes and in fuel cells. Changing the preparation and pretreatment conditions of Nafion membranes allows for the optimization of their properties. In this work, a Nafion-NMP membrane with a higher conductivity than the commercial Nafion^® 212 membrane (11.5 and 8.7 mS∙cm^-1 in contact with water at t = 30 °C) and a comparable hydrogen permeability was obtained by casting from a Nafion dispersion in N-methyl-2-pyrrolidone. Since the ion-exchange capacity and the water uptake of these membranes are similar, it can be assumed that the increase in conductivity is the result of optimizing the Nafion-NMP microstructure by improving the connectivity of the pores and channels system. This leads to a 27% increase in the capacity of the membrane electrode assembly with the Nafion-NMP membrane compared to the Nafion^® 212 membrane. Thus, the method of obtaining a Nafion membrane has a great influence on its properties and performance of fuel cells based on them.

From Local Covalent Bonding to Extended Electric Field Interactions in Proton Hydration

Seemingly simple yet surprisingly difficult to probe, excess protons in water constitute complex quantum objects with strong interactions with the extended and dynamically changing hydrogen-bonding network of the liquid. Proton hydration plays pivotal roles in energy transport in hydrogen fuel cells and signal transduction in transmembrane proteins. While geometries and stoichiometry have been widely addressed in both experiment and theory, the electronic structure of these specific hydrated proton complexes has remained elusive. Here we show, layer by layer, how utilizing novel flatjet technology for accurate x-ray spectroscopic measurements and combining infrared spectral analysis and calculations, we find orbital-specific markers that distinguish two main electronic-structure effects: Local orbital interactions determine covalent bonding between the proton and neigbouring water molecules, while orbital-energy shifts measure the strength of the extended electric field of the proton.

Vibrational Dynamics of the Intramolecular H-Bond in Acetylacetone Investigated with Transient and 2D IR Spectroscopy

Acetylacetone (AcAc) has proven to be a fruitful but highly challenging model system for the experimental and computational interrogation of strong intramolecular hydrogen bonds. Key questions remain, however, regarding the identity of the minimum-energy structure of AcAc and the dynamics of intramolecular proton transfer. Here, we investigate the OH/OD stretch and bend regions of the enol tautomer of AcAc and its deuterated isotopologue with transient absorption and 2D IR spectroscopy. The OH bend region reveals a single dominant diagonal transition near 1625 cm^-1 with intense cross peaks to lower-frequency modes, demonstrating highly mixed fingerprint transitions that contain OH bend character. The anharmonic coupling of the OH bend results in a highly elongated OH bend excited-state absorption transition that indicates a large manifold of OH bend overtone/combination bands in the OH stretch region that leads to strong bend-stretch Fermi resonance interactions. The OH and OD stretch regions consist of broad ground-state bleach signals, but there is no clear evidence of ω₂₁ excited-state absorptions due to rapid population relaxation arising from strong intramolecular coupling to bending, fingerprint, and low-frequency H-bond modes. Orientational relaxation dynamics persist for timescales longer than the vibrational lifetimes, with polarization anisotropy components decaying within approximately 2 and 10 periods of the O-O oscillation for the OH and OD stretch, respectively. The significant isotopic dependence of the orientational dynamics is discussed in the context of intramolecular mode coupling, diffusional processes, and contributions from proton/deuteron transfer dynamics.

Proton Transport and Related Chemical Processes of Ice

Excess protons play a key role in the chemical reactions of ice because of their exceptional mobility, even when the diffusion of atoms and molecules is suppressed in ice at low temperatures. This article reviews the current state of knowledge on the properties of excess protons in ice, with a focus on the involvement of protons in chemical reactions. The mechanism of efficient proton transport in ice, which involves a proton-hopping relay along the hydrogen-bond ice network and the reorientation of water, is discussed and compared with the inefficient transport of hydroxide in ice. Distinctly different properties of protons residing in the ice interior and on the ice surface are emphasized. Recent observations of the spontaneous occurrence of reactions in ice at low temperatures, which include the dissociation of protic acids and the hydrolysis of acidic oxides, are discussed with regard to the kinetic and thermodynamic effects of mobile protons on the promotion of unique chemical processes of ice.

Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy

The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H⁺(aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H⁺⋅(H₂O)_n cluster ions isolated in the gas phase shed light on the chemical nature of H⁺(aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared-infrared (IR-IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H⁺(aq).

Recent Progress in the Manipulation of Molecules with DC Electric Fields

The structure and reactivity of a molecule in the condensed phase are governed by its intermolecular interactions with the surrounding environment. The multipole expansion of each molecule in the condensed phase indicates that the intermolecular interactions are essentially electrostatic (e.g., ion-dipole, dipole-dipole, dipole-quadrupole, dipole-induced dipole). The electrostatic field is a fundamental language of intermolecular communications. Therefore, understanding the influence of the electrostatic field on a molecule, that is, the mechanisms by which an electrostatic field manipulates a molecule, from the perspective of molecular structure, energy states, and dynamics is indispensable for illustrating and, by extension, controlling the chemistry in molecular systems.In this Account, we describe the recent progress made in manipulation of molecular processes using an external DC electrostatic field. An electrostatic field with unprecedentedly high strength (≤4 × 10⁸ V/m) was applied in a controlled manner across a molecular film sample using the ice film nanocapacitor method. This field strength is comparable in magnitude to that of weak intermolecular interactions such as van der Waals interactions in the condensed phases. The samples were prepared using a thin film growing technique in vacuum to obtain the desired chemically tailored molecular systems. The examples of prepared systems included small molecules and molecular clusters isolated in cryogenic Ar matrices, frozen molecular films in amorphous or crystalline phase, and interfaces of multilayered molecular films. The response of the molecules to the external field was monitored by reflection-absorption infrared spectroscopy. This approach allowed us to investigate a variety of molecular systems with various intermolecular strength and environments under the influence of strong electrostatic fields. The range of observed molecular behaviors includes the manipulation of molecular orientation, intramolecular dynamics, and proton transfer reactions as an example of stereodynamic control of chemical reactivity. These observations improve our understanding of molecular behaviors in strong electric fields and broaden our perspective on electrostatic manipulation of molecules. This information is also relevant to a variety of research topics in physical and biological sciences where electric fields play a role in molecular and biological functions.

Green Synthesis, SC-XRD, Non-Covalent Interactive Potential and Electronic Communication via DFT Exploration of Pyridine-Based Hydrazone

Ultrasound-based synthesis at room temperature produces valuable compounds greener and safer than most other methods. This study presents the sonochemical fabrication and characterization of a pyridine-based halogenated hydrazone, (E)-2-((6-chloropyridin-2-yl)oxy)-N′-(2-hydroxybenzylidene) acetohydrazide (HBPAH). The NMR spectroscopic technique was used to determine the structure, while SC-XRD confirmed its crystalline nature. Our structural studies revealed that strong, inter-molecular attractive forces stabilize this crystalline organic compound. Moreover, the compound was optimized at the B3LYP/6-311G(d,p) level using the Crystallographic Information File (CIF). Natural bonding orbital (NBO) and natural population analysis (NPA) were performed at the same level using optimized geometry. Time-dependent density functional theory (DFT) was performed at the B3LYP/6-311G (d,p) method to calculate the frontier molecular orbitals (FMOs) and molecular electrostatic potential (MEP). The global reactivity descriptors were determined using HOMO and LUMO energy gaps. Theoretical calculations based on the Quantum Theory of Atoms in Molecules (QT-AIM) and Hirshfeld analyses identified the non-covalent and covalent interactions of the HBPAH compound. Consequently, QT-AIM and Hirshfeld analyses agree with experimental results.

Féry Infrared Spectrometer for Single-Shot Analysis of Protein Dynamics

Current submillisecond time-resolved broad-band infrared spectroscopy, one of the most frequently used techniques for studying structure-function relationships in life sciences, is typically limited to fast-cycling reactions that can be repeated thousands of times with high frequency. Notably, a majority of chemical and biological processes do not comply with this requirement. For example, the activation of vertebrate rhodopsin, a prototype of many protein receptors in biological organisms that mediate basic functions of life, including vision, smell, and taste, is irreversible. Here we present a dispersive single-shot Féry spectrometer setup that extends such spectroscopy to irreversible and slow-cycling systems by exploiting the unique properties of brilliant synchrotron infrared light combined with an advanced focal plane detector array embedded in a dispersive optical concept. We demonstrate our single-shot method on microbial actinorhodopsin with a slow photocycle and on vertebrate rhodopsin with irreversible activation.

Exploring Nature and Predicting Strength of Hydrogen Bonds: A Correlation Analysis Between Atoms-in-Molecules Descriptors, Binding Energies, and Energy Components of Symmetry-Adapted Perturbation Theory

This work studies the underlying nature of H-bonds (HBs) of different types and strengths and tries to predict binding energies (BEs) based on the properties derived from wave function analysis. A total of 42 HB complexes constructed from 28 neutral and 14 charged monomers were considered. This set was designed to sample a wide range of HB strengths to obtain a complete view about HBs. BEs were derived with the accurate coupled cluster singles and doubles with perturbative triples correction (CCSD(T))(T) method and the physical components of the BE were investigated by symmetry-adapted perturbation theory (SAPT). Quantum theory of atoms-in-molecules (QTAIM) descriptors and other HB indices were calculated based on high-quality density functional theory wave functions. We propose a new and rigorous classification of H-bonds (HBs) based on the SAPT decomposition. Neutral complexes are either classified as "very weak" HBs with a BE ≥ -2.5 kcal/mol that are mainly dominated by both dispersion and electrostatic interactions or as "weak-to-medium" HBs with a BE varying between -2.5 and -14.0 kcal/mol that are only dominated by electrostatic interactions. On the other hand, charged complexes are divided into "medium" HBs with a BE in the range of -11.0 to -15.0 kcal/mol, which are mainly dominated by electrostatic interactions, or into "strong" HBs whose BE is more negative than -15.0 kcal/mol, which are mainly dominated by electrostatic together with induction interactions. Among various explored correlations between BEs and wave function-based HB descriptors, a fairly satisfactory correlation was found for the electron density at the bond critical point (BCP; ρ_BCP ) of HBs. The fitted equation for neutral complexes is BE/kcal/mol =  - 223.08 × ρ_BCP /a. u. + 0.7423 with a mean absolute percentage error (MAPE) of 14.7%, while that for charged complexes is BE/kcal/mol =  - 332.34 × ρ_BCP /a. u. - 1.0661 with a MAPE of 10.0%. In practice, these equations may be used for a quick estimation of HB BEs, for example, for intramolecular HBs or large HB networks in biomolecules. © 2019 Wiley Periodicals, Inc.

Decomposition of the Experimental Raman and Infrared Spectra of Acidic Water into Proton, Special Pair, and Counterion Contributions

Textbooks describe excess protons in liquid water as hydronium (H₃O⁺) ions, although their true structure remains lively debated. To address this question, we have combined Raman and infrared (IR) multivariate curve resolution spectroscopy with ab initio molecular dynamics and anharmonic vibrational spectroscopic calculations. Our results are used to resolve, for the first time, the vibrational spectra of hydrated protons and counterions and reveal that there is little ion-pairing below 2 M. Moreover, we find that isolated excess protons are strongly IR active and nearly Raman inactive (with vibrational frequencies of ∼1500 ± 500 cm^-1), while flanking water OH vibrations are both IR and Raman active (with higher frequencies of ∼2500 ± 500 cm^-1). The emerging picture is consistent with Georg Zundel's seminal work, as well as recent ultrafast dynamics studies, leading to the conclusion that protons in liquid water are primarily hydrated by two flanking water molecules, with a broad range of proton hydrogen bond lengths and asymmetries.

Molecular modeling and assignment of IR spectra of the hydrated excess proton in isotopically dilute water

Infrared (IR) spectroscopy of the water O-H stretch has been widely used to probe both the local hydrogen-bonding structure and dynamics of aqueous systems. Although of significant interest, the IR spectroscopy of excess protons in water remains difficult to assign as a result of extensive and strong intermolecular interactions in hydrated proton complexes. As an alternate approach, we develop a mixed quantum-classical model for the vibrational spectroscopy of the excess proton in isotopically dilute water that draws on frozen proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multi-state empirical valence bond proton model (MS-EVB 3.2). A semi-empirical single oscillator spectroscopic map for the instantaneous transition frequency and transition dipole moment is constructed using potential energy surfaces for the O-H stretch coordinate of the excess proton using electronic structure calculations. Calculated spectra are compared with experimental spectra of dilute H⁺ in D₂O obtained from double-difference FTIR to demonstrate the validity of the map. The model is also used to decompose IR spectra into contributions from different aqueous proton configurations. We find that the O-H transition frequency continuously decreases as the oxygen-oxygen length for a special pair proton decreases, shifting from Eigen- to Zundel-like configurations. The same shift is accompanied by a shift of the flanking water stretches of the Zundel complex to higher frequency than the hydronium O-H vibrations.

Role of Presolvation and Anharmonicity in Aqueous Phase Hydrated Proton Solvation and Transport

Results from condensed phase ab initio molecular dynamics (AIMD) simulations suggest a proton transfer reaction is facilitated by "presolvation" in which the hydronium is transiently solvated by four water molecules, similar to the typical solvation structure of water, by accepting a weak hydrogen bond from the fourth water molecule. A new version 3.2 multistate empirical valence bond (MS-EVB 3.2) model for the hydrated excess proton incorporating this presolvation behavior is therefore developed. The classical MS-EVB simulations show similar structural properties as those of the previous model but with significantly improved diffusive behavior. The inclusion of nuclear quantum effects in the MS-EVB also provides an even better description of the proton diffusion rate. To quantify the influence of anharmonicity, a second model (aMS-EVB 3.2) is developed using the anharmonic aSPC/Fw water model, which provides similar structural properties but improved spectroscopic responses at high frequencies.

Molecular Modes and Dynamics of HCl and DCl Guests of Gas Clathrate Hydrates

Recent years have yielded advances in the placement of unusual molecules as guests within clathrate hydrates (CHs) without severe distortion of the classic lattice structures. Reports describing systems for which observable but limited distortion does occur are available for methanol, ammonia, acetone, and small ether molecules. In these particular examples, the large-cage molecules often participate as non-classical guests H-bonded to the cage walls. Here, we expand the list of such components to include HCl/DCl and HBr as small-cage guests. Based on FTIR spectra of nanocrystalline CHs from two distinct preparative methods combined with critical insights derived from on-the-fly molecular dynamics and ab initio computational data, a coherent argument emerges that these strong acids serve as a source of molecular small-cage guests, ions, and orientational defects. Depending on the HCl/DCl content the ions, defects and molecular guests determine the CH structures, some of which form in sub-seconds via an all-vapor preparative method.

Ions at hydrophobic interfaces

We review the present understanding of the behavior of ions at the air-water and oil-water interfaces. We argue that while the alkali metal cations remain strongly hydrated and are repelled from the hydrophobic surfaces, the anions must be classified into kosmotropes and chaotropes. The kosmotropes remain strongly hydrated in the vicinity of a hydrophobic surface, while the chaotropes lose their hydration shell and can become adsorbed to the interface. The mechanism of adsorption is still a subject of debate. Here, we argue that there are two driving forces for anionic adsorption: the hydrophobic cavitational energy and the interfacial electrostatic surface potential of water. While the cavitational contribution to ionic adsorption is now well accepted, the role of the electrostatic surface potential is much less clear. The difficulty is that even the sign of this potential is a subject of debate, with the ab initio and the classical force field simulations predicting electrostatic surface potentials of opposite sign. In this paper, we will argue that the strong anionic adsorption found in the polarizable force field simulations is the result of the artificial electrostatic surface potential present in the classical water models. We will show that if the adsorption of anions were as large as predicted by the polarizable force field simulations, the excess surface tension of the NaI solution would be strongly negative, contrary to the experimental measurements. While the large polarizability of heavy halides is a fundamental property and must be included in realistic modeling of the electrolyte solutions, we argue that the point charge water models, studied so far, are incompatible with the polarizable ionic force fields when the translational symmetry is broken. The goal for the future should be the development of water models with very low electrostatic surface potential. We believe that such water models will be compatible with the polarizable force fields, which can then be used to study the interaction of ions with hydrophobic surfaces and proteins.

NMR and FT-IR Studies on the Association of Derivatives of Thymidine, Adenosine, and 6-N-Methyl-Adenosine in Aprotic Solvents

Associates of 3',5'-O-TBDMS protected derivatives of the 2'-deoxy forms of the nucleosides adenosine, 6-N-methyl-adenosine and thymidine (henceforward simply addressed by their parents' names) and further model systems in dichloromethane and Freon (CDClF2/CDF3) solutions are studied at low temperatures by 1H NMR and FT-IR spectroscopy. N⋯N distances in hydrogen bonds are estimated from chemical shifts of protons in hydrogen bonds employing geometric and spectroscopic hydrogen bond correlations. These distances are in turn employed to derive N–H stretching frequencies from IR spectroscopic hydrogen bond correlations which may be compared to corresponding experimental results. Three isomeric hydrogen bonded dimers of thymidine are characterized in Freon solution at 120 K. Binary associates of thymidine and a series of pyridines are studied; estimated N⋯N distances in the range of 3.08 to 2.85 Å are qualitatively correlated to shifts of N–H stretching bands where in all cases considerable contributions are found in the spectral region below 3000 cm-1. For adenosine, three isomeric binary associates with 4-nitrophenol are found allowing for an assessment of site-specific acceptor capabilities. In associates of thymidine and adenosine, Watson-Crick and Hoogsteen type 1:1 associates (estimated N⋯N distances of 2.85 and 2.90 Å) as well as 2:1 associates bearing only marginally longer H-bonds could be characterized. Two 1:1 associates between thymidine and 6-N-methyl-adenosine are described that are exclusively bonded via N–H⋯N bridges of about 2.97 and 3.08 Å for Watson-Crick and Hoogsteen sites, respectively, which leads to the conclusion that cooperative effects among coupled N–H⋯O and N–H⋯N hydrogen bonds in A-T base pairs are significant as formation of the N–H⋯O bond induces a contraction of around 0.15 Å in the neighboring N–H⋯N bond.

Probing Small Protonated Water Clusters in Acetonitrile Solutions by 1H NMR

In a previous publication by Kalish et al. (J. Phys. Chem. A 115 (2011) 4063) the existence of well defined small protonated water clusters in acetonitrile has been established by IR spectroscopy. Here we report on a 1H NMR study of triflic acid, CF3SO3H, in acetonitrile-water solutions. Using NMR we are able to corroborate the general solvation scheme we have proposed for the hydrated proton in acetonitrile as a function of the molar ratio between the strong mineral acid and water, n = [H2O]/[acid]. According to this scheme, backed now by both IR absorption spectroscopy and NMR measurements, the very strong triflic acid completely dissociates in acetonitrile/water solutions to yield the aqueous proton and the triflate anion when n > 1. Furthermore, increasing n results in the proton solvated in increasingly larger water clusters formed within the acetonitrile solution.

Clearly distinguishable by NMR are the smallest protonated water clusters, the protonated water monomer, H3 +O, and the protonated water dimer, H5 +O2, which dominate the solution for n = 1,2,3. For larger n the NMR study indicates the gradual increase of the average protonated water cluster size as a function of n while the proton inner solvation core more closely retaining the characteristics of a deformed protonated water dimer, (H2O-H+⋯OH2) s than that of the protonated water monomer (H3 +O) s .

Hydrated alizarin complexes: hydrogen bonding and proton transfer

We investigated the hydrogen bonding structures and proton transfer for the hydration complexes of alizarin (Az) produced in a supersonic jet using fluorescence excitation (FE), dispersed laser induced fluorescence (LIF), visible-visible hole burning (HB), and fluorescence detected infrared (FDIR) spectroscopy. The FDIR spectrum of bare Az with two O-H groups exhibits two vibrational bands at 3092 and 3579 cm(-1), which, respectively, correspond to the stretching vibration of O1-H1 that forms a strong intramolecular hydrogen bond with the C9=O9 carbonyl group and the stretching vibration of O2-H2 that is weakly hydrogen-bonded to O1-H1. For the 1:1 hydration complex Az(H(2)O)(1), we identified three conformers. In the most stable conformer, the water molecule forms hydrogen bonds with the O1-H1 and O2-H2 groups of Az as a proton donor and proton acceptor, respectively. In the other conformers, the water binds to the C10=O10 group in two nearly isoenergetic configurations. In contrast to the sharp vibronic peaks in the FE spectra of Az and Az(H(2)O)(1), only broad, structureless absorption was observed for Az(H(2)O)(n) (n≥ 2), indicating a facile decay process, possibly due to proton transfer in the electronic excited state. The FDIR spectrum with the wavelength of the probe laser fixed at the broad band exhibited a broad vibrational band near the O2-H2 stretching vibration frequency of the most stable conformer of Az(H(2)O)(1). With the help of theoretical calculations, we suggest that the broad vibrational band may represent the occurrence of proton transfer by tunnelling in the electronic ground state of Az(H(2)O)(n) (n≥ 2) upon excitation of the O2-H2 vibration.

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.

Probing the dynamics of solvation and structure of the OH- ion in aqueous solution from picosecond transient absorption measurements

The reaction of intracomplex proton transfer (44BPY^-....HO-H) → 44BPYH^. + OH^- that follows the photoreduction of 4,4’-bipyridine (44BPY) into its anion radical 44BPY^- in the presence of 1,4-diazabicyclo[2.2.2]octane (DABCO) is investigated in acetonitrile-water mixtures by using picosecond transient absorption. The dependence of the appearance kinetics of the 44BPYH^. radical on the water content reveals a highly diffusional proton transfer process that is controlled by the dynamics of solvation of the released hydroxide ion. The results are interpreted on the basis of a two-step mechanism where an intermediate solvation complex (44BPYH^.)OH^-(H₂O)₃ is formed first before evolving toward a final four-water hydration structure OH^-(H₂O)₄.

Direct ab initio dynamics calculations for rates and the kinetic isotope effects of multiproton transfer in ClONO2 + HCl --> HNO3 + Cl2 reactions with water clusters: breakdown of the rule of the geometric mean

We performed high-level quantum mechanical calculations and direct ab initio reaction dynamics calculations for multiple proton transfers in ClONO(2)+HCl-->HNO(3)+Cl(2) with water clusters containing one to two water molecules, which can be used as a model of the reactions occurring on ice surface in stratospheric clouds. The energy barriers of these reactions depend on the number of water molecules involved. Two and three protons in these reactions with one and two water molecules, respectively, were transferred concertedly and asynchronously. The potential energy barrier at the MP2/6-311++(3df,3pd)//MP2/6-31G(d,p) level was 4.8 kcal/mol for the triple proton transfer involving two water molecules with a rate constant of 1.6x10(3) s(-1) at 197 K. The potential energy curve near the saddle points was very flat and the tunneling effect on the proton transfer was negligible. The primary HH/DH kinetic isotope effect for the double proton transfer involving one water molecule was lower than unity due to the enhanced force constant at the transition state. The rule of the geometric mean for the concerted proton transfer does not hold in these reactions because the zero-point energy changes of each proton in flight at the transition state are not the same in the highly asynchronous processes.