Determination of kinetic parameters of probe migration in micelles using simultaneous analysis of the fluorescence decay surface

Monte Carlo simulation of mixed lennard-jones nonionic surfactant adsorption at the liquid/vapor interface

New Monte Carlo simulations are presented for nonionic surfactant adsorption at the liquid/vapor interface of a monatomic solvent specifically investigating the roles of tail attraction and binary mixtures of different tail lengths. Surfactant molecules consist of an amphiphilic chain with a solvophilic head and a solvophobic tail. All molecules in the system, solvent and surfactant, are characterized by the Lennard-Jones (LJ) potential. Adjacent atoms along the surfactant chain are connected by finitely extensible harmonic springs. Solvent molecules move via the Metropolis random-walk algorithm, whereas surfactant molecules move according to the continuum configurational bias Monte Carlo (CBMC) method. We generate thermodynamic adsorption and surface-tension isotherms and compare results quantitatively to single-surfactant adsorption (Langmuir, 2007, 23, 1835). Surfactant tail groups with attractive interaction lead to cooperative adsorption at high surface coverage and higher maximum adsorption at the interface than those without. Moreover, adsorption and surface-tension isotherms with and without tail attraction are identical at low concentrations, deviating only near maximum coverage. Simulated binary mixtures of surfactants with differing lengths give intermediate behavior between that of the corresponding single-surfactant adsorption and surface-tension isotherms both with and without tail attraction. We successfully predict simulated mixture results with the thermodynamically consistent ideal adsorbed solution (IAS) theory for binary mixtures of unequal-sized surfactants using only the simulations from the single surfactants. Ultimately, we establish that a coarse-grained LJ surfactant system is useful for understanding actual surfactant systems when tail attraction is important and for unequal-sized mixtures of amphiphiles.

Compartmental analysis in photophysics: Kinetics and identifiability of models for quenching of fluorescent probes in micelles

The parameters describing the kinetics of excited-state processes can possibly be recovered by analysis of the fluorescence decay surface measured as a function of the experimental variables. The identifiability analysis of a photophysical model assuming errorless time-resolved fluorescence data can verify whether the model parameters can be determined. In this work, we have used the methods of similarity transformation and Taylor series to investigate the identifiability of two models utilized to describe the time-resolved fluorescence quenching of stationary probes in micelles. The first model assumes that exchange of the quencher between micelles is much slower than the fluorescence decay of the unquenched probe (the 'immobile' quencher model). The second model assumes that quenchers exchange between the aqueous and micellar phases (the 'mobile' quencher model). For the 'immobile' quencher model, the rate constants for deactivation (k(0)) and quenching (k(q)) of the excited probe are uniquely identified together with the average number of quencher molecules per micelle. For the 'mobile' quencher model, the rate constants k(0) and k(q) are uniquely identified, as are the rate constants for entry (k(+)) and exit (k(-)) of one quencher molecule into and from a micelle, and the micellar aggregation number. The concomitant rate equations describing the time-resolved fluorescence are solved using z-transforms.

Effect of surfactant structure on the residual fluorescence of micelle-based fluorescent probes

In the present paper we have investigated some photo-physical characteristics of different micellar-based fluorescent probes containing a fluorophore (pyrene) and a quencher unit (dodecyl-dioxo 2,3,2). The fluorescent response of the probe in the presence of Cu(II) ions was studied using different micellar substrates, and it was found that the pH at which the On-Off jump occurs is not influenced by the chemical structure of surfactant. In addition, the experimental residual fluorescence is not proportionally affected by microviscosity or by the size of the micellar aggregates. The signal of the native fluorescence of pyrene was observed even when the quencher's occupancy number was greater than one. Moreover, we observed discrepancies between experimental values and calculated residual fluorescence using Laplace data. These results were interpreted suggesting that the residual fluorescence has two main components, one that seems to be independent on micellar properties, while the other is directly related to location of molecules inside the surfactant aggregates that serve as substrate.

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.

Quantifying the Hydrophobic Effect. 1. A Computer Simulation−Molecular-Thermodynamic Model for the Self-Assembly of Hydrophobic and Amphiphilic Solutes in Aqueous Solution

Surfactant micellization and micellar solubilization in aqueous solution can be modeled using a molecular-thermodynamic (MT) theoretical approach; however, the implementation of MT theory requires an accurate identification of the portions of solutes (surfactants and solubilizates) that are hydrated and unhydrated in the micellar state. For simple solutes, such identification is comparatively straightforward using simple rules of thumb or group-contribution methods, but for more complex solutes, the hydration states in the micellar environment are unclear. Recently, a hybrid method was reported by these authors in which hydrated and unhydrated states are identified by atomistic simulation, with the resulting information being used to make MT predictions of micellization and micellar solubilization behavior. Although this hybrid method improves the accuracy of the MT approach for complex solutes with a minimum of computational expense, the limitation remains that individual atoms are modeled as being in only one of two states-head or tail-whereas in reality, there is a continuous spectrum of hydration states between these two limits. In the case of hydrophobic or amphiphilic solutes possessing more complex chemical structures, a new modeling approach is needed to (i) obtain quantitative information about changes in hydration that occur upon aggregate formation, (ii) quantify the hydrophobic driving force for self-assembly, and (iii) make predictions of micellization and micellar solubilization behavior. This article is the first in a series of articles introducing a new computer simulation-molecular thermodynamic (CS-MT) model that accomplishes objectives (i)-(iii) and enables prediction of micellization and micellar solubilization behaviors, which are infeasible to model directly using atomistic simulation. In this article (article 1 of the series), the CS-MT model is introduced and implemented to model simple oil aggregates of various shapes and sizes, and its predictions are compared to those of the traditional MT model. The CS-MT model is formulated to allow the prediction of the free-energy change associated with aggregate formation (gform) of solute aggregates of any shape and size by performing only two computer simulations-one of the solute in bulk water and the other of the solute in an aggregate of arbitrary shape and size. For the 15 oil systems modeled in this article, the average discrepancy between the predictions of the CS-MT model and those of the traditional MT model for gform is only 1.04%. In article 2, the CS-MT modeling approach is implemented to predict the micellization behavior of nonionic surfactants; in article 3, it is used to predict the micellization behavior of ionic and zwitterionic surfactants.

Molecular-thermodynamic theory of micellization of pH-sensitive surfactants

A predictive, molecular-thermodynamic theory is developed to model the micellization of pH-sensitive surfactants. The theory combines a molecular-thermodynamic description of micellization in binary surfactant mixtures with the protonation equilibrium of the surfactant monomers. The thermodynamic component of the theory models the pH-mediated equilibrium between micelles, surfactant monomers, and counterions. These counterions may originate from the surfactant or from added salt, acid, or base. The molecular component of the theory models the various contributions to the free energy of micellization, which corresponds to the free-energy change associated with forming a mixed micelle from the protonated and deprotonated forms of the surfactant and from the bound counterions. The free energy of micellization includes hydrophobic, interfacial, packing, steric, electrostatic, and entropic contributions, which are all calculated molecularly. The theory also requires knowledge of the surfactant molecular structure and the solution conditions, including the temperature and the amount of any added salt, acid, or base. To account for the pH sensitivity of the surfactant, the theory requires knowledge of the surfactant monomer equilibrium deprotonation constant (pK1), which may be obtained from experimental titration data obtained below the critical micelle concentration (cmc). The theory can be utilized to predict the equilibrium micelle and solution properties, including the cmc, the micelle composition, the micelle shape and aggregation number, the solution pH, and the micelle deprotonation equilibrium constant (pKm). Theoretical predictions of the cmc, the micelle aggregation number, and the pKm compare favorably with the available experimental data for alkyldimethylamine oxide surfactants. This class of pH-sensitive surfactants exhibits a form of self-synergy, which has previously been attributed to hydrogen-bond formation at the micelle interface. Instead, we show that much of the observed synergy is related to the electrostatic contribution to the free energy of micellization. Although we do not explicitly include hydrogen bonding in the molecular model of micellization, we briefly discuss how it may be incorporated and its anticipated effect on the predicted micellization behavior.

Complementary use of simulations and molecular-thermodynamic theory to model micellization

Molecular-thermodynamic descriptions of micellization in aqueous media can be utilized to model the self-assembly of surfactants possessing relatively simple chemical structures, where it is possible to identify a priori what equilibrium position they will adopt in the resulting micellar aggregate. For such chemical structures, the portion of the surfactant molecule that is expected to be exposed to water upon aggregate self-assembly can be identified and used as an input to the molecular-thermodynamic description. Unfortunately, for many surfactants possessing more complex chemical structures, it is not clear a priori how they will orient themselves within a micellar aggregate. In this paper, we present a computational approach to identify what portions of a surfactant molecule are hydrated in a micellar environment through the use of molecular dynamics simulations of such molecules at an oil/water interface (modeling the micelle core/water interface). The local environment of each surfactant segment is determined by counting the number of contacts of each segment with the water and oil molecules. After identifying the hydrated and the unhydrated segments of the surfactant molecule, molecular-thermodynamic modeling can be performed to predict: (i) the free-energy change associated with forming a micellar aggregate, (ii) the critical micelle concentration (CMC), and (iii) the optimal shape and size of the micellar aggregate. The computer simulation results were found to be sensitive to the atomic charge parameters utilized during the simulation runs. Two different methods of assigning atomic charges were tested, and the computer simulation and molecular-thermodynamic modeling results obtained using both sets of atomic charges are presented and compared. The combined computer simulation/molecular-thermodynamic modeling approach presented here is validated first by implementing it in the case of anionic (sodium dodecyl sulfate, SDS), cationic (cetyltrimethylammonium bromide, CTAB), zwitterionic (dodecylphosphocholine, DPC), and nonionic (dodecyl poly(ethylene oxide), C12E8) surfactants possessing relatively simple chemical structures and verifying that good predictions of CMCs and micelle aggregation numbers are obtained. In the case of C12E8, the challenges and limitations associated with simulating a single, polymeric E8 moiety at the oil/water interface to model its behavior at the micelle/water interface are discussed. Subsequently, the combined modeling approach is implemented in the case of the anionic surfactant 3-hydroxy sulfonate (AOS) and of the nonionic surfactant decanoyl-n-methylglucamide (MEGA-10), which possess significantly more complex chemical structures. The good predictions obtained for these two surfactants indicate that the combined computer simulation/molecular-thermodynamic modeling approach presented here extends the range of applicability of molecular-thermodynamic theory to allow modeling of the micellization behavior of surfactants possessing more complex chemical structures.

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.

Identifiability of Models for Fluorescence Quenching in Aqueous Micellar Systems

The first deterministic identifiability analysis is presented for four commonly used kinetic models of fluorescence quenching of an excited probe in aqueous micelles: A) model with immobile quenchers, B) model with mobile quenchers, C) an extension of model B in which exchange of quenchers both via the aqueous phase and during micelle collisions is taken into account, and D) model with probe migration. It is shown that these specific models for fluorescence decay of an excited probe solubilized in a micelle and quenched by molecules or ions that are Poisson-distributed over the micelles, resulting in the generalized four-parameter equation f(t)=A(1) exp{-A(2)t-A(3)[1-A(4)t]}, are uniquely identifiable in terms of four descriptive A parameters. Moreover, each model also can be uniquely identified in terms of the underlying rate constants and micellar concentration or mean micellar aggregation number. This means that these parameters can be extracted in a unique way from time-resolved fluorescence quenching experiments on a probe in micelles. For each model the recommended analysis approach is given.

Modeling counterion binding in ionic-nonionic and ionic-zwitterionic binary surfactant mixtures

A predictive molecular-thermodynamic theory is developed to model the effect of counterion binding on micellar solution properties of binary surfactant mixtures of ionic and nonionic (or zwitterionic) surfactants. The theory combines a molecular-thermodynamic description of micellization in binary surfactant mixtures with a recently developed model of counterion binding to single-component ionic surfactant micelles. The thermodynamic component of the theory models the equilibrium between the surfactant monomers, the counterions, and the mixed micelles. The molecular component of the theory models the various contributions to the free-energy change associated with forming a mixed micelle from ionic surfactants, nonionic (or zwitterionic) surfactants, and bound counterions (referred to as the free energy of mixed micellization). Specifically, the various molecular contributions to the free energy of mixed micellization model the underlying physics associated with the assembly of, and the interactions between, the surfactant polar heads, the surfactant nonpolar tails, and the bound counterions. Utilizing known structural characteristics of the surfactants and the counterions, along with the solution conditions, the free energy of mixed micellization is minimized to predict various optimal micelle characteristics, including the degree of counterion binding, the micelle composition, and the micelle shape and size. These predicted optimal micelle characteristics are then used to predict the critical micelle concentration (cmc) and the average micelle aggregation number. Our predictions of the degree of counterion binding, the cmc, and the average micelle aggregation number show good agreement with available experimental results from the literature for several binary surfactant mixtures. In addition, the theory is used to shed light on the relationship between the micelle composition, counterion binding and ion condensation, and the micelle shape transition.

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.

Aggregation of Perfluoropolyether Carboxylic Salts in Aqueous Solutions. Fluorescence Probe Study

Aqueous solutions of anionic surfactants Cl(C3F6O)nCF2COOX, consisting of n = 2 and 3 perfluoroisopropoxy units and the counterion X = Na+ or NH4+, were studied by the method of fluorescence quenching with the use of (1-pyrenylbutyl)trimethylammonium bromide as a luminophore, and 1,1'-dimethyl-4,4'bipyridinium dichloride (methyl viologen) as a quencher. From the kinetics of fluorescence decay (time-resolved experiments) micellar aggregation numbers, N, and rate constants of the intramicellar quenching were determined for a wide range of surfactant concentrations, on the basis of the model developed by Infelta and Tachiya. The results are discussed in terms of the shape of the aggregates and the degree of counterion binding. The most important conclusions include: (i) a significant increase of N with increasing surfactant concentration suggests that spherical micelles formed at critical micellar concentration (CMC) transform into ellipsoidal aggregates, (ii) the degree of counterion binding to micelles is higher for NH4+ than for Na+, leading to higher N values in the case of the ammonium salt (n = 2), and (iii) at concentrations close to CMC the longer chain surfactant (n = 3) forms loose aggregates suggesting significant permeation with water molecules. An additional finding of this study is that the micelle-bound luminophore and quencher can form a ground-state complex, and for this reason the N values cannot be evaluated properly from the steady-state fluorescence intensity data using the equation proposed by Turro and Yekta.

Cholesterol solubilization in aqueous micellar solutions of quillaja saponin, bile salts, or nonionic surfactants

Quillaja saponin in aqueous solution enhanced cholesterol solubility by as much as a factor of 10(3) at room temperature. Increased temperature and [NaCl] increased cholesterol solubility, whereas solubility was greatest at an aqueous pH of 4.6 at 298 K. Although various saponin sources were observed to differ in their abilities to solubilize cholesterol, trends in their solubilization properties with changing aqueous phase parameters were consistent. Surfactant molecules containing fused-ring structures as their hydrophobic portion, such as sodium cholate, sodium deoxycholate, and quillaja saponin, solubilized cholesterol significantly better than the linear hydrocarbon chain surfactants Tween 20 and Triton X-100. Mixtures of surfactants studied were found to exhibit synergistic effects: they formed micelles at lower concentrations than did those formed by the individual surfactants themselves, and they had a better ability to solubilize cholesterol. The knowledge obtained from these studies improves our understanding of cholesterol association with saponin and other types of surfactants and enhances the potential for using saponins for the solubilization and extraction of hydrophobic solutes in various pharmacological and industrial applications.