Assessing long-range contributions to the charge asymmetry of ion adsorption at the air–water interface

Anions generally associate more favorably with the air–water interface than cations. In addition to solute size and polarizability, the intrinsic structure of the unperturbed interface has been discussed as an important contributor to this bias. Here we assess quantitatively the role that intrinsic charge asymmetry of water''s surface plays in ion adsorption, using computer simulations to compare model solutes of various size and charge. In doing so, we also evaluate the degree to which linear response theory for solvent polarization is a reasonable approach for comparing the thermodynamics of bulk and interfacial ion solvation. Consistent with previous works on bulk ion solvation, we find that the average electrostatic potential at the center of a neutral, sub-nanometer solute at the air–water interface depends sensitively on its radius, and that this potential changes quite nonlinearly as the solute''s charge is introduced. The nonlinear response closely resembles that of the bulk. As a result, the net nonlinearity of ion adsorption is weaker than in bulk, but still substantial, comparable to the apparent magnitude of macroscopically nonlocal contributions from the undisturbed interface. For the simple-point-charge model of water we study, these results argue distinctly against rationalizing ion adsorption in terms of surface potentials inherent to molecular structure of the liquid''s boundary.

Cations and anions have different affinities for the air-water interface. The intrinsic orientation of surface molecules suggests such an asymmetry, but the bias is dominated by solvent response that is spatially local and significantly nonlinear.

Counter to expectations from conventional theories of solvation, there is a large body of both computational and experimental evidence indicating that small ions can adsorb to the air–water interface.^1–9 Implications across the biological, atmospheric and physical sciences have inspired efforts to understand the microscopic driving forces for ions associating with hydrophobic interfaces in general.^10–21 A particular emphasis has been placed on understanding ion specificity, i.e., why some ions exhibit strong interfacial affinity while others do not. Empirical trends indicate that ion size and polarizability are important factors, as could be anticipated from conventional theory. More surprisingly, the sign of a solute''s charge can effect a significant bias, with anions tending to adsorb more favorably than cations.Here we examine the microscopic origin of this charge asymmetry in interfacial ion adsorption. We specifically assess whether the thermodynamic preference can be simply and generally understood in terms of long-range biases that are intrinsic to an aqueous system surrounded by vapor. By “long-range” and “nonlocal” we refer to macroscopically large scales, i.e., collective forces that are felt at arbitrarily long distance. Such a macroscopically long-range bias is expected from the air–water interface due to its average polarization, and by some measures the bias is quite strong. By contrast, “local” contributions comprise the entire influence of a solute''s microscopic environment, including electrostatic forces from molecules that are many solvation shells away – any influence that decays over a sub-macroscopic length scale.The importance of macroscopically nonlocal contributions has been discussed extensively in the context of ion solvation in bulk liquid water, which we review in Section 1 as a backdrop for interfacial solvation. The notion that such contributions strongly influence charge asymmetry of solvation at the air–water interface has informed theoretical approaches and inspired criticism of widely used force fields for molecular simulation.^22,23 A full understanding of their role in interfacial adsorption, however, is lacking.In the course of this study, we will also evaluate the suitability of dielectric continuum theory (DCT) to describe the adsorption process. DCT has provided an essential conceptual framework for rationalizing water''s response to electrostatic perturbations. But a more precise understanding of its applicability is needed, particularly for the construction of more elaborate models (e.g., with heterogeneous polarizability near interfaces^24–26) and for the application of DCT to evermore complex (e.g., nanoconfined^27,28) environments.