No abstract

A simple rigid three-site model for methanol compatible with the simple point charge ͑SPC͒ water and the GROMOS96 force field is parametrized and tested. The influence of different force-field parameters, such as the methanol geometry and the charge distribution on several properties calculated by molecular dynamics is investigated. In particular an attempt was made to obtain good agreement with experimental data for the static dielectric constant and the mixing enthalpy with water. The model is compared to other methanol models from the literature in terms of the ability to reproduce a range of experimental properties.

We report a molecular dynamics study of the interface between water and (macroscopically) water-immiscible room-temperature ionic liquids "ILs", composed of PF6(-) anions and butyl- versus octyl-substituted methylimidazolium+ cations (noted BMI+ and OMI+). Because the parameters used to simulate the pure ILs were found to exaggerate the water/IL mixing, they have been modified by scaling down the atomic charges, leading to better agreement with the experiment. The comparison of [OMI][PF6] versus [BMI][PF6] ILs demonstrates the importance of the N-alkyl substituent on the extent of solvent mixing and on the nature of the interface. With the most hydrophobic [OMI][PF6] liquid, the "bulk" IL phase is dryer than with the [BMI][PF6] liquid. At the interface, the OMI+ cations retain direct contacts with the bulk IL, whereas the more hydrophilic PF6(-) anions gradually dilute in the local water micro-environment and are thus isolated from the "bulk" IL. The interfacial OMI+ cations are ordered with their imidazolium moiety pointing toward the aqueous side and their octyl chains toward the IL side of the interface. With the [BMI][PF6] liquid, the system gradually evolves from an IL-rich to a water-rich medium, leading to an ill-defined interfacial domain with high intersolvent mixing. As a result, the BMI+ cations are isotropically oriented "at the interface". Because the imidazolium cations are more hydrophobic than the PF6(-) anions, the charge distribution at the interface is heterogeneous, leading to a positive electrostatic potential at the interface with the two studied ILs. Mixing-demixing simulations on [BMI][PF6]/water mixtures are also reported, comparing Ewald versus reaction field treatments of electrostatics. Phase separation is very slow (at least 30 ns), in marked contrast with mixtures involving classical organic liquids, which separate in less than 0.5 ns at the microscopic level. The results allow us to better understand the specificity of the aqueous interfaces with hydrophobic ionic liquids, compared with classical organic solvents, which has important implications as far as the mechanism of liquid-liquid ion extraction is concerned.

Density-functional-based Car-Parrinello molecular dynamics (CPMD) simulations have been performed for the ionic liquid 1,3-dimethylimidazolium chloride, [dmim]Cl, at 438 K. The local structure of the liquid is described in terms of various partial radial distribution functions and anisotropic spatial distributions, which reveal a significant extent of hydrogen bonding. The cation-anion distribution simulated with the BP86 functional is in qualitative agreement with the structural model derived from neutron diffraction data for the liquid, whereas the theoretical cation-cation distribution shows less satisfactory accord. Population analyses indicate noticeable charge transfer from anions to cations, and specific CH...Cl hydrogen bonds are characterized in terms of donor-acceptor interactions between lone pairs on Cl and antibonding sigma(CH) orbitals.

No abstract

We report molecular dynamics (MD) simulations of the aqueous interface of the hydrophobic [BMI][Tf2N] ionic liquid (IL), composed of 1-butyl-3-methylimidazolium cations (BMI+) and bis(trifluoromethylsulfonyl)imide anions (Tf2N-). The questions of water/IL phase separation and properties of the neat interface are addressed, comparing different liquid models (TIP3P vs TIP5P water and +1.0/-1.0 vs +0.9/-0.9 charged IL ions), the Ewald vs the reaction field treatments of the long range electrostatics, and different starting conditions. With the different models, the "randomly" mixed liquids separate much more slowly (in 20 to 40 ns) than classical water-oil mixtures do (typically, in less than 1 ns), finally leading to distinct nanoscopic phases separated by an interface, as in simulations which started with a preformed interface, but the IL phase is more humid. The final state of water in the IL thus depends on the protocol and relates to IL heterogeneities and viscosity. Water mainly fluctuates in hydrophilic basins (rich in O(Tf2N) and aromatic CH(BMI) groups), separated by more hydrophobic domains (rich in CF3(Tf2N) and alkyl(BMI) groups), in the form of monomers and dimers in the weakly humid IL phase, and as higher aggregates when the IL phase is more humid. There is more water in the IL than IL in water, to different extents, depending on the model. The interface is sharper and narrower (approximately 10 A) than with the less hydrophobic [BMI][PF6] IL and is overall neutral, with isotropically oriented molecules, as in the bulk phases. The results allow us to better understand the analogies and differences of aqueous interfaces with hydrophobic (but hygroscopic) ILs, compared to classical organic liquids.

No abstract

We report a series of molecular dynamics simulations on the demixing of “homogeneous” binary water−chloroform mixtures containing species involved in the assisted ion extraction process. We consider an ionophore L (L = 1,3-alternate calix4arene-crown6), uncomplexed salts of Cs+ and the LCs+ and LNa+ cation complexes with a lipophilic (Pic-) and a hydrophilic (Cl-) counterion, respectively, as being solutes. In all cases, the liquids separate rapidly, leading to two solvent slabs separated by a well-defined interface. However, the final state is very different, depending on the hydrophilic/hydrophobic balance of the solutes:  the Cs+ and NO3 - ions of the CsNO3 salt are completely immersed in the aqueous phase, whereas Pic- anions display a strong adsorption at the interface. The LCs+ complex and the free ligand L, although more soluble in the organic phase than in water, also display a surfactant like behavior. Similar conclusions are obtained when L, LCs+, Cs+ Pic-, and Cs+ NO3 - ions are simultaneously present in the solution. On the basis of free energy perturbation calculations on LM+ complexes, we calculate a marked Cs+/Na+ recognition by L at the interface. These results have important implications concerning the mechanism of ionophore assisted liquid−liquid ion extraction and recognition processes at the interface.