Using Microemulsions: Formulation Based on Knowledge of Their Mesostructure.

Microemulsions, as thermodynamically stable mixtures of oil, water, and surfactant, are known and have been studied for more than 70 years. However, even today there are still quite a number of unclear aspects, and more recent research work has modified and extended our picture. This review gives a short overview of how the understanding of microemulsions has developed, the current view on their properties and structural features, and in particular, how they are related to applications. We also discuss more recent developments regarding nonclassical microemulsions such as surfactant-free (ultraflexible) microemulsions or ones containing uncommon solvents or amphiphiles (like antagonistic salts). These new findings challenge to some extent our previous understanding of microemulsions, which therefore has to be extended to look at the different types of microemulsions in a unified way. In particular, the flexibility of the amphiphilic film is the key property to classify different microemulsion types and their properties in this review. Such a classification of microemulsions requires a thorough determination of their structural properties, and therefore, the experimental methods to determine microemulsion structure and dynamics are reviewed briefly, with a particular emphasis on recent developments in the field of direct imaging by means of electron microscopy. Based on this classification of microemulsions, we then discuss their applications, where the application demands have to be met by the properties of the microemulsion, which in turn are controlled by the flexibility of their amphiphilic interface. Another frequently important aspect for applications is the control of the rheological properties. Normally, microemulsions are low viscous and therefore enhancing viscosity has to be achieved by either having high concentrations (often not wished for) or additives, which do not significantly interfere with the microemulsion. Accordingly, this review gives a comprehensive account of the properties of microemulsions, including most recent developments and bringing them together from a united viewpoint, with an emphasis on how this affects the way of formulating microemulsions for a given application with desired properties.

Contribution of Molecular Dynamics in pNMR for the Structural Determination of AnV and AnVI Complexes in Solution.

In this study, we propose to use classical molecular dynamics (MD) coupled with 1H NMR spectroscopy to study the conformations of different actinyl AnVI (An = U, Np, and Pu) and AnV (An = Np) complexes with tetra-ethyl dyglicolamide (TEDGA) ligands in order to have a better representation of such complexes in solution. Molecular dynamics simulations showed its effectiveness in interpreting the experiments by the calculation of geometric factors needed for the determination of magnetic properties of these complexes. We demonstrated that different conformations of the AnV and AnVI complexes with TEDGA exist in solution with different coordination modes, which is experimentally confirmed by 1H NMR and EXAFS spectroscopies. Furthermore, MD simulations provide additional insights into the structures of complexes in solution since conformations with fast exchanges, which are not accessible from NMR experiments, have been observed by MD simulations.

Why local and non-local terms are essential for second harmonic generation simulation?

Second Harmonic Generation (SHG) today represents one of the most powerful techniques to selectively probe all types of interfaces. However, the origin of the SHG signal at a molecular level is still debated since the local dipole contribution, which is strongly correlated to the molecular orientation can be counterbalanced by non-local quadrupole contributions. Here, we propose a method to simulate the SHG signal of a model water/air interface from the molecular response of each contribution. This method includes both local and non-local terms, which are represented, respectively, by the dependency of the polarisability and hyperpolarisability upon the chemical environment of the molecule and by the bulk quadrupole response. The importance of both terms for the sound simulation of the SHG signals and their interpretation is assessed. We demonstrate that the sole dipole term is unable to simulate a SHG signal, even if the dependency of the hyperpolarisability on the local environment is considered. The inclusion of the bulk quadrupole contribution, which largely dominates the dipole contribution, is essential to predict the SHG response, although the accuracy of the prediction is increased when the dependency upon the local environment is considered.

Liquid/liquid interface in periodic boundary condition.

We study how surface phenomena can change the interface geometry in liquid-liquid two-phase systems with periodic boundary conditions. Without any curvature effect on surface tension, planar (slab), cylindrical, and spherical structures are successively obtained as a function of the total composition and elongation of the box, in accordance with molecular dynamics simulations for a water/heptane system. The curvature effects described by Tolman relationship desymmetrize the phase diagram by stabilizing a concavity but it leads to inconsistencies with high curvature. Helfrich model partially resolves this and predicts the possible presence of shells reflecting a frustrated system.

UO22+ structure in solvent extraction phases resolved at molecular and supramolecular scales: a combined molecular dynamics, EXAFS and SWAXS approach.

A new polarizable force field for describing the solvation of the uranyl (UO22+) cation in solvent extraction phases has been developed for molecular dynamics simulations. The validity of the polarizable force field has been established by comparison with EXAFS and SWAXS experiments. This new force field allows for describing both the UO22+ hydration and solvation properties in good agreement with the experiments. In aqueous phases we demonstrated that the UO22+ force field has been improved from the previous one we developed. Indeed, the UO22+ structural and dynamics properties, i.e., the dynamics of the water molecules in the vicinity of the uranyl cation, calculated from molecular dynamics are in better agreement with the EXAFS experiments. Furthermore, the transferability of the UO22+ force field proposed here has been validated on typical solvent extraction phases containing uranyl nitrate salts with extractant molecules, namely DMDOHEMA molecules, in n-heptane. The good agreements observed between the theoretical (MD simulations) and experimental UO22+ structures at the molecular (EXAFS) and supramolecular (SWAXS) scales prove the accuracy of the UO22+ force field developed and proposed in the present paper.

The role of curvature effects in liquid-liquid extraction: assessing organic phase mesoscopic properties from MD simulations.

The bending rigidity of small reverse aggregates involved in liquid-liquid extraction processes has been investigated by molecular dynamics simulations. Simulations of a common extractant (DMDOHEMA) with four hydrophobic chains in explicit solvent (n-heptane) and in vacuum have been performed to determine the effect of solvent penetration on film stiffness. Elastic film bending energy that is needed for mesoscopic modelling of transfer of species between complex fluids is harmonic in terms of curvature (Helfrich formalism) and the packing parameter only if the solvent is explicitly taken into account. In terms of the packing parameter of the real molecular film constituting the reverse water in oil aggregates and taking into account molecular volume, area and film thickness (that is in agreement with Tanford's model), the bending rigidity is calculated to be about 16 kBT per extractant molecule (about 40 kJ mol-1), which is smaller than the free energy of transfer from an isolated "monomer" molecule to a weak aggregate, but of the order of magnitude of the free energy of transfer used in liquid-liquid extraction processes.

A Combined Spectroscopic/Molecular Dynamic Study for Investigating a Methyl-Carboxylated PEI as a Potential Uranium Decorporation Agent.

Natural uranium has a very limited radioactive dose impact, but its chemical toxicity due to chronic exposure is still a matter of debate. Once inside the human body, the soluble uranium, under its uranyl form (U(VI)), is quickly removed from the blood system, partially excreted from the body, and partially retained in targeted organs, that is, the kidneys and bone matrix essentially. It is then crucial to remove or prevent the incorporation of uranium in these organs to limit the long-term chronic exposure. A lot of small chelating agents such as aminocarboxylates, catecholamides, and hydroxypyridonates have been developed so far. However, they suffer from poor selectivity and targeting abilities. Macromolecules and polymers are known to present a passive accumulation (size related), that is, the so-called enhanced permeability and retention effect, toward the main organs, which can be used as indirect targeting. Very interestingly, the methyl carboxylated polyethylenimine (PEI-MC) derivative has been described as a potent sequestering agent for heavy metals. It would be therefore an interesting candidate to evaluate as a new class of decorporation agents with passive targeting capabilities matching uranium preferential sequestering sites. In the present work, we explored the ability of a highly functionalized (89% rate) PEI-MC to uptake U(VI) close to physiological pH using a combination of analytical and spectroscopic techniques (inductively coupled plasma optical emission spectrometry (ICP-OES); extended X-ray absorption fine structure (EXAFS); and Fourier transformed infrared (FT-IR)) together with molecular dynamics (MD) simulation. A maximum loading of 0.47 mg U(VI) per milligram of PEI-MC was determined by ICP-OES measurements. From FT-IR data, a majority of monodentate coordination of the carboxylate functions of the PEI-MC seems to occur. From EXAFS and MD, a mix of mono and bidentate coordination mode was observed. Note that agreement between the EXAFS metrical parameters and MD radial distribution functions is remarkable. To the best of our knowledge, this is the first comprehensive structural study of a macromolecular PEI-based agent considered for uranium decorporation purposes.

Stability of reverse micelles in rare-earth separation: a chemical model based on a molecular approach.

Molecular complexes formed in the organic phase during solvent extraction may self-assemble as reverse micelles, and therefore induce a supramolecular organization of this phase. In most of the cases, water molecules play an essential role in the organization of this non polar medium. The aim of this work is to investigate the speciation of the aggregates formed in the organic phase during solvent extraction, and especially to assess their stability as a function of the number of water molecules included in their polar core. We have focused on malonamide extractants that have already been investigated experimentally. Different stoichiometries of reverse micelles in the organic phase have been studied by means of classical molecular dynamics simulations. Furthermore, umbrella-sampling molecular dynamics simulations have been used to calculate the equilibrium constant (K°) representing the association/dissociation pathways of water molecules in the aggregates and the corresponding reaction free energies (ΔrG°).

Multi-scale modelling of uranyl chloride solutions.

Classical molecular dynamics simulations with explicit polarization have been successfully used to determine the structural and thermodynamic properties of binary aqueous solutions of uranyl chloride (UO2Cl2). Concentrated aqueous solutions of uranyl chloride have been studied to determine the hydration properties and the ion-ion interactions. The bond distances and the coordination number of the hydrated uranyl are in good agreement with available experimental data. Two stable positions of chloride in the second hydration shell of uranyl have been identified. The UO2(2+)-Cl(-) association constants have also been calculated using a multi-scale approach. First, the ion-ion potential averaged over the solvent configurations at infinite dilution (McMillan-Mayer potential) was calculated to establish the dissociation/association processes of UO2 (2+)-Cl(-) ion pairs in aqueous solution. Then, the association constant was calculated from this potential. The value we obtained for the association constant is in good agreement with the experimental result (KUO2Cl(+) = 1.48 l mol(-1)), but the resulting activity coefficient appears to be too low at molar concentration.

Predicting for thermodynamic instabilities in water/oil/surfactant microemulsions: a mesoscopic modelling approach.

The thermodynamics and structural properties of flexible and rigid nonionic water/oil/surfactant microemulsions have been investigated using a two level-cut Gaussian random field method based on the Helfrich formalism. Ternary stability diagrams and scattering spectra have been calculated for different surfactant rigidities and spontaneous curvatures. A more important contribution of the Gaussian elastic constants compared to the bending one is observed on the ternary stability diagrams. Furthermore, influence of the spontaneous curvature of the surfactant points out a displacement of the instability domains which corresponds to the difference between the spontaneous and effective curvatures. We enlighten that a continuous transition from a connected water in oil droplets to a frustrated locally lamellar (oil in water in oil droplets) microstructure is found to occur when increasing the temperature for an oil-rich microemulsion. This continuous transition translated in a shift in the scattering functions, points out that the phase inversion phenomenon occurs by a coalescence of the water droplets.

Modeling Lanthanide Ions in Solution: A Versatile Force Field in Aqueous and Organic Solvents.

In this paper, we propose a new nonpolarizable force field for describing the Ln3+ (Ln = lanthanide) series based on a 12-6-4 Lennard-Jones potential. The development of the force field was performed in pure water by adjusting both the ion-oxygen distance and the hydration free energy. This force field accurately reproduces the Ln3+ hydration properties through the series, especially the coordination number that is hardly accessible using a nonpolarizable force field. Then, the validity and the transferability of the current force field were evaluated for two different systems containing Ln3+ in various solvents, namely, 0.1 mol L-1 La(NO3)3 salts in methanol and Eu(NO3)3 salts in solvent organic phases composed of DMDOHEMA molecules in n-heptane. The good agreement between our simulations and the data available in the literature confirms the accuracy of the force field for describing the lanthanide cations in both aqueous and nonaqueous media.

Mesoscopic modelling of frustration in microemulsions.

The swelling behaviour of water-oil microemulsions - considering a surfactant layer between oil and water - has been studied using a two level-cuts Gaussian random field approach based on the Helfrich formalism. Microstructures and scattering properties of microemulsions have been calculated for different amounts of oil (and water) for flexible and rigid microemulsions. When the stiffness, the spontaneous curvature of the interfacial film, and the surface to volume ratio of the immiscible fluids are varied, the microemulsion topology and morphology change in order to minimize the microemulsion free energy. Our simulations point out a change in the microemulsion morphology as a function of the surfactant film rigidity and the composition of oil, water and the surfactant. Locally lamellar structures are found for rigid microemulsions, whereas for more flexible ones, the connected-droplet and/or bicontinuous structures are preferred. Furthermore, we show that the microemulsion swelling versus the volume fraction gives a specific signature of the microemulsion microstructure. This allows for discriminating between different types of microemulsions: flexible, frustrated and unfrustrated (close to bi-liquid foams), and connected structures as molten hexagonal and cubic phases. The universal swelling behaviour is compared to different analytic expressions of Disordered Open Connected (DOC) models for the microemulsion swelling versus the volume fraction.

Exploring the Aggregation Behavior of Extractant Molecules in Ionic Liquids: A Coupled Polarizable Molecular Dynamics and SAXS Study.

This study presents a comprehensive investigation of the aggregation behavior of a malonamide extractant molecule (N,N'-dimethyl,N,N'-dioctylhexylethoxymalonamide (DMDOHEMA)) in three different solvents, including two piperidinium- and (trifluoromethylsulfonyl)imide-based ionic liquids (1-ethyl-1-butylpiperidinium bis(trifluoromethylsulfonyl)imide ([EBPip+][NTf2-]) and 1-ethyl-1-octylpiperidinium bis(trifluoromethylsulfonyl)imide ([EOPip+][NTf2-])) and n-dodecane. By combining polarizable molecular dynamics simulations and small-angle X-ray scattering experiments, we extensively investigated the arrangement of supramolecular assemblies of the extractant molecules. Our results showed that the insertion of the alkyl chains of the extractant molecules into the apolar domain of [EOPip+][NTf2-] has a significant impact on the aggregation behavior of the extractant molecules, leading to the formation of smaller aggregates having a higher dispersion compared to other solvents. These findings provide new insights into the physicochemical properties of this type of system and are crucial in designing more effective solvents for rare earth metal extraction.

Long-Range Organization Study of Piperidinium-Based Ionic Liquids by Polarizable Molecular Dynamics Simulations.

The nanoscale organization of some classes of ionic liquids is responsible for their singular properties. In this paper, we use polarizable molecular dynamics simulations and small-angle X-ray scattering to probe the structure of two piperidinium- and (trifluoromethylsulfonyl)imide-based ionic liquids ([EBPip+][NTf2-] and [EOPip+][NTf2-]) that differ in the alkyl chain length of their cation. The X-ray scattering intensities calculated numerically, from the radial distribution functions, are in excellent agreement with the experimental data. The analysis of the different contributions of the X-ray scattering data allowed us to highlight the correlations responsible for the low q peak observed for the long-chain alkyl cations. New angular analyses showed that anions were more likely to align with alkyl chains as their size increased, inducing angular correlation between anions at larger distances. They also showed that the long alkyl chains of the cations aligned more with each other than the short ones. These more aligned alkyl chains induce a smaller volume of the apolar microdomains compared to the well-studied imidazolium-based ionic liquids, leading to the smaller correlation distance for piperidinium-based ionic liquids.

Revised ionic radii of lanthanoid(III) ions in aqueous solution.

A new set of ionic radii in aqueous solution has been derived for lanthanoid(III) cations starting from a very accurate experimental determination of the ion-water distances obtained from extended X-ray absorption fine structure (EXAFS) data. At variance with previous results, a very regular trend has been obtained, as expected for this series of elements. A general procedure to compute ionic radii in solution by combining the EXAFS technique and molecular dynamics (MD) structural data has been developed. This method can be applied to other ions allowing one to determine ionic radii in solution with an accuracy comparable to that of the Shannon crystal ionic radii.

Complexation of lanthanides(III), americium(III), and uranium(VI) with bitopic N,O ligands: an experimental and theoretical study.

New functionalized terpyridine-diamide ligands were recently developed for the group actinide separation by solvent extraction. In order to acquire a better understanding of their coordination mode in solution, protonation and complexation of lanthanides(III), americium(III), and uranium(VI) with these bitopic N,O-bearing ligands were studied in homogeneous methanol/water conditions by experimental and theoretical approaches. UV-visible spectrophotometry was used to determine the protonation and stability constants of te-tpyda and dedp-tpyda. The conformations of free and protonated forms of te-tpyda were investigated using NMR and theoretical calculations. The introduction of amide functional groups on the terpyridine moiety improved the extracting properties of these new ligands by lowering their basicity and enhancing the stability of the corresponding 1:1 complexes with lanthanides(III). Coordination of these ligands was studied by density functional theory and molecular dynamics calculations, especially to evaluate potential participation of hard oxygen and soft nitrogen atoms in actinide coordination and to correlate with their affinity and selectivity. Two predominant inner-sphere coordination modes were found from the calculations: one mode where the cation is coordinated by the nitrogen atoms of the cavity and by the amide oxygen atoms and the other mode where the cation is only coordinated by the two amide oxygen atoms and by solvent molecules. Further simulations and analysis of UV-visible spectra using both coordination modes indicate that inner-sphere coordination with direct complexation of the three nitrogen and two oxygen atoms to the cation leads to the most likely species in a methanol/water solution.

Understanding the nitrate coordination to Eu3+ ions in solution by potential of mean force calculations.

Coordination of nitrate anions with lanthanoid cations (Ln(3+)) in water, methanol and octanol-1 has been studied by means of molecular dynamics simulations with explicit polarization. Potential of mean force (PMF) profiles have been calculated for a mono-complex of lanthanoid nitrate (Ln(NO(3))(2+)) in these solvents using umbrella-sampling molecular dynamics. In pure water, no difference in the nitrato coordination to lanthanoids (Nd(3+), Eu(3+) and Dy(3+)) is observed, i.e. the nitrate anion prefers the monodentate coordination, which promotes the salt dissociation. Then, the influence of the nature of the solvating molecules on the nitrato coordination to Eu(3+) has been investigated. PMF profiles point out that both monodenate and bidentate coordinations are stable in neat methanol, while in neat octanol, only the bidentate one is. MD simulations of Eu(NO(3))(3) in water-octanol mixtures with different concentrations of water have been then performed and confirm the importance of the water molecules' presence on the nitrate ion's coordination mode.

Polarizable interaction potential for molecular dynamics simulations of actinoids(III) in liquid water.

In this work, we have developed a polarizable classical interaction potential to study actinoids(III) in liquid water. This potential has the same analytical form as was recently used for lanthanoid(III) hydration [M. Duvail, P. Vitorge, and R. Spezia, J. Chem. Phys. 130, 104501 (2009)]. The hydration structure obtained with this potential is in good agreement with the experimentally measured ion-water distances and coordination numbers for the first half of the actinoid series. In particular, the almost linearly decreasing water-ion distance found experimentally is replicated within the calculations, in agreement with the actinoid contraction behavior. We also studied the hydration of the last part of the series, for which no structural experimental data are available, which allows us to provide some predictive insights on these ions. In particular we found that the ion-water distance decreases almost linearly across the series with a smooth decrease of coordination number from nine to eight at the end.

Thermodynamics of Malonamide Aggregation Deduced from Molecular Dynamics Simulations.

The aggregation of malonamide extractants diluted in an aliphatic solvent phase has been studied in the presence of water by molecular dynamics simulation. Using association criteria based on distances between molecules and graphs theory, the aggregate distribution has been computed and the corresponding Gibbs energy of aggregates and mass action law constants have been determined. Finally, a model allowing us to the compute critical micelle concentration and osmotic data for a variable concentration of extractants, with or without a correction of the organic phase activity, was developed. It appears however that the accurate depiction of the aggregation allows modeling the thermodynamics of the solution even without an explicit calculation of the activity: both models give results in good agreement with the experiments.

Deciphering second harmonic generation signals.

Second harmonic generation (SHG) has emerged as one of the most powerful techniques used to selectively monitor surface dynamics and reactions for all types of interfaces as well as for imaging non-centrosymmetric structures, although the molecular origin of the SHG signal is still poorly understood. Here, we present a breakthrough approach to predict and interpret the SHG signal at the atomic level, which is freed from the hyperpolarisability concept and self-consistently considers the non-locality and the coupling with the environment. The direct ab initio method developed here shows that a bulk quadrupole contribution significantly overwhelms the interface dipole term in the purely interfacial induced second-order polarisation for water/air interfaces. The obtained simulated SHG responses are in unprecedented agreement with the experimental signal. This work not only paves the road for the prediction of SHG response from more complex interfaces of all types, but also suggests new insights in the interpretation of the SHG signal at a molecular level. In particular, it highlights the modest influence of the molecular orientation and the high significance of the bulk quadrupole contribution, which does not depend on the interface, in the total experimental response.