A multiscale model for kinetics of formation and disintegration of spherical micelles.

Dynamics of self-assembly and structural transitions in surfactant systems often involve a large span of length and time scales. A comprehensive understanding of these processes requires development of models connecting phenomena taking place on different scales. In this paper, we develop a multiscale model for formation and disintegration of spherical nonionic micelles. The study is performed under the assumption that the dominant mechanism of micelle formation (disintegration) is a stepwise addition (removal) of single monomers to (from) a surfactant aggregate. Different scales of these processes are investigated using a combination of coarse-grained molecular dynamics simulations, analytical and numerical solution of stochastic differential equations, and a numerical solution of kinetic equations. The removal of a surfactant from an aggregate is modeled by a Langevin equation for a single reaction coordinate, the distance between the centers of mass of the surfactant and the aggregate, with parameters obtained from a series of constrained molecular dynamics simulations. We demonstrate that the reverse process of addition of a surfactant molecule to an aggregate involves at least two additional degrees of freedom, orientation of the surfactant molecule and micellar microstructure. These additional degrees of freedom play an active role in the monomer addition process and neglecting their contribution leads to qualitative discrepancies in predicted surfactant addition rates. We propose a stochastic model for the monomer addition which takes the two additional degrees of freedom into account and extracts the model parameters from molecular dynamics simulations. The surfactant addition rates are determined from Brownian dynamics simulations of this model. The obtained addition and removal rates are then incorporated into the kinetic model of micellar formation and disintegration.

[1]  H. Kramers.  Brownian motion in a field of force and the diffusion model of chemical reactions , 1940 .

[2]  Jason E. Butler,et al.  Simulations of concentrated suspensions of rigid fibers: relationship between short-time diffusivities and the long-time rotational diffusion. , 2005, The Journal of chemical physics.

[3]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[4]  Kenneth S. Schmitz,et al.  Role of orientation constraints and rotational diffusion in bimolecular solution kinetics , 1972 .

[5]  Eric S. G. Shaqfeh,et al.  Dynamic simulations of the inhomogeneous sedimentation of rigid fibres , 2002, Journal of Fluid Mechanics.

[6]  A. Halperin,et al.  Polymeric micelles: their relaxation kinetics , 1989 .

[7]  Richard L. Rowley,et al.  Diffusion and viscosity equations of state for a Lennard-Jones fluid obtained from molecular dynamics simulations , 1997 .

[8]  U. Kaatze,et al.  ON THE KINETICS OF THE FORMATION OF SMALL MICELLES. 2. EXTENSION OF THE MODEL OF STEPWISE ASSOCIATION , 1997 .

[9]  Walter H. Stockmayer,et al.  Kinetics of Diffusion‐Controlled Reaction between Chemically Asymmetric Molecules. I. General Theory , 1971 .

[10]  Alan E. Mark,et al.  Molecular dynamics simulation of the kinetics of spontaneous micelle formation , 2000 .

[11]  D. Shah,et al.  Importance of micellar kinetics in relation to technological processes. , 2002, Journal of colloid and interface science.

[12]  Peter A. J. Hilbers,et al.  Computer simulations of surfactant self-assembly. , 1993 .

[13]  A. D. Mackie,et al.  Aggregation Behavior of a Lattice Model for Amphiphiles , 1997 .

[14]  H. Risken.  The Fokker-Planck equation : methods of solution and applications , 1985 .

[15]  G. Voth,et al.  A theory for the activated barrier crossing rate constant in systems influenced by space and time dependent friction , 1994 .

[16]  T. A. Hatton,et al.  Dynamics of self-assembled surfactant systems , 1998 .

[17]  Jun Liu,et al.  Simulations of Micelle Self-Assembly in Surfactant Solutions , 1996 .

[18]  J. Ávalos,et al.  Prediction of the critical micelle concentration in a lattice model for amphiphiles using a single-chain mean-field theory. , 2005, The Journal of chemical physics.

[19]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[20]  I. Szleifer,et al.  SELF-ASSEMBLY OF MODEL NONIONIC AMPHIPHILIC MOLECULES , 1999 .

[21]  H. Hasimoto.  On the periodic fundamental solutions of the Stokes equations and their application to viscous flow past a cubic array of spheres , 1959, Journal of Fluid Mechanics.

[22]  Jan D. Miller,et al.  The effect of cosurfactants on sodium dodecyl sulfate micellar structures at a graphite surface , 2006 .

[23]  W. Briels,et al.  THE CALCULATION OF FREE-ENERGY DIFFERENCES BY CONSTRAINED MOLECULAR-DYNAMICS SIMULATIONS , 1998 .

[24]  S. Kwon,et al.  Topological transition in aqueous nonionic micellar solutions. , 2002, Physical review letters.

[25]  C. W. Gardiner,et al.  Handbook of stochastic methods - for physics, chemistry and the natural sciences, Second Edition , 1986, Springer series in synergetics.

[26]  J. Mccammon,et al.  Brownian dynamics simulation of diffusion‐influenced bimolecular reactions , 1984 .

[27]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[28]  D. Ermak,et al.  Brownian dynamics with hydrodynamic interactions , 1978 .

[29]  P. Bolhuis,et al.  Accurate free energies of micelle formation. , 2005, The journal of physical chemistry. B.

[30]  I. Kevrekidis,et al.  Coarse molecular dynamics of a peptide fragment: Free energy, kinetics, and long-time dynamics computations , 2002, physics/0212108.

[31]  M. Teubner,et al.  Theory of ultrasonic absorption in micellar solutions , 1979 .

[32]  A. Mark,et al.  Coarse grained model for semiquantitative lipid simulations , 2004 .

[33]  R. C. Reeder,et al.  A Coarse Grain Model for Phospholipid Simulations , 2001 .

[34]  A. Szabó,et al.  Diffusion-controlled bimolecular reaction rates. The effect of rotational diffusion and orientation constraints. , 1981, Biophysical journal.

[35]  J. Lang,et al.  Theory of the kinetics of micellar equilibria and quantitative interpretation of chemical relaxation studies of micellar solutions of ionic surfactants , 1976 .

[36]  Noriyuki Yoshii,et al.  A molecular dynamics study of free energy of micelle formation for sodium dodecyl sulfate in water and its size distribution. , 2006, The Journal of chemical physics.

[37]  Herman J. C. Berendsen,et al.  Simulation of Water Transport through a Lipid Membrane , 1994 .

[38]  S. Edwards,et al.  The Theory of Polymer Dynamics , 1986 .

[39]  P. Bolhuis,et al.  Prediction of an autocatalytic replication mechanism for micelle formation. , 2006, Physical review letters.

[40]  Håkan Wennerström,et al.  The Colloidal Domain: Where Physics, Chemistry, Biology and Technology Meet , 1994 .

[41]  Michiel Sprik,et al.  Free energy from constrained molecular dynamics , 1998 .

[42]  D. Baker,et al.  A general expression for bimolecular association rates with orientational constraints , 2002, physics/0401127.

[43]  Benoît Roux,et al.  Ion transport in a gramicidin-like channel: dynamics and mobility , 1991 .

[44]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[45]  E. A. G. Aniansson,et al.  Kinetics of step-wise micelle association , 1974 .

[46]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[47]  Koert N. J. Burger,et al.  Greasing Membrane Fusion and Fission Machineries , 2000, Traffic.

[48]  Ioannis G Kevrekidis,et al.  Coarse-grained kinetic computations for rare events: application to micelle formation. , 2005, The Journal of chemical physics.

[49]  E. Ruckenstein,et al.  Theory of surfactant self-assembly : a predictive molecular thermodynamic approach , 1991 .

[50]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[51]  J. Israelachvili.  Intermolecular and surface forces , 1985 .

[52]  Three eras of micellization. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[53]  M. Teubner,et al.  On the kinetics of micellization in aqueous solutions , 1980 .

[54]  Reinhard Lipowsky,et al.  Computer simulations of bilayer membranes - self-assembly and interfacial tension. , 1998 .