What Coordinate Best Describes the Affinity of the Hydrated Excess Proton for the Air-Water Interface?

Ab Initio Characterization of the CO2-Water Interface Using Unsupervised Machine Learning for Dimensionality Reduction

Precise characterization of the supercritical CO2-water interface under high pressure and temperature conditions is crucial for the geological storage of carbon dioxide (CO2) in deep saline aquifers. Molecular dynamics (MD) simulations offer a valuable approach to gaining insight into the CO2-water interface at a microscopic level. However, no attempt has been made to characterize the CO2-water interface with the accuracy afforded by ab initio calculations. In this study, we performed ab initio MD (AIMD) simulations to investigate the structural and dynamical properties of the CO2-water interface, comparing the results with those obtained from classical force-field MD (FF-MD) simulations. Molecular orientation at the interface was well reproduced in both AIMD and FF-MD simulations. Characteristic structural fluctuations of water at the interface were unveiled by applying multidimensional scaling and time-dependent principal component analysis to the AIMD trajectories; however, they were not prominent in the FF-MD simulations. Furthermore, our study demonstrated a marked difference in the residence time of molecules in the interface region between AIMD and FF-MD simulations, indicating that time-dependent properties of the CO2-water interface strongly depend on the description of the intermolecular forces.

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system's bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

Influence of H+, OH- and salts on hydrophobic self-assembly

In spite of the ubiquity of acid/base ions and salts in biological systems, their influence on hydrophobic self-assembly remains an open question. Here we use a combined experimental and theoretical strategy to quantify the influence of H+ and OH-, as well as salts containing Li+, Na+, Cl- and Br-, on the hydrophobic self-assembly of micelles composed of neutral oily 1,2-hexanediol surfactants. The distributions of aggregate sizes, both below and above the critical micelle concentration (CMC), are determined using Raman multivariate curve resolution (Raman-MCR) spectroscopy to quantify the multi-aggregation chemical potential surface (MCPS) that drives self-assembly. The results reveal that ions have little influence on the formation of hydrophobic contact dimers but can significantly drive high-order self assembly. Moreover, the hydration-shells of oily solutes are found to expel the above salt ions and OH-, but to attract H+, with wide-ranging implications.

Tuning the Microenvironment of Water Confined in Ti3C2Tx MXene by Cation Intercalation

The local microenvironment has recently been found to play a major role in the electrocatalytic activity of nanomaterials. Modulating the microenvironment by adding alkali metal cations into the electrolyte can be used to either suppress hydrogen or oxygen evolution, thereby extending the electrochemical window of energy storage systems, or to tune the selectivity of electrocatalysts. MXenes are a large family of two-dimensional transition metal carbides, nitrides, and carbonitrides that have shown potential for use in electrochemical energy storage applications. Due to their negatively charged surfaces, MXenes can accommodate cations and water molecules between the layers. Nevertheless, the nature of the aqueous microenvironment in the MXene interlayer space is poorly understood. Here, we apply Fourier transform infrared spectroscopy (FTIR) to probe the hydrogen bonding of intercalated water in Ti3C2Tx as a function of intercalated cation and relative humidity. Substantial changes in the FTIR spectra after cation exchange demonstrate that the hydrogen bonding of water molecules confined between the MXene layers is strongly cation-dependent. Furthermore, the IR absorbance of the confined water correlates with resistivity estimated by 4-point probe measurements and interlayer distance calculated from XRD patterns. This work demonstrates that cation intercalation strongly modulates the confined microenvironment, which can be used to tune the activity or selectivity of electrochemical reactions in the interlayer space of MXenes in the future.

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.

Acidic Conditions Impact Hydrophobe Transfer across the Oil-Water Interface in Unusual Ways

Molecular dynamics simulation and enhanced free energy sampling are used to study hydrophobic solute transfer across the water-oil interface with explicit consideration of the effect of different electrolytes: hydronium cation (hydrated excess proton) and sodium cation, both with chloride counterions (i.e., dissociated acid and salt, HCl and NaCl). With the Multistate Empirical Valence Bond (MS-EVB) methodology, we find that, surprisingly, hydronium can to a certain degree stabilize the hydrophobic solute, neopentane, in the aqueous phase and including at the oil-water interface. At the same time, the sodium cation tends to "salt out" the hydrophobic solute in the expected fashion. When it comes to the solvation structure of the hydrophobic solute in the acidic conditions, hydronium shows an affinity to the hydrophobic solute, as suggested by the radial distribution functions (RDFs). Upon consideration of this interfacial effect, we find that the solvation structure of the hydrophobic solute varies at different distances from the oil-liquid interface due to a competition between the bulk oil phase and the hydrophobic solute phase. Together with an observed orientational preference of the hydroniums and the lifetime of water molecules in the first solvation shell of neopentane, we conclude that hydronium stabilizes to a certain degree the dispersal of neopentane in the aqueous phase and eliminates any salting out effect in the acid solution; i.e., the hydronium acts like a surfactant. The present molecular dynamics study provides new insight into the hydrophobic solute transfer across the water-oil interface process, including for acid and salt solutions.

Regarding the Influence of Additives and Additional Plasma-Induced Chemical Ionization on Adduct Formation in ESI/IMS/MS

Ion mobility spectrometers (IMS) separate ions based on their ion mobility, which depends mainly on collision cross-section, mass, and charge of the ions. However, the performance is often hampered in electrospray ionization (ESI) by the appearance of multiple ion mobility peaks in the spectrum for the same analyte due to clustering and additional sodium adducts. In this work, we investigate the influence of solvents and buffer additives on the detected ion mobility peaks using ESI. Additionally, we investigate the effects of an additional chemical ionization (CI) induced by plasma ionization on the ions formed by electrospray. For this purpose, we coupled our high-resolution IMS with a resolving power of R_p = 100 to a time-of-flight mass spectrometer. Depending on the analyte and the chosen additives, the ionization process can be influenced during the electrospray process. For the herbicide isoproturon, the addition of 5 mM sodium acetate results in the formation of the sodium adduct [M + Na]⁺, which is reflected in the ion mobility K₀ of 1.22 cm²/(V·s). In contrast, the addition of 5 mM ammonium acetate yields the protonated species [M + H]⁺ and a correspondingly higher K₀ of 1.29 cm²/(V·s). In some cases, as with the herbicide pyrimethanil, the addition of sodium acetate can completely suppress ionizations. By carefully choosing the solvent additive for ESI-IMS or additional CI, the formation of different ion mobility peaks can be observed. This can facilitate the assignment of ions to ion mobility peaks using IMS as a compact, stand-alone instrument, e.g., for on-site analysis.

Highly Altered State of Proton Transport in Acid Pools in Charged Reverse Micelles

Transport mechanisms of solvated protons of 1 M HCl acid pools, confined within reverse micelles (RMs) containing the negatively charged surfactant sodium bis(2-ethylhexyl) sulfosuccinate (NaAOT) or the positively charged cetyltrimethylammonium bromide (CTABr), are analyzed with reactive force field simulations to interpret dynamical signatures from TeraHertz absorption and dielectric relaxation spectroscopy. We find that the forward proton hopping events for NaAOT are further suppressed compared to a nonionic RM, while the Grotthuss mechanism ceases altogether for CTABr. We attribute the sluggish proton dynamics for both charged RMs as due to headgroup and counterion charges that expel hydronium and chloride ions from the interface and into the bulk interior, thereby increasing the pH of the acid pools relative to the nonionic RM. For charged NaAOT and CTABr RMs, the localization of hydronium near a counterion or conjugate base reduces the Eigen and Zundel configurations that enable forward hopping. Thus, localized oscillatory hopping dominates, an effect that is most extreme for CTABr in which the proton residence time increases dramatically such that even oscillatory hopping is slow.

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.