Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Polarization Relaxation in an Ionic Liquid Confined between Electrified Walls

The response of a room temperature molten salt to an external electric field when it is confined to a nanoslit is studied by molecular dynamics simulations. The fluid is confined between two parallel and oppositely charged walls, emulating two electrified solid-liquid interfaces. Attention is focused on structural, electrostatic, and dynamical properties, which are compared with those of the nonpolarized fluid. It is found that the relaxation of the electrostatic potential, after switching the electric field off, occurs in two stages. A first, subpicosecond process accounts for 80% of the decay and is followed by a second subdiffusive process with a time constant of 8 ps. Diffusion is not involved in the relaxation, which is mostly driven by small anion translations. The relaxation of the polarization in the confined system is discussed in terms of the spectrum of charge density fluctuations in the bulk.

Structure and dynamics of a confined ionic liquid. topics of relevance to dye-sensitized solar cells

The behavior of a model ionic liquid (IL) confined between two flat parallel walls was studied at various interwall distances using computer simulations. The results focus both on structural and dynamical properties. Mass and charge density along the confinement axis reveal a structure of layers parallel to the walls that leads to an oscillatory profile in the electrostatic potential. Orientational correlation functions indicate that cations at the interface orient tilted with respect to the surface and that any other orientational order is lost thereafter. The diffusion coefficients of the ions exhibit a maximum as a function of the confinement distance, a behavior that results from a combination of the structure of the liquid as a whole and a faster molecular motion in the vicinity of the walls. We discuss the relevance of the present results and elaborate on topics that need further attention regarding the effects of ILs in the functioning of IL-based dye-sensitized solar cells.

Mutual Solubilities of Water and Hydrophobic Ionic Liquids

The ionic nature of ionic liquids (ILs) results in a unique combination of intrinsic properties that produces increasing interest in the research of these fluids as environmentally friendly "neoteric" solvents. One of the main research fields is their exploitation as solvents for liquid-liquid extractions, but although ILs cannot vaporize leading to air pollution, they present non-negligible miscibility with water that may be the cause of some environmental aquatic risks. It is thus important to know the mutual solubilities between ILs and water before their industrial applications. In this work, the mutual solubilities of hydrophobic yet hygroscopic imidazolium-, pyridinium-, pyrrolidinium-, and piperidinium-based ILs in combination with the anions bis(trifluoromethylsulfonyl)imide, hexafluorophosphate, and tricyanomethane with water were measured between 288.15 and 318.15 K. The effect of the ILs structural combinations, as well as the influence of several factors, namely cation side alkyl chain length, the number of cation substitutions, the cation family, and the anion identity in these mutual solubilities are analyzed and discussed. The hydrophobicity of the anions increases in the order [C(CN)3] <[PF6] <[Tf2N] while the hydrophobicity of the cations increases from [Cnmim] <[Cnmpy] [Cnmpyr] <[Cnmpip] and with the alkyl chain length increase. From experimental measurements of the temperature dependence of ionic liquid solubilities in water, the thermodynamic molar functions of solution, such as Gibbs energy, enthalpy, and entropy at infinite dilution were determined, showing that the solubility of these ILs in water is entropically driven and that the anion solvation at the IL-rich phase controls their solubilities in water. The COSMO-RS, a predictive method based on unimolecular quantum chemistry calculations, was also evaluated for the description of the water-IL binary systems studied, where it showed to be capable of providing an acceptable qualitative agreement with the experimental data.

Mutual Solubilities of Water and the [Cnmim][Tf2N] Hydrophobic Ionic Liquids

Ionic liquids (ILs) have recently garnered increased attention because of their potential environmental benefits as "green" replacements over conventional volatile organic solvents. While ILs cannot significantly volatilize and contribute to air pollution, even the most hydrophobic ones present some miscibility with water posing environmental risks to the aquatic ecosystems. Thus, the knowledge of ILs toxicity and their water solubility must be assessed before an accurate judgment of their environmental benefits and prior to their industrial applications. In this work, the mutual solubilities for [C2-C8mim][Tf2N] (n-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) and water between 288.15 and 318.15 K at atmospheric pressure were measured. Although these are among the most hydrophobic ionic liquids known, the solubility of water in these compounds is surprisingly large, ranging from 0.17 to 0.36 in mole fraction, while the solubility of these ILs in water is much lower ranging from 3.2 × 10-5 to 1.1 × 10-3 in mole fraction, in the temperature and pressure conditions studied. From the experimental data, the molar thermodynamic functions of solution and solvation such as Gibbs energy, enthalpy, and entropy at infinite dilution were estimated, showing that the solubility of these ILs in water is entropically driven. The predictive capability of COSMO-RS, a model based on unimolecular quantum chemistry calculations, was evaluated for the description of the binary systems investigated providing an acceptable agreement between the model predictions and the experimental data both with the temperature dependence and with the ILs structural variations.

Modulating the selectivity for CO and butane oxidation over heterogeneous catalysis through amorphous catalyst coatings

The preparation of porous films directly deposited onto the surface of catalyst particles is attracting increasing attention. We report here for the first time a method that can be carried out at ambient pressure for the preparation of porous films deposited over 3 mm diameter catalyst particles of silica-supported Pt-Fe. Characterization of the sample prepared at ambient pressure (i.e., open air, OA) and its main structural differences as compared with a Na-A (LTA) coated catalyst made using an autoclave-based method are presented. The OA-coated material predominantly exhibited an amorphous film over the catalyst surface with between 4 and 13% of crystallinity as compared with fully crystallized LTA zeolite crystals. This coated sample was highly selective for CO oxidation in the presence of butane with no butane oxidation observed up to 350 degrees C. This indicates, for the first time, that the presence of a crystalline membrane is not necessary for the difference in light off temperature between CO and butane to be achieved and that amorphous films may also produce this effect. An examination of the space velocity dependence and adsorption of Na+ on the catalysts indicates that the variation in CO and butane oxidation activity is not caused by site blocking predominantly, although the Pt activity was lowered by contact with this alkali.

Thermodynamics, structure, and dynamics in room temperature ionic liquids: The case of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6])

A detailed investigation of the phase diagram of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6]) is presented on the basis of a wide set of experimental data accessing thermodynamic, structural, and dynamical properties of this important room temperature ionic liquid (RTIL). The combination of quasi adiabatic, continuous calorimetry, wide angle neutron and X-ray diffraction, and quasi elastic neutron scattering allows the exploration of many novel features of this material. Thermodynamic and microscopic structural information is derived on both glassy and crystalline states and compared with results that recently appeared in the literature allowing direct information to be obtained on the existence of two crystalline phases that were not previously characterized and confirming the view that RTILs show a substantial degree of order (even in their amorphous states), which resembles the crystalline order. We highlight a strong connection between structure and dynamics, showing the existence of three temperature ranges in the glassy state across which both the spatial correlation and the dynamics change. The complex crystalline polymorphism in [bmim][PF6] also is investigated; we compare our findings with the corresponding findings for similar RTILs. These results provide a strong experimental basis for the exploration of the features of the phase diagram of RTILs and for the further study of longer alkyl chain salts.

Relaxation processes in room temperature ionic liquids: The case of 1-butyl-3-methyl imidazolium hexafluorophosphate

A detailed investigation on the nature of the relaxation processes occurring in a typical room temperature ionic liquid (RTIL), namely, 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6]), is reported. The study was conducted using both elastic and inelastic neutron scattering over a wide temperature range from 10 to 400 K, accessing the dynamic features of both the liquid and glassy amorphous states. In this study, the inelastic fixed energy scan technique has been applied for the first time to this class of materials. Using this technique, the existence of two relaxation processes below the glass transition and a further diffusive process occurring above the glass-liquid transition are observed. The low temperature processes are associated with methyl group rotation and butyl chain relaxation in the glassy state and have been modeled in terms of two Debye-like, Arrhenius activated processes. The high temperature process has been modeled in terms of a Kohlraush-Williams-Watts relaxation, with a distinct Vogel-Fulcher-Tamman temperature dependence. These results provide novel information that will be useful in rationalizing the observed structural and dynamical behavior of RTILs in the amorphous state.

A structural and electrochemical investigation of 1-alkyl-3-methylimidazollum salts of the nitratodioxouranate(VI) anions [{UO2(NO3)(2)}(2)(mu(4)-C2O4)](2-), [UO2(NO3)(3)](-), and [UO2(NO3)(4)](2-)

The properties of the 1-butyl-3-methylimidazolium salt of the dinuclear mu(4)-(O,O,O',O'-ethane-1,2-dioato)bis[bis-(nitrato-O,O)dioxouranate(VI)] anion have been investigated using electrochemistry, single-crystal X-ray crystallography, and extended X-ray absorbance fine structure spectroscopy: the anion structures from these last two techniques are in excellent agreement with each other. Electrochemical reduction of the complex leads to the a two-electron metal-centered reduction of U(VI) to U(IV), and the production Of UO2, or a complex containing UO2. Under normal conditions, this leads to the coating of the electrode with a passivating film. The presence of volatile organic compounds in the ionic liquids 1-alkyl-3-methylimidazolium nitrate (where the 1-alkyl chain was methyl, ethyl, propyl, butyl, pentyl, hexyl, dodecyl, hexadecyl, or octadecyl) during the oxidative dissolution of uranium(IV) oxide led to the formation of a yellow precipitate. To understand the effect of the cation upon the composition and structure of the precipitates, 1-alkyl-3-methylimidazolium salts of a number of nitratodioxouranate(VI) complexes were synthesized and then analyzed using X-ray crystallography. It was demonstrated that the length of the 1-alkyl chain played an important role, not only in the composition of the complex salt, but also in the synthesis of dinuclear anions containing the bridging mu(4)-(O,O,O',O'-ethane-1,2-dioato), or oxalato, ligand, by protecting it from further oxidation.

Voltammetry of oxygen in the room-temperature ionic liquids 1-ethyl-3-methylimidazolium bis((trifluoromethyl)sulfonyl)imide and hexyltriethylammonium bis((trifluoromethyl)sulfonyl)imide: One-electron reduction to form superoxide. Steady-state and transient behavior in the same cyclic voltammogram resulting from widely different diffusion coefficients of oxygen and superoxide

The electrochemical reduction of oxygen in two different room-temperature ionic liquids, 1-ethyl-3-methylimidazolium bis((trifluoromethyl)sulfonyl)imide ([EMIM][N(Tf)(2)]) and hexyltriethylammonium bis((trifluoromethyl)sulfonyl)imide ([N-6222][N(Tf)(2)]) was investigated by cyclic voltammetry at a gold microdisk electrode. Chronoamperometric measurements were made to determine the diffusion coefficient, D, and concentration, c, of the electroactive oxygen dissolved in the ionic liquid by fitting experimental transients to the Aoki model. [Aoki, K.; et al. J. Electroanal. Chem. 1981, 122, 19]. A theory and simulation designed for cyclic voltammetry at microdisk electrodes was then employed to determine the diffusion coefficient of the electrogenerated superoxide species, O-2(.-), as well as compute theoretical voltammograms to confirm the values of D and c for neutral oxygen obtained from the transients. As expected, the diffusion coefficient of the superoxide species was found to be smaller than that of the oxygen in both ionic liquids. The diffusion coefficients of O-2 and O-2(.-) in [N-6222][N(Tf)(2)], however, differ by more than a factor of 30 (D-O2 = 1.48 x 10(-10) m(2) s(-1), DO2.- = 4.66 x 10(-12) m(2) s(-1)), whereas they fall within the same order of magnitude in [EMIM][N(Tf)(2)] (D-O2 = 7.3 x 10(-10) m(2) s(-1), DO2.- = 2.7 x 10(-10) m(2) s(-1)). This difference in [N-6222][N(Tf)(2)] causes pronounced asymmetry in the concentration distributions of oxygen and superoxide, resulting in significant differences in the heights of the forward and back peaks in the cyclic voltammograms for the reduction of oxygen. This observation is most likely a result of the higher viscosity of [N-6222][N(Tf)(2)] in comparison to [EMIM][N(Tf)(2)], due to the structural differences in cationic component.

Temperature Dependence of the Primary Relaxation in 1-Hexyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}imide

We present results from complementary characterizations of the primary relaxation rate of a room temperature ionic liquid (RTIL), 1-hexyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl} imide, [C(6)mim][Tf2N], over a wide temperature range. This extensive data set is successfully merged with existing literature data for conductivity, viscosity, and NMR diffusion coefficients thus providing, for the case of RTILs, a unique description of the primary process relaxation map over more than 12 decades in relaxation rate and between 185 and 430 K. This unique data set allows a detailed characterization of the VTF parameters for the primary process, that are: B = 890 K, T-0 = 155.2 K, leading to a fragility index m = 71, corresponding to an intermediate fragility. For the first time neutron spin echo data from a fully deuteriated sample of RTIL at the two main interference peaks, Q = 0.76 and 1.4 angstrom(-1) are presented. At high temperature (T > 250 K), the collective structural relaxation rate follows the viscosity behavior; however at lower temperatures it deviates from the viscosity behavior, indicating the existence of a faster process.

Choline saccharinate and choline acesulfamate: Ionic liquids with low toxicities

Choline saccharinate and choline acesulfamate are two examples of hydrophilic ionic liquids, which can be prepared from easily available starting materials (choline chloride and a non-nutritive sweetener). The (eco)toxicity of these ionic liquids in aqueous solution is very low in comparison to other types of ionic liquids. A general method for the synthesis and purification of hydrophilic ionic liquids is presented. The method consists of a silver-free metathesis reaction, followed by purification of the ionic liquid by ion-exchange chromatography. The crystal structures show a marked difference in hydrogen bonding between the two ionic liquids, although the saccharinate and the acesulfamate anions show structural similarities. The optimized structures, the energetics, and the charge distribution of cation-anion pairs in the ionic liquids were studied by density functional theory (DFT) and second-order (Moller-Plesset) perturbation theory calculations. The occupation of the non-Lewis orbitals was considered to obtain a qualitative picture of the Lewis structures. The calculated interaction energies and the dipole moments for the ion pairs in the gas phase were discussed.

Electroconvection in a homeotropic bent-rod nematic liquid crystal beyond the dielectric inversion frequency

We characterize the structural transitions in an initially homeotropic bent-rod nematic liquid crystal excited by ac fields of frequency f well above the dielectric inversion point f(i). From the measured principal dielectric constants and electrical conductivities of the compound, the Carr-Helfrich conduction regime is anticipated to extend into the sub-megahertz region. Periodic patterned states occur through secondary bifurcations from the Freedericksz distorted state. An anchoring transition between the bend Freedericksz (1317) and degenerate planar (DP) states is detected. The BF state is metastable well above the Freedericksz threshold and gives way to the DP state, which persists in the field-off condition for several hours. Numerous +1 and -1 umbilics form at the onset of BF distortion, the former being largely of the chiral type. They survive in the DP configuration as linear defects, nonsingular in the core. In the BF regime, not far from fi, periodic Williams-like domains form around the umbilics; they drift along the director easy axis right from their onset. With increasing f, the wave vector of the periodic domains switches from parallel to normal disposition with respect to the c vector. Well above fi, a broadband instability is found.

Computational Verification of Two Universal Relations for Simple Ionic Liquids. Kinetic Properties of a Model 2:1 Molten Salt

Two semianalytical relations [Nature, 1996, 381, 137 and Phys. Rev. Lett. 2001, 87, 245901] predicting dynamical coefficients of simple liquids on the basis of structural properties have been tested by extensive molecular dynamics simulations for an idealized 2:1 model molten salt. In agreement with previous simulation studies, our results support the validity of the relation expressing the self-diffusion coefficient as a Function of the radial distribution functions for all thermodynamic conditions such that the system is in the ionic (ie., fully dissociated) liquid state. Deviations are apparent for high-density samples in the amorphous state and in the low-density, low-temperature range, when ions condense into AB(2) molecules. A similar relation predicting the ionic conductivity is only partially validated by our data. The simulation results, covering 210 distinct thermodynamic states, represent an extended database to tune and validate semianalytical theories of dynamical properties and provide a baseline for the interpretation of properties of more complex systems such as the room-temperature ionic liquids.