Collective vibrations of water-solvated hydroxide ions investigated with broadband 2DIR spectroscopy

Dressing a Nonpolarizable Force Field for OH- in TIP4P/2005 Aqueous Solutions with Corrected Hirshfeld Charges

We present a rigid model for the OH- ion parametrized for binary mixtures with TIP4P/2005-type water molecules. Li+, Na+ and K+ were selected as counterions, hence mimicking the important and widely used solutions of soluble alkaline hydroxides. The optimized atomic charge distributions were obtained by scaling in a factor of 0.85 those derived from the atomic dipole corrected Hirshfeld approach. The agreement between experimental and Molecular Dynamics simulation results is remarkable for a set of properties, namely, the dependence of the density of the solutions on the hydroxide concentration and on temperature, the structure (i.e., positions of the atom-to-atom radial distribution functions and coordination numbers), the viscosity coefficients, the surface tension, or the freezing point depression. The proposed optimized potential parameters for OH- thus enlarge the set of models comprised within the Madrid-2019 force field and widen the potential applicability of the TIP4P/2005 water model in basic environments.

Surface stratification determines the interfacial water structure of simple electrolyte solutions

The distribution of ions at the air/water interface plays a decisive role in many natural processes. Several studies have reported that larger ions tend to be surface-active, implying ions are located on top of the water surface, thereby inducing electric fields that determine the interfacial water structure. Here we challenge this view by combining surface-specific heterodyne-detected vibrational sum-frequency generation with neural network-assisted ab initio molecular dynamics simulations. Our results show that ions in typical electrolyte solutions are, in fact, located in a subsurface region, leading to a stratification of such interfaces into two distinctive water layers. The outermost surface is ion-depleted, and the subsurface layer is ion-enriched. This surface stratification is a key element in explaining the ion-induced water reorganization at the outermost air/water interface.

Microscopic Mechanism of Proton Transfer in Pure Water under Ambient Conditions

Water molecules and the associated proton transfer (PT) are prevalent in chemical and biological systems and have been a hot research topic. Spectroscopic characterization and ab initio molecular dynamics (AIMD) simulations have previously revealed insights into acidic and basic liquids. Presumably, the situation in the acidic/basic solution is not necessarily the same as in pure water; in addition, the autoionization constant for water is only 10^-14 under ambient conditions, making the study of PT in pure water challenging. To overcome this issue, we modeled periodic water box systems containing 1000 molecules for tens of nanoseconds based on a neural network potential (NNP) with quantum mechanical accuracy. The NNP was generated by training a dataset containing the energies and atomic forces of 17 075 configurations of periodic water box systems, and these data points were calculated at the MP2 level that considers electron correlation effects. We found that the size of the system and the duration of the simulation have a significant impact on the convergence of the results. With these factors considered, our simulations showed that hydronium (H₃O⁺) and hydroxide (OH^-) ions in water have distinct hydration structures, thermodynamic and kinetic properties, e.g., the longer-lasting and more stable hydrated structure of OH^- ions than that of H₃O⁺, as well as a significantly higher free energy barrier for the OH^--associated PT than that of H₃O⁺, leading the two to exhibit completely different PT behaviors. Given these characteristics, we further found that PT via OH^- ions tends not to occur multiple times or between many molecules. In contrast, PT via H₃O⁺ can synergistically occur among multiple molecules and prefers to adopt a cyclic pattern among three water molecules, while it occurs mostly in a chain pattern when more water molecules are involved. Therefore, our studies provide a detailed and solid microscopic explanation for the PT process in pure water.

Free energy difference to create the M-OH* intermediate of the oxygen evolution reaction by time-resolved optical spectroscopy

Theoretical descriptors differentiate the catalytic activity of materials for the oxygen evolution reaction by the strength of oxygen binding in the reactive intermediate created upon electron transfer. Recently, time-resolved spectroscopy of a photo-electrochemically driven oxygen evolution reaction followed the vibrational and optical spectra of this intermediate, denoted M-OH^*. However, these inherently kinetic experiments have not been connected to the relevant thermodynamic quantities. Here we discover that picosecond optical spectra of the Ti-OH^* population on lightly doped SrTiO₃ are ordered by the surface hydroxylation. A Langmuir isotherm as a function of pH extracts an effective equilibrium constant relatable to the free energy difference of the first oxygen evolution reaction step. Thus, time-resolved spectroscopy of the catalytic surface reveals both kinetic and energetic information of elementary reaction steps, which provides a critical new connection between theory and experiment by which to tailor the pathway of water oxidation and other surface reactions.

Mode-Selective Vibrational Energy Transfer Dynamics in 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX) Thin Films

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. Herein, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N-N and axial NO₂ stretching modes with the long-lived, excited phonon bath.

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.

Hydration interactions beyond the first solvation shell in aqueous phenolate solution

We investigate the orientational dynamics of water molecules solvating phenolate ions using ultrafast vibrational spectroscopy and density functional theory-based molecular dynamics simulations.

Generation and Applications of the Hydroxide Trihydrate Anion, [OH(OH2 )3 ]- , Stabilized by a Weakly Coordinating Cation

The reaction of a strongly basic phosphazene (Schwesinger base) with water afforded the corresponding metastable hydroxide trihydrate [OH(OH₂)₃]⁻ salt. This is the first hydroxide solvate that is not in contact with a cation and furthermore one of rare known water‐stabilized hydroxide anions. Thermolysis in vacuum results in the decomposition of the hydroxide salt and quantitative liberation of the free phosphazene base. This approach was used to synthesize the Schwesinger base from its hydrochloride salt after anion exchange in excellent yields of over 97 %. This deprotonation method can also be used for the phosphazene‐base‐catalyzed preparation of the Ruppert–Prakash reagent Me₃SiCF₃ using fluoroform (HCF₃) as the trifluoromethyl building block and sodium hydroxide as the formal deprotonation agent.

Unprecedented O:⇔:O compression and H↔H fragilization in Lewis solutions

Charge injection in terms of protons, lone pairs, cations and anions by acid and base solvation mediates the HB network and properties of Lewis solutions through H↔H fragilization, O:⇔:O compression and polarization, ionic polarization and hydrating H₂O dipolar screen shielding, anion–anion repulsion, compressed solvent H–O bond elongation and undercoordinated solute H–O bond contraction.

IR spectral assignments for the hydrated excess proton in liquid water

The local environmental sensitivity of infrared (IR) spectroscopy to a hydrogen-bonding structure makes it a powerful tool for investigating the structure and dynamics of excess protons in water. Although of significant interest, the line broadening that results from the ultrafast evolution of different solvated proton-water structures makes the assignment of liquid-phase IR spectra a challenging task. In this work, we apply a normal mode analysis using density functional theory of thousands of proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multistate empirical valence bond proton model (MS-EVB 3.2). These calculations are used to obtain a vibrational density of states and IR spectral density, which are decomposed on the basis of solvated proton structure and the frequency dependent mode character. Decompositions are presented on the basis of the proton sharing parameter δ, often used to distinguish Eigen and Zundel species, the stretch and bend character of the modes, the mode delocalization, and the vibrational mode symmetry. We find there is a wide distribution of vibrational frequencies spanning 1200-3000 cm^-1 for every local proton configuration, with the region 2000-2600 cm^-1 being mostly governed by the distorted Eigen-like configuration. We find a continuous red shift of the special-pair O⋯H⁺⋯O stretching frequency, and an increase in the flanking water bending intensity with decreasing δ. Also, we find that the flanking water stretch mode of the Zundel-like species is strongly mixed with the flanking water bend, and the special pair proton oscillation band is strongly coupled with the bend modes of the central H₅O2+moiety.

Ultrafast Two-Dimensional Infrared Spectroscopy of a Quasifree Rotor: J Scrambling and Perfectly Anticorrelated Cross Peaks

Ultrafast two-dimensional infrared (2DIR) spectra of the N_{2}O ν_{3} mode in moderately dense SF_{6} gas exhibit complex line shapes with diagonal and antidiagonal features in contrast to condensed phase vibrational 2DIR spectroscopy. Observed spectra for this quasifree rotor system are well captured by a model that includes all 36 possible rovibrational pathways and treats P (ΔJ=-1) and R (ΔJ=+1) branch resonances as distinct Kubo line shape features. Transition frequency correlation decay is due to J scrambling within one to two gas collisions at each density. Studies of supercritical solvation and relaxation at high pressure and temperature are enabled by this methodology.

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.

Dimension of discrete variable representation for mixed quantum/classical computation of three lowest vibrational states of OH stretching in liquid water

Three low-lying vibrational states of molecular systems are responsible for the signals of linear and third-order nonlinear vibrational spectroscopies. Theoretical studies based on mixed quantum/classical calculations provide a powerful way to analyze those experiments. A statistically meaningful result can be obtained from the calculations by solving the vibrational Schrödinger equation over many numbers of molecular configurations. The discrete variable representation (DVR) method is a useful technique to calculate vibrational eigenstates subject to an arbitrary anharmonic potential surface. Considering the large number of molecular configurations over which the DVR calculations are repeated, the calculations are desired to be optimized in balance between the cost and accuracy. We determine a dimension of the DVR method which appears to be optimum for the calculations of the three states of molecular vibrations with anharmonic strengths often found in realistic molecular systems. We apply the numerical technique to calculate the local OH stretching frequencies of liquid water, which are well known to be widely distributed due to the inhomogeneity in molecular configuration, and found that the frequencies of the 0-1 and 1-2 transitions are highly correlated. An empirical relation between the two frequencies is suggested and compared with the experimental data of nonlinear IR spectroscopies.

How the Zundel (H5O2+) Potential Can Be Used to Predict the Proton Stretch and Bend Frequencies of Larger Protonated Water Clusters

From a series of seminal experiments on the IR spectra of protonated water clusters and associated theoretical analyses, it is clear that the energies and spectral features of the proton stretch and bend modes are very sensitive functions of the cluster size. Here we show that this dynamic range can be understood by examining the sensitivity of these modes in the potential of the Zundel cation, H₅O₂⁺, as the separation of the two water monomers is varied. As this distance increases, the proton increasingly localizes on a monomer, and this is encoded in the IR spectrum of the proton vibrational modes. The quantitative predictions from this simple correlation are verified for the H₇O₃⁺ and H₉O₄⁺ clusters, for which new benchmark harmonic frequencies are reported. The predictions are also in good accord with trends seen experimentally and previous calculations for these and five other clusters, including H⁺(H₂O)₂₁.

Characterizing the local solvation environment of OH(-) in water clusters with AIMD

In this work, we use ab initio molecular dynamics coupled with metadynamics to explore and characterize the glassy potential energy landscape of the OH(-) in a 20 and 48 water cluster. The structural, energetic, and topological properties of OH(-) are characterized for both clusters and the molecular origins of the IR signatures are examined. We find that in both the small and large clusters, the OH(-) can donate or accept a varying number of hydrogen bonds confirming that the amphiphilic character does not depend on cluster size. However, we highlight some important differences found between the energetic and topological properties of both families of clusters which may have implications on understanding the changes in the solvation structure of OH(-) between bulk and interfacial environments. By studying the IR spectra of smaller subsets of molecules within the 20 water molecule cluster, we find that the IR spectrum of the bare OH(-) as well as the water molecule donating a strong hydrogen bond to it exhibits characteristic absorption along the amphiphilic band between 1500 and 3000 cm(-1) at positions very similar to those found for the entire hydroxide cluster. The results presented here will be useful in the calibration and improvement of both ab initio and semi-empirical methods to model this complex anion.

Protons and Hydroxide Ions in Aqueous Systems

Understanding the structure and dynamics of water's constituent ions, proton and hydroxide, has been a subject of numerous experimental and theoretical studies over the last century. Besides their obvious importance in acid-base chemistry, these ions play an important role in numerous applications ranging from enzyme catalysis to environmental chemistry. Despite a long history of research, many fundamental issues regarding their properties continue to be an active area of research. Here, we provide a review of the experimental and theoretical advances made in the last several decades in understanding the structure, dynamics, and transport of the proton and hydroxide ions in different aqueous environments, ranging from water clusters to the bulk liquid and its interfaces with hydrophobic surfaces. The propensity of these ions to accumulate at hydrophobic surfaces has been a subject of intense debate, and we highlight the open issues and challenges in this area. Biological applications reviewed include proton transport along the hydration layer of various membranes and through channel proteins, problems that are at the core of cellular bioenergetics.

Fourier transform two-dimensional electronic-vibrational spectroscopy using an octave-spanning mid-IR probe

The development of coherent Fourier transform two-dimensional electronic-vibrational (2D EV) spectroscopy with acousto-optic pulse-shaper-generated near-UV pump pulses and an octave-spanning broadband mid-IR probe pulse is detailed. A 2D EV spectrum of a silicon wafer demonstrates the full experimental capability of this experiment, and a 2D EV spectrum of dissolved hexacyanoferrate establishes the viability of our 2D EV experiment for studying condensed phase molecular ensembles.

Vibrational dynamics of aqueous hydroxide solutions probed using broadband 2DIR spectroscopy

We employed ultrafast transient absorption and broadband 2DIR spectroscopy to study the vibrational dynamics of aqueous hydroxide solutions by exciting the O-H stretch vibrations of the strongly hydrogen-bonded hydroxide solvation shell water and probing the continuum absorption of the solvated ion between 1500 and 3800 cm(-1). We observe rapid vibrational relaxation processes on 150-250 fs time scales across the entire probed spectral region as well as slower vibrational dynamics on 1-2 ps time scales. Furthermore, the O-H stretch excitation loses its frequency memory in 180 fs, and vibrational energy exchange between bulk-like water vibrations and hydroxide-associated water vibrations occurs in ∼200 fs. The fast dynamics in this system originate in strong nonlinear coupling between intra- and intermolecular vibrations and are explained in terms of non-adiabatic vibrational relaxation. These measurements indicate that the vibrational dynamics of the aqueous hydroxide complex are faster than the time scales reported for long-range transport of protons in aqueous hydroxide solutions.

Two-dimensional vibrational-electronic spectroscopy

Two-dimensional vibrational-electronic (2D VE) spectroscopy is a femtosecond Fourier transform (FT) third-order nonlinear technique that creates a link between existing 2D FT spectroscopies in the vibrational and electronic regions of the spectrum. 2D VE spectroscopy enables a direct measurement of infrared (IR) and electronic dipole moment cross terms by utilizing mid-IR pump and optical probe fields that are resonant with vibrational and electronic transitions, respectively, in a sample of interest. We detail this newly developed 2D VE spectroscopy experiment and outline the information contained in a 2D VE spectrum. We then use this technique and its single-pump counterpart (1D VE) to probe the vibrational-electronic couplings between high frequency cyanide stretching vibrations (νCN) and either a ligand-to-metal charge transfer transition ([Fe(III)(CN)6](3-) dissolved in formamide) or a metal-to-metal charge transfer (MMCT) transition ([(CN)5Fe(II)CNRu(III)(NH3)5](-) dissolved in formamide). The 2D VE spectra of both molecules reveal peaks resulting from coupled high- and low-frequency vibrational modes to the charge transfer transition. The time-evolving amplitudes and positions of the peaks in the 2D VE spectra report on coherent and incoherent vibrational energy transfer dynamics among the coupled vibrational modes and the charge transfer transition. The selectivity of 2D VE spectroscopy to vibronic processes is evidenced from the selective coupling of specific νCN modes to the MMCT transition in the mixed valence complex. The lineshapes in 2D VE spectra report on the correlation of the frequency fluctuations between the coupled vibrational and electronic frequencies in the mixed valence complex which has a time scale of 1 ps. The details and results of this study confirm the versatility of 2D VE spectroscopy and its applicability to probe how vibrations modulate charge and energy transfer in a wide range of complex molecular, material, and biological systems.

Solvation dynamics of an ionic probe in choline chloride-based deep eutectic solvents

Study of the solvation dynamics of an ionic probe in different choline-based deep eutectic solvents shows that the process is controlled by the motions of the choline ions within the pseudo lattice formed by the solvent.

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.

Insights into hydrogen bonding via ice interfaces and isolated water

Water in a confined environment has a combination of fewer available configurations and restricted mobility. Both affect the spectroscopic signature. In this work, the spectroscopic signature of water in confined environments is discussed in the context of competing models for condensed water: (1) as a system of intramolecular coupled molecules or (2) as a network with intermolecular dipole-dipole coupled O-H stretches. Two distinct environments are used: the confined asymmetric environment at the ice surface and the near-isolated environment of water in an infrared transparent matrix. Both the spectroscopy and the environment are described followed by a perspective discussion of implications for the two competing models. Despite being a small molecule, water is relatively complex; perhaps not surprisingly the results support a model that blends inter- and intramolecular coupling. The frequency, and therefore the hydrogen-bond strength, appears to be a function of donor-acceptor interaction and of longer-range dipole-dipole alignment in the hydrogen-bonded network. The O-H dipole direction depends on the local environment and reflects intramolecular O-H stretch coupling.