RESUMO
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].
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.
Assuntos
Solventes Eutéticos Profundos/química , Biomassa , Indústria Farmacêutica , Eletroquímica , Aprendizado de Máquina , Simulação de Dinâmica Molecular , Nanotecnologia , Teoria QuânticaRESUMO
Polarizable force fields are gradually becoming a common choice for ionic soft matter, in particular, for molecular dynamics (MD) simulations of ionic liquids (ILs) and deep eutectic solvents (DESs). The CL&Pol force field introduced in 2019 is the first general, transferable, and polarizable force field for MD simulations of different types of DESs. The original formulation contains, however, some problems that appear in simulations of ethaline and may also have a broader impact. First, the originally proposed atomic diameter parameters are unbalanced, resulting in too weak interactions between the chlorides and the hydroxyl groups of the ethylene glycol molecules. This, in turn, causes an artificial phase separation in long simulations. Second, there is an overpolarization of chlorides due to strong induced dipoles that give rise to the presence of peaks and antipeaks at very low q-vector values (2.4 nm-1) in the partial components of the structure factors. In physical terms, this is equivalent to overestimated spatial nanoscale heterogeneity. To correct these problems, we adjusted the chloride-hydroxyl radial distribution functions against ab initio data and then extended the use of the Tang-Toennis damping function for the chlorides' induced dipoles. These adjustments correct the problems without losing the robustness of the CL&Pol force field. The results were also compared with the nonpolarizable version, the CL&P force field. We expect that the corrections will facilitate reliable use of the CL&Pol force field for other types of DESs.
Assuntos
Cloretos , Líquidos Iônicos , Cloretos/química , Solventes Eutéticos Profundos , Líquidos Iônicos/química , Simulação de Dinâmica Molecular , Solventes/químicaRESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.