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

Classical Models of Hydroxide for Proton Hopping Simulations

Hydronium (H3O+) and hydroxide (OH-) ions perform structural diffusion in water via sequential proton transfers ("Grotthuss hopping"). This phenomenon can be accounted for by interspersing stochastic proton transfer events in classical molecular dynamics simulations. The implementation of OH--mediated proton hopping is particularly challenging because classical force fields are known to produce overcoordinated solvation structures around the OH- ion. Here, we first explore the ability of two-particle point-charge models to reproduce both the solvation free energy and coordination number in TIP3P water. We find that this is possible only with unphysical changes in the nonbonded parameters which create problems in proton hopping simulations. We then construct a classical OH- model with the charge of oxygen distributed among three auxiliary particles. This model favors a lower coordination number by accepting three hydrogen bonds and weakly donating one. The model was implemented in the MOBHY module of the CHARMM program and was fit to reproduce the experimental aqueous diffusion coefficient of OH-. This parameterization gave reasonable electrophoretic mobilities and the expected accelerated transport under nanoconfinement.

How liquids charge the superhydrophobic surfaces

Liquid-solid contact electrification (CE) is essential to diverse applications. Exploiting its full implementation requires an in-depth understanding and fine-grained control of charge carriers (electrons and/or ions) during CE. Here, we decouple the electrons and ions during liquid-solid CE by designing binary superhydrophobic surfaces that eliminate liquid and ion residues on the surfaces and simultaneously enable us to regulate surface properties, namely work function, to control electron transfers. We find the existence of a linear relationship between the work function of superhydrophobic surfaces and the as-generated charges in liquids, implying that liquid-solid CE arises from electron transfer due to the work function difference between two contacting surfaces. We also rule out the possibility of ion transfer during CE occurring on superhydrophobic surfaces by proving the absence of ions on superhydrophobic surfaces after contact with ion-enriched acidic, alkaline, and salt liquids. Our findings stand in contrast to existing liquid-solid CE studies, and the new insights learned offer the potential to explore more applications.

Overcoming Low C2+ Yield in Acidic CO2 Electroreduction: Modulating Local Hydrophobicity for Enhanced Performance

Operating electrochemical CO2 reduction reaction (CO2RR) in acidic media has garnered considerable attention due to its sustainable electrolyte cycling and stable performance. Nevertheless, the severe parasitic hydrogen evolution reaction (HER) and decayed multi-carbon species (C2+) yield still hampers efficient CO2RR in acid. Here, this work investigates the influence of local hydrophobicity on the acidic CO2RR. By employing direct electrodeposition, the hydrophobicity of the catalyst layer can be finely tuned over a wide range without additive. It is revealed that the hydrophobic microenvironment significantly suppressed HER, improved CO2RR performance and boosted C2+ yield. A Faradaic efficiency (FE) of ≈74% for C2+ is achieved in pH = 2 on electrodeposited copper with a highly hydrophobic environment. Moreover, this phenomenon can be extended to industrial application. An ≈81% total FE for the CO2RR, along with a ≈62% FE for C2+ species, is achieved even with commercial copper. Remarkably, the system exhibited stable operation for a continuous period exceeding 50 h at an industrially applied current density of 300 mA cm-2. This work highlights the crucial role of interface hydrophobicity in acidic CO2RR and proposes a facile and universally applicable method for achieving efficient and stable CO2RR to high-value products in acidic media.

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.

OH-···Au Hydrogen Bond and Its Effect on the Oxygen Reduction Reaction on Au(100) in Alkaline Media

Using ab initio molecular dynamics simulations with fully solvated ions, we demonstrate that solvated OH^- forms a stable hydrogen bond with Au(100). Unlike the hydrogen bond between H₂O and Au reported previously, which is more favorable for negatively charged Au, the OH^-···Au interaction is stabilized when a small positive charge is added to the metal slab. For electro-catalysis, this means that while OH₂···Au plays a significant role in the hydrogen evolution reaction, OH^-···Au could be a significant factor in the oxygen reduction reaction in alkaline media. It also points to a fundamental difference in the mechanism of oxygen reduction between gold and platinum electrodes.

Unveiling the structure of aqueous magnesium nitrate solutions by combining X-ray diffraction and theoretical calculations

The structure of aqueous magnesium nitrate solution is gaining significant interest among researchers, especially whether contact ion pairs exist in concentrated solutions. Here, combining X-ray diffraction experiments, quantum chemical calculations and ab initio molecular dynamics simulations, we report that the [Mg(NO₃)₂] molecular structure in solution from the coexistence of a free [Mg(H₂O)₆]²⁺ octahedral supramolecular structure with a free [NO₃(H₂O)_n]^- (n = 11-13) supramolecular structure to an [Mg²⁺(H₂O)_n(NO₃^-)_m] (n = 3, 4, 5; m = 3, 2, 1) associated structure with increasing concentration. Interestingly, two hydration modes of NO₃^--the nearest neighbor hydration with a hydration distance less than 3.9 Å and the next nearest neighbor hydration with hydration distance ranging from 3.9 to 4.3 Å-were distinguished. With an increase in the solution concentration, the hydrated NO₃^- ions lost outer layer water molecules, and the hexagonal octahedral hydration structure of [Mg(H₂O)₆²⁺] was destroyed, resulting in direct contact between Mg²⁺ and NO₃^- ions in a monodentate way. As the concentration of the solution further increased, NO₃^- ions replaced water molecules in the hydration layer of Mg²⁺ to form three-ion clusters and even more complex chains or linear ion clusters.

Controlling Hydronium Diffusivity in Model Proton Exchange Membranes

Fuel-cell-based proton exchange membranes (PEMs) show great potential as cost-effective and clean energy conversion devices. In our recent work, we found that for the low-hydrated model PEMs with a inhomogeneous water distribution and a sulfonate anionic functional end group (SO₃^-), the H₃O⁺ reacts with SO₃^- according to SO₃^- + H₃O⁺ ↔ SO₃H + H₂O, indicating that the anions in PEMs become active participants in the hydronium diffusion. In this work, we use fully atomistic ab initio molecular dynamics simulations to elucidate the optimal conditions that would promote the participation of SO₃^- in the hydronium diffusion mechanism by increasing the H₃O⁺/SO₃^- reactivity, thus increasing the hydronium diffusivity along the cell. The results presented in this work allow us to suggest a set of design rules for creating novel, highly conductive PEMs operating at high temperatures under a nonuniform water distribution using a linker/anion with a relatively high pK_a such as (CH₂)₂SO₃. We expect that the discovery of these key design principles will play an important role in the synthesis of high-performing materials for emerging PEM-based fuel cell technologies.

Structure and solvation dynamics of the hydroxide ion in ice-like water clusters: a CCSD(T) and car-parrinello molecular dynamics study

Using MP2, CCSD(T) electronic structure theory and ab initio molecular dynamics simulations, we explore the structure, solvation dynamics and vibrational spectra of OH^-(H₂O)_n clusters. Our study reports new cubic and fused cubic global minima structures of OH^-(H₂O)_n for n = 8-26 with surface and interior solvation arrangements. In the case of OH^-(H₂O)₂₆, we show that MP2 and CCSD(T) calculations predict global minima structures with the hydroxide ion occupying the interior region of a densely packed cubic cluster that is secured by ionic hydrogen bonds. More importantly, results from ab initio molecular dynamics simulations of OH^-(H₂O)₂₆ demonstrate that the hydroxide ion remains in the cluster interior and hexa-coordinated, irrespective of the temperature, up to around 175 K, then incrementally transitions from a surface-exposed penta- (170-200 K), to a tetra- (225 K) to a tri-coordinated OH^-(H₂O)₃ structure at 300 K. Building on our temperature-dependent vibrational power spectra, we are also able to disentangle structure and temperature effects on individual spectral contributions arising from water molecules located in the inner and outer shell of OH^-(H₂O)₂₆. Some of these theoretical results provide valuable guidance for the interpretation of IRMPD spectra of small hydroxide-water clusters, but there are also several intriguing implications of these results, in particular, for the solvation of the OH^- ion at the surface of water nanodroplets and aqueous interfaces.

ABEEM/MM OH- Models for OH-(H2O)n Clusters and Aqueous OH-: Structures, Charge Distributions, and Binding Energies

Based on the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), two fluctuating charge models of OH^--water system were proposed. The difference between these two models is whether there is charge transfer between OH^- and its first-shell water molecules. The structures, charge distributions, charge transfer, and binding energies of the OH^-(H₂O)_n (n = 1-8, 10, 15, 23) clusters were studied by these two ABEEM/MM models, the OPLS/AA force field, the OPLS-SMOOTH/AA force field, and the QM methods. The results demonstrate that two ABEEM/MM models can search out all stable structures just as the QM methods, and the structures and charge distributions agree well with those from the QM calculations. The structures, the charge transfer, and the strength of hydrogen bonds in the first hydration shell are closely related to the coordination number of OH^-. Molecular dynamics simulations on the aqueous OH^- solution are performed at 298 and 278 K using ABEEM/MM-I model. The MD results show that the populations of three-, four-, and five-coordinated OH^- are 29.6%, 67.1%, and 3.4% at 298 K, respectively, and those of two-, three-, four-, and five-coordinated OH^- are 10.8%, 44.9%, 39.2%, and 4.9% at 278 K, respectively; the average hydrogen bond lengths and the hydrogen bond angle in the first shell increase with the temperature decreasing.

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

Charge injection in terms of lone pairs ':', protons, and ions upon acid and base solvation mediates the hydrogen bonding network and properties of Lewis solutions, and is ubiquitously important in many subject areas of Chemical Physics. This work features the recent progress and future trends in this aspect with a focus on the solute-solvent interactions and hydrogen bond (O:H-O or HB) transition from the vibration mode of ordinary water to the hydrating states. A combination of the O:H-O bond cooperativity notion, differential phonon spectrometrics, calorimetric detection, and quantum computations clarified the solute capabilities of O:H-O bond transition in HX and YOH (X = Cl, Br, I and Y = Li, Na, K) solutions. The H+ and the lone pair do not stay alone to move or shuttle freely between adjacent H2O molecules, but they are attached to a H2O molecule to form (H3O+ and OH-)·4H2O tetrahedral motifs, which transits an O:H-O bond into the H↔H anti-HB point breaker in acidic solutions and into the O:⇔:O super-HB compressor and polarizer in basic solutions, respectively. H↔H disrupts the solvent network and surface stress, having the same effect of liquid heating on HB bond relaxation and thermal fluctuation on surface stress. The O:⇔:O compression lengthens and weakens the solute H-O bond, which heats up the solution during solvation. The H-O bonds due to H3O+ contract by 3% and due to OH- shrink by 10%. The Y+ and X- ions perform in the same manner as they do in salt solutions to form hydration shells through electrostatic polarization and hydrating H2O dipolar screen shielding. Focusing more on the O:H-O bond transition would be even more promising and revealing than on the manner and mobility of lone pair and proton transportation.

Effect of Basicity on the Hydrolysis of the Bi(III) Aqua Ion in Solution: An Ab Initio Molecular Dynamics Study

Hydrolysis of the Bi(III) aqua ion under a range of solution conditions has been studied by means of ab initio molecular dynamics simulations. While the Bi(III) aqua ion is stable in pure water, there is an increasing degree of hydrolysis with the number of hydroxide anions in the medium. This is accompanied by a monotonic decrease of the total coordination number to an asymptotic value of ∼6, reached under extreme basicity conditions. Comparison of the simulated Bi(III) hydrolyzed species with the experimental species distribution at different degrees of basicity suggests that, at the PBE/DFT level of theory here employed, liquid water shows an overly acidic character. Predictions of theoretical EXAFS and XANES spectra were generated from the AIMD trajectories for different Bi hydrolyzed species, [Bi(HO)_m(H₂O)_n]^3-m+, m = 0-3 and n = 7-2. Comparison with available experimental spectra is presented. Spectral features joined to the degree of hydrolysis and hydration are analyzed.

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.