Room temperature ionic liquids (IL) have been used in numerous applications in chemistry. Addition of water alters many of their properties making it possible to custom design solvents for specific applications. Along with experiments, computational studies using various approaches have provided key insights into the structure and dynamics of IL systems, as well as aggregate formation and phase behavior of the IL/water mixtures. These systems provide computational challenges since ILs and IL/water mixtures are viscous liquids with intrinsically slow processes and structural organization over surprisingly large length scales, which push the limits of applicability of the available techniques. Recent developments in the studies of IL/water mixtures using computational methodologies are reviewed and the future prospects for the field are briefly discussed.

Kinetics and equilibrium are studied on the hydrothermal decarbonylation and decarboxylation of formic acid, the intermediate of the water-gas-shift (WGS) reaction, in hot water at temperatures of 170-330 degrees C, to understand and control the hydrothermal WGS reaction. (1)H and (13)C NMR spectroscopy is applied to analyze as a function of time the quenched reaction mixtures in both the liquid and gas phases. Only the decarbonylation is catalyzed by HCl, and the reaction is first-order with respect to both [H(+)] and [HCOOH]. Consequently, the reaction without HCl is first and a half (1.5) order due to the unsuppressed ionization of formic acid. The HCl-accelerated decarbonylation path can thus be separated in time from the decarboxylation. The rate and equilibrium constants for the decarbonylation are determined separately by using the Henry constant (gas solubility data) for carbon monoxide in hot water. The rate constant for the decarbonylation is 1.5 x 10(-5), 2.0 x 10(-4), 3.7 x 10(-3), and 6.3 x 10(-2) mol(-1) kg s(-1), respectively, at 170, 200, 240, and 280 degrees C on the liquid branch of the saturation curve. The Arrhenius plot of the decarbonylation is linear and gives the activation energy as 146 +/- 3 kJ mol(-1). The equilibrium constant K(CO) = [CO]/[HCOOH] is 0.15, 0.33, 0.80, and 4.2, respectively, at 170, 200, 240, and 280 degrees C. The van't Hoff plot results in the enthalpy change of DeltaH = 58 +/- 6 kJ mol(-1). The decarboxylation rate is also measured at 240-330 degrees C in both acidic and basic conditions. The rate is weakly dependent on the solution pH and is of the order of 10(-4) mol kg(-1) s(-1) at 330 degrees C. Furthermore, the equilibrium constant K(CO2) = [CO(2)][H(2)]/[HCOOH] is estimated to be 1.0 x10(2) mol kg(-1) at 330 degrees C.

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.

The equilibrium for the reversible decomposition of formic acid into carbon dioxide and hydrogen is studied in the ionic liquid (IL) 1,3-dipropyl-2-methylimidazolium formate. The equilibrium is strongly favored to the formic acid side because of the strong solvation of formic acid in the IL through the strong Coulombic solute-solvent interactions. The comparison of the equilibrium constants in the IL and water has shown that the pressures required to transform hydrogen and carbon dioxide into formic acid can be reduced by a factor of approximately 100 by using the IL instead of water. The hydrogen transformation in such mild conditions can be a chemical basis for the hydrogen storage and transportation using formic acid.

Translational diffusion coefficients of diphenylcyclopropenone (DPCP), diphenylacetylene (DPA), and carbon monoxide (CO) in 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([BMIm][NTf2]) and 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([EMIm][NTf2]) were determined by the transient grating (TG) spectroscopy under pressure from 0.1 to 200 MPa at 298 K and from 298 to 373 K under 0.1 MPa. Diffusion coefficients of these molecules at high temperatures in tributylmethylphosphonium bis(trifluoromethanesulfonyl)imide ([P4441][NTf2]), and tetraoctylphosphonium bis(trifluoromethanesulfonyl)imide ([P8888][NTf2]), and also in the mixtures of [BMIm][NTf2], N-methyl-N-propylpiperidinium bis(trifluoromethanesulfonyl)imide ([Pp13][NTf2]), and trihexyltetradecylphosphonium bis(trifluoromethanesulfonyl)imide ([P66614][NTf2]) with ethanol or chloroform have been determined. Diffusion coefficients except in ILs of phosphonium cations were well scaled by the power law of T/η, i.e., (T/η)(P), where T and η are the absolute temperature and the viscosity, irrespective of the solvent species, pressure and temperature, and the compositions of mixtures. The values of the exponent P were smaller for the smaller size of the molecules. On the other hand, the diffusion coefficients in ILs of phosphonium cations with longer alkyl chains were larger than the values expected from the correlation obtained by other ILs and conventional liquids. The deviation becomes larger with increasing the number of carbon atoms of alkyl-chain of cation, and with decreasing the molecular size of diffusing molecules. The molecular size dependence of the diffusion coefficient was correlated by the ratio of the volume of the solute to that of the solvent as demonstrated by the preceding work (Kaintz et al., J. Phys. Chem. B 2013 , 117 , 11697 ). Diffusion coefficients have been well correlated with the power laws of both T/η and the relative volume of the solute to the solvent.

We have investigated the assembly of a two-dimensional coordination polymer, Nd(2)(C(6)H(2)N(2)O(4))(2)(C(2)O(4))(H(2)O)(2), that has been prepared from the hydrothermal reaction of Nd(NO(3))(3)·6H(2)O and 2,3-pyrazinedicarboxylic acid (H(2)pzdc). In situ oxalate formation as observed in this system has been been investigated using (1)H and (13)C nuclear magnetic resonance spectroscopy, and a pathway for C(2)O(4)(2-) anion formation under hydrothermal conditions has been elucidated. The oxalate ligands found in Nd(2)(C(6)H(2)N(2)O(4))(2)(C(2)O(4))(H(2)O)(2) result from the oxidation of H(2)pzdc, which proceeds through intermediates, such as 2-pyrazinecarboxylic acid (2-pzca), 2-hydroxyacetamide, 3-amino-2-hydroxy-3-oxopropanoic acid, 2-hydroxymalonic acid, 2-oxoacetic acid (glyoxylic acid), and glycolic acid. The species are generated through a ring-opening that occurs via cleavage of the C-N bond of the pyrazine ring, followed by hydrolysis/oxidation of the resulting species.

The rotational dynamics of benzene and water in the ionic liquid (IL) 1-butyl-3-methylimidazolium chloride are studied using molecular dynamics (MD) simulation and NMR T(1) measurements. MD trajectories based on an effective potential are used to calculate the (2)H NMR relaxation time, T(1) via Fourier transform of the relevant rotational time correlation function, C(2R)(t). To compensate for the lack of polarization in the standard fixed-charge modeling of the IL, an effective ionic charge, which is smaller than the elementary charge is employed. The simulation results are in closest agreement with NMR experiments with respect to the temperature and Larmor frequency dependencies of T(1) when an effective charge of ±0.5e is used for the anion and the cation, respectively. The computed C(2R)(t) of both solutes shows a bi-modal nature, comprised of an initial non-diffusive ps relaxation plus a long-time ns tail extending to the diffusive regime. Due to the latter component, the solute dynamics is not under the motional narrowing condition with respect to the prevalent Larmor frequency. It is shown that the diffusive tail of the C(2R)(t) is most important to understand frequency and temperature dependencies of T(1) in ILs. On the other hand, the effect of the initial ps relaxation is an increase of T(1) by a constant factor. This is equivalent to an "effective" reduction of the quadrupolar coupling constant (QCC). Thus, in the NMR T(1) analysis, the rotational time correlation function can be modeled analytically in the form of aexp (-t/τ) (Lipari-Szabo model), where the constant a, the Lipari-Szabo factor, contains the integrated contribution of the short-time relaxation and τ represents the relaxation time of the exponential (diffusive) tail. The Debye model is a special case of the Lipari-Szabo model with a = 1, and turns out to be inappropriate to represent benzene and water dynamics in ILs since a is as small as 0.1. The use of the Debye model would result in an underestimation of the QCC by a factor of 2-3 as a compensation for the neglect of the Lipari-Szabo factor.

The H/D exchange reaction and the rotational dynamics of heavy water (D2O) are studied at 50 degrees C in the ionic liquid, 1-butyl-3-methylimidazolium chloride ([bmim][Cl]), in the [D2O] range of 3-55 M. The initial H/D exchange rates are observed as 1.0 x 10(-7), 4.5 x 10(-6), 1.0 x 10(-5), 4.1 x 10(-5), 1.1 x 10(-4), and 3.7 x 10(-4) s(-1), respectively, at [D2O] of 2.8, 7.1, 8.1, 11, 15, and 25 M. The rate is very slow and less than 10(-5) s(-1) at [D2O] below approximately 7 M. It steeply increases to the order of 10(-4)s(-1) for 7 M < [D2O] < 10 M, and linearly increases with [D2O] in the more water-rich region. The intercept of the linear region at [D2O] = approximately 9 M is interpreted by considering that each chloride anion deactivates 1.6 equiv water molecules due to the strong solvation. Correspondingly, the rotational correlation time of D2O at [D2O] < 7 M is 1 order of magnitude larger than that in water-rich conditions.