Phase transitions of choline dihydrogen phosphate: A vibrational spectroscopy and periodic DFT study.

Choline dihydrogen phosphate, [Chol][H2PO4], is a proton-conducting ionic plastic crystal exhibiting a complicated sequence of phase transitions. Here, we address the argument in the literature around the thermal properties of [Chol][H2PO4] using Raman and infrared microspectroscopy. The known structure of the low-temperature crystal, which contains the anti-conformer of [Chol]+ and hydrogen-bonded dimers of anions, was used to do periodic density functional theory calculations of the vibrational frequencies. Raman spectra indicate that the solid-solid transition at 20 °C is linked to a conformational change to the gauche [Chol] conformer with a concurrent local rearrangement of the anions. The distinct bands of lattice modes in the low-frequency range of the Raman spectra vanish at the 20 °C transition. Given the ease with which metastable crystals can be produced, Raman mappings demonstrate that a sample of [Chol][H2PO4] at ambient temperature can contain a combination of anti- and gauche conformers. Heating to 120 °C causes continuous changes in the local environment of anions rather than melting as suggested by a recent calorimetric investigation of [Chol][H2PO4]. The monotonic change in vibrational spectra is consistent with earlier observations of a very small entropy of fusion and no abrupt jump in the temperature dependence of ionic conductivity along the phase transitions of [Chol][H2PO4].

DFT and ab initio molecular dynamics simulation study of the infrared spectrum of the protic ionic liquid 2-hydroxyethylammonium formate.

Protic ionic liquids (PILs) typically show a complex band shape in their infrared (IR) spectra in the high-frequency range due to the hydrogen stretching vibrations of functional groups forming rather strong hydrogen bonds (H-bonds). In the low-frequency range, the intermolecular stretching mode of the H-bond leaves a mark in the far-IR spectrum of PILs. In this study, the IR spectrum of the PIL 2-hydroxyethylammonium formate, [HOCH2CH2NH3][HCOO], is investigated in order to identify the different modes that contribute to the high-frequency band shape, i.e. the cation ν(NH), ν(OH), and ν(CH) modes, and the anion ν(CH) mode, as well as the intermolecular mode of the strongest H-bond in the far-IR spectrum. The assignment is validated by quantum chemistry calculations of clusters at the density functional theory (DFT) level for four ionic pairs and by ab initio molecular dynamics (AIMD) simulations of ten ionic pairs. There is good agreement between the vibrational frequencies obtained from DFT and AIMD simulations for both the high- and low-frequency ranges. Based on the calculations, the strong H-bond interaction between the cation -NH3 group and [HCOO]- gives a broad band envelope associated with the ν(NH) mode in the high-frequency range of the IR spectrum on which there are narrower peaks corresponding to the ν(OH) and ν(CH) modes. In the far-IR (FIR) spectrum, the anions' rattling motion gives a broad feature with a maximum at 160 cm-1, while the H-bond's intermolecular NH⋯O stretching mode appears as a peak at 255 cm-1.

Confined ionic liquids films under shear: The importance of the chemical nature of the solid surface.

Ionic liquids have generated interest in applications as lubricants and as additives to conventional lubricants due to their unique physical properties. In these applications, the liquid thin film can be subjected simultaneously to extremely high shear and loads in addition to nanoconfinement effects. Here, we use molecular dynamics simulations with a coarse-grained model to study a nanometric film of an ionic liquid confined between two planar solid surfaces both at equilibrium and at several shear rates. The strength of the interaction between the solid surface and the ions was changed by simulating three different surfaces with enhanced interactions with different ions. The increase in the interaction with either the cation or the anion leads to the formation of a solid-like layer that moves alongside the substrates; however, this layer can exhibit different structures and stability. An increase in the interaction with the high symmetry anion produces a more regular structure that is more resistant to the effects of shear and viscous heating. Two definitions were proposed and used for the calculation of the viscosity: a local definition based on the microscopic characteristics of the liquid and an engineering definition based on the forces measured at the solid surfaces, with the former displaying a correlation with the layered structure induced by the surfaces. Because of the shear thinning behavior of the ionic liquids as well as the temperature rise brought on by viscous heating, both the engineering and the local viscosities decrease as the shear rate increases.

Role of density and electrostatic interactions in the viscosity and non-newtonian behavior of ionic liquids - a molecular dynamics study.

Strong ionic interactions, as well as the consequent correlations between cation and anion dynamics, give ionic liquids various physical features that set them apart from ordinary organic solvents. In particular, they result in larger viscosities and larger densities than mixtures of neutral compounds with similar molecular structures. However, both the direct effect of electrostatic interactions and the increase of liquid density contribute to the high viscosity and so far no experimental or computational work enabled a clear quantification of those effects. Also, the effects over the shear thinning behavior, which may have important consequences for application as lubricants, were not considered yet. Here, these questions were tackled by performing non-equilibrium molecular dynamics (NEMD) simulations changing both the strength of ionic interactions and liquid density at several shear rates using a coarse grained model. The relative dielectric constant was adjusted to reproduce viscosity data from all-atoms simulations on both zero shear and high shear conditions. Elimination of ionic interactions results in a reduction of density and zero shear viscosity and also delays the beginning of shear thinning to higher shear rates. Restoring density to the ionic liquid's value only partially reverses the alterations. Correlations of the non-newtonian behavior and changes in the intermolecular structure and contact lifetimes were also explored.

Relating the structure and dynamics of ionic liquids under shear by means of reverse non-equilibrium molecular dynamics simulations.

The effect of the shear rate on the viscosity and the structure of 1-ethyl-3-methylimidazolium based ionic liquids with three different anions (tetrafluoroborate, dicyanamide, and bis(trifluoromethylsulfonyl)imide) was studied by means of reverse non-equilibrium molecular dynamics (RNEMD) simulations using a polarizable force field. The three liquids display a Newtonian plateau followed by a shear thinning regime at shear rates of the order of GHz. Even though the main features of the liquid structure remains under shear, systematic changes were noticed at the GHz rates, with coordination shells becoming more diffuse as noticed by the reduction in the difference between consecutive maxima and minima in the radial distribution function. Interestingly, these structural changes with the shear rate can be precisely fitted using the Carreau equation, which is a well-known expression for the shear rate dependence of the viscosity. The fitting parameters for different distributions can be used to explain qualitatively the shear thinning behavior of these liquids. In the GHz range, the cations and, in a minor extension, some anions, tend to assume preferentially a parallel orientation with the flux, which contributes to the shear thinning behavior and may have consequences for adhesion in applications as lubricants.

An inelastic neutron scattering, Raman, far-infrared, and molecular dynamics study of the intermolecular dynamics of two ionic liquids.

The intermolecular dynamics in the THz frequency range of the ionic liquids n-butyl-trimethylammonium bis(trifluoromethanesulfonyl)imide, [N1114][NTf2], and methyl-tributylammonium bis(trifluoromethanesulfonyl)imide, [N1444][NTf2], were investigated by a combined usage of inelastic neutron scattering (INS), Raman, and far-infrared (FIR) spectroscopies and the power spectrum calculated by molecular dynamics (MD) simulations. The collective dynamics of the simulated systems is also discussed by the calculation of time correlation functions of charge and mass currents that are projected onto acoustic- and optic-like motions. The INS and Raman measurements have been performed as a function of temperature in the glassy, crystalline, and liquid phases. The excess in the vibrational density of states over the expectation of the Debye theory, the so-called boson peak, is found in the INS and Raman spectra as a peak at ∼2 meV (∼16 cm-1) and also in the direct measurement of heat capacity at very low temperatures (4-20 K). This low-frequency vibration is incorporated into the curve fits of Raman, FIR, and MD data at room temperature. Fits of spectra from these different sources in the range below 100 cm-1 are consistently achieved with three components at ca. 25, 50, and 80 cm-1, but with distinct relative intensities among the different techniques. It is proposed as the collective nature of the lowest-frequency component and the anion-cation intermolecular vibration nature of the highest-frequency component. The MD results indicate that there is no clear distinction between acoustic and optic vibrations in the spectral range investigated in this work for the ionic liquids [N1114][NTf2] and [N1444][NTf2]. The analysis carried out here agrees in part, but not entirely, with other propositions in the literature, mainly from optical Kerr effect (OKE) and FIR spectroscopies, concerning the intermolecular dynamics of ionic liquids.

Effect of alkyl-group flexibility on the melting point of imidazolium-based ionic liquids.

The low melting point of room temperature ionic liquids is usually explained in terms of the presence of bulky, low-symmetry, and flexible ions, with the first two factors related to the lattice energy while an entropic effect is attributed to the latter. By means of molecular dynamics simulations, the melting points of 1-ethyl-3-methyl-imidazolium hexafluorophosphate and 1-decyl-3-methyl-imidazolium hexafluorophosphate were determined, and the effect of the molecular flexibility over the melting point was explicitly computed by restraining the rotation of dihedral angles in both the solid and the liquid phases. The rotational flexibility over the bond between the ring and the alkyl chain affects the relative ordering of the anions around the cations and results in substantial effects over both the enthalpy and the entropy of melting. For the other dihedral angles of the alkyl group, the contributions are predominantly entropic and an alternating behavior was found. The flexibility of some dihedral angles has negligible effects on the melting point, while others can lead to differences in the melting point as large as 20 K. This alternating behavior is rationalized by the different probabilities of conformation defects in the crystal.

Ion pair free energy surface as a probe of ionic liquid structure.

Numerous combinations of cations and anions are possible for the production of ionic liquids with fine-tuned properties once the correlation with the molecular structure is known. In this sense, computer simulations are useful tools to explain and even predict the properties of ionic liquids. However, quantum mechanical methods are usually restricted to either small clusters or short time scales so that parameterized force fields are required to study the bulk liquids. In this work, a method is proposed to enable a comparison between the quantum mechanical system and both polarizable and nonpolarizable force fields by means of the calculation of free energy surfaces for the translational motion of the anion around the cation in gas phase. This method was tested for imidazolium-based cations with 3 different anions, [BF4]-, [N(CN)2]-, and [NTf2]-. Better agreement was found with the density functional theory calculations when polarizability is introduced in the force field. In addition, the ion pair free energy surfaces reproduced the main structural patterns observed in the first coordination shell in molecular dynamics simulations of the bulk liquid, proving to be useful probes for the liquid phase structure that can be computed with higher level methods and the comparison with forcefields can indicate further improvements in their parameterization.

Vibrational Spectroscopy of Ionic Liquids.

Vibrational spectroscopy has continued use as a powerful tool to characterize ionic liquids since the literature on room temperature molten salts experienced the rapid increase in number of publications in the 1990's. In the past years, infrared (IR) and Raman spectroscopies have provided insights on ionic interactions and the resulting liquid structure in ionic liquids. A large body of information is now available concerning vibrational spectra of ionic liquids made of many different combinations of anions and cations, but reviews on this literature are scarce. This review is an attempt at filling this gap. Some basic care needed while recording IR or Raman spectra of ionic liquids is explained. We have reviewed the conceptual basis of theoretical frameworks which have been used to interpret vibrational spectra of ionic liquids, helping the reader to distinguish the scope of application of different methods of calculation. Vibrational frequencies observed in IR and Raman spectra of ionic liquids based on different anions and cations are discussed and eventual disagreements between different sources are critically reviewed. The aim is that the reader can use this information while assigning vibrational spectra of an ionic liquid containing another particular combination of anions and cations. Different applications of IR and Raman spectroscopies are given for both pure ionic liquids and solutions. Further issues addressed in this review are the intermolecular vibrations that are more directly probed by the low-frequency range of IR and Raman spectra and the applications of vibrational spectroscopy in studying phase transitions of ionic liquids.

Low-frequency Raman spectra of a glass-forming ionic liquid at low temperature and high pressure.

The frequency range below ∼100 cm-1 of the Raman spectrum of a glass-forming liquid exhibits two features that characterize the short-time (THz) dynamics: the quasi-elastic scattering (QES) tail and the boson peak (BP). In this work, we follow temperature and pressure effects on the intermolecular dynamics of a typical ionic liquid, 1-butyl-1-methylpiperidinium bis(trifluoromethanesulfonyl)imide, [Pip14][[NTf2]. The glass transition temperature of [Pip14][[NTf2] at atmospheric pressure is Tg = 198 K, and the pressure of glass transition at room temperature is Pg = 1.1 GPa. Raman spectra obtained while cooling the liquid or heating the glass exhibit hysteresis in QES and BP intensities, IQES and IBP. The dependence of IQES, IBP, and the BP frequency, ωBP, with pressure up to the glass transition is steeper than the temperature dependence due to the stronger pressure effect on density within the GPa range. The temperature and pressure behaviors of the parameters IQES, IBP, and ωBP obtained here for [Pip14][[NTf2] are discussed in light of known results for other glass-formers.

Spatial and thermal signatures of α and ß relaxations in glassy and glacial aliphatic ionic liquids.

The competition between Coulomb and van der Waals interactions brings forth unique dynamic features and broad applications to ionic liquids. Herein, we present a combined calorimetric, X-ray diffraction, incoherent elastic, and quasi-elastic neutron scattering study, over a wide temperature range (180-340 K), of the relaxational dynamics of the liquid, supercooled liquid, crystalline, glassy, and glacial states of two model ionic liquids: tributylmethylammonium (a good glass-former) and butyltrimethylammonium (a good crystal-former) cations and the bis(trifluoromethanesulfonyl)imide anion. In both systems, we observed two distinct relaxation processes. The Q-dependence of the respective relaxation time shows that the α-process is diffusive, while the ß-process is modulated by the structure of the liquids.

Structure and Liquid Fragility in Sodium Carbonate.

The relationship between local structure and dynamics is explored for molten sodium carbonate. A flexible fluctuating-charge model, which allows for changes in the shape and charge distribution of the carbonate molecular anion, is developed. The system shows the evolution of highly temperature-dependent complex low-dimensional structures which control the dynamics (and hence the liquid fragility). By varying the molecular anion charge distribution, the key interactions responsible for the formation of these structures can be identified and rationalized. An increase in the mean charge separation within the carbonate ions increases the connectivity of the emerging structures and leads to an increase in the system fragility.

Mechanical heterogeneity in ionic liquids.

Molecular dynamics (MD) simulations of five ionic liquids based on 1-alkyl-3-methylimidazolium cations, [C n C1im]+, have been performed in order to calculate high-frequency elastic moduli and to evaluate heterogeneity of local elastic moduli. The MD simulations of [C n C1im][NO3], n = 2, 4, 6, and 8, assessed the effect of domain segregation when the alkyl chain length increases, and [C8C1im][PF6] assessed the effect of strength of anion-cation interaction. Dispersion curves of excitation energies of longitudinal and transverse acoustic, LA and TA, modes were obtained from time correlation functions of mass currents at different wavevectors. High-frequency sound velocity of LA modes depends on the alkyl chain length, but sound velocity for TA modes does not. High-frequency bulk and shear moduli, K ∞ and G ∞ , depend on the alkyl chain length because of a density effect. Both K ∞ and G ∞ are strongly dependent on the anion. The calculation of local bulk and shear moduli was accomplished by performing bulk and shear deformations of the systems cooled to 0 K. The simulations showed a clear connection between structural and elastic modulus heterogeneities. The development of nano-heterogeneous structure with increasing length of the alkyl chain in [C n C1im][NO3] implies lower values for local bulk and shear moduli in the non-polar domains. The mean value and the standard deviations of distributions of local elastic moduli decrease when [NO3]- is replaced by the less coordinating [PF6]- anion.

Communication: Glass transition and melting lines of an ionic liquid.

The phase diagram of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethanesufonyl)imide, [Pyrr1,4][NTf2], was explored by synchroton X-ray diffraction and Raman scattering measurements as a function of temperature and pressure. Glass transition Tg(p) and melting Tm(p) temperatures were obtained from atmospheric pressure up to ca. 2.0 GPa. We found that both the Tg(p) and Tm(p) curves follow essentially the same pressure dependence. The similarity of pressure coefficients, dTg/dp ≈ dTm/dp, is explained within the non-equilibrium thermodynamics approach for the glass transition by assuming that one of the Ehrenfest equations is appropriated for Tg(p), whereas Tm(p) follows the Clausius-Clapeyron equation valid for the first-order transitions. The results highlight that ionic liquids are excellent model systems to address fundamental questions related to the glass transition.

The structure of liquid alkali nitrates and nitrites.

High energy X-ray diffraction has been combined with containerless techniques to determine the structure of a series of alkali and ammonium nitrate and nitrite liquids. The systems have been modelled using molecular dynamics simulation which allows for the flexibility of, and movement of charge within, the molecular anions. The model reproduces the experimentally-determined scattering functions in both the low- and high-Q regimes reflecting the inter- and intra-molecular length-scales. For ammonium nitrate the best fit to the diffraction data is obtained by assuming the NH4+ cation to have a radius closer to that for Cs+ rather than a smaller cation such as Rb+ as often previously assumed. The alkali nitrites show an emergent length scale, attributed to the nitrogen-nitrogen spatial correlations, that depends on both temperature and the identity of the alkali cation. The corresponding nitrates show a more subtle effect in the nitrogen-nitrogen correlations. As a result, the nature of this N-N length-scale appears different for the respective nitrites and nitrates.

Comparing two tetraalkylammonium ionic liquids. I. Liquid phase structure.

X-ray scattering experiments at room temperature were performed for the ionic liquids n-butyl-trimethylammonium bis(trifluoromethanesulfonyl)imide, [N1114][NTf2], and methyl-tributylammonium bis(trifluoromethanesulfonyl)imide, [N1444][NTf2]. The peak in the diffraction data characteristic of charge ordering in [N1444][NTf2] is shifted to longer distances in comparison to [N1114][NTf2], but the peak characteristic of short-range correlations is shifted in [N1444][NTf2] to shorter distances. Molecular dynamics (MD) simulations were performed for these ionic liquids using force fields available from the literature, although with new sets of partial charges for [N1114](+) and [N1444](+) proposed in this work. The shifting of charge and adjacency peaks to opposite directions in these ionic liquids was found in the static structure factor, S(k), calculated by MD simulations. Despite differences in cation sizes, the MD simulations unravel that anions are allowed as close to [N1444](+) as to [N1114](+) because anions are located in between the angle formed by the butyl chains. The more asymmetric molecular structure of the [N1114](+) cation implies differences in partial structure factors calculated for atoms belonging to polar or non-polar parts of [N1114][NTf2], whereas polar and non-polar structure factors are essentially the same in [N1444][NTf2]. Results of this work shed light on controversies in the literature on the liquid structure of tetraalkylammonium based ionic liquids.

Structure of cyano-anion ionic liquids: X-ray scattering and simulations.

Ionic liquids with cyano anions have long been used because of their unique combination of low-melting temperatures, reduced viscosities, and increased conductivities. Recently we have shown that cyano anions in ionic liquids are particularly interesting for their potential use as electron donors to excited state photo-acceptors [B. Wu et al., J. Phys. Chem. B 119, 14790-14799 (2015)]. Here we report on bulk structural and quantum mechanical results for a series of ionic liquids based on the 1-ethyl-3-methylimidazolium cation, paired with the following five cyano anions: SeCN(-), SCN(-), N(CN)2 (-), C(CN)3 (-), and B(CN)4 (-). By combining molecular dynamics simulations, high-energy X-ray scattering measurements, and periodic boundary condition DFT calculations, we are able to obtain a comprehensive description of the liquid landscape as well as the nature of the HOMO-LUMO states for these ionic liquids in the condensed phase. Features in the structure functions for these ionic liquids are somewhat different than the commonly observed adjacency, charge-charge, and polarity peaks, especially for the bulkiest B(CN)4 (-) anion. While the other four cyano-anion ionic liquids present an anionic HOMO, the one for Im2,1 (+)/B(CN)4 (-) is cationic.

Comparing two tetraalkylammonium ionic liquids. II. Phase transitions.

Phase transitions of the ionic liquids n-butyl-trimethylammonium bis(trifluoromethanesulfonyl)imide, [N1114][NTf2], and methyl-tributylammonium bis(trifluoromethanesulfonyl)imide, [N1444][NTf2], were investigated by differential scanning calorimetry (DSC), X-ray diffraction (XRD) measurements, and Raman spectroscopy. XRD and Raman spectra were obtained as a function of temperature at atmospheric pressure, and also under high pressure at room temperature using a diamond anvil cell (DAC). [N1444][NTf2] experiences glass transition at low temperature, whereas [N1114][NTf2] crystallizes or not depending on the cooling rate. Both the ionic liquids exhibit glass transition under high pressure. XRD and low-frequency Raman spectra provide a consistent physical picture of structural ordering-disordering accompanying the thermal events of crystallization, glass transition, cold crystallization, pre-melting, and melting. Raman spectra in the high-frequency range of some specific cation and anion normal modes reveal conformational changes of the molecular structures along phase transitions.

Two phosphonium ionic liquids with high Li(+) transport number.

This work presents the physicochemical characterization of two ionic liquids (ILs) with small phosphonium cations, triethylpenthylphosphonium bis(trifluoromethanesulfonyl)imide ([P2225][Tf2N]) and (2-methoxyethyl)trimethylphosphonium bis(trifluoromethanesulfonyl)imide ([P222(201)][Tf2N]), and their mixtures with Li(+). Properties such as the electrochemical window, density, viscosity and ionic conductivity are presented. The diffusion coefficient was obtained using two different techniques, PGSE-NMR and Li electrodeposition with microelectrodes. In addition, the Li(+) transport number was calculated using the PGSE-NMR technique and an electrochemical approach. The use of these three techniques showed that the PGSE-NMR technique underestimates the diffusion coefficient for charged species. The Li(+) transport number was found to be as high as 0.54. Raman spectroscopy and molecular dynamics simulations were used to evaluate the short-range structure of the liquids. These experiments suggested that the interaction between the Li(+) and the Tf2N(-) anion is similar to that seen with other ILs containing the same anion. However, the MD simulations also showed that the Li(+) ions interact differently with the cation containing an alkyl ether chain. The results found in this work suggest that these Li(+) mixtures have promising potential to be applied as electrolytes in batteries.

Glass transition of ionic liquids under high pressure.

The glass transition pressure at room temperature, pg, of six ionic liquids based on 1-alkyl-3-methylimidazolium cations and the anions [BF4](-), [PF6](-), and bis(trifluromethanesulfonyl)imide, [NTf2](-), has been obtained from the pressure dependence of the bandwidth of the ruby fluorescence line in diamond anvil cells. Molar volume, Vm(pg), has been estimated by a group contribution model (GCM) developed for the ionic liquids. A density scaling relation, TV(γ), has been considered for the states Vm(pg, 295 K) and Vm(Tg, 0.1 MPa) using the simplifying condition that the viscosity at the glass transition is the same at pg at room temperature and at atmospheric pressure at Tg. Assuming a constant γ over this range of density, a reasonable agreement has been found for the γ determined herein and that of a previous density scaling analysis of ionic liquids viscosities under moderate conditions. Further support for the appropriateness of extrapolating the GCM equation of state to the GPa pressure range is provided by comparing the GCM and an equation of state previously derived in the power law density-scaling regime.