Water nanodroplets charged with excess protons play a central role during electrospray ionization (ESI). In the current study molecular dynamics (MD) simulations were used for gaining insights into the nanodroplet behavior based on classical mechanics. The SPC/E water model was modified to permit the inclusion of protons as highly mobile point charges at minimum computational cost. A spherical trapping potential was assigned to every SPC/E oxygen, thereby allowing the formation of protonated water molecules. Within a tightly packed nanodroplet the individual potential wells merge to form a three-dimensional energy landscape that facilitates rapid proton hopping between water molecules. This approach requires short-range modifications to the standard Coulomb potential for modeling electrostatic proton-water interactions. Simulations on nanodroplets consisting of 1248 water molecules and 10 protons (radius, ca. 21 A) result in a proton diffusion coefficient that is in agreement with the value measured in bulk solution. Radial proton distributions extracted from 1 ns MD runs exhibit a large peak around 14 A, in addition to substantial population density closer to the droplet center. Similar radial distributions were found for nanodroplets charged with Na+ ions. This behavior is dramatically different from that expected on the basis of continuum electrostatic theory, which predicts that excess charge should be confined to a thin layer on the droplet surface. One important contributor to this effect seems to be the ordering of water molecules at the liquid/vacuum interface. This ordering results in an electrical double layer, generating a potential gradient that tends to pull positive charge carriers (such as protons, but also others such as Na+ ions) toward the droplet interior. This deviation from the widely assumed surface charge paradigm could have implications for the mechanism by which protonated analyte ions are formed during ESI.
