Solvation of sodium chloride in the 1-butyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide ionic liquid: a molecular dynamics study

We report molecular dynamics studies on the solvation of sodium chloride in the 1-butyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide ionic liquid ([BMI][Tf2N] IL). We first consider the potential of mean force for dissociating a single Na+Cl- ion pair, showing that the latter prefers to be undissociated rather than dissociated (by ca. 9 kcal/mol), with a free energy barrier of ca. 5 kcal/mol (at d approximately 5.2 A) for the association process. The preference for Na+Cl- association is also observed from a 100 ns molecular dynamics simulation of a concentrated solution, where the Na+Cl- ions tend to form oligomers and microcrystals in the IL. Conversely, the simulation of Na13Cl14- and Na14Cl13+ cubic microcrystals (with, respectively, Cl- and Na+ at the vertices) does not lead to dissolution in the IL. Among these, Na14Cl13+ is found to be better solvated than Na13Cl14-, mainly due to the stronger Na+...Tf2N- interactions as compared to the Cl-...BMI+ interactions at the vertices of the cube. We finally consider the solid/liquid interface between the 100 face of NaCl and the IL, revealing that, in spite of its polar nature, the crystal surface is solvated by the less polar IL components (CF3(Tf2N) and butyl(BMI) groups) rather than by the polar ones (O(Tf2N) and imidazolium(BMI) ring). Specific ordering at the interface is described for both Tf2N- anions and BMI+ cations. In the first IL layer, the ions are rather parallel to the surface, whereas in the second "layer" they are more perpendicular. A similar IL structure is found at the surface of the all-neutral Na0Cl0 solid analogue, confirming that the solvation of the crystal is rather "apolar", due to the mismatch between the IL and the crystal ions. Several comparisons with water, methanol, or different BMI+-based ILs as solvents are presented, allowing us to better understand the specificity of the ionic liquid-NaCl interactions.     

Many-body potentials for aqueous Li(+), Na(+), Mg(2+), and Al(3+): comparison of effective three-body potentials and polarizable models

Many-body potentials for the aqueous Li(+), Na(+), Mg(2+), and Al(3+) ions have been constructed from ab initio cluster calculations. Pure pair, effective pair, effective three-body, and effective polarizable models were created and used in subsequent molecular dynamics simulations. The structures of the first and second solvation shells were studied using radial distribution functions and angular-radial distribution functions. The effective three-body and polarizable potentials yield similar first-shell structures, while the contraction of the O-O distances between the first and second solvation shells is more pronounced with the polarizable potentials. The definition of the tilt angle of the water molecules around the ions is discussed. When a proper definition is used, it is found that for Li(+), Mg(2+), and Al(3+) the water molecules prefer a trigonal orientation, but for Na(+) a tetrahedral orientation (ion in lone-pair direction) is preferred. The self-diffusion coefficients for the water molecules and the ions were calculated; the ionic values follow the order obtained from experiment, although the simulated absolute values are smaller than experiment for Mg(2+) and Al(3+).     

Halide anion recognition in water by a hexaprotonated octaaza-cryptand: a molecular dynamics investigation

Based on molecular dynamics simulations, we describe the F- versus Cl- complexation by an hexaprotonated cryptand L6+ in aqueous solution, in order to elucidate their structures, solvation properties and the status of external halide counterions. In water, F- and Cl- simulated inclusive complexes adopt a structure somewhat different from the solid state structure of the F- complex: The anion binding involves two diammonium bridges only, and the accompanying counterions are dissociated from the +5 charged complex. A remarkable result is obtained for the dissociated L6+,3F- ,3Cl- system, where spontaneous complexation of F- (the anion which forms the most stable complex with L6+) takes place during the dynamics. The resulting complex is of facial type; this suggests that the equilibrium involves multiple binding modes and structures in aqueous solution. The question of F-/Cl- binding selectivity is investigated by free energy perturbations simulations which nicely reproduce the spectacular preference for F- over Cl-. Two different methodologies used for the treatment of electrostatics (standard versus Ewald calculations) yield similar conclusions.     

Free energy of solvation of simple ions: molecular-dynamics study of solvation of Cl- and Na+ in the ice/water interface

Molecular-dynamics simulations of Cl(-) and Na(+) ions are performed to calculate ionic solvation free energies in both bulk simple point-charge/extended water and ice 1 h at several different temperatures, and at the basal ice 1 h/water interface. For the interface we calculate the free energy of "transfer" of the ions across the ice/water interface. For the ions in bulk water in the NPT ensemble at 298 K and 1 atm, results are found to be in good agreement with experiments, and with other simulation results. Simulations performed in the NVT ensemble are shown to give equivalent solvation free energies, and this ensemble is used for the interfacial simulations. Solvation free energies of Cl(-) and Na(+) ions in ice at 150 K are found to be approximately 30 and approximately 20 kcal mol(-1), respectively, less favorable than for water at room temperature. Near the melting point of the model the solvation of the ions in water is the same (within statistical error) as that measured at room temperature, and in the ice is equivalent and approximately 10 kcal mol(-1) less favorable than the liquid. The free energy of transfer for each ion across ice/water interface is calculated and is in good agreement with the bulk observations for the Cl(-) ion. However, for the model of Na(+) the long-range electrostatic contribution to the free energy was more negative in the ice than the liquid, in contrast with the results observed in the bulk calculations.     

Molecular dynamics study of the coordination sphere of trivalent lanthanum in a highly concentrated LiCl aqueous solution: a combined classical and ab initio approach

The first coordination sphere of trivalent lanthanum in a highly concentrated (14 M) lithium chloride solution is studied with a combination of classical molecular dynamics and density functional theory based first principle molecular dynamics. This method enables us to obtain a solvation shell of La3+ containing 2 chloride ions and 6 water molecules. After refinement using first principle molecular dynamics, the resulting cation-water and cation-anion distances are in very good agreement with experiment. The 2Cl- and the 6 water molecules arrange in a square antiprism around La3+. Exchange of water molecules was also observed in the first-principle simulation, with an intermediate structure comprising 7 water molecules stable for 2.5 ps. Finally, evaluation of dipole moments using maximally localized Wannier functions shows a substantial polarization of the choride anions and the water molecules in the first solvation shell of trivalent lanthanum.     

Indolocarbazole-based foldamers capable of binding halides in water

A series of indolocarbazole-based foldamers have been prepared which can fold into a helical array with an internal cavity encircled by multiple indole NHs, thus allowing for binding anions by hydrogen bonds. The helical folding has been confirmed by computer modeling, 1H NMR spectroscopy, 2D ROESY experiments, and binding studies. A water-soluble derivative binds small, hydrophilic anions in the order Cl- > F- > Br- but negligibly with large, diffuse anions such as I- and ClO4-. Interestingly, the relative binding affinities of the fluoride and chloride ions in water are opposite to each other in an organic medium 4:1 (v/v) DMSO/MeOH, possibly due to the difference in the competing solvation energy.     

Hydration of sodium, potassium, and chloride ions in solution and the concept of structure maker/breaker

Neutron diffraction data with hydrogen isotope substitution on aqueous solutions of NaCl and KCl at concentrations ranging from high dilution to near-saturation are analyzed using the Empirical Potential Structure Refinement technique. Information on both the ion hydration shells and the microscopic structure of the solvent is extracted. Apart from obvious effects due to the different radii of the three ions investigated, it is found that water molecules in the hydration shell of K+ are orientationally more disordered than those hydrating a Na+ ion and are inclined to orient their dipole moments tangentially to the hydration sphere. Cl- ions form instead hydrogen-bonded bridges with water molecules and are readily accommodated into the H-bond network of water. The results are used to show that concepts such as structure maker/breaker, largely based on thermodynamic data, are not helpful in understanding how these ions interact with water at the molecular level.     

Hofmeister salt effects on surface tension arise from partitioning of anions and cations between bulk water and the air-water interface

We apply a recently developed surface-bulk partitioning model to interpret the effects of individual Hofmeister cations and anions on the surface tension of water. The most surface-excluded salt (Na2SO4) provides a minimum estimate for the number of water molecules per unit area of the surface region of 0.2 H2O A-2. This corresponds to a lower bound thickness of the surface region of approximately 6 A, which we assume is a property of this region and not of the salt investigated. At salt concentrations < or = 1 m, single-ion partition coefficients Kp,i, defined relative to Kp,Na+ = Kp,SO42- = 0, are found to be independent of bulk salt concentration and additive for different salt ions. Semiquantitative agreement with surface-sensitive spectroscopy data and molecular dynamics simulations is attained. In most cases, the rank orders of Kp,i for both anions and cations follow the conventional Hofmeister series, qualitative rankings of ions based on their effects on protein processes (folding, precipitation, assembly). Most anions that favor processes that expose protein surface to water (e.g., SCN-), and hence must interact favorably with (i.e., accumulate at) protein surface, are also accumulated at the air-water interface (Kp >1, e.g., Kp,SCN- =1.6). Most anions that favor processes that remove protein surface from water (e.g., F-), and hence are excluded from protein surface, are also excluded from the air-water interface (Kp,F- = 0.5). The guanidinium cation, a strong protein denaturant and therefore accumulated at the protein surface exposed in unfolding, is somewhat excluded from the air-water surface (Kp,GuH+ = 0.7), but is much less excluded than alkali metal cations (e.g., Kp,Na+ identical with 0, Kp,K+ = 0.1). Hence, cation Kp values for the air-water surface appear shifted (toward exclusion) as compared with values inferred for interactions of these cations with protein surface.     

Solvation of cations in the LiNO3-Ca(NO3)2-H2O system at 25°C. A molecular dynamics simulation study

A series of solutions of the ternary LiNO₃-Ca(NO₃)₂-H₂O system is simulated under standard conditions by means of classical molecular dynamics (MD). Several variants of the potentials of interparticle interactions predicting different structures and compositions of the lithium cation hydration shell are used in the calculations. The coordination numbers of ions, the mean residence times of water molecules and nitrate anions in the solvation shells of cations, and the self-diffusion coefficients of solution components are estimated for all of the investigated systems. The structural features of the solvation shells of cations in multicomponent aqueous solutions are described on the basis of the obtained data.     

Simulation of Ca2+ and Mg2+ solvation using polarizable atomic multipole potential

The alkaline earth metals calcium and magnesium are critically involved in many biomolecular processes. To understand the hydration thermodynamics of these ions, we have performed molecular dynamics simulations using a polarizable potential. Particle-mesh Ewald for point multipoles has been applied to the calculation of electrostatic interactions. The parameters in this model have been determined from an ab initio quantum mechanical calculation of dimer interactions between ions and water. Two methods for ion solvation free energy calculation, free energy perturbation, and the Bennett acceptance ratio have been compared. Both predict results consistent with other theoretical estimations while the Bennett approach leads to a much smaller statistical error. Based on the Born theory and the ion-oxygen radial distribution functions, we estimate the effective size of the ions in solution, concluding that K(+) > Na(+) congruent with Ca(2+) > Mg(2+). There appears to be much stronger perturbation in water structure, dynamics, and dipole moment around the divalent cations than the monovalent K(+) and Na(+). The average water coordination numbers for Ca(2+) and Mg(2+) are 7.3 and 6, respectively. The lifetime of water molecules in the first solvation shell of Mg(2+) is on the order of hundreds of picoseconds, in contrast to only few picoseconds for Ca(2+), K(+), or Na(+).     

Alkali cation extraction by calix[4]crown-6 to room-temperature ionic liquids. The effect of solvent anion and humidity investigated by molecular dynamics simulations

We report a molecular dynamics study on the solvation of M+ (Na+ to Cs+) alkali cations and of their LM+ complexes with a calix[4]arene host (L = 1,3-dimethoxy-calix[4]arene-crown-6 in the 1,3-alternate conformation) in the [BMI][PF6] and [BMI][Tf2N] room-temperature ionic liquids "ILs" based on the BMI+ (1-butyl-3-methylimidazolium) cation. The comparison of the two liquids and the dry versus humid form of the former one (with a 1:1 ratio of H2O and BMI+PF6- species) reveals the importance of humidity: in [BMI][PF6]-dry as in the [BMI][Tf2N] liquid, the first solvation shell of the "naked" M+ ions is composed of solvent anions only (four PF6- anions, and from four to five Tf2N- anions, respectively, quasi-neutralized by a surrounding cage of BMI+ cations), while in the [BMI][PF6]-humid IL, it comprises from one to three solvent anions and about four H2O molecules. In the LM+ complexes, the cation is shielded from solvent, but still somewhat interacts with a solvent anion in the dry ILs and with water in the humid IL. We also report tests on M+ interactions with solvent anions PF6- and Tf2N- in the gas phase, showing that the AMBER results are in satisfactory agreement with QM results obtained at different levels of theory. The question of ion recognition by L is then examined by free energy perturbation studies in the three liquids, predicting a high Cs+/Na+ selectivity upon liquid extraction from an aqueous phase, in agreement with experimental results on a parent calixarene host. A similar Cs+/Na+ selectivity is predicted upon complexation in a homogeneous IL phase, mainly due to the desolvation energy of the free cations. Thus, despite their polar character, ionic liquids qualitatively behave as classical weakly polar organic liquids (e.g., choroform) as far as liquid-liquid extraction is concerned but more like polar liquids (water, alcohols) as far as complexation in a single phase is concerned.     

Theoretical study of the solvation of fluorine and chlorine anions by water

The solvation of fluoride and chloride anions (F(-) and Cl(-), respectively) by water has been studied using effective fragment potentials (EFPs) for the water molecules and ab initio quantum mechanics for the anions. In particular, the number of water molecules required to fully surround each anion has been investigated. Monte Carlo calculations have been used in an attempt to find the solvated system X(-)(H(2)O)(n) (X = F, Cl) with the lowest energy for each value of n. It is predicted that 18 water molecules are required to form a complete solvation shell around a Cl(-) anion, where "complete solvation" is interpreted as an ion that is completely surrounded by solvent molecules. Although fewer water molecules may fully solvate the Cl(-) anion, such structures are higher in energy than partially solvated molecules, up to n > or = 18. Calculations on the F(-) anion suggest that 15 water molecules are required for a complete solvation shell. The EFP predictions are in good agreement with the relative energies predicted by ab initio energy calculations at the EFP geometries.     

Mechanism of unassisted ion transport across membrane bilayers

To establish how charged species move from water to the nonpolar membrane interior and to determine the energetic and structural effects accompanying this process, we performed molecular dynamics simulations of the transport of Na+ and Cl- across a lipid bilayer located between two water lamellae. The total length of molecular dynamics trajectories generated for each ion was 10 ns. Our simulations demonstrate that permeation of ions into the membrane is accompanied by the formation of deep, asymmetric thinning defects in the bilayer, whereby polar lipid head groups and water penetrate the nonpolar membrane interior. Once the ion crosses the midplane of the bilayer the deformation "switches sides"; the initial defect slowly relaxes, and a defect forms in the outgoing side of the bilayer. As a result, the ion remains well solvated during the process; the total number of oxygen atoms from water and lipid head groups in the first solvation shell remains constant. A similar membrane deformation is formed when the ion is instantaneously inserted into the interior of the bilayer. The formation of defects considerably lowers the free energy barrier to transfer of the ion across the bilayer and, consequently, increases the permeabilities of the membrane to ions, compared to the rigid, planar structure, by approximately 14 orders of magnitude. Our results have implications for drug delivery using liposomes and peptide insertion into membranes.     

A case against large cluster-induced nonequilibrium solvent effects in supercritical water

Using a partially compressible continuum solvation model, we have shown that solvent compression in just the first two solvation shells (or thereabouts) is all that is required to gain the bulk of the compression-induced enhancement to the solvation energy of ions in supercritical water. This result is found to hold even when the direct, equilibrium solvent-solute cluster involves well over a hundred solvent molecules. We argue that, for charge variation reactions in supercritical water, the observed short-range behavior of the compression-induced solvation free energy precludes the existence of any anomalously large nonequilibrium solvent effects which might be expected on the basis of the very large size of the equilibrium clusters. Received: 8 January 1997 / Accepted: 17 January 1997     

Formation and interaction of hydrated alkali metal ions at the graphite-water interface

Ion hydration at a solid surface ubiquitously exists in nature and plays important roles in many natural processes and technological applications. Aiming at obtaining a microscopic insight into the formation of such systems and interactions therein, we have investigated the hydration of alkali metal ions at a prototype surface-graphite (0001), using first-principles molecular dynamics simulations. At low water coverage, the alkali metal ions form two-dimensional hydration shells accommodating at most four (Li, Na) and three (K, Rb, Cs) waters in the first shell. These two-dimensional shells generally evolve into three-dimensional structures at higher water coverage, due to the competition between hydration and ion-surface interactions. Exceptionally K was found to reside at the graphite-water interface for water coverages up to bulk water limit, where it forms an "umbrellalike" surface hydration shell with an average water-ion-surface angle of 115 degrees . Interactions between the hydrated K and Na ions at the interface have also been studied. Water molecules seem to mediate an effective ion-ion interaction, which favors the aggregation of Na ions but prevents nucleation of K. These results agree with experimental observations in electron energy loss spectroscopy, desorption spectroscopy, and work function measurement. In addition, the sensitive dependence of charge transfer on dynamical structure evolution during the hydration process, implies the necessity to describe surface ion hydration from electronic structure calculations.     

Solvation of uranyl(II), europium(III) and europium(II) cations in "basic" room-temperature ionic liquids: a theoretical study

We report a molecular dynamics study of the solvation of UO2(2+), Eu3+ and Eu2+ ions in two "basic" (Lewis acidity) room-temperature ionic liquids (IL) composed of the 1-ethyl-3-methylimidazolium cation (EMI+) and a mixture of AlCl4- and Cl- anions, in which the Cl-/AlCl4- ratio is about 1 and 3, respectively. The study reveals the importance of the [UO2Cl4]2- species, which spontaneously form during most simulations, and that the first solvation shell of europium is filled with Cl- and AlCl4- ions embedded in a cationic EMI+ shell. The stability of the [UO2Cl4]2- and [Eu(III)Cl6]3- complexes is supported by quantum mechanical calculations, according to which the uranyl and europium cations intrinsically prefer Cl- to the AlCl4- ion. In the gas phase, however, [Eu(III)Cl6]3- and [Eu(II)Cl6]4- complexes are predicted to be metastable and to lose two to three Cl- ions. This contrasts with the results of simulations of complexes in ILs, in which the "solvation" of the europium complexes increases with the number of coordinated chlorides, leading to an equilibrium between different chloro species. The behavior of the hydrated [Eu(OH2)8]3+ complex is considered in the basic liquids; the complex exchanges H2O molecules with Cl- ions to form mixed [EuCl3(OH2)4] and [EuCl4(OH2)3]- complexes. The results of the simulations allow us to better understand the microscopic nature and solvation of lanthanide and actinide complexes in "basic" ionic liquids.     

Ab initio molecular dynamics study of carbon dioxide and bicarbonate hydration and the nucleophilic attack of hydroxide on CO2

We apply ab initio molecular dynamics (AIMD) to study the hydration structures of the carbon dioxide molecule and the bicarbonate and carbonate anions in liquid water. We also compute the free energy change associated with the nucleophilic attack of the hydroxide ion on carbon dioxide. CO2 behaves like a hydrophobic species and exhibits weak interactions with water molecules. The bicarbonate and carbonate ions are strongly hydrated and coordinate to an average of 6.9 and 8.7 water molecules, respectively. The energetics for the reaction in the gas phase are investigated using density functional theory and second-order M?ller-Plesset perturbation theory (MP2) in conjunction with high-quality basis sets. Using umbrella sampling techniques, we compute the standard state, aqueous phase free energy difference associated with the reaction CO2+OH--->HCO3- after correcting AIMD energies with MP2 results. Our predictions are in good agreement with experiments. The hydration structures along the reaction coordinate, which give rise to a predicted 9.7 kcal/mol standard state free energy barrier, are further analyzed.