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.

Liquid structure of the ionic liquid, 1-methyl-4-cyanopyridinium bis{(trifluoromethyl)sulfonyl}imide determined from neutron scattering and molecular dynamics simulations

The liquid structure of 1-methyl-4-cyanopyridinium bis {(trifluoromethyl)sulfonyl}imide, a prototypical ionic liquid containing an electron-withdrawing group on the cation, has been investigated at 368 K. Experimental neutron scattering combined with empirical potential structure refinement analysis of the data and classical molecular dynamics simulations have been used to probe the liquid structure in detail. Both techniques generated highly consistent results that provide valuable validation of the force fields and refinement approaches. A significant degree of apparent charge ordering is found in the liquid structure, although the nonspherical shape of the ions results in interpenetration of cations into the first shell of adjacent cations, with much shorter closest contact distances than the averaged center-of-mass cation-cation and cation-anion separations.

Small angle neutron scattering from 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids ([C(n)mim][PF6], n=4, 6, and 8)

The presence of local anisotropy in the bulk, isotropic, and ionic liquid phases-leading to local mesoscopic inhomogeneity-with nanoscale segregation and expanding nonpolar domains on increasing the length of the cation alkyl-substituents has been proposed on the basis of molecular dynamics (MD) simulations. However, there has been little conclusive experimental evidence for the existence of intermediate mesoscopic structure between the first/second shell correlations shown by neutron scattering on short chain length based materials and the mesophase structure of the long chain length ionic liquid crystals. Herein, small angle neutron scattering measurements have been performed on selectively H/D-isotopically substituted 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids with butyl, hexyl, and octyl substituents. The data show the unambiguous existence of a diffraction peak in the low-Q region for all three liquids which moves to longer distances (lower Q), sharpens, and increases in intensity with increasing length of the alkyl substituent. It is notable, however, that this peak occurs at lower values of Q (longer length scale) than predicted in any of the previously published MD simulations of ionic liquids, and that the magnitude of the scattering from this peak is comparable with that from the remainder of the amorphous ionic liquid. This strongly suggests that the peak arises from the second coordination shells of the ions along the vector of alkyl-chain substituents as a consequence of increasing the anisotropy of the cation, and that there is little or no long-range correlated nanostructure in these ionic liquids.

Study of structural change in Wyodak coal in high pressure CO2 by small-angle neutron scattering

Small angle neutron scattering (SANS) has been applied to examine the effect of high pressure CO2 on the structure of Wyodak coal. Significant decrease in the scattering intensities upon exposure of the coal to high pressure CO2 showed that high pressure CO2 rapidly adsorbs on the coal and reaches to all pores in the structure. This is confirmed by strong and steep exothermic peaks observed on DSC scans during coal/ CO2 interactions. In situ small angle neutron scattering on coal at high pressure CO2 atmosphere showed an increase in scattering intensities with time suggesting that after adsorption, high pressure CO2 immediately begins to diffuse into the coal matrix, changes the macromolecular structure of the coal, swells the matrix and probably creates microporosity in coal structure by extraction of volatile components from coal. Significant decrease in the glass transition temperature of coal caused by high pressure CO2 also confirms that CO2 at elevated pressures dissolve in the coal matrix, results in significant plasticization and physical rearrangement of the coal’s macromolecular structure.

Association and liquid structure of pyridine-acetic acid mixtures determined from neutron scattering using a ‘free proton’ EPSR simulation model

The liquid structure of pyridine-acetic acid mixtures have been investigated using neutron scattering at various mole fractions of acetic acid, χHOAc = 0.33, 0.50, and 0.67, and compared to the structures of neat pyridine and acetic acid. Data has been modelled using Empirical Potential Structure Refinement (EPSR) with a ‘free proton’ reference model, which has no prejudicial weighting towards either the existence of molecular or ionised species. Analysis of the neutron scattering results shows the existence of hydrogen-bonded acetic acid chains with pyridine inclusions, rather than the formation of an ionic liquid by proton transfer.

Neutron Scattering of Aromatic and Aliphatic Liquids

Organic solvents, such as cyclohexane, cyclohexene, methylcyclohexane, benzene and toluene, are widely used as both reagents and solvents in industrial processes. Despite the ubiquity of these liquids, the local structures that govern the chemical properties have not been studied extensively. Herein, we report neutron diffraction measurements on liquid cyclohexane, cyclohexene, methylcyclohexane, benzene and toluene at 298 K to obtain a detailed description of the local structure in these compounds. The radial distribution functions of the centres of the molecules, as well as the partial distribution functions for the double bond for cyclohexene and methyl group for methylcyclohexane and toluene have been calculated. Additionally, probability density functions and angular radial distribution functions were extracted to provide a full description of the local structure within the chosen liquids. Structural motifs are discussed and compared for all liquids, referring specifically to the functional group and aromaticity present in the different liquids.

Evolution of elastic x-ray scattering in laser-shocked warm dense lithium

We have studied the dynamics of warm dense Li with near-elastic x-ray scattering. Li foils were heated and compressed using shock waves driven by 4-ns-long laser pulses. Separate 1-ns-long laser pulses were used to generate a bright source of 2.96 keV Cl Ly-alpha photons for x-ray scattering, and the spectrum of scattered photons was recorded at a scattering angle of 120 degrees using a highly oriented pyrolytic graphite crystal operated in the von Hamos geometry. A variable delay between the heater and backlighter laser beams measured the scattering time evolution. Comparison with radiation-hydrodynamics simulations shows that the plasma is highly coupled during the first several nanoseconds, then relaxes to a moderate coupling state at later times. Near-elastic scattering amplitudes have been successfully simulated using the screened one-component plasma model. Our main finding is that the near-elastic scattering amplitudes are quite sensitive to the mean ionization state Z and by extension to the choice of ionization model in the radiation-hydrodynamics simulations used to predict plasma properties within the shocked Li.

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.

Neutron diffraction, NMR and molecular dynamics study of glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate

beta-D-glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate in a 6 : 1 molar ratio (ionic liquid : glucose) has been studied by neutron scattering, NMR and molecular dynamics simulations. Good agreement was found between simulated neutron scattering profiles generated for isotopically substituted liquid systems and those experimentally determined as well as between simulated and experimental diffusion coefficients obtained by Pulsed Field Gradient NMR spectroscopy. The overriding glucose-ionic liquid interactions in the liquid are hydrogen-bonding between acetate oxygens and sugar hydroxyl groups. The ionic liquid cation was found to play only a minor role in the solvation of the sugar and does not participate in hydrogen-bonding with the sugar to any significant degree. NOESY experiments lend further evidence that there is no direct interaction between sugar hydroxyl groups and acidic hydrogens on the ionic liquid cation.

Local order and frustration in the geometrically frustrated spinels Cd1-xZnxV2O4

Orbitally degenerate frustrated spinels, Cd1-xZnxV2O4, with 0 <= x <= 1 were investigated using elastic and inelastic neutron scattering techniques. In the end members with x=0 and 1, a tetragonal distortion (c < a) has been observed upon cooling mediated by a Jahn-Teller distortion that gives rise to orbital ordering. This leads to the formation of spin chains in the ab-plane that upon further cooling, Neel ordering is established due to interchain coupling. In the doped compositions, however, the bulk susceptibility, chi, shows that the macroscopic transitions to cooperative orbital ordering and long-range antiferromagnetic ordering are suppressed. However, the inelastic neutron scattering measurements suggest that the dynamic spin correlations at low temperatures have similar one-dimensional characteristics as those observed in the pure samples. The pair density function analysis of neutron diffraction data shows that the local atomic structure does not become random with doping but rather consists of two distinct environments corresponding to ZnV2O4 and CdV2O4. This indicates that short-range orbital ordering is present which leads to the one-dimensional character of the spin correlations even in the low temperature cubic phase of the doped compositions.

Determination of toluene hydrogenation kinetics with neutron diffraction

Total neutron scattering has been used to follow the hydrogenation of toluene-d₈ to methylcyclohexane-d₁₄ over 3 wt% platinum supported on highly ordered mesoporous silica (MCM-41) at 298 K and under 150 mbar D₂ pressure. The detailed kinetic information so revealed indicates that liquid reorganisation inside pores is the slowest step of the whole process. Additionally, the results were compared with the reaction performed under 250 mbar D₂ pressure as well as with toluene-h₈ hydrogenation using D₂ at 150 mbar.

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.

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.