Quantifying the hydrophobic effect. 2. A computer simulation-molecular-thermodynamic model for the micellization of nonionic surfactants in aqueous solution

In this article, the validity and accuracy of the CS-MT model is evaluated by using it to model the micellization behavior of seven nonionic surfactants in aqueous solution. Detailed information about the changes in hydration that occur upon the self-assembly of the surfactants into micelles was obtained through molecular dynamics simulation and subsequently used to compute the hydrophobic driving force for micelle formation. This information has also been used to test, for the first time, approximations made in traditional molecular-thermodynamic modeling. In the CS-MT model, two separate free-energy contributions to the hydrophobic driving force are computed. The first contribution, gdehydr, is the free-energy change associated with the dehydration of each surfactant group upon micelle formation. The second contribution, ghydr, is the change in the hydration free energy of each surfactant group upon micelle formation. To enable the straightforward estimation of gdehydr and ghydr in the case of nonionic surfactants, a number of simplifying approximations were made. Although the CS-MT model can be used to predict a variety of micellar solution properties including the micelle shape, size, and composition, the critical micelle concentration (CMC) was selected for prediction and comparison with experimental CMC data because it depends exponentially on the free energy of micelle formation, and as such, it provides a stringent quantitative test with which to evaluate the predictive accuracy of the CS-MT model. Reasonable agreement between the CMCs predicted by the CS-MT model and the experimental CMCs was obtained for octyl glucoside (OG), dodecyl maltoside (DM), octyl sulfinyl ethanol (OSE), decyl methyl sulfoxide (C10SO), decyl dimethyl phosphine oxide (C10PO), and decanoyl-n-methylglucamide (MEGA-10). For five of these surfactants, the CMCs predicted using the CS-MT model were closer to the experimental CMCs than the CMCs predicted using the traditional molecular-thermodynamic (MT) model. In addition, CMCs predicted for mixtures of C10PO and C10SO using the CS-MT model were significantly closer to the experimental CMCs than those predicted using the traditional MT model. For dodecyl octa(ethylene oxide) (C12E8), the CMC predicted by the CS-MT model was not in good agreement with the experimental CMC and with the CMC predicted by the traditional MT model, because the simplifying approximations made to estimate gdehydr and ghydr in this case were not sufficiently accurate. Consequently, we recommend that these simplifying approximations only be used for nonionic surfactants possessing relatively small, non-polymeric heads. For MEGA-10, which is the most structurally complex of the seven nonionic surfactants modeled, the CMC predicted by the CS-MT model (6.55 mM) was found to be in much closer agreement with the experimental CMC (5 mM) than the CMC predicted by the traditional MT model (43.3 mM). Our results suggest that, for complex, small-head nonionic surfactants where it is difficult to accurately quantify the hydrophobic driving force for micelle formation using the traditional MT modeling approach, the CS-MT model is capable of making reasonable predictions of aqueous micellization behavior.     

Quantifying the hydrophobic effect. 3. A computer simulation-molecular-thermodynamic model for the micellization of ionic and zwitterionic surfactants in aqueous solution

In this article, the validity and accuracy of the CS-MT model introduced in article 1 for oil aggregates and in article 2 for nonionic surfactants is further evaluated by using it to model the micellization behavior of ionic and zwitterionic surfactants in aqueous solution. In the CS-MT model, two separate free-energy contributions to the hydrophobic driving force for micelle formation are computed using hydration data obtained from computer simulation: gdehydr, the free-energy change associated with dehydration, and ghydr, the change in the hydration free energy. To enable straightforward estimation of gdehydr and ghydr for ionic and zwitterionic surfactants, a number of simplifying approximations were made. Reasonable agreement between the CMCs predicted using the CS-MT model and the experimental CMCs was obtained for sodium dodecyl sulfate (SDS), dodecylphophocholine (DPC), cetyltrimethylammonium bromide (CTAB), two 3-hydroxy sulfonate surfactants (AOS-12 and AOS-16), and a homologous series of four DCNA bromide surfactants with a dimethylammonium head attached to a dodecyl alkyl tail and to an alkyl side chain of length CN, having the chemical formula C12H25CNH2N+1N(CH3)2Br, with N = 1 (DC1AB), 2 (DC2AB), 4 (DC4AB), and 6 (DC6AB). For six of these nine surfactants, the CMCs predicted using the CS-MT model are closer to the experimental CMCs than the CMCs predicted using the traditional molecular-thermodynamic (MT) model. For DC2AB, DC4AB, and DC6AB, which are the most structurally complex of the ionic surfactants modeled, the CMCs predicted using the CS-MT model are in remarkably good agreement with the experimental CMCs, and the CMCs predicted using the traditional MT model are quite inaccurate. Our results suggest that the CS-MT model accurately quantifies the hydrophobic driving force for micelle formation for ionic and zwitterionic surfactants in aqueous solution. For complex ionic and zwitterionic surfactants where it is difficult to accurately quantify the hydrophobic driving force for micelle formation using the traditional MT modeling approach, the CS-MT model represents a very promising alternative.     

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.     

Micellar solubilization thermodynamics of acetophenone in surfactant mixtures: Case of benzyldimethyltetradecylammonium chloride and trimethyltetradecylammonium chloride

The thermodynamics of micellar solubilization of acetophenone in mixtures of two cationic surfactants [benzyldimethyltetradecylammonium chloride +trimethyltetradecylammonium chloride] has been derived from calorimetric measurements at controlled solute activity. The partition coefficient between micelles and water as well as the standard enthalpy and entropy of transfer between micelles and water were calculated. The results were compared to the case of benzylalcohol in the same cationic mixtures. For acetophenone, the variation of all thermodynamic transfer functions with micellar composition may be described by the regular solution formalism. The same conclusion has been achieved for most polar solutes in various surfactant mixtures: favorable interaction between unlike surfactants induces an unfavorable micellar solubilization. Exceptions should be found with the cases where solute solubilization induces profound micellar changes. It seems to be the case with some alcohols in the cationic surfactant mixtures studied.     

Prediction of conformational characteristics and micellar solution properties of fluorocarbon surfactants

A molecular-thermodynamic theory is developed to model the micellization of fluorocarbon surfactants in aqueous solutions, by combining a molecular model that evaluates the free energy of micellization of fluorocarbon surfactant micelles with a previously developed thermodynamic framework describing the free energy of the micellar solution. In the molecular model of micellization developed, a single-chain mean-field theory is combined with an appropriate rotational isomeric state model of fluorocarbon chains to describe the packing of the fluorocarbon surfactant tails inside the micelle core. Utilizing this single-chain mean-field theory, the packing free energies of fluorocarbon surfactants are evaluated and compared with those of their hydrocarbon analogues. We find that the greater rigidity of the fluorocarbon chain promotes its packing in micellar aggregates of low curvatures, such as bilayers. In addition, the mean-field approach is utilized to predict the average conformational characteristics (specifically, the bond order parameters) of fluorocarbon and hydrocarbon surfactant tails within the micelle core, and the predictions are found to agree well with the available experimental results. The electrostatic effects in fluorocarbon ionic surfactant micelles are modeled by allowing for counterion binding onto the charged micelle surface, which accounts explicitly for the effect of the counterion type on the micellar solution properties. In addition, a theoretical formulation is developed to evaluate the free energy of micellization and the size distribution of finite disklike micelles, which often form in the case of fluorocarbon surfactants. We find that, compared to their hydrocarbon analogues, fluorocarbon surfactants exhibit a greater tendency to form cylindrical or disklike micelles, as a result of their larger molecular volume as well as due to the greater conformational rigidity of the fluorocarbon tails. The molecular-thermodynamic theory developed is then applied to several ionic fluorocarbon surfactant-electrolyte systems, including perfluoroalkanoates and perfluorosulfonates with added LiCl or NH(4)Cl, and various micellar solution properties, including critical micelle concentrations (cmc's), optimal micelle shapes, and average micelle aggregation numbers, are predicted. The predicted micellar solution properties agree reasonably well with the available experimental results.     

Organization of amphiphiles: XII. Evidence in favor of formation of hydrophobic complexes in aqueous solution

The effects of nonionic surfactants OP-10 and OP-30 (polyoxyethylated octyl phenols with 10 and 30 oxyethylene groups, respectively) in surfactant mixtures with ionic surfactants hexadecyltrimethylammonium bromide (CTAB) and sodium dodecyl sulphate (SDS) have been investigated by a conductometric method in conjunction with fluorescence, surface tension, zeta potential, and DLS measurements. The interactions are found to be antagonistic in nature for each of the systems; i.e., micellization of CTAB as well as SDS is hindered on addition of the nonionic surfactants. The antagonism is found to be more prominent in the presence of OP-10 compared to that of OP-30. Two types of mechanistic paths, path A operating below the critical micellar concentration and path B operating beyond the critical micellar concentration of nonionic surfactants, have been suggested. In path A, the retardation in micellization has been attributed to a decrease in monomeric concentration of the ionic surfactants from solution as a result of the formation of a hydrophobic complex between nonionic and ionic surfactants. In path B, the decrease in monomer concentration is due to the solubilization of the ionic surfactant in micelles of the nonionic surfactants in a 1:1 stoichiometric ratio. A theoretical treatment to the interaction in each ionic-nonionic pair yields a positive value of the interaction parameter supporting the concept of antagonism. The formation of the hydrophobic complex is supported by fluorescence and surface tension measurements. A schematic representation of the stabilization of these hydrophobic complexes has been suggested. The association of ionic surfactants by nonionic micelles is suggested by zeta potential and DLS studies.     

Dimeric and oligomeric surfactants. Behavior at interfaces and in aqueous solution: a review

Dimeric and oligomeric surfactants are novel surfactants that are presently attracting considerable interest in the academic and industrial communities working on surfactants. This paper first presents a number of chemical structures that have been reported for ionic, amphoteric and nonionic dimeric and oligomeric surfactants. The following aspects of these surfactants are then successively reviewed the state of dimeric and oligomeric surfactants in aqueous solutions at concentration below the critical micellization concentration (cmc); their behavior at the air/solution and solid/solution interfaces; their solubility in water, cmc and thermodynamics of micellization; the properties of the aqueous micelles of dimeric and oligomeric surfactants (ionization degree, size, shape, micropolarity and microviscosity, solution microstructure, solution rheology, micelle dynamics, micellar solubilization, interaction between dimeric surfactants and water-soluble polymers); the mixed micellization of dimeric surfactants with various conventional surfactants; the phase behavior of dimeric surfactants and the applications of these novel surfactants.     

Capillary electrophoresis study on the micellization and critical micelle concentration of sodium dodecyl sulfate. Influence of solubilized solutes

The influence of solubilized solutes on the micellization and critical micelle concentration (CMC) of sodium dodecyl sulfate (SDS) were investigated by means of capillary electrophoresis (CE). Three different structural types of test solutes, including chloropyridines. chlorophenols and cephalosporins with different binding strength to SDS micelles, were selected in this study. The variations of the effective electrophoretic mobility of these solutes as a function of SDS concentration in the premicellar and micellar regions were analyzed. Interestingly, the results indicate that, in the presence of these solubilized solutes, the micellization of SDS may occur over a range of SDS concentration, with the aggregate size increasing over this range. Depending on the nature of solubilized solutes and the extent of the interactions between solubilized solutes and SDS micelles, the CMC value of SDS may vary significantly. The incorporation of solubilized solutes into SDS micelles to form mixed micelles is proposed to interpret the migration behavior of solubilized solutes in CE.     

Separation of beta-lactam antibiotics by micellar electrokinetic chromatography

The retention behaviour of beta-lactam antibiotics in micellar electrokinetic chromatography (EKC) was investigated. Sodium dodecyl sulphate (SDS) and sodium N-lauroyl-N-methyltaurate were used an anionic surfactants at concentrations of 0.05-0.3 M. It was found that the retention of ionic substances in micellar EKC is determined by the following three factors: the electrophoretic migration of the ionic substances, the interaction between the ionic substances and ionic surfactants and solubilization of the solute by the micellar phase. A difference in the retention behaviours of cationic substances was observed between the two anionic surfactants, which have different groups neighbouring the charge-bearing groups. The effect of an ion-pairing reagent was also investigated to make the effect of the micelle clearer. All test solutes were successfully separated by micellar EKC at SDS concentrations above 0.1 M, with theoretical plate numbers ranging from 70,000 to 260,000.     

Micellization,mixed micellization and solubilization: The role of interfacial interactions

The importance of interfacial Interactions in governing micellization, mixed micellization, polymer-micene complexation and solubilization is examined in this review. A common thermodynamic approach is used to treat these different phenomena involving surfactant self-assembly. In all the cases, the free energy of self-assembly can be decomposed into bulk and interfacial components. The interfacial component arises from two competing contributions. One is due to the free energy of formation of the micellar core-solvent interface while the other is due to the steric and electrostatic interactions among the head groups at the micellar surface. The competition between these two contributions is shown here as determining all the fundamental features of self-assembly. Specifically, we discuss in this review the influence of interfacial interactions on the cooperativity of self-assembly, the critical micelle concentration, the size and size distribution of micelles, the transition from spherical to cylindrical micelles, the non-ideal behavior in mixed surfactants, the complexation or non-complexation of micelles with polymers, the solubilization of aliphatic and aromatic hydrocarbons and the selective and synergistic solubilization of hydrocarbon mixtures.     

Application of a chemical equilibrium model to thermodynamic functions of transfer of alcohols from water to aqueous surfactants

In ternary aqueous solutions, hydrophobic solutes such as alcohols tend to aggregate with surfactants to form mixed micelles. These systems can be studied by meas of the functions of transfer of hydrophobic solutes from water to aqueous solutions of surfactant. These thermodynamic functions often go through extrema in the critical micellar concentration (CMC) region of the surfactant. A simple model based on interactions between surfactant and hydrophobic solute monomers, on the distribution of the hydrophobic solute between water and the micelles and on the shift in the CMC induced by the hydrophobic solute, can simulate the magnitude and trends of the transfer functions using parameters which are mostly derived from the binary systems. In order to check the model more quantitatively, volumes and heat capacities of transfer of alcohols from water to aqueous solutions of a nonionic surfactant, octyldimethylamine oxide, were measured. A quantitative agreement was achieved with three adjustable parameters. Good fits are also obtained for the transfers to the ionic surfactants, octylamine hydrobromide and sodium dodecylsulfate. When the equilibrium displacement contribution is small, the distribution constants and the partial molar properties of the alcohols in the micellar phase agree well with the parameters obtained with similar models.     

Micellar solutions of ionic surfactants and their mixtures with nonionic surfactants: Theoretical modeling vs. Experiment

Here, we review two recent theoretical models in the field of ionic surfactant micelles and discuss the comparison of their predictions with experimental data. The first approach is based on the analysis of the stepwise thinning (stratification) of liquid films formed from micellar solutions. From the experimental step-wise dependence of the film thickness on time, it is possible to determine the micelle aggregation number and charge. The second approach is based on a complete system of equations (a generalized phase separation model), which describes the chemical and mechanical equilibrium of ionic micelles, including the effects of electrostatic and non-electrostatic interactions, and counterion binding. The parameters of this model can be determined by fitting a given set of experimental data, for example, the dependence of the critical micellization concentration on the salt concentration. The model is generalized to mixed solutions of ionic and nonionic surfactants. It quantitatively describes the dependencies of the critical micellization concentration on the composition of the surfactant mixture and on the electrolyte concentration, and predicts the concentrations of the monomers that are in equilibrium with the micelles, as well as the solution’s electrolytic conductivity; the micelle composition, aggregation number, ionization degree and surface electric potential. These predictions are in very good agreement with experimental data, including data from stratifying films. The model can find applications for the analysis and quantitative interpretation of the properties of various micellar solutions of ionic surfactants and mixed solutions of ionic and nonionic surfactants.     

Phenyl-ring-bearing cationic surfactants: effect of ring location on the micellar structure

A series of isomeric cationic surfactants (S1-S5) bearing a long alkyl chain that carries a 1,4-phenylene unit and a trimethyl ammonium headgroup was synthesized; the location of the phenyl ring within the alkyl tail was varied in an effort to understand its influence on the amphiphilic properties of the surfactants. The cmc's of the surfactants were estimated using ionic conductivity measurements and isothermal calorimetric titrations (ITC); the values obtained by the two methods were found to be in excellent agreement. The ITC measurements provided additional insight into the various thermodynamic parameters associated with the micellization process. Although all five surfactants have exactly the same molecular formula, their micellar properties were seen to vary dramatically depending on the location of the phenyl ring; the cmc was seen to decrease by almost an order of magnitude when the phenyl ring was moved from the tail end (cmc of S1 is 23 mM) to the headgroup region (cmc of S5 is 3 mM). In all cases, the enthalpy of micellization was negative but the entropy of micellization was positive, suggesting that in all of these systems the formation of micelles is both enthalpically and entropically favored. As expected, the decrease in cmc values upon moving the phenyl ring from the tail end to the headgroup region is accompanied by an increase in the thermodynamic driving force (ΔG) for micellization. To understand further the differences in the micellar structure of these surfactants, small-angle neutron scattering (SANS) measurements were carried out; these measurements reveal that the aggregation number of the micelles increases as the cmc decreases. This increase in the aggregation number is also accompanied by an increase in the asphericity of the micellar aggregate and a decrease in the fractional charge. Geometric packing arguments are presented to account for these changes in aggregation behavior as a function of phenyl ring location.     

Thermal signature of hydrophobic hydration dynamics

Hydrophobic hydration, the perturbation of the aqueous solvent near an apolar solute or interface, is a fundamental ingredient in many chemical and biological processes. Both bulk water and aqueous solutions of apolar solutes behave anomalously at low temperatures for reasons that are not fully understood. Here, we use (2)H NMR relaxation to characterize the rotational dynamics in hydrophobic hydration shells over a wide temperature range, extending down to 243 K. We examine four partly hydrophobic solutes: the peptides N-acetyl-glycine-N'-methylamide and N-acetyl-leucine-N'-methylamide, and the osmolytes trimethylamine N-oxide and tetramethylurea. For all four solutes, we find that water rotates with lower activation energy in the hydration shell than in bulk water below 255 +/- 2 K. At still lower temperatures, water rotation is predicted to be faster in the shell than in bulk. We rationalize this behavior in terms of the geometric constraints imposed by the solute. These findings reverse the classical "iceberg" view of hydrophobic hydration by indicating that hydrophobic hydration water is less ice-like than bulk water. Our results also challenge the "structural temperature" concept. The two investigated osmolytes have opposite effects on protein stability but have virtually the same effect on water dynamics, suggesting that they do not act indirectly via solvent perturbations. The NMR-derived picture of hydrophobic hydration dynamics differs substantially from views emerging from recent quasielastic neutron scattering and pump-probe infrared spectroscopy studies of the same solutes. We discuss the possible reasons for these discrepancies.     

Experimental and theoretical investigation of the micellar-assisted solubilization of ibuprofen in aqueous media

Surfactants can be used to increase the solubility of poorly soluble drugs in water and to increase drug bioavailability. In this article, the aqueous solubilization of the nonsteroidal, antiinflammatory drug ibuprofen is studied experimentally and theoretically in micellar solutions of anionic (sodium dodecyl sulfate, SDS), cationic (dodecyltrimethylammonium bromide, DTAB), and nonionic (dodecyl octa(ethylene oxide), C12E8) surfactants possessing the same hydrocarbon "tail" length but differing in their hydrophilic headgroups. We find that, for these three surfactants, the aqueous solubility of ibuprofen increases linearly with increasing surfactant concentration. In particular, we observed a 16-fold increase in the solubility of ibuprofen relative to that in the aqueous buffer upon the addition of 80 mM DTAB and 80 mM C12E8 but only a 5.5-fold solubility increase upon the addition of 80 mM SDS. The highest value of the molar solubilization capacity (chi) was obtained for DTAB (chi = 0.97), followed by C12E8 (chi = 0.72) and finally by SDS (chi = 0.23). A recently developed computer simulation/molecular-thermodynamic modeling approach was extended to predict theoretically the solubilization behavior of the three ibuprofen/surfactant mixtures considered. In this modeling approach, molecular-dynamics (MD) simulations were used to identify which portions of ibuprofen are exposed to water (hydrated) in a micellar environment by simulating a single ibuprofen molecule at an oil/water interface (modeling the micelle core/water interface). On the basis of this input, molecular-thermodynamic modeling was then implemented to predict (i) the micellar composition as a function of surfactant concentration, (ii) the aqueous solubility of ibuprofen as a function of surfactant concentration, and (iii) the molar solubilization capacity (chi). Our theoretical results on the solubility of ibuprofen in aqueous SDS and C12E8 surfactant solutions are in good agreement with the experimental data. The ibuprofen solubility in aqueous DTAB solutions was somewhat overpredicted because of challenges associated with accurately modeling the strong electrostatic interactions between the anionic ibuprofen and the cationic DTAB. Our results indicate that computer simulations of ibuprofen at a flat oil/water interface can be used to obtain accurate information about the hydrated and the unhydrated portions of ibuprofen in a micellar environment. This information can then be used as input to a molecular-thermodynamic model of self-assembly to successfully predict the aqueous solubilization behavior of ibuprofen in the three surfactant systems studied.     

On the mechanism of hydrophobic association of nanoscopic solutes

The hydration behavior of two planar nanoscopic hydrophobic solutes in liquid water at normal temperature and pressure is investigated by calculating the potential of mean force between them at constant pressure as a function of the solute-solvent interaction potential. The importance of the effect of weak attractive interactions between the solute atoms and the solvent on the hydration behavior is clearly demonstrated. We focus on the underlying mechanism behind the contrasting results obtained in various recent experimental and computational studies on water near hydrophobic solutes. The length scale where crossover from a solvent separated state to the contact pair state occurs is shown to depend on the solute sizes as well as on details of the solute-solvent interaction. We find the mechanism for attractive mean forces between the plates is very different depending on the nature of the solute-solvent interaction which has implications for the mechanism of the hydrophobic effect for biomolecules.     

Effect of stiffness on the micellization behavior of model H4T4 surfactant chains

The micellization behavior of a series of model surfactants, all with four head and tail groups (H4T4) but with different degrees of chain stiffness, was studied using grand canonical Monte Carlo simulations on a cubic lattice. The critical micelle concentration, micellar size, and thermodynamics of micellization were examined. In all cases investigated, the critical micelle concentration was found to increase with increasing temperature as observed for nonionic surfactants in apolar or slightly polar solvents. At a fixed reduced temperature and increasing chain stiffness, in agreement with previous observations, it was found that the critical micelle concentration decreased and the average micelle size increased. This behavior is qualitatively consistent with that experimentally observed when comparing hydrogenated and homologous fluorinated surfactants. Thermodynamic considerations based on the analysis of the temperature dependence of the critical micelle concentration indicated that both effects could be interpreted as arising from an increased number of heterocontacts between hydrophobic portions of stiff surfactants and a lower entropic cost on packing rigid chains. Structural analysis that was also based on considering the inner micellar radial dependence of the surfactant head and tail site fraction distributions suggested that, on stiffening the molecular backbone, the resulting micellar aggregates grew, without appreciable deviations from spherical symmetry. Stiffer surfactants led to a slightly denser micellar core because of better packing.