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.

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.

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.

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.

Thermodynamical studies of irreversible sorption of CO2 by Wyodak coal

Differential scanning calorimetry (DSC), temperature programmed desorption mass spectrometry (TPD-MS) and small angle neutron scattering (SANS) were used to investigate CO2 uptake by the Wyodak coal. The adsorption of carbon dioxide on Wyodak coal was studied by DSC. The exotherms evident at low temperatures are associated with the uptake of CO2 suggesting that carbon dioxide interacts strongly with the coal surface. The reduction in the value of the exotherms between the first and second runs for the Wyodak coal suggests that some CO2 is irreversibly bound to the structure even after heating to 200 °C DSC results also showed that adsorption of CO2 on the coal surface is an activated process and presumably at the temperature of the exotherms there is enough thermal energy to overcome the activation energy for adsorption. The adsorption process is instantly pursued by much slower diffusion of the gas molecules into the coal matrix (absorption). Structural rearrangement in coal by CO2 is examined by change in the glass transition temperature of coal after CO2 uptake at different pressures. The amount of gas dissolved in the coal increases with increasing CO2 pressure. TPD-MS showed that CO2 desorption from the Wyodak coal follows a first order kinetic model. Increase in the activation energy for desorption with pre-adsorbed CO2 pressure suggests that higher pressures facilitate the transport of CO2 molecules through the barriers therefore the amount of CO2 uptake by the coal is greater at higher pressures and more attempts are required to desorb CO2 molecules sorbed at elevated pressures. These conclusions were further confirmed by examining the Wyodak coal structure in high pressure CO 2 by SANS.

Fluid Shear Induces Conformation Change in Human Blood Protein von Willebrand Factor in Solution

Many of the physiological functions of von Willebrand Factor (VWF), including its binding interaction with blood platelets, are regulated by the magnitude of applied fluid/hydrodynamic stress. We applied two complementary strategies to study the effect of fluid forces on the solution structure of VWF. First, small-angle neutron scattering was used to measure protein conformation changes in response to laminar shear rates (G) up to 3000/s. Here, purified VWF was sheared in a quartz Couette cell and protein conformation was measured in real time over length scales from 2-140 nm. Second, changes in VWF structure up to 9600/s were quantified by measuring the binding of a fluorescent probe 1,1'-bis(anilino)-4-,4'-bis(naphtalene)-8,8'-disulfonate (bis-ANS) to hydrophobic pockets exposed in the sheared protein. Small angle neutron scattering studies, coupled with quantitative modeling, showed that VWF undergoes structural changes at G < 3000/s. These changes were most prominent at length scales <10 nm (scattering vector (q) range >0.6/nm). A mathematical model attributes these changes to the rearrangement of domain level features within the globular section of the protein. Studies with bis-ANS demonstrated marked increase in bis-ANS binding at G > 2300/s. Together, the data suggest that local rearrangements at the domain level may precede changes at larger-length scales that accompany exposure of protein hydrophobic pockets. Changes in VWF conformation reported here likely regulate protein function in response to fluid shear.