Molecular Dynamics Simulations of Sodium Dodecyl Sulfate Micelles in Water—The Effect of the Force Field

Molecular dynamic (MD) simulations of preassembled sodium dodecyl sulfate (SDS) micelles are carried out using three versions of GROMOS, as well as CHARMM36, OPLS-AA, and OPLS-UA force fields at different aggregation numbers and box sizes. The differences among force fields have little effect on the overall micelle structure of small aggregates of size 60 or 100, but for micelles of an aggregation number of 300 or higher, bicelle structures with ordered tails, rather than the more realistic rodlike or cylindrical micelles with disordered tails, occur when using versions of GROMOS45A3 or the OPLS-AA force fields that are adapted to model the sulfate head group atoms using methods given in the literature. We find that the Lennard-Jones (L-J) parameters for the sodium ions and the ionic oxygens of the SDS head group, as well as the water model, control the transition to bicelles, regardless of other L-J parameters. A closer binding of the sodium ions to the head group ionic oxygens screens the electrostatic repulsions more strongly, resulting in condensation of SDS head groups, leading to unphysical bicelles for GROMOS45A3 or the OPLS-AA force fields, when the aggregation number is large. A telltale sign that the sodium–oxygen interaction is too strong shows up in high nearest neighbor peaks (height >8 and height >20 for micelles with 60 and 100 surfactants, respectively) in the radial distribution functions (RDFs) of sodium ions to ionic oxygens. In the 100-surfactant micelles, the high RDF peak is accompanied by “crystal-like” layering of sodium ions onto the surface of the micelle. The distance between the sodium ions and micelle also depends on the number of waters binding to sodium ions in the presence of surfactant head groups, which depends on both the sodium ion and water models, and for the same sodium model increases as the water model is changed in the order: TIP4P, SPC/E, SPC, and TIP3P.