Concerted hydrogen-bond dynamics in the transport mechanism of the hydrated proton: a first-principles molecular dynamics study

Reduced rovibrational coupling Cartesian dynamics for semiclassical calculations: Application to the spectrum of the Zundel cation

We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from the work of Huang et al. [J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows for the physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid initially exciting lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of ∼29 cm^-1 with respect to available multiconfiguration time-dependent Hartree calculations and MAE ∼ 14 cm^-1 when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster 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.

Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide Salt on the Deprotonation of Formic Acid in Aqueous Solution

The deprotonation of formic acid is investigated using metadynamics in tandem with Born-Oppenheimer molecular dynamics simulations. We compare our findings for formic acid in pure water with previous studies before examining formic acid in aqueous solutions of lithium bromide. We carefully consider different definitions for the collective variable(s) used to drive the metadynamics, emphasizing that the variables used must include all of the possible reactive atoms in the system, in this case carboxylate oxygens and water hydrogens. This ensures that all the various possible proton exchange events can be accommodated and the collective variable(s) can distinguish the protonated and deprotonated states, even over rather long ab initio simulation runs (ca. 200-300 ps). Our findings show that the formic acid deprotonation barrier and the free energy of the deprotonated state are higher in concentrated lithium bromide, in agreement with the available experimental data for acids in salt solution. We show that the presence of Br^- in proximity to the formic acid hydroxyl group effectively inhibits deprotonation. Our study extends previous work on acid deprotonation in pure water and at air-water interfaces to more complex multicomponent systems of importance in atmospheric and marine chemistry.

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.

Proton transfer in hydrogen-bonded degenerate systems of water and ammonia in metal-organic frameworks

Porous crystalline metal-organic frameworks (MOFs) or porous coordination polymers (PCPs) are emerging as a new class of proton conductors with numerous investigations. Some of the MOFs exhibit an excellent proton-conducting performance (higher than 10-2 S cm-1) originating from the interesting hydrogen(H)-bonding networks with guest molecules, where the conducting medium plays a crucial role. In the overwhelming majority of MOFs, the conducting medium is H2O because of its degenerate conjugate acid-base system (H3O+ + H2O ⇔ H2O + H3O+ or OH- + H2O ⇔ H2O + OH-) and the efficient H-bonding ability through two proton donor and two acceptor sites with a tetrahedral geometry. Considering the systematic molecular similarity to water, ammonia (NH3; NH4 + + NH3 ⇔ NH3 + NH4 +) is promising as the next proton-conducting medium. In addition, there are few reports on NH3-mediated proton conductivity in MOFs. In this perspective, we provide overviews of the degenerate water (hydronium or hydroxide)- or ammonia (ammonium)-mediated proton conduction system, the design strategies for proton-conductive MOFs, and the conduction mechanisms.

Correlated dynamics in aqueous proton diffusion

The aqueous proton displays an anomalously large diffusion coefficient that is up to 7 times that of similarly sized cations. There is general consensus that the proton achieves its high diffusion through the Grotthuss mechanism, whereby protons hop from one molecule to the next. A main assumption concerning the extraction of the timescale of the Grotthuss mechanism from experimental results has been that, on average, there is an equal probability for the proton to hop to any of its neighboring water molecules. Herein, we present ab initio simulations that show this assumption is not generally valid. Specifically, we observe that there is an increased probability for the proton to revert back to its previous location. These correlations indicate that the interpretation of the experimental results need to be re-examined and suggest that the timescale of the Grotthuss mechanism is significantly shorter than was previously thought.

Unusual Proton Transfer Kinetics in Water at the Temperature of Maximum Density

Water exhibits numerous anomalous properties, many of which remain poorly understood. One of its intriguing behaviors is that it exhibits a temperature of maximum density (TMD) at 4 °C. We provide here new experimental evidence for hitherto unknown abrupt changes in proton transfer kinetics at the TMD. In particular, we show that the lifetime of OH^{-} ions has a maximum at this temperature, in contrast to hydronium ions. Furthermore, base-catalyzed proton transfer shows a sharp local minimum at this temperature, and activation energies change abruptly as well. The measured lifetimes agree with earlier theoretical predictions as the temperature approaches the TMD. Similar results are also found for heavy water at its own TMD. These findings point to a high propensity of forming fourfold coordinated OH^{-} solvation complexes at the TMD, underlining the asymmetry between hydroxide and hydronium transport. These results could help to further elucidate the unusual properties of water and related liquids.

Broadband 2D IR spectroscopy reveals dominant asymmetric H5O2+ proton hydration structures in acid solutions

Given the critical role of the aqueous excess proton in redox chemistry, determining its structure and the mechanism of its transport in water are intense areas of experimental and theoretical research. The ultrafast dynamics of the proton's hydration structure has made it extremely challenging to study experimentally. Using ultrafast broadband two-dimensional infrared spectroscopy, we show that the vibrational spectrum of the aqueous proton is fully consistent with a protonated water complex broadly defined as a Zundel-like H₅O₂⁺ motif. Analysis of the inhomogeneously broadened proton stretch two-dimensional lineshape indicates an intrinsically asymmetric, low-barrier O-H⁺-O potential that exhibits surprisingly persistent distributions in both its asymmetry and O-O distance. This structural characterization has direct implications for the extent of delocalization exhibited by a proton's excess charge and for the possible mechanisms of proton transport in water.

Decoding the spectroscopic features and time scales of aqueous proton defects

Acid solutions exhibit a variety of complex structural and dynamical features arising from the presence of multiple interacting reactive proton defects and counterions. However, disentangling the transient structural motifs of proton defects in the water hydrogen bond network and the mechanisms for their interconversion remains a formidable challenge. Here, we use simulations treating the quantum nature of both the electrons and nuclei to show how the experimentally observed spectroscopic features and relaxation time scales can be elucidated using a physically transparent coordinate that encodes the overall asymmetry of the solvation environment of the proton defect. We demonstrate that this coordinate can be used both to discriminate the extremities of the features observed in the linear vibrational spectrum and to explain the molecular motions that give rise to the interconversion time scales observed in recent nonlinear experiments. This analysis provides a unified condensed-phase picture of the proton structure and dynamics that, at its extrema, encompasses proton sharing and spectroscopic features resembling the limiting Eigen [H₃O(H₂O)₃]⁺ and Zundel [H(H₂O)₂]⁺ gas-phase structures, while also describing the rich variety of interconverting environments in the liquid phase.

Spectroscopy of Proton Coordination with Ethylenediamine

Protonated ethylenediamine monomer, dimer, and trimer were produced in the gas phase by an electrical discharge/supersonic expansion of argon seeded with ethylenediamine (C₂H₈N₂, en) vapor. Infrared spectra of these ions were measured in the region from 1000 to 4000 cm^-1 using laser photodissociation and argon tagging. Computations at the CBS-QB3 level were performed to explore possible isomers and understand the infrared spectra. The protonated monomer exhibits a gauche conformation and an intramolecular hydrogen bond. Its parallel shared proton vibration occurs as a broad band around 2785 cm^-1, despite the formally equivalent proton affinities of the two amino groups involved, which usually leads to low frequency bands. The barrier to intramolecular proton transfer is 2.2 kcal mol^-1 and does not vanish upon addition of the zero-point energy, unlike the related protonated ammonia dimer. The structure of the dimer is formed by chelation of the monomer's NH₃⁺ group, thereby localizing the excess proton and increasing the frequency of the intramolecular shared proton vibration to 3157 cm^-1. Other highly fluxional dimer structures with facile intermolecular proton transfer and concomitant structural reorganization were computed to lie within 2 kcal mol^-1 of the experimentally observed structure. The spectrum of the trimer is rather diffuse, and a clear assignment is not possible. However, an isomer with an intramolecular proton transfer like that of the monomer is most consistent with the experimental spectrum.

From classical to quantum and back: Hamiltonian adaptive resolution path integral, ring polymer, and centroid molecular dynamics

Path integral-based methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical simulations. To reduce this numerical effort, we recently proposed a method, based on a rigorous Hamiltonian formulation, which restricts the quantum modeling to a small but relevant spatial region within a larger reservoir where particles are treated classically. In this work, we extend this idea and show how it can be implemented along with state-of-the-art path integral simulation techniques, including path-integral molecular dynamics, which allows for the calculation of quantum statistical properties, and ring-polymer and centroid molecular dynamics, which allow the calculation of approximate quantum dynamical properties. To this end, we derive a new integration algorithm that also makes use of multiple time-stepping. The scheme is validated via adaptive classical-path-integral simulations of liquid water. Potential applications of the proposed multiresolution method are diverse and include efficient quantum simulations of interfaces as well as complex biomolecular systems such as membranes and proteins.

Picosecond Proton Transfer Kinetics in Water Revealed with Ultrafast IR Spectroscopy

Aqueous proton transport involves the ultrafast interconversion of hydrated proton species that are closely linked to the hydrogen bond dynamics of water, which has been a long-standing challenge to experiments. In this study, we use ultrafast IR spectroscopy to investigate the distinct vibrational transition centered at 1750 cm^-1 in strong acid solutions, which arises from bending vibrations of the hydrated proton complex. Broadband ultrafast two-dimensional IR spectroscopy and transient absorption are used to measure vibrational relaxation, spectral diffusion, and orientational relaxation dynamics. The hydrated proton bend displays fast vibrational relaxation and spectral diffusion timescales of 200-300 fs; however, the transient absorption anisotropy decays on a remarkably long 2.5 ps timescale, which matches the timescale for hydrogen bond reorganization in liquid water. These observations are indications that the bending vibration of the aqueous proton complex is relatively localized, with an orientation that is insensitive to fast hydrogen bonding fluctuations and dependent on collective structural relaxation of the liquid to reorient. We conclude that the orientational relaxation is a result of proton transfer between configurations that are well described by a Zundel-like proton shared between two flanking water molecules.

Interfacial water molecules at biological membranes: Structural features and role for lateral proton diffusion

Proton transport at water/membrane interfaces plays a fundamental role for a myriad of bioenergetic processes. Here we have performed ab initio molecular dynamics simulations of proton transfer along two phosphatidylcholine bilayers. As found in previous theoretical studies, the excess proton is preferably located at the water/membrane interface. Further, our simulations indicate that it interacts not only with phosphate head groups, but also with water molecules at the interfaces. Interfacial water molecules turn out to be oriented relative to the lipid bilayers, consistently with experimental evidence. Hence, the specific water-proton interaction may help explain the proton mobility experimentally observed at the membrane interface.

Mechanisms of photoexcitation and photoionization in small water clusters

The S₀ → S₁ excitation leads to strong polarization and formation of [(H₂O)₂]⁺˙ from which both photoexcited and photoionized products are generated.

Towards a dissociative SPC-like water model II. The impact of Lennard-Jones and Buckingham non-coulombic forces

The dissociative water potential by Garofalini and coworkers has been re-formulated in the framework of the widely employed Lennard-Jones and Buckingham potentials, enhancing the transferability of the model to third party simulation programs.

Spectral signatures of proton delocalization in H+(H2O)n=1−4 ions

Vibrational couplings in protonated water clusters are described by harmonic analysis, vibrational perturbation theory (VPT2) and diffusion Monte Carlo (DMC) approaches.

Classical Molecular Dynamics with Mobile Protons

An important limitation of standard classical molecular dynamics simulations is the inability to make or break chemical bonds. This restricts severely our ability to study processes that involve even the simplest of chemical reactions, the transfer of a proton. Existing approaches for allowing proton transfer in the context of classical mechanics are rather cumbersome and have not achieved widespread use and routine status. Here we reconsider the combination of molecular dynamics with periodic stochastic proton hops. To ensure computational efficiency, we propose a non-Boltzmann acceptance criterion that is heuristically adjusted to maintain the correct or desirable thermodynamic equilibria between different protonation states and proton transfer rates. Parameters are proposed for hydronium, Asp, Glu, and His. The algorithm is implemented in the program CHARMM and tested on proton diffusion in bulk water and carbon nanotubes and on proton conductance in the gramicidin A channel. Using hopping parameters determined from proton diffusion in bulk water, the model reproduces the enhanced proton diffusivity in carbon nanotubes and gives a reasonable estimate of the proton conductance in gramicidin A.

Predicting the Ionic Product of Water

We present a first-principles calculation and mechanistic characterization of the ion product of liquid water (K _W ), based on Quantum Cluster Equilibrium (QCE) theory with a variety of ab initio and density functional methods. The QCE method is based on T-dependent Boltzmann weighting of different-sized clusters and consequently enables the observation of thermodynamically less favored and therefore low populated species such as hydronium and hydroxide ions in water. We find that common quantum chemical methods achieve semi-quantitative accuracy in predicting K _W and its T-dependence. Dominant ion-pair water clusters of the QCE equilibrium distribution are found to exhibit stable 2-coordinate buttress-type motifs, all with maximally Grotthus-ordered H-bond patterns that successfully prevent recombination of hydronium and hydroxide ions at 3-coordinate bridgehead sites. We employ standard quantum chemistry techniques to describe kinetic and mechanistic aspects of ion-pair formation, and we obtain NBO-based bonding indices to characterize other electronic, structural, spectroscopic, and reactive properties of cluster-mediated ionic dissociation.

Spectroscopic snapshots of the proton-transfer mechanism in water

The Grotthuss mechanism explains the anomalously high proton mobility in water as a sequence of proton transfers along a hydrogen-bonded (H-bonded) network. However, the vibrational spectroscopic signatures of this process are masked by the diffuse nature of the key bands in bulk water. Here we report how the much simpler vibrational spectra of cold, composition-selected heavy water clusters, D⁺(D₂O)_n, can be exploited to capture clear markers that encode the collective reaction coordinate along the proton-transfer event. By complexing the solvated hydronium "Eigen" cluster [D₃O⁺(D₂O)₃] with increasingly strong H-bond acceptor molecules (D₂, N₂, CO, and D₂O), we are able to track the frequency of every O-D stretch vibration in the complex as the transferring hydron is incrementally pulled from the central hydronium to a neighboring water molecule.

Nuclear Quantum Effects in Water Reorientation and Hydrogen-Bond Dynamics

We combine classical and ring polymer molecular dynamics simulations with the molecular jump model to provide a molecular description of the nuclear quantum effects (NQEs) on water reorientation and hydrogen-bond dynamics in liquid H₂O and D₂O. We show that while the net NQE is negligible in D₂O, it leads to a ∼13% acceleration in H₂O dynamics compared to a classical description. Large angular jumps-exchanging hydrogen-bond partners-are the dominant reorientation pathway (just as in a classical description); the faster reorientation dynamics arise from the increased jump rate constant. NQEs do not change the jump amplitude distribution, and no significant tunneling is found. The faster jump dynamics are quantitatively related to decreased structuring of the OO radial distribution function when NQEs are included. This is explained, via a jump model analysis, by competition between the effects of water's librational and OH stretch mode zero-point energies on the hydrogen-bond strength.

Facile Anhydrous Proton Transport on Hydroxyl Functionalized Graphane

Graphane functionalized with hydroxyl groups is shown to rapidly conduct protons under anhydrous conditions through a contiguous network of hydrogen bonds. Density functional theory calculations predict remarkably low barriers to diffusion of protons along a 1D chain of surface hydroxyls. Diffusion is controlled by the local rotation of hydroxyl groups, a mechanism that is very different from that found in 1D water wires in confined nanopores or in bulk water. The proton mean square displacement in the 1D chain was observed to follow Fickian diffusion rather than the expected single-file mobility. A charge analysis reveals that the charge on the proton is essentially equally shared by all hydrogens bound to oxygens, effectively delocalizing the proton.

Proton-Transfer Mechanisms at the Water-ZnO Interface: The Role of Presolvation

The dissociation of water is an important step in many chemical processes at solid surfaces. In particular, water often spontaneously dissociates near metal oxide surfaces, resulting in a mixture of H₂O, H⁺, and OH^- at the interface. Ubiquitous proton-transfer (PT) reactions cause these species to dynamically interconvert, but the underlying mechanisms are poorly understood. Here, we develop and use a reactive high-dimensional neural-network potential based on density functional theory data to elucidate the structural and dynamical properties of the interfacial species at the liquid-water-metal-oxide interface, using the nonpolar ZnO(101̅0) surface as a prototypical case. Molecular dynamics simulations reveal that water dissociation and recombination proceed via two types of PT reactions: (i) to and from surface oxide and hydroxide anions ("surface-PT") and (ii) to and from neighboring adsorbed hydroxide ions and water molecules ("adlayer-PT"). We find that the adlayer-PT rate is significantly higher than the surface-PT rate. Water dissociation is, for both types of PT, governed by a predominant presolvation mechanism, i.e., thermal fluctuations that cause the adsorbed water molecules to occasionally accept a hydrogen bond, resulting in a decreased PT barrier and an increased dissociation rate as compared to when no hydrogen bond is present. Consequently, we are able to show that hydrogen bond fluctuations govern PT events at the water-metal-oxide interface in a way similar to that in acidic and basic aqueous bulk solutions.

Structure of a model TiO2 photocatalytic interface

The interaction of water with TiO₂ is crucial to many of its practical applications, including photocatalytic water splitting. Following the first demonstration of this phenomenon 40 years ago there have been numerous studies of the rutile single-crystal TiO₂(110) interface with water. This has provided an atomic-level understanding of the water-TiO₂ interaction. However, nearly all of the previous studies of water/TiO₂ interfaces involve water in the vapour phase. Here, we explore the interfacial structure between liquid water and a rutile TiO₂(110) surface pre-characterized at the atomic level. Scanning tunnelling microscopy and surface X-ray diffraction are used to determine the structure, which is comprised of an ordered array of hydroxyl molecules with molecular water in the second layer. Static and dynamic density functional theory calculations suggest that a possible mechanism for formation of the hydroxyl overlayer involves the mixed adsorption of O₂ and H₂O on a partially defected surface. The quantitative structural properties derived here provide a basis with which to explore the atomistic properties and hence mechanisms involved in TiO₂ photocatalysis.

Reduced Acid Dissociation of Amino-Acids at the Surface of Water

We use surface-specific intensity vibrational sum-frequency generation and attenuated total reflection spectroscopy to probe the ionization state of the amino-acids l-alanine and l-proline at the air/water surface and in the bulk. The ionization state is determined by probing the vibrational signatures of the carboxylic acid group, representing the nondissociated acid form, and the carboxylate anion group, representing the dissociated form, over a wide range of pH values. We find that the carboxylic acid group deprotonates at a significantly higher pH at the surface than in the bulk.

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.

Collective proton transfer in ordinary ice: local environments, temperature dependence and deuteration effects

Ordinary ice at low temperature: what about collective nuclear quantum effects in its chiral six rings?

Towards a dissociative SPC-like water model – probing the impact of intramolecular Coulombic contributions

A modification of the dissociative Garofalini water model towards an SPC-like Coulombic formulation proved to enhance accuracy and transferability of this successful force field approach.

Effect of Ions on H-Bond Structure and Dynamics at the Quartz(101)-Water Interface

We use ab initio molecular dynamics simulations to study the effect of ions on the structure and dynamics of the quartz(101)-water interface. We study several IA (Na⁺, Rb⁺) and IIA (Mg²⁺, Sr²⁺) cations, with Cl^- as counterion, adsorbed onto acidic, neutral, and basic surface configurations at 300 and 373 K. We find that both cations and anions can bond directly to the surface and perturb the local H-bond network. The adsorbed ions promote the formation of intrasurface H-bonds, as shown by vibrational density of states and orientations of the surface silanols. Both local and global structural correlations of the interfacial H-bond network are studied using a structural definition of the H-bond and a network correlation function. We find the ions' effect on the solvent structure exhibits a complex dependence on specific surface interactions. The structure-making properties of ions are enhanced at the quartz surface, particularly for ions adsorbed without a complete hydration shell, and the structuring effect extends beyond the first solvation shell. The ions have a lesser effect on solvent structure in solution, especially in the presence of counterions. In fact, cations that are the greatest "structure makers" at the surface are the greatest "structure breakers" when in solution with a counterion. Therefore, we find the ions cannot be simply classified as "structure making" or "structure breaking". We discuss the implications of these findings for the effect of ions on the dissolution rate, surface charge, and solvent structure.

Proton transfer aiding phase transitions in oxalic acid dihydrate under pressure

Oxalic acid dihydrate, an important molecular solid in crystal chemistry, ecology and physiology, has been studied for nearly 100 years now. The most debated issues regarding its proton dynamics have arisen due to an unusually short hydrogen bond between the acid and water molecules. Using combined in situ spectroscopic studies and first-principles simulations at high pressures, we show that the structural modification associated with this hydrogen bond is much more significant than ever assumed. Initially, under pressure, proton migration takes place along this strong hydrogen bond at a very low pressure of 2 GPa. This results in the protonation of water with systematic formation of dianionic oxalate and hydronium ion motifs, thus reversing the hydrogen bond hierarchy in the high pressure phase II. The resulting hydrogen bond between a hydronium ion and a carboxylic group shows remarkable strengthening under pressure, even in the pure ionic phase III. The loss of cooperativity of hydrogen bonds leads to another phase transition at ∼ 9 GPa through reorientation of other hydrogen bonds. The high pressure phase IV is stabilized by a strong hydrogen bond between the dominant CO2 and H2O groups of oxalate and hydronium ions, respectively. These findings suggest that oxalate systems may provide useful insights into proton transfer reactions and assembly of simple molecules under extreme conditions.

Advanced Potential Energy Surfaces for Molecular Simulation

Advanced potential energy surfaces are defined as theoretical models that explicitly include many-body effects that transcend the standard fixed-charge, pairwise-additive paradigm typically used in molecular simulation. However, several factors relating to their software implementation have precluded their widespread use in condensed-phase simulations: the computational cost of the theoretical models, a paucity of approximate models and algorithmic improvements that can ameliorate their cost, underdeveloped interfaces and limited dissemination in computational code bases that are widely used in the computational chemistry community, and software implementations that have not kept pace with modern high-performance computing (HPC) architectures, such as multicore CPUs and modern graphics processing units (GPUs). In this Feature Article we review recent progress made in these areas, including well-defined polarization approximations and new multipole electrostatic formulations, novel methods for solving the mutual polarization equations and increasing the MD time step, combining linear-scaling electronic structure methods with new QM/MM methods that account for mutual polarization between the two regions, and the greatly improved software deployment of these models and methods onto GPU and CPU hardware platforms. We have now approached an era where multipole-based polarizable force fields can be routinely used to obtain computational results comparable to state-of-the-art density functional theory while reaching sampling statistics that are acceptable when compared to that obtained from simpler fixed partial charge force fields.

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.