Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Molecular electrostatic properties of ions in an ionic liquid

We have analysed the electronic wave functions from an ab initio simulation of the ionic liquid (room temperature molten salt) dimethyl imidazolium chloride ([dmim][Cl] or [C1mim][Cl]) using localized Wannier orbitals. This allows us to assign electron density to individual ions. The probability distributions of the ionic dipole moments for an isolated ion and for ions in solution are compared. The liquid environment is found to polarize the cation by about 0.7 D and to increase the amplitude of the fluctuations in the dipole moments of both cation and anion. The relative changes in nuclear and electronic contributions are shown. The implications for classical force fields are discussed.

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.

Identical extraction behavior and coordination of trivalent or hexavalent f-element cations using ionic liquid and molecular solvents

The extraction of both UO22+ and trivalent lanthanide and actinide ions (Am3+, Nd3+, Eu3+) by dialkylphosphoric or dialkylphosphinic acids from aqueous solutions into the ionic liquid, 1-decyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide has been studied and compared to extractions into dodecane. Radiotracer partitioning measurements show comparable patterns of distribution ratios for both the ionic liquid/aqueous and dodecane/aqueous systems, and the limiting slopes at low acidity indicate the partitioning of neutral complexes in both solvent systems. The metal ion coordination environment, elucidated from EXAFS and UV-visible spectroscopy measurements, is equivalent in the ionic liquid and dodecane solutions with coordination of the uranyl cation by two hydrogen-bonded extractant dimers, and of the trivalent cations by three extractant dimers. This is the first definitive report of a system where both the biphasic extraction equilibria and metal coordination environment are the same in an ionic liquid and a molecular organic solvent.

Non-ideal behaviour of a room temperature ionic liquid in an alkoxyethanol or poly ethers at T = (298.15 to 318.15) K

Non-ideal behaviour of 1-butyl-3-methylimidazolium hexafluorophosphate [bmim][PF6] in ethylene glycol monomethyl ether; CH3OCH2CH2OH (EGMME), ethylene glycol dimethyl ether; CH3OCH2CH2OCH3 (EGDME) and diethylene glycol dimethyl ether; CH3(OCH2CH2)2OCH3 (DEGDME) have been investigated over the whole composition range at T = (298.15 to 318.15) K. To gain insight into the mixing behaviour, results of density measurements were used to estimate excess molar volumes, image, apparent molar volumes, Vphi,i, partial molar volumes, image, excess partial molar volumes, image, and their limiting values at infinite dilution, image, image, and image, respectively. Volumetric results have been analyzed in the light of Prigogine–Flory–Patterson (PFP) statistical mechanical theory. Measurements of refractive indices n were also performed for all the binary mixtures over whole composition range at T = 298.15 K. Deviations in refractive indices ?phin and the deviation of molar refraction ?xR have been calculated from experimental data. Refractive indices results have been correlated with volumetric results and have been interpreted in terms of molecular interactions. Excess properties are fitted to the Redlich–Kister polynomial equation to obtain the binary coefficients and the standard errors.

Application of static charge transfer within an ionic-liquid force field and its effect on structure and dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3-dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are +/- 0.5 e. While the three-dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly-with the positive result that for ionic charges of +/- 0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial-charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic-liquid models in which the ionic charges are assumed to be +/- 1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled-charge sets may be of benefit in the simulation of these systems-especially when looking at properties beyond liquid structure-thus providing on alternative to computationally expensive polarisable force fields.

Density-viscosity product of small-volume ionic liquid samples using quartz crystal impedance analysis

Quartz crystal impedance analysis has been developed as a technique to assess whether room-temperature ionic liquids are Newtonian fluids and as a small-volume method for determining the values of their viscosity-density product, rho eta. Changes in the impedance spectrum of a 5-MHz fundamental frequency quartz crystal induced by a water-miscible room-temperature ionic liquid, 1-butyl-3-methylimidazolium. trifluoromethylsulfonate ([C(4)mim][OTf]), were measured. From coupled frequency shift and bandwidth changes as the concentration was varied from 0 to 100% ionic liquid, it was determined that this liquid provided a Newtonian response. A second water-immiscible ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C(4)mim][NTf2], with concentration varied using methanol, was tested and also found to provide a Newtonian response. In both cases, the values of the square root of the viscosity-density product deduced from the small-volume quartz crystal technique were consistent with those measured using a viscometer and density meter. The third harmonic of the crystal was found to provide the closest agreement between the two measurement methods; the pure ionic liquids had the largest difference of similar to 10%. In addition, 18 pure ionic liquids were tested, and for 11 of these, good-quality frequency shift and bandwidth data were obtained; these 12 all had a Newtonian response. The frequency shift of the third harmonic was found to vary linearly with square root of viscosity-density product of the pure ionic liquids up to a value of root(rho eta) approximate to 18 kg m(-2) s(-1/2), but with a slope 10% smaller than that predicted by the Kanazawa and Gordon equation. It is envisaged that the quartz crystal technique could be used in a high-throughput microfluidic system for characterizing ionic liquids.

Electrode kinetics and mechanism of iodine reduction in the room-temperature ionic liquid [C(4)mim][NTf2]

The fast electrochemical reduction of iodine in the RTIL 1-butyl-3-methylimidazolium bis(trifluoromethyl-sulfonyl)imide, [C(4)mim][NTf2], is reported and the kinetics and mechanism of the process elucidated. Two reduction peaks were observed. The first reduction peak is assigned to the process