The Electrochemical Surface Potential Due to Classical Point Charge Models Drives Anion Adsorption to the Air-Water Interface

The Electrochemical Surface Potential Due to Classical Point Charge Models Drives Anion Adsorption to the Air-Water Interface

Marcel D. Baer Chemical and Materials Science Division, Pacific Northwest National Laboratory, Richland    Abraham C. Stern [    Yan Levin Instituto de Física, UFRGS Porto Alegre, RS, 91501-970, Brazil    Douglas J. Tobias [    Christopher J. Mundy Chemical and Materials Science Division, Pacific Northwest National Laboratory, Richland
July 12, 2019
July 12, 2019

We demonstrate that the driving forces for ion adsorption to the air-water interface for point charge models results from both cavitation and a term that is of the form of a negative electrochemical surface potential. We carefully characterize the role of the free energy due to the electrochemical surface potential computed from simple empirical models and its role in ionic adsorption within the context of dielectric continuum theory. Our research suggests that the electrochemical surface potential due to point charge models provides anions with a significant driving force to the air-water interface. This is contrary to the results of ab initio simulations that indicate that the average electrostatic surface potential should favor the desorption of anions at the air-water interface. The results have profound implications for the studies of ionic distributions in the vicinity of hydrophobic surfaces and proteins.

University of California, Irvine] Department of Chemistry, University of California, Irvine, Irvine, California 92697-2025, United States University of California, Irvine] Department of Chemistry, University of California, Irvine, Irvine, California 92697-2025, United States {scheme}

Keywords: Statistical Mechanics, Dielectric Continuum Theory, Interfaces, Free-energy, Ion solvation

Wolfgang Pauli used to say “God made the bulk; the surface was invented by the devil.” 1 Pauli’s frustration referred to solid surfaces, which are relatively simple compared to liquid interfaces. Besides the broken symmetry associated with the solid surface, liquid-liquid and liquid-air interfaces exhibit capillary fluctuations, making their study even more difficult. If one of the liquid phases contains a mixture — solvent plus solute — the level of complexity deepens. We can only speculate about what Pauli would say about such a confounding geometry.

The first clues to the behavior of ions near the air-water interface came from the measurements of excess surface tension of electrolyte solutions. 2 In particular, it was observed that halide salts increase the surface tension of the air-solution interface. For this to be the case, thermodynamics requires that ions must be depleted from the interfacial region. 3 On the other hand, electrochemical measurements of Frumkin showed that the electrolyte solutions possess an electrostatic potential gradient across their surface. 4 Frumkin’s work clearly indicated that cations and anions can behave very differently near the surface by demonstrating that anions are closer to the surface than cations. Frumkin’s work was also consistent with another more than 100 year old mystery: the Hofmeister series of electrolyte solution. 4 In 1888 Hofmeister observed that various monovalent electrolytes have very different effects on protein solutions. 5 While some electrolytes help solubilize proteins, denaturing them in the process, others lead to protein precipitation. The interaction of ions with surfaces, such as water-protein or air-water interfaces, has remained an outstanding puzzle of physical chemistry.

Indications to the mechanism of the ion-surface interaction started to appear when computational methods, hardware and software became sufficiently powerful to allow realistic microscopic simulations. The work of Berkowitz on small water clusters in the early 1990s 6 and of Jungwirth and Tobias (JT) on extended air-water interface a decade ago brought to the attention of the physical chemistry community a new phenomenon associated with aqueous electrolytes, 7 namely, that soft, polarizable anions can adsorb to the air-water interface.

Levin and co-workers have recently shown that it is possible to quantitatively explain the surface tension measurements and the surface adsorption of large polarizable anions in terms of a dielectric continuum theory that explicitly takes into account the finite size and a full polarization response of the ions comprising the electrolyte to the presence of the interface. 8, 9, 10 By using the accepted literature values of anionic radii and polarizabilities, the radius of the hydrated sodium counter ion was fit to quantitatively reproduce the surface tension data of the alkali halide and other oxyanion salts. The results of this polarizable anion dielectric continuum theory (PA-DCT) were qualitatively consistent with the original work of JT.

An important difference between PA-DCT theory and JT is that the predicted amount of adsorption for the most polarizable iodide anion is significantly less in the former. State-of-the-art direct surface measurements suggest that there is indeed adsorption of large polarizable anions to the air-water interface. 11, 12, 13 However, precise measurements on the degree of anion adsorption at low concentrations usually rely on using the Gibbs’ adsorption isotherm (GAI) that relates surface excess of anions to the electrolyte concentration dependence of the surface tension. 14 For the case of iodide, the large adsorption free energies predicted by JT would produce concentration dependent surface tensions that resemble surfactants rather than electrolytes.

Although quantitatively different, both JT and the PA-DCT approaches point to the fundamental role of ionic polarizability for ion adsorption. Recently, a potential of mean force (PMF) of iodide at the air-water interface was computed using molecular interactions based on the density functional theory (DFT) (a so called ab initio simulation that naturally takes into account the full polarization response of the anion to the air-water interface) and yielded results in almost quantitative agreement with the PA-DCT. 15

In a recent study, Horinek and co-workers have fit a non-polarizable point charge, soft-sphere interaction potentials for the alkali-halide salts in a fixed charge non-polarizable water model. The resulting PMF obtained from the simulation was fitted to reproduce the experimentally measured surface tensions as a function of electrolyte concentration. 16 Interestingly, the PMF of a single iodide anion from the ab initio calculation and the PA-DCT compare well to the work of Horinek and co-workers: all of the PMFs contain a shallow minimum that is less than , where is Boltzmann’s constant and K. This is to be contrasted to the PMFs of iodide obtained with a polarizable force field similar to those used in the original study of JT that predicts the surface adsorption of roughly 17, 15 The results of the fitted non-polarizable force fields of Horinek and co-workers 18, 19 provide additional validation of the depth of the PMF minimum at the surface that is necessary to be consistent with surface tension measurements. However, other questions regarding the role of polarization in the driving force of anions to interfaces remain unanswered.

The purpose of this paper is to show that the driving force for adsorption of large halide anions when using non-polarizable point charge models for both the anion and water is due to very different physics than the ab initio MD simulation and PA-DCT. We do not dispute the success of point charge models (e.g. SPC/E) for describing a range of both bulk and interfacial properties of neat water. 20, 21 However, clues to a possible model dependence for the driving force of anions to interfaces may lie in the surface potential of the neat air-water interface that could strongly influence ionic adsorption. 19 The surface potential of the air-water interface computed utilizing a classical point charge model, differs in both sign and magnitude to the ab initio result of 22, 23, 24 In this study we will make a distinction between the electrochemical surface potential — the potential of mean force that an ion experiences due to the presence of an interface — and the average electrostatic potential due to the presence of the interface, which is the surface potential.

The PA-DCT is based on the assumption that there is no significant free-energy contribution from the electrochemical surface potential across the neat air-water interface. This ansatz is corroborated by the good agreement of the computed electrochemical surface potentials with PA-DCT and the measurements of Frumkin on electrolyte solutions which are indeed small (on the order of ). 4 An important conclusion that arises from the comparison PA-DCT and ab initio simulations is that ions in the DFT approximation, for reasons that are not completely clear at this time, feel an electrochemical surface potential of the air-water interface that is nearly zero. 13, 25, 15

In this letter we will isolate the effect of the electrochemical surface potential on driving ions to interfaces by examining the solvation of a hard-sphere ion in both an SPC/E model of water and a Stockmayer fluid, which by symmetry has no surface potential. Through the choice of a hard-sphere solute we can cleanly isolate the volume dependent cavitational energy from the dielectric response of the charging process. Thus, we will demonstrate that the models presented herein provide a straightforward decomposition of the ion adsorption propensity into cavitation and self-energy. The tug-of-war between the cavitation and self-energy penalty usually favors the bulk ionic solvation for non-polarizable models of anions. However, we will show that for fixed point charge models of ions and water there is an additional contribution to the free energy due to the electrochemical surface potential of the air-water interface. From these results two important findings will emerge: first, it will be shown that the free energy due to the electrochemical surface potential that a charged hard sphere feels as it approaches the interface is not the average surface potential of the air-water interface as is normally computed 26, 27, 24, 19. Second, we will demonstrate that the principles of anion adsorption are model dependent. The fundamental role of the electrochemical surface potential as the driving force of ions to interfaces found in point charge models is inconsistent with the mechanism that arises when a quantum mechanical description of the charge density is utilized.

An extended interfacial system containing 215 solvent molecules and a single iodide anion was modeled in slab geometry within a supercell of 151575 Å. This choice of system size has been shown to produce a stable bulk liquid in the center of the slab and is able to quantitatively reproduce structural properties of larger system sizes that have been used elsewhere. 21

It is clear from our results (see the supporting information (SI)) that we have a stable liquid for the Stockmayer fluid that can be used to embed a charged hard sphere cavity. Moreover, the surface potentials for both the Stockmayer and SPC/E fluid can be computed using the simple formula first presented by Wilson et al28 The resulting surface potential for SPC/E water is in excellent agreement with all studies performed on this system to date yielding a value of roughly V and validating our choice of system size. As expected, due to symmetry, the Stockmayer fluid yields a value of approximately V. In order to make contact with a dielectric continuum theory we estimate the dielectric constant of both SPC/E water and the Stockmayer fluid using the same number of solvent molecules in a cubic box at the estimated bulk density. 29 The dielectric constant for the SPC/E (65) is in good agreement with the reported values for 216 water molecules by van der Spoel et al. (69), 29 as our estimated bulk density is a bit lower.

Thus, we are positioned to compare the statistical mechanics of ion adsorption of a dielectric medium with zero and finite values of the surface potential. The dielectric continuum theory (DCT) presented herein, is built by considering various contributions to the solvation free energy of a hard-sphere ion (which differs from the PA-DCT discussed above). We first note that to transfer an ion into water requires the creation of a cavity. For small cavities, less than 4 Å, the cavitational energy is proportional to the volume of the cavity. When the ion of charge moves across the interface, the cavitation energy decreases proportionally to the volume exposed to the vapor phase. The cavitation energy for an ion of radius located at distance from the Gibbs dividing surface (GDS) 30 (positive is towards water) is given by:


where is in units of and the prefactor is in units of Å8 When the ion is in the bulk water, its electrostatic field is screened by the surrounding water molecules so that the electrostatic self-energy is , where the Bjerrum length is . When the ion is at distance from the interface, besides its interaction with the water molecules, it also interacts with the induced charge on the dielectric interface. The electric field produced by the induced charge is exactly the same as the field produced by a point image charge located across the interface. Since the dielectric contrast between water and air is so large, the image charge has the same sign and magnitude as the ion, and is located at . The electrostatic self energy of an ion located at can then be easily calculated to be . When the ion crosses the interface and is at , it sees an image in water, which has the charge , and the self energy is therefore . Unfortunately, there is no simple way to calculate the self energy for an ion that has only partially penetrated the interface; the energy will depend on the precise model that is use to treat the interior of an ion and even simple models result in very complicated mathematical expressions. 31 Within the mixed dielectric model of Tamashiro and Constantino, 31 we can make an estimate of the electrostatic self-energy of an ion located at , , for the SPC/E water with . We can now interpolate between the different limits to obtain on the interval ,


Finally, we relate the surface potential of SPC/E water to a free energy difference across the air-water interface. If we approximate this free energy as a discontinuous drop across a point dipole layer, we obtain the following model for free energy change due to the surface potential as the ion moves across the GDS,


Here will be determined by the value of the electrochemical surface potential. Combining all these contributions, we obtain the total potential for the ion in the SPC/E water, .

For the Stockmayer liquid there is no free energy due to the electrochemical surface potential contribution to the solvation free energy. The cavitational energy is more difficult to define because the air-water interface of the Stockmayer liquid is very diffuse. However if we adjust the energy scale so that the ion is fully solvated at one ionic radius from the GDS, the cavitation energy can be written as


The electrostatic self-energy of an ion inside a Stockmayer liquid with is


and the total energy is the sum of these two contributions.

In order to compare the free energy of ion adsorption from the aforementioned classical MD simulations to the DCT outlined above, we have used restrained molecular dynamics runs where the histogram of the coordinate of the ion is reconstructed using the weighted histogram analysis method (WHAM). 32, 33 Specifically, 32 windows spanning the range from  Å and  Å relative to the GDS for the SPC/E and Stockmayer fluid, respectively. The ion position was restrained with a harmonic force constant of 42.82 kcal/(mol Å) and each window was simulated for at least 2 ns.

Figure 1: The potential of mean force of transferring a hard-sphere anion with varying charge across the liquid-air interface (negative bulk, positive air) by molecular simulations in a Stockmayer (b) and SPC/E (d) fluid compared to DCT results for ions with fixed charge with and , respectively. Each curve represents the charge state of the hard-sphere. The excellent qualitative agreement with DCT is apparent in both cases indicating that DCT theory presented here for hard charged sphere adsorption does contain the essential physics to model the adsorption seen in our molecular simulations. The Stockmayer fluid is modeled with an electrochemical surface potential of and the SPC/E water employs the value in  3. Without the electrochemical surface potential the cavity should have the lowest free energy of adsorption. For the SPC/E, it is clear this is not the case, and the subtle effects obtained from molecular simulation are within the error bars of the WHAM method 32, 33 (see Figure 3 in the SI). For completeness, we have presented the PMF for (a fully charged monovalent ion) for comparison with the DCT model in the SI where similar agreement is obtained. Thus, it is clear that we can model the statistical mechanics of the full range of monovalent charges for hard-sphere adsorption with the procedure presented in this letter.

The main result of our study is presented in 1. Here one can see excellent qualitative agreement with simulation and DCT for the case of a fractionally charged hard sphere in SPC/E water and the Stockmayer fluid. The comparison between the simulation and the DCT in 1 in different panels has been made to establish an important point. The results in 1 are not a result of a fit to one another, but are the results of two different modeling approaches. Due to its technical nature, a detailed explanation of the quantitative differences between the simulations and the DCT is given in the SI.

In the case of SPC/E water interacting with a hard sphere in the limit of zero charge, the cavity potential for our simulations is in excellent agreement with the studies of Rajamani et al., which were obtained using Widom insertion techniques. 34 It should also be pointed out that in the DCT study the cavitation potential given in Eq.(1) is fit to the simulation results of Rajamani et al.34 not the simulation results presented here (see the SI for an explanation of these differences). Without the inclusion of the electrochemical surface potential, the DCT for fixed charge ions dictates that the uncharged cavity PMF should produce the lowest free energy in the charged cavity series. In the case of the SPC/E simulations, this is clearly not the case. Our study shows that indeed small charges adsorb more than the uncharged cavity. These, small but real effects can only be accounted for through DCT by the inclusion of an electrochemical surface potential of corresponding to in Eq.(3), or , which is roughly half of the computed surface potential for SPC/E water.

It is important to recall that DCT theory assumes, in general, that the two terms that contribute to the solvation free energy are associated with the formation of a cavity followed by the charging process. The latter is assumed to be completely determined by dielectric response. From the point of view of hard-spheres in a dielectric, the asymmetry of cations and anions in the adsorption at the air-water interface is due only to the presence of the electrochemical surface potential. Recent work by Arslanargin and Beck has described this term in the context of molecular simulation as coming from both near and far-field electrostatic contributions. 19 The near-field is due to the local arrangements of water molecules around the ion, and the far-field to the electrostatic potential (surface potential) of the air-water interface. Both of these contributions are contained in in  3, which gives a free energy gain of , due to the electrochemical surface potential of , in excellent agreement with the value obtained by Arslanargin and Beck. 19

When comparing the DCT model to the simulation of a charged hard sphere in a Stockmayer fluid a very different picture emerges. Again, we see excellent agreement between the DCT and our simulations. However, we find that the uncharged cavity indeed has the most negative adsorption free energy, as expected. This is in stark contrast to the PMFs obtained for the SPC/E water model, in which the uncharged cavity does not exhibit the lowest free energy. Moreover, we have verified that cations and anions give identical free energy profiles as dictated from the symmetry of the Stockmayer fluid. Thus, we have provided direct evidence through molecular simulation that the principles for anion adsorption are different for the Stockmayer fluid. Only fully charged monovalent anions of unrealistic size will be able to adsorb to the air-water interface.

Horinek and co-workers have fit a point charge, soft-sphere interaction potentials to yield PMF that reproduce the experimentally measured surface tensions as a function of electrolyte concentration. 16 This, however comes with a price since other thermodynamic and structural properties will not be well described by modifications of the interaction potential. 16 On the other hand, it has been shown using ab initio based interaction potentials that we recover the experimentally determined local solvation structure around the anion. 35 Therefore, our present work suggests that the small adsorption predicted by Horinek and co-workers is due to the electrochemical surface potential of SPC/E water. 19 In other words, the electrochemical surface potential is implicitly contained in their simulated PMF, while in our treatment it appears explicitly, i.e., it is separated from the cavity and electrostatic self-energy contributions to the free energy. Thus, we conclude that the driving force for adsorption of large halide anions to the interface in point charge models of ions and water arises from different principles than in our ab initio MD simulation and PA-DCT theory. For example, we assert that fully charged, non-polarizable anions of realistic size will only adsorb if they couple to free energy at the air-water interface corresponding to a electrochemical surface potential of 19 In contrast, DFT and PA-DCT suggest it is only cavitation, aided by polarizability, that drives large anions to the air-water interface.

It is clear that future research needs to be focused on the role of the ab initio derived surface potential and how it influences the free energy profile for an ion in the vicinity of the vapor-liquid interface. There have been other independent studies that have speculated as to why an ion does not feel the full surface potential given by point charge models 36, 24, 19. These studies give hints as to why it could be possible for an ion in a quantum mechanical study to feel a very small and possibly negligible surface potential. Verification of these ideas with ab initio calculations is a subject for future research where, once again, proper simulation protocols must be developed in order to isolate the different contributions to the free energy of adsorption of ions to the air-water interface.

We have studied the driving forces for the adsorption of ions to the air-water interface using both molecular simulation and DCT for non-polarizable ions. Within the class of partial charge water models (e.g., SPC/E), the free energy of adsorption for non-polarizable, hard sphere ions has two contributions: cavitation, and the negative electrochemical surface potential of water. The microscopic origin of the cavitation energy is the perturbation to the hydrogen bond network produced when the ion is inside water. When the ion moves across the interface, the perturbation to the hydrogen bond network decreases, and so does the cavitation energy. In addition, the classical water models give rise to a negative free energy of adsorption that is associated with the electrochemical surface potential, which has been shown elsewhere to have both near and far-field contributions. 19 Our calculations suggests that the free energy due to the surface potential that an ion experiences is approximately half of the average electrostatic surface potential of SPC/E water, that amounts to a significant driving force for moving the ion across the interface. However, moving a hard (non-polarizable) ion completely into the low dielectric of air results in a huge electrostatic self-energy penalty that overwhelms the favorable contributions from the cavitational energy and the electrostatic surface potential. Simulations of the Stockmayer fluid, which by symmetry does not have a contribution from the electrochemical surface potential, have further emphasized the fundamental role played by the electrostatic surface potential in determining the driving force for ion adsorption within classical water models.

We anticipate that the effect of the electrochemical surface potential of empirical water models utilizing polarizable force fields will yield similar effects on ionic adsorption. It has been shown that the computed electrochemical surface potential for water obtained from simulations using a polarizable force field is similar to those obtained using non-polarizable point charge models. 37 Additionally, as argued in many recent studies, the currently available polarizable force fields tend to over-polarize even under bulk conditions. 26, 38, 35 There have been attempts to remedy this problem using different methods for screening the multipole electrostatics that yield dipole moments in good agreement with density functional theory studies. 37 However, even these modified models will predict too much adsorption (), leading to an incorrect, negative surface tension increment for solutions of large halides.

At this point, it is important to stress that the surface potential of classical water models is likely due to the nature of a point charge model in the vicinity of a broken symmetry. This has been argued in earlier work 28, 22, 23, 24 where it is shown that, as long as the width of the electron density is larger than that of the positive nuclear charge density, then the mean-inner potential, or surface potential, will be positive. Thus, the real average electrostatic potential across the air-water interface as calculated using the ab initio simulations is (again, note the difference in sign). In Ref. 24 it was determined that this average potential can be probed by high energy electron holography measurements, but does not translate directly into the free energy felt by an anion. Our previous ab initio study suggests ions feel a free energy due to the electrochemical surface potential that is small and possibly negligible.15

Insofar as charged hard-spheres in dielectric media are concerned, it would seem that we have a way of decomposing the interactions that drive ions to interfaces using molecular simulation and DCT. Unfortunately, this simplified view does not capture the degree of complexity that exists in real systems modeled either by sophisticated interaction potentials derived empirically or by quantum mechanics. Moreover, in the case of soft-spheres (e.g. Lennard-Jones interactions) it is much more difficult to make contact with DCT because one cannot easily separate the cavitation free energy from the dielectric response. Despite this complication, similar conclusions concerning the role of the electrochemical surface potential in driving ions to interfaces were obtained in a recent study of ion adsorption using the soft sphere potentials of Horinek et al19

Clearly, if the electrochemical surface potential is directly related to the surface potential of water (say, , as computed by DFT) it should strongly favor anion repulsion from the surface, and not adsorption as predicted by the classical water models. The relationship between the surface potential (which can be probed by high energy electrons) to the electrochemical surface potential that an ion would feel is still not apparent. In the present study, we have shown that we can isolate the free energy contribution due to an electrochemical surface potential in the PMF for the adsorption of a hard sphere ion to the air-water interface. The neglect of a free energy due to the electrochemical surface potential in the PA-DCT leads to excellent agreement with surface tension measurements and the PMF for iodide adsorption obtained using DFT. 15 These interesting results suggest that the real electrochemical surface potential that an ion feels is model dependent and is likely much smaller () than that obtained in simulations employing classical water models. Given that we can achieve near quantitative agreement between molecular simulation and DCT theory for fixed charge ions for the models employed in the present investigation, we are confident that similar connections can be made with ab initio simulation and PA-DCT.


It is a pleasure to thank Prof. Tom Beck for sharing preprints of his work and for candid discussions on this subject. We also acknowledge insightful conversations with Greg Schenter and Shawn Kathmann. This work was supported by National Science Foundation grant CHE-0431312 and by the U.S. Department of Energy‘s (DOE) Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated for the Department of Energy by Battelle. The calculations presented herein were performed with PNNLs Institutional Computing. MDB is grateful for the support of the Linus Pauling Distinguished Postdoctoral Fellowship Program at PNNL. Y.L. would like to acknowledge partial support by the CNPq, FAPERGS, INCT-FCx, and by the US-AFOSR under the grant FA9550-09-1-0283.



  • Jamtveit and Meakin 1999 Jamtveit, B., Meakin, P., Eds. Growth, Dissolution and Pattern Formation in Geosystems; Springer, 1999.
  • Heydweiller 1910 Heydweiller, A. Concerning the physical characteristics of solutions in correlation. II. Surface tension and elcetronic conductivity of watery salt solutions. Annalen Der Physik 1910, 33, 145–185.
  • Onsager and Samaras 1934 Onsager, L.; Samaras, N. N. T. The Surface Tension of Debye-Hückel Electrolytes. J. Chem. Phys. 1934, 2, 528–537.
  • Frumkin 1924 Frumkin, A. Phase interface powers and adsorption on the segregative surface air - Solution of anorganic electrolyte. Zeitschrift für Physikalische Chemie–Stochiometrie und Verwandtschaftslehre 1924, 109, 34–48.
  • Kunz et al. 2004 Kunz, W.; Henle, J.; Ninham, B. W. ’Zur Lehre von der Wirkung der Salze’ (About the Science of the Effect of Salts): Franz Hofmeister’s Historical Papers. Curr. Opin. Colloid In. 2004, 9, 19–37.
  • Perera and Berkowitz 1993 Perera, L.; Berkowitz, M. Stabilization Energies of Cl, Br, and I Ions in Water Clusters. J. Chem. Phys. 1993, 99, 4222–4224.
  • Jungwirth and Tobias 2001 Jungwirth, P.; Tobias, D. J. Molecular Structure of Salt Solutions: A New View of the Interface with Implications for Heterogeneous Atmospheric Chemistry. J. Phys. Chem. B 2001, 105, 10468–10472.
  • Levin 2009 Levin, Y. Polarizable Ions at Interfaces. Phys. Rev. Lett. 2009, 102, 147803.
  • Levin et al. 2009 Levin, Y.; dos Santos, A. P.; Diehl, A. Ions at the Air-Water Interface: An End to a Hundred-Year-Old Mystery? Phys. Rev. Lett. 2009, 103, 257802.
  • dos Santos et al. 2010 dos Santos, A. P.; Diehl, A.; Levin, Y. Surface Tensions, Surface Potentials, and the Hofmeister Series of Electrolyte Solutions. Langmuir 2010, 26, 10778–10783.
  • Ghosal et al. 2005 Ghosal, S.; Hemminger, J. C.; Bluhm, H.; Mun, B. S.; Hebenstreit, E. L. D.; Ketteler, G.; Ogletree, D. F.; Requejo, F. G.; Salmeron, M. Electron Spectroscopy of Aqueous Solution Interfaces Reveals Surface Enhancement of Halides. Science 2005, 307, 563–566.
  • Petersen and Saykally 2006 Petersen, P. B.; Saykally, R. J. Probing the Interfacial Structure of Aqueous Electrolytes with Femtosecond Second Harmonic Generation Spectroscopy. J. Phys. Chem. B 2006, 110, 14060–14073.
  • Baer et al. 2009 Baer, M. D.; Kuo, I.-F. W.; Bluhm, H.; Ghosal, S. Interfacial Behavior of Perchlorate versus Chloride Ions in Aqueous Solutions. J. Phys. Chem. B 2009, 113, 15843–15850.
  • Ho et al. 2003 Ho, C.; Tsao, H.; Sheng, Y. Interfacial tension of a salty droplet: Monte Carlo study. J. Chem. Phys. 2003, 119, 2369–2375.
  • Baer and Mundy 2011 Baer, M. D.; Mundy, C. J. Toward an Understanding of the Specific Ion Effect Using Density Functional Theory. J. Phys. Chem. lett. 2011, 2, 1088–1093.
  • Horinek et al. 2009 Horinek, D.; Herz, A.; Vrbka, L.; Sedlmeier, F.; Mamatkulov, S. I.; Netz, R. R. Specific Ion Adsorption at the Air/Water Interface: The Role of Hydrophobic Solvation. Chem. Phys. Lett. 2009, 479, 173–183.
  • Dang 2002 Dang, L. X. Computational Study of Ion Binding to the Liquid Interface of Water. J. Phys. Chem. B 2002, 106, 10388–10394.
  • Netz and Horinek 2012 Netz, R. R.; Horinek, D. Progress in Modelling of Ion Effects at the Vapor/Water Interface. Annu. Rev. Phys Chem. 2012, 63, 401–418.
  • Arslanargin and Beck 2012 Arslanargin, A.; Beck, T. L. Free energy partitioning analysis of the driving forces that determine ion density profiles near the water liquid-vapor interface. J. Chem. Phys. 2012, 136, 104503–12.
  • Paschek and Garcia 2004 Paschek, D.; Garcia, A. Reversible temperature and pressure denaturation of a protein fragment: A replica exchange molecular dynamics simulation study. Phys. Rev. Lett. 2004, 93, 238105.
  • Kuo et al. 2006 Kuo, I.-F. W.; Mundy, C. J.; Eggimann, B. L.; McGrath, M. J.; Siepmann, J. I.; Chen, B.; Vieceli, J.; Tobias, D. J. Structure and Dynamics of the Aqueous Liquid-Vapor Interface: A Comprehensive Particle-Based Simulation Study. J. Phys. Chem. B 2006, 110, 3738–3746.
  • Kathmann et al. 2008 Kathmann, S. M.; Kuo, I.-F. W.; Mundy, C. J. Electronic Effects on the Surface Potential at the Vapor-Liquid Interface of Water. J. Am. Chem. Soc. 2008, 130, 16556–16561.
  • Kathmann et al. 2009 Kathmann, S. M.; Kuo, I.-F. W.; Mundy, C. J. Electronic Effects on the Surface Potential at the Vapor-Liquid Interface of Water (vol 130, 16556, 2008). J. Am. Chem. Soc. 2009, 131, 17522–17522.
  • Kathmann et al. 2011 Kathmann, S. M.; Kuo, I.-F. W.; Mundy, C. J.; Schenter, G. K. Understanding the Surface Potential of Water. J. Phys. Chem. B 2011, 115, 4369–4377.
  • Mundy et al. 2009 Mundy, C. J.; Kuo, I.-F. W.; Tuckerman, M. E.; Lee, H.-S.; Tobias, D. J. Hydroxide Anion at the Air-Water Interface. Chem. Phys. Lett. 2009, 481, 2–8.
  • Wick 2009 Wick, C. D. Electrostatic Dampening Dampens the Anion Propensity for the Air-Water Interface. J. Chem. Phys. 2009, 131, 084715.
  • Leung 2010 Leung, K. Surface Potential at the Air-Water Interface Computed Using Density Functional Theory. J. Phys. Chem. Lett. 2010, 1, 496–499.
  • Wilson et al. 1988 Wilson, M.; Pohorille, A.; Pratt, L. Surface-potential of the Water Liquid Vapor Interface. J. Chem. Phys. 1988, 88, 3281–3285.
  • van der Spoel et al. 1998 van der Spoel, D.; van Maaren, P. J.; Berendsen, H. J. C. A systematic study of water models for molecular simulation: Derivation of water models optimized for use with a reaction field. J. Chem. Phys. 1998, 108, 10220–10230.
  • Jungwirth and Tobias 2006 Jungwirth, P.; Tobias, D. J. Specific Ion Effects at the Air/Water Interface. Chem. Rev. 2006, 106, 1259–1281.
  • Tamashiro and Constantino 2010 Tamashiro, M. N.; Constantino, M. A. Ions at the Water-Vapor Interface. J. Phys. Chem. B 2010, 114, 3583–3591.
  • Kumar et al. 1992 Kumar, S.; Bouzida, D.; Swendsen, R. H.; Kollman, P. A.; Rosenberg, J. M. The Weighted Histogram Analysis Method for Free-Energy Calculations on Biomolecules. 1. the Method. J. Comput. Chem. 1992, 13, 1011–1021.
  • Roux 1995 Roux, B. The Calculation of the Potential of Mean Force using Computer Simulations. Comput. Phys. Commun. 1995, 91, 275–282.
  • Rajamani et al. 2005 Rajamani, S.; Truskett, T.; Garde, S. Hydrophobic hydration from small to large lengthscales: Understanding and manipulating the crossover. Proc. Natl. Acad. Sci. USA 2005, 102, 9475–9480.
  • Fulton et al. 2010 Fulton, J. L.; Schenter, G. K.; Baer, M. D.; Mundy, C. J.; Dang, L. X.; Balasubramanian, M. Probing the Hydration Structure of Polarizable Halides: A Multiedge XAFS and Molecular Dynamics Study of the Iodide Anion. J. Phys. Chem. B 2010, 114, 12926–12937.
  • Harder and Roux 2008 Harder, E.; Roux, B. On the origin of the electrostatic potential difference at a liquid-vacuum interface. J. Chem. Phys. 2008, 129, 234706.
  • Wick et al. 2006 Wick, C. D.; Dang, L. X.; Jungwirth, P. Simulated surface potentials at the vapor-water interface for the KCl aqueous electrolyte solution. J. Chem. Phys. 2006, 125, 024706.
  • Rogers and Beck 2010 Rogers, D. M.; Beck, T. L. Quasichemical and Structural Analysis of Polarizable Anion Hydration. J. Chem. Phys. 2010, 132, 014505.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description