Investigating Intermolecular Interactions in a DME-Based Hybrid Ionic Liquid Electrolyte by HOESY NMR

The intermolecular interactions in a hybrid electrolyte based on various compositions of the ionic liquid N-methyl-N-propyl pyrrolidinium bis-fluorosulfonylimide (C3mpyrFSI), LiFSI salt and an ether-based additive, 1,2-dimethoxy ethane (DME), have been investigated using the HOESY (Heteronuclear Overhauser Effect SpectroscopY) NMR experiment. This NMR technique allows a quantification of the intermolecular interactions in ionic liquids (ILs) by measuring the cross-relaxation rate (σ) between different pairs of nuclei. Thereby, we compare the cross-relaxation rates between the cations, anions and DME in these hybrid electrolyte systems using 1H-7Li and 1H-19F HOESY experiments, and interpret the measured parameters in terms of ionic and molecular associations. The results give insights into the local coordination environment of the Li+ cations and their solvation by the FSI anions and DME.


INTRODUCTION
Ionic liquids (ILs), typically consisting of an organic cation and an inorganic anion, have gained a lot of attention from researchers over the last two decades (Weingärtner, 2013). ILs are generally liquids below 100 • C and have shown a lot of promise as electrolytes in batteries and other energy storage applications due to their versatile properties including low volatility, high thermal and electrochemical stability as well as high ionic conductivity (Macfarlane et al., 2014). Lithium-ion batteries play a very crucial role among current energy storage technologies due to their high energy and power densities (Bruce et al., 2012). They are commonly used in consumer electronics and increasingly in electric vehicles, and hence are very important in addressing global energy concerns. However, current lithium ion battery technologies still pose a safety risk primarily due to the volatile and flammable nature of the commonly used organic electrolytes (Wongittharom et al., 2014;Zhou et al., 2017). IL electrolytes are promising materials to replace these battery electrolytes in particular because of their low volatility and non-flammability. However, these properties typically coincide with lower ionic conductivities compared to traditional electrolyte systems.
In order to advance the design of ILs with improved ionic conductivities for battery applications, their molecular interactions need to be understood (Damodaran, 2016). The location and local environment of Li + relative to the other ions in the electrolyte will play a very important role in the charge transport mechanisms and other processes occurring in the electrolyte (Martin et al., 2018b), which will in turn be major factors in determining the overall performance of the device. It is therefore vital to understand the molecular level structure and dynamics of ILs in order to improve their properties by tailoring their structure and composition.
Nuclear magnetic resonance (NMR) spectroscopy is a useful tool that has helped researchers to probe the structure and dynamics of ionic liquids. Pulsed field gradients (Pope et al., 2016;Hilder et al., 2017), relaxation measurements (Han et al., 2012;Filippov et al., 2015), and dynamic nuclear polarization (Sani et al., 2018) are some NMR techniques that have been used to investigate the structure and dynamics of ILs. Nuclear Overhauser Effect (NOE) based experiments are particularly useful to probe nanostructured arrangements in ILs through the transfer of nuclear spin polarization by crossrelaxation (Damodaran, 2016). In particular, the Heteronuclear Overhauser Effect SpectroscopY (HOESY) pulse sequence (Lingscheid et al., 2012;Bai et al., 2013;Chiappe et al., 2013;Castiglione et al., 2014Castiglione et al., , 2015Giernoth et al., 2014;Tripathi and Saha, 2014) has been used to probe the various molecular interactions in ILs and their interactions with other solvents.
In our recent work on HOESY (Martin et al., 2018b), we demonstrated the importance of taking into account both the longitudinal relaxation times and self-diffusion coefficients of the nuclear spins when fitting HOESY build-up curves, as well as the normalization of the measured signals to the equilibrium magnetization. This analysis allows the quantitative comparison of NOEs (quantified as cross-relaxation rates) among different ILs, concentrations or temperatures. There have however been controversies regarding the interpretation of the intermolecular NOEs measured using this technique, owing to the fact that the internuclear distances are dependent on both the rotational and translational dynamics of the ions. Gabl et al. (2013) proposed that intermolecular NOEs measured between nuclear spins with similar frequencies, such as 1 H and 19 F, are sensitive to longer distances than nuclear spins whose frequencies are far apart. Our very recent work corroborates this prediction (Martin et al., 2018a), demonstrating the validity of interpreting intermolecular NOEs in terms of distances provided that such effects are taken into consideration.
Herein, we apply the HOESY NMR technique to a class of hybrid ionic liquid electrolytes recently reported (Pal et al., 2018) by our group that shows promise for applications in lithium electrochemistry with low polysulfide dissolution that is of importance for Li-S batteries. The addition of DME in these hybrid electrolytes has been shown to increase the relative mobility of the lithium ion through the formation of a chelate compound with lithium. Further understanding of the structure and dynamics in these electrolytes is of interest in terms of their future development. Herein, we quantitatively measure the Overhauser cross-relaxation rates in this DME-based hybrid ionic liquid electrolyte system for the intermolecular heteronuclear spin pairs 1 H-19 F and 1 H-7 Li. The results provide insights into the ionic associations and lithium cation solvation environment via the interactions between the Li + /FSI − and the different hydrogen environments on the pyrrolidinium cations and DME molecules.

Sample Preparation
The ionic liquid N-methyl-N-propyl pyrrolidinium bisfluorosulfonylimide (C 3 mpyrFSI) was purchased from Solvionic (France) with 99.5% purity, 1,2-dimethoxy ethane (DME) was purchased from Sigma-Aldrich with 99.9% purity (anhydrous), and the Li salt of bis-fluorosulfonylimide (LiFSI) was purchased from Nippon Shokubai (Japan) with 99.5% purity and used without further purification. The structures of these compounds are shown in Figure 1. These components were mixed together at different molar ratios as specified in Table 1. A neat sample of the C 3 mpyrFSI IL was also studied. All the materials were stored and packed in an Argon filled glove box. For the NMR experiments, the samples were packed in 3 mm capillary tubes inside the glove box to avoid contamination with water. The samples were then flame-sealed externally and inserted into 5 mm NMR tubes with D 2 O as outer layer for field lock and to minimize thermal convection effects.

Diffusion Coefficient Measurements
Self-diffusion measurements for 1 H, 7 Li, and 19 F were performed at 20 • C using a pulsed field gradient stimulated echo NMR pulse sequence with a Bruker Avance III 7.05 T spectrometer equipped with a 5 mm Bruker Diff50 probe with a maximum gradient strength of 3,000 G/cm. The gradient pulse length used in the experiment was 2 ms and the diffusion time was 25 ms. To extract the diffusion coefficients, which are displayed in Table 2 below, the data were fitted in the Bruker Topspin software using the Stejskal-Tanner equation: where, I is the observed signal intensity, I o is the maximum signal intensity, γ is the gyromagnetic ratio of the nucleus under observation, g is the gradient strength, δ is the gradient pulse duration and is the diffusion time.

Longitudinal Relaxation (T 1 ) Measurements
Longitudinal relaxation times (T 1 ) for 1 H, 7 Li, and 19 F were measured at 20 • C at a magnetic field strength of 11.7 T (500 MHz 1 H frequency) using an inversion recovery pulse sequence and a Bruker Avance III spectrometer equipped with a 5 mm HX solution state probe. The data were fitted to an exponential curve using home-written Maple code.

Cross-Relaxation Time Measurements
Cross-relaxation rates for 7 Li-1 H and 19 F-1 H for each sample were measured at 20 • C using a 1 H-detected HOESY pulse sequence, following an identical procedure to that previously reported (Martin et al., 2018a) and with the same NMR spectrometer and probe used for the longitudinal relaxation measurements. Mixing times for the HOESY build-up curves were varied over 22 values ranging from 3.7 ms to 4.5 s and FIGURE 1 | Structures of (A) the C 3 mpyr + cation and (B) the DME molecule, with the eight groups of distinct protons labeled, (C) the FSI − anion, and (D) the 1 H NMR spectrum of LiDME-0.5 with peak assignment.
440 scans were recorded for each mixing time. A recycle delay between 5.0 and 6.0 s was used to ensure quantitative measurements. The HOESY build-up curves were fitted to the following expression with σ (the cross-relaxation rate) as the sole variable: where M(τ ) is the observed 1 H signal, M 0 is the thermal equilibrium magnetization of the 7 Li or 19 F spin, R I/S are the longitudinal relaxation rates of the two spins, D I/S are the diffusion coefficients of the ions containing the spins, γ I/S are the nuclear gyromagnetic ratios, g is the gradient pulse strength, δ the gradient pulse length, the delay time between the first and last gradient pulses, and: A complete description of the signal normalization and HOESY curve fitting procedure can be found in our previous publication (Martin et al., 2018b).

RESULTS AND DISCUSSION 7 Li-1 H Cross-Relaxation Rates
The single-pulse excitation 1 H NMR spectrum of the LiDME-0.5 sample is shown in Figure 1D with the peaks assigned according to the structures in Figures 1A,B. All of the peaks are wellresolved with the exception of sites 5 (the ring CH 2 protons closest to the nitrogen on the C 3 mpyr cation) and 7 (the DME methyl protons), which show partial overlap. In the 1 H-7 Li HOESY experiments, the observed signals were much weaker, and the HOESY signal from site 5 was found to be generally less intense than that from site 7 in most cases. This made the extraction of the 1 H-7 Li cross-relaxation rates of the protons of site 5 from the DME-containing samples extremely difficult. The peak at around 3.8 ppm in the HOESY spectra was therefore fitted to a single cross-relaxation rate assumed to reflect only that of the site 7 protons on the DME molecule. For the 1 H-19 F HOESY experiments, higher signal intensities were obtained and these peaks were therefore able to be deconvoluted with separate HOESY build-up curves fitted for sites 5 and 7. The 7 Li-1 H cross-relaxation rates (×10 −4 s −1 ) measured at 20 • C from all three samples are presented in Table 3, and example fitted HOESY build-up curves from which these parameters were extracted are shown in Figure 2. The larger the value of the cross-relaxation rate, the stronger the NOE interaction between the two spins and this is typically interpreted as a closer average distance between the two nuclei. However, it should be noted here that variations in dynamics (specifically of the internuclear vector) will also play a role in determining the σ-value.
Considering first the interaction between the Li + and the C 3 mpyr cation in sample LiMDE-0, the largest NOE is observed between the 7 Li and the methyl proton on the carbon directly attached to the nitrogen site. This is consistent with a previous DNP NMR and molecular dynamics (MD) study of a glassystate LiFSI-containing C 3 mpyrFSI sample in which measured C-Li distances suggested that the lithium cations are preferentially located on the side of the pyrrolidinium ring with this methyl   Uncertainties in these values are estimated at ± 10%. * Values for peak 5 could not be obtained from the samples containing DME due to low intensity and overlap with peak 7.
group (Sani et al., 2018). The smallest NOEs are observed for the protons on the propyl chain. Interestingly, upon addition of DME to this system, the NOEs between the lithium and the various sites on the C 3 mpyr cation remain roughly the same, with the notable exception of the terminal methyl protons on the propyl chain (site 1), which show a significant increase in the cross relaxation rate.
With regards to the interactions between the Li + cations and DME molecules, the σ-values in Table 3 show Li-DME NOEs that are larger by an order of magnitude or more than the Li-C 3 mpyr NOEs. This could be expected given the solvation of the Li + cations by the ether oxygen groups of the DME molecules, as well as the fact that DME is a neutral molecule and hence will not show any charge repulsion from the Li. LiDME-1 shows the highest Li-DME NOEs and the lowest Li-C 3 mpyr NOEs (with the exception of the propyl CH 3 group). It can be noted here that the diffusion coefficients of the Li + and FSI − ions are very similar in LiDME-0 ( Table 2) indicating strong solvation of the lithium cations by the anions, but with the addition of DME the diffusion coefficient FIGURE 2 | Fits of the 1 H-7 Li HOESY build-up curves from all of the 1 H peaks of sample LiDME-0.5 at 20 • C (site numbers refer to the structures/peaks shown in Figure 1). of the Li + is closer to that of the DME molecules while the FSI anions diffuse faster. In combination with the HOESY results, this suggests that the DME solvates the Li + by displacing (at least partially) the FSI − ions. The solvation of the Li + by the DME molecules and the displacement of the FSI anions around the Li + means that there will be more repulsion between the Li and the C 3 mpyr (e.g., a Li-DME-FSI cluster will have a net neutral charge compared to an Li(FSI) − 2 cluster that will be more attracted to the C 3 mpyr), and this can explain the reduced Li-C 3 mpyr NOEs as the DME concentration increases. It was not possible to obtain reliable 7 Li-1 H cross-relaxation rates for the CH 2 on the IL cation ring near N + because it strongly overlaps with the terminal DME proton peak. Finally, the Li-DME NOEs are stronger for the CH 2 protons on the DME molecule than for the CH 3 protons. This may reflect the conformation of the DME molecule when solvating the lithium, however, the different dynamics of the two groups of protons (e.g., the rotation of the methyl group) may also play a role in the different NOE efficiencies.

F-1 H Cross-Relaxation Rates
The cross-relaxation rates extracted from 19 F-1 H HOESY buildup curves are displayed in Table 4, including values obtained from the neat C 3 mpyrFSI IL (without LiFSI) for comparison. These values give insights into the interactions between the FSI anions and the C 3 mpyr + cation/DME species, and once again we make the assumption that larger NOEs reflect closer average Neat IL LiDME-0 LiDME-0.5 LiDME-1 distances between these anions and the different proton sites.
In the neat IL, the 19 F-1 H cross-relaxation rate is largest for the N-CH 3 group and the values decrease as the proton sites become more distant from this central nitrogen site. This can be interpreted in terms of the localized positive charge at this region of the cation which acts to attract the FSI anions. Upon addition of the LiFSI salt (i.e., sample LiDME-0), the σ-values increase for all proton sites. Rather than a net increase in the FSI-C 3 mpyr interaction strength, we interpret this as being due to a change in the overall dynamics of the system, as it is well known that adding lithium salts to ionic liquids results in a significant increase in viscosity (Bose et al., 2018). This phenomenon is confirmed by the decrease in diffusion coefficients of both ions by around an order of magnitude upon the addition of LiFSI salt in LiDME-0 as shown in Table 2. The LiDME-0 sample also shows the largest 19 F-1 H NOE for the N-CH 3 protons which is again consistent with the picture established by DNP NMR and MD simulations (Sani et al., 2018) of a negatively charged Li(FSI) n cluster existing close to this part of the cation. Upon addition of DME to this system, the σ-values can be seen to decrease overall, and once again this can be attributed at least in part to a reduction in the sample viscosity and an increase in the overall dynamics of the liquid (reflected also by an increase in the diffusion coefficients in Table 2). With regards to the FSI-C 3 mpyr interactions, the overall trends remain similar to the LiDME-0 sample with the strongest interaction still being between the FSI and the N-CH 3 group. The fact that the interaction between the FSI anions and the terminal methyl group on the propyl chain of the cation does not appear to increase upon DME addition suggests that the increased interactions between the lithium cations and this particular group may be due to the presence of Li(DME) n clusters. However, the much larger 19 F-1 H NOE values observed for the protons on the DME molecules strongly indicate the presence of lithium cations that are simultaneously solvated by both FSI and DME. MD simulations are currently being carried out on this system and should help with the further interpretation of these NOE values.

CONCLUSIONS
We have measured 7 Li-1 H and 19 F-1 H Overhauser crossrelaxation rates from several pyrrolidinium-based ionic liquid electrolytes containing lithium salt and the 1,2-dimethoxy ethane additive in different ratios. These parameters allow a quantitative comparison of the interactions of the Li + and FSI − ions with the IL cation and DME molecules in the different electrolytes. Both FSI-C 3 mpyr and Li-C 3 mpyr NOEs decrease overall as the DME content increases. The Li-DME and FSI-DME NOEs were found to be the largest, indicating the solvation of Li + by both of these species. Strong FSI-DME NOEs indicate simultaneous solvation of a Li + by both DME and FSI. Li-C 3 mpyr and FSI-C 3 mpyr interactions were strongest for the N-CH 3 protons which is consistent with a previous study that concluded that Li(FSI) − n clusters preferentially locate close to that region of the cation (Sani et al., 2018). The measured diffusion coefficients assisted greatly in interpreting the NOE values, providing additional evidence for ion associations while also reflecting composition-dependent changes in the overall dynamics that affect the sizes of the measured NOE values. These results help to quantify the ionic interactions and the role that the additive 1,2dimethoxy ethane plays in the ionic liquid N-methyl-N-propyl pyrrolidinium bis-fluorosulfonylimide. Further measurements as a function of temperature, along with the correlation to molecular dynamics simulation data, will allow a more detailed interpretation of these cross-relaxation rates, and such work is currently underway.

AUTHOR CONTRIBUTIONS
DG carried out the experiments, data analysis and interpretation, and wrote the paper. P-AM, UP, MD, MF, and LO contributed to the interpretation of the results and edited the paper.