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.

Coalescence Initiation in Liquid-Liquid Systems: A Stochastic Model.

Coalescence is a complex phenomenon leading to the merging of deformable particles of fluid. The complexity stems largely from the simultaneous occurrence of phenomena of a different nature (hydrodynamic, electrostatic, physicochemical) acting at different scales. The stochastic effects controlling the formation of the liquid bridge between two droplets of the same liquid, immersed in another nonmiscible liquid, are studied through a series of molecular dynamics simulations. The case of heptane droplets in water, relevant to solvent extraction, a key process of the circular economy, is considered. From this series of simulations, we have confirmed that the probability function of coalescence of two identical droplets in contact follows a Poisson distribution. We moreover propose a criterion for the initiation of coalescence based on nucleation theory. A complete description of the stochastic initiation of coalescence is hence provided, opening many perspectives for the simulation of coalescence in continuous approaches used in fluid mechanics and chemical engineering. The methodology can be generalized to droplets of different size and composition, immersed in gas, or to bubbles, i.e., to other physical problems whose kinetics is influenced by the molecular scale.

Cation Adsorption in TiO2 Nanotubes: Implication for Water Decontamination.

TiO2 nanotubes constitute very promising nanomaterials for water decontamination by the removal of cations. We combined a range of experimental techniques from structural analyses to measurements of the properties of aqueous suspensions of nanotubes, with (i) continuous solvent modeling and (ii) quantum DFT-based simulations to assess the adsorption of Cs+ on TiO2 nanotubes and to predict the separation of metal ions. The methodology is set to be operable under realistic conditions, which, in this case, include the presence of CO2 that needs to be treated as a substantial contaminant, both in experiments and in models. The mesoscopic model, based on the Poisson-Boltzmann equation and surface adsorption equilibrium, predicts that H+ ions are the charge-determining species, while Cs+ ions are in the diffuse layer of the outer surface with a significant contribution only at high concentrations and high pH. The effect of the size of nanotubes in terms of the polydispersity and the distribution of the inner and outer radii is shown to be a third-order effect that is very small when the nanotube layer is not very thick (ranging from 1 to 2 nm). Besides, DFT-based molecular dynamics simulations demonstrate that, for protonation, the one-site and successive association assumption is correct, while, for Cs+ adsorption, the size of the cation is important and the adsorption sites should be carefully defined.

When Ions Defy Electrostatics: The Case of Superchaotropic Nanoion Adsorption.

Nanometer-sized anions, like polyoxometalates and borate clusters, bind to nonionic hydrated matter driven by the chaotropic effect, which arises from the favorable dehydration of the ions. Herein, we evaluate the adsorption and activity coefficient of the superchaotropic Keggin polyoxometalate SiW12O404- (SiW) on nonionic surfactant (C8E4) micelles by modeling small-angle X-ray and neutron-scattering spectra. Neither hard sphere nor electrostatic repulsion models reproduce the experimental activity coefficient of adsorbed SiW ions on the micelles. However, the activity and binding of SiW on the micelles is well-described by a Langmuir adsorption isotherm. These results imply that adsorbed SiW ions are non-interacting and "create" around themselves adsorption sites on the micelle. The temperature dependence of the adsorption constant showed that the SiW adsorption is enthalpically driven and entropically unfavorable, in line with the typical chaotropic thermochemical signature. The adsorption enthalpy can be split into an electrostatic term and a water-recovery term to evaluate and qualitatively predict the superchaotropicity of a nanoion.

Liquid-Liquid Flow at Nanoscale: Slip and Hydrodynamic Boundary Conditions.

Nonequilibrium molecular dynamics (NEMD) simulations have been performed to describe the flow of a fluid nanolayer confined by another fluid. The results show that the behavior of liquids can still be described by the Navier-Stokes equation at the nanoscale, i.e., when only few molecular layers are involved. NEMD furthermore gives additional knowledge on flow. Indeed, while a very small slip is evidenced for a solid-liquid interface as, e.g., in lubrication, the slip lengths are significantly larger at the liquid-liquid interface, as encountered, e.g., in droplet coalescence. The slip lengths of the two fluids are linked. The increase in hydrodynamic slip for liquid-liquid interfaces is attributed to the enhancement of fluid diffusion, which reduces friction.

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.

How Temperature Rise Can Induce Phase Separation in Aqueous Biphasic Solutions.

Ionic-liquid-based acidic aqueous biphasic solutions (AcABSs) recently offered a breakthrough in the field of metal recycling. The particular mixture of tributyltetradecylphosphonium chloride ([P4,4,4,14]Cl), acid, and water presents the unusual characteristic of a lower solution critical temperature (LCST), leading to phase separation upon a temperature rise of typically a few tens of degrees. We address here the microscopic mechanisms driving the phase separation. Using small-angle neutron scattering, we characterized the spherical micelle formation in a binary ionic liquid/water solution and the micelle aggregation upon the addition of acid due to the screening of electrostatic repulsion. The increase in both the acid concentration and the temperature eventually leads to micelle flocculation and phase separation. This last step is achieved through chloride ion adsorption at the surface of the micelle. This exothermic adsorption compensates for the entropic cost, leading to a counterintuitive behavior, and may be generalized to a number of molecular systems with an LCST.

Experimentally probing ionic solutions in single-digit nanoconfinement.

Understanding ionic solutions in single-digit nanoconfinement is crucial to explain the behavioral transition of confined solutions. This is particularly the case when the system length scale crosses the classical key length scales describing energetics and equilibrium of ionic solutions next to surfaces. Experimentally probing nanoconfinement would open large perspectives to test modelling or theory predictions. Here, using a new test vehicle that consists in 3 and 5 nm-height silica nanochannels associated with an original characterization technique based on the interface hard X-ray reflectivity analysis, we directly probed the transport of solutions containing cations having increasing kosmotropic properties (XCl2 with X: Ba < Ca < Mg) and obtained their distributions inside the nanochannels. We observed that cation adsorption decreases with the size of the confinement and that small cation adsorption is favored. In addition, nanochannel clogging occurs when ions tend to form ion pairs. These ion pairs may play the role of nano-sized prenucleation clusters leading to phase precipitation. These results evidence the specific ion effect in single-digit nanoconfinement that may result in dramatic changes of solution properties. In this line, our new method opens new perspectives for the characterization of ionic solutions and of interfaces in single-digit nanoconfinement.

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.

Theory of Ternary Fluids under Centrifugal Fields.

In this work, we developed a general theoretical description of ternary solutions of small molecules under a centrifugal field, from which we obtained the centrifugation map (CMap) as a general tool to understand observations or to predict composition profiles in centrifugal fields of arbitrary strength. The theoretical formalism is based on the classical density functional theory with established models for liquid mixtures. Thermodynamics also yields a general criterion for apparent aggregation. The strength of the CMap approach is illustrated for a ternary model system where ethanol is a co-solvent.

Molecular Forces in Liquid-Liquid Extraction.

The phase transfer of ions is driven by gradients of chemical potentials rather than concentrations alone (i.e., by both the molecular forces and entropy). Extraction is a combination of high-energy interactions that correspond to short-range forces in the first solvation shell such as ion pairing or complexation forces, with supramolecular and nanoscale organization. While the latter are similar to the long-range solvent-averaged interactions in the colloidal world, in solvent extraction they are associated with lower characteristic lengths of the nanometric domain. Modeling of such complex systems is especially complicated because the two domains are coupled, whereas the resulting free energy of extraction is around kBT to guarantee the reversibility of the practical process. Nevertheless, quantification is possible by considering a partitioning of space among the polar cores, interfacial film, and solvent. The resulting free energy of transfer can be rationalized by utilizing a combination of terms which represent strong complexation energies, counterbalanced by various entropic effects and the confinement of polar solutes in nanodomains dispersed in the diluent, together with interfacial extractant terms. We describe here this ienaics approach in the context of solvent extraction systems; it can also be applied to further complex ionic systems, such as membranes and biological interfaces.

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.

Phase separation of binary mixtures induced by soft centrifugal fields.

We use the model system ethanol-dodecane to demonstrate that giant critical fluctuations induced by easily accessible weak centrifugal fields as low as 2000g can be observed above the miscibility gap of a binary liquid mixture. Moreover, several degrees above the phase transition, i.e. in the one-phase region, strong gradients of ethanol concentration occur upon centrifugation. In this case, the standard interpretation of sedimentation equilibrium in the analytical ultracentrifuge (AUC) yields an apparent molar mass of ethanol three orders of magnitude higher than the real value. Notably, these composition gradients have no influence on the distribution gradient of solutes such as dyes like Nile red. The thick opaque interphase formed upon centrifugation does not appear as the commonly observed sharp meniscus, but as a turbidity zone, similar to critical opalescence. This layer is a few millimeters thick and separates two fluids with low compositional gradients. All these effects can be qualitatively understood and explained using the Flory-Huggins solution model coupled to classical density functional theory (DFT). In this domain hetero-phase fluctuations can be triggered by gravity even far from the critical point. Taking into account Jean Perrin's approach to external fields in colloids, a self-consistent definition of the Flory effective volume and an explicit calculation of the total free energy per unit volume is possible.

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.

How Adsorption of Pheromones on Aerosols Controls Their Transport.

We propose a general transport theory for pheromone molecules in an atmosphere containing aerosols. Many pheromones are hydrophobic molecules containing polar groups. They are low volatile and have some properties similar to those of hydrotropes. They therefore form a nonsoluble film at the water-air interface of aerosols. The fate of a small pheromone puff in air is computed through reaction-diffusion equations. Partitioning of pheromones between the gas and the aerosol surface over time is studied for various climate conditions (available aerosol surface) and adsorption affinities (energy of adsorption). We show that, for adsorption energy above 30 k B T per molecule, transport of pheromones on aerosols dominates over molecular transport typically 10 s after pheromone emission, even when few adsorbing aerosols are present. This new communication path for airborne chemicals leads to distinctive features including enhanced signal sensibility and increased persistence of pheromone concentration in the air due to slow diffusion of aerosols. Each aerosol droplet has the ability to adsorb thousands of pheromones to the surface, keeping a "history" of the atmospheric content between emission and reception. This new mechanism of pheromone transport leads to dramatic consequences on insect sensing revisiting the way we figure the capture of chemical signals.

Deformation of a Liquid Near an AFM Tip: Molecular Dynamics Approach.

The interaction between an atomic force microscopy (AFM) probe and a thin film of water deposited over a flat substrate is studied using molecular dynamics (MD). The effects of the film thickness and the probe radius on both the deformation height of the liquid interface and the distance of the jump to contact at which the liquid comes in direct contact with the probe are investigated. The dynamics of the surface deformation and the role of interface fluctuations are studied in detail. The systems considered belong to the thin-film regime described in a semianalytical model previously established by Ledesma-Alonso et al. (Langmuir 2013, 29, 7749-7757). MD simulations predict that for shallow films, both the distance at which the jump to contact occurs and the surface maximal deformation height increase steadily with the layer thickness regardless of the probe radius, which is in agreement with the previously proposed theoretical model. The deformation of the surface was shown to be unstable because of the strong effect of thermal fluctuations. For each of the considered systems, the film thickness was such that interface fluctuations induced the jump to contact. The comparison of the deformation obtained in MD with the profiles predicted by the continuous model points out the complementarity between the two approaches. The results of the molecular approach not only are consistent with those of the continuous model but also provide more information on the description of nanoscale phenomena. In particular, MD results point out the importance of fluctuations when it comes to the description of the particular dynamics of nanosystems involving soft interfaces. This shows the need to improve continuous models by complementing them with a molecular approach for a better accuracy.

Synergistic Solvent Extraction Is Driven by Entropy.

In solvent extraction, the self-assembly of amphiphilic molecules into an organized structure is the phenomenon responsible for the transfer of the metal ion from the aqueous phase to the organic solvent. Despite their significance for chemical engineering and separation science, the forces driving the solute transfer are not fully understood. Instead of assuming the simple complexation reaction with predefined stoichiometry, we model synergistic extraction systems by a colloidal approach that explicitly takes into account the self-assembly resulting from the amphiphilic nature of the extractants. Contrary to the current paradigm of simple stoichiometry behind liquid-liquid extraction, there is a severe polydispersity of aggregates completely different in compositions, but similar in the free energy. This variety of structures on the nanoscale is responsible for the synergistic transfer of ions to the organic phase. Synergy can be understood as a reciprocal effect of chelation: it enhances extraction because it increases the configurational entropy of an extracted ion. The global overview of the complex nature of a synergistic mixture shows different regimes in self-assembly, and thus in the extraction efficiency, which can be tuned with respect to the green chemistry aspect.

Investigation of the Structure of Concentrated NaOH Aqueous Solutions by Combining Molecular Dynamics and Wide-Angle X-ray Scattering.

Classical molecular dynamics has been performed with explicit polarization on NaOH aqueous solutions from 0.5 mol L-1 up to 9.7 mol L-1. We adapted a force field of OH- for polarizable simulation in order to reproduce the NaOH structural and thermodynamics properties in aqueous solutions. A good agreement between theoretical and experimental results has been found. Wide-angle X-ray scattering (WAXS) intensities issued from molecular dynamics are compared to experimental ones measured on Synchrotron facilities. The structure of the first coordination shell of Na+ has been studied to determine the variation of the oxygen number and hydroxide oxygen around the cation. In addition, Na+-OH- McMillan-Mayer potential issued from molecular dynamics simulations has been calculated and allows for calculating Na+-OH- pair association constant of 0.1 L mol-1, which is in good agreement with the experiments.

Colloidal Model for the Prediction of the Extraction of Rare Earths Assisted by the Acidic Extractant.

We propose the statistical thermodynamic model for the prediction of the liquid-liquid extraction efficiency in the case of rare-earth metal cations using the common bis(2-ethyl-hexyl)phosphoric acid (HDEHP) extractant. In this soft matter-based approach, the solutes are modeled as colloids. The leading terms in free-energy representation account for: the complexation, the formation of a highly curved extractant film, lateral interactions between the different extractant head groups in the film, configurational entropy of ions and water molecules, the dimerization, and the acidity of the HDEHP extractant. We provided a full framework for the multicomponent study of extraction systems. By taking into account these different contributions, we are able to establish the relation between the extraction and general complexation at any pH in the system. This further allowed us to rationalize the well-defined optimum in the extraction engineering design. Calculations show that there are multiple extraction regimes even in the case of lanthanide/acid system only. Each of these regimes is controlled by the formation of different species in the solvent phase, ranging from multiple metal cation-filled aggregates (at the low acid concentrations in the aqueous phase), to the pure acid-filled aggregates (at the high acid concentrations in the aqueous phase). These results are contrary to a long-standing opinion that liquid-liquid extraction can be modeled with only a few species. Therefore, a traditional multiple equilibria approach is abandoned in favor of polydisperse spherical aggregate formations, which are in dynamic equilibrium.