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

Interfacial chemical reactivity enhancement

Interfacial enhancements of chemical reaction equilibria and rates in liquid droplets are predicted using a combined theoretical and experimental analysis strategy. Self-consistent solutions of reaction and adsorption equilibria indicate that interfacial reactivity enhancement is driven primarily by the adsorption free energy of the product (or activated complex). Reactant surface activity has a smaller indirect influence on reactivity due to compensating reactant interfacial concentration and adsorption free energy changes, as well as adsorption-induced depletion of the droplet core. Experimental air-water interfacial adsorption free energies and critical micelle concentration correlations provide quantitative surface activity estimates as a function of molecular structure, predicting an increase in interfacial reactivity with increasing product size and decreasing product polarity, aromaticity, and charge (but less so for anions than cations). Reactions with small, neutral, or charged products are predicted to have little reactivity enhancement at an air-water interface unless the product is rendered sufficiently surface active by, for example, interactions with interfacial water dangling OH groups, charge transfer, or voltage fluctuations.

Open-[60]fullerenols with water adsorbed both inside and outside

The water affinity on [60]fullerenols was found to be governed by surface electrostatic potential while water aggregation is initiated by the hydroxy groups attached on the carbon surface. The molecular water adsorption at the internal sphere caused a significnat inhibition of water adsorption at the external carbon surface.

SHINERS Study of Chloride Order-Disorder Phase Transition and Solvation of Cu(100)

Shell-isolated nanoparticle enhanced Raman spectroscopy (SHINERS) and density functional theory (DFT) are used to probe Cl- adsorption and the order-disorder phase transition associated with the c(2 × 2) Cl- adlayer on Cu(100) in acid media. A two-component ν(Cu-Cl) vibrational band centered near 260 ± 1 cm-1 is used to track the potential dependence of Cl- adsorption. The potential dependence of the dominant 260 cm-1 component tracks the coverage of the fluctional c(2 × 2) Cl- phase on terraces in good agreement with the normalized intensity of the c(2 × 2) superstructure rods in prior surface X-ray diffraction (SXRD) studies. As the c(2 × 2) Cl- coverage approaches saturation, a second ν(Cu-Cl) component mode emerges between 290 and 300 cm-1 that coincides with the onset and stiffening of step faceting where Cl- occupies the threefold hollow sites to stabilize the metal kink saturated Cu <100> step edge. The formation of the c(2 × 2) Cl- adlayer is accompanied by the strengthening of ν(O-H) stretching modes in the adjacent non-hydrogen-bonded water at 3600 cm-1 and an increase in hydronium concentration evident in the flanking H2O modes at 3100 cm-1. The polarization of the water molecules and enrichment of hydronium arise from the combination of Cl- anionic character and lateral templating provided by the c(2 × 2) adlayer, consistent with SXRD studies. At negative potentials, Cl- desorption occurs followed by development of a sulfate νs(S═O) band. Below -1.1 V vs Hg/HgSO4, a new 200 cm-1 mode emerges congruent with hydride formation and surface reconstruction reported in electrochemical scanning tunneling microscopy studies.

Thermal dependence of the hydrated proton and optimal proton transfer in the protonated water hexamer

Water is a key ingredient for life and plays a central role as solvent in many biochemical reactions. However, the intrinsically quantum nature of the hydrogen nucleus, revealing itself in a large variety of physical manifestations, including proton transfer, gives rise to unexpected phenomena whose description is still elusive. Here we study, by a combination of state-of-the-art quantum Monte Carlo methods and path-integral molecular dynamics, the structure and hydrogen-bond dynamics of the protonated water hexamer, the fundamental unit for the hydrated proton. We report a remarkably low thermal expansion of the hydrogen bond from zero temperature up to 300 K, owing to the presence of short-Zundel configurations, characterised by proton delocalisation and favoured by the synergy of nuclear quantum effects and thermal activation. The hydrogen bond strength progressively weakens above 300 K, when localised Eigen-like configurations become relevant. Our analysis, supported by the instanton statistics of shuttling protons, reveals that the near-room-temperature range from 250 K to 300 K is optimal for proton transfer in the protonated water hexamer.

Calculating Vibrational Excited State Absorptions with Excited State Constrained Minimized Energy Surfaces

The modeling and interpretation of vibrational spectra are crucial for studying reaction dynamics using vibrational spectroscopy. Most prior theoretical developments focused on describing fundamental vibrational transitions while fewer developments focused on vibrational excited state absorptions. In this study, we present a new method that uses excited state constrained minimized energy surfaces (CMESs) to describe vibrational excited state absorptions. The excited state CMESs are obtained similarly to the previous ground state CMES development in our group but with additional wave function orthogonality constraints. Using a series of model systems, including the harmonic oscillator, Morse potential, double-well potential, quartic potential, and two-dimensional anharmonic potential, we demonstrate that this new procedure provides good estimations of the transition frequencies for vibrational excited state absorptions. These results are significantly better than those obtained from harmonic approximations using conventional potential energy surfaces, demonstrating the promise of excited state CMES-based methods for calculating vibrational excited state absorptions in real systems.

Ion-Size Dependent Adsorption Crossover on the Surface of a Water Droplet

The adsorption of ionic and neutral spherical solutes on the surface of a liquid water droplet are investigated using molecular dynamics simulations and theoretical analyses. The results reveal a crossover in the sign of the adsorption free energy as a function of ion size, with ions larger than iodide predicted to be increasingly surface active. Adsorption free energies are decomposed into competing energetic and entropic contributions arising from direct solute-water interaction energy and its fluctuations. The entropically driven surface activity of large ions is predicted to increase with ion size, while small ions are typically driven away from the interface by a more delicate balance of energetic and entropic contributions, with a nonmonotonic ion size dependence linked to the ion's hydration-shell structure and stability. The physical interpretation of the results is illuminated by comparisons with dielectric linear response and cavity formation predictions and implications to interfacial acidity and enhanced chemical reactivity are discussed.

Proton-transfer spectroscopy beyond the normal-mode scenario

A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single [Formula: see text] cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.

Hydrogen Bond Thermodynamics in Aqueous Acid Solutions: A Combined DFT and Classical Force-Field Approach

The thermodynamics of hydrogen bonds in aqueous and acidic solutions significantly impacts the kinetics and thermodynamics of acid reaction chemistry. We utilize in this work a multiscale approach, combining density functional theory (DFT) with classical molecular dynamics (MD) to model hydrogen bond thermodynamics in an acidic solution. Using thermodynamic cycles, we split the solution phase free energy into its gas phase counterpart plus solvation free energies. We validate this DFT/MD approach by calculating the aqueous phase hydrogen bond free energy between two water molecules (H₂O-···-H₂O), the free energy to transform an H₃O⁺ cation into an H₅O₂⁺ cation, and the hydrogen bond free energy of protonated water clusters (H₃O⁺-···-H₂O and H₅O₂⁺-···-H₂O). The computed equilibrium hydrogen bond free energy of H₂O-···-H₂O is remarkably accurate, especially considering the large individual contributions to the thermodynamic cycle. Turning to cations, we find the ion to be more stable than H₃O⁺ by roughly 1-2 k_BT. This small free energy difference allows for thermal fluctuation between the two idealized motifs, consistent with spectroscopic and simulation studies. Lastly, hydrogen bonding free energies between either H⁺ cation and H₂O in solution were found to be stronger than between two H₂O, though much less so than in vacuum because of dielectric screening in solution. Altogether, our results suggest the DFT/MD approach is promising for application in modeling hydrogen bonding and proton transfer thermodynamics in condensed phases.

Importance of Hydrogen Bonding in Crowded Environments: A Physical Chemistry Perspective

Cells are heterogeneous on every length and time scale; cytosol contains thousands of proteins, lipids, nucleic acids, and small molecules, and molecular interactions within this crowded environment determine the structure, dynamics, and stability of biomolecules. For decades, the effects of crowding at the atomistic scale have been overlooked in favor of more tractable models largely based on thermodynamics. Crowding can affect the conformations and stability of biomolecules by modulating water structure and dynamics within the cell, and these effects are nonlocal and environment dependent. Thus, characterizing water's hydrogen-bond (H-bond) networks is a critical step toward a complete microscopic crowding model. This perspective provides an overview of molecular crowding and describes recent time-resolved spectroscopy approaches investigating H-bond networks and dynamics in crowded or otherwise complex aqueous environments. Ultrafast spectroscopy combined with atomistic simulations has emerged as a powerful combination for studying H-bond structure and dynamics in heterogeneous multicomponent systems. We discuss the ongoing challenges toward developing a complete atomistic description of macromolecular crowding from an experimental as well as a theoretical perspective.

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200–300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm⁻¹, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm⁻¹ and agrees in location and width with previous experimental mid-infrared spectra.

The spectroscopic signatures of excess protons in HCl solutions are studied by ab initio simulations and THz experiments. Two contributions beyond the normal-mode scenario are identified that reflect proton-waiting and proton-transfer processes.

Infrared Spectra and Hydrogen-Bond Configurations of Water Molecules at the Interface of Water-Insoluble Polymers under Humidified Conditions

Elucidating the state of interfacial water, especially the hydrogen-bond configurations, is considered to be key for a better understanding of the functions of polymers that are exhibited in the presence of water. Here, an analysis in this direction is conducted for two water-insoluble biocompatible polymers, poly(2-methoxyethyl acrylate) and cyclic(poly(2-methoxyethyl acrylate)), and a non-biocompatible polymer, poly(n-butyl acrylate), by measuring their IR spectra under humidified conditions and by carrying out theoretical calculations on model complex systems. It is found that the OH stretching bands of water are decomposed into four components, and while the higher-frequency components (with peaks at ∼3610 and ∼3540 cm^-1) behave in parallel with the C═O and C-O-C stretching and CH deformation bands of the polymers, the lower-frequency components (with peaks at ∼3430 and ∼3260 cm^-1) become pronounced to a greater extent with increasing humidity. From the theoretical calculations, it is shown that the OH stretching frequency that is distributed from ∼3650 to ∼3200 cm^-1 is correlated to the hydrogen-bond configurations and is mainly controlled by the electric field that is sensed by the vibrating H atom. By combining these observed and calculated results, the configurations of water at the interface of the polymers are discussed.

Microfluidic In-Situ Spectrophotometric Approaches to Tackle Actinides Analysis in Multiple Oxidation States

The study and development of present and future processes for the treatment/recycling of spent nuclear fuels require many steps, from design in the laboratory to setting up on an industrial scale. In all of these steps, analysis and instrumentation are key points. For scientific reasons (small-scale studies, control of phenomena, etc.) but also with regard to minimizing costs, risks, and waste, such developments are increasingly carried out on milli- or microfluidic devices. The logic is the same for the chemical analyses associated with their follow-up and interpretation. Due to this, over the last few years, opto-microfluidic analysis devices adapted to the monitoring of different processes (dissolution, liquid-liquid extraction, precipitation, etc.) have been increasingly designed and developed. In this work, we prove that photonic lab-on-a-chip (PhLoC) technology is fully suitable for all actinides concentration monitoring along the plutonium uranium refining extraction (plutonium, uranium, reduction, extraction, or Purex) process. Several PhLoC microfluidic platforms were specifically designed and used in different nuclear research and development (R&D) laboratories, to tackle actinides analysis in multiple oxidation states even in mixtures. The detection limits reached (tens of µmol·L^-1) are fully compliant with on-line process monitoring, whereas a range of analyzable concentrations of three orders of magnitude can be covered with less than 150 µL of analyte. Finally, this work confirms the possibility and the potential of coupling Raman and ultraviolet-visible (UV-Vis) spectroscopies at the microfluidic scale, opening the perspective of measuring very complex mixtures.

Hydration Dynamics of Model Peptides with Different Hydrophobic Character

The multi-scale dynamics of aqueous solutions of the hydrophilic peptide N-acetyl-glycine-methylamide (NAGMA) have been investigated through extended frequency-range depolarized light scattering (EDLS), which enables the broad-band detection of collective polarizability anisotropy fluctuations. The results have been compared to those obtained for N-acetyl-leucinemethylamide (NALMA), an amphiphilic peptide which shares with NAGMA the same polar backbone, but also contains an apolar group. Our study indicates that the two model peptides induce similar effects on the fast translational dynamics of surrounding water. Both systems slow down the mobility of solvating water molecules by a factor 6–8, with respect to the bulk. Moreover, the two peptides cause a comparable far-reaching spatial perturbation extending to more than two hydration layers in diluted conditions. The observed concentration dependence of the hydration number is explained considering the random superposition of different hydration shells, while no indication of solute aggregation phenomena has been found. The results indicate that the effect on the dynamics of water solvating the amphiphilic peptide is dominated by the hydrophilic backbone. The minor impact of the hydrophobic moiety on hydration features is consistent with structural findings derived by Fourier transform infrared (FTIR) measurements, performed in attenuated total reflectance (ATR) configuration. Additionally, we give evidence that, for both systems, the relaxation mode in the GHz frequency range probed by EDLS is related to solute rotational dynamics. The rotation of NALMA occurs at higher timescales, with respect to the rotation of NAGMA; both processes are significantly slower than the structural dynamics of hydration water, suggesting that solute and solvent motions are uncoupled. Finally, our results do not indicate the presence of super-slow water (relaxation times in the order of tens of picoseconds) around the peptides investigated.

Quantum mechanical effects in acid–base chemistry

Acid–base chemistry has immense importance for explaining and predicting the chemical products formed by an acid and a base when mixed together. However, the traditional chemistry theories used to describe acid–base reactions do not take into account the effect arising from the quantum mechanical nature of the acidic hydrogen shuttling potential and its dependence on the acid base distance. Here, infrared and NMR spectroscopies, in combination with first principles simulations, are performed to demonstrate that quantum mechanical effects, including electronic and nuclear quantum effects, play an essential role in defining the acid–base chemistry when 1-methylimidazole and acetic acid are mixed together. In particular, it is observed that the acid and the base interact to form a complex containing a strong hydrogen bond, in which the acidic hydrogen atom is neither close to the acid nor to the base, but delocalized between them. In addition, the delocalization of the acidic hydrogen atom in the complex leads to characteristic IR and NMR signatures. The presence of a hydrogen delocalized state in this simple system challenges the conventional knowledge of acid–base chemistry and opens up new avenues for designing materials in which specific properties produced by the hydrogen delocalized state can be harvested.

Acid-based theories do not consider the quantum mechanical nature of the acidic hydrogen shuttling potential. Here, it is demonstrated that this particularity is needed to explain the formation acid-base complex with a delocalized acidic hydrogen.

The Proton in Biochemistry: Impacts on Bioenergetics, Biophysical Chemistry, and Bioorganic Chemistry

The proton is the smallest atomic particle, and in aqueous solution it is the smallest hydrated ion, having only two waters in its first hydration shell. In this article we survey key aspects of the proton in chemistry and biochemistry, starting with the definitions of pH and pK _a and their application inside biological cells. This includes an exploration of pH in nanoscale spaces, distinguishing between bulk and interfacial phases. We survey the Eigen and Zundel models of the structure of the hydrated proton, and how these can be used to explain: a) the behavior of protons at the water-hydrophobic interface, and b) the extraordinarily high mobility of protons in bulk water via Grotthuss hopping, and inside proteins via proton wires. Lastly, we survey key aspects of the effect of proton concentration and proton transfer on biochemical reactions including ligand binding and enzyme catalysis, as well as pH effects on biochemical thermodynamics, including the Chemiosmotic Theory. We find, for example, that the spontaneity of ATP hydrolysis at pH ≥ 7 is not due to any inherent property of ATP (or ADP or phosphate), but rather to the low concentration of H⁺. Additionally, we show that acidification due to fermentation does not derive from the organic acid waste products, but rather from the proton produced by ATP hydrolysis.

Water as a probe for pH measurement in individual particles using micro-Raman spectroscopy

Atmospheric aerosol acidity impacts numerous physicochemical processes, but the determination of particle pH remains a significant challenge due to the nonconservative nature of the H⁺ concentration ([H⁺]). Traditional measurements have difficulty in describing the practical state of an aerosol because they comprise chemical components or hypotheses that change the nature of the particles. In this work, we present a direct pH measurement that uses water as a general probe to detect [H⁺] in individual particles by micro-Raman spectroscopy. Containing the vibrational bands of ions and water influenced by ions, the spectra of hydrated ion were decomposed from the solution spectra as standard spectra by multivariate curve resolution analysis. Meanwhile, ratios of hydrated ions were calculated between the Raman spectra and standard spectra to evaluate concentration profiles of each ion. It demonstrated that good quantitative models between the ratio and concentration for all ions including H⁺ can be built with correlation coefficients (R²) higher than 0.95 for the solutions. The method was further applied to individual particle pH measurement. The pH value of sulfate aerosol particles was calculated, and the standard error was 0.09 using pH values calculated from the [HSO₄^-]/[SO₄^2-] as a reference. Furthermore, the applicability of the method was proven by detecting the pH value of chloride particles. Therefore, utilizing water, the most common substance, as the spectroscopic probe to measure [H⁺] without restriction of the ion system, this method has potential to measure the pH value of atmospheric particles with various compounds, although more work needs to be done to improve the sensitivity of the method.

Resolving the Structural Debate for the Hydrated Excess Proton in Water

It has long been proposed that the hydrated excess proton in water (aka the solvated "hydronium" cation) likely has two limiting forms, that of the Eigen cation (H₉O₄⁺) and that of the Zundel cation (H₅O₂⁺). There has been debate over which of these two is the more dominant species and/or whether intermediate (or "distorted") structures between these two limits are the more realistic representation. Spectroscopy experiments have recently provided further results regarding the excess proton. These experiments show that the hydrated proton has an anisotropy reorientation time scale on the order of 1-2 ps. This time scale has been suggested to possibly contradict the picture of the more rapid "special pair dance" phenomenon for the hydrated excess proton, which is a signature of a distorted Eigen cation. The special pair dance was predicted from prior computational studies in which the hydrated central core hydronium structure continually switches (O-H···O)* special pair hydrogen-bond partners with the closest three water molecules, yielding on average a distorted Eigen cation with three equivalent and dynamically exchanging distortions. Through state-of-art simulations it is shown here that anisotropy reorientation time scales of the same magnitude are obtained that also include structural reorientations associated with the special pair dance, leading to a reinterpretation of the experimental results. These results and additional analyses point to a distorted and dynamic Eigen cation as the most prevalent hydrated proton species in aqueous acid solutions of dilute to moderate concentration, as opposed to a stabilized or a distorted (but not "dancing") Zundel cation.

Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters

Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H⁺(H₂O)₄ and H⁺(H₂O)₂₁ agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H⁺(H₂O)₅, which can adopt several low-energy conformations.

Protein Hydration in a Bioprotecting Mixture

We combined broad-band depolarized light scattering and infrared spectroscopies to study the properties of hydration water in a lysozyme-trehalose aqueous solution, where trehalose is present above the concentration threshold (30% in weight) relevant for biopreservation. The joint use of the two different techniques, which were sensitive to inter-and intra-molecular degrees of freedom, shed new light on the molecular mechanism underlying the interaction between the three species in the mixture. Thanks to the comparison with the binary solution cases, we were able to show that, under the investigated conditions, the protein, through preferential hydration, remains strongly hydrated even in the ternary mixture. This supported the water entrapment scenario, for which a certain amount of water between protein and sugar protects the biomolecule from damage caused by external agents.

Signature of the surface hydrated proton and associated restructuring of water at model membrane interfaces: a vibrational sum frequency generation study

Hydrated proton at membrane interfaces plays an important role in the bioenergetic process of almost all organisms. Herein, the signature of the hydrated proton at membrane interfaces has been investigated by measuring the vibrational sum frequency generated (VSFG) spectra of negatively charged and zwitterionic lipids in the presence of different concentrations of acids. The addition of acids decreases the intensity of the OH stretch of the VSFG signal of water present at the negatively charged and zwitterionic lipids along with the enhanced intensity of the broad VSFG signal in the range of 2500-2800 cm-1. The enhanced intensity of the broad continuum observed in the range of 2500-2800 cm-1 has been assigned to the signature of the hydrated proton at the lipid interfaces. The decrease in the VSFG signal of the OH stretch of water along with the appearance of the broad signal suggests that the hydrated proton exists in the vicinity of the lipid interfaces and restructures the interaction between the interfacial water molecules.

The hopping mechanism of the hydrated excess proton and its contribution to proton diffusion in water

In this work, a series of analyses are performed on ab initio molecular dynamics simulations of a hydrated excess proton in water to quantify the relative occurrence of concerted hopping events and "rattling" events and thus to further elucidate the hopping mechanism of proton transport in water. Contrary to results reported in certain earlier papers, the new analysis finds that concerted hopping events do occur in all simulations but that the majority of events are the product of proton rattling, where the excess proton will rattle between two or more waters. The results are consistent with the proposed "special-pair dance" model of the hydrated excess proton wherein the acceptor water molecule for the proton transfer will quickly change (resonate between three equivalent special pairs) until a decisive proton hop occurs. To remove the misleading effect of simple rattling, a filter was applied to the trajectory such that hopping events that were followed by back hops to the original water are not counted. A steep reduction in the number of multiple hopping events is found when the filter is applied, suggesting that many multiple hopping events that occur in the unfiltered trajectory are largely the product of rattling, contrary to prior suggestions. Comparing the continuous correlation function of the filtered and unfiltered trajectories, we find agreement with experimental values for the proton hopping time and Eigen-Zundel interconversion time, respectively.

Nuclear Quantum Effects Largely Influence Molecular Dissociation and Proton Transfer in Liquid Water under an Electric Field

Proton transfer in liquid water controls acid-base chemistry, crucial enzyme reactions, and the functioning of fuel cells. Externally applied static electric fields in water are capable of dissociating molecules and transferring protons across the H-bond network. However, the impact of nuclear quantum effects (NQEs) on these fundamental field-induced phenomena has not yet been reported. By comparing state-of-the-art ab initio molecular dynamics (AIMD) and path integral AIMD simulations of water under electric fields, I show that quantum delocalization of the proton lowers the molecular ionization threshold to approximately one-third. Moreover, also the water behavior as a protonic semiconductor is considerably modified by the inclusion of NQEs. In fact, when the quantum nature of the nuclei is taken into account, the proton conductivity is ∼50% larger. This work proves that NQEs sizably affect the protolysis phenomenon and proton transfer in room-temperature liquid water.

Molecular Origins of the Barriers to Proton Transport in Acidic Aqueous Solutions

The self-consistent iterative multistate empirical valence bond (SCI-MS-EVB) method is used to analyze the structure, thermodynamics, and dynamics of hydrochloric acid solutions. The reorientation time scales of irreversible proton transport are elucidated by simulating 0.43, 0.85, 1.68, and 3.26 M HCl solutions at 270, 285, 300, 315, and 330 K. The results indicate increased counterion pairing with increasing concentration, which manifests itself via a reduced hydronium oxygen-chloride (O*-Cl) structuring in the radial distribution functions. Increasing ionic concentration also reduces the diffusion of the hydrated excess protons, principally by reducing the contribution of the Grotthuss proton hopping (shuttling) mechanism to the overall diffusion process. In agreement with prior experimental findings, a decrease in the activation energy of reorientation time scales was also observed, which is explicitly explained by using activated rate theory and an energy-entropy decomposition of the state-averaged radial distribution functions. These results provide atomistic verification of suggestions from recent two-dimensional infrared spectroscopy experiments that chloride anions (as opposed to hydrated excess protons) create entropic barriers to proton transport.

Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations

The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H₅O₂ ⁺ and Eigen H₉O₄ ⁺ ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure-spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm^-1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. Finally, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.

Slowing Down of the Molecular Reorientation of Water in Concentrated Alkaline Solutions

It is generally accepted that the hydroxide ion (OH^-) is a strong hydrogen bond acceptor and that its anomalously high diffusion constant in water results from a Grotthuss-like structural diffusion mechanism. However, the spatial extent over which OH^- ions influence the dynamics of the hydrogen-bond network of water remained largely unclear. Here, we measure the ultrafast dynamics of OH groups of HDO molecules interacting with the deuterated hydroxide ion OD^-. For solutions with OD^- concentrations up to 4 M, we find that HDO molecules that are not directly interacting with the ions have a reorientation time constant of ∼2.7 ps, similar to that of pure liquid water. When the concentration of OD^- ions is increased, the reorientation time constant increases, indicating a strong slowing down of the structural dynamics of the solution.

Composition-Dependent Hydrogen-Bonding Motifs and Dynamics in Brønsted Acid-Base Mixtures

In recent years the interaction of organophosphates and imines, which is at the core of Brønsted acid organocatalysis, has been established to be based on strong ionic hydrogen bonds. Yet, besides the formation of homodimers consisting of two acid molecules and heterodimers consisting of one acid and one base, also multimeric molecular aggregates are formed in solution. These multimeric aggregates consist of one base and several acid molecules. The details of the intermolecular bonding in such aggregates, however, have remained elusive. To characterize composition-dependent bonding and bonding dynamics in these aggregates, we use linear and nonlinear infrared (IR) spectroscopy at varying molar ratios of diphenyl phosphoric acid and quinaldine. We identify the individual aggregate species, giving rise to the structured, strong, and very broad infrared absorptions, which span more than 1000 cm^–1. Linear infrared spectra and density functional theory calculations of the proton transfer potential show that doubly ionic intermolecular hydrogen bonds between the acid and the base lead to absorptions which peak at ∼2040 cm^–1. The contribution of singly ionic hydrogen bonds between an acid anion and an acid molecule is observed at higher frequencies. As common to such strong hydrogen bonds, ultrafast IR spectroscopy reveals rapid, ∼ 100 fs, dissipation of energy from the proton transfer coordinate. Yet, the full dissipation of the excess energy occurs on a ∼0.8–1.1 ps time scale, which becomes longer when multimers dominate. Our results thus demonstrate the coupling and collectivity of the hydrogen bonds within these complexes, which enable efficient energy transfer.

Slow Proton Transfer in Nanoconfined Water

The transport of protons in nanoconfined environments, such as in nanochannels of biological or artificial proton conductive membranes, is essential to chemistry, biology, and nanotechnology. In water, proton diffusion occurs by hopping of protons between water molecules. This process involves the rearrangement of many hydrogen bonds and as such can be strongly affected by nanoconfinement. We study the vibrational and structural dynamics of hydrated protons in water nanodroplets stabilized by a cationic surfactant using polarization-resolved femtosecond infrared transient absorption spectroscopy. We determine the time scale of proton hopping in the center of the water nanodroplets from the dynamics of the anisotropy of the transient absorption signals. We find that in small nanodroplets with a diameter <4 nm, proton hopping is more than 10 times slower than in bulk water. Even in relatively large nanodroplets with a diameter of ∼7 nm, we find that the rate of proton hopping is slowed by ∼4 times compared with bulk water.

Understanding and Tracking the Excess Proton in Ab Initio Simulations; Insights from IR Spectra

Proton transport in aqueous media is ubiquitously important in chemical and biological processes. Although ab initio molecular dynamics (AIMD) simulations have made great progress in characterizing proton transport, there has been a long-standing challenge in defining and tracking the excess proton, or more properly, the center of excess charge (CEC) created when a hydrogen nucleus distorts the electron distributions of water molecules in a delocalized and highly dynamic nature. Yet, defining (and biasing) such a CEC is essential when combining AIMD with enhanced sampling methods to calculate the relevant macroscopic properties via free-energy landscapes, which is the standard practice for most processes of interest. Several CEC formulas have been proposed and used, but none have yet been systematically tested or rigorously derived. In this paper, we show that the CEC can be used as a computational tool to disentangle IR features of the solvated excess proton from its surrounding solvent, and in turn, how correlating the features in the excess charge spectrum with the behavior of CEC in simulations enables a systematic evaluation of various CEC definitions. We present a new definition of CEC and show how it overcomes the limitations of those currently available both from a spectroscopic point of view and from a practical perspective of performance in enhanced sampling simulations.

Identifying Eigen-like hydrated protons at negatively charged interfaces

Despite the importance of the hydrogen ion in a wide range of biological, chemical, and physical processes, its molecular structure in solution remains lively debated. Progress has been primarily hampered by the extreme diffuse nature of the vibrational signatures of hydrated protons in bulk solution. Using the inherently surface-specific vibrational sum frequency spectroscopy technique, we show that at selected negatively charged interfaces, a resolved spectral feature directly linked to the H₃O⁺ core in an Eigen-like species can be readily identified in a biologically compatible pH range. Centered at ~2540 cm⁻¹, the band is seen to shift to ~1875 cm⁻¹ when forming D₃O⁺ upon isotopic substitution. The results offer the possibility of tracking and understanding from a molecular perspective the behavior of hydrated protons at charged interfaces.

Hydrated protons are always present in aqueous solution, but their molecular structure remains under debate. Here the authors use vibrational sum frequency spectroscopy to show that at negatively charged liquid–vapor interfaces, protons adopt a specific configuration characteristic of Eigen-like species.

Nature of Excess Hydrated Proton at the Water-Air Interface

Understanding the interfacial molecular structure of acidic aqueous solutions is important in the context of, e.g., atmospheric chemistry, biophysics, and electrochemistry. The hydration of the interfacial proton is necessarily different from that in the bulk, given the lower effective density of water at the interface, but has not yet been elucidated. Here, using surface-specific vibrational spectroscopy, we probe the response of interfacial protons at the water-air interface and reveal the interfacial proton continuum. Combined with spectral calculations based on ab initio molecular dynamics simulations, the proton at the water-air interface is shown to be well-hydrated, despite the limited availability of hydration water, with both Eigen and Zundel structures coexisting at the interface. Notwithstanding the interfacial hydrated proton exhibiting bulk-like structures, a substantial interfacial stabilization by -1.3 ± 0.2 kcal/mol is observed experimentally, in good agreement with our free energy calculations. The surface propensity of the proton can be attributed to the interaction between the hydrated proton and its counterion.

Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk

The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H₂O)₁₀₀H⁺, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.

Hydration of ferric chloride and nitrate in aqueous solutions: water-mediated ion pairing revealed by Raman spectroscopy

Iron is the most abundant transition metal in the earth's crust and is important for the proper functioning of many technological and natural processes. Despite the importance, a complete microscopic understanding of the hydration of ferric ions and water mediated ion pairing has not been realized. Hydrated Fe(iii) is difficult to study due to the process of complexation to the anion and hydrolysis of the hydrating water molecules leading to a heterogeneous solution with diverse speciation. Here, ferric chloride and nitrate aqueous solutions are studied using polarized Raman spectroscopy as a function of concentration and referenced to their respective sodium salt or mineral acid. Perturbed water spectra (PWS) were generated using multivariate curve resolution-alternating least squares (MCR-ALS) to show the residual spectral response uniquely attributable to the hydration of ferric speciation. The hydrogen bonding network associated with the hydrating water molecules in ferric chloride solutions are found to be more similar to hydrochloric acid solutions, whereas in ferric nitrate solutions, the network behaves more similar to sodium nitrate, despite increased acidity. Thus, in the FeNO3 and FeCl3 solutions, ion pairing and coordination, respectively, are significantly influencing the hydration spectra signature. These results further reveal concentration dependent changes to the hydrogen bonding network, hydrating water symmetry, and changes to the relative abundance of solvent shared ion pairs that are governed primarily by the ferric salt identity.

Electrically induced liquid-liquid phase transition in water at room temperature

In this work we expand on findings previously reported [Wexler et al., Phys. Chem. Chem. Phys., 2016, 18, 16281] on the experimental observation of a phase transition in a hydrogen bonded liquid manifesting in long range dipole-dipole interactions. The studied system, liquid water stressed by an electric field, exhibits collective oscillations brought about through spontaneous breakdown of symmetry. Raman spectroscopy identifies the primary excitation of the emergent phase as transverse optically active phonon-like sidebands that appear on the hydrogen bonded asymmetric stretch mode. The phase transition is observed throughout the entire volume of liquid. The system also exhibits a self-similarity relation between the scattered Raman intensity and the electric field strength which further supports the conclusion that collective behavior persists against thermal disruption. The experimental findings are discussed in terms of a quantum field theory for macroscopic quantum systems.

High-Level VSCF/VCI Calculations Decode the Vibrational Spectrum of the Aqueous Proton

The hydrated excess proton is a common species in aqueous chemistry, which complexes with water in a variety of structures. The infrared spectrum of the aqueous proton is particularly sensitive to this array of structures, which manifests as continuous IR absorption from 1000 to 3000 cm^-1 known as the "proton continuum". Because of the extreme breadth of the continuum and strong anharmonicity of the involved vibrational modes, this spectrum has eluded straightforward interpretation and simulation. Using protonated water hexamer clusters from reactive molecular dynamics trajectories, and focusing on their central H⁺(H₂O)₂ structures' spectral contribution, we reproduce the linear IR spectrum of the aqueous proton with a high-level local monomer quantum method and highly accurate many-body potential energy surface. The accuracy of this approach is first verified in the vibrational spectra of the two isomers of the protonated water hexamer in the gas phase. We then apply this approach to 800 H⁺(H₂O)₆ clusters, also written as [H⁺(H₂O)₂](H₂O)₄, drawn from multistate empirical valence bond simulations of the bulk liquid to calculate the infrared spectrum of the aqueous proton complex. Incorporation of anharmonic effects to the vibrational potential and quantum mechanical treatment of the proton produces a better agreement to the infrared spectrum compared to that of the double-harmonic approximation. We assess the correlation of the proton stretching mode with different atomistic coordinates, finding the best correlation with ⟨R_OH⟩, the expectation value of the proton-oxygen distance R_OH. We also decompose the IR spectrum based on normal mode vibrations and ⟨R_OH⟩ to provide insight on how different frequency regions in the continuum report on different configurations, vibrational modes, and mode couplings.

Tracking Aqueous Proton Transfer by Two-Dimensional Infrared Spectroscopy and ab Initio Molecular Dynamics Simulations

Proton transfer in water is ubiquitous and a critical elementary event that, via proton hopping between water molecules, enables protons to diffuse much faster than other ions. The problem of the anomalous nature of proton transport in water was first identified by Grotthuss over 200 years ago. In spite of a vast amount of modern research effort, there are still many unanswered questions about proton transport in water. An experimental determination of the proton hopping time has remained elusive due to its ultrafast nature and the lack of direct experimental observables. Here, we use two-dimensional infrared spectroscopy to extract the chemical exchange rates between hydronium and water in acid solutions using a vibrational probe, methyl thiocyanate. Ab initio molecular dynamics (AIMD) simulations demonstrate that the chemical exchange is dominated by proton hopping. The observed experimental and simulated acid concentration dependence then allow us to extrapolate the measured single step proton hopping time to the dilute limit, which, within error, gives the same value as inferred from measurements of the proton mobility and NMR line width analysis. In addition to obtaining the proton hopping time in the dilute limit from direct measurements and AIMD simulations, the results indicate that proton hopping in dilute acid solutions is induced by the concerted multi-water molecule hydrogen bond rearrangement that occurs in pure water. This proposition on the dynamics that drive proton hopping is confirmed by a combination of experimental results from the literature.

Entropic barriers in the kinetics of aqueous proton transfer

Aqueous proton transport is uniquely rapid among aqueous processes, mediated by fluctuating hydrogen bond reorganization in liquid water. In a process known as Grotthuss diffusion, the excess charge diffuses primarily by sequential proton transfers between water molecules rather than standard Brownian motion, which explains the anomalously high electrical conductivity of acidic solutions. Employing ultrafast IR spectroscopy, we use the orientational anisotropy decay of the bending vibrations of the hydrated proton complex to study the picosecond aqueous proton transfer kinetics as a function of temperature, concentration, and counterion. We find that the orientational anisotropy decay exhibits Arrhenius behavior, with an apparent activation energy of 2.4 kcal/mol in 1M and 2M HCl. Interestingly, acidic solutions at high concentration with longer proton transfer time scales display corresponding decreases in activation energy. We interpret this counterintuitive trend by considering the entropic and enthalpic contributions to the activation free energy for proton transfer. Halide counteranions at high concentrations impose entropic barriers to proton transfer in the form of constraints on the solution's collective H-bond fluctuations and obstruction of potential proton transfer pathways. The corresponding proton transfer barrier decreases due to weaker water-halide H-bonds in close proximity to the excess proton, but the entropic effects dominate and result in a net reduction in the proton transfer rate. We estimate the activation free energy for proton transfer as ∼1.0 kcal/mol at 280 K.

Hydration-Shell Vibrational Spectroscopy

Hydration-shell vibrational spectroscopy provides an experimental window into solute-induced water structure changes that mediate aqueous folding, binding, and self-assembly. Decomposition of measured Raman and infrared (IR) spectra of aqueous solutions using multivariate curve resolution (MCR) and related methods may be used to obtain solute-correlated spectra revealing solute-induced perturbations of water structure, such as changes in water hydrogen-bond strength, tetrahedral order, and the presence of dangling (non-hydrogen-bonded) OH groups. More generally, vibrational-MCR may be applied to both aqueous and nonaqueous solutions, including multicomponent mixtures, to quantify solvent-mediated interactions between oily, polar, and ionic solutes, in both dilute and crowded fluids. Combining vibrational-MCR with emerging theoretical modeling strategies promises synergetic advances in the predictive understanding of multiscale self-assembly processes of both biological and technological interest.

Hydrated Excess Protons in Acetonitrile/Water Mixtures: Solvation Species and Ultrafast Proton Motions

The solvation structure of protons in aqueous media is highly relevant to electric properties and to proton transport in liquids and membranes. At ambient temperature, polar liquids display structural fluctuations on femto- to picosecond time scales with a direct impact on proton solvation. We use two-dimensional infrared (2D-IR) spectroscopy to follow proton dynamics in acetonitrile/water mixtures with the Zundel cation H₅O₂⁺ prepared in neat acetonitrile as a benchmark. The 2D-IR spectra of the proton transfer mode of H₅O₂⁺ demonstrate stochastic large-amplitude motions in the double-minimum proton potential, driven by fluctuating electric fields. In all cases, the excess proton is embedded in a water dimer, forming an H₅O₂⁺ complex as the major solvation species. This observation is rationalized by quantum mechanics/molecular mechanics molecular dynamics simulations including up to four water molecules embedded in acetonitrile. The Zundel motif interacts with its closest water neighbor in an H₇O₃⁺ unit without persistent proton localization.

Neat Water-Vapor Interface: Proton Continuum and the Nonresonant Background

Whether the surface of neat water is "acidic" or "basic" remains an active and controversial field of research. Most of the experimental evidence supporting the preferential adsorption of H₃O⁺ ions stems from nonlinear optical spectroscopy methods typically carried out at extreme pH conditions (pH < 1). Here, we use vibrational sum frequency spectroscopy (VSFS) to target the "proton continuum", an unexplored frequency range characteristic of hydrated protons and hydroxide ions. The VSFS spectra of neat water show a broad and nonzero signal intensity between 1700 and 3000 cm^-1 in the three different polarization combinations examined. By comparing the SF response of water with that from dilute HCl and NaOH aqueous solutions, we conclude the intensity does not originate from either adsorbed H₃O⁺ or OH^- ions. Contributions from the nonresonant background are then critically considered by comparing the experimental results with many-body molecular dynamics (MB-MD) simulated spectra.

Energy Relaxation and Structural Dynamics of Protons in Water/DMSO Mixtures

We investigate the structure and dynamics of proton solvation structures in mixed water/dimethyl sulfoxide (DMSO) solvents using two-color mid-infrared femtosecond pump–probe spectroscopy. At a water fraction below 20%, protons are mainly solvated as (DMSO-H)⁺ and (DMSO-H)⁺-H₂O structures. We find that excitation of the OH-stretch vibration of the proton in (DMSO-H)⁺-H₂O structures leads to an ultrafast contraction of the hydrogen bond between (DMSO-H)⁺ and H₂O. This excited state relaxes rapidly with T₁ = 95 ± 10 fs and leads in part to a strong local heating effect and in part to predissociation of the protonated cluster into (DMSO-H)⁺ and water monomers.