Sustainable and Printable Nanocellulose-Based Ionogels as Gel Polymer Electrolytes for Supercapacitors

A new gel polymer electrolyte (GPE) based supercapacitor with an ionic conductivity up to 0.32–0.94 mS cm−2 has been synthesized from a mixture of an ionic liquid (IL) with nanocellulose (NC). The new NC-ionogel was prepared by combining the IL 1-ethyl-3-methylimidazolium dimethyl phosphate (EMIMP) with carboxymethylated cellulose nanofibers (CNFc) at different ratios (CNFc ratio from 1 to 4). The addition of CNFc improved the ionogel properties to become easily printable onto the electrode surface. The new GPE based supercapacitor cell showed good electrochemical performance with specific capacitance of 160 F g−1 and an equivalent series resistance (ESR) of 10.2 Ω cm−2 at a current density of 1 mA cm−2. The accessibility to the full capacitance of the device is demonstrated after the addition of CNFc in EMIMP compared to the pristine EMIMP (99 F g−1 and 14.7 Ω cm−2).


Introduction
From all the different types of energy storage systems (ESS), the interest in supercapacitors (SCs) has grown thanks to their energy storage mechanism by electrical double layer, which increases cycle life and allows the use of safer ionic conductors and electrode materials than in traditional batteries [1]. Their promising properties such as high-power density, fast charge-discharge rate and long-cycle stability, have turned SCs into a potentially useful tool for application in flexible and wearable electronics [2][3][4]. In this sense, the use of these quasi-solid electrolytes could play an important role in the production of printed supercapacitors providing flexibility, transparency, stability, and safety to the devices [5].
Different types of electrolytes can be used for ESS, including liquid, solid state electrolytes (SSEs) [6] and gel polymer electrolytes (GPEs). GPEs consist basically of a liquid electrolyte, that can be an ionic liquid (IL), or an ionic salt trapped inside a threedimensional (3D) polymer network [7,8]. ILs, generally a room temperature molten salt, chemical characterizations (capacity, ionic conductivity and cyclability) were performed to investigate its properties and understand the functionality of the novel nanocellulose-based ionogel electrolytes.

Materials
Carboxymethylated cellulose nanofibers were produced by RISE Research Institutes of Sweden from a commercial sulfite softwood pulp (Domsjö Dissolving plus; Domsjö Fabriker AB, Domsjö, Sweden). Ethanol (Rectapur ® ) was purchased from VWR (VWR AB International, Spånga, Sweden). Monochloroacetic acid (99%, ACS reagent, ClCH 2 COOH), acetic acid (ACS reagent), iso-propanol (ACS reagent), sodium hydroxide (ACS reagent), sodium hydrogen carbonate (ACS reagent) and methanol (ACS reagent) were purchased from Sigma Aldrich (Stockholm, Sweden). 1-ethyl-3-methylimidazolium dimethyl phosphate (>98% grade) (EMIMP) was purchased at Iolitec (Heilbronn, Germany) and used as received. Active Carbon YP50F (Kuraray, Kurashiki, Japan), conductive carbon TIMCAL C 65 (Cabot), 1-Methyl-2-pyrrolidone (NMP, Scharlau, Sant Feliu del Llobregat, Spain) and 5130 PVDF Solef (Solvay, Barcelona, Spain) were used in the formulation of the active layers of the supercapacitors. The fabricated electrodes were assembled in CR2032 coin cells from MTI by using cellulose separator Kodoshi (Nippon Kodoshi Corporation, Kochi, Japan) and 1-ethyl-3-methylimidazolium dimethyl phosphate (EMIMP, >98%, Iolitec) or EMIMP:CNFc mixtures (96:4 to 99:1 in ratio) as electrolyte. 5 samples electrolytes are reported including EMIMP pure electrolyte and 4 mixtures adding 1%, 2%, 3% and 4% of CNF in EMIMP ionic liquid. It has not been possible to produced jelly electrolyte with higher rate of CNF due to the loss of the mechanical properties not allowing to print and process the component. The preparation of the electrode is identic for each symmetric coin cell and described in Section 2.3. Each electrode is prepared by coating of the jelly electrolyte directly on the top of the dried carbon electrode and then the CR2032 coin cell is assembled with staking 2 symmetric electrode separated by Kodoshi nanocellulose separator (thickness 50 µm) without further addition of electrolyte or solvent.

Carboxymethylated Cellulose Nanofibers (CNFc) Fabrication
CNFc fabrication was done according to the general methodology described in the literature [26]. Most commonly, CNFc is prepared with a total charge of~600 µeq/g, corresponding to a carboxmethyl degree of substitution (DS) of 0.1 [27]. In this work, a higher charge DS 0.3 CNFc was produced as has been recently described [28]. Generally, resulted pulp is solvent exchanged to ethanol followed by carboxymethylation reaction, which begins by impregnating the pulp with a solution of monochloric acetic acid in isopropanol, followed by transfer of the pulp to a heated methanolic solution of NaOH in isopropanol and reaction under reflux for 1 h. Compared to the reaction conditions for DS 0.1, to achieve DS 0.3, the dosage of monochloroacetic acid is increased 4.4-fold and NaOH 2.7-fold. After reaction, the carboxymethylated pulp is washed with NaHCO 3 to neutralize, and finally washed with water to remove excess reagents, until the wash has a conductivity below 10 µS/cm. To generate CNFc, the carboxymethylated pulp is microfluidized (Microfluidizer M-110EH, Microfluidics Corp., Newton, MA, USA) at 1700 bar (1 pass) to give a suspension of DS 0.3 CNFc at 1 wt%. The CNFc was then subjected to lyophilization and morturation before use in Leitat (Barcelona, Spain).

Supercapacitor Preparation
The active ink used in the supercapacitor application was fabricated by mixing NMP, active carbon YP50F, conductive carbon TiMCALC65 and binder PVDF at dry weight ratio of 80:10:10 by using a Dipermat (VMA-GETZMANN GmbH, Reischsof, Germany) with a 3D printed head reproducing the shape of an industrial planetary mixer but adapted for lower volume of ink. The prepared ink was printed onto a 20 μm aluminum foil by bare coating (gap of 100 μm) and dried for 2 h at 120 °C in vacuum obtaining a dry coating of 50 μm thickness and active material loading of 1.8 mg cm −2 . For the CNFc-ionogel printing, heating to 80-100 °C allowed direct printing using a Doctor Blade technique (gap of 150 μm bar) onto a carbon electrode surface as shown in Figure 2.

Characterization
X-ray diffraction (XRD) patterns were recorded for each sample on a diffractometer (Malvern PANalytical X'pert Pro MPD, X-Ray Diffraction Facility, Institut Català de Nanociència i Nanotecnologia) at room temperature from 10 to 50 2Ɵ. Fourier transform infrared spectroscopy were recorded with a Shimadzu Iraffinity-1S spectrometer in ATR mode. The specimens were measured directly with a scan range from 400 to 4000 cm −1 . The thermogravimetric analysis (TGA) of the samples (5-10 mg) were carried out by using a TGA/DSC (Mettler-Toledo, University of Alicante, Alicante, Spain), equipped with autosampler and gas-controlled unit, from room temperature to 900 °C with heating rate of 20 °C/min under N2 conditions. Scanning Electron Microscopy (SEM) were performed in

Supercapacitor Preparation
The active ink used in the supercapacitor application was fabricated by mixing NMP, active carbon YP50F, conductive carbon TiMCALC65 and binder PVDF at dry weight ratio of 80:10:10 by using a Dipermat (VMA-GETZMANN GmbH, Reischsof, Germany) with a 3D printed head reproducing the shape of an industrial planetary mixer but adapted for lower volume of ink. The prepared ink was printed onto a 20 µm aluminum foil by bare coating (gap of 100 µm) and dried for 2 h at 120 • C in vacuum obtaining a dry coating of 50 µm thickness and active material loading of 1.8 mg cm −2 . For the CNFc-ionogel printing, heating to 80-100 • C allowed direct printing using a Doctor Blade technique (gap of 150 µm bar) onto a carbon electrode surface as shown in Figure 2. different ratios of EMIMP (p/p, EMIMP/CNFc: 99:1, 98:2, 97:3 and 96:4, referred as EM-IMP/CNFc 99:1, EMIMP/CNFc 98:2, EMIMP/CNFc 97:3 and EMIMP/CNFc 96:4) with heating (100 °C) for 16 h to obtain a yellowish viscous gel (final appearance pictures are depicted in Figure S1 from Supplementary Materials). The as prepared ionogel could be used without any additional step in ambient air, showing long-term stability since no destabilization phenomena, like phase-separation, have been observed after approx. 2 years.

Supercapacitor Preparation
The active ink used in the supercapacitor application was fabricated by mixing NMP, active carbon YP50F, conductive carbon TiMCALC65 and binder PVDF at dry weight ratio of 80:10:10 by using a Dipermat (VMA-GETZMANN GmbH, Reischsof, Germany) with a 3D printed head reproducing the shape of an industrial planetary mixer but adapted for lower volume of ink. The prepared ink was printed onto a 20 μm aluminum foil by bare coating (gap of 100 μm) and dried for 2 h at 120 °C in vacuum obtaining a dry coating of 50 μm thickness and active material loading of 1.8 mg cm −2 . For the CNFc-ionogel printing, heating to 80-100 °C allowed direct printing using a Doctor Blade technique (gap of 150 μm bar) onto a carbon electrode surface as shown in Figure 2.

Characterization
X-ray diffraction (XRD) patterns were recorded for each sample on a diffractometer (Malvern PANalytical X'pert Pro MPD, X-Ray Diffraction Facility, Institut Català de Nanociència i Nanotecnologia) at room temperature from 10 to 50 2Ɵ. Fourier transform infrared spectroscopy were recorded with a Shimadzu Iraffinity-1S spectrometer in ATR mode. The specimens were measured directly with a scan range from 400 to 4000 cm −1 . The thermogravimetric analysis (TGA) of the samples (5-10 mg) were carried out by using a TGA/DSC (Mettler-Toledo, University of Alicante, Alicante, Spain), equipped with autosampler and gas-controlled unit, from room temperature to 900 °C with heating rate of 20 °C/min under N2 conditions. Scanning Electron Microscopy (SEM) were performed in

Characterization
X-ray diffraction (XRD) patterns were recorded for each sample on a diffractometer (Malvern PANalytical X'pert Pro MPD, X-ray Diffraction Facility, Institut Català de Nanociència i Nanotecnologia) at room temperature from 10 to 50 2θ. Fourier transform infrared spectroscopy were recorded with a Shimadzu Iraffinity-1S spectrometer in ATR mode. The specimens were measured directly with a scan range from 400 to 4000 cm −1 . The thermogravimetric analysis (TGA) of the samples (5-10 mg) were carried out by using a TGA/DSC (Mettler-Toledo, University of Alicante, Alicante, Spain), equipped with autosampler and gas-controlled unit, from room temperature to 900 • C with heating rate of 20 • C/min under N 2 conditions. Scanning Electron Microscopy (SEM) were performed in a Zeiss Merlin FE-SEM (Carl Zeiss Microscopy, Autonomous University of Barcelona, Bellaterra, Spain) with samples mounted directly on SEM stubs with adhesive carbon tape and metalized with Au thin film (E5000 Sputter Coater. Polaroid Equipment Limited. Wat-ford). Images were recorded at different steps of magnification. The rheologic properties of EMIMP and all fabricated electrolytes were analyzed by using the rheometer Bohlin CVO 100-901 in the share rate range from 0.0517 to 888 s −1 .
The electrochemical testing was performed with a multichannel potentiostat VMP3 (Bio-Logic) run with the EC-Lab software. Cyclic voltammetry was performed between 0 < V < 2.0, and with a scan rate between 10 < mV s −1 < 50, to investigate side reactions, evaluate the potential window of stability and the device behavior. Galvanostatic measurements at a constant areal capacity of mA·cm −2 and between 0 < V < 2.0 V) were used to calculate the specific capacity (mAh·g −1 of AC, energy (unit), power (unit), and equivalent series resistance (ESR) is calculated from the ohmic drop between of the developed supercapacitors. Electrochemical impedance spectrum measurements were carried out with frequency range of 50 mHz to 1 MHz. All these analyses were carried out by VMP3 from Biologic. Finally, a long cycling at 1 mA·s −1 for 10,000 cycles were carried out by Neware cycler. The electrochemical CR2032 coin-cell type used for the tests is built from the stack of two symmetric electrodes separated by cellulose membrane separators and with 40 µL of EMIMP or 150 µm (wet) bar-coated EMIMP:CNFc gel as electrolyte ensuring it excess.
The cell capacitance (C cell ) of a symmetrical capacitor can be calculated according to Equation (1) [29]: where C e = C + = C − are the specific capacitance of the two electrodes in mAh·g −1 of activated carbon, and the cell capacitance of a supercapacitor can be calculated from the charge-discharge experiment Equation (2): This can be then used, when combined with C cell equation to obtain the electrode capacity, and subsequently obtain information on the capacity of an investigated material Equation (3): where I (A) is the discharge current, m (g) the total mass of active materials of two electrodes and dV/dt the slope of the discharge curve between 80% and 40% of the cut-off voltage (V). From the total specific cell capacitance of two-electrode system (C cell ), the maximum energy (E max. ) and maximum power (P max. ) of a supercapacitor cell can be calculated according to Equations (4) and (5): where, V (V) is the cell voltage minus the ohmic drop, and ESR (Ω) is the total equivalent series resistance of the supercapacitor. From the EIS, the ionic conductivity of the different gels has been calculated using the following formula Equation (6): where σ is the ionic conductivity (S m −1 ), l (m) is the thickness of the membrane, A (m 2 ) is the contact area between the electrode and the electrolyte and R 1 (Ω) is the ionic resistance at high frequency corresponding to the first interception of the circle in the Nyquist plots.

Results and Discussion
An extensive characterization of the prepared CNFc-ionogels have been realized to study their composition, homogeneity, thermal and rheological properties. XRD patterns for each IL/CNFc ratio, pristine CNFc (lyophilized) and EMIMP are shown in Figure 3a. Broad signals are obtained in all cases, obtaining major intensities at 12 • and 25 • for all the EMIMP/CNFc mixtures, which are related to the ionic liquid EMIMP (triangle signals), with a minor signal present as a broad band at 15 • corresponding to (002) nanocellulose crystalline plane (star signals). The major component of ionic liquid ensures a good ionic conductivity and a good gel flexibility due to their mainly amorphous nature. Fourier transform infrared spectroscopy (FTIR) of each CNFc-ionogel were also recorded and are depicted in Figure 3b, where the IR profiles of pure CNFc and EMIMP are also added. As in case of previous XRD studies, the EMIMP has the major intensity, showing its typical bands of -CH from the imidazolium ring (3060 and 1570 cm −1 ), the -CH 2 st and -CH 3 st (2920-2970 cm −1 ), -CH 3 bend (1450 cm −1 ) and -CH 2 bend (780 cm −1 ) from the ethyl chains, the -C-N st at 1240 cm −1 and finally, the -P-O-C (930-1050 cm −1 ) and -P=O (1090-1200 cm −1 ) vibrations of the diethylphosphate group. CNFc (highlighted in blue) shows an -OH band at 3400 cm −1 , which increases with the higher ratio of CNFc in the mixtures, a -C=O peak related with the carboxylated groups at 1600 cm −1 , which appears slightly shifted to 1650 cm −1 for all ionogel mixtures, and a -C-O-C streaching at 1000 cm −1 from CNFc glucose units, which can not be appreciated in the mixtures due to overlaping with EMIMP P=O and P-O-C signals.

Results and Discussion
An extensive characterization of the prepared CNFc-ionogels have been realized to study their composition, homogeneity, thermal and rheological properties. XRD patterns for each IL/CNFc ratio, pristine CNFc (lyophilized) and EMIMP are shown in Figure 3a. Broad signals are obtained in all cases, obtaining major intensities at 12° and 25° for all the EMIMP/CNFc mixtures, which are related to the ionic liquid EMIMP (triangle signals), with a minor signal present as a broad band at 15° corresponding to (002) nanocellulose crystalline plane (star signals). The major component of ionic liquid ensures a good ionic conductivity and a good gel flexibility due to their mainly amorphous nature. Fourier transform infrared spectroscopy (FTIR) of each CNFc-ionogel were also recorded and are depicted in Figure 3b, where the IR profiles of pure CNFc and EMIMP are also added. As in case of previous XRD studies, the EMIMP has the major intensity, showing its typical bands of -CH from the imidazolium ring (3060 and 1570 cm −1 ), the -CH2 st and -CH3 st (2920-2970 cm −1 ), -CH3 bend (1450 cm −1 ) and -CH2 bend (780 cm −1 ) from the ethyl chains, the -C-N st at 1240 cm −1 and finally, the -P-O-C (930-1050 cm −1 ) and -P=O (1090-1200 cm −1 ) vibrations of the diethylphosphate group. CNFc (highlighted in blue) shows an -OH band at 3400 cm −1 , which increases with the higher ratio of CNFc in the mixtures, a -C=O peak related with the carboxylated groups at 1600 cm −1 , which appears slightly shifted to 1650 cm −1 for all ionogel mixtures, and a -C-O-C streaching at 1000 cm −1 from CNFc glucose units, which can not be appreciated in the mixtures due to overlaping with EMIMP P=O and P-O-C signals. Thermal properties (TGA) were studied for each CNFc-ionogel and are depicted in Figure 3c. The TGA profiles shown in all cases a two-step weight loss at temperatures below 100 • C and above 250 • C, which can be reasonably attributed to evaporation of adsorbed water (less than a 15% wt) and the degradation of CNFc and the EMIMP, respectively. This data shows the high thermal stability of the new gel electrolyte. DSC experiments have been also studied to determine the degradation temperatures for each ionogel in comparison to pure EMIMP ( Figure S2 and Table S1, Supplementary Materials), obtaining similar results between them (290 • C for EMIMP, to 293 • C for EMIMP/CNFc 96:4). The addition of smalls amounts of nanocellulose has a minor affection on the decomposition temperatures, improving them slightly. ILs are mostly Newtonian fluids and more viscous than most common molecular solvents [30]. Their ionogels show similar viscoelastic properties under stress, providing a higher-quality electrode/electrolyte contact, in comparison with other rigid solid-state electrolytes [10,31]. In the case of this work, viscosity measurements ( Figure S3, Supplementary Materials) reveal that almost all the prepared NC-ionogels have similar viscosity characteristics, except for the mixture EMIMP/CNFc 96:4, that reveals a thixotropic behavior, typical of NC. Regarding the apparent viscosity dependence of the applied shear stress (Figure 3d), all NC-ionogels, and CNFc (1% wt, gel), show a reduction in the viscosity, which is more pronounced in the ionogels with a higher CNFc content (EMIMP/CNFc 96:4 and EMIMP/CNFc 97:3), with no changes in case of pure EMIMP. This characteristic, named shear-thinning [32], is common of pseudoplastic fluids and consents on a more effective particles flow through the liquid phase, reducing the viscosity of the solution. It is intrinsic of the nanocellulose and its derivatives [33,34], and it is responsible for the good printability of the presented ionogels, while the desired mechanical properties were maintained.
High resolution scanning electronic microscopy (FE-SEM) was done for the ionogel EMIMP/CNFc 96:4 to study its morphology, showing a smooth and homogeneous surface. Moreover, an EDS mapping of nitrogen, phosphorous and sodium elements that are present in the ionic liquid (N and P) and the nanocellulose (Na) was also performed to confirm the homogeneity of the GPE. (See Figure 4).  In the same way, to use these NC-based ionogels as electrolytes for supercapacitors, their electrochemical performance and properties were first evaluated. Symmetric cells were assembled in CR2032 coin cell to test the electrochemical performance of the EMIMP electrolyte, alone and with the different CNFc ratios, using a cellulosic Kodoshi separator. In the same way, to use these NC-based ionogels as electrolytes for supercapacitors, their electrochemical performance and properties were first evaluated. Symmetric cells were assembled in CR2032 coin cell to test the electrochemical performance of the EMIMP electrolyte, alone and with the different CNFc ratios, using a cellulosic Kodoshi separator. The cyclic voltammogram at 10 mV·s −1 (See Figure 5a) of the supercapacitor with EMIMP electrolyte shows a near rectangle shape, approximated to the ideal situation of double layer capacitor (or supercapacitor). There is no visible redox peak from Faradaic current over the potential region. However, a peak onset appears as 2 V is approached. At high scan rates (See Figure 5b), the shape of the CV curve was distorted due to an increase of the internal resistance and the peak onset near 2 V was shifted to higher voltage, due to the increase in cell resistance at higher scan rate. In order to clarify if the peak onset observed in the cyclic voltammetry was caused by the electrolyte degradation or because of its own behavior (as has been seen before in different IL electrolytes), [35] a floating test was carried out. [36] The supercapacitors with the different studied electrolytes were charged at 2V in a stepwise ramp and the voltage was maintained for a period of 24 h, as shown Figure S4 from Supplementary Materials. The leakage current obtained at each time was recorded and could be used to directly assess the importance of the electrochemical degradation of the electrolytes in the fully charged state. The results showed that the leakage current of the EMIMP alone was 6 times higher than the one obtained by adding CNFc (i.e., 120 μA vs 20 μA). The data obtained for the different EMIMP/CNFc electrolytes were low enough to determine that the peak onset showed in the cyclic voltammogram was not caused by electrolyte degradation. [37] Moreover, these results demonstrated that the addition of CNFc increase the EMIMP voltage resistance. The capacity and resistance of the supercapacitors with EMIMP and with the gel electrolytes (by the introduction of different rates of CNFc in the EMIMP) were evaluated by galvanostatic charge-discharge tests. All these parameters were calculated using the equations described in Section 2 and the results are showed in Table S1 in Supplementary Materials and Figure 6. In the case of capacities calculated from cyclic voltammetry, as shown Table S4 from Supplementary Materials, the capacitance tendency in all cases is decreasing at higher scan rate. Higher specific capacity is measured after addition of CNFc, which is online with galvanostatic charge-discharge results. The supercapacitors with CNFc in the electrolyte exhibit a smaller ESR in all intensity ranges than the supercapacitors with only EMIMP (4.0-10.2 Ω for different rates of EMIMP:CNFc and 14.7 Ω for EMIMP alone). These results also indicate that EMIMP/CNFc gel electrolytes (in all the studied ranges) have higher ionic conductivity than EMIMP alone, and this has been confirmed by the electrochemical impedance spectroscopy showed in Figure 7. In addition, When CNFc is added to the EMIMP ionic liquid, the capacity of the device increases with the addition of CNFc, the peak onset near 2 V in the oxidation zone becomes more evident and a new one between 0 and 0.5 V in the reduction zone appears (Figure 5a). At high scan rates (Figure 5b), the shape of the CV curve was distorted due to an increase of the internal resistance but to lesser extent than with EMIMP alone and the peaks onset near 2 V and between 0 and 0.5 V disappear. However, no clear differences in capacity and resistance between the electrolytes with different ratios of EMIMP/CNFc can be appreciated from the cyclic voltammetry studies.
In order to clarify if the peak onset observed in the cyclic voltammetry was caused by the electrolyte degradation or because of its own behavior (as has been seen before in different IL electrolytes) [35], a floating test was carried out [36]. The supercapacitors with the different studied electrolytes were charged at 2V in a stepwise ramp and the voltage was maintained for a period of 24 h, as shown Figure S4 from Supplementary Materials. The leakage current obtained at each time was recorded and could be used to directly assess the importance of the electrochemical degradation of the electrolytes in the fully charged state. The results showed that the leakage current of the EMIMP alone was 6 times higher than the one obtained by adding CNFc (i.e., 120 µA vs 20 µA). The data obtained for the different EMIMP/CNFc electrolytes were low enough to determine that the peak onset showed in the cyclic voltammogram was not caused by electrolyte degradation [37]. Moreover, these results demonstrated that the addition of CNFc increase the EMIMP voltage resistance. The capacity and resistance of the supercapacitors with EMIMP and with the gel electrolytes (by the introduction of different rates of CNFc in the EMIMP) were evaluated by galvanostatic charge-discharge tests. All these parameters were calculated using the equations described in Section 2 and the results are showed in Table S1 in Supplementary Materials and Figure 6. In the case of capacities calculated from cyclic voltammetry, as shown Table S4 from Supplementary Materials, the capacitance tendency in all cases is decreasing at higher scan rate. Higher specific capacity is measured after addition of CNFc, which is online with galvanostatic charge-discharge results. The supercapacitors with CNFc in the electrolyte exhibit a smaller ESR in all intensity ranges than the supercapacitors with only EMIMP (4.0-10.2 Ω for different rates of EMIMP:CNFc and 14.7 Ω for EMIMP alone). These results also indicate that EMIMP/CNFc gel electrolytes (in all the studied ranges) have higher ionic conductivity than EMIMP alone, and this has been confirmed by the electrochemical impedance spectroscopy showed in Figure 7. In addition, according to Equations (2) and (3), the capacity for the electrolyte EMIMP and EMIMP/CNFc with different ratios were calculated to be 99 and 132-160 F g −1 of YP50F active material at current density of 1 mA cm −2 . The introduction of CNF in the electrolyte significantly increases the capacity of the electrode, and this increase on capacity and decrease on ESR is also significant between the four gel electrolytes with different CNFc ratios. Maximum of capacity is reached for EMIMP/CNFc 97:3 ionogel but presenting similar capacity than the EMIMP/CNFc 96:4 ratio. However, the ESR of EMIMP/CNFc 96:4 is clearly lower than the other gel electrolytes that present similar values. In addition, we previously demonstrated that this increase on capacity does not come from the electrolyte (or CNFc) degradation by the floating test, rather this increment on capacity suggests that part of the CNFc present in the electrolyte is displaced to the electrode surface (due to his high polarity), creating an extra porous structure where more ions are allocated during the charge process, and therefore, increasing the capacity of the device. It has not been tested above 96:4 due to the poor mechanical of gel above this range. (or CNFc) degradation by the floating test, rather this increment on capacity suggests that part of the CNFc present in the electrolyte is displaced to the electrode surface (due to his high polarity), creating an extra porous structure where more ions are allocated during the charge process, and therefore, increasing the capacity of the device. It has not been tested above 96:4 due to the poor mechanical of gel above this range. The capacity variation with the cycles is shown in Figure S5. The supercapacitors with EMIMP show a good capacity retention of 80% after 10,000 cycles at 1 mA cm −2 between 0 < V <2. The supercapacitors with CNFc had a quick capacity decay in the first 2000 cycles (~0.020 mAh g −1 cycle −1 of capacity decay) and then present a moderate degradation (~0.001 mAh g −1 cycle −1 of capacity decay) reaching the 10,000 cycles with a capacity retention below 50%. It is possible that during the cycles, the CNFc cannot be retained in the matrix of the electrolyte and that creates a larger SEI in the first cycle, starting to cover part of the carbon porous and therefore limiting the capacity of the electrodes.
To rule out a possible dissolution of CNFc to carboxymethyl cellulose (CMC) units, an ionogel with a 96:4 ratio EMIMP/CMC was prepared and assessed. Hence, in the case of use CMC instead CNFc, the electrochemical performance of these supercapacitors,  The capacity variation with the cycles is shown in Figure S5. The supercapacitors with EMIMP show a good capacity retention of 80% after 10,000 cycles at 1 mA cm −2 between 0 < V <2. The supercapacitors with CNFc had a quick capacity decay in the first 2000 cycles (~0.020 mAh g −1 cycle −1 of capacity decay) and then present a moderate degradation (~0.001 mAh g −1 cycle −1 of capacity decay) reaching the 10,000 cycles with a capacity retention below 50%. It is possible that during the cycles, the CNFc cannot be retained in the matrix of the electrolyte and that creates a larger SEI in the first cycle, starting to cover part of the carbon porous and therefore limiting the capacity of the electrodes. Finally, electrochemical impedance spectroscopy (EIS) for the supercapacitors with EMIMP/CNFc gel electrolytes and EMIMP are shown in Figure 7a. The equivalent circuit was made with EC-Lab (Figure 7b). The Nyquist plots show a semi-circle at high frequencies region and a line perpendicular to the X-axis at low frequencies that correspond the ionic diffusion and the double layer formation at the surface of the conductive solid phase. R1 is mostly corresponding to the ionic conductivity in the dielectric, both gel and membrane separator. R2 is referring to the charge transfer both ionic and electronic in the bulk of the electrode. Q1 and Q2 are constant phase element simulating the real behaviour of a pure capacitor. The first thing to be considered is a charge transfer resistance (Rct, that includes the ionic resistance), which could be intensively investigated in the semi-circle loop appearing in the middle frequency region. A decrease in resistances is shown in supercapacitors with CNFc compared to supercapacitors with only EMIMP. This decrease in resistances has been observed above all in the value of the charge transfer resistance (Rct). These results are in line with the ones obtained in the CV and galvanostatic tests. The real and imaginary capacitance versus the frequency is given in supporting in Figure S7 and using the definition from P.L. Taberna and al. [40] The real part of the capacitance shows the variation of the available stored energy with the frequency. The imaginary part of the capacitance represents the losses that occur during charge storage. In these graphs, can observe that the maximum capacitance is reached faster at higher frequencies in the case of the blended electrolyte containing CNFc versus EMIMP electrolyte. No real difference is observed within the different percentage of CNF used in this studied. No degradation of the electrolyte is demonstrated by the floating test, also the same active material and electrode systems have been used for all electrolyte testing, then we assume no intrinsic supplementary capacity is reached by adding CNFc to EMIMP but a faster access to this maximal capacity is demonstrated.
All the values of ionic conductivity (σ, mS cm −1 ) and resistances (R1 and R2, Ω) are included in Table 1. The increase in conductivity with the addition of CNFc was not expected because of the important increment on the electrolyte viscosity (almost 0 to 60 Pa ·s), Figure 3d and Figure S3 from Supplementary Materials). However, as can be seen in Table 1, adding nanocellulose has helped to increase the ionic conductivity of the GPEs, To rule out a possible dissolution of CNFc to carboxymethyl cellulose (CMC) units, an ionogel with a 96:4 ratio EMIMP/CMC was prepared and assessed. Hence, in the case of use CMC instead CNFc, the electrochemical performance of these supercapacitors, shown in Table S3, provides specific capacitance values of 123 F g −1 , being like EMIMP/CNFc 99:1 ionogel values. Moreover, it has been observed the formation of a solid precipitates by mixing CMC to EMIMP, presumably due to a lack of interaction between the CMC and IL because of the diminution of available -OH groups. As it has been mentioned, the presence of these hydroxyl groups plays an important role in the dissolution of cellulose in ILs, due to they are the responsible of the formation of new hydrogen bonds between the NC and the ILs anions [38,39].
Finally, electrochemical impedance spectroscopy (EIS) for the supercapacitors with EMIMP/CNFc gel electrolytes and EMIMP are shown in Figure 7a. The equivalent circuit was made with EC-Lab (Figure 7b). The Nyquist plots show a semi-circle at high frequencies region and a line perpendicular to the X-axis at low frequencies that correspond the ionic diffusion and the double layer formation at the surface of the conductive solid phase. R1 is mostly corresponding to the ionic conductivity in the dielectric, both gel and membrane separator. R2 is referring to the charge transfer both ionic and electronic in the bulk of the electrode. Q1 and Q2 are constant phase element simulating the real behaviour of a pure capacitor. The first thing to be considered is a charge transfer resistance (R ct, that includes the ionic resistance), which could be intensively investigated in the semi-circle loop appearing in the middle frequency region. A decrease in resistances is shown in supercapacitors with CNFc compared to supercapacitors with only EMIMP. This decrease in resistances has been observed above all in the value of the charge transfer resistance (R ct ). These results are in line with the ones obtained in the CV and galvanostatic tests. The real and imaginary capacitance versus the frequency is given in supporting in Figure S7 and using the definition from P.L. Taberna and et al. [40]. The real part of the capacitance shows the variation of the available stored energy with the frequency. The imaginary part of the capacitance represents the losses that occur during charge storage. In these graphs, can observe that the maximum capacitance is reached faster at higher frequencies in the case of the blended electrolyte containing CNFc versus EMIMP electrolyte. No real difference is observed within the different percentage of CNF used in this studied. No degradation of the electrolyte is demonstrated by the floating test, also the same active material and electrode systems have been used for all electrolyte testing, then we assume no intrinsic supplementary capacity is reached by adding CNFc to EMIMP but a faster access to this maximal capacity is demonstrated.
All the values of ionic conductivity (σ, mS cm −1 ) and resistances (R1 and R2, Ω) are included in Table 1. The increase in conductivity with the addition of CNFc was not expected because of the important increment on the electrolyte viscosity (almost 0 to 60 Pa ·s), Figure 3d and Figure S3 from Supplementary Materials). However, as can be seen in Table 1, adding nanocellulose has helped to increase the ionic conductivity of the GPEs, causing a decrease in the internal resistance of the system (0.32-0.94 mS cm −1 for EMIMP/CNFc ionogels against 0.26 mS cm −1 for EMIMP). Plenty of polar functional groups in the system helped to increase the dielectric constant of the gel electrolytes, which is proportional to the capacity assuming all electrode surface area and distances between the two electrodes are the same for each cell assembled. This also increases the interaction between ions and polymer chains and helps to dissociate these ion pairs into separate ions. As a result, more free ions are generated and consequently, the ionic conductivity is enhanced. Moreover, the abundant polar functional groups supplied by cellulose offer an ionic conduction channel by aiding facile ion hopping onto nearby polar groups However no tendency could be measured with the increased ratio of CNFc.

Conclusions
Novel GPEs have been prepared by the addition of different CNFc ratios (from 1 to 4) to EMIMP, showing excellent thermal and mechanical properties. Contrary to the usual liquid electrolytes, these new NC based GPEs not only improve safety issues, but also show a good printability, very interesting for printed electronics applications. Moreover, enhanced electrochemical properties have also been measured. In case of specific capacitance, EMIMP/CNFc ionogels have shown a significant increase (160 F g −1 ) compared to the pristine IL (99 F g −1 ), and its internal resistances have considerable decreased from 14.7 Ω cm −2 to 10.2 Ω cm −2 at a current density of 1 mA cm −2 . Although the viscosity (Figure 3d) of the electrolyte increases after CNF addition, which should promote lower ionic conductivity, all electrochemical measurements show both higher ionic conductivity and capacity with CNF addition. These data confirm that the addition of polar groups coming from nanocellulose play a beneficial role. However, one of the points to be improved in future experiments is to maintain the capacitance over the cycles (50 % decrease in capacitance in the first 2000 cycles). A drastic loss of the capacity is observed for all samples containing CNF after long cycling (Supplementary Materials Figure S5), although the leakage current test showed no electrolyte degradation ( Figure S4). With these results we assume the addition of CNF is beneficial at device level, reducing the initial internal resistance and increasing the quality of the interface between the electrolyte and the electrode. But this effect disappears within the cycling, where samples with and without CNFc show similar capacity after 2.000 cycles A possible passivation of the electrode that it is reducing both capacity and ionic conductivity is not discarded and needs to be further studied. Moreover, investigations about the use of CNFc vs CMC have been also performed, showing that the use of the nanofibrils of cellulose form (CNFc) provides better solubility in EMIMP and higher electrochemical performances. With all this results we can conclude that these new GPEs have shown great potential for applications in energy storage devices.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.