Molecular dynamics simulation of self-assembly of n-decyltrimethylammonium bromide micelles

In this paper, a molecular dynamics simulation of surfactant self-assembly using realistic atomistic models is presented. The simulations are long enough to enable the observation of several processes leading to equilibrium, such as monomer addition and detachment, micelle dissolution, and micelle fusion. The self-assembly of DeTAB surfactants takes place in three stages: fast aggregation of monomers to form small disordered oligomers; ripening process by which larger aggregates grow at the expense of smaller ones; slower stage involving collisions between large micelles. The first two stages were described well by a simple kinetic model with a size-independent rate constant estimated from the self-diffusion coefficient and collision radius of an isolated monomer. The average cluster size, area per headgroup, degree of counterion dissociation, and critical micelle concentration estimated from the simulation are in reasonable agreement with experimental values. An all-atom and united-atom surfactant model were compared, and the results were seen to be almost independent of the choice of model. DeTAB micelles are spheroidal, with a hydrophobic core composed of tail atoms surrounded by a hydrophilic corona of head atoms. A Stern layer composed of bromide counterions was also identified. Water molecules solvate the counterions and the head atoms, penetrating into the micelle up to the location of the atom connecting the head to the aliphatic tail, in agreement with recent experimental observations.     

Decay of standing foams: drainage,coalescence and collapse

A summary of recent theoretical work on the decay of foams is presented. In a series of papers, we have proposed models for the drainage, coalescence and collapse of foams with time. Each of our papers dealt with a different aspect of foam decay and involved several assumptions. The fundamental equations, the assumptions involved and the results obtained are discussed in detail and presented within a unified framework.Film drainage is modeled using the Reynolds equation for flow between parallel circular disks and film rupture is assumed to occur when the film thickness falls below a certain critical thickness which corresponds to the maximum disjoining pressure. Fluid flow in the Plateau border channels is modeled using a Hagen-Poiseuille type flow in ducts with triangular cross-section.The foam is assumed to be composed of pentagonal dodecahedral bubbles and global conservation equations for the liquid, the gas and the surfactant are solved to obtain information about the state of the decaying foam as a function of time. Homogeneous foams produced by mixing and foams produced by bubbling (pneumatic foams) are considered. It is shown that a draining foam eventually arrives at a mechanical equilibrium when the opposing forces due to gravity and the Plateau-border suction gradient balance each other. The properties of the foam in this equilibrium state can be predicted from the surfactant and salt concentration in the foaming solution, the density of the liquid and the bubble radius.For homogeneous foams, it is possible to have conditions under which there is no drainage of liquid from the foam. There are three possible scenarios at equilibrium: separation of a single phase (separation of the continuous phase liquid by drainage or separation of the dispersed phase gas via collapse), separation of both phases (drainage and collapse occurs) or no phase separation (neither drainage nor collapse occurs). It is shown that the phase behavior depends on a single dimensionless group which is a measure of the relative magnitudes of the gravitational and capillary forces. A generalized phase diagram is presented which can be used to determine the phase behavior.For pneumatic foams, the effects of various system parameters such as the superficial gas velocity, the bubble size and the surfactant and salt concentrations on the rate of foam collapse and the evolution of liquid fraction profile are discussed. The steady state height attained by pneumatic foams when collapse occurs during generation is also evaluated.Bubble coalescence is assumed to occur due to the non-uniformity in the sizes of the films which constitute the faces of the polyhedral bubbles. This leads to a non-uniformity of film-drainage rates and hence of film thicknesses within any volume element in the foam. Smaller films drain faster and rupture earlier, causing the bubbles containing them to coalesce. This leads to a bubble size distribution in the foam, with the bubbles being larger in regions where greater coalescence has occurred.The formation of very stable Newton black films at high salt and surfactant concentrations is also explained.     

Micelle–monomer equilibria in solutions of ionic surfactants and in ionic–nonionic mixtures: A generalized phase separation model

On the basis of a detailed physicochemical model, a complete system of equations is formulated that describes the equilibrium between micelles and monomers in solutions of ionic surfactants and their mixtures with nonionic surfactants. The equations of the system express mass balances, chemical and mechanical equilibria. Each nonionic surfactant is characterized by a single thermodynamic parameter — its micellization constant. Each ionic surfactant is characterized by three parameters, including the Stern constant that quantifies the counterion binding. In the case of mixed micelles, each pair of surfactants is characterized with an interaction parameter, β, in terms of the regular solution theory. The comparison of the model with experimental data for surfactant binary mixtures shows that β is constant — independent of the micelle composition and electrolyte concentration. The solution of the system of equations gives the concentrations of all monomeric species, the micelle composition, ionization degree, surface potential and mean area per head group. Upon additional assumptions for the micelle shape, the mean aggregation number can be also estimated. The model gives quantitative theoretical interpretation of the dependence of the critical micellization concentration (CMC) of ionic surfactants on the ionic strength; of the CMC of mixed surfactant solutions, and of the electrolytic conductivity of micellar solutions. It turns out, that in the absence of added salt the conductivity is completely dominated by the contribution of the small ions: monomers and counterions. The theoretical predictions are in good agreement with experimental data.     

Investigation on properties of aqueous foams stabilized by aliphatic alcohols and polypropylene glycol

Foams stabilized by nonionic surfactants are usually moderately stable due to high drainage rate and intense bubble coalescence and coarsening. This study aimed to investigate comparatively the foam properties of aliphatic alcohols (methyl isobutyl carbinol (MIBC) and 2-octanol) and polypropylene glycol (PPG400). Experiments were conducted using the FoamScan method at various surfactant concentrations and gas flow rates where the foam volume, liquid content of foam and foam half-life were determined. The results showed that both foamability and foam stability of surfactant solution increased with increasing gas flow rate and surfactant concentration for all tested surfactants. PPG400 was an unusually strong surfactant having the largest surface activity compared with MIBC and 2-octanol, which exhibited the maximum foaming performance and foam stability at all tested gas flow rates and concentrations. The present study suggested that foam properties depended primarily on the type of surfactant and its concentration and secondarily on the gas flow rate. In addition, properties of interface are closely related to that of foam, which is a significant point if one wants to produce foams for specific applications.     

Effect of counterions on surface and foaming properties of dodecyl sulfate

The influence of counterions of surfactant on interfacial properties is studied by measuring foamability, foam stability, equilibrium and dynamic surface tension, and surface viscosity. The surfactant chosen is anionic dodecyl sulfate with various counterions, Li(+), Na(+), Cs(+), and Mg(++). Surface tension measurements show a decrease in the following order: LiDS > NaDS > CsDS > Mg(DS)(2). Foamability done using shaking method shows similar order as surface tension, i.e., LiDS > NaDS > CsDS > Mg(DS)(2). This has been explained in terms of the differences in micellar stability and diffusion of monomers. This is further confirmed by our dynamic surface tension results, which show the same order as equilibrium surface tension (i.e., LiDS > NaDS > CsDS > Mg(DS)(2)) at low bubble frequencies but the order is LiDS > NaDS = Mg(DS)(2) > CsDS at high bubble frequencies. Foam stability measurements were done at concentrations below and above cmc to elucidate the role of micelles. It was found that there is no significant change in foam stability when counterions are changed for surfactant concentration values below the cmc, but at concentration above cmc the foam stability of CsDS and Mg(DS)(2) are much greater than LiDS and NaDS indicating presence of stable micelles are essential to high foam stabilities. Surface viscosity measurements correlated well with the foam stability trends and gave the following order LiDS < NaDS < CsDS < Mg(DS)(2), indicating that the molecules of CsDS and Mg(DS)(2) are tightly packed at the air/water interface.     

Dynamic surface tension of micellar solutions studied by the maximum bubble pressure method

A theoretical model for the dynamic surface tension of an air bubble expanding in micellar surfactant solution is proposed. The model accounts for the effect of expansion of the bubble surface during the adsorption of surfactant molecules (monomers) and the effect of disintegration of polydisperse micelles on the surfactant diffusion. Assuming small deviations from equilibrium and constant rate of expansion analytical expression for the surface tension and the subsurface concentration of monomers as a function of time is derived. The characteristic time of micellization is computed from the experimental data for two surfactants (sodium dodecyl sulfate and nonylphenol polyglycol ether) obtained by the maximum bubble pressure method.     

Counterion condensation and release in micellar solutions

Counterion condensation and release in micellar solutions are investigated by direct measurement of counterion concentration with ion-selective electrode. Monte Carlo simulations based on the cell model are also performed to analyze the experimental results. The degree of counterion condensation is indicated by the concentration ratio of counterions in the bulk to the total ionic surfactant added, alpha< or =1. The ionic surfactant is completely dissociated below the critical micelle concentration (cmc). However, as cmc is exceeded, the free counterion ratio alpha declines with increasing the surfactant concentration and approaches an asymptotic value owing to counterion condensation to the surface of the highly charged micelles. Micelle formation leads to much stronger electrostatic attraction between the counterion and the highly charged sphere in comparison to the attraction of single surfactant ion with its counterion. A simple model is developed to obtain the true degree of ionization, which agrees with our Monte Carlo results. Upon addition of neutral polymer or monovalent salts, some of the surfactant counterions are released to the bulk. The former is due to the decrease of the intrinsic charge (smaller aggregation number) and the degree of ionization is increased. The latter is attributed to competitive counterion condensation, which follows the Hefmeister series. This consequence indicates that the specific ion effect plays an important role next to the electrostatic attraction.     

Surfactant accumulation within the top foam layer due to rupture of external foam films

The analysis of processes taking place in a steady pneumatic (dynamic) foam shows the possibility of different modes of surfactant accumulation within the top layers of bubbles due to rupture of external foam films. An increasing surfactant concentration within the top layers promotes the stabilisation of bubbles and the foam as a whole. Considering the balance of surfactant and water during the bursting of films it is possible to estimate the accumulated surfactant loss caused by a downwards flow through the Plateau borders of the subsurface bubble layer. This effect depends on the particular conditions, especially on the surfactant activity and concentration of the surfactant, water volume fraction in the foam and size of foam bubbles. The process of surfactant accumulation in the top foam bubble layer can be complicated due to the removal of part of the accumulated surfactant through transport with droplets spread out during bubble bursting.     

Counterion binding on mixed anionic/cationic micelles

Measurements of counterion binding in mixtures of surfactant aqueous solutions have been performed to study the structure of the anionic/cationic mixed micelle/solution interface. The mixtures studied were SDS/DDAC and STS/TDPC. The binding of chloride and sodium ions to mixed anionic/cationic micelles was measured using ion-specific electrodes. Counterion binding was found to be strongly dependent on the molar ratio of surfactants present. The mixed micelle/solution interface includes the headgroups of both surfactants and counterions of surfactant in excess. The addition of oppositely charged surfactant caused an increasing dissociation of counterions.     

Control of Ostwald ripening by using surfactants with high surface modulus

We describe results from systematic measurements of the rate of bubble Ostwald ripening in foams with air volume fraction of 90%. Several surfactant systems, with high and low surface modulus, were used to clarify the effect of the surfactant adsorption layer on the gas permeability across the foam films. In one series of experiments, glycerol was added to the foaming solutions to clarify how changes in the composition of the aqueous phase affect the rate of bubble coarsening. The experimental results are interpreted by a new theoretical model, which allowed us to determine the overall gas permeability of the foam films in the systems studied, and to decompose the film permeability into contributions coming from the surfactant adsorption layers and from the aqueous core of the films. For verification of the theoretical model, the gas permeability determined from the experiments with bulk foams are compared with values, determined in an independent set of measurements with the diminishing bubble method (single bubble attached at large air-water interface) and reasonably good agreement between the results obtained by the two methods is found. The analysis of the experimental data showed that the rate of bubble Ostwald ripening in the studied foams depends on (1) type of used surfactant-surfactants with high surface modulus lead to much slower rate of Ostwald ripening, which is explained by the reduced gas permeability of the adsorption layers in these systems; (2) presence of glycerol which reduces the gas solubility and diffusivity in the aqueous core of the foam film (without affecting the permeability of the adsorption layers), thus also leading to slower Ostwald ripening. Direct measurements showed that the foam films in the studied systems had very similar thicknesses, thus ruling out the possible explanation that the observed differences in the Ostwald ripening are due to different film thicknesses. Experiments with the Langmuir trough were used to demonstrate that the possible differences in the surface tensions of the shrinking and expanding bubbles in a given foam are too small to strongly affect the rate of Ostwald ripening in the specific systems studied here, despite the fact that some of the surfactant solutions have rather high surface modulus. The main reason for the latter observation is that the rate of surface deformation of the coarsening bubbles is extremely low, on the order of 10(-4) s(-1), so that the relaxation of the surface tension (though also slow for the high surface modulus systems) is still able to reduce the surface tension variations down to several mN/m. Thus, we conclude that the main reason for the reduced rate of bubble Ostwald ripening in the systems with high surface modulus is the low solubility and diffusivity of the gas molecules in the respective condensed adsorption layers (which have solid rather than fluid molecular packing).     

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.     

Styrène solubilization in micelles of dodecyltrimethylammonium hydroxide

 The solubilization of styrene molecules in aqueous dodeciltrimethylammonium Hydroxide (DTAOH) solution was studied by UV-Vis spectroscopy. In short, fully ionized DTAOH aggregates the styrene molecules in the micelle double layer, oriented with their vinyl group to the micelle core and the aryl ring to the intermicellar solution. At increased surfactant concentration, when the aggregates incorporate counterions in their Stern layer, the orientation is maintained, but styrene molecules gradually penetrate into the micelle core as the micelle size and degree of counterion union increased, until they were completely immersed in the hydrocarbon core of rod-like micelles. Received: 6 November 1996 Accepted: 10 February 1997     

Small-angle neutron scattering of mixed ionic perfluoropolyether micellar solutions

Aqueous mixed micellar solutions of perfluoropolyether carboxylic salts with ammonium counterions have been studied by small-angle neutron scattering. Two surfactants differing in the tail length were mixed in proportions n2/n3 = 60/40 w/w, where n2 and n3 are the surfactants with two and three perfluoroisopropoxy units in the tail, respectively. The tails are chlorine-terminated. The mixed micellar solutions, in the concentration range 0.1-0.2 M and thermal interval 20-40 degrees C, show structural characteristics of the interfacial shell that are very similar to ammonium n2 micellar solutions previously investigated; thus, the physics of the interfacial region is dominated by the polar head and counterion. The shape and dimensions of the micelles are influenced by the presence of the n3 surfactant, whose chain length in the micelle is 2 A longer than that of the n2 surfactant. The n3 surfactant favors the ellipsoidal shape in the concentration range 0.1-0.2 M with a 1/2 ionization degree of n2 micelles. The very low surface charge of the mixed micelles is attributed to the increase in hydrophobic interactions between the surfactant tails, due to the longer n3 surfactant molecules in micelles. The closer packing of the tails decreases the micellar curvature and the repulsions between the polar heads, by surface charge neutralization of counterions migrating from the Gouy-Chapman diffuse layer, leading to micellar growth in ellipsoids with greater axial ratios.     

Phase behavior in mixtures of cationic and anionic surfactants in aqueous solutions

Phase behavior of cationic/anionic surfactant mixtures of the same chain length (n=10, 12 or 14) strongly depends on the molar ratio and actual concentration of the surfactants. Precipitation of catanionic surfactant and mixed micelles formation are observed over the concentration range investigated. Coacervate and liquid crystals are found to coexist in the transition region from crystalline catanionic surfactant to mixed micelles.The addition of oppositely charged surfactant diminishes the surface charge density at the mixed micelle/solution interface and enhances the apparent degree of counterion dissociation from mixed micelles. Cationic surfactants have a greater tendency to be incorporated in mixed micelles than anionic ones.     

Thermodynamic parameters and counterion binding to the micelle in binary anionic surfactant systems

Competitive counterion binding of sodium and calcium to micelles, and mixed micellization have been investigated in the systems sodium dodecylsulfate (NaDS)/sodium decylsulfate (NaDeS) and NaDS/sodium 4-octylbenzenesulfonate (NaOBS) in order to accurately model the activity of the relevant species in solution. The critical micelle concentration (CMC) and equilibrium micelle compositions of mixtures of these anionic surfactants, which is necessary for determining fractional counterion binding measurements, is thermodynamically modeled by regular solution theory. The mixed micelle is ideal (the regular solution parameter β(M)=0) for the NaDS/NaOBS system, while the mixed micelle for NaDS/NaDeS has β(M)=-1.05 indicating a slight synergistic interaction. Counterion binding of sodium to the micelle is influenced by the calcium ion concentration, and vice versa. However, the total degree of counterion binding is essentially constant at approximately 0.65 charge negation at the micelle's surface. The counterion binding coefficients can be quantitatively modeled using a simple equilibrium model relating concentrations of bound and unbound counterions.     

Aggregation and adsorption of sodium dioctylsulfosuccinate in aqueous ammonium chloride solution: role of mixed counterions

The critical micelle concentration (cmc) of sodium dioctylsulfosuccinate (AOT) was determined at 25 °C from surface tension and fluorescence methods in aqueous NH(4)Cl solution for assessing the influence of mixed counterions on the special counterion binding behavior (SCB) of AOT. The SCB of AOT refers to a sudden twofold increase in the value of the counterion binding constant (β) in aqueous medium when the concentration (c(*)) of the added 1:1 sodium salt is about 0.015 mol kg(-1), and it has been tested so far for sodium ion only. In the presence of sodium and ammonium mixed counterions also the SCB of AOT exist, but with lower c(*) (0.009 mol kg(-1) NH(4)Cl). Synergism in the cmc occurs due to mixed counterions. In the case of inorganic counterions, unlike the case with organic counterions, the cmc is dependent on the total counterion concentration in solution and negligibly on the specific type of counterion. Na(+) and NH(4)(+) bind almost equally to the micelle in the region of low β (below c(*)), but in the region of high β (above c(*)) NH(4)(+) binds predominantly. It has been shown that the theoretical expression for the surface excess of ionic surfactant+electrolyte system containing a single counterion can also be used to evaluate the surface excess in the presence of mixed counterions if the two counterions are considered to undergo Henry-type adsorption at the air-solution interface.     

Presence or absence of counterion specificity in the interaction of alkylammonium surfactants with alkylacrylamide gels

Patterns in the interaction of cationic surfactants with nonionic polymer gels, which were inferred from a recent study from our laboratory, are confirmed by measurements of a series of alkylammonium surfactants with different counterions with a series of alkyl acrylamide gels of increasing hydrophobicity. Two swelling patterns were observed: Either the swelling continued above the surfactant critical micelle concentration (cmc) and the maximum swelling differed for different counterions and increased in the order of Br-相似文献