On the density scaling of pVT data and transport properties for molecular and ionic liquids

In this work, a general equation of state (EOS) recently derived by Grzybowski et al. [Phys. Rev. E 83, 041505 (2011)] is applied to 51 molecular and ionic liquids in order to perform density scaling of pVT data employing the scaling exponent γ(EOS). It is found that the scaling is excellent in most cases examined. γ(EOS) values range from 6.1 for ammonia to 13.3 for the ionic liquid [C(4)C(1)im][BF(4)]. These γ(EOS) values are compared with results recently reported by us [E. R. López, A. S. Pensado, M. J. P. Comu?as, A. A. H. Pádua, J. Fernández, and K. R. Harris, J. Chem. Phys. 134, 144507 (2011)] for the scaling exponent γ obtained for several different transport properties, namely, the viscosity, self-diffusion coefficient, and electrical conductivity. For the majority of the compounds examined, γ(EOS) > γ, but for hexane, heptane, octane, cyclopentane, cyclohexane, CCl(4), dimethyl carbonate, m-xylene, and decalin, γ(EOS) < γ. In addition, we find that the γ(EOS) values are very much higher than those of γ for alcohols, pentaerythritol esters, and ionic liquids. For viscosities and the self-diffusion coefficient-temperature ratio, we have tested the relation linking EOS and dynamic scaling parameters, proposed by Paluch et al. [J. Phys. Chem. Lett. 1, 987-992 (2010)] and Grzybowski et al. [J. Chem. Phys. 133, 161101 (2010); Phys. Rev. E 82, 013501 (2010)], that is, γ = (γ(EOS)/φ) + γ(G), where φ is the stretching parameter of the modified Avramov relation for the density scaling of a transport property, and γ(G) is the Gru?neisen constant. This relationship is based on data for structural relaxation times near the glass transition temperature for seven molecular liquids, including glass formers, and a single ionic liquid. For all the compounds examined in our much larger database the ratio (γ(EOS)/φ) is actually higher than γ, with the only exceptions of propylene carbonate and 1-methylnaphthalene. Therefore, it seems the relation proposed by Paluch et al. applies only in certain cases, and is really not generally applicable to liquid transport properties such as viscosities, self-diffusion coefficients or electrical conductivities when examined over broad ranges of temperature and pressure.     

Effect of anions on static orientational correlations, hydrogen bonds, and dynamics in ionic liquids: a simulational study

Three different ionic liquids are investigated via atomistic molecular dynamics simulations using the force field of Lopes and PAdua (J. Phys. Chem. B 2006, 110, 19586). In particular, the 1-ethyl-3-methylimidazolium cation EMIM+ is studied in the presence of three different anions, namely, chloride Cl-, tetrafluoroborate BF(4)(-), and bis((trifluoromethyl)sulfonyly)imide TF2N-. In the focus of the present study are the static distributions of anions and cations around a cation as a function of anion size. It is found that the preferred positions of the anions change from being close to the imidazolium hydrogens to being above and below the imidazolium rings. Lifetimes of hydrogen bonds are calculated and found to be of the same order of magnitude as those of pure liquid water and of some small primary alcohols. Three kinds of short-range cation-cation orderings are studied, among which the offset stacking dominates in all of the investigated ionic liquids. The offset stacking becomes weaker from [EMIM][Cl] to [EMIM][BF4] to [EMIM][TF2N]. Further investigation of the dynamical behavior reveals that cations in [EMIM][TF2N] have a slower tumbling motion compared with those in [EMIM][Cl] and [EMIM][BF4] and that pure diffusive behavior can be observed after 1.5 ns for all three systems at temperatures 90 K above the corresponding melting temperatures.     

Nucleophilicity in ionic liquids. 3. Anion effects on halide nucleophilicity in a series of 1-butyl-3-methylimidazolium ionic liquids

We have continued the study of halide nucleophilicity in ionic liquids, concentrating on the effect of changing the anion ([BF(4)](-), [PF(6)](-), [SbF(6)](-), [OTf](-), and [N(Tf)(2)](-)) when the cation is [bmim](+) (where bmim = 1-butyl-3-methylimidazolium). It was found that the nucleophilicities of all the halides were lower in all of the ionic liquids than in dichloromethane. Changing the anion affected the order of halide nucleophilicity, e.g., in [bmim][BF(4)] the order of nucleophilicity was Cl(-)>Br(-)>I(-) while in [bmim][N(Tf)(2)] the order was Cl(-)相似文献   

Simulation of the surface structure of butylmethylimidazolium ionic liquids

Molecular dynamics simulations of the liquid/vacuum surfaces of the room temperature ionic liquids [bmim][PF(6)], [bmim][BF(4)] and [bmim][Cl] have been carried out at various temperatures. The surfaces are structured with a top monolayer containing oriented cations and anions. The butyl side chains tend to face the vacuum and the methyl side chains the liquid. However, as the butyl chains are not densely packed, both anions and rings are visible from the vacuum phase. The effects of temperature and the anion on the degree of cation orientation is small, but the potential drop from the vacuum to the interior of the liquid is greater for liquids with smaller anions. We compare the simulation results with a range of experimental observations and suggest that neutron reflection from samples with protiated butyl groups would be a sensitive probe of the structure.     

Effect of the structures of ionic liquids and alkylbenzene-derived amphiphiles on the inhibition of asphaltene precipitation from CO2-injected reservoir oils

The inhibition of asphaltene precipitation from high-pressure, CO(2)-injected reservoir oils by ionic and nonionic amphiphiles, the ionic liquids based on p-alkylpyridinium ([C(n)()py](+)) and N-butylisoquinolinium ([C(4)iql](+)) cations, and the alkylbenzene-derived amphiphiles p-alkylphenol (C(n)()phol), p-alkylbenzenesulfonic acid (C(n)()bsa), and sodium p-alkylbenzenesulfonate (C(n)()bsNa) was investigated for the first time. The influences of the structures of these compounds and the effect of the combination of their cations and anions were studied. The results show that the inhibition abilities of the alkylbenzene-derived amphiphiles first increase when n = 2-8 and then remain almost constant when n >/=8 and that the effectiveness follows the order C(n)()phol < C(n)()bsa approximately C(n)()bsNa. The inverse trend is observed for the ionic liquids [C(n)()py][Cl]; that is, their inhibition abilities decrease as n increases from 4 to 8 to 12. [C(4)iql][Cl] is more effective than [C(4)py][Cl], but [C(n)()py][BF(4)] and [C(n)()py][PF(6)] have almost no effect on the stabilization of asphaltenes. It was found that the effectiveness of an alkylbenzene-derived amphiphile on the inhibition of asphaltene precipitation from reservoir oils relies on its ability to form a stable steric-stabilization layer around asphaltenes, which is controlled by the polarity of its headgroup and the length of its alkyl tail. The novel mechanism of inhibiting asphaltene precipitation using the ionic liquids [C(n)()py] ([Cl], [BF(4)], and [PF(6)]) and [C(4)iql][Cl] was proposed. The mechanism states that the ionic liquids can effectively prevent asphaltene precipitation from the reservoir oils by breaking the asphaltene associations, which are due to the local nonneutrality of the charge densities of the cation and the anion. The ionic liquids that are based on an anion with high charge density, in connection with cations with sufficiently low charge densities, can effectively inhibit asphaltene precipitation from the reservoir oils. This mechanism is also important for studying the thermodynamic properties and phase behavior of the ionic liquids.     

分子模拟研究气体在室温离子液体中的溶解度

在作者先前建立的分子力场基础上, 采用Widom粒子插入法预测了CO₂、N₂、O₂、Ar及CH₄等5种气体在多种咪唑类离子液体中的溶解度, 包括2种侧链长度的阳离子和3种阴离子. 首先考察了256个离子对组成的体系中溶质分子插入次数对计算结果的影响, 在此基础上计算了不同温度下气体在1-丁基-3-甲基咪唑的四氟化硼盐([bmim][BF₄])和六氟化磷盐([bmim][PF₆])中的溶解度. 计算结果正确反映了CO₂气体溶解度的变化趋势, 在[bmim][BF₄]中溶解度的模拟结果与实验值符合很好, 且明显优于Pádua等的模拟结果；在[bmim][PF₆]中的溶解度较实验值偏高, 精度与文献模拟结果相当；并预测了较高温度下CO₂气体在[bmim][BF₄]和[bmim][PF₆]中的溶解度. 计算结果也正确地反映了5种气体在[bmim][PF₆]中溶解度实验值的相对大小. 另外考察了常温下几种气体在不同室温离子液体中的溶解度, 模拟结果表明气体在含有较长碳链和双-三氟甲基磺酰胺阴离子(Tf₂N^－)的离子液体中溶解度较大.     

Computation of methodology-independent single-ion solvation properties from molecular simulations. IV. Optimized Lennard-Jones interaction parameter sets for the alkali and halide ions in water

The raw single-ion solvation free energies computed from atomistic (explicit-solvent) simulations are extremely sensitive to the boundary conditions and treatment of electrostatic interactions used during these simulations. However, as shown recently [M. A. Kastenholz and P. H. Hu?nenberger, J. Chem. Phys. 124, 224501 (2006); M. M. Reif and P. H. Hu?nenberger, J. Chem. Phys. 134, 144103 (2010)], the application of appropriate correction terms permits to obtain methodology-independent results. The corrected values are then exclusively characteristic of the underlying molecular model including in particular the ion-solvent van der Waals interaction parameters, determining the effective ion size and the magnitude of its dispersion interactions. In the present study, the comparison of calculated (corrected) hydration free energies with experimental data (along with the consideration of ionic polarizabilities) is used to calibrate new sets of ion-solvent van der Waals (Lennard-Jones) interaction parameters for the alkali (Li(+), Na(+), K(+), Rb(+), Cs(+)) and halide (F(-), Cl(-), Br(-), I(-)) ions along with either the SPC or the SPC/E water models. The experimental dataset is defined by conventional single-ion hydration free energies [Tissandier et al., J. Phys. Chem. A 102, 7787 (1998); Fawcett, J. Phys. Chem. B 103, 11181] along with three plausible choices for the (experimentally elusive) value of the absolute (intrinsic) hydration free energy of the proton, namely, ΔG(hyd)(?)[H(+)] = -1100, -1075 or -1050 kJ mol(-1), resulting in three sets L, M, and H for the SPC water model and three sets L(E), M(E), and H(E) for the SPC/E water model (alternative sets can easily be interpolated to intermediate ΔG(hyd)(?)[H(+)] values). The residual sensitivity of the calculated (corrected) hydration free energies on the volume-pressure boundary conditions and on the effective ionic radius entering into the calculation of the correction terms is also evaluated and found to be very limited. Ultimately, it is expected that comparison with other experimental ionic properties (e.g., derivative single-ion solvation properties, as well as data concerning ionic crystals, melts, solutions at finite concentrations, or nonaqueous solutions) will permit to validate one specific set and thus, the associated ΔG(hyd)(?)[H(+)] value (atomistic consistency assumption). Preliminary results (first-peak positions in the ion-water radial distribution functions, partial molar volumes of ionic salts in water, and structural properties of ionic crystals) support a value of ΔG(hyd)(?)[H(+)] close to -1100 kJ·mol(-1).     

Communication: Thermodynamic scaling of the Debye process in primary alcohols

The molecular dynamics of hydrogen-bonded liquids usually does not satisfy the thermodynamic scaling. However, very recently, two opposite conclusions about validity of thermodynamical scaling in monohydroxy alcohol, 2-ethyl-1-hexanol, were presented by Reiser et al. [J. Chem. Phys. 132, 181101 (2010)] and Fragiadakis et al. [J. Chem. Phys. 132, 144505 (2010)]. In this communication we present new experimental results that can explain this ostensible contradiction.     

Molecular dynamics simulation studies of the influence of imidazolium structure on the properties of imidazolium/azide ionic liquids

Atomistic molecular dynamics simulations were performed on 1-butyl-3-methyl-imidazolium azide [bmim][N(3)], 1-butyl-2,3-dimethylimidazolium azide [bmmim][N(3)], and 1-butynyl-3-methyl-imidazolium azide [bumim][N(3)] ionic liquids. The many-body polarizable APPLE&P force field was augmented with parameters for the azide anion and the bumim cation. Good agreement between the experimentally determined and simulated crystal structure of [bumim][N(3)] as well as the liquid-state density and ionic conductivity of [bmmim][N(3)] were found. Methylation of bmim (yielding bmmim) resulted in dramatic changes in ion structuring in the liquid and slowing of ion motion. Conversely, replacing the butyl group of bmim with the smaller 2-butynyl group resulted in an increase of ion dynamics.     

The free energy of the metastable supersaturated vapor via restricted ensemble simulations

Pressure, excess chemical potential, and excess free energy, with respect to ideal gas data at different densities of the supersaturated Lennard-Jones particle vapor at the reduced temperature 0.7 are obtained by the restricted canonical ensemble Monte Carlo simulation method [D. S. Corti and P. Debenedetti, Chem. Eng. Sci. 49, 2717 (1994)]. The excess free energy values depend upon the constraints imposed on the system with local minima exhibited for densities below the spinodal density and monotonic variation for densities larger than the spinodal density. The results are compared with a molecular dynamics simulation study [A. Linharton et al., J. Chem. Phys. 122, 144506 (2005)] on the same system. The current study verifies the conclusion drawn by the simulation work that clustering of Lennard-Jones atoms exists even in the vicinity of spinodal. Our method gives an alternative to molecular dynamic simulations for the determination of equilibrium properties of a metastable fluid, especially close to the spinodal, and does not require a very large system to carry out the simulation.     

A new force field model for the simulation of transport properties of imidazolium-based ionic liquids

A new, non-polarizable force field model (FFM) for imidazolium-based, room-temperature ionic liquids (RTILs), 1-ethyl-3-methyl-imidazolium tetrafluoroborate and 1-butyl-3-methyl-imidazolium tetrafluoroborate, has been developed. Modifying the FFM originally designed by Liu et al. (J. Phys. Chem. B, 2004, 108, 12978-12989), the electrostatic charges on interacting sites are refined according to partial charges calculated by explicit-ion density functional theory. The refined FFM reproduces experimental heats of vaporization, diffusion coefficients, ionic conductivities, and shear viscosities of RTILs, which is a significant improvement over the original model (Zh. Liu, Sh. Huang and W. Wang, J. Phys. Chem. B, 2004, 108, 12978-12989). The advantages of the proposed procedure include clarity, simplicity, and flexibility. Expanding the functionality of our FFM conveniently only requires modification of the electrostatic charges. Our FFM can be extended to other classes of RTILs as well as condensed matter systems in which the ionic interaction requires an account of polarization effects.     

Diffusion in linear porous media with periodic entropy barriers: A tube formed by contacting spheres

The problem of transport in quasi-one-dimensional periodic structures has been studied recently by several groups [D. Reguera et al., Phys. Rev. Lett.96, 130603 (2006); P. S. Burada et al., Phys. Rev. E75, 051111 (2007); B. Q. Ai and L. G. Liu, ibid.74, 051114 (2006); B. Q. Ai et al., ibid.75, 061126 (2007); B. Q. Ai and L. G. Liu, J. Chem. Phys.126, 204706 (2007); 128, 024706 (2008); E. Yariv and K. D. Dorfman, Phys. Fluids19, 037101 (2007); N. Laachi et al., Europhys. Lett.80, 50009 (2007); A. M. Berezhkovskii et al., J. Chem. Phys.118, 7146 (2003); 119, 6991 (2003)]. Using the concept of "entropy barrier" [R. Zwanzig, J. Phys. Chem.96, 3926 (1992)] one can classify such structures based on the height of the entropy barrier. Structures with high barriers are formed by chambers, which are weakly connected with each other because they are connected by small apertures. To escape from such a chamber a diffusing particle has to climb a high entropy barrier to find an exit that takes a lot of time [I. V. Grigoriev et al., J. Chem. Phys.116, 9574 (2002)]. As a consequence, the particle intrachamber lifetime tau(esc) is much larger than its intrachamber equilibration time, tau(rel), tau(esc)>tau(rel). When the aperture is not small enough, the intrachamber escape and relaxation times are of the same order and the hierarchy fails. This is the case of low entropy barriers. Transport in this case is analyzed in the works of Schmid and co-workers, Liu and co-workers, and Dorfman and co-workers, while the work of Berezhkovskii et al. is devoted to diffusion in the case of high entropy barriers.     

The glass transition and the distribution of voids in room-temperature ionic liquids: a molecular dynamics study

The glass transition in prototypical room temperature ionic liquids has been investigated by molecular dynamics simulations based on an Amber-like empirical force field. Samples of [C(4)mim][PF(6)], [C(4)mim][Tf(2)N], and [C(3)mim][Tf(2)N] have been quenched from the liquid phase at T = 500 to a glassy state at T ～ 0 K in discontinuous steps of 20 K every 1.2 ns. The glass temperature estimated by simulation (T(g) = 209 K for [C(4)mim][PF(6)], T(g) = 204 K for [C(4)mim][Tf(2)N], and T(g) = 196 K for [C(3)mim][Tf(2)N]) agrees semi-quantitatively with the experimental values (T(g) = 193÷196 K for [C(4)mim][PF(6)], T(g) = 186÷189 K for [C(4)mim][Tf(2)N], and T(g) = 183 K for [C(3)mim][Tf(2)N]). A model electron density is introduced to identify voids in the system. The temperature dependence of the size distribution of voids provided by simulation reproduce well the experimental results of positron annihilation lifetime spectroscopy reported in G. Dlubek, Y. Yu, R. Krause-Rehberg, W. Beichel, S. Bulut, N. Pogodina, I. Krossing, and Ch. Friedrich, J. Chem. Phys. 133, 124502 (2010), with only one free parameter needed to fit the experimental data.     

Shear viscosity of the ionic liquid 1-n-butyl 3-methylimidazolium hexafluorophosphate [bmim][pf6] computed by reverse nonequilibrium molecular dynamics

Reverse nonequilibrium molecular dynamics and equilibrium molecular dynamics simulations were carried out to compute the shear viscosity of the pure ionic liquid system [bmim][PF 6] at 300 K. The two methods yielded consistent results which were also compared to experiments. The results showed that the reverse nonequilibrium molecular dynamics (RNEMD) methodology can successfully be applied to computation of highly viscous ionic liquids. Moreover, this study provides a validation of the atomistic force-field developed by Bhargava and Balasubramanian ( J. Chem. Phys. 2007, 127, 114510 ) for dynamic properties.     

Correlation between polarity parameters and dielectric properties of [Na][TOTO]--a sodium ionic liquid

The ionic liquid (IL) [Na][TOTO], with sodium as the cation and an oligoethercarboxylate as the anion, shows properties that differ significantly from conventional ionic liquids, like imidazolium salts. Its polarity, determined in the temperature range of (293-333) K from measurements of Reichardt's E value and the Kamlet-Taft parameters is extraordinarily low and matches the dielectric constant extrapolated from frequency dependent complex permittivity measurements. The dielectric spectra also reveal split dynamics with the dominating slow mode probably associated with the reorientation of -COO(-)...Na(+) ion pairs and the fast mode arising from the flexibility of the oligoethylene moieties of the anion. This assignment, supported by the derived effective dipole moments, argues against our previous hypothesis [O. Zech et al., Chem.-Eur. J., 2009, 15, 1341-1345] that salts, such as [Na][TOTO], have low melting points because of cation complexation in a pseudo crown ether-like fashion and resulting "intra-molecular" charge neutralization. The present data rule out a rigid chelate-like complex as the dominating species. Considering the present findings together with the viscosities of [TOTO](-) salts, a crosslinked structure of the ionic liquid with strong -COO(-)...Na(+) interactions and relatively weaker forces between cations and ether oxygen groups of the anions appears to be plausible.     

Rotational dynamics of water and benzene controlled by anion field in ionic liquids: 1-butyl-3-methylimidazolium chloride and hexafluorophosphate

The rotational correlation time (tau(2R)) is determined for D(2)O (polar) and C(6)D(6) (apolar) in 1-butyl-3-methylimidazolium chloride ([bmim][Cl]) and hexafluorophosphate ([bmim][PF(6)]) by measuring (2)H (D) nuclear magnetic resonance spin-lattice relaxation time (T(1)) in the temperature range from -20 to 110 degrees C. The tau(2R) ratio of water to benzene (tau(WB)) was used as a measure of solute-solvent attraction. tau(WB) is 0.73 and 0.52 in [bmim][Cl] and [bmim][PF(6)], respectively, whereas the molecular volume ratio is as small as 0.11. The slowdown of the water dynamics compared to the benzene dynamics in ionic liquids is interpreted by the Coulombic attractive interaction between the polar water molecule and the anion. As for the anion effect, the rotational dynamics of water solvated by Cl(-) is slower than that solvated by PF(6) (-), whereas the rotational dynamics of benzene is similar in the two ionic liquids. This is interpreted as an indication of the stronger solvation by the anion with a larger surface charge density. The slowdown of the water dynamics via Coulombic solvation is actually significant only at water concentrations lower than approximately 9 mol dm(-3) at room temperature, and it is indistinguishable at temperatures above approximately 100 degrees C. The quadrupolar coupling constants determined for D(2)O and C(6)D(6) in the ionic liquids were smaller by a factor of 2-3 than those in the pure liquid state.     

Molecular dynamics studies of melting and some liquid-state properties of 1-ethyl-3-methylimidazolium hexafluorophosphate [emim][PF6

Molecular dynamics simulations are used to study the liquid-state properties and melting of 1-ethyl-3-methylimidazolium hexafluorosphosphate [emim][PF6] using the force field of Canongia Lopes et al. [J. Phys. Chem. B 108, 2038 (2004)] and geometric constants from crystallographic data. The structures of the solid and liquid states are characterized by carbon-carbon, carbon-phosphorous, and phosphorous-phosphorous radial distribution functions. Spatial correlations among the ions are strong in the liquid state. The cohesive energy density and the temperature dependences of the molar volume and density of the liquid have been computed. The melting point is determined by equilibrating the solid-state supercells in which void defects have been introduced to eliminate the free-energy barrier for the formation of a solid-liquid interface. The computed melting point is 375+/-10 K, which is approximately 10% higher than the experimental value of 333 K.     

Single particle dynamics in ionic liquids of 1-alkyl-3-methylimidazolium cations

Ionic dynamics in room temperature molten salts (ionic liquids) containing 1-alkyl-3-methylimidazolium cations is investigated by molecular-dynamics simulations. Calculations were performed with united atom models, which were used in a previous detailed study of the equilibrium structure of ionic liquids [S. M. Urahata and M. C. C. Ribeiro, J. Chem. Phys. 120, 1855 (2004)]. The models were used in a systematic study of the dependency of several single particle time correlation functions on anion size (F-, Cl-, Br-, and PF6-) and alkyl chain length (1-methyl-, 1-ethyl-, 1-butyl-, and 1-octyl-). Despite of large mass and size of imidazolium cations, they exhibit larger mean-square displacement than anions. A further detailed picture of ionic motions is obtained by using appropriate projections of displacements along the plane or perpendicular to the plane of the imidazolium ring. A clear anisotropy in ionic displacement is revealed, the motion on the ring plane and almost perpendicular to the 1-alkyl chain being the less hindered one. Similar projections were performed on velocity correlation functions, whose spectra were used to relate short time ionic rattling with the corresponding long time diffusive regime. Time correlation functions of cation reorientation and dihedral angles of the alkyl chains are discussed, the latter decaying much faster than the former. A comparative physical picture of time scales for distinct dynamical processes in ionic liquids is provided.     

Revealing long-range density fluctuations in dialkylimidazolium chloride ionic liquids by dynamic light scattering

In this study, the structures and dynamics of ionic liquids of 1-allyl-3-methylimidazolium chloride ([AMIM][Cl]) and 1-n-butyl-3-methylimidazolium chloride ([BMIM][Cl]) were studied by dynamic light scattering with polarized and depolarized geometries in the temperature range from 300 to 400 K. The temperature range covered supercooled and liquid states for [BMIM][Cl] and covered the liquid state for [AMIM][Cl]. The results show that for these ionic liquids at all chosen temperatures only one ultraslow relaxation is observed in the polarized component of dynamic light scattering, however, the ultraslow relaxation is not observed in the depolarized component. The ultraslow relaxation exhibited several typical features of the "cluster" mode generally found in glass-forming liquids and polymer melts, such as diffusive, strongly scattering-vector-dependent, and nearly exponential characters, which thus corresponded to long-range density fluctuations. The physical origin for long-range density fluctuations was the existence of heterogeneities with large characteristic length scales in these ionic liquids. It was further considered that molecules of these ionic liquids not only tended to aggregate to form dynamic clusters but also possibly formed dynamic networks in the supercooled state and the heterogeneities could exist even at temperatures higher than the melting points.     

Protein structure and dynamics in ionic liquids. Insights from molecular dynamics simulation studies

We present in this work the first molecular simulation study of an enzyme, the serine protease cutinase from Fusarium solani pisi, in two ionic liquids (ILs): 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]) and 1-butyl-3-methylimidazolium nitrate ([BMIM][NO(3)]). We tested different water contents in these ILs at room temperature (298 K) and high temperature (343 K), and we observe that the enzyme structure is highly dependent on the amount of water present in the IL media. We show that the enzyme is preferentially stabilized in [BMIM][PF6] at 5-10% (w/w) (weight of water over protein) water content at room temperature. [BMIM][PF6] renders a more nativelike enzyme structure at the same water content of 5-10% (w/w) as previously found for hexane, and the system displays a similar bell-shape-like dependence with the water content in the IL media. [BMIM][PF6] is shown to increase significantly the protein thermostability at high temperatures, especially at low hydration. Our analysis indicates that the enzyme is less stabilized in [BMIM][NO(3)] relative to [BMIM][PF6] at both temperatures, most likely due to the strong affinity of the [NO(3)]- anion toward the protein main chain. These findings are in accordance with the experimental knowledge for these two ionic liquids. We also show that these ILs "strip off" most of the water from the enzyme surface in a degree similar to that found for polar organic solvents such as acetonitrile, and that the remaining waters at the enzyme surface are organized in many small clusters.