Monte carlo simulations of micellization in model ionic surfactants: application to sodium dodecyl sulfate

A lattice model for ionic surfactants with explicit counterions is proposed for which the micellization behavior can be accurately determined from grand canonical Monte Carlo simulations. The model is characterized by a few parameters that can be adjusted to represent various linear surfactants with ionic headgroups. The model parameters have a clear physical interpretation and can be obtained from experimental data unrelated to micellization, namely, geometric information and solubilities of tail segments. As a specific example, parameter values for sodium dodecyl sulfate were obtained by optimizing for the solubility of hydrocarbons in water and the structural properties of dodecane. The critical micelle concentration (cmc), average aggregation number, degree of counterion binding, and their dependence on temperature were determined from histogram reweighting grand canonical Monte Carlo simulations and were compared to experimental results. The model gives the correct trend and order of magnitude for all quantities but underpredicts the cmc and aggregation number. We suggest ways to modify the model that may improve agreement with experimental values.     

A simplified molecular-thermodynamic model of a microphase-separated ionic gel swollen in salt solution

A molecular-thermodynamic model [A.I. Victorov, C.J. Radke, J.M. Prausnitz, Mol. Phys., 2004, submitted for publication] of diblock copolymer ionic gels swollen in brine is simplified by deriving asymptotic expressions for electrostatic terms. This model derived recently from the self-consistent field theory in the limit of strong segregation gives thermodynamic functions for gels of lamellar, bicontinuous, cylindrical and spherical morphologies and details the gel structure including equilibrium microdomain spacing, distribution of mobile ions, polymer segments, and the electric potential across the microdomains. The model reflects the copolymer chain composition, length, rigidity, ionization degree and the effective polymer–polymer and polymer–solvent interactions. Several asymptotic regimes are considered that lead to simple formulae for the solvent chemical potential. Applicability of asymptotic relations is tested. Equilibrium uptakes of salt are calculated for gels of varying ionic charge over a wide range of solution salinity.     

Aggregation work at polydisperse micellization: ideal solution and "dressed micelle" models comparing to molecular dynamics simulations

General thermodynamic relations for the work of polydisperse micelle formation in the model of ideal solution of molecular aggregates in nonionic surfactant solution and the model of "dressed micelles" in ionic solution have been considered. In particular, the dependence of the aggregation work on the total concentration of nonionic surfactant has been analyzed. The analogous dependence for the work of formation of ionic aggregates has been examined with regard to existence of two variables of a state of an ionic aggregate, the aggregation numbers of surface active ions and counterions. To verify the thermodynamic models, the molecular dynamics simulations of micellization in nonionic and ionic surfactant solutions at two total surfactant concentrations have been performed. It was shown that for nonionic surfactants, even at relatively high total surfactant concentrations, the shape and behavior of the work of polydisperse micelle formation found within the model of the ideal solution at different total surfactant concentrations agrees fairly well with the numerical experiment. For ionic surfactant solutions, the numerical results indicate a strong screening of ionic aggregates by the bound counterions. This fact as well as independence of the coefficient in the law of mass action for ionic aggregates on total surfactant concentration and predictable behavior of the "waterfall" lines of surfaces of the aggregation work upholds the model of "dressed" ionic aggregates.     

The use of multidimensional Franck-Condon simulations to assess model chemistries: a case study on phenol

Multidimensional Franck-Condon simulations of the dispersed fluorescence spectra of phenol generated with geometries obtained from the highly correlated post-Hartree-Fock methods CASSCF, MRCI, and SACCI are presented. While the simulations based on CASSCF and MRCI optimized geometries are very similar to each other and fail to reproduce the experimentally measured intensities faithfully, the simulations obtained from SACCI optimized geometries are very close to the experimental spectra. The code developed for the multidimensional Franck-Condon simulations is described. It is shown that the integral storage problem common to the evaluation of multidimensional Franck-Condon integrals can be overcome by saving all quantities needed to disk. This strategy allows the code to run on computers with limited resources and is very well suited for application to molecules with a very large number of vibrational modes.     

Macroion solutions in the cell model studied by field theory and Monte Carlo simulations

Aqueous solutions of charged spherical macroions with variable dielectric permittivity and their associated counterions are examined within the cell model using a field theory and Monte Carlo simulations. The field theory is based on separation of fields into short- and long-wavelength terms, which are subjected to different statistical-mechanical treatments. The simulations were performed by using a new, accurate, and fast algorithm for numerical evaluation of the electrostatic polarization interaction. The field theory provides counterion distributions outside a macroion in good agreement with the simulation results over the full range from weak to strong electrostatic coupling. A low-dielectric macroion leads to a displacement of the counterions away from the macroion.     

Thermodynamic properties of model solids with short-ranged potentials from Monte Carlo simulations and perturbation theory

Monte Carlo simulations have been performed to determine the excess energy and the equation of state of fcc solids with Sutherland potentials for wide ranges of temperatures, densities, and effective potential ranges. The same quantities have been determined within a perturbative scheme by means of two procedures: (i) Monte Carlo simulations performed on the reference hard-sphere system and (ii) second-order Barker-Henderson perturbation theory. The aim was twofold: on the one hand, to test the capability of the "exact" MC-perturbation theory of reproducing the direct MC simulations and, on the other hand, the reliability of the Barker-Henderson perturbation theory, as compared with direct MC simulations and MC-perturbation theory, to determine the thermodynamic properties of these solids depending on temperature, density, and potential range. We have found that the simulation data for the excess energy obtained from the two procedures are in close agreement with each other. For the equation of state, the results from the MC-perturbation procedure also agree well with the direct MC simulations except for very low temperatures and extremely short-ranged potentials. Regarding the Barker-Henderson perturbation theory, we have found that in general the second-order approximation does not provide significant improvement over the first-order one.     

Refining of nonionic surfactant micellization theory based on the law of mass action

Two approaches to determining critical micelle concentration (CMC) are assessed, i.e., from the inflection point in the curve for the concentration dependence of the degree of micellization and as K^1/(1–n), where K is the constant of the law of mass action and n is the aggregation number. The latter approach makes the theory simpler, while the former explicitly expresses the critical degree of micellization via the aggregation number. The concentrations of monomers and micelles are analyzed as functions of the overall concentration of a surfactant in a micellar solution. These functions look much simpler in the graphical form as compared with their complex exact analytical representation. This has resulted in derivation of simple analytical approximations for these functions, with these approximations being useful for calculations. The concentration dependence of the surfactant diffusion coefficient has been considered based on these approximations. It turned out that this dependence not only provides the known method for determining the diffusion coefficient of micelles, but also gives the possibility in principle to determine the aggregation number from the slope of the dependence of the diffusion coefficient on the inverse concentration (counted from the CMC in the CMC units). This new method for determining the aggregation number has been tested using the literature data on the diffusion coefficient of penta(ethylene glycol)-1-hexyl ether in an aqueous solution.     

Complementary use of inter-laboratory certification study data

The agreed purpose of an inter-laboratory certification study is to characterise a certified reference material (CRM). In addition to this, any (successful) certification study may also be utilised as a laboratory/method performance study by the participants. The intention of this note is to emphasise the complementary use of certification study data by the participants of the certification study and to explain the options involved. On this occasion, an apparent paradox arising in the interpretation of certification study data is resolved. The procedures for bias correction and estimation of bias-related uncertainty contributions discussed on this occasion are, however, generally applicable to normal CRM use. Received: 4 February 1999 / Revised: 16 July 1999 / Accepted: 21 July 1999     

Complementary use of inter-laboratory certification study data

The agreed purpose of an inter-laboratory certification study is to characterise a certified reference material (CRM). In addition to this, any (successful) certification study may also be utilised as a laboratory/method performance study by the participants. The intention of this note is to emphasise the complementary use of certification study data by the participants of the certification study and to explain the options involved. On this occasion, an apparent paradox arising in the interpretation of certification study data is resolved. The procedures for bias correction and estimation of bias-related uncertainty contributions discussed on this occasion are, however, generally applicable to normal CRM use. Received: 4 February 1999 / Revised: 16 July 1999 / Accepted: 21 July 1999     

Using Monte Carlo simulations and self-consistent field theory to design interfacially active copolymers

We use both Monte Carlo computer simulations and numerical self-consistent field lattice calculations to determine the behavior of copolymers at penetrable and impenetrable interfaces. These computational techniques are useful as “design tools”: they allow us to systemically vary the copolymer architecture, determine optimal structures for specific applications, and establish guidelines for fabricating copolymers that yield the desired interfacial properties. We illustrate this principle with three different examples. In the first study, we combine the techniques to design copolymer compatibilizers that enhance the strength of immiscible polymer blends. These copolymers contain teeth that associate across the penetrable interface between the phase-separated regions and form a “molecular velcro” that effectively binds the regions together. In the case of impenetrable interfaces, we determine how the copolymer sequence distribution affects the structure of a layer of copolymers grafted onto a solid surface. The results indicate how to control the morphology of the layer and the surface properties of the substrate, by varying the microstructure of the grafted copolymers. Finally, we design a polymer channel that “opens” and “closes” in response to changes in the pH and quality of the surrounding solvent. The channel is formed from polyacid chains that are anchored onto a solid surface. Due to these properties, the system can be used for controlled release or sensor devices.     

The use of Monte Carlo simulations to evaluate kinetic data and analytic approximations

Experimental kineticists are always faced with the problem of reducing kinetic data to extract physically meaningful information. A particularly vexing problem arises when different models reproduce the data but yield different values for the physical parameters. For over forty-five years Monte Carlo simulation techniques have been used to study the statistical behavior of parameters extracted from data. Not only do these simulations provide realistic uncertainties, correlation coefficients, and confidence envelopes, but they also provide insight into the nature of the model. These insights may be obtained by viewing two-dimensional scatter plots of the fractional changes of the parameters and one-dimensional histograms of the distributions of the changes in the parameters. Monte Carlo simulations are illustrated with examples from OH+CH₄ → CH₃+H₂O and the high-pressure rate coefficient for methyl-methyl association. A more complex problem involves models for pressure-dependent rate coefficients in the falloff region. We have modeled methyl-methyl association with five of the most current analytic approximations for behavior in the falloff region. All of these reproduce the data to within their uncertainties. However, when Monte Carlo techniques are applied the correlations between the parameters and the nonlinear nature of their behavior become evident. We postulate that the statistical behavior of the parameters of a model may be used to distinguish one model from another and, thereby, identify those analytic approximations that hold promise for further investigation and utilization. Finally, the recent advent of high-speed workstations implies that Monte Carlo simulations should become a routine part of the analysis of kinetic data. © 1997 John Wiley & Sons, Inc. Int J Chem Kinet 29: 803–817, 1997     

Kinetics of micellization: its significance to technological processes

The association of many classes of surface active molecules into micellar aggregates is a well-known phenomenon. Micelles are often drawn as static structures of spherical aggregates of oriented molecules. However, micelles are in dynamic equilibrium with surfactant monomers in the bulk solution constantly being exchanged with the surfactant molecules in the micelles. Additionally, the micelles themselves are continuously disintegrating and reforming. The first process is a fast relaxation process typically referred to as τ₁. The latter is a slow relaxation process with relaxation time τ₂. Thus, τ₂ represents the entire process of the formation or disintegration of a micelle. The slow relaxation time is directly correlated with the average lifetime of a micelle, and hence the molecular packing in the micelle, which in turn relates to the stability of a micelle. It was shown earlier by Shah and coworkers that the stability of sodium dodecyl sulfate (SDS) micelles plays an important role in various technological processes involving an increase in interfacial area, such as foaming, wetting, emulsification, solubilization and detergency. The slow relaxation time of SDS micelles, as measured by pressure-jump and temperature-jump techniques was in the range of 10⁻⁴–10¹ s depending on the surfactant concentration. A maximum relaxation time and thus a maximum micellar stability was found at 200 mM SDS, corresponding to the least foaming, largest bubble size, longest wetting time of textile, largest emulsion droplet size and the most rapid solubilization of oil. These results are explained in terms of the flux of surfactant monomers from the bulk to the interface, which determines the dynamic surface tension. The more stable micelles lead to less monomer flux and hence to a higher dynamic surface tension. As the SDS concentration increases, the micelles become more rigid and stable as a result of the decrease in intermicellar distance. The smaller the intermicellar distance, the larger the Coulombic repulsive forces between the micelles leading to enhanced stability of micelles (presumably by increased counterion binding to the micelles). The Center for Surface Science & Engineering at the University of Florida has developed methods using stopped-flow and pressure-jump with optical detection to determine the slow relaxation time of micelles of nonionic surfactants. The results show relaxation times τ₂ in the range of seconds for Triton X-100 to minutes for polyoxyethylene alkyl ethers. The slow relaxation times are much longer for nonionic surfactants than for ionic surfactants, because of the absence of ionic repulsion between the head groups. The observed relaxation time τ₂ was related to dynamic surface tension and foaming experiments. A slow break-up of micelles, (i.e. a long relaxation time τ₂) corresponds to a high dynamic surface tension and low foamability, whereas a fast break-up of micelles, leads to a lower dynamic surface tension and higher foamability. In conclusion, micellar stability and thus the micellar break-up time is a key factor in controlling technological processes involving a rapid increase in interfacial area, such as foaming, wetting, emulsification and oil solubilization. First, the available monomers adsorb onto the freshly created interface. Then, additional monomers must be provided by the break-up of micelles. Especially when the free monomer concentration is low, as indicated by a low CMC, the micellar break-up time is a rate limiting step in the supply of monomers, which is the case for many nonionic surfactant solutions. Therefore, relaxation time data of surfactant solutions enables us to predict the performance of a given surfactant solution. Moreover, the results suggest that one can design appropriate micelles with specific stability or τ₂ by controlling the surfactant structure, concentration and physico-chemical conditions, as well as by mixing anionic/cationic or ionic/nonionic surfactants for a desired technological application.     

Molecular-thermodynamic theory of micellization of multicomponent surfactant mixtures: 1. Conventional (pH-Insensitive) surfactants

A molecular-thermodynamic (MT) theory was developed to model the micellization of mixtures containing an arbitrary number of conventional (pH-insensitive) surfactants. The theory was validated by comparing predicted and experimental cmc's of ternary surfactant mixtures, yielding results that were comparable to, and sometimes better than, the cmc's determined using regular solution theory. The theory was also used to model a commercial nonionic surfactant (Genapol UD-079), which was modeled as a mixture of 16 surfactant components. The predicted cmc agreed well with the experimental cmc, and the monomer concentration was predicted to increase significantly above the cmc. In addition, the monomer and the micelle compositions were predicted to vary significantly with surfactant concentration. These composition variations were rationalized in terms of competing steric and entropic effects and a micelle shape transition near the cmc. To understand the packing constraints imposed on ternary surfactant mixtures better, the maximum micelle radius was also examined theoretically. The MT theory presented here represents the first molecular-based theory of the micellization behavior of mixtures of three or more conventional surfactants. In article 2 of this series, the MT theory will be extended to model the micellization of mixtures of conventional and pH-sensitive surfactants.     

Lattice Monte Carlo simulations of phase separation and micellization in supercritical CO2/surfactant systems: effect of CO2 density

Lattice Monte Carlo simulations are used to study the effect of nonionic surfactant concentration and CO2 density on the micellization and phase equilibria of supercritical CO2/surfactant systems. The interaction parameter for carbon dioxide is obtained by matching the critical temperature of the model fluid with the experimental critical temperature. Various properties such as the critical micelle concentration and the size, shape, and structure ofmicelles are calculated, and the phase diagram in the surfactant concentration-CO2 density space is constructed. On increasing the CO2 density, we find an increase in the critical micelle concentration and a decrease in the micellar size; this is consistent with existing experimental results. The variation of the micellar shape and structure with CO2 density shows that the micelles are spherical and that the extension of the micellar core increases with increasing micellar size, while the extension of the micellar corona increases with increasing CO2 density. The predicted phase diagram is in qualitative agreement with experimental phase diagrams for nonionic surfactants in carbon dioxide.     

Some aspects of theory and simulations of polymers in steady flows

We present a study of various properties of bead spring chains in steady flows. The Langevin equation of the normal modes of the chain is solved by Fourier transformation. From the resulting power spectrum, the autocorrelation functions of all configuration-dependent quantities can be calculated. In equilibrium, the influence of the bead masses on the short-time dynamics is discussed. The influence of different flow fields (shear, elongational and Kramers potential flow) on the mean-square chain dimension is calculated. A comparison with results obtained from non-equilibrium molecular dynamics and Monte Carlo calculations is made. Finally, the influence of shear flow on the configurational and rheological properties of cyclic polymers and on the excluded volume behavior of linear chains is examined.     

Correlations in anticrossing spectra and scattering theory: Numerical simulations

For the interpretation of experimental results of anticrossing spectroscopy in molecules, the results of numerical simulations are presented. We have focussed our attention on the visibility of the correlation hole. In the formulation of the model, we combine scattering theory and a random-matrix model.     

Binary hairy nanoparticles: Recent progress in theory and simulations

Binary polymer brushes, including mixed homopolymer brushes and diblock copolymer brushes, are an attractive class of environmentally responsive nanostructured materials. Owing to microphase separation of the two chemically distinct components in the brush, multifaceted nanomaterials with functionalized and patterned surfaces can be obtained. This review summarizes recent progress on the theory and simulations related to binary polymer brushes grafted to flat, spherical, and cylindrical substrates, with a focus on patterned morphologies of multifaceted hairy nanoparticles, an intriguing class of hybrid nanostructured particles (e.g., nanospheres and nanorods). In particular, powerful field theory and particle-based simulations suitable for revealing novel structures on these patterned surfaces, including self-consistent field theory and dissipative particle dynamics simulations, are emphasized. The unsolved yet critical issues in this research field, such as dynamic response of binary polymer brushes to environmental stimuli and the hierarchical self-assembly of binary hairy nanoparticles, are briefly discussed. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2014 , 52, 1583–1599