Molecular dynamics and metadynamics simulations of electrosprayed water nanodroplets including sodium bis(2-ethylhexyl)sulfosuccinate micelles

The behavior of aqueous solutions of sodium bis(2-ethylhexyl)sulfosuccinate (AOTNa) under conditions of electrospray ionization (ESI) has been investigated by molecular dynamics (MD) and well-tempered metadynamics (WTM) simulations at 300 K and 400 K. We have examined water droplets with initial fixed numbers of water molecules (1000) and AOT- anions (100), and with sodium cations in the range of 70-130. At 300 K, all charged droplets show the water evaporation rate increasing with the absolute value of the initial droplet charge state (Z), accompanied by ejection of an increasing number of solvated sodium ions or by expulsion of AOT- anions depending on the sign of Z and by fragmentation in the case of high |Z|. At 400 K, the water evaporation becomes more rapid and the fission process more extensive. In all cases, the AOTNa molecules, arranged as a direct micelle inside the aqueous system, undergo a rapid inversion in vacuo so that the hydrophilic heads and sodium ions surrounded by water molecules move toward the droplet interior. At the end of the 100-ns MD simulations, some water molecules remain within the aggregates at both temperatures. The subsequent metadynamics simulations accelerate the droplet evolution and show that all systems become anhydrous, in agreement with the experimental results of ESI mass spectrometry. This complete water loss is accompanied by sodium counterion emission for positively charged aggregates at 300 K. The analysis shows how the temperature and droplet charge state affect the populations of the generated surfactant aggregates, providing information potentially useful in designing future ESI experimental conditions.