Kinetic description of the relaxation of surfactant solutions containing spherical and cylindrical micelles

Monotonically decaying relaxation of a materially isolated nonionic surfactant solution containing spherical and cylindrical micelles at the arbitrary heights of the first and second potential barriers of aggregation work is kinetically substantiated. The realistic situation, where the height of second potential barrier is at least slightly higher (by the relative value) than that of the first barrier, is studied. Analytical expressions for two relaxation times of materially isolated surfactant solution are calculated. The shortest of these times corresponds to the relatively fast establishment of the mutual quasi-equilibrium of spherical and cylindrical micelles, beginning with relatively small cylindrical micelles. The longest of relaxation times corresponds to the relatively slow establishment of the total equilibrium of surfactant solution. It is shown that this time (the only significant for the establishment of the final equilibrium of materially isolated surfactant solution) is determined by the height of the first potential barrier of aggregation work and is by no means dependent on the height of the second potential barrier about which not much is known. Variations (with time) of the total concentrations of spherical and cylindrical micelles, surfactant monomer concentration, and the total amount of the substance in cylindrical micelles in the approach of solution to the final equilibrium state are described analytically. It is shown that theoretically admitted small relative deviations of the concentrations of spherical and cylindrical micelles from their values in the final equilibrium state are fully measurable in experiment. Calculated relaxation time of surfactant solution can also be measured experimentally together with the aforementioned values. It is elucidated that this time is approximately proportional to the overall solution concentration, if the second critical micellization concentration (CMC2) by the order of magnitude exceeds the first critical micellization concentration (CMC1), and is virtually independent of the overall solution concentration, if the CMC2 exceeds the CMC1 by two orders of magnitude. The characteristic time of the establishment of quasi-equilibrium distribution of cylindrical micelles throughout the region of their sizes is estimated, thus allowing us to establish the lower limit of the height of the first barrier of aggregation work.Translated from Kolloidnyi Zhurnal, Vol. 67, No. 1, 2005, pp. 47–56.Original Russian Text Copyright © 2005 by Kuni, Shchekin, Grinin, Rusanov.     

Boltzmann distributions and slow relaxation in systems with spherical and cylindrical micelles

Equilibrium and nonequilibrium distributions of molecular aggregates in a solution of a nonionic surfactant are investigated at the total surfactant concentration above the second critical micelle concentration (CMC2). The investigation is not limited by the choice of a specific micellar model. Expressions for the direct and reverse fluxes of molecular aggregates over the potential humps of the aggregation work are derived. These aggregation work humps set up activation barriers for the formation of spherical and cylindrical micelles. With the aid of the expressions for molecular aggregate fluxes, a set of two kinetic equations of micellization is derived. This set, along with the material balance equation, describes the molecular mechanism of the slow relaxation of micellar solution above the CMC2. A realistic situation has been analyzed when the CMC2 exceeds the first critical micelle concentration, CMC1, by an order of magnitude, and the total surfactant concentration varies within the range lying markedly above the CMC2 but not by more than 2 orders of magnitude. For such conditions, an equation relating the parameters of the aggregation work of a cylindrical micelle to the observable ratio of the total surfactant concentration and the monomer concentration is found for an equilibrium solution. For the same conditions, but in the nonequilibrium state of materially isolated surfactant solution, a closed set of linearized relaxation equations for total concentrations of spherical and cylindrical micelles is derived. These equations determine the time development of two modes of slow relaxation in micellar solutions markedly above the CMC2. Solving the set of equations yields two rates and two times of slow relaxation.     

System of relaxation equations for materially isolated surfactant solution with spherical and cylindrical micelles

Analytical expressions for the direct and reverse fluxes of molecular aggregates over the first and second potential barriers of the aggregation work in the presence of spherical and cylindrical micelles in non-ionic surfactant solution were derived. Expressions for the sum (entering into kinetic equations of micellization) of direct and reverse fluxes of molecular aggregates over the first and second potential barriers of the aggregation work in the vicinity of the final equilibrium state of materially isolated surfactant solution were linearized. In the experimentally important range of the values of overall surfactant concentration in solution where the predominant contribution to the total surfactant amount is introduced by cylindrical micelles, we derived a closed system of two linearized relaxation equations determining the buildup (with time) of experimentally observed total concentrations of spherical and cylindrical micelles in the vicinity of the final equilibrium state of materially isolated surfactant solution. The case of the solutions of such surfactants, for which the spherical shape of a micelle appeared to be unrealizable due to the structure and packing conditions of molecules, was considered separately.Translated from Kolloidnyi Zhurnal, Vol. 67, No. 1, 2005, pp. 38–46. Original Russian Text Copyright © 2005 by Kuni, Shchekin, Rusanov, Grinin.     

Micellization kinetics with allowance for fussion and fission of spherical and cylindrical micelles: 1. Set of nonlinear equations describing slow relaxation

Based on the general kinetic equation that describes the aggregation and fragmentation of surfactant molecular aggregates, a closed set of nonlinear equations is derived for the slow relaxation of surfactant monomer concentration and the total concentrations of coexisting spherical and cylindrical micelles to the equilibrium state of a micellar solution. Both the transitions accompanied by the emission and capture of surfactant monomers by micelles and the transitions resulting from the fussion and fission of micelles, are taken into account. The derived set of equations describes all stages of the slow relaxation from the initial perturbance to the final equilibrium state of a micellar solution.     

Concentrations of Monomers and Cylindrical Micelles above the Second CMC

Based on thermodynamically substantiated linear dependence of the work of cylindrical micelle formation on the aggregation number within a wide range of aggregation numbers where the cylindrical micelles are accumulated in a surfactant solution, the second critical micellization concentration (CMC) is introduced as an overall surfactant concentration at which the ratio of the total amount of substance in cylindrical micelles to the amount of substance in monomers is equal to 0.1, i.e., it is already noticeable. It is shown that this ratio increases rather rapidly with a monomer concentration. The coefficient of the linear dependence of the work of cylindrical micelle formation on the aggregation number in the important practical situation where the ratios of the total concentration of cylindrical micelles and total amount of substance in these micelles to the monomer concentration are equal by the order of magnitude to 1 and 10⁵, respectively, while disc micelles and extended bilayers are still not appeared. In the same situation, the ratios of the total concentration of spherical micelles and total amount of substance in these micelles to the monomer concentration are equal by the order of magnitude to 1 and 10², respectively. The relationship between the overall surfactant concentration and monomer concentration is found. It is shown that the second CMC exceeds by two orders of magnitude the first CMC corresponding to the onset of the noticeable accumulation of surfactant in spherical micelles. The distribution of cylindrical micelles over the aggregation numbers is analyzed. It is demonstrated that, in agreement with the experiment, the distribution is almost uniform in the considerable part of the wide range of aggregation numbers and drops exponentially in the remaining (right-hand) part of this range. Experimental result is confirmed that the total concentration of cylindrical micelles, the mean value, and the mean statistical scatter of aggregation numbers in a cylindrical micelle is proportional to the square root of the overall surfactant concentration. The balance equation of surfactant amount in the vicinity of the final equilibrium state of a materially isolated solution is linearized. This linearization makes it possible to express the deviations of monomer and aggregate concentrations from their equilibrium values at the lower boundary of the region of the linear dependence of the work of cylindrical micelle formation on the aggregation numbers via the deviations of experimentally observed total concentrations of spherical and cylindrical micelles from their equilibrium values. The case of the solutions of such surfactants, for which spherical shape appeared to be unrealizable due to their molecular structure and packing conditions, is considered separately.     

Kinetics of Aggregation and Relaxation in Micellar Surfactant Solutions

Theoretical results published in the last 17 years on the kinetics of aggregation and relaxation in micellar surfactant solutions have been reviewed. The results obtained by the analytical and direct numerical solution of the Becker–Döring kinetic equations and the Smoluchowski generalized equations, which describe different possible mechanisms of aggregation and relaxation on all time scales from ultrafast relaxation while reaching the quasi-equilibrium in the region of subcritical molecular aggregates to the last stage of slow relaxation of micelles to the final aggregated state, have been considered in detail. The droplet model and the model linear with respect to aggregation numbers have been used for the work of aggregation to describe the dynamics of the rearrangement of micellar systems consisting of only spherical, only cylindrical, and coexisting spherical and cylindrical aggregates, with the dynamics being both linear and nonlinear with respect to deviations from equilibrium. The results of molecular simulation of the rearrangement kinetics of micellar systems subjected to initial disturbance have been reviewed.     

Thermodynamic and Kinetic Foundations of the Micellization Theory: 4. Kinetics of Establishment of Equilibrium in a Micellar Solution

A system of the kinetic equations of the material balance for the concentrations of surfactant monomers and micelles in a micellar nonionic surfactant solution was formulated. The equilibrium state of a materially isolated micellar solution was analyzed. The system of the kinetic equations of the material balance of a micellar solution was solved. The total time of the establishment of equilibrium in a micellar solution was determined. It was shown that this time increases or (typically) decreases with an increase in micelle concentration, depending on the degree of micellization.     

Thermodynamic and Kinetic Foundations of the Micellization Theory: 5. Hierarchy of Kinetic Times

The characteristic kinetic times of micellization in the solution of a nonionic surfactant: the times of establishment of quasi-equilibrium concentrations of molecular aggregates in micellar, subcritical, and overcritical regions, times of establishment of quasi-equilibrium concentrations of molecular aggregates in the near-critical region of their sizes, the average time between two successive acts of emission of surfactant monomers by a micelle, the average value of micelle lifetime, the time of establishment of quasi-stationary mode of matter exchange between the solution and molecular aggregate, as well as the times of fast and slow relaxation in a solution were analyzed. The hierarchy of these times disclosing complex multistage kinetic process of micelle formation and decomposition and the establishment of equilibrium in the micellar solution was revealed. It was shown that this hierarchy is provided by the small parameters of the kinetic theory. The inverse problem of micellization kinetics was discussed; this problem allows us to find the characteristics of the formation work for micellar aggregate from the experimental data on the relaxation time of micellar solution.     

Power-law stage of slow relaxation in solutions with spherical micelles

General (independent of models selected for surfactant molecular aggregates) analytical relations are derived to describe the initial stage of slow relaxation in micellar solutions with spherical micelles. This stage precedes the final stage of the relaxation occurring via an exponential decay of disturbances with time. The relations obtained are applicable throughout the interval of micellar solution concentrations from the first to the second critical micellization concentration. It is shown that the initial stage is characterized by power laws of variations in the concentrations of monomers and micelles with time, these laws being different for the relaxation processes proceeding from above and below toward equilibrium values of micellar solution parameters. Relations are derived for the duration of this stage, and the effect of initial conditions is studied. Characteristic times of the power-law stage are determined and compared with the characteristic time of the final exponent-law relaxation stage. The behavior of these times is investigated at surfactant solution concentrations in the vicinity of, and noticeably above, the first critical micellization concentration. On the basis of the droplet and quasi-droplet thermodynamic models of surfactant molecular aggregates, numerical solutions are found for nonlinearized equations of slow relaxation for the time dependence of surfactant monomer concentrations at all stages of the slow relaxation. Numerical results obtained from the models are compared with the results of a general analytical study.     

Viscometric detection of sphere to cylinder transition and polydispersity in aqueous micellar solutions

Viscosity measurements of aqueous surfactant solutions using a Cannon-Fenske capillary viscometer have been made to examine the conditions under which spherical micelles grow into cylindrical micelles. Surfactants with different polar head groups and hydrocarbon chain lengths have been studied at various solution conditions. The observed transition from spherical to cylindrical micelles is explained in terms of the attractive and repulsive forces associated with the micellization process. Further the viscosity of SDS micellar solutions have been coraputed assuming the size distribution data generated by light scattering measurements. It is found that the relative viscosity calculated for a polydisperse solution is close to that calculated assuming monodispersed, rigid rods having size equal to the weight average aggregation numbers. Also the calculated viscosities compare well with the experimental viscometric data; thus demonstrating the compatibility of micellar size; determined from light scattering and viscometric measurements. Finally, the relative roles of solution polydispersity and micelle flexibility in the interpretation of experimental viscometric data are evaluated.     

Structure and dynamics of concentrated micellar solutions of sodium dodecyl sulfate

In micellar solutions of sodium dodecyl sulfate, as the concentration of surfactants increases, the spheroid shape of the micelles changes from almost spherical to ellipsoidal with increasing ratio of half-axes ratio, and further the transition to cylindrical micelles occurs. The micelles in an aqueous solution can directly contact (compact aggregates) or be separated from one another by layers of intermicellar medium (periodical colloid structures). In the latter case, the thickness of the layer can significantly exceed the micelle size, and then no mutual correlation in micelle arrangement is observed. According to the data of small-angle X-ray scattering, the relationship between the surfactant concentration and formation of “quasi-crystalline” micellar structure is nonlinear, which can be due to both micelle aggregation processes and nonuniformity of their structure. The possible influence of ordered micellar structures on the diffusion mobility of micelles is shown.     

The “fine structure” of the slow micellar relaxation mode and the aggregation rates in the range between a potential hump and well in the work of aggregation

An analytical expression has been derived for the quasi-stationary size distribution of surfactant aggregates in a micellar system approaching the final equilibrium state. In contrast to previously known relations, the derived expression takes into account variations in the concentration of monomers during the slow relaxation and enables one to determine the previously unknown fine structure of the linearized mode of slow relaxation, i.e., its dependence on the aggregation numbers in the range between the maximum and minimum of the work of aggregation. This dependence has been reliably confirmed by the numerical solution of the set of linearized Becker–Döering difference equations, which describe the molecular mechanism of the kinetics of micellization and micellar relaxation. In turn, the expression found for the relaxation mode makes it possible to refine the analogous “fine structure” of aggregation rates at different points of the same range between the maximum and minimum of the work of aggregation, in which the aggregation rates appear to be low but exhibit a nonmonotonic behavior. This behavior is also confirmed by the numerical solution of the Becker–Döering difference kinetic equations.     

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.     

Length control in rigid cylindrical nanoassembly by tuning molecular interactions in aqueous solutions

We have studied aqueous micellar solutions of nonionic surfactant (pentaethylene glycol mono-n-dodecyl ether, C12E5) doped by cationic surfactant (dodecyl trimethylamoniumbromide, DTAB) as a function of doping level, using small angle neutron scattering. At a doping level of at least 6 mol %, rigid cylindrical micelles formed and the local cylindrical structure of the doped micelles showed no variation across the range of doping levels covered in this study (0-10 mol %). However, the total micellar length decreased rapidly as doping level increased, following well the prediction of micellar aggregation number based on molecular-thermodynamic theory. There was no synergistic interaction between surfactants, leading to monotonically decreasing the micellar aggregation number (shortening of the micellar length).     

Toward a Theory of Diffusion of a Nonionic Surfactant with Variable Aggregation Number in a Micellar System

Since the aggregation number of micelles always grows with concentration, and, in some cases this dependence is noticeable even for spherical micelles, there is a need to revise the theory of micellization, in which the aggregation number is assumed to be constant. This work reformulates the theory of diffusion of nonionic surfactants in micellar solutions with regard to the variability of the aggregation number. A new formula, which expresses the diffusion coefficient of a surfactant via the diffusion coefficients of monomers and micelles, contains an additional factor capable of increasing the diffusion coefficient with the surfactant concentration. However, this factor is not overly strong, and the “old” part of the formula acts in the opposite direction; as a result, the conventional decrease in the diffusion coefficient of a nonionic surfactant remains prevailing. The analytical consideration has been supplemented with numerical calculations, the results of which are presented in the tables.