Journal
SOFT MATTER
Volume 8, Issue 8, Pages 2385-2397Publisher
ROYAL SOC CHEMISTRY
DOI: 10.1039/c1sm06787g
Keywords
-
Categories
Funding
- National Science Foundation [CHE 09-46358, CNS 09-58854]
- Longhorn at TACC [TG-MCA93S020, TG-ASC090088]
- Princeton Center for Complex Materials [NSF-DMR 0819860]
- Department of Energy, Office of Basic Energy Sciences [DE-SC0002128]
- Direct For Mathematical & Physical Scien
- Division Of Materials Research [819860, 1120901] Funding Source: National Science Foundation
- Division Of Chemistry
- Direct For Mathematical & Physical Scien [0946358] Funding Source: National Science Foundation
- Division Of Computer and Network Systems
- Direct For Computer & Info Scie & Enginr [0958854] Funding Source: National Science Foundation
Ask authors/readers for more resources
Due to the relatively long time scales inherent to ionic surfactant self-assembly (>mu s), an aggressive computational approach is needed to obtain converged data on micellar solutions. This work presents a study of micellization using a coarse-grained (CG) model of aqueous ionic surfactants in replicated molecular dynamics (MD) simulations run on graphics processing unit hardware. The performance of our implementation of the CG model with electrostatics into the HOOMD-Blue GPU-accelerated MD software package is comparable to that of a modest sized cluster running a highly optimized parallel CPU code. From 0.36 ms of cumulative trajectory data, we are able to predict equilibrium thermodynamic and morphological properties of ionic surfactant micellar solutions. Estimating the critical micelle concentrations (CMC) from the free monomer (rho(1)) and premicellar concentrations obtained from simulations of sodium hexyl sulfate (S6S, CMC of 460 +/- 6 mM) at high (1 M) concentration, a value in good agreement with experimental results is obtained; however, the same method applied to simulations of sodium nonyl sulfate (S9S, rho(1) of 2.4 +/- 0.01 mM) and sodium dodecyl sulfate (SDS, rho(1) of 0.02 +/- 0.01 mM) at the same total concentration systematically underestimates the CMCs. An alternative method for calculating the CMC is presented, where the free monomer concentration computed from high concentration CG-MD data is used as the input to a simple theoretical model which can be used to extrapolate to a more accurate prediction of the CMC. Better agreement between the empirical and predicted CMC is obtained from this theory for S9S (28.7 +/- 0.3 mM) and SDS (3.32 +/- 0.04 mM), though the CMC for S6S is slightly underestimated (304 +/- 3 mM). We also present statistically converged morphological data, including aggregation number distributions and the principal components of the gyration tensor. This data suggest a transition from spherical micelles to rod-like at a specific aggregation number, which increases with increasing hydrocarbon length.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available