Molecular dynamics simulations of electrosprayed water nanodroplets: internal potential gradients, location of excess charge centers, and "hopping" protons.
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.
Ahadi
Elias
E
Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada.
Konermann
Lars
L
D004058:Diffusion; D008958:Models, Molecular; D008968:Molecular Conformation; D053758:Nanoparticles / Q000737:chemistry*; D011522:Protons; D015203:Reproducibility of Results; D012997:Solvents / Q000737:chemistry; D021241:Spectrometry, Mass, Electrospray Ionization; D055672:Static Electricity; D014867:Water / Q000737:chemistry*
D011522:Protons; D012997:Solvents; D014867:Water