Molecular dynamics studies of cation aggregation in the room temperature ionic liquid [C10mim][Br] in aqueous solution

The structure of an aqueous 1-n-decyl-3-methylimidazolium bromide solution and its vapor-liquid interface has been studied using molecular dynamics (MD) simulations. Starting from an isotropic solution, spontaneous self-assembly of cations into small micellar aggregates has been observed. The decyl chains are buried inside the micelle to avoid unfavorable interactions with water, leaving the polar headgroups exposed to water. The cation aggregation numbers, ranging from 15 to 24 compare favorably with experimental estimates. Results are presented for the organization of solvent around the cations. The structure of the aggregates as determined from the present MD simulations does not support the staircase model proposed on the basis of nuclear magnetic resonance studies on similar aqueous ionic-liquid solutions. The distribution of ions in bulk solutions and at an air/water interface is also discussed.     

Effect of salts on local mobility and rheological properties of micelles of new long-chain surfactant

The effects of organic (sodium salicylate) and inorganic (KCl) salts on the rheological properties of micellar solutions and the local characteristics (local mobility and ordering) of micelle cores is studied for a cationic surfactant containing a long (C₁₈) unsaturated alkyl radical. The polar head of the surfactant contains two hydroxyl groups. The local characteristics are determined employing spin probe ESR spectroscopy. It is shown that the incorporation of a salt into a micellar solution reduces the local mobility of radicals of surfactant molecules in micelle cores and increases their local Lordering and the viscosity of the solution. Sodium salicylate has a stronger influence on the solution viscosity and the local characteristics of micelle cores than KCl does. Variations in the local characteristics of micelle cores under the action of the salts are in close correlation with variations in the rheological properties of the micellar solutions.     

Computer Simulations of the Formation of Bile Salt Micelles and Bile Salt/DPPC Mixed Micelles in Aqueous Solutions

Brownian dynamics simulations for a coarse-grained model have been performed to study the formation of micelles from bile salts and mixed micelles with dipalmitoyl-phosphatidylcholine (DPPC) in aqueous solutions. The particular association behavior of bile salts as facial surfactants was shown to be caused by their special molecular architecture with a hydrophilic and a hydrophobic side. The experimentally observed smooth transition into the micellar region with increasing concentration is reproduced. Micelle size distributions have been evaluated at different bile salt concentrations. Typical structures of pure bile salt micelles could be identified. The composition and the structure of mixed micelles have been studied in their dependence on the bile salt/lipid concentration ratio in the aqueous solution. We have found that the bile salt fraction in the mixed micelles increases considerably with increasing bile salt/lipid concentration ratio and decreasing micelle size. The structural and thermodynamic features of micelle formation in the aqueous bile salt solutions with DPPC, which we have studied with the coarse-grained model, are in good qualitative agreement with experimental findings.     

A diffusive anomaly of water in aqueous sodium chloride solutions at low temperatures

Molecular dynamics simulations are presented for the self-diffusion coefficient of water in aqueous sodium chloride solutions. At temperatures above the freezing point of pure water, the self-diffusion coefficient is a monotonically decreasing function of salt concentration. Below the freezing point of pure water, however, the self-diffusion coefficient is a non-monotonic function of salt concentration, showing a maximum at approximately one molal salt. This suggests that sodium chloride, which is considered a structure-making salt at room temperature, becomes a structure-breaking salt at low temperatures. A qualitative understanding of this effect can be obtained by considering the effect of ions on the residence time of water molecules near other water molecules. A consideration of the freezing point depression of aqueous sodium chloride solutions suggests that the self-diffusion coefficient of water in supercooled sodium chloride solutions is always higher than that in pure (supercooled) water at the same temperature.     

Effect of sodium salicylate and other electrolytes on the partition of 1-pentanol between sodium dodecylsulfate micelles and water: A gas-chromatographic study

The salt effect of sodium salicylate (NaS) on the micellization and micellar solubilization of sodium dodecylsulfate (NaDS) has been studied. The experimental and theoretical conditions for the determination of the thermodynamic partition coefficient P of 1-pentanol between the micellar pseudo-phase and water in presence of added salt is discussed in the case of a precise gas-chromatographic method. In Particular, it is shown that P decreases with solute concentration in aqueous NaDS and sodium perfluorooctanoate surfactant solutions in opposition to the classical behavior in water-organic immiscible phases. As a reference salt effect, it is shown that P is constant with added NaCl in a large salt concentration domain where NaDS micelles are known to undergo dranatic structural changes. In the case of added NaS, P decreases slightly at very high salt concentration. It is suggested that this behavior might be the consequence of partial mixed micelle formation between the salicylate ion and NaDS micelles.     

Binary interactions and salt-induced coalescence of spherical micelles of cationic surfactants from molecular dynamics simulations

A direct estimation of salt-mediated potential of mean force (PMF) between spherical micelles of cationic surfactants is obtained for the first time using molecular dynamics (MD) simulations. Coarse-grained (CG) potentials benchmarked in an earlier study [Langmuir, 2011, 27(11), 6628-6638] are used to model a binary system of cetyltrimethylammonium chloride (CTAC) surfactant micelles at varying concentrations of sodium chloride (NaCl) or sodium salicylate (NaSal). The shape and structure of micelles are not subject to external constraints. NaSal is significantly more efficient in screening the intermicelle repulsive interactions shown by the PMF compared to NaCl due to a stronger binding of salicylate counterions to the micelle corona. Upon contact with each other, the micelles coalesce in the presence of NaSal to form a cylindrical structure which is stabilized by the adsorbed salicylate anions. Comparison of the PMF with Derjaguin-Landau-Verwey-Overbeek (DLVO) potentials shows qualitative agreement, while the magnitude of PMF is significantly greater than that of the DLVO potentials. To understand this discrepancy, PMF is evaluated by turning off (a) long-ranged electrostatic interactions and (b) solvent polarizability. The above effects are shown to play an important role in determining the solvent-mediated and ion-correlated interactions between the two micelles, which are not explicitly captured by mean-field double layer theories such as DLVO.     

Solvation and aggregation of n,n'-dialkylimidazolium ionic liquids: a multinuclear NMR spectroscopy and molecular dynamics simulation study

The solvation and aggregation of the ionic liquid (IL) 1-n-butyl-3-methylimidazolium chloride ([C4mim]Cl) in water and dimethylsulfoxide (DMSO) were examined by analysis of (1)H and (35/37)Cl chemical shift perturbations and molecular dynamics (MD) simulations. Evidence of aggregation of the IL n-butyl chains in aqueous environments at IL concentrations of 75-80 wt% was observed both in the NMR experiments and in the MD simulations. The studies also show that [C4mim]Cl behaves as a typical electrolyte in water, with both ions completely solvated at low concentrations. On the other hand, the data reveal that the interactions between the [C4mim](+) and Cl(-) ions strengthen as the DMSO content of the solutions increases, and IL-rich clusters persist in this solvent even at concentrations below 10 wt%. These results provide an experimentally supported atomistic explanation of the effects that these two solvents have on some of the macroscopic properties of [C4mim]Cl. The implications that these findings could have on the design of IL-based solvent systems are briefly discussed.     

Effect of sodium salicylate, sodium oxalate, and sodium chloride on the micellization and adsorption of sodium deoxycholate in aqueous solutions

The salicylate ion increases the rate of bile flow (choleretic effect) and bile salts are known to affect the colonic absorption of oxalate. Owing to this physiological relevance of salicylate and oxalate ions, critical micelle concentration (cmc) values of sodium deoxycholate (NaDC) were determined in aqueous sodium oxalate, sodium salicylate, and sodium chloride solutions by using surface tension, fluorescence, and EMF methods. The results indicate, besides a counterion effect, the influence of coanions on the cmc. In the range from 25 to 40 °C, cmc increases almost linearly with temperature. In the temperature range from 30 to 40 °C, the counterion binding constant β of NaDC micelles has the same value (0.17±0.01) in the presence of sodium chloride and sodium salicylate. On the other hand, in sodium oxalate solution β=0.05±0.02 when oxalate concentration is less than or equal to c* and β=0.48±0.04 above c*, where c*≈0.038 mol kg(-1). EMF measurements also supported this type of counterion binding to NaDC micelles in sodium oxalate solutions. In sodium oxalate solution, at c* a change in the shape of deoxycholate micelles is expected to take place. Salicylate, oxalate, and chloride coanions have a similar effect on the adsorption of NaDC. This study reveals that the choleretic effect of salicylate is not due to the influence of salicylate ions on the micellization of NaDC.     

Thermodynamic studies on the interaction of some ureas and salts with micellar triton-X-100 in aqueous solutions

We have determined the enthalpies of solution in the micellar state (ΔH_s) for Triton-X-100 in 1 m aqueous solutions of urea, 1,3-dimethyl urea, tetramethyl urea, sodium chloride and calcium chloride at 298.15 K and 308.15 K. These results were used to evaluate the heat capacities of solution (ΔC_p,s) for Triton-X-100 micelles in these solvent systems. It has been observed that ΔC_p,s values of micellar Triton-X-100 decreases drastically upon transfer from water to these solutions but is positive in all cases. Thus, the heat capacities of transfer of Triton-X-100 micelles (ΔC_p,tr) are negative in all the systems studied. A comparison of the effect of non-electrolytes (ureas) and electrolytes (salts) on the micelle has been presented. The results have been discussed in terms of the relative water-structure-disrupting tendencies of the ureas and the salts.     

Effect of concentration on surfactant micelle shapes——A molecular dynamics study

Above the critical micellar concentration (cmc), micelles with a wide variety of structures and shapes are formed with the increase of surfactant concentra-tion in surfactant-water or surfactant-water-oil systems, such as spherical micelles, rodlike micelles, and bilay-ers. The viewpoint that micelle should be in spheres of constant size was first proposed by Hartley[1]. Later experiments by light scatter indicated that most mi-celles were indeed spherical, and their aggregation numbers were c…     

Revised implicit solvent model for the simulation of surfactants in aqueous solutions

The implicit solvent model (ISM) proposed previously for the simulation of surfactant aqueous solutions, in which no water molecules of the solvent are treated explicitly, but the effects are incorporated using the solvent-averaged interactions between the surfactant segments in water at infinite dilution, has been revised to represent the surfactant aggregates more appropriately. In the revised model (ISM-2), the interactions between the hydrophobic sites of the surfactants are varied depending on their surroundings, namely, the local hydrocarbon density. The ISM-2 has been applied to the molecular dynamics simulations of (i) the single n-hexane droplets of different sizes in water and (ii) the single micelle composed of 30 n-decyltrimethylammonium chloride (C10TAC) cationic surfactants. As a result, it was found that the ISM-2 can mimic the n-hexane/water interface and represent the fluidity of the hydrocarbon interior of the surfactant micelle that the original ISM fails to do. The results will be compared to those from experiments and atomistic model simulations.     

Effect of concentration on surfactant micelle shapes —A molecular dynamics study

Many aspects of the behavior of surfactants have not been well understood due to the coupling of many different mechanisms. Computer simulation is, therefore, attractive in the sense that it can explore the effect of different mechanisms separately. In this paper, the shapes, structures and sizes of sodium dodecylbenzenesulfonate (SDBS) micelles under different concentrations in an oil/water mixture were studied via molecular dynamics (MD) simulations using a simplified atomistic model which basically maintains the hydrophile and lipophile properties of the surfactant molecules. Above the critical micellar concentration (cmc), surfactant molecules aggregate spontaneously to form a wide variety of assemblies, from spherical to rodlike, wormlike and bilayer micelles. Changes in their ratios of the principle moments of inertia (g₁/g₃, g₂/g₃) indicated the transition of micelle shapes at different concentrations. The aggregation number of micelle is found to have a power-law dependence on surfactant concentration.     

Multifold curdlan gel formation by dialysis into aqueous solutions of metal salts

We have found that the dialysis of curdlan dissolved in alkaline solution into aqueous solutions of metal salts yielded multifold gel structures. Aqueous sodium chloride and potassium chloride as well as pure water induced isotropic gels. Aqueous calcium salts induced liquid crystalline gel with refractive index gradient/amorphous gel alternative structure. Aqueous salts of trivalent aluminum and ferric cations induced a rigid liquid crystalline gel, which shrank above a threshold concentration of each salt. On the other hand, Liesegang ring-like pattern was observed with aqueous solutions of mixed salts of calcium chloride and magnesium chloride. The patterns have been classified to discuss the mechanism of forming the variety of structures.     

From crystal to micelle: A new approach to the micellar structure

The structure of sodium and rubidium deoxycholate micellar aggregates in aqueous solutions was found to be helical and to be stabilized mainly by polar interactions. Astonishingly, the lateral surface of the helix is covered by nonpolar groups and the interior part is filled with cations surrounded by water molecules, as in the case of an inverted micelle. This helical model was inferred from the crystal structures of sodium and rubidium deoxycholates and proved by spectroscopic and diffractometric experimental data. The strategy of the approach to the determination of the micellar structure and the comparison with another model, previously proposed for the bile salt micelles, are reported. On the basis of some results obtained for sodium tauro- and glyco-deoxycholates, micellar models are suggested which could account for the biological function of these important conjugated bile salts.     

Molecular dynamics simulation of solution structure and dynamics of aqueous sodium chloride solutions from dilute to supersaturated concentration

Constant-temperature and constant-pressure (NpT) molecular dynamics simulations were performed to study the effects of salt concentration ranging from dilute to supersaturated concentrations on solution structure and dynamical properties of aqueous sodium chloride solutions at 298 K. The rigid SPC/E model was used for water molecules, and sodium and chloride ions were modeled as charged Lennard–Jones particles. Na⁺–Cl⁻ radial distribution functions showed the presence of contact ion pairs and solvent separated ion pairs. The coordination numbers of Na⁺–Cl⁻ ion pairs increased with salt concentration up to saturated concentration, although the number of contact ion pairs was almost constant in supersaturated regions. The tracer diffusion coefficients of both ions decreased with salt concentration up to saturated concentration, while that of sodium ion was almost constant in supersaturated regions. The tracer diffusion coefficients of both ions were therefore quite close to each other. The constant number of the contact ion pairs and the almost equality of the tracer diffusion coefficients of both ions would lead to the formation of clusters in supersaturated solutions.     

Asymmetry of lipid bilayers induced by monovalent salt: atomistic molecular-dynamics study

Interactions between salt ions and lipid components of biological membranes are essential for the structure, stability, and functions of the membranes. The specific ionic composition of aqueous buffers inside and outside of the cell is known to differ considerably. To model such a situation we perform atomistic molecular-dynamics (MD) simulations of a single-component phosphatidylcholine lipid bilayer which separates two aqueous reservoirs with and without NaCl salt. To implement the difference in electrolyte composition near two membrane sides, a double bilayer setup (i.e., two bilayers in a simulation box) is employed. It turns out that monovalent salt, being in contact with one leaflet only, induces a pronounced asymmetry in the structural, electrostatic, and dynamical properties of bilayer leaflets after 50 ns of MD simulations. Binding of sodium ions to the carbonyl region of the leaflet which is in contact with salt results in the formation of "Na-lipids" complexes and, correspondingly, reduces mobility of lipids of this leaflet. In turn, attractive interactions of chloride ions (mainly located in the aqueous phase close to the water-lipid interface) with choline lipid groups lead to a substantial (more vertical) reorientation of postphatidylcholine headgroups of the leaflet adjoined to salt. The difference in headgroup orientation on two sides of a bilayer, being coupled with salt-induced reorientation of water dipoles, leads to a notable asymmetry in the charge-density profiles and electrostatic potentials of bilayer constitutes of the two leaflets. Although the overall charge density of the bilayer is found to be almost insensitive to the presence of salt, a slight asymmetry in the charge distribution between the two bilayer leaflets results in a nonzero potential difference of about 85 mV between the two water phases. Thus, a transmembrane potential of the order of the membrane potential in a cell can arise without ionic charge imbalance between two aqueous compartments.     

Local order in aqueous NaCl solutions and pure water: X-ray scattering and molecular dynamics simulations study

The microstructures of pure water and aqueous NaCl solutions over a wide range of salt concentrations (0-4 m) under ambient conditions are characterized by X-ray scattering and molecular dynamics (MD) simulations. MD simulations are performed with the rigid SPC water model as a solvent, while the ions are treated as charged Lennard-Jones particles. Simulated data show that the first peaks in the O...O and O...H pair correlation functions clearly decrease in height with increasing salt concentration. Simultaneously, the location of the second O...O peak, the signature of the so-called tetrahedral structure of water, gradually disappears. Consequently, the degree of hydrogen bonding in liquid water decreases when compared to pure fluid. MD results also show that the hydration number around the cation decreases as the salt concentration increases, which is most likely because some water molecules in the first hydration shell are occasionally substituted by chlorine. In addition, the fraction of contact ion pairs increases and that of solvent-separated ion pairs decreases. Experimental data are analyzed to deduce the structure factors and the pair correlation functions of each system. X-ray results clearly show a perturbation of the association structure of the solvent and highlight the appearance of new interactions between ions and water. A model of intermolecular arrangement via MD results is then proposed to describe the local order in each system, as deduced from X-ray scattering data.     

The phenomenon of water penetration into sodium octanoate micelles studied by molecular dynamics computer simulation

 In this publication we have studied the penetration process of water molecules into the hydrophobic core of a sodium octanoate micelle. The analysis of this phenomenon was based on a molecular dynamics computer simulation. We calculated the probability to find water molecules within a specific sphere which was adjusted in the center of the micelle. It turned out that the position of the micellar mass and geometry center was not too different, so that this reference point was well characterized. Water penetration was observed within the whole aggregate but if the radius is smaller than 300 pm, polar solvent molecules are very rarely observed. The results of our computer simulations suggest that significant water diffusion into the micelle occurs at larger distances from the micellar center with a lower threshold value of about 400 pm. In addition to these calculations, we used the Connolly algorithm in order to determine the solvent accessible surfaces of different micellar equilibrium structures. We observed large dynamical fluctuations with the formation of pores and channels. These structures are occasionally filled with water molecules. Received: 29 April 1998 Accepted: 27 May 1998     

Solvent effects and hydration of a tripeptide in sodium halide aqueous solutions: an in silico study

In this work we are trying to gain insight into the mechanisms of ion-protein interactions in aqueous media at the molecular scale through fully atomistic molecular dynamics simulations. We present a systematic molecular simulation study of interactions of sodium and halide ions with a trialanine peptide in aqueous sodium halide solutions with different salts concentrations (0.20, 0.50, 1.0 and 2.0 M). Each simulation covers more than fifty nanoseconds to ensure the convergence of the results and to enable a proper determination of the tripeptide-ion interactions through the potentials of mean force. Changes in ion densities in the vicinity of different peptide groups are analysed and implications for the tripeptide conformations are discussed.     

Sphere-to-rod transitions of nonionic surfactant micelles in aqueous solution modeled by molecular dynamics simulations

Control of the size and agglomeration of micellar systems is important for pharmaceutical applications such as drug delivery. Although shape-related transitions in surfactant solutions are studied experimentally, their molecular mechanisms are still not well understood. In this study, we use coarse-grained molecular dynamics simulations to describe micellar assemblies of pentaethylene glycol monododecyl ether (C(12)E(5)) in aqueous solution at different concentrations. The obtained size and aggregation numbers of the aggregates formed are in very good agreement with the available experimental data. Importantly, increase of the concentration leads to a second critical micelle concentration where a transition to rod-like aggregates is observed. This transition is quantified in terms of shape anisotropy, together with a detailed structural analysis of the micelles as a function of aggregation number.