Parametric studies of interaction strengths in polymer/CO2 systems: discontinuous molecular dynamics simulations

Discontinuous molecular dynamics simulations are performed on homopolymer/solvent and surfactant/solvent systems. The homopolymer and surfactant molecules are modeled as freely jointed square-well chains. Solvent molecules are modeled as both hard spheres and square-well spheres. We explore how the various interaction parameters affect the types of phase behavior and micellization observed in the homopolymer/solvent and surfactant/solvent systems. Increasing the packing fraction of homopolymers in both hard-sphere solvents and square-well solvents increases the solvent's ability to dissolve homopolymers only when the segment-solvent interaction strength exceeds a critical value. Although only upper critical solution temperature (UCST) behavior is observed for homopolymers in hard-sphere solvents, both UCST and lower critical solution temperature (LCST) behavior are observed for homopolymers in square-well solvents, depending upon the interaction strengths and chain length. This indicates that it is necessary to account for the solvent-solvent attraction to model LCST behavior in supercritical CO2. Our simulation results on surfactants in hard-sphere solvents show that it is necessary to account for the interactions experienced by both the head and tail blocks in order to capture the essential features of surfactant/supercritical CO2 systems.     

Atomistic structure of a micelle in solution determined by wide Q-range neutron diffraction

The accepted picture of the structure of a micelle in solution arises from the idea that the surfactant molecules self-assemble into a spherical aggregate, driven by the conflicting affinity of their head and tail groups with the solvent. It is also assumed that the micelle's size and shape can be explained by simple arguments involving volumetric packing parameters and electrostatic interactions. By using wide Q-range neutron diffraction measurements of H/D isotopically substituted solutions of decyltrimethylammonimum bromide (C(10)TAB) surfactants, we are able to determine the complete, atomistic structure of a micelle and its surroundings in solution. The properties of the micelle we extract are in agreement with previous experimental studies. We find that ~45 surfactant molecules aggregate to form a spherical micelle with a radius of gyration of 14.2 ? and that the larger micelles are more ellipsoidal. The surfactant tail groups are hidden away from the solvent to form a central dry hydrophobic core. This is surrounded by a disordered corona containing the surfactant headgroups, counterions, water, and some alkyl groups from the hydrophobic tails. We find a Stern layer of 0.7 bromide counterion per surfactant molecule, in which the bromide counterions maintain their hydration shells. The atomistic resolution of this technique provides us with unprecedented detail of the physicochemical properties of the micelle in its solvent.     

Monte Carlo simulations of Lennard-Jones nonionic surfactant adsorption at the liquid/vapor interface

We present Monte Carlo simulations of nonionic surfactant adsorption at the liquid/vapor interface of a monatomic solvent. All molecules in the system, solvent and surfactant, are characterized by the Lennard-Jones (LJ) potential using differing interaction parameters. Surfactant molecules consist of an amphiphilic chain with a solvophilic head and a solvophobic tail. 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 quantitative thermodynamic adsorption and surface tension isotherms in addition to surfactant radius of gyration, tilt angles, and potentials of mean force. Surface tension simulations compared to those calculated from the simulated adsorbed amounts and the Gibbs adsorption isotherm agree confirming equilibrium in our simulations. We find that the classical Langmuir isotherm is obeyed for our LJ surfactants over the range of head and tail lengths studied. Although simulated surfactant chains in the bulk solution exhibit random orientations, surfactant chains at the interface orient roughly perpendicular and the tails elongate compared to bulk chains even in the submonolayer adsorption regime. At a critical surfactant concentration, designated as the critical aggregation concentration (CAC), we find aggregates in the solution away from the interface. At higher concentrations, simulated surface tensions remain practically constant. Using the simulated potential of mean force in the submonolayer regime and an estimate of the surfactant footprint at the CAC, we predict a priori the Langmuir adsorption constant, KL, and the maximum monolayer adsorption, Gammam. Adsorption is driven not by proclivity of the surfactant for the interface, but by the dislike of the surfactant tails for the solvent, that is by a "solvophobic" effect. Accordingly, we establish that a coarse-grained LJ surfactant system mimics well the expected equilibrium behavior of aqueous nonionic surfactants adsorbing at the air/water interface.     

Molecular dynamics simulation of amphiphilic dimers at a liquid-vapor interface

Molecular dynamics simulations are utilized to simulate a model liquid-vapor-amphiphile system. Amphiphilic surfactant molecules are modeled as dimers composed of a hydrophilic head and a hydrophobic tail. Three dimer models with three different head sizes and two different head-to-tail size ratios are studied. The surfactant molecules distribute preferentially at the interfaces at low concentrations and form micelles in the bulk liquid phase as the concentration increases. We find that the surface tension decreases as molecular concentration increases, with a reduction in the rate of decrease after micellization occurs. The extent to which a surfactant can reduce the surface tension at a given concentration is found to depend on the head size. Furthermore, the head size and concentration dependence of the surfactant tilt-angle distribution is studied and compared to experimental data.     

Phase behavior in model homopolymer/CO2 and surfactant/CO2 systems: discontinuous molecular dynamics simulations

Discontinuous molecular dynamics simulations are performed on surfactant (HmTn)/solvent systems modeled as a mixture of single-sphere solvent molecules and freely jointed surfactant chains composed of m slightly solvent-philic head spheres (H) and n solvent-philic tail spheres (T), all of the same size. We use a square-well potential to account for the head-head, head-solvent, tail-tail, and tail-solvent interactions and a hard-sphere potential for the head-tail and solvent-solvent interactions. We first simulate homopolymer/supercritical CO2 (scCO2) systems to establish the appropriate interaction parameters for a surfactant/scCO2 system. Next, we simulate surfactant/scCO2 systems and explore the effect of the surfactant volume fraction, packing fraction, and temperature on the phase behavior. The transition from the two-phase region to the one-phase region is located by monitoring the contrast structure factor of the equilibrated surfactant/scCO2 system, and the micelle to unimer transition is located by monitoring the aggregate size distribution of the equilibrated surfactant/scCO2 system. We find a two-phase region, a micelle phase, and a unimer phase with increasing packing fraction at fixed temperature or with increasing temperature at fixed packing fraction. The phase diagram for the surfactant/scCO2 system in the surfactant volume fraction-packing fraction plane and the density dependence of the critical micelle concentration are in qualitative agreement with experimental observations. The phase behavior of a surfactant/scCO2 system can be directly related to the solubilities of the corresponding homopolymers that serve as the head and tail blocks for the surfactant. The influence of surfactant structure (head and tail lengths) on the phase transitions is explored.     

Structures of micelles formed by synthetic alkyl glycosides with unsaturated alkyl chains

Three new alkyl glycosides with similar molecular structures (oleyl and oleoyl alkyl chains and various head groups: disaccharide, trisaccharide and disaccharide with an additional amidoethoxy spacer) were synthesized and their supramolecular structure in aqueous solution was investigated. Small angle neutron scattering, surface tension measurement and the contact preparation method were applied to get molecular structure-property relationships. Although the chemical structures differ only in small details, their CMC values, lyotropic phase behaviour, surface area per surfactant molecule in the micelle and at the liquid-air interface, and the size and shape of the micelles are very different. We have found three different types of aggregates: spherical, cylindrical and polymer-like micelles in dilute solutions.     

Effect of short-chain alcohols on surfactant-mediated reversed-phase liquid chromatographic systems

The behaviour of β-blockers in a reversed-phase liquid chromatographic (RPLC) column with mobile phases containing a short-chain alcohol (methanol, ethanol or 1-propanol), with and without the surfactant sodium dodecyl sulphate (SDS), was explored. Two surfactant-mediated RPLC modes were studied, where the mobile phases contained either micelles or only surfactant monomers at high concentration. Acetonitrile was also considered for comparison purposes. A correlation was found between the effects of the organic solvent on micelle formation (monitored by the drop weight procedure) and on the nature of the chromatographic system (as revealed by the retention, elution strength and peak shape of β-blockers). When SDS is added to the mobile phase, the free surfactant monomers bind the C18 bonded chains on the stationary phase, forming an anionic layer, which attracts strongly the cationic β-blockers. The retention is modified as a consequence of the solving power of the organic solvent, micelles and surfactant monomers. The molecules of organic solvent bind the micelles, modify their shape, and may avoid their formation. They also bind the monomers of surfactant, desorbing them from the stationary phase, which affects the retention. The remaining surfactant covers the free silanols on the siliceous support, avoiding the interaction with the cationic solutes. The retention of β-blockers results from a combination of electrostatic and hydrophobic interactions, the latter being weaker compared to the hydro-organic system. The peak efficiencies and asymmetries are excellent tools to probe the surfactant layer on the stationary phase in an SDS/organic solvent system. The peaks will be nearly symmetrical wherever enough surfactant coats the stationary phase (up to 60% methanol, 40% ethanol, 35% 1-propanol, and 50% acetonitrile).     

三种不同分子结构阴离子表面活性剂胶束微结构的NMR研究

用核磁共振测定自旋-晶格、自旋-自旋弛豫时间(t₁,t₂)、自扩散系数(D),用2DNOESY技术对正十四烷基硫酸钠、β-戊基壬烷基硫酸钠和β-戊基壬烷基聚氧乙烯醚(4)硫酸钠三类阴离子表面活性剂水溶液进行了观测,烷烃链各基团的t₂/t₁值给出了这三类分子形成各自胶束的水合层位点信息以及烷烃链在胶束内核中堆积程度的比较,自扩散系数结果表明,β-戊基壬烷基硫酸钠比正十四烷基硫酸钠形成的胶束的水合动力学半径小,但β-戊基壬烷基聚氧乙烯醚(4)硫酸钠形成的胶束水合动力学半径明显大于其它两类表面活性剂胶束,2DNOESY谱图提供了β-戊基壬烷基聚氧乙烯醚(4)硫酸钠分子中聚氧乙烯基键在胶束外层卷曲排列的信息.     

Removal of organic dyes from water by liquid-liquid extraction using reverse micelles

In the present work, solvent extraction using reverse micelles is proposed for the removal of organic dyes from water. In this approach, the dye is solubilized in the aqueous core of the reverse micelles, which are present in the organic phase. The organic phase is subsequently separated from the aqueous phase leading to signifi-cant removal of dye. Experimental results reveal that the electrostatic interaction between the oppositely charged surfactant head group present in the reverse micelles and the dye molecule plays a key role in the separation. The removal of the anionic methyl orange dye from water is carried out in the presence of cationic hexadecyltrimethyl ammonium bromide surfactant, whereas the removal of the cationic methylene blue dye is carried out in the presence of anionic sodium dodecylbenzene sulfonate surfactant. Amyl alcohol is used as the solvent. The influence of parameters such as dye concentrations, surfactant concentrations, pH, and KCl and NaBr concentrations on the percentage removal of dye was studied. The percentage removal of dye is decreased with the increase in dye concentration in the feed. The increase in surfactant concentration resulted in higher dye removal, because more reverse micelles could be hosted in the organic phase. The increase in aqueous phase pH resulted in enhanced removal of methyl orange from water, while in the case of methylene blue the percentage removal decreased. The increase in KCl and NaBr concentrations resulted in decreased percentage removal of methylene blue, whereas the percentage removal of methyl orange was increased. The effect of pH and salt concentration is explained based on charge transfer mechanism and electrostatic interactions and dye-surfactant complex formation.     

The thermodynamics of micelles in surfactant solution: Lattice mean field theory

The lattice self-consistent mean field theory is applied to study the thermodynamics of micelles in surfactant solution. The model surfactants used are H ₄ T ₄ and H ₂ T ₄. The formation of spherical micelles is considered. The effect of the head length on the thermodynamic stability of the micellar solution is examined. The critical micelle concentration is studied at different lengths of head, fractional charge, solvent quality, and univalent salt concentration. A lower critical micelle concentration is associated with a larger aggregation number, while the smallest micelles are found at the lowest univalent salt concentration. The text was submitted by the authors in English.     

Soluble complexes of sodium poly(isoprene-b-methacrylate) micelles with cationic surfactants in aqueous media

Water-soluble complexes between sodium poly(isoprene-b-methacrylate) (NaIMA) amphiphilic block copolymer micelles and two cationic surfactants with different hydrophobic tail lengths, namely, dodecyltrimethylammonium bromide (DTMAB) and octyltrimethylammonium bromide (OTMAB), were prepared by mixing individual aqueous solutions of block copolymers and surfactants. The complexes were characterized in terms of size, overall charge, and micropolarity by dynamic light scattering, zeta-potential measurements, and fluorescence spectroscopy. Properties of the systems were investigated as a function of surfactant concentration and surfactant type and state in the initial solutions, as well as temperature. Experiments reveal surfactant complexation at the coronal sodium poly(methacrylate) (NaMA) chains, followed by an increase in mass and a decrease in size of the micelles. Complexation of individual surfactant micelles was observed when the DTMAB concentration in the starting solutions was higher than the surfactant cmc. The complexes show a temperature dependence of their dimension due to the hydrophobic effect.     

Simulation and model development for the equation of state of self assembling non-additive hard chain–hard monomer mixture

Molecular dynamics simulations for a short hard chain composed of a head and three tail groups interacting with non-additive size interactions with a hard sphere solvent were performed. Different densities and non-additivities were used. The equation of state for this mixture was investigated and models based on the first-order thermodynamic perturbation theory (TPT1) and the polymeric analog of the Percus–Yevick approximation (PPY) were developed to predict the compressibility factor of the mixture. The models predicted the compressibilities of the mixtures accurately at zero and negative non-additivities. However, at positive non-additivities, the models overpredicted the compressibilities especially at high densities. The TPT1 model was generally more accurate in predicting the compressibility factor than the PPY model. Microphase separation was observed at high densities and positive non-additivities.     

Aggregate morphologies of amphiphilic graft copolymers in dilute solution studied by self-consistent field theory

Microstructures assembled by amphiphilic graft copolymers in a selective solvent (poor for the backbone chain and good for graft chains or poor for graft chains and good for the backbone chain) were investigated on the basis of a real-space algorithm of self-consistent field theory in two-dimensions. Circle-like micelles, line-like micelles, large compound micelles, and vesicles are obtained by tailoring the architectural parameters and interaction parameter between the graft blocks and solvents. The aggregate morphology stability regions of graft copolymers as functions of the position of first graft point and the number of branches are constructed. It is found that the architectural parameters play a remarkable role in the complex microstructure formation. The interaction between the graft blocks and solvents is also shown to exert an effect on the morphology stability regions. The distributions of the free end and inner blocks of the backbone are found to be different in various aggregate structures. For the circle-like micelles assembled by graft copolymers with a hydrophobic backbone and vesicles assembled by graft copolymers with a hydrophilic backbone, the free end and inner blocks segregate and localize in different parts of the aggregates depending on their length. However, with respect to the large compound micelles and vesicles assembled by graft copolymers with a hydrophobic backbone, the free end and inner blocks uniformly mix in the clusters.     

Model calculations of the spontaneous curvature, mean and Gaussian bending constants for a thermodynamically open surfactant film

The spontaneous curvature (H(0)), mean and Gaussian bending constants (k(c) and k (c)), as defined in the well-known Helfrich expression, have been calculated from a detailed model for a thermodynamically open surfactant layer. The effect of head group cross-section area, surfactant tail length and electrolyte concentration for monovalent ionic surfactants have been investigated. Geometrical packing constraints subjected to the aggregated hydrocarbon tails and electrostatics are found to be the dominant contributions to H(0), k(c) and k (c). In addition, the transition from spherocylindrical micelles to vesicles were investigated in terms of the three parameters and the following simple expressions were derived as criteria for coexistence between micelles and vesicles H(0)=1/4 xi and N(ves)/N(mic)=exp[4 pi(k(c)+k (c))/kT], where xi is the thickness of the hydrocarbon part of the film and N(mic) and N(ves) the average aggregation numbers of micelles and vesicles, respectively. However, it is found that the ratio N(ves)/N(mic) is order of magnitudes too large for vesicles to form at all in charged single-surfactant systems where the surfactant head is of moderate size.     

表面活性剂对污泥脱水性能影响的机理研究:Gemini表面活性剂与聚电解质相互作用的分子动力学模拟

表面活性剂可以与污泥表面的胞外聚合物(EPS)吸附形成胶束,释放出自由水和结合水,从而达到改善污泥脱水性能的目的.本文采用粗粒化的分子动力学模拟方法,研究了Gemini表面活性剂与EPS形成复合物的过程和结构.聚电解质链的亲疏水性对吸附过程有显著影响,亲水聚电解质链与Gemini表面活性剂吸附的主要驱动力为静电吸引,Gemini表面活性剂头基吸附在链上,尾链朝向溶剂;疏水聚电解质链与Gemini表面活性剂吸附过程由静电作用与疏水作用共同促进,Gemini表面活性剂以平行于聚电解质链的构型存在.Gemini表面活性剂联结基团长度对吸附过程的影响甚微;聚电解质链的电荷密度对亲水聚电解质链的吸附产生协同作用,对疏水聚电解质链的吸附不产生作用.     

Structure and dynamics of surfactin studied by NMR in micellar media

The NMR structure of the cyclic lipopeptide surfactin from Bacillus subtilis was determined in sodium dodecyl sulfate (SDS) micellar solution. The two negatively charged side chains of surfactin form a polar head opposite to most hydrophobic side chains, accounting for its amphiphilic nature and its strong surfactant properties. Disorder was observed around the fatty acid chain, and 15N relaxation studies were performed to investigate whether it originates from a dynamic phenomenon. A very large exchange contribution to transverse relaxation rate R(2) was effectively observed in this region, indicating slow conformational exchange. Temperature variation and Carr-Purcell-Meiboom-Gill (CPMG) delay variation relaxation studies provided an estimation of the apparent activation energy around 35-43 kJ x mol(-1) and an exchange rate of about 200 ms(-1) for this conformational exchange. 15N relaxation parameters were also recorded in dodecylphosphocholine (DPC) micelles and DMSO. Similar chemical exchange around the fatty acid was found in DPC but not in DMSO, which demonstrates that this phenomenon only occurs in micellar media. Consequently, it may either reflect the disorder observed in our structures determined in SDS or originate from an interaction of the lipopeptide with the detergent, which would be qualitatively similar with an anionic (SDS) or a zwitterionic (DPC) detergent. These structural and dynamics results on surfactin are the first NMR characterization of a lipopeptide incorporated in micelles. Moreover, they provide a model of surfactin determined in a more biomimetic environment than an organic solvent, which could be useful for understanding the molecular mechanism of its biological activity.     

Modeling and simulation of equation of state for nonadditive hard disk mixtures

We present exact results for mixtures of nonadditive hard disks and use some of them to derive a consistent model for the equation of state. We also performed molecular dynamics simulation for hard disks over a wide range of size ratios. Comparison of the model to the data shows that the model is accurate for all densities in the case of additive and slightly nonadditive (nonadditivity parameter within ±0.1) mixtures. For large nonadditivity, the model is accurate for low to moderate densities only, and starts to deteriorate at high densities.     

Fluorescence spectroscopy of fulvic acids’ interaction with surfactants

Fluorescence spectra of two fulvic acid (FA) samples, FA0 from underground water and FA1 from forest soil, were recorded in various surfactant solutions. Alkyltrimethylammonium ions with different alkyl chain lengths induced a decrease in the fluorescence intensity for both FAs at concentrations below the critical micelle concentration (cmc) and an enhancement above the cmc. The intensity minimum thus obtained at the cmc was deeper for surfactants with longer alkyl chains. This effect was attributable to the formation of insoluble FA–surfactant complexes below the cmc and to the solubilization of the complex into micelles above the cmc. Dodecylpyridinium chloride caused a monotonic decrease in the FA fluorescence even far above the cmc. This was attributable to the quenching of FA fluorescence by the positioning of the pyridinium head group near the FA fluorophore. Anionic and nonionic surfactants showed little to no effect on the FA fluorescence.     

Enantioseparations by capillary electrophoresis using chiral glycosidic surfactants. I. Evaluation of cyclohexyl-pentyl-beta-D-maltoside surfactant.

Chiral cyclohexyl-pentyl-beta-D-maltoside (CYMAL-5) surfactant was evaluated in the enantioseparation of charged racemic species by capillary electrophoresis. CYMAL-5 is a glycosidic surfactant (GS) with a chiral maltose polar head group and a cyclohexyl-pentyl hydrophobic tail. At concentrations above its critical micellar concentration (CMC), CYMAL-5 produces neutral micelles in aqueous media. The neutral micelles migrate at the velocity of the electroosmotic flow (EOF). As expected, the CYMAL-5 system was only useful for the enantioseparation of charged chiral solutes. The enantioresolution of the CYMAL-5 can be manipulated over a wide range of electrolyte composition, e.g., pH, ionic strength and surfactant concentration. In the presence of EOF, and in all cases, there is an optimum surfactant concentration for maximum enantioresolution, which is located at low surfactant concentration for strongly hydrophobic solutes and at high surfactant concentration for relatively hydrophilic solutes. The presence of an optimum surfactant concentration for maximum enantioresolution is attributed to the EOF. At low pH values where the EOF is negligible, enantioresolution increased with increasing surfactant concentration in the useful concentration range in a way similar to chromatography.     

Salt effect on the complex formation between polyelectrolyte and oppositely charged surfactant in aqueous solution

The complex formation between sodium carboxymethylcellulose (NaCMC) and dodecyltrimethylammonium bromide (DTAB) at various sodium bromide concentrations (C(NaBr)) has been studied by microcalorimetry, turbidimetric titration, steady-state fluorescence measurements, and the fluorescence polarization technique. The addition of salt is found to influence the formation of NaCMC/DTAB complexes markedly. At C(NaBr) = 0.00, 0.01, 0.02, 0.10, and 0.20 M, DTAB monomers form micelle-like aggregates on NaCMC chains to form NaCMC/DTAB complexes above the critical surfactant concentration (C1). At C(NaBr) = 0.23 M, DTAB molecules first form micelles above a 2.46 mM DTAB concentration prompted by the added salt, and then, above C1 = 4.40 mM, these micelles can aggregate with NaCMC chains to form NaCMC/DTAB complexes. However, at C(NaBr) = 0.25 M, there is no NaCMC/DTAB complex formation because of the complete salt screening of the electrostatic attraction between DTAB micelles and NaCMC chains. It is also surprisingly found that the addition of NaBr can bring out a decrease in C1 at C(NaBr) < 0.20 M. Moreover, the addition of NaBr to a mixture of 0.01 g/L NaCMC and 3.6 mM DTAB can directly induce the formation of NaCMC/DTAB complexes. This salt-enhancing effect on the complex formation is explained as the result of competition between the screening of interaction of polyelectrolyte with surfactant and the increasing of polyelectrolyte/surfactant interaction owing to the growth of micelles by added salt. When the increasing of polyelectrolyte/surfactant interaction exceeds the screening of interaction, the complex formation can be enhanced.