Propensity of soft ions for the air/water interface

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Effects of Charged Surfactants on Interfacial Water Structure and Macroscopic Properties of the Air-Water Interface

Surfactants are used to control the macroscopic properties of the air-water interface. However, the link between the surfactant molecular structure and the macroscopic properties remains unclear. Using sum-frequency generation spectroscopy and molecular dynamics simulations, two ionic surfactants (dodecyl trimethylammonium bromide, DTAB, and sodium dodecyl sulphate, SDS) with the same carbon chain lengths and charge magnitude (but different signs) of head groups interact and reorient interfacial water molecules differently. DTAB forms a thicker but sparser interfacial layer than SDS. It is due to the deep penetration into the adsorption zone of Br- counterions compared to smaller Na+ ones, and also due to the flip-flop orientation of water molecules. SDS alters two distinctive interfacial water layers into a layer where H+ points to the air, forming strong hydrogen bonding with the sulphate headgroup. In contrast, only weaker dipole-dipole interactions with the DTAB headgroup are formed as they reorient water molecules with H+ point down to the aqueous phase. Hence, with more molecules adsorbed at the interface, SDS builds up a higher interfacial pressure than DTAB, producing lower surface tension and higher foam stability at a similar bulk concentration. Our findings offer improved knowledge for understanding various processes in the industry and nature.

Does the Sign of Charge Affect the Surface Affinity of Simple Ions?

The role the charge sign of simple ions plays in determining their surface affinity in aqueous solutions is investigated by computer simulation methods. For this purpose, the free surface of aqueous solutions of fictitious salts is simulated at finite concentration both with nonpolarizable point-charge and polarizable Gaussian-charge potential models. The salts consist of monovalent cations and anions that are, apart from the sign of their charge, identical to each other. In particular, we consider the small Na⁺ and the large I^- ions together with their charge-inverted counterparts. In an attempt to avoid the interference even between the behavior of cations and anions, we also simulate systems containing only one of the above ions, and determine the free energy profile of these ions across the liquid-vapor interface of water at infinite dilution by potential of mean force (PMF) calculations. The obtained results reveal that, in the case of small ions, the anion is hydrated considerably stronger than the cation due to the close approach of water H atoms, bearing a positive fractional charge. As a consequence, the surface affinity of a small anion is even smaller than that of its cationic counterpart. However, considering that small ions are effectively repelled from the water surface, the importance of this difference is negligible. Further, a change in the hydration energy trends of the two oppositely charged ions is observed with their increasing size. This change is largely attributed to the fact that, with increasing ion size, the factor of 2 in the magnitude of the fractional charge of the closely approaching water atoms (i.e., O around cations and H around anions) outweighs the closer approach of the H than the O atom in the hydration energy. Thus, for large ions, being already surface active themselves, the surface affinity of the anion is larger than that of its positively charged counterpart. Further, such a difference is seen even in the case when the sign of the surface potential favors the adsorption of cations.

Surface Affinity of Tetramethylammonium Iodide in Aqueous Solutions: A Combined Experimental and Computer Simulation Study

The surface affinity of tetramethylammonium iodide (TMAI) in aqueous solutions is investigated by surface tension measurements and molecular dynamics computer simulations. Experiments, performed in the entire composition range of solubility using the pendant drop method with two different setups, clearly reveal that TMAI is a weakly capillary active salt. Computer simulations performed with the AMBER force field reproduce the experimental data very well, while two other major force fields (i.e., CHARMM and OPLS) can still reproduce the experimental trend qualitatively; however, even qualitative reproduction of the experimental trend requires scaling down the ion charges according to the Leontyev-Stuchebrukhov correction. On the other hand, the GROMOS force field fails in reproducing the experimentally confirmed capillary activity of TMAI. Molecular dynamics simulation results show that, among the two ions, iodide has a clearly larger surface affinity than tetramethylammonium (TMA⁺). Further, the adsorption of the I^- anions is strictly limited to the first molecular layer beneath the liquid-vapor interface, which is followed by several layers of their depletion. On the other hand, the net negative charge of the surface layer, caused by the excess amount of I^- with respect to TMA⁺, is compensated by a diffuse layer of adsorbed TMA⁺ cations, extending to or beyond the fourth molecular layer beneath the liquid surface.

Mechanism, Kinetics, and Potential of Mean Force of Evaporation of Water from Aqueous Sodium Chloride Solutions of Varying Concentrations

The mechanism, kinetics, and potential of mean force of evaporation of water from aqueous NaCl solutions are investigated through both unbiased molecular dynamics simulations and also biased simulations using the umbrella sampling method. The results are obtained for aqueous solutions of three different NaCl concentrations ranging from 0.6 to 6.0 m and also for pure water. The rate of evaporation is found to decrease in the presence of ions. It is found that the process of evaporation of a surface water molecule from ionic solutions can be triggered through its collision with another water or chloride ion. Such collisions provide the additional kinetic energy that is required for evaporation. However, when the collision takes place with a Cl^- ion, the evaporation of the escaping water also involves a collision with water in the vicinity of the ion at the same time along with the ion-water collision. These two collisions together provide the required kinetic energy for escape of the evaporating water molecule. Thus, the mechanism of evaporation process of ionic solutions can be more complex than that of pure water. The potential of mean force (PMF) of evaporation is found to be positive and it increases with increasing ion concentration. Also, no barrier in the PMF is found to be present for the condensation of water from vapor phase to the surfaces of the solutions. A detailed analysis of the unsuccessful evaporation attempts by surface water molecules is also made in the current study.

Chemical Feedback in the Self-Assembly and Function of Air-Liquid Interfaces: Insight into the Bottlenecks of CO2 Direct Air Capture

As fossil fuels remain a major source of energy throughout the world, developing efficient negative emission technologies, such as direct air capture (DAC), which remove carbon dioxide (CO₂) from the air, becomes critical for mitigating climate change. Although all DAC processes involve CO₂ transport from air into a sorbent/solvent, through an air-solid or air-liquid interface, the fundamental roles the interfaces play in DAC remain poorly understood. Herein, we study the interfacial behavior of amino acid (AA) solvents used in DAC through a combination of vibrational sum frequency generation spectroscopy and molecular dynamics simulations. This study revealed that the absorption of atmospheric CO₂ has antagonistic effects on subsequent capture events that are driven by changes in bulk pH and specific ion effects that feedback on surface organization and interactions. Among the three AAs (leucine, valine, and phenylalanine) studied, we identify and separate behaviors from CO₂ loading, chemical changes, variations in pH, and specific ion effects that tune structural and chemical degrees of freedom at the air-aqueous interface. The fundamental mechanistic findings described here are anticipated to enable new approaches to DAC based on exploiting interfaces as a tool to address climate change.

Orcinol and resorcinol induce local ordering of water molecules near the liquid-vapor interface

Resorcinol and orcinol are simple members of the family of phenolic compounds present in particulate matter in the atmosphere; they are amphiphilic in nature and thus surface active in aqueous solution. Here, we used X-ray photoelectron spectroscopy to probe the concentration of resorcinol (benzene-1,3-diol) and orcinol (5-methylbenzene-1,3-diol) at the liquid-vapor interface of aqueous solutions. Qualitatively consistent surface propensity and preferential orientation was obtained by molecular dynamics simulations. Auger electron yield near-edge X-ray absorption fine structure (NEXAFS) spectroscopy was used to probe the hydrogen bonding (HB) structure, indicating that the local structure of water molecules near the surface of the resorcinol and orcinol solutions tends towards a larger fraction of tetrahedrally coordinated molecules than observed at the liquid-vapor interface of pure water. The order parameter obtained from the molecular dynamics simulations confirm these observations. This effect is being discussed in terms of the formation of an ordered structure of these molecules at the surface leading to patterns of hydrated OH groups with distances among them that are relatively close to those in ice. These results suggest that the self-assembly of phenolic species at the aqueous solution-air interface could induce freezing similar to the case of fatty alcohol monolayers and, thus, be of relevance for ice nucleation in the atmosphere. We also attempted at looking at the changes of the O 1b₁, 3a₂ and 1b₂ molecular orbitals of liquid water, which are known to be sensitive to the HB structure as well, in response to the presence of resorcinol and orcinol. However, these changes remained negligible within uncertainty for both experimentally obtained valence spectra and theoretically calculated density of states.

Liquid-Liquid Extraction of the Eu(III) Cation by BTP Ligands into Ionic Liquids: Interfacial Features and Extraction Mechanisms Investigated by MD Simulations

What happens at the ionic-liquid (IL)/water interface when the Eu³⁺ cation is complexed and extracted by bis(dimethyltriazinyl) pyridine "BTP" ligands has been investigated by molecular dynamics and potential of mean force simulations on the interface crossing by key species: neutral BTP, its protonated BTPH⁺ form, Eu³⁺, and the Eu(BTP)₃³⁺ complex. At both the [BMI][Tf₂N]/water and [OMI][Tf₂N]/water interfaces, neither BTP nor Eu(BTP)₃³⁺ are found to adsorb. The distribution of Eu(BTP)₂³⁺ and Eu(BTP)³⁺ precursors of Eu(BTP)₃³⁺, and of their nitrate adducts, implies the occurrence of a stepwise complexation process in the interfacial domain, however. The analysis of the ionic content of the bulk phases and of their interface before and after extraction highlights the role of charge buffering by interfacial IL cations and anions, by different amounts depending on the IL. Comparison of ILs with octanol as the oil phase reveals striking differences regarding the extraction efficiency, the affinity of Eu(BTP)₃³⁺ for the interface, the effects of added nitric acid and of counterions (NO₃^- vs Tf₂N^-), charge neutralization mechanisms, and the extent of "oil" heterogeneity. Extraction into octanol is suggested to proceed via adsorption at the surface of water pools, nanoemulsions, or droplets, with marked counterion effects.

Computer Simulation of the Surface of Aqueous Ionic and Surfactant Solutions

The surface of aqueous solutions of simple salts was not the main focus of scientific attention for a long while. Considerable interest in studying such systems has only emerged in the past two decades, following the pioneering finding that large halide ions, such as I^-, exhibit considerable surface affinity. Since then, a number of issues have been clarified; however, there are still several unresolved points (e.g., the effect of various salts on lateral water diffusion at the surface) in this respect. Computer simulation studies of the field have largely benefited from the appearance of intrinsic surface analysis methods, by which the particles staying right at the boundary of the two phases can be unambiguously identified. Considering complex ions instead of simple ones opens a number of interesting questions, both from the theoretical point of view and from that of the applications. Besides reviewing the state-of-the-art of intrinsic surface analysis methods as well as the most important advances and open questions concerning the surface of simple ionic solutions, we focus on two such systems in this Perspective, namely, the surface of aqueous mixtures of room temperature ionic liquids and that of ionic surfactants. In the case of the former systems, for which computer simulation studies have still scarcely been reported, we summarize the theoretical advances that could trigger such investigations, which might well be of importance also from the point of view of industrial applications. Computer simulation methods are, on the other hand, widely used in studies of the surface of surfactant solutions. Here we review the most important theoretical advances and issues to be addressed and discuss two areas of applications, namely, the inclusion of information gathered from such simulations in large scale atmospheric models and the better understanding of the airborne transmission of viruses, such as SARS-CoV-2.

Expulsion of Hydroxide Ions from Methyl Hydration Shells

The affinity of hydroxide ions for methyl hydration shells is assessed using a combined experimental and theoretical analysis of tert-butyl alcohol (TBA) dissolved in pure water and aqueous NaOH and NaI. The experimental results are obtained using Raman multivariate curve resolution (Raman-MCR) and a new three-component total least squares (Raman-TLS) spectral decomposition strategy used to highlight vibrational perturbations resulting from interactions between TBA and aqueous ions. The experiments are interpreted and extended with the aid of effective fragment potential molecular dynamics (EFP-MD) simulations, as well as Kirkwood-Buff calculations and octanol/water partition measurements, to relate TBA-ion distribution functions to TBA solubility changes. The combined experimental and simulation results reveal that methyl group hydration shells more strongly expel hydroxide than iodide anions, whose populations near the methyl groups of TBA are predicted to be correlated with sodium counterion localization near the TBA hydroxyl group.

Valorisation of agricultural waste derived biochars in aquaculture to remove organic micropollutants from water - experimental study and molecular dynamics simulations

In this work, we evaluated the valorisation of agricultural waste materials by transforming coconut husks and shells, corncobs and rice straw into biochar for water treatment in aquaculture. We compared the biochars' suitability for removal of organic micropollutants (acetaminophen, oxytetracycline, tetracycline, enrofloxacin, atrazine, diuron and diclofenac) from surface water needed for aquaculture. The biochars were prepared by three methods ranging from inexpensive drum kilns (200 °C) to pyrolysis with biogasfication (350-750 °C). Overall, antibiotics tetracycline and enrofloxacin were the most strongly sorbed micropollutants, and coconut husk biochar prepared at 750 °C was the best sorbent material. Molecular Dynamics simulations indicated that the major sorption mechanism is via π-π stacking interactions and there is a possibility of multilayer sorption for some of the micropollutants. We observed, a strong impact of ionic strength (salinity), which is an important consideration in coastal aquaculture applications. High salinity decreased the sorption for antibiotics oxytetracycline, tetracycline and enrofloxacin but increased diclofenac, atrazine and diuron sorption. We considered coconut husk biochar produced in drum kilns the most practical option for biochar applications in small-scale coastal aquacultures in South Asia. Pilot trials of canal water filtration at an aquaculture farm revealed that micropollutant sorption by coconut husk biochar under real-world conditions might be 10-500 times less than observed in the laboratory studies. Even so, biochar amendment of sand enhanced the micropollutant retention, which may facilitate subsequent biodegradation and improve the quality of brackish surface water used for food production in coastal aquaculture.

Effects of Salts on the Solvation of Hydrophobic Objects in Water

The solvation of large, hydrophobic objects in water is facilitated by the formation of a low-density region surrounding the solute that is separated from the bulk liquid by an interface, which has a structure that resembles that between a liquid and its vapor. We study the effect of dissolved sodium chloride on the thermodynamics of solvation and on the solvent structure surrounding hydrophobic solutes in the size regime where this interface is not yet fully formed. Using biased Molecular Dynamics computer simulations, we calculate solvation free energies and orientational distributions of water molecules at different salt concentrations and solute sizes. We find that while the effects of sodium chloride on thermodynamic properties are small, the ions' response to the presence of a hydrophobic solute differs significantly from that of the water. Our findings provide mechanistic insight into how our understanding of hydrophobic solvation in water can be extended to electrolyte solutions.

Contribution of Different Molecules and Moieties to the Surface Tension in Aqueous Surfactant Solutions. II: Role of the Size and Charge Sign of the Counterions

Understanding the role of the counterion species in surfactant solutions is a complicated task, made harder by the fact that, experimentally, it is not possible to vary independently bulk and surface quantities. Here, we perform molecular dynamics simulations at constant surface coverage of the liquid/vapor interface of lithium, sodium, potassium, rubidium, and cesium dodecyl sulfate aqueous solutions. We investigate the effect of counterion type and charge sign on the surface tension of the solution, analyzing the contribution of different species and moieties to the lateral pressure profile. The observed trends are qualitatively compatible with the Hofmeister series, with the notable exception of sodium. We point out a possible shortcoming of what is at the moment, in our experience, the most realistic nonpolarizable force field (CHARMM36) that includes the parametrization for the whole series of alkali counterions. In the artificial system where the counterion and surfactant charges are inverted in sign, the counterions become considerably harder. This charge inversion changes considerably the surface tension contributions of the counterions, surfactant headgroups, and water molecules, stressing the key role of the hardness of the counterions in this respect. However, the hydration free energy gain of the counterions, occurring upon charge inversion, is compensated by the concomitant free energy loss of the headgroups and water molecules, leading to a negligible change in the surface tension of the entire system.

Insight into the Ionizing Surface Potential Method and Aqueous Sodium Halide Surfaces

Complementing the microscopic picture of the surface structure of electrolyte solutions set out by previous theoretical and experimental studies, the ionizing surface potential technique offers a unique approach to quantifying the impact of aqueous inorganic ions upon the interfacial electric field of the air-aqueous interface. In this Feature Article, we review the vulnerability of theoretical and empirically derived χ_water values as a normative reference for aqueous ion surface potentials. Instead, we recognize and evaluate aqueous ion surface potentials relative to well-known ionic surfactants cetyltrimethylammonium bromide (CTAB) and sodium dodecyl sulfate (SDS). Additionally, we also explore factors that impact the magnitude of the measured surface potentials using the ionizing method, particularly in the type of reference electrode and ionizing gas environment. With potential measurements of sodium halide solutions, we show that iodide has a dominant effect on the air-aqueous electric field. Compared to chloride and bromide, iodide is directly observed with a net negatively charged surface electric field at all salt concentrations measured (0.2 to 3.0 mol/kg water). Also, above the 2 M region, bromide is observed with a net negatively charged surface. Although several scenarios contribute to this effect, it is most likely due to the surface enrichment of bromide and iodide. While the results of this study are pertinent to determining the specific interfacial reactivity of aqueous halides, these anions seldom transpire as single-halide systems in the natural environment. Therefore, we also provide an outlook on future research concerning surface potential methods and more complex aqueous electrolyte systems.

Relationship between liquid nanoscale structure in solvents and the strength of the Hofmeister effect in extraction experiments

In this study, we used mixtures of carboxylic acids and amines as solvents for the liquid-liquid extraction of copper salts with various anions from aqueous phase, and systematically varied the acid/amine ratio to determine its influence on extraction efficiency. The organic phases resulting from these extraction experiments were studied using small-angle X-ray scattering (SAXS), establishing a connection between the extraction process and the liquid structure. A relationship is found between the extent of extraction for the metal salt, the strength of the Hofmeister effect of the anions of the salt, and the characteristic lengthscale of the observed liquid nanoscale structure before and after extraction.

Single-Particle Dynamics at the Intrinsic Surface of Aqueous Alkali Halide Solutions

The distribution of ions in the proximity of the liquid-vapor interface of their aqueous solution has been the subject of an intense debate during the last decade. The effects of ionic polarizability have been one of its salient aspects. Much less has been said about the corresponding dynamical properties, which are substantially unexplored. Here, we investigate the single-particle dynamics at the liquid-vapor interface of several alkali halide solutions, using molecular dynamics simulations with polarizable and nonpolarizable force fields and intrinsic surface analysis. We analyze the diffusion coefficient, residence time, and velocity autocorrelation function of water and ions and investigate how these properties depend on the molecular layer where they reside. While anions are found in the first molecular layer for relatively long times, cations are only making quick excursions into it, thanks to thermal fluctuations. The in-layer residence time of ions and their molar fraction in the layer turned out to be linearly dependent on each other. We interpret this unexpected result using a simple two-state model. In addition, we found that, unlike water and other neat molecular liquids that show a different diffusion mechanism at the surface than in the bulk of their liquid phase, ions do not enjoy enhanced mobility in the surface layer of their aqueous solution. This result indicates that ions in the surface layer are shielded by their nearest water neighbors from being exposed to the vapor phase as much as possible. Such positions are available for the ions at the negatively curved troughs of the molecularly rugged liquid surface.

Interfacial and bulk properties of concentrated solutions of ammonium nitrate

We conducted molecular dynamics (MD) simulations to calculate the density and surface tension of concentrated ammonium nitrate (AN) solutions up to the solubility limit of ammonium nitrate in water, by combining the SPC/E, SPCE/F and TIP4P/2005 water models with OPLS model for ammonium and nitrate ions. This is the first time that the properties of concentrated solutions of nitrates, especially AN, have been studied by molecular dynamics. We effectively account for the polarisation effects by the electronic continuum correction (ECC), practically realised via rescaling of the ionic charges. We found that, the full-charge force field MD simulations overestimate the experimental results, as the ions experience repulsion from the interface and prefer to remain in the subsurface layer and the bulk solution. In contrast, reducing the ionic charges results in the behaviour that fits well with the experimental data. The nitrate anions display a greater propensity for the interface than the ammonium cations. We accurately predict both the density and the rise in the surface tension of concentrated solutions of AN, recommending TIP4P/2005 for water and the scaled-charge OPLS model (OPLS/ECC) for the ions in the solutions. We observe that, the adsorption of anions to the interface accompanies their depletion in the subsurface layer, which is preferentially occupied by cations, resulting in the formation of the electric double layer. We demonstrate the ion deficiency for up to 3 Å below the surface and establish the requirement to include the polarisability effects in the OPLS model for AN. While these results confirmed the findings of the previous studies for dilute solutions, they are new in the solubility limit. Concentrated solutions exhibit a strong effect of the abundance of solute on the coordination numbers of ions and on the degree of ion pairing. Surprisingly, ion pairing decreases significantly at the interface compared with the bulk. The present study identifies OPLS/ECC, along with TIP4P/2005, to yield accurate predictions of physical properties of concentrated AN, with precision required for industrial applications, such as a formulation of emulsion and fuel-oil explosives that now predominate the civilian use of AN. An application of this model will allow one to predict the surface properties of supersaturated solutions of AN which fall outside the capability of the present laboratory experiments but are important industrially.

Surface Affinity of Alkali and Halide Ions in Their Aqueous Solution: Insight from Intrinsic Density Analysis

The surface tension of all aqueous alkali halide solutions is higher than that of pure water. According to the Gibbs adsorption equation, this indicates a net depletion of these ions in the interfacial region. However, simulations and experiments show that large, soft ions, such as I^-, can accumulate at the liquid/vapor interface. The presence of a loose hydration shell is usually considered to be the reason for this behavior. In this work, we perform computer simulations to characterize the liquid-vapor interface of aqueous alkali chloride and sodium halide solutions systematically, considering all ions from Li⁺ to Cs⁺ and from F^- to I^-. Using computational methods for the removal of surface fluctuations, we analyze the structure of the interface at a dramatically enhanced resolution, showing that the positive excess originates in the very first molecular layer and that the next 3-4 layers account for the net negative excess. With the help of a fictitious system with charge-inverted ion pairs, we also show that it is not possible to rationalize the surface affinity of ions in solutions in terms of the properties of anions and cations separately. Moreover, the surface excess is generally dominated by the smaller of the two ions.

Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling

Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.

Specific counter-cation effect on the molecular orientation of thiocyanate anions at the aqueous solution interface

Understanding the interfacial structure of aqueous electrolyte solutions is important and relevant to a wide range of systems, ranging from atmospheric aerosols to electrochemistry, and biological environments. Though significant efforts have been made to unravel the interfacial structure of water molecules, the structure and dynamics of ions at the interface have not yet been fully elucidated. Here, the interfacial structure of the aqueous solution was investigated directly by monitoring the thiocyanate (SCN-) anions using surface-specific sum frequency generation (SFG) vibrational spectroscopy. The molecular orientation of the SCN- anions and their adsorption behavior at the air/water interface were systematically determined by quantitative polarization analysis. The transition dipole of the CN stretching of the SCN- anion is oriented around 44° from the surface normal of the NaSCN aqueous solution surface and remained unchanged with the bulk concentration varying from 1 mol kg-1 to 13 mol kg-1. The free energy of adsorption of SCN- anions at the air/water interface was determined to be -1.53 ± 0.04 kcal mol-1. Furthermore, a new SFG peak positioned at 2080 cm-1 in the ppp polarization combination was observed at the air/15.0 mol kg-1 NaSCN aqueous solution interface for the first time. Concentration-dependent SFG analysis and density functional theory (DFT) calculation further revealed that the SCN- anions form an ion clustering structure at the air/water interface. The subtle and specific Na+ and K+ counter-cation effects on the interfacial structure of the SCN- anions at the aqueous solution interface were also observed, which showed that ion cooperativity plays an important role in affecting the interfacial structure of ions at the air/water interface. The results are expected to yield significant insights into the understanding of the structure of aqueous solution surfaces and the molecular level mechanism of the cationic Hofmeister effect.

Surface Propensity of Aqueous Atmospheric Bromine at the Liquid-Gas Interface

Multiphase reactions of halide ions in aqueous solutions exposed to the atmosphere initiate the formation of molecular halogen compounds in the gas phase. Their photolysis leads to halogen atoms, which are catalytic sinks for ozone, making these processes relevant for the regional and global tropospheric ozone budget. The affinity of halide ions in aqueous solution for the liquid-gas interface, which may influence their reactivity with gaseous species, has been debated. Our study focuses on the surface properties of the bromide ion and its oxidation products. In situ X-ray photoelectron spectroscopy carried out on a liquid jet combined with classical and first-principles molecular dynamics calculations was used to investigate the interfacial depth profile of bromide, hypobromite, hypobromous acid, and bromate. The simulated core electron binding energies support the experimentally observed values, which follow a correlation with bromine oxidation state for the anion series. Bromide ions are homogeneously distributed in the solution. Hypobromous acid, a key species in the multiphase cycling of bromine, is the only species showing surface propensity, which suggests a more important role of the interface in multiphase bromine chemistry than thought so far.

Gating of Hydrophobic Nanopores with Large Anions

Understanding ion transport in nanoporous materials is critical to a wide variety of energy and environmental technologies, ranging from ion-selective membranes, drug delivery, and biosensing, to ion batteries and supercapacitors. While nanoscale transport is often described by continuum models that rely on a point charge description for ions and a homogeneous dielectric medium for the solvent, here, we show that transport of aqueous solutions at a hydrophobic interface can be highly dependent on the size and hydration strength of the solvated ions. Specifically, measurements of ion current through single silicon nitride nanopores that contain a hydrophobic-hydrophilic junction show that transport properties are dependent not only on applied voltage but also on the type of anion. We find that in Cl^--containing solutions the nanopores only conducted ionic current above a negative voltage threshold. On the other hand, introduction of large polarizable anions, such as Br^- and I^-, facilitated the pore wetting, making the pore conductive at all examined voltages. Molecular dynamics simulations revealed that the large anions, Br^- and I^-, have a weaker solvation shell compared to that of Cl^- and consequently were prone to migrate from the aqueous solution to the hydrophobic surface, leading to the anion accumulation responsible for pore wetting. The results are essential for designing nanoporous systems that are selective to ions of the same charge, for realization of ion-induced wetting in hydrophobic pores, as well as for a fundamental understanding on the role of ion hydration shell on the properties of solid/liquid interfaces.

Role of the Counterions in the Surface Tension of Aqueous Surfactant Solutions. A Computer Simulation Study of Alkali Dodecyl Sulfate Systems

We have investigated the surface tension contributions of the counterions, surfactant headgroups and tails, and water molecules in aqueous alkali dodecyl sulfate (DS) solutions close to the saturated surface concentration by analyzing the lateral pressure profile contribution of these components using molecular dynamics simulations. For this purpose, we have used the combination of two popular force fields, namely KBFF for the counterions and GROMOS96 for the surfactant, which are both parameterized for the SPC/E water model. Except for the system containing Na+ counterions, the surface tension of the surfactant solutions has turned out to be larger rather than smaller than that of neat water, showing a severe shortcoming of the combination of the two force fields. We have traced back this failure of the potential model combination to the unphysically strong attraction of the KBFF counterions, except for Na+, to the anionic head of the surfactants. Despite this failure of the model, we have observed a clear relation between the soft/hard character (in the sense of the Hofmeister series) and the surface tension contribution of the counterions, which, given the above limitations of the model, can only be regarded as an indicative result. We emphasize that the obtained results, although in a twisted way, clearly stress the crucial role the counterions of ionic surfactants play in determining the surface tension of the aqueous surfactant solutions.

Adsorption of ionic surfactants at the air-water interface: The gap between theory and experiment

We review the experimental and theoretical results for the adsorption and structure of ionic surfactants at the air-liquid interface. The results show that ionic surfactants form thick adsorption layers at the interfacial region. We also review several adsorption models for ionic surfactants, which become increasingly complex as they capture the many features of adsorption layers. However, the adsorption layer structures determined by experiments and the structures predicted by models do not match because most models assume very thin adsorption layers. We show the discrepancies between measured and predicted surface properties and provide several explanations. We conclude that the mismatch in the adsorption layer structure provided by experiments and the structure provided by adsorption models is the main reason for the discrepancies in the surface excess and the surface potential.

Do soft anions promote protein denaturation through binding interactions? A case study using ribonuclease A

It has long been known that large soft anions like bromide, iodide and thiocyanate are protein denaturing agents, but their mechanism of action is still unclear. In this work we have investigated the protein denaturing properties of these anions using Ribonuclease A (RNase A) as a model protein system. Salt-induced perturbations to the protein folding free energy were determined using differential scanning calorimetry and the results demonstrate that the addition of sodium iodide and sodium thiocyanate significantly decreases the melting temperature of the protein. In order to account for this reduction in protein stability, we show that the introduction of salts that contain soft anions to the aqueous solvent perturbs the protein unfolding free energy through three mechanisms: (a) screening Coulomb interactions that exist between charged protein residues, (b) Hofmeister effects, and (c) specific anion binding to CH and CH2 moieties in the protein polypeptide backbone. Using the micellization of 1,2-hexanediol as a ruler for hydrophobicity, we have devised a practical methodology that separates the Coulomb and Hofmeister contributions of salts to the protein unfolding free energy. This allowing us to isolate the contribution of soft anion binding interactions to the unfolding process. The analysis shows that binding contributions have the largest magnitude, confirming that it is the binding of soft anions to the polypeptide backbone that is the main promoter of protein unfolding.

Permeating disciplines: Overcoming barriers between molecular simulations and classical structure-function approaches in biological ion transport

Ion translocation across biological barriers is a fundamental requirement for life. In many cases, controlling this process-for example with neuroactive drugs-demands an understanding of rapid and reversible structural changes in membrane-embedded proteins, including ion channels and transporters. Classical approaches to electrophysiology and structural biology have provided valuable insights into several such proteins over macroscopic, often discontinuous scales of space and time. Integrating these observations into meaningful mechanistic models now relies increasingly on computational methods, particularly molecular dynamics simulations, while surfacing important challenges in data management and conceptual alignment. Here, we seek to provide contemporary context, concrete examples, and a look to the future for bridging disciplinary gaps in biological ion transport. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

A comparison of sodium and hydrogen halides at the air-water interface

New molecular models, parameterized to ab initio calculations, were developed to describe HBr and HI at the air-water interface. These were used to compare how the air-water interface influenced dissociation of NaX and HX, with X being Cl, Br, or I, and also their propensity for the interface. The polarizable multistate empirical valence bond method, which explicitly describes proton sharing, was used to model HX. Results showed that the air-water interface suppressed HX dissociation from a contact ion pair to a solvent separated to a greater degree than NaX dissociation. Furthermore, HX had a greater propensity for the interface than NaX, which was a consequence of the hydronium ion having a greatest interfacial activity of all species studied. As a consequence of this, the average configuration of dissociated HX, while in both contact ion and solvent separated ion pairs near the air-water interface, is with the dissociated hydrogen oriented more towards the air than the X atom.

Salt Promotes Protonation of Amine Groups at Air/Water Interface

Interfacial water reorientation caused by charged Langmuir monolayers consisting of primary fatty amine (ODA) and cationic lipid having quaternary amine headgroup (DPTAP) were investigated by interface-selective vibrational sum-frequency generation spectroscopy. For DPTAP monolayer, initially large sum-frequency intensity from interfacial water OH band decreased steadily by increasing monovalent salt (NaCl, NaI) concentration due to counterion adsorption. On the other hand, ODA/water exhibited significantly smaller sum-frequency intensity than DPTAP/water, implying only small portion of protonated amine group (-NH₃⁺) initially existed. By increasing the ionic strength, however, SF intensity of water OH band was enhanced markedly up to ∼1 mM, and then decreased in both NaCl and NaI solutions. By measuring the phase of the sum-frequency spectra, it was found that water dipoles under the ODA headgroup point downward, indicating that the surfaces were always positively charged. This demonstrated that increasing ionic strength facilitates protonation of primary amine headgroups. A simple model based on Poisson-Boltzmann (PB) theory explained this protonation behavior of primary amines.

Ambient pressure photoelectron spectroscopy: Practical considerations and experimental frontiers

Over the past several decades, ambient pressure x-ray photoelectron spectroscopy (APXPS) has emerged as a powerful technique for in situ and operando investigations of chemical reactions under relevant ambient atmospheres far from ultra-high vacuum conditions. This review focuses on exemplary cases of APXPS experiments, giving special consideration to experimental techniques, challenges, and limitations specific to distinct condensed matter interfaces. We discuss APXPS experiments on solid/vapor interfaces, including the special case of 2D films of graphene and hexagonal boron nitride on metal substrates with intercalated gas molecules, liquid/vapor interfaces, and liquid/solid interfaces, which are a relatively new class of interfaces being probed by APXPS. We also provide a critical evaluation of the persistent limitations and challenges of APXPS, as well as the current experimental frontiers.

Observation of isolated ionic liquid cations and water molecules in an inert solvent

1-Octyl-3-methyl imidazolium halides ([OMIM]I and [OMIM]Cl) were loaded on top of CCl₄, and an in situ inclusion process was monitored from the CCl₄ phase as time elapses by infrared absorption spectroscopy. Absorption from the bands corresponding to the C(2)-H and C(4,5)-H stretch modes in the imidazolium cation was reduced significantly compared to the bulk IL spectra. This indicates that (1) the [OMIM] cation exists in CCl₄ as a monomer, dissociated from the anion and other cations, and (2) hydrogen bonding between the anion and the cation increases the dipole strength of the CH moieties in the imidazolium ring. [OMIM]I was found to transfer into the CCl₄ phase much faster than [OMIM]Cl, and this instigated us to compare the transfer of aqueous solutions of ionic liquids, 1-butyl-3-methyl imidazolium halides ([BMIM]I and [BMIM]Cl) into the CCl₄ matrix. Not only [BMIM]I but also water molecules transferred faster compared to those in the [BMIM]Cl aqueous solution. Water molecules in CCl₄ were shown to form clusters in [BMIM]I; presumably, I^- anions work as nucleation centers of water clusters.

Liquid–liquid extraction of alkali cations by 18-crown-6: complexation and interface crossing studied by MD and PMF simulations

The 18C6/M⁺Pic⁻complexes form and adsorb “right at the nano-interface” where 18C6 prefers the K⁺guest.

Surface evolution of manganese chloride aqueous droplets resulting in self-suppressed evaporation

The exchange kinetics of liquid water, which are of fundamental interest and have potential applications, remain unclear. A fantastic and extraordinary phenomenon was observed during the evaporation of a water droplet doped with manganese chloride. As observed from the evolution of this type of droplet, a thin film was formed on the surface with an exothermic phase transition, resulting in self-suppressed evaporation. The MnCl₂-doped water droplets were maintained in a relative humidity (RH) of 50% at 40 °C for more than a week and for longer than two months at a temperature of 25 °C. In contrast, a pure water droplet can only be sustained for a few minutes. The self-suppressed evaporation of doped water may be due to the special hydration of the accumulated manganese and chloride ions at the surface, decreasing the surface tension.

β-Cyclodextrin at the Water/1-Bromobutane Interface: Molecular Insight into Reverse Phase Transfer Catalysis

Molecular insight into the role of β-cyclodextrin (βCD) as a phase transfer catalyst at the liquid/liquid interface is obtained by molecular dynamics simulations of the structure and dynamics of βCD adsorbed at the interface between water and 1-bromobutane. In particular, we consider the structure and dynamics of the water and bromobutane molecules inside the βCD cavity and compare them with the behavior when βCD is dissolved in bulk water. βCD is preferentially oriented at the interface, with the cavity opening along the interface normal. While in bulk water the cavity includes 6-8 water molecules that are relatively mobile with short residence time, at the interface the cavity is mostly dehydrated and includes a single bromobutane molecule. This inclusion complex is stable in bulk water. The implication of this behavior for reverse phase transfer catalysis is discussed.