Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 02 May 2019
Sec. Green and Sustainable Chemistry
This article is part of the Research Topic Ionic Liquids: Properties and Applications View all 30 articles

Microscopic Structural and Dynamic Features in Triphilic Room Temperature Ionic Liquids

\nFabrizio Lo CelsoFabrizio Lo Celso1Giovanni B. AppetecchiGiovanni B. Appetecchi2Elisabetta SimonettiElisabetta Simonetti2Man ZhaoMan Zhao3Edward W. Castner Jr.Edward W. Castner Jr.3Uwe KeiderlingUwe Keiderling4Lorenzo GontraniLorenzo Gontrani5Alessandro Triolo*Alessandro Triolo6*Olga Russina*Olga Russina5*
  • 1Dipartimento di Fisica e Chimica, Università di Palermo, Palermo, Italy
  • 2ENEA, Laboratory SSPT-PROMAS-MATPRO, Rome, Italy
  • 3Department of Chemistry and Chemical Biology, Rutgers University, The State University of New Jersey, Newark, NJ, United States
  • 4Soft Matter and Functional Materials, Helmholtz-Zentrum für Materialien und Energie GmbH, Berlin, Germany
  • 5Department of Chemistry, University of Rome Sapienza, Rome, Italy
  • 6Laboratorio Liquidi Ionici, Istituto Struttura della Materia, Consiglio Nazionale delle Ricerche (ISM-CNR), Rome, Italy

Here we report a thorough investigation of the microscopic and mesoscopic structural organization in a series of triphilic fluorinated room temperature ionic liquids, namely [1-alkyl,3-methylimidazolium][(trifluoromethanesulfonyl)(nonafluorobutylsulfonyl)imide], with alkyl=ethyl, butyl, octyl ([Cnmim][IM14], n = 2, 4, 8), based on the synergic exploitation of X-ray and Neutron Scattering and Molecular Dynamics simulations. This study reveals the strong complementarity between X-ray/neutron scattering in detecting the complex segregated morphology in these systems at mesoscopic spatial scales. The use of MD simulations delivering a very good agreement with experimental data allows us to gain a robust understanding of the segregated morphology. The structural scenario is completed with determination of dynamic properties accessing the diffusive behavior and a relaxation map is provided for [C2mim][IM14] and [C8mim][IM14], highlighting their natures as fragile glass formers.

Introduction

Ionic liquids (ILs) constitute an interesting class of compounds, composed solely of ionic species and with a melting point below 100°C. Although some debate exists on their environmental compatibility, they are attracting a great deal of attention as they are considered potential replacements for noxious solvents and flammable electrolytes. Their fascinating chemical and physical properties can be widely modulated upon slight chemical modification of constituent ionic species. In the framework of their potential application in a variety of fields, the possibility of introducing a fluorous moiety into either the cation or the anion (e.g., by replacing an alkyl with an equivalent length perfluoroalkyl chain) is quite attractive, as it allows accessing further modulation of properties and performances. In these cases, the already appealing properties of conventional ILs, including non-flammability, negligible vapor pressure, high thermal, chemical and electrochemical stability, are maintained and can be further fine-tuned by introduction of fluorous moieties leading to high hydrophobicity, enhanced surface activity, interesting gas uptake properties etc. Fluorinated Ionic Liquids (FILs) show properties arising from both their fluorous and ionic natures, thus introducing interesting variations into the more conventional IL landscape. As both fluorinated compounds and ILs are considered neoteric solvents with the potential toimpact sustainable processes and green chemistry, it can be envisaged that their merging into the FILs technology will soon be of high impact (Merrigan et al., 2000; van den Broeke et al., 2002; Kim et al., 2004; Xue and Shreeve, 2005; Xue et al., 2006; Almantariotis et al., 2010, 2017; Smith et al., 2010; Kunze et al., 2011; Yoshida and Saito, 2011; Jeremias et al., 2013; Pereiro et al., 2013a, 2018; Weber et al., 2013; Hollóczki et al., 2015; Suarez et al., 2015; Wu et al., 2016; Rauber et al., 2017; Bastos et al., 2018).

Among other peculiar features in FILs, their mesoscopic structural organization is quite interesting. It is well-known that alkyl chain bearing ILs are characterized by a distinct degree of nm-scale order, associated to the spatial segregation of apolar domains (viz. alkyl chains) with respect to a polar matrix (viz. anion and cation heads) (Urahata and Ribeiro, 2004; Wang and Voth, 2005, 2006; Canongia Lopes and Padua, 2006; Canongia Lopes et al., 2006; Pádua et al., 2007; Triolo et al., 2007; Almantariotis et al., 2010, 2017; Smith et al., 2010). Such a structural organization can be experimentally detected by means of X-ray and/or neutron scattering techniques, as it is fingerprinted by a distinct low momentum transfer (Q) peak (that is centered at a position reflecting the characteristic size of the segregated domain), as it was both experimentally and computationally found (Triolo et al., 2007, 2008, 2009; Atkin and Warr, 2008; Russina et al., 2009, 2012, 2017a,b; Hayes et al., 2011; Zheng et al., 2011; Li et al., 2012; Macchiagodena et al., 2012; Russina and Triolo, 2012, 2017; Song et al., 2012; Rocha et al., 2013). When introducing a fluorous moiety, as it is incompatible with both polar moieties and alkyl tails, it will tend to segregate into a third further kind of domain, where fluorous tails alone tend to cluster. This will lead to a characteristic mesoscopically organized morphology where three different, macroscopically incompatible, domains are forced to co-exist at the nm spatial scale: this led some of us to term such a kind of compounds as triphilic IL (Russina et al., 2013; Hollóczki et al., 2015). Such a scenario (Shen et al., 2012; Greaves et al., 2013; Pereiro et al., 2013b; Hettige et al., 2014; Brehm et al., 2015; Vieira et al., 2015; Ferreira et al., 2017; Shimizu et al., 2017) opens the way to simultaneous solvation of completely incompatible components into the homogeneous FIL, as each compound (e.g., a salt, an oil, and a fluorous compound) can be simultaneously dissolved into the same mesoscopically organized medium (Weiss et al., 2017; Tsurumaki and Ohno, 2018). In this contribution, we report on the structural properties of a series of 1-alkyl-3-methyimidazolium ([Cnmim], where n indicates the number of carbon atoms along the side alkyl chain, with n = 2, 4, and 8) (trifluoromethanesulfonyl)(nonafluorobutylsulfonyl) imide ([IM14], 1 and 4 represent the number of carbon atoms on the two sides of the imide anion) [see e.g., (Johansson et al., 2010; Montanino et al., 2012; Jeremias et al., 2013; Russina et al., 2013, 2017a; Castiglione et al., 2014; Paolone et al., 2018)], at ambient conditions. Scheme 1 represents the chemical structure of the probed FILs.

SCHEME 1
www.frontiersin.org

Scheme 1. Description of chemical nature of the investigated ionic liquids. On the (Left), the anion's structure is shown, while the (Right) part of the Scheme shows the different cations [1-alkyl-3-methylimidazolium, with alkyl = ethyl (Top), butyl (Middle), and octyl (Bottom)].

Experimental Details

Samples Synthesis

The imidazolium-based, (trifluoromethanesulfonyl)(nonafluorobutylsulfonyl)imide ionic liquids were synthesized through a procedure developed at ENEA and described in detail elsewhere (Jeremias et al., 2013; Montanino et al., 2013). The chemicals, i.e., 1-methylimidazole (99 wt.%), bromoethane (99 wt.%), 1-bromobutane (99 wt.%), 1-bromooctane (99 wt.%), were purchased by Sigma-Aldrich and used as received. The sorbent materials, i.e., activated carbon (Darco-G60, Sigma-Aldrich) and alumina (acidic, Brockmann I, Sigma-Aldrich), were previously rinsed in deionized water according to a route reported elsewhere (Palumbo et al., 2017). Acidic (trifluoromethanesulfonyl)(nonafluorobutylsulfonyl)imide, HIM14 (59 wt.% solution in water), was purchased by 3M and used as received. Deionized water, used as the processing solvent, was obtained with a Millipore ion-exchange resin deionizer.

The imidazolium precursor was synthesized by reacting the 1-methylimidazole, previously dissolved in deionized water, with the appropriate amount of bromoethane (or 1-bromobutane or 1-bromooctane). A concentrated aqueous solution of imidazolium bromide was obtained, which was then purified through activated carbon and acidic alumina. The liquid fraction (clear and colorless) was separated from the solid one (sorbent materials) by vacuum filtering and separately collected. Successively, the solid phase was rinsed with deionized water to recover the precursor trapped through the purifying materials, obtaining a further aqueous solution (clear and colorless) which was unified to the previous one.

The imidazolium IM14 ionic liquids were obtained by reacting the precursor water solutions and the suitable amount of HIM14. The reaction (anion exchange) led to fast formation of hydrophobic ionic liquid and hydrophilic HBr (in water). After 1 h the mixing was interrupted and phase separation was allowed for a few minutes. The upper phase was mostly composed of water, HBr and HIM14 excess whereas the lower one was constituted of the ionic liquid with traces of water and HBr. After removal of the aqueous phase, the imidazolium IM14 ionic liquid was rinsed several times with deionized water to remove the HBr and HIM14 impurities. The ionic liquids were then housed in a vacuum rotary evaporator at 90°C to remove most of the water and, finally, vacuum dried at 120°C. The ionic liquid materials were stored in sealed glass tubes within a dry room.

The Li+ and Br content was checked by atomic absorption analysis (SpetcrAA mod. 220 atomic absorption Spectrometer) and X-ray fluorescence spectrometry (Shimadzu energy-dispersion EDX-720 spectrometer), respectively. Less than 2 ppm of lithium and bromide were detected in the ionic liquid samples.

X-Ray Scattering

Small-Medium Angle X-ray scattering data were collected at the 11-ID-C beamline at the Argonne Photon Source (Ren, 2012). A monochromatic beam at a wavelength of 0.11165 A was used, with a beam size of 0.7 × 0.7 mm2. Samples were contained in a 2 mm diameter quartz capillary purchased from WJM-Berlin. Thirty consecutive data acquisitions were collected for a total acquisition time of 200 sec. Before merging the different data sets, care was taken that no radiation damage occurred to samples. Raw 2D data sets were integrated and transformed into intensities using the FIT2D software (Hammersley et al., 1996). The scattering structure factor S(Q) was calculated using the PDFgetX2 software (Qiu et al., 2004), correcting for background contribution, absorption, Compton scattering, and multiple scattering.

Additionally Small-Medium-Wide Angle X-ray scattering data were collected at beamline BL04B2 SPring-8 (Japan Synchrotron Radiation Research Institute, JASRI, Japan) (Kohara et al., 2001; Ohmura et al., 2009; Aoun et al., 2010). Measurements were conducted at ambient conditions. A monochromatized X-ray (61.6 keV) was obtained using a Si(220) monochromator. The observed scattering intensity of the X-ray was corrected for absorption, polarization, and incoherent scatterings (Cromer, 1969) to obtain coherent scattering intensities.

Small Angle Neutron Scattering

Small Angle Neutron Scattering experiments (SANS) measurements on [Cnmim][IM14] (n = 2, 4, 8) were performed at the small-angle scattering instrument V4 which is placed in the cold neutron guide of Helmholtz-Zentrum Berlin (HZB). The magnitude of the scattering vector is defined as Q = (4π/λ)(sinθ) with λ being the wavelength and 2θ the scattering angle. The measured neutron flux of the V4 is ~106 cm−2.s−1 for the wavelength used, λ = 4.5 Å (Keiderling and Wiedenmann, 1995). The scattering data were obtained at one sample detector distance of 1.0 m, which yields a total momentum transfer range of 0.5 nm−1 < Q < 8 nm−1. For further information regarding the V4 instrument and its resolution the reader is referred to Keiderling and Wiedenmann (1995)and Gilles et al. (2000). The sample was placed into a circular quartz cuvette with inner spacing of 1 mm and placed in the beam for measurement. A Cd aperture of 13 mm was used for the scattering measurements. The 2D scattering data were reduced to a scattering curve [dΣ/dΩ(Q) vs. Q; hereinafter indicated as I(Q) vs. Q] by means of the BerSANS software. The raw data is then corrected for transmission, the quartz cell background scattering subtracted and converted to absolute units taking into account the scattering from water (Keiderling, 2002).

Wide Angle Neutron Scattering

Neutron diffraction measurements on [C2mim][IM14] were performed on the SANDALS diffractometer at ISIS. The neutron wavelength range is 0.05–4.95 Å. Data were collected over the momentum transfer (Q) range between ca. 0.3 and 50 Å−1. The sample was contained in chemically inert, null scattering Ti0.68Zr0.32 flat cans (with size 3.5 × 3.5 × 1 mm3) sealed with Teflon O-rings. Measurements were also collected on the empty cell, so to properly subtract empty cell contribution. Furthermore, a measurement on a vanadium standard sample was also collected for data normalization purposes. Diffraction experiments were conducted at 298 K under vacuum. The diffraction pattern was measured for approximately 6 h at conditions of fully operating source. Data analysis was carried out using the GUDRUN software available at the facility that allows the application of (a) normalization to the incident flux, absorption and multiple scattering corrections, empty can subtraction, and normalization to absolute units by dividing the measured differential cross section by the scattering of a vanadium standard and (b) corrections for single atom scattering and hydrogen inelasticity effects.

Diffusivity Measurements From PG-SE NMR

Self-diffusion coefficients for the anions and cations were measured using pulsed gradient spin-echo (PG-SE) NMR methods. Anionic diffusivities were measured from 19F NMR signals, and cationic diffusivities were obtained from 1H signals. The Diffusion-ordered bipolar pulse pair stimulated echo (DBPPSTE) pulse sequence was used (Wu et al., 1995). NMR samples were prepared from ionic liquids dried for 48 h on a Schlenk vacuum line, and prepared in an argon glovebox with water and oxygen levels below 0.1 and 0.4 ppm, respectively. Samples were prepared in 3 mm O.D. NMR tubes, capped in the glovebox, and flame-sealed immediately after removal from the glovebox. A Doty Scientific diffusion probe was installed on a 400 MHz Varian DirectDrive spectrometer for the PG-SE NMR measurements. PG-SE experiments had temperature uncertainties of ±1 K.

Additional details of our NMR diffusion protocols have been described previously (Liang et al., 2015; Wu et al., 2015, 2018).

Viscosity Measurements

Temperature-dependent viscosity measurements of the ionic liquid samples were made using a Cambridge Viscosity ViscoLab 4100 instrument, as described previously (Funston et al., 2007; Mariani et al., 2017). Temperatures for the viscosity measurements were controlled to ±0.1 K by water flow from a Lauda Brinkmann RMT-6 recirculating chiller/heater.

Dielectric Spectroscopy Measurements

Dielectric spectroscopy data were collected between 160 and 280 K in the range 2 10−2 – 3 106 Hz, applying the gain-phase analysis technique with a Solartron SI-1260 analyser and a Novo-Control BDS 4000 spectrometer. The sample was sandwiched between two gold plated flat electrodes, with diameter of 20 mm and with a sample thickness of 50 μm. The temperature was controlled by a Quatro cryosystem.

Computational Details

Molecular dynamic simulations were performed using the GROMACS 5.1.1 package (Van Der Spoel et al., 2005; Hess et al., 2008). Interactions were described using an all-atoms potential (Lopes and Padua, 2004; Shimizu et al., 2010). The simulations for [Cnmim][IM14] ILs were performed using a cubic box of 900, 800, and 700 ion pairs for n = 2, 4, and 8, respectively; periodic boundary conditions were applied. Force field parameter files and initial configuration were created by DLPGEN software (Bernardes and Joseph, 2015); initial density was fixed 10% higher than the experimental one. A procedure of energy minimization has been performed starting from the initial configuration prior to the first NPT equilibration. The steepest descent algorithm implemented in Gromacs has been used and convergence has been obtained when every force on each atom would not exceed 1,000 kJ/mol/nm. Cutoff for the short range electrostatics and Van der Waals interaction were set to 1.0 nm. The equilibration procedure was done in several steps, starting from a series of NPT simulation at high temperatures and scaled partial charges, followed by lowering progressively the temperature and increasing the charges to their final value at 298,15 K and 1 bar, after a 15 ns run. This procedure was repeated two further times until an equilibrated system was obtained. After the equilibration phase, the system was run for a total of 50 ns for a production run, and then the trajectory of the last 5 ns was saved at a frequency of 1 ps, for calculation of the structural properties. The simulations were always checked vs. the experimental density and the energy profile. During the production runs for the temperature coupling, we used a velocity rescaling thermostat (Bussi et al., 2007) (with a time coupling constant of 0.1 ps), while for the pressure coupling, we used a Parrinello–Rahman barostat (Parrinello and Rahman, 1981) (1 ps for the relaxation constant). The Leap-Frog algorithm with a 1 fs time step was used for integrating the equations of motion. Cut-offs for the Lennard- Jones and real space part of the Coulombic interactions were set to 15 Å. For the electrostatic interactions, the Particle Mesh Ewald (PME) summation method (Darden et al., 1993; Essmann et al., 1995) was used, with an interpolation order of 6 and 0.08 nm of FFT grid spacing. Selected graphs were done using Matplotlib (Hunter, 2007). Weighted and partial structure factors were computed by using in-house developed software, accordingly to text book formulas as highlighted in Margulis's work (Kashyap et al., 2012), while selected pair correlation function, angular distribution function were obtained by TRAVIS (Brehm and Kirchner, 2011; Hollóczki et al., 2015). Ion aggregation analysis was performed using the AGGREGATES software (Bernardes, 2017).

Results and Discussion

The chemical structure of [Cnmim][IM14] samples with n = 2, 4, 8 that have been studied in this report is shown in Scheme 1. These samples differ only on the length of the side alkyl chain in the cation, the anion remaining unaltered. [C2mim][IM14] has been previously reported to be an homogeneous liquid at ambient conditions (Appetecchi et al., 2011). Analogously, [C4mim][IM14] and [C8mim][IM14] are liquid at ambient conditions and their synthesis and properties are here reported for the first time.

Viscosity

The values for the viscosities for [C2mim][IM14] measured between 1.5 and 90°C spanned the range from 520 down to 12 cP. For [C8mim][IM14], the viscosities ranged from 1,011 to 28 cP over the temperature range from 3.7 to 70°C. As with all of our previous recorded viscosities for ionic liquids the data are best fit to the well-known Vogel-Fulcher-Tammann (VFT) equation in logarithmic form, given by: ln (η(T), cP) = ln (ηVFT)+B/(T−To); where ηVFT, B, and To are fitting parameters, representing infinite temperature viscosity, an equivalent activation energy and a characteristic temperature, respectively.

The primary significance of these fits is that they provide an excellent representation of the data, permitting us to calculate an extrapolated viscosity for any arbitrary temperature in order to match the viscosity/temperature point for comparison with the diffusivities measured at different temperatures via PG-SE NMR. While the viscosity data does not fit as well to the Arrhenius equation as it does to the VFT equation, the Arrhenius fits are adequate for a qualitative discussion. The viscosity activation energies were found to be 37.2 kJ/mol for [C2mim][IM14] and 44.1 kJ/mol for [C8mim][IM14].

The fits of the viscosity data to the VFT and the Arrhenius equations are given in Figures S1, S2.

Anion and Cation Diffusivities From PG-SE NMR

PG-SE NMR measurements were used to obtain the anionic and cationic diffusivities over the temperature range from 5 to 75°C for both [C2mim][IM14] and [C8mim][IM14]. Cationic diffusivities were observed to always larger than anionic diffusivities for both [C2mim][IM14] and [C8mim][IM14], though the difference is more pronounced for the former case. For the [C2mim][IM14] IL, the anion is significantly larger than the cation, while the opposite is true for [C8mim][IM14]. It is worth noting that the mass of the [IM14] anion (430 g/mol) is roughly four times that of the [C2mim] cation (Fw = 111 g/mol) and twice that of the [C8mim] cation (Fw = 195 g/mol). For [C2mim][IM14], the diffusivities range from 3.4 × 10−12 m2s−1 for the [IM14] anion at 5°C up to 1.1 × 10−10 m2s−1. For comparison, we note that the diffusivity of neat water is 1.26 × 10−9 m2s−1 at 4°C.

As is the case for viscosities, the ionic self-diffusivities can also be fit by both Arrhenius and VFT laws, with the latter always providing somewhat superior fit statistics, albeit at the cost of one additional fit parameter. Considering the activation energies from the Arrhenius fits, one finds that the Ea values for [C2mim][IM14] are slightly less than that for the viscosity: 32.5 and 35.4 kJ/mol for the cation and anion, respectively, as compared with 37.2 kJ/mol for the viscosity. The Ea values for [C8mim][IM14] self-diffusion coefficients are larger than [C2mim][IM14], but again, less than the Ea value for the viscosity: Ea(cation) = 35.9 kJ/mol and Ea(anion) = 37.7 kJ/mol, vs. 44.1 kJ/mol (viscosity).

Harris showed that hydrodynamic scaling is not exact for any molecular fluids considered, including ionic liquids (Harris, 2009). Specifically, the Stokes-Einstein prediction that the diffusivities should be proportional to the ratio of temperature over viscosity does not hold exactly. Rather, a scaled or fractional Stokes-Einstein law does hold, where D = (T/η)α, where α is an exponent between 0.9 and 1.0. By using the VFT fit parameters, we can calculate precise viscosities in order to make graphs of the ionic diffusivities vs. the ratio of T/η These plots of D vs. T/η are shown in Figure S3.

Structure

In Figure 1, the X-ray diffraction patterns from the three samples are shown: in the reported Q range they are characterized by the presence of three peaks (QI−III, in order of increasing Q position), that fingerprint the existence of different kinds of correlations (QI is better discernible in the case of [C8mim][IM14]). In particular, in order of decreasing Q position, these peaks have been shown to reflect the existence of: (a) adjacency, (b) charge, and (c) polar correlations (Annapureddy et al., 2010; Araque et al., 2015). It appears that, similarly to other cases, the peak that is most affected by increase in chain length is QI that reflects the establishment of polar-apolar correlations (Triolo et al., 2007, 2009; Russina and Triolo, 2012; Russina et al., 2012): this peak has negligible amplitude for short chain length and it grows in amplitude and shifts to lower Q value upon increasing n. Such a behavior has been reported for a multitude of ILs and reflects the progressive increase in size of the apolar domain embedded into the polar matrix that can be qualitatively quantified as DI~2π/QI. In the inset of Figure 1, the characteristic domain sizes obtained for the case of [Cnmim][IM11] (where [IM11] is the commonly used anion bistriflamide, [Tf2N]), reported by some of us in the past (Russina et al., 2009), are shown as a function of n. It can be observed that the corresponding quantity for the case of [C8mim][IM14] (for which QI peak position can be easily determined) is quite similar to the value reported for [C8mim][IM11], thus suggesting that X-ray scattering is a good reporter of the alkyl domain heterogeneities, whose size is only marginally affected by the different anion's nature. On the other hand, Figure 2 reports Small Angle Neutron Scattering (SANS) data from the same [Cnmim][IM14] samples at ambient conditions, over the Q range up to 8 nm−1 (thus covering only the Q range where QI typically occurs, see Figure 1 for comparison). It clearly appears that while [C8mim][IM14] is characterized by a peak centered at ca. 3 nm−1 and [C4mim][IM14] shows a broad peaks superposition between 1 and 6 nm−1, [C2mim][IM14] is instead characterized by a distinct peak centered at ca. 4.5 nm−1. The latter peak has no analogous counterpart in the X-ray diffraction pattern shown in Figure 1. Moreover, the inset of Figure 2 shows a comparison between the SANS data from [C2mim][IM14] and [C2mim][IM11]: it is clear that while [C2mim][IM11] is essentially featureless in the probed Q range, on the other hand [C2mim][IM14] shows a distinct peak. The featureless pattern from [C2mim][IM11] is rationalized considering that the side alkyl chain in the cation (an ethyl moiety) is too short to deliver a clustering of chains segregating from the polar moieties, accordingly no polar-apolar separation occurs in this system. On the other hand the clear existence of low Q peak in the case of [C2mim][IM14] fingerprints the existence of mesoscopic structural heterogeneities that, due to contrast reasons, neutrons can detect. As recently reported in a series of papers (Shen et al., 2012; Russina et al., 2013; Hettige et al., 2014; Hollóczki et al., 2015; Lo Celso et al., 2017, 2018a,b), neutron scattering (and occasionally, X-ray scattering, vide infra) succeed in detecting the occurrence of fluorous tails clusters embedded into the charged matrix: the results presented in Figures 1, 2 confirm this observation, indicating that the perfluorobutyl chains in the [IM14] anions tend to associate into small fluorous domains and in the rest of the manuscript we will further explore this structural feature in the probed FILs.

FIGURE 1
www.frontiersin.org

Figure 1. Experimental X-ray diffraction patterns from [Cnmim][IM14] (n = 2, 4, 8) at ambient conditions. In the inset, the characteristic sizes of structural heterogeneities, obtained from 2π/QI, where QI is the low Q peak in the diffraction pattern, are shown for [Cnmim][IM11] (with n = 6–10) [from Russina et al. (2009)] (black symbols), together with the corresponding value from [C8mim][IM14] (present work, red symbol).

FIGURE 2
www.frontiersin.org

Figure 2. Experimental small angle neutron scattering data from [Cnmim][IM14] (n = 2, 4, 8) at ambient conditions. In the inset a comparison between SANS data from [C2mim][IM14] and [C2mim][IM11] is shown.

In order to better address the issue of microscopic and mesoscopic structural organization in [Cnmim][IM14] FILs, we undertook a series of detailed Molecular Dynamics simulations on these compounds, aiming to extract information at atomistic level. In the experimental part, we described the procedure and other details used for these simulations. In the case of [C2mim][IM14], the computed neutron and X-ray diffraction patterns have been directly compared with experimentally derived data sets. Figures 3A,B show such a comparison, while Figures 3C,D describe the comparison between the experimental SANS data and the corresponding computed quantities for the case of [Cnmim][IM14], with n = 4 and 8: we consider the overall agreement between experimental and computational diffraction pattern very satisfactory. As a matter of fact, the simulations (along with the chosen interatomic potentials) succeed in accounting for all the structural features fingerprinted in the experimental diffraction patterns and this makes us confident that we can reliably interrogate the MD simulations to extract structural details at atomistic level.

FIGURE 3
www.frontiersin.org

Figure 3. Comparison between (A) X-ray and (B) neutron diffraction data and the output from present MD simulations for the case of [C2mim][IM14] at ambient conditions. Also, comparison between experimental (black line) and un-normalized MD-computed (red line) SANS data from [Cnmim][IM14], with (C) n = 4 and (D) n = 8, respectively, at ambient conditions.

Figure 4 show representative snapshots of simulation boxes for the three different samples [Cnmim][IM14] (n = 2, 4, 8), where polar moieties (imidazolium ring and SO2-N-SO2 moieties, belonging to cation and anion, respectively), alkyl tails and fluorous moieties are identified. A distinct mutual spatial segregation can be observed for these moieties, depending on the alkyl side chain length and we will better clarify this issue later on.

FIGURE 4
www.frontiersin.org

Figure 4. Representative snapshots derived from MD simulations from [Cnmim][IM14], with n = 2, 4, 8 (left, center, and right, respectively), where color distinction of charged (red), alkyl (blue), and fluorous (yellow) moieties has been applied.

The characteristic alternation of oppositely charged layers around a chosen reference ion that is encountered in ILs, can be observed also in the present choice of compounds. Figure 5 shows the representative case of [C2mim][IM14] (different values of n lead to qualitatively similar results), where pair distribution functions (pdf) for: (a) center of mass of the imidazolium ring (ir) – ir; (b) ir – Nanion and (c) Nanion – Nanion, Nanion being the nitrogen atom in the anion, are reported. It appears that cation-anion correlations are the shorter ones, occurring at ca. 4 Å, in agreement with the nature of coulombic interactions. Spatial correlations between similarly charged ions occur at a larger distance (ca. 10 Å) and such correlations are out of phase with opposite charges correlations (their maxima fall where opposite charge pdf has minima).

FIGURE 5
www.frontiersin.org

Figure 5. (Left) MD-derived pair distribution functions for ir-ir (red line), ir-Nanion (blue line), and Nanion-Nanion (green line), where ir is the imidazolium ring center of mass and Nanion is the anion's nitrogen atom, from MD simulation of [C2mim][IM14] at ambient conditions. (Right) combined distribution function with abscissa the ir-ir distance and ordinate the angle between vector ir-ir and the vector normal to the imidazolium ring, from MD simulation of [C2mim][IM14] at ambient conditions. The plot explains the two maxima features occurring in the ir-ir pdf between 7 and 11 Å.

We notice that the cation-cation correlation (ir-ir) is characterized by a short distance correlation (red line in Figure 5 left) with two maxima at ca. 7 and 10 Å. In order to better rationalize this correlation, the combined distribution function (cdf), containing in abscissa the distance between neighbor imidazolium rings and in ordinate the angle (°) formed by vectors perpendicular to the imidazolium rings, has been evaluated and shown in Figure 5 right. We can then envisage that the observed splitting of the first solvation shell of imidazolium rings around a reference ring is the consequence of different packing of neighbor rings: the hottest lobes in Figure 5 right are found at ca. 8.2 Å for parallel or anti-parallel neighbor rings; however, non-negligible population is detected for perpendicular rings at distances of order of 7 and 10 Å.

Following an analysis approach introduced for the case of ILs by Margulis et al. (Annapureddy et al., 2010; Kashyap et al., 2012, 2013; Hettige et al., 2014; Araque et al., 2015), here we decompose the computed X-ray/neutron diffraction patterns in their components arising from structural correlations between polar and apolar moieties in the samples. While several examples exist of such a decomposition for the case of conventional ILs, where charged moieties and long alkyl tails represent the polar and apolar moieties, respectively, in the present case, similarly to recent analogous situations (Hettige et al., 2014), two different kinds of apolar moieties have been identified: namely, the alkyl and the fluorous tails. Such an approach stems from the quantitative interpretation of the visual situation depicted in Figure 4, where polar moieties as well as apolar alkyl and fluorous moieties tend to mutually exclude each other, thus introducing a complex triphilic structural scenario, with three different classes of nm-scale domains, whose size and mutual distribution determines the experimentally determined diffraction patterns. In the simple case of polar (P) and one (e.g., the alkyl moieties) apolar (A) domains, Margulis' group has shown that a characteristic decomposition can be achieved for the low Q peak, in terms of P-P, A-A, and A-P contributions, the two former terms manifesting into a positive amplitude peak and the latter into a negative amplitude one at the experimental peak position (Annapureddy et al., 2010; Kashyap et al., 2012, 2013; Hettige et al., 2014; Araque et al., 2015). In the case of two different apolar (Aa and Af, where the subscript a and f stand for alkyl and fluorous, respectively) domains coexisting with the polar (P) one, all the different self and distinct contributions should be considered to contribute to the overall experimental pattern (Hettige et al., 2014). Figure 6 report such decompositions of computed X-ray and neutron diffraction patterns; in the figures, the self-terms P-P (indicated therein as P), Aa-Aa (indicated therein as Aa), and Af-Af (indicated therein as Af) are shown together with the summation of the different distinct terms P-Aa, P-Af, Aa-Af (the combination of the latter terms being indicated as “cross”): the combination of these contributions leads to the total diffraction pattern. Overall, due to contrast reasons (contrast acts as a scale factor for each contribution and depends on either electron densities or scattering length densities for X-ray and neutron, respectively), the contribution of P/Aa alternation appears evident only in the case of [C8mim][IM14], where a positive amplitude peak is found for Aa-Aa and a negative amplitude one is found for P-Aa, both centered at a position QP/Aa ~3.5 nm−1, fingerprinting the existence of the P/Aa alternation and its characteristic size of the order of 2π/QP/Aa~18 Å. The decomposition of both X-ray and neutron scattering allows detecting in a very clear way the structural alternation P/Af, between polar and fluorous domains. For example in Figure 7, the contributions P (Polar-Polar), Af (Af-Af), and cross (P-Af) to X-ray and neutron scattering are shown for [C4mim][IM14]. It appears that both in the case of X-ray and neutron scattering, peaks and anti-peaks related to both P-P and Af-Af and P-Af alternations, exist, respectively and fall at a position QP/Af ~4 nm−1. It is noteworthy that while in the case of X-ray scattering (Figure 7 left), the combination of these three contributions essentially vanishes, thus leading to no net contribution of the P/Af alternation to the total diffraction pattern, on the other hand, in the case of neutron scattering, the combination of the three terms leads to an overall non-zero contribution, that results into the low Q peak observed in the experimental neutron scattering data set (see Figures 1, 2).

FIGURE 6
www.frontiersin.org

Figure 6. (Top) X-ray MD-computed diffraction patterns for [Cnmim][IM14], with n = 2, 4, 8 (left, center, and right, respectively) and their decomposition into polar-polar (P), alkyl-alkyl (Aa), and fluorous-fluorous (Af-Af), total and cross terms; (Bottom) neutron MD-computed diffraction patterns for [Cnmim][IM14], with n = 2, 4, 8 (left, center, and right, respectively) and their decomposition into polar-polar (P), alkyl-alkyl (Aa), and fluorous-fluorous (Af), total and cross terms.

FIGURE 7
www.frontiersin.org

Figure 7. X-ray (Left) and neutron (Right) MD-computed diffraction patterns decomposition into polar-polar (P), fluorous-fluorous (Af), and polar-fluorous (cross) terms, for [C4mim][IM14] at ambient conditions.

An estimation of the degree of clustering of both alkyl and perfluoroalkyl tails can be also be obtained by monitoring the sizes of the clusters of such moieties. Using an approach developed by Bernardes et al. (2011, 2014), Shimizu et al. (2014) we analyzed our MD trajectories using the AGGREGATES software made available by Bernardes (2017) We considered CH2/CH3 moieties as belonging to an alkyl domain when neighbor carbon atoms were 5 Å apart, while we considered CF2/CF3 groups belonging to a fluorous domain when neighbor carbon atoms were 6 Å apart (Shimizu et al., 2017). Under these conditions, we monitored the tendency of both alkyl and fluorous chains to segregate into clusters. Figure 8 shows the trends observed for the dependence of (alkyl or fluorous) aggregates distribution sizes on the alkyl side chain length. It can be noticed that alkyl tails tend to form limited size clusters for the case of [Cnmim][IM14] with n = 2 and 4 (clusters containing up to ~10 and ~30 methyl groups, respectively), while in [C8mim][IM14], octyl chains tend to percolate, as the size of aggregates reaches ~700 (representing the total number of anions present in the simulation box). Accordingly a gradual transition from finite size to percolating clusters is observed. In the case of anion's perfluorobutyl chains, [Cnmim][IM14] with n = 2 and 4 show percolating nature of the chains clusters, while, in the case of [C8mim][IM14] a finite size distribution with maximum size of ~25 members is found, indicating also in this case the existence of a gradual transition from finite size to percolating clusters, though, with opposite trend with respect to n. When alkyl tails form finite size clusters, perfluorobutyl chain succeed in percolating the simulating box and, vice versa, when the latter chains form only finite clusters, it will be the alkyl chains that will percolate across the simulation box. Recently [Cnmim]-based FILs with a fluorous anion have been simulated by Canongia- Lopes et al. (Pereiro et al., 2013b; Vieira et al., 2015; Ferreira et al., 2017; Shimizu et al., 2017). They found a behavior that is nice agreement with the present observations: while monitoring the aggregates detected in [Cnmim][C4F9SO3], they also found that by increasing the side alkyl chain length of the cation, the alkyl domains tend to grow in size until they eventually percolate and this process is inverse in trend to the one observed for the fluorous domains sizes that instead decrease from the percolating regime to finite size clusters. Herein, on the basis of experimental results, supported by MD simulations, we further (Lo Celso et al., 2017, 2018b) validate this proposal.

FIGURE 8
www.frontiersin.org

Figure 8. Normalized distribution functions for the alkyl (Left) and fluorous (Right) domain sizes from MD simulations of [Cnmim][IM14] (n = 2, 4, 8, black, red, and blue, respectively) at ambient conditions. The highest value of the abscissa corresponds to the number of ion pairs contained in the simulation boxes.

Main Relaxation Process

We aim now at providing a preliminary description of relaxation phenomena occurring in selected [Cnmim][IM14] ILs, on the basis of experimental data. The nature of relaxation processes occurring in both [C2mim][IM14] and [C8mim][IM14] has been explored using a range of techniques. Figure 9 reports dielectric spectroscopy data collected for [C2mim][IM14] in the form of isochronal temperature dependence of the imaginary part of the electrical modulus, M′′. This formalism is commonly used when dealing with conducting samples, as it allows suppressing the contribution from conductivity and detecting the main α process. The latter process is generally associated to diffusive dynamics in glass formers and has equivalent signatures when monitored with other complementary techniques, such as viscosity or conductivity. In this case it is customary to account for the overall temperature dependence of the process by using the time/temperature superposition (tTS) that would provide a complete picture of the relaxation map of the chosen material. In Figure 10, we then report the relaxation map built up by the merging of different experimental data sets obtained for the case of [C2mim][IM14] as well as [C8mim][IM14]. While a more extended interpretation of these experimental results will be presented elsewhere, here we compare data arising from the above mentioned dielectric spectroscopy experiment, from previously published (Appetecchi et al., 2011) and new (this work) viscosity experiments. Using the tTS approach, data referring to the viscosities of the two ILs are combined with the characteristic times extracted from the M′′ data sets (once properly vertically scaled) and are jointly modeled in terms of a Vogel-Fulcher-Tamman (VFT) temperature dependence (Fulcher, 1925; Tammann and Hesse, 1926): log x = log xo+B/(T−To), where x refers to a given quantity (e.g., viscosity, relaxation time etc.), xo, B, and To are fitting parameters, accounting for the temperature dependence of the observed property. The VFT formalism nicely accounts for the divergence of diffusion-related quantities (viscosity and characteristic time for the main α-process in the dielectric spectrum), when approaching the glass transition. The dielectric data sets indicate that characteristic times of 100 s are reached at 187 and 190 K for [C2mim][IM14] and [C8mim][IM14], respectively. These are then considered to be the dielectric glass transitions (Tg,diel) for these compounds. On the basis of the VFT fitting parameters and the obtained Tg,diel's, one can derive the corresponding fragility indexes (Angell, 1985, 1991; Böhmer and Angell, 1992; Böhmer et al., 1993), after: m = B/ln(10) Tg/(Tg−To)2: we obtain then values of 84 and 76 for [C2mim][IM14]'s and [C8mim][IM14]'s fragilities, respectively. These are relatively large values for the fragility index and reflect some degree of fragility in this class of compounds [the larger m, the more fragile is the compound, according to the definition proposed by Angell (1985, 1991), Böhmer and Angell (1992), Böhmer et al. (1993). The trend toward decreasing fragility (m) upon increasing alkyl chain length in ILs has been shown elsewhere, e.g., Sippel et al. (2015), Leys et al. (2008), Tao et al. (2017), where m values for a series of [Cnmim]-based ILs with [BF4] and [PF6] anions have been reported and its relevance for energy applications highlighted.

FIGURE 9
www.frontiersin.org

Figure 9. Imaginary part of the electric modulus (M′′) at isochronous conditions as a function of temperature, for [C2mim][IM14]. The temperature/frequency region associated to the main α process is shown.

FIGURE 10
www.frontiersin.org

Figure 10. Relaxation map for the primary α-process for [C2mim][IM14] and [C8mim][IM14]. Viscosity data obtained in the present work and from the literature [GBA refers to Appetecchi et al. (2011)] and characteristic times from the imaginary part of the electrical modulus, M′′, are reported. The latter data sets have been vertically shifted in the logarithmic scale, taking advantage of the tTS that applies when dealing with different techniques probing the same relaxation process. Lines refer to the VFT fitting of the experimental data sets.

Conclusion

Here we have reported a detailed characterization of the triphilic nature of morphology in a series of imidazolium based ILs bearing a fluorinated anion. In particular samples with an alkyl chain in the cation and a perfluorobutyl moiety in the anion have been investigated. The synergic exploitation of complementary scattering techniques and Molecular Dynamics simulations allows extracting precise structural information at atomistic level. The complementarity of X-ray and neutron scattering is crucial to assess a specific feature of this class of FILs: It is clear that due to contrast reasons, X-ray scattering can detect only the existence of segregated clusters formed by alkyl chains. Neutron scattering instead allows detection of the existence of both alkyl and perfluoro alkyl segregated domains. Our careful MD simulations were able to accurately reproduce this behavior as well as the whole experimental data sets. Complementary approaches based on either diffraction pattern decomposition or aggregates analysis allowed to quantify the extent of (alkyl or perfluoroalkyl) chains clustering: upon increasing the alkyl chain length a complex transition from small to percolating alkyl chain clusters is observed; such a transition develops simultaneously with the inverse transition from percolating to small perfluoroalkyl chains clusters.

We also accounted for experimental determination of the diffusion related relaxation process for the case of [C2mim][IM14] and [C8mim][IM14], by determining their relaxation map and the fragility index, indicating their fragile behavior.

Author Contributions

FLC developed MD simulations. GA and ES synthetized and characterized samples. MZ and EC executed viscosity and NMR measurements. UK, AT, and OR executed SANS experiments. AT and OR executed x-ray/neutron experiments. FLC, GA, EC, AT, and OR discussed the results. FLC, AT, and OR wrote the manuscript. All the authors revised the manuscript.

Funding

This work has been supported by the University of Rome Sapienza Project: Microscopic and mesoscopic organization in ionic liquid-based systems. (RG11715C7CC660BE). MZ and EC gratefully acknowledge support from the U.S. NSF through grant CHE-1664809.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor declared a shared affiliation, though no other collaboration, with one of the authors FLC.

Acknowledgments

AT and OR gratefully acknowledge access to the Japan Synchrotron Radiation Research Institute (JASRI-Spring8) and the kind and skilful support from Dr. S. Kohara in exploiting beamline BL04B2. This research also used resources of the Advanced Photon Source, a U.S. Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357. AT and OR acknowledge the competent support from Dr. B. Aoun in exploiting beamline 11-ID-C at APS. We acknowledge Helmholtz Zentrum Berlin (HZB) for the allocation of neutron beamtime. Access to the dielectric spectroscopy setup has been kindly provided by Dr. M. Russina (Institute for Soft Matter and Functional Materials at the HZB). This work has also been performed within the Agreement No 0018318 (02/06/2014) between STFC and CNR, concerning collaboration in scientific research at the spallation neutron source ISIS and with partial financial support of CNR. Beamtime awarded by ISIS under proposal RB 1700047 and the support from Dr. T. G. A. Youngs are gratefully acknowledged.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2019.00285/full#supplementary-material

References

Almantariotis, D., Gefflaut, T., Pádua, A. A. H., Coxam, J.-Y., and Costa Gomes, M. F. (2010). Effect of fluorination and size of the Alkyl side-chain on the solubility of carbon dioxide in 1-alkyl-3-methylimidazolium Bis(Trifluoromethylsulfonyl)Amide ionic liquids. J. Phys. Chem. B 114, 3608–3617. doi: 10.1021/jp912176n

PubMed Abstract | CrossRef Full Text | Google Scholar

Almantariotis, D., Pensado, A. S., Gunaratne, H. Q. N. N., Hardacre, C., Pádua, A. A. H. H., Coxam, J.-Y., et al. (2017). Influence of fluorination on the solubilities of carbon dioxide, ethane, and nitrogen in 1- n -fluoro-alkyl-3-methylimidazolium Bis( n -Fluoroalkylsulfonyl)Amide ionic liquids. J. Phys. Chem. B 121, 426–436. doi: 10.1021/acs.jpcb.6b10301

PubMed Abstract | CrossRef Full Text | Google Scholar

Angell, C. (1991). Relaxation in liquids, polymers and plastic crystals — strong/fragile patterns and problems. J. Non. Cryst. Solids 131–133, 13–31. doi: 10.1016/0022-3093(91)90266-9

CrossRef Full Text | Google Scholar

Angell, C. A. (1985). Spectroscopy simulation and scattering, and the medium range order problem in glass. J. Non. Cryst. Solids 73, 1–17. doi: 10.1016/0022-3093(85)90334-5

CrossRef Full Text | Google Scholar

Annapureddy, H. V. R., Kashyap, H. K., De Biase, P. M., and Margulis, C. J. (2010). What is the origin of the prepeak in the X-ray scattering of imidazolium-based room-temperature ionic liquids? J. Phys. Chem. B 114, 16838–16846. doi: 10.1021/jp108545z

PubMed Abstract | CrossRef Full Text | Google Scholar

Aoun, B., Goldbach, A., Kohara, S., Wax, J.-F., Gonzalez, M. A., and Saboungi, M.-L. (2010). Structure of a prototypic ionic liquid: ethyl-methylimidazolium bromide. J. Phys. Chem. B 114, 12623–12628. doi: 10.1021/jp1070715

PubMed Abstract | CrossRef Full Text | Google Scholar

Appetecchi, G. B., Montanino, M., Carewska, M., Moreno, M., Alessandrini, F., and Passerini, S. (2011). Chemical–physical properties of Bis(Perfluoroalkylsulfonyl)Imide-based ionic liquids. Electrochim. Acta 56, 1300–1307. doi: 10.1016/j.electacta.2010.10.023

CrossRef Full Text | Google Scholar

Araque, J. C., Hettige, J. J., and Margulis, C. J. (2015). Modern room temperature ionic liquids a simple guide to understanding their structure and how it may relate to dynamics. J. Phys. Chem. B 119, 12727–12740. doi: 10.1021/acs.jpcb.5b05506

PubMed Abstract | CrossRef Full Text | Google Scholar

Atkin, R., and Warr, G. G. (2008). The smallest amphiphiles: nanostructure in protic room-temperature ionic liquids with short alkyl groups. J. Phys. Chem. B 112, 4164–4166. doi: 10.1021/jp801190u

PubMed Abstract | CrossRef Full Text | Google Scholar

Bastos, J. C., Carvalho, S. F., Welton, T., Canongia Lopes, J. N., Rebelo, L. P. N., Shimizu, K., et al. (2018). Design of task-specific fluorinated ionic liquids: nanosegregation: versus hydrogen-bonding ability in aqueous solutions. Chem. Commun. 54, 3524–3527. doi: 10.1039/C8CC00361K

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardes, C. E. S. (2017). AGGREGATES: finding structures in simulation results of solutions. J. Comput. Chem. 38, 753–765. doi: 10.1002/jcc.24735

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardes, C. E. S., and Joseph, A. (2015). Evaluation of the OPLS-AA force field for the study of structural and energetic aspects of molecular organic crystals. J. Phys. Chem. A 119, 3023–3034. doi: 10.1021/jp512349r

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardes, C. E. S., Minas da Piedade, M. E., and Canongia Lopes, J. N. (2011). The structure of aqueous solutions of a hydrophilic ionic liquid: the full concentration range of 1-ethyl-3-methylimidazolium ethylsulfate and water. J. Phys. Chem. B 115, 2067–2074. doi: 10.1021/jp1113202

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardes, C. E. S., Shimizu, K., Lobo Ferreira, A. I. M. C., Santos, L. M. N. B. F., and Canongia Lopes, J. N. (2014). Structure and aggregation in the 1,3-dialkyl-imidazolium Bis(Trifluoromethylsulfonyl)Imide ionic liquid family: 2. From single to double long alkyl side chains. J. Phys. Chem. B 118, 6885–6895. doi: 10.1021/jp502968u

PubMed Abstract | CrossRef Full Text | Google Scholar

Böhmer, R., and Angell, C. A. (1992). Correlations of the nonexponentiality and state dependence of mechanical relaxations with bond connectivity in Ge-As-Se supercooled liquids. Phys. Rev. B 45, 10091–10094. doi: 10.1103/PhysRevB.45.10091

PubMed Abstract | CrossRef Full Text | Google Scholar

Böhmer, R., Ngai, K. L., Angell, C. A., and Plazek, D. J. (1993). Nonexponential relaxations in strong and fragile glass formers. J. Chem. Phys. 99, 4201–4209. doi: 10.1063/1.466117

CrossRef Full Text | Google Scholar

Brehm, M., and Kirchner, B. (2011). TRAVIS - a free analyzer and visualizer for monte carlo and molecular dynamics trajectories. J. Chem. Inf. Model. 51, 2007–2023. doi: 10.1021/ci200217w

PubMed Abstract | CrossRef Full Text | Google Scholar

Brehm, M., Weber, H., Thomas, M., Hollóczki, O., and Kirchner, B. (2015). Domain analysis in nanostructured liquids: a post-molecular dynamics study at the example of ionic liquids. Chem. Phys. Chem. 16, 3271–3277. doi: 10.1002/cphc.201500471

PubMed Abstract | CrossRef Full Text | Google Scholar

Bussi, G., Donadio, D., and Parrinello, M. (2007). Canonical sampling through velocity rescaling. J. Chem. Phys. 126:014101. doi: 10.1063/1.2408420

PubMed Abstract | CrossRef Full Text | Google Scholar

Canongia Lopes, J. N., Costa Gomes, M. F., and Pádua, A. A. H. (2006). Nonpolar, polar, and associating solutes in ionic liquids. J. Phys. Chem. B 110, 16816–16818. doi: 10.1021/jp063603r

PubMed Abstract | CrossRef Full Text | Google Scholar

Canongia Lopes, J. N., and Padua, A. A. H. (2006). Nanostructural organization in ionic liquids. J. Phys. Chem. B 110:3330–3335. doi: 10.1021/jp056006y

PubMed Abstract | CrossRef Full Text | Google Scholar

Castiglione, F., Famulari, A., Raos, G., Meille, S. V., Mele, A., Appetecchi, G. B., et al. (2014). Pyrrolidinium-based ionic liquids doped with lithium salts: how does Li+coordination affect its diffusivity? J. Phys. Chem. B 118, 13679–13688. doi: 10.1021/jp509387r

PubMed Abstract | CrossRef Full Text | Google Scholar

Cromer, D. T. (1969). Compton scattering factors for aspherical free atoms. J. Chem. Phys. 50, 4857–4859. doi: 10.1063/1.1670980

CrossRef Full Text | Google Scholar

Darden, T., York, D., and Pedersen, L. (1993). Particle mesh ewald: an N·log(N) method for ewald sums in large systems. J. Chem. Phys. 98, 10089–10092. doi: 10.1063/1.464397

CrossRef Full Text | Google Scholar

Essmann, U., Perera, L., Berkowitz, M. L., Darden, T., Lee, H., and Pedersen, L. G. (1995). A smooth particle mesh ewald method. J. Chem. Phys. 103, 8577–8593. doi: 10.1063/1.470117

CrossRef Full Text | Google Scholar

Ferreira, M. L., Pastoriza-Gallego, M. J., Araújo, J. M. M., Canongia Lopes, J. N., Rebelo, L. P. N. M., Piñeiro, M., et al. (2017). Influence of nanosegregation on the phase behavior of fluorinated ionic liquids. J. Phys. Chem. C 121, 5415–5427. doi: 10.1021/acs.jpcc.7b00516

CrossRef Full Text | Google Scholar

Fulcher, G. S. (1925). Analysis of recent measurements of the viscosity of glasses. J. Am. Ceram. Soc. 8, 339–355. doi: 10.1111/j.1151-2916.1925.tb16731.x

CrossRef Full Text | Google Scholar

Funston, A. M., Fadeeva, T. A., Wishart, J. F., and Castner, E. W. (2007). Fluorescence probing of temperature-dependent dynamics and friction in ionic liquid local environments. J. Phys. Chem. B 111, 4963–4977. doi: 10.1021/jp068298o

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilles, R., Keiderling, U., Strunz, P., Wiedenmann, A., and Fuess, H. (2000). Silver behenate as a standard for instrumental resolution and wavelength calibration for small angle neutron scattering. Mater. Sci. Forum 321–323, 264–269. doi: 10.4028/www.scientific.net/MSF.321-324.264

CrossRef Full Text | Google Scholar

Greaves, T. L., Kennedy, D. F., Shen, Y., Hawley, A., Song, G., and Drummond, C. J. (2013). Fluorous protic ionic liquids exhibit discrete segregated nano-scale solvent domains and form new populations of nano-scale objects upon primary alcohol addition. Phys. Chem. Chem. Phys. 15, 7592–7598. doi: 10.1039/c3cp44589e

PubMed Abstract | CrossRef Full Text | Google Scholar

Hammersley, A. P., Svensson, S. O., Hanfland, M., Fitch, A. N., and Hausermann, D. (1996). Two-dimensional detector software: from real detector to idealised image or two-theta scan. High Press. Res. 14, 235–248. doi: 10.1080/08957959608201408

CrossRef Full Text | Google Scholar

Harris, K. R. (2009). The fractional stokes–einstein equation: application to lennard-jones, molecular, and ionic liquids. J. Chem. Phys. 131:054503. doi: 10.1063/1.3183951

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayes, R., Imberti, S., Warr, G. G., and Atkin, R. (2011). Amphiphilicity determines nanostructure in protic ionic liquids. Phys. Chem. Chem. Phys. 13, 3237–3247. doi: 10.1039/C0CP01137A

PubMed Abstract | CrossRef Full Text | Google Scholar

Hess, B., Kutzner, C., Van Der Spoel, D., and Lindahl, E. (2008). GRGMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4, 435–447. doi: 10.1021/ct700301q

PubMed Abstract | CrossRef Full Text | Google Scholar

Hettige, J. J., Araque, J. C., and Margulis, C. J. (2014). Bicontinuity and multiple length scale ordering in triphilic hydrogen-bonding ionic liquids. J. Phys. Chem. B 118, 12706–12716. doi: 10.1021/jp5068457

PubMed Abstract | CrossRef Full Text | Google Scholar

Hollóczki, O., Macchiagodena, M., Weber, H., Thomas, M., Brehm, M., Stark, A., et al. (2015). Triphilic ionic-liquid mixtures: fluorinated and non-fluorinated aprotic ionic-liquid mixtures. Chem. Phys. Chem. 16, 3325–3333. doi: 10.1002/cphc.201500473

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunter, J. D. (2007). Matplotlib: a 2D graphics environment. Comput. Sci. Eng. 9, 90–95. doi: 10.1109/MCSE.2007.55

CrossRef Full Text | Google Scholar

Jeremias, S., Carewska, M., Conte, L., Passerini, S., and Appetecchi, G. B. (2013). Asymmetry effect of novel per(Fluoroalkylsulfonyl)Imide anions in pyrrolidinium ionic liquids. RSC Adv. 3:17755. doi: 10.1039/c3ra42980f

CrossRef Full Text | Google Scholar

Johansson, P., Fast, L. E., Matic, A., Appetecchi, G. B., and Passerini, S. (2010). The conductivity of pyrrolidinium and sulfonylimide-based ionic liquids: a combined experimental and computational study. J. Power Sources 195, 2074–2076. doi: 10.1016/j.jpowsour.2009.10.029

CrossRef Full Text | Google Scholar

Kashyap, H. K., Hettige, J. J., Annapureddy, H. V. R., and Margulis, C. J. (2012). SAXS anti-peaks reveal the length-scales of dual positive-negative and polar-apolar ordering in room-temperature ionic liquids. Chem. Commun. 48, 5103–5105. doi: 10.1039/c2cc30609c

PubMed Abstract | CrossRef Full Text | Google Scholar

Kashyap, H. K., Santos, C. S., Daly, R. P., Hettige, J. J., Murthy, N. S., Shirota, H., et al. (2013). How does the ionic liquid organizational landscape change when nonpolar cationic alkyl groups are replaced by polar isoelectronic diethers? J. Phys. Chem. B 117, 1130–1135. doi: 10.1021/jp311032p

PubMed Abstract | CrossRef Full Text | Google Scholar

Keiderling, U. (2002). The new “BerSANS-PC” software for reduction and treatment of small angle neutron scattering data. Appl. Phys. A Mater. Sci. Process. 74, s1455–s1457. doi: 10.1007/s003390201561

CrossRef Full Text | Google Scholar

Keiderling, U., and Wiedenmann, A. (1995). New SANS instrument at the BER II reactor in Berlin, Germany. Phys. B 213–214, 895–897. doi: 10.1016/0921-4526(95)00316-2

CrossRef Full Text | Google Scholar

Kim, J., Singh, R. P., and Shreeve, J. M. (2004). Low melting inorganic salts of alkyl-, fluoroalkyl-, alkyl ether-, and fluoroalkyl ether-substituted oxazolidine and morpholine. Inorg. Chem. 43, 2960–2966. doi: 10.1021/ic049916n

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohara, S., Suzuya, K., Kashihara, Y., Matsumoto, N., Umesaki, N., and Sakai, I. (2001). A horizontal two-axis diffractometer for high-energy X-ray diffraction using synchrotron radiation on bending magnet beamline BL04B2 at SPring-8. Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 467–468, 1030–1033. doi: 10.1016/S0168-9002(01)00630-1

CrossRef Full Text | Google Scholar

Kunze, M., Paillard, E., Jeong, S., Appetecchi, G. B., Sch, M., Winter, M., et al. (2011). Inhibition of self-aggregation in ionic liquid electrolytes for high-energy electrochemical devices. J. Phys. Chem. C 115, 19431–19436. doi: 10.1021/jp2055969

CrossRef Full Text | Google Scholar

Leys, J., Wübbenhorst, M., Preethy Menon, C., Rajesh, R., Thoen, J., Glorieux, C., et al. (2008). Temperature dependence of the electrical conductivity of imidazolium ionic liquids. J. Chem. Phys. 128:064509. doi: 10.1063/1.2827462

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Bañuelos, J. L., Guo, J., Anovitz, L., Rother, G., Shaw, R. W., et al. (2012). Alkyl chain length and temperature effects on structural properties of pyrrolidinium-based ionic liquids: a combined atomistic simulation and small-angle X-ray scattering study. J. Phys. Chem. Lett. 3, 125–130. doi: 10.1021/jz2013209

CrossRef Full Text | Google Scholar

Liang, M., Khatun, S., and Castner, E. W. (2015). Communication: unusual structure and transport in ionic liquid-hexane mixtures. J. Chem. Phys. 142, 12–16. doi: 10.1063/1.4916388

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo Celso, F., Appetecchi, G. B., Jafta, C. J., Gontrani, L., Canongia Lopes, J. N., Triolo, A., et al. (2018a). Nanoscale organization in the fluorinated room temperature ionic liquid: tetraethyl ammonium (Trifluoromethanesulfonyl)(Nonafluorobutylsulfonyl)Imide. J. Chem. Phys. 148:193816. doi: 10.1063/1.5016236

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo Celso, F., Yoshida, Y., Castiglione, F., Ferro, M., Mele, A., Jafta, C. J., et al. (2017). Direct experimental observation of mesoscopic fluorous domains in fluorinated room temperature ionic liquids. Phys. Chem. Chem. Phys. 19, 13101–13110. doi: 10.1039/C7CP01971H

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo Celso, F., Yoshida, Y., Lombardo, R., Jafta, C., Gontrani, L., Triolo, A., et al. (2018b). Mesoscopic structural organization in fluorinated room temperature ionic liquids. Comptes Rendus Chim. 21, 757–770. doi: 10.1016/j.crci.2018.02.001

CrossRef Full Text | Google Scholar

Lopes, J. N. C., and Padua, A. A. H. (2004). Molecular force field for ionic liquids composed of triflate or bistriflyimide anions. J. Phys. Chem. B 108, 16893–16898. doi: 10.1021/jp0476545

CrossRef Full Text | Google Scholar

Macchiagodena, M., Ramondo, F., Triolo, A., Gontrani, L., and Caminiti, R. (2012). Liquid structure of 1-ethyl-3-methylimidazolium alkyl sulfates by X-ray scattering and molecular dynamics. J. Phys. Chem. B 116, 13448–13458. doi: 10.1021/jp306982e

PubMed Abstract | CrossRef Full Text | Google Scholar

Mariani, A., Bonomo, M., Wu, B., Centrella, B., Dini, D., Castner, E. W., et al. (2017). Intriguing transport dynamics of ethylammonium nitrate-acetonitrile binary mixtures arising from nano-inhomogeneity. Phys. Chem. Chem. Phys. 19, 27212–27220. doi: 10.1039/C7CP04592A

PubMed Abstract | CrossRef Full Text | Google Scholar

Merrigan, T. L., Bates, E. D., Dorman, S. C., and Davis, J. H. (2000). New fluorous ionic liquids function as surfactants in conventional room-temperature ionic liquids. Chem. Commun. 2051–2052. doi: 10.1039/b005418f

CrossRef Full Text | Google Scholar

Montanino, M., Alessandrini, F., Passerini, S., and Appetecchi, G. B. (2013). Water-based synthesis of hydrophobic ionic liquids for high-energy electrochemical devices. Electrochim. Acta 96, 124–133. doi: 10.1016/j.electacta.2013.02.082

CrossRef Full Text | Google Scholar

Montanino, M., Moreno, M., Alessandrini, F., Appetecchi, G. B., Passerini, S., Zhou, Q., et al. (2012). Physical and electrochemical properties of binary ionic liquid mixtures: (1–x) PYR14TFSI–(x) PYR14IM14. Electrochim. Acta 60, 163–169. doi: 10.1016/j.electacta.2011.11.030

CrossRef Full Text | Google Scholar

Ohmura, A., Sato, K., Hamaya, N., Isshiki, M., and Ohishi, Y. (2009). Structure of pressure-induced amorphous form of SnI4 at high pressure. Phys. Rev. B 80:054201. doi: 10.1103/PhysRevB.80.054201

CrossRef Full Text | Google Scholar

Pádua, A. A. H., Costa Gomes, M. F., and Canongia Lopes, J. N. A. (2007). Molecular Solutes in Ionic Liquids: A Structural Perspective. Acc. Chem. Res. 40, 1087–1096. doi: 10.1021/ar700050q

PubMed Abstract | CrossRef Full Text | Google Scholar

Palumbo, O., Trequattrini, F., Appetecchi, G. B., and Paolone, A. (2017). A study of the conformers of the (Nonafluorobutanesulfonyl)Imide ion by means of infrared spectroscopy and Density Functional Theory (DFT) calculations. Challenges 8:7. doi: 10.3390/challe8010007

CrossRef Full Text | Google Scholar

Paolone, A., Palumbo, O., Trequattrini, F., and Appetecchi, G. B. (2018). Relaxational dynamics in the PYR14-IM14 ionic liquid by mechanical spectroscopy. Mater. Res. 21, 1–6. doi: 10.1590/1980-5373-mr-2017-0870

CrossRef Full Text | Google Scholar

Parrinello, M., and Rahman, A. (1981). Polymorphic transitions in single crystals: a new molecular dynamics method. J. Appl. Phys. 52, 7182–7190. doi: 10.1063/1.328693

CrossRef Full Text | Google Scholar

Pereiro, A. B., Araújo, J. M. M., Esperança, J. M. S. S., and Rebelo, L. P. N. (2018). “Chapter 4, Surfactant fluorinated ionic liquids,” in RSC Smart Materials, eds A. Eftekhari, Vol. 2018 (The Royal Society of Chemistry) 79–102. doi: 10.1039/9781788011839-00079

CrossRef Full Text | Google Scholar

Pereiro, A. B., Araújo, J. M. M., Martinho, S., Alves, F., Nunes, S., Matias, A., et al. (2013a). Fluorinated ionic liquids: properties and applications. ACS Sustain. Chem. Eng. 1, 427–439. doi: 10.1021/sc300163n

CrossRef Full Text | Google Scholar

Pereiro, A. B., Pastoriza-Gallego, M. J., Shimizu, K., Marrucho, I. M., Lopes, J. N. C., Piñeiro, M. M., et al. (2013b). On the formation of a third, nanostructured domain in ionic liquids. J. Phys. Chem. B 117, 10826–10833. doi: 10.1021/jp402300c

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, X., Thompson, J. W., and Billinge, S. J. L. (2004). PDFgetX2 : a GUI-driven program to obtain the pair distribution function from X-ray powder diffraction data. J. Appl. Crystallogr. 37, 678–678. doi: 10.1107/S0021889804011744

CrossRef Full Text | Google Scholar

Rauber, D., Zhang, P., Hempelmann, R., and Huch, V. (2017). Lamellar structures in fluorinated phosphonium ionic liquids : the roles of fluorination and chain length. Phys. Chem. Chem. Phys. 19, 27251–27258. doi: 10.1039/C7CP04814A

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, Y. (2012). High-energy synchrotron X-ray diffraction and its application to in situ structural phase-transition studies in complex sample environments. JOM 64, 140–149. doi: 10.1007/s11837-011-0218-8

CrossRef Full Text | Google Scholar

Rocha, M. A. A., Neves, C. M. S. S., Freire, M. G., Russina, O., Triolo, A., Coutinho, J. A. P., et al. (2013). Alkylimidazolium based ionic liquids: impact of cation symmetry on their nanoscale structural organization. J. Phys. Chem. B 117, 10889–10897. doi: 10.1021/jp406374a

PubMed Abstract | CrossRef Full Text | Google Scholar

Russina, O., Lo Celso, F., Di Michiel, M., Passerini, S., Appetecchi, G. B., Castiglione, F., et al. (2013). Mesoscopic structural organization in triphilic room temperature ionic liquids. Faraday Discuss. 167, 499–513. doi: 10.1039/c3fd00056g

PubMed Abstract | CrossRef Full Text | Google Scholar

Russina, O., Lo Celso, F., Plechkova, N., Jafta, C. J., Appetecchi, G. B., and Triolo, A. (2017a). Mesoscopic organization in ionic liquids. Top. Curr. Chem. 375:58. doi: 10.1007/s41061-017-0147-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Russina, O., Lo Celso, F., Plechkova, N. V., and Triolo, A. (2017b). Emerging evidences of mesoscopic-scale complexity in neat ionic liquids and their mixtures. J. Phys. Chem. Lett. 8, 1197–1204. doi: 10.1021/acs.jpclett.6b02811

PubMed Abstract | CrossRef Full Text | Google Scholar

Russina, O., and Triolo, A. (2012). New experimental evidence supporting the mesoscopic segregation model in room temperature ionic liquids. Faraday Discuss. 154, 97–109. doi: 10.1039/C1FD00073J

PubMed Abstract | CrossRef Full Text | Google Scholar

Russina, O., and Triolo, A. (2017). “Ionic liquids and neutron scattering,” in Experimental Methods in the Physical Sciences, Vol. 49, eds F. Fernandez-Alonso and D. L. Price (Academic Press), 213–278. doi: 10.1016/B978-0-12-805324-9.00004-2

CrossRef Full Text | Google Scholar

Russina, O., Triolo, A., Gontrani, L., and Caminiti, R. (2012). Mesoscopic structural heterogeneities in room-temperature ionic liquids. J. Phys. Chem. Lett. 3, 27–33. doi: 10.1021/jz201349z

CrossRef Full Text | Google Scholar

Russina, O., Triolo, A., Gontrani, L., Caminiti, R., Xiao, D., Hines, L. G. Jr, et al. (2009). Morphology and intermolecular dynamics of 1-alkyl-3-methylimidazolium Bis{(Trifluoromethane)Sulfonyl}amide ionic liquids: structural and dynamic evidence of nanoscale segregation. J. Phys. Condens. Matter 21:424121. doi: 10.1088/0953-8984/21/42/424121

CrossRef Full Text | Google Scholar

Shen, Y., Kennedy, D. F., Greaves, T. L., Weerawardena, A., Mulder, R. J., Kirby, N., et al. (2012). Protic ionic liquids with fluorous anions: physicochemical properties and self-assembly nanostructure. Phys. Chem. Chem. Phys. 14, 7981–7992. doi: 10.1039/c2cp40463j

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimizu, K., Almantariotis, D., Gomes, M. F. C., Pádua, A. A. H., and Canongia Lopes, J. N. (2010). Molecular force field for ionic liquids V: hydroxyethylimidazolium, dimethoxy-2- methylimidazolium, and fluoroalkylimidazolium cations and Bis(Fluorosulfonyl)Amide, perfluoroalkanesulfonylamide, and fluoroalkylfluorophosphate anions. J. Phys. Chem. B 114, 3592–3600. doi: 10.1021/jp9120468

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimizu, K., Bernardes, C. E. S., and Canongia Lopes, J. N. (2014). Structure and aggregation in the 1-alkyl-3-methylimidazolium Bis(Trifluoromethylsulfonyl)Imide ionic liquid homologous series. J. Phys. Chem. B 118, 567–576. doi: 10.1021/jp409987d

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimizu, K., Freitas, A. A., and Canongia Lopes, J. N. (2017). Structural characterization of the [CnC1im][C4F9SO3] ionic liquid series: alkyl versus perfluoroalkyl side chains. J. Mol. Liq. 226, 28–34. doi: 10.1016/j.molliq.2016.08.014

CrossRef Full Text | Google Scholar

Sippel, P., Lunkenheimer, P., Krohns, S., Thoms, E., and Loidl, A. (2015). Importance of liquid fragility for energy applications of ionic liquids. Sci. Rep. 5:13922. doi: 10.1038/srep13922

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, G. D., Borodin, O., Magda, J. J., Boyd, R. H., Wang, Y., Bara, J. E., et al. (2010). A comparison of fluoroalkyl-derivatized imidazolium:TFSI and Alkyl-derivatized imidazolium:TFSI ionic liquids: a molecular dynamics simulation study. Phys. Chem. Chem. Phys. 12, 7064–7076. doi: 10.1039/c001387k

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, X., Hamano, H., Minofar, B., Kanzaki, R., Fujii, K., Kameda, Y., et al. (2012). Structural heterogeneity and unique distorted hydrogen bonding in primary ammonium nitrate ionic liquids studied by high-energy X-ray diffraction experiments and MD simulations. J. Phys. Chem. B 116, 2801–2813. doi: 10.1021/jp209561t

PubMed Abstract | CrossRef Full Text | Google Scholar

Suarez, S. N., Rúa, A., Cuffari, D., Pilar, K., Hatcher, J. L., and Ramati, S. F. (2015). Do TFSA anions slither? Pressure exposes the role of TFSA conformational exchange in self-diffusion. J. Phys. Chem. B 119, 14756–14765. doi: 10.1021/acs.jpcb.5b08658

PubMed Abstract | CrossRef Full Text | Google Scholar

Tammann, G., and Hesse, W. (1926). Die abhängigkeit der viscosität von der temperatur bie unterkühlten flüssigkeiten. Zeitschrift Anorg. Allg. Chemie 156, 245–257. doi: 10.1002/zaac.19261560121

CrossRef Full Text | Google Scholar

Tao, R., Gurung, E., Cetin, M. M., Mayer, M. F., Quitevis, E. L., and Simon, S. L. (2017). Fragility of ionic liquids measured by flash differential scanning calorimetry. Thermochim. Acta 654, 121–129. doi: 10.1016/j.tca.2017.05.008

CrossRef Full Text | Google Scholar

Triolo, A., Russina, O., Bleif, H.-J., and Di Cola, E. (2007). Nanoscale segregation in room temperature ionic liquids. J. Phys. Chem. B 111, 4641–4644. doi: 10.1021/jp067705t

PubMed Abstract | CrossRef Full Text | Google Scholar

Triolo, A., Russina, O., Fazio, B., Appetecchi, G. B., Carewska, M., and Passerini, S. (2009). Nanoscale organization in piperidinium-based room temperature ionic liquids. J. Chem. Phys. 130:164521. doi: 10.1063/1.3119977

PubMed Abstract | CrossRef Full Text | Google Scholar

Triolo, A., Russina, O., Fazio, B., Triolo, R., and Di Cola, E. (2008). Morphology of 1-alkyl-3-methylimidazolium hexafluorophosphate room temperature ionic liquids. Chem. Phys. Lett. 457, 362–365. doi: 10.1016/j.cplett.2008.04.027

CrossRef Full Text | Google Scholar

Tsurumaki, A., and Ohno, H. (2018). Dissolution of Oligo(Tetrafluoroethylene) and preparation of Poly(Tetrafluoroethylene)-based composites by using fluorinated ionic liquids. Chem. Commun. 54, 409–412. doi: 10.1039/C7CC08449H

PubMed Abstract | CrossRef Full Text | Google Scholar

Urahata, S. M., and Ribeiro, M. C. C. (2004). Structure of ionic liquids of 1-alkyl-3-methylimidazolium cations: a systematic computer simulation study. J. Chem. Phys. 120, 1855–1863. doi: 10.1063/1.1635356

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Broeke, J., Winter, F., Deelman, B., and van Koten, G. (2002). A highly fluorous room-temperature ionic liquid exhibiting fluorous biphasic behavior and its use in catalyst recycling. Org. Lett. 4, 3851–3854. doi: 10.1021/ol026700l

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Der Spoel, D., Lindahl, E., Hess, B., Groenhof, G., Mark, A. E., and Berendsen, H. J. C. (2005). GROMACS: fast, flexible, and free. J. Comput. Chem. 26, 1701–1718. doi: 10.1002/jcc.20291

PubMed Abstract | CrossRef Full Text | Google Scholar

Vieira, N. S. M., Reis, P. M., Shimizu, K., Cortes, O. A., Marrucho, I. M., Araújo, J. M. M., et al. (2015). A thermophysical and structural characterization of ionic liquids with alkyl and perfluoroalkyl side chains. RSC Adv. 5, 65337–65350. doi: 10.1039/C5RA13869H

CrossRef Full Text | Google Scholar

Wang, Y., and Voth, G. A. (2005). Unique spatial heterogeneity in ionic liquids. J. Am. Chem. Soc. 127, 12192–12193. doi: 10.1021/ja053796g

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., and Voth, G. A. (2006). Tail aggregation and domain diffusion in ionic liquids. J. Phys. Chem. B 110, 18601–18608. doi: 10.1021/jp063199w

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, H., Hollóczki, O., Pensado, A. S., and Kirchner, B. (2013). Side chain fluorination and anion effect on the structure of 1-butyl-3-methylimidazolium ionic liquids. J. Chem. Phys. 139:084502. doi: 10.1063/1.4818540

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiss, H., Mars, J., Li, H., Kircher, G., Ivanova, O., Feoktystov, A., et al. (2017). Mesoscopic correlation functions in heterogeneous ionic liquids. J. Phys. Chem. B 121, 620–629. doi: 10.1021/acs.jpcb.6b11220

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, B., Liang, M., Maroncelli, M., and Castner, E. W. (2015). Photoinduced bimolecular electron transfer from cyano anions in ionic liquids. J. Phys. Chem. B 119, 14790–14799. doi: 10.1021/acs.jpcb.5b09216

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, B., Liang, M., Zmich, N., Hatcher, J., Lall-Ramnarine, S. I., Wishart, J. F., et al. (2018). Photoinduced bimolecular electron transfer in ionic liquids: cationic electron donors. J. Phys. Chem. B 122, 2379–2388. doi: 10.1021/acs.jpcb.7b12542

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, B., Yamashita, Y., Endo, T., Takahashi, K., and Castner, E. W. (2016). Structure and dynamics of ionic liquids: trimethylsilylpropyl-substituted cations and Bis(Sulfonyl)Amide anions. J. Chem. Phys. 145:244506. doi: 10.1063/1.4972410

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, D. H., Chen, A., and Johnson, C. S. (1995). An improved diffusion-ordered spectroscopy experiment incorporating bipolar-gradient pulses. J. Magn. Reson. Ser. A 115, 260–264. doi: 10.1006/jmra.1995.1176

CrossRef Full Text | Google Scholar

Xue, H., and Shreeve, J. M. (2005). Ionic liquids with fluorine-containing cations. Eur. J. Inorg. Chem. 13, 2573–2580. doi: 10.1002/ejic.200500129

CrossRef Full Text | Google Scholar

Xue, H., Verma, R., and Shreeve, J. M. (2006). Review of ionic liquids with fluorine-containing anions. J. Fluor. Chem. 127, 159–176. doi: 10.1016/j.jfluchem.2005.11.007

CrossRef Full Text | Google Scholar

Yoshida, Y., and Saito, G. (2011). Ionic liquids based on Diethylmethyl(2-Methoxyethyl)Ammonium cations and Bis(Perfluoroalkanesulfonyl)Amide anions: influence of anion structure on liquid properties. Phys. Chem. Chem. Phys. 13, 20302–20310. doi: 10.1039/c1cp21783f

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, W., Mohammed, A., Hines, L. G., Xiao, D., Martinez, O. J., Bartsch, R. A., et al. (2011). Effect of cation symmetry on the morphology and physicochemical properties of imidazolium ionic liquids. J. Phys. Chem. B 115, 6572–6584. doi: 10.1021/jp1115614

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: fluorous tail, triphilic, ionic liquid, neutron scattering, molecular dynamics (MD)

Citation: Lo Celso F, Appetecchi GB, Simonetti E, Zhao M, Castner EW Jr, Keiderling U, Gontrani L, Triolo A and Russina O (2019) Microscopic Structural and Dynamic Features in Triphilic Room Temperature Ionic Liquids. Front. Chem. 7:285. doi: 10.3389/fchem.2019.00285

Received: 30 November 2018; Accepted: 08 April 2019;
Published: 02 May 2019.

Edited by:

Francesca D'Anna, Università degli Studi di Palermo, Italy

Reviewed by:

Barbara Kirchner, University of Bonn, Germany
Bhabani S. Mallik, Indian Institute of Technology Hyderabad, India
Adilson Alves De Freitas, Universidade de Lisboa, Portugal

Copyright © 2019 Lo Celso, Appetecchi, Simonetti, Zhao, Castner, Keiderling, Gontrani, Triolo and Russina. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Alessandro Triolo, triolo@ism.cnr.it
Olga Russina, olga.russina@uniroma1.it

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.