Reduced Graphene Oxide Derived from Low-Grade Coal for High-Performance Flexible Supercapacitors with Ultrahigh Cyclability

Preparation of reduced graphene oxide (RGO) from abundant and inexpensive low-grade coal is regarded as one of the most promising methods for utilizing this resource in a high-value and environmentally sustainable manner. As the main precursor for the fabrication of RGO, graphene oxide (GO) can be extracted from low-grade coal such as lignite, but its size is just in the range of tens to hundreds of nanometers, which limits its practical application. Herein, we demonstrate that large-size RGO sheets can be prepared in large quantities by the pretreatment of lignite using the high temperature–high pressure (HTHP) method. The RGO electrode after the reduction reaction by 50 mM NaBH4 at 105 °C features porosity and high conductivity, which can facilitate high electrochemical reaction efficiency. Thus, we also demonstrate the use of lignite-derived RGO for supercapacitor electrode materials with high performance. The lignite-derived RGO supercapacitor can deliver outstanding volumetric capacitance (30.6 F cm−3), high energy density (4.2 mW h cm−3), excellent flexibility (79.5% retention of the initial capacitance at 180° bending), and a long lifespan (112.3% retention of the initial capacitance after 20,000 cycles). It is believed that the proposed large-size RGO based on reasonable optimization of inferior lignite will offer a new prospect for next-generation energy storage applications.


Introduction
As epidermal and intelligent electronics become more prevalent in our daily lives, there is a high demand for flexible and high-performance energy storage devices [1][2][3]. Flexible supercapacitors are expected to become an emerging focus among many energy storage devices, owing to their outstanding characteristics, including fast charge-discharge rate, high output power density, high cyclic stability, especially in flexibility and easy integration [4][5][6]. Recently, a lot of works have been devoted to the development of flexible high-performance supercapacitors using novel carbon nanomaterials and their composites [7,8]. Among the various carbon nanomaterials, graphene is a significant type of important functional material which has attracted intensive attention for its promising applications in supercapacitor electrode materials, primarily because of its exceptional electrical conductivity, robust mechanical flexibility, notably large specific surface area, and exceptional chemical inertness to the surrounding environment [9]. Several recent studies have shown that graphene has been explored as an outstanding electrode material in symmetric supercapacitors and shows superior specific capacitance and rate performance [10,11]. However, the mass manufacturing of graphene in a low-cost way is a significant barrier to its widespread industrial application.
Over the past few years, many endeavors have been made toward the preparation of graphene. Mechanical exfoliation using the Scotch-tape method and chemical vapor deposition (CVD) can usually produce good-quality graphene, but are inappropriate for synthetic graphite powder (1.0 g) and NaNO 3 (0.5 g) were dispersed in concentrated H 2 SO 4 (23 mL) and placed in an ice bath environment with stirring. Under constant agitation, KMnO 4 (3.0 g) was slowly added to prevent the reaction temperature from rising above 20 • C. Successively, the mixed suspension was transferred to a 35 • C oil bath for 1 h and the reaction temperature was further increased to 85 • C before attenuation and washing by 100 mL of deionized water and 15 mL of H 2 O 2 (30%). The mixture was filtered and then washed with HCl aqueous solution (5%), ethyl alcohol, and deionized water in sequence until the pH value became around 7. After being purified by dialysis for one week, the GO powder was collected after freeze-drying. For comparison, GO (diameter: 0.5-5 µm, thickness: 1-3 nm) was also acquired from Nanjing XFNANO Materials Tech Co., Ltd., located in China.

Reduction of GO and Preparation of Electrodes
To confirm the effect of the NaBH 4 reducing agent, 10 mg of GO powder was immersed in NaBH 4 solutions (20 mL) with varying concentrations (15,50, and 100 mM) at room temperature (RT) and kept for 4 h. The products were labeled as RGO-RT-15, RGO-RT-50, and RGO-RT-100 corresponding to reduction concentrations of 15, 50, and 100 mM, respectively. To deepen the reduction of GO, GO powder was also dipped in aqueous NaBH 4 (50 mM) at 105 • C for 4 h, and the final product was referred to as RGO-105-50. For the comparison of energy storage performance of RGO from commercial GO and lignite-derived GO, the purchased GO was reduced at the same condition of RGO-105-50, and the as-prepared product was named as c-RGO-105-50. To fabricate film electrodes, the post-reduction GO dispersions were vacuum filtered through a porous alumina membrane (Whatman, 20 nm pore size, 47 mm diameter) to make a homogenous thin film, which was then vacuum-dried at 70 • C for 24 h.

Preparation of Supercapacitor Devices
For the three-electrode system, a platinum sheet and Ag/AgCl electrode were used as the counter electrode and reference electrode, respectively. All the electrodes were prepared by mixing 80% as-prepared RGO, 10% PVDF, and 10% carbon black in NMP to form a homogeneous slurry. After that, the slurry was homogeneously blade-coated onto carbon cloth, and the as-prepared electrodes were dried in a vacuum oven at 80 • C for 10 h. The mass of active material on the coating area (1 cm 2 ) was 2-3 mg.
In the assembly of the flexible supercapacitor, the electrode film was transferred on the gold-coated PET film, before porous alumina was etched by 3 M NaOH solution. The electrodes were placed into a vacuum oven for 12 h at 60 • C before the copper tapes were adhered on one edge of each Au-coated PET piece with silver plastic. The following steps were taken to make the PVA/H 3 PO 4 electrolyte [29]: 3 g of PVA was dispersed in 30 mL of deionized water and heated to 85 • C under continuous stirring until a clear solution was obtained. The PVA/H 3 PO 4 electrolyte was then spread on the RGO films. After drying to get rid of any extra water, two pieces of Au-coated PETs were pressed together to create a device with a classic sandwich structure. The electrode volumes of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, RGO-105-50, and c-RGO-105-50 were determined to be around 2.05 × 10 −4 cm −3 , 2.15 × 10 −4 cm −3 , 2.4 × 10 −4 cm −3 , 2.15 × 10 −4 cm −3 , 2.25 × 10 −4 cm −3 , 4.85 × 10 −4 cm −3 , and 5.2 × 10 −4 cm −3 .

Characterization
X-ray diffraction (XRD) was carried out on X'Pert3 diffractometer with a Cu-Ka target. The field emission scanning electron microscopy (Carl Zeiss Ultra 55, Oberkochen, Germany) and transmission electron microscope (JEOL-2100PLUS, Tokyo, Japan) were used to characterize the morphology and nanostructure of the samples, respectively. Raman spectra (HORIBA LabRAM HR Evolution, Kyoto, Japan) were obtained through the use of a 532 nm laser as the light source. The XPS measurements were carried out on an Axis Ultra spectrometer (ESCALAB 250Xi, Thermo Fisher, Waltham, MA, USA), and the spectra were calibrated to the C 1s peak, which was located at 284.5 eV. Fourier transform infrared spectrometer (Thermo Scientific Nicolet iS20, Waltham, MA, USA) was utilized in order to acquire the surface information of GO and RGO. The resistance of the RGO films was measured using the four-point probe method (Keithley 2400, Beaverton, OR, USA).

Electrochemical Measurements
All of the electrochemical measurements were evaluated using a CHI-660E electrochemical workstation at room temperature. These measurements included cyclic voltammetry (CV), galvanostatic charging-discharging (GCD), and electrochemical impedance spectroscopy (EIS). The CV curves were recorded using different scan rates ranging from 5 mV s −1 to 10 V s −1 in the potential range of 0.0 to 1.0 V. Different current densities of 0.03 to 0.2 A cm −3 were used to conduct GCD tests. EIS was used to characterize frequencydependent behaviors in the frequency range of 100 kHz to 0.01 Hz at open circuit potential with a 5 mV ac perturbation.
For three-electrode system, the gravimetric capacitance (C g ) calculated by galvanostatic discharge curves was obtained by the following formula [30,31].
The specific volumetric capacitance values (C V , F cm −3 ) were evaluated from the CV curves by the formula [29].
The specific gravimetric capacitance values (C g , F g −1 ) were determined based on the CV curves by the formula [32].
The following formulas were used to define the specific energy density (E, mWh cm −3 ) and power density (P, W cm −3 ) under the dynamic conditions of cyclic voltammetry [29].
where I is the applied current (A), i(V) is the response current (A), V (cm 3 ) and m (g) are the total volume and the total mass of active electrode materials in the measuring system, ν is the potential scan rate (V s −1 ), ∆V (V) is the voltage range, i.e., 1 V during the charge-discharge process, ∆t (s) is the discharging time, respectively, C V (F cm −3 ) and C g (F g −1 ) are the volumetric and gravimetric capacitance of the device, E (mWh cm −3 ) and P (W cm −3 ) are the energy and power densities, respectively.

Results and Discussions
Following the HTHP method, we have recently successfully fabricated highly ordered graphite using natural lignite, which is regarded as an inferior coal resource that is difficult to graphitize compared with the commonly high carbon content anthracite [27]. The enhanced graphitization degree of lignite under HTHP can be characterized by XRD and Raman, as shown in Figure S1a,b. The XRD of raw lignite has a relatively broad (002) peak because there is a large amount of disordered structures and amorphous carbon in lignite. After the HTHP process, it is noted that the (002) peak of the sample becomes strong and sharp, located at 26.4 • (corresponding to a d-spacing of~3.36 Å), and a small (110) peak appears, implying the feature of lignite-based synthetic graphite with high orientation. The Raman spectra of samples after the HTHP process and raw lignite show the G-band and the D-band are around 1590 cm −1 and 1350 cm −1 , respectively [33]. Lignite-based synthetic graphite exhibits a prominent G-band and a weak D-band, and the ratio of D-band and G-band intensity is I D /I G = 0.17, further confirming the formation of high-quality graphite. The morphology and nanostructure analysis of lignite after HTHP pyrolysis were investigated by SEM and TEM. It can be seen that after the HTHP process the sample displays a sheet-like structure ( Figure S1c). Figure S1d provides a typical lowmagnification TEM image of lignite-based synthetic graphite. According to the pattern of the corresponding selected area electron diffraction (SAED), a typical diffraction ring with some bright spots of synthetic graphite is clarified, indicating a partial amorphization structure in the lignite-based synthetic graphite. Figure 1a shows a representative preparation procedure of RGO nanomaterials. Firstly, high-quality graphite was prepared by HTHP. Then the Hummers method was successfully utilized to obtain GO. Successively, the RGO was converted by chemical reduction using various molar concentrations of NaBH 4 and different treatment temperatures to improve its reducing capacity [34]. Systematic characterization analyses including XRD, Raman, FTIR, and XRS were used to investigate the effects of molar concentrations and reduction temperatures on the degree of reduction. The XRD patterns of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50 are shown in Figure 1b. The calculated interlayer distance from the (002) peak of GO is 8.27 Å (2θ = 10.7 • ), which is larger than the value for synthetic graphite (3.38 Å) [35]. The generation of oxygen-containing functional groups (hydroxyl, epoxy, and carboxyl) in the interlamellar gap of graphite is presumed to be the reason for the expanded interlayer distance. Under the reduction of NaBH 4 , the interlayer distance of RGO-RT-15 was further expanded to 8.9 Å (2θ = 9.8 • ). With further improving the molar concentration of the reducing agent, the d-spacing of RGO was slightly expanded from 9.6 Å (2θ = 9.1 • ) to 9.8 Å (2θ = 8.9 • ) for RGO-RT-50 and RGO-RT-100. When GO was reduced by 50 mM of reductant at 105 • C, the steep peak of RGO-105-50 vanished entirely, and a broad peak near 3.7 Å (2θ = 23.5 • ) occurred instead. This implies that high reduction temperature is more favorable for promoting deoxygenation of GO even under medium concentration of reducing agent.
Raman spectroscopy was adopted to explore the structural evolution under different reducing conditions. All of the spectra have significant D-band and G-band ( Figure 1c). The peak position of the G band for GO and RGO samples basically remains unchanged, while the peak of the D band is red-shifted as the degree of reduction increases. Moreover, the ratio of the intensities of the D and G bands (I D /I G ) in the Raman spectra can be used to monitor the formation of defects [36]. As depicted in Figure S2, the I D /I G ratio increases significantly even at low reductant concentrations in comparison with GO, whereas only a slight increase is observed in the samples of RGO-105-50, RGO-RT-100, and RGO-105-50. A higher I D /I G ratio can be assigned to more structural defects or partially disordered structure in the graphene sheets. This observation shows the removal of oxygen atoms would promote the conversion of sp 3 bonds to sp 2 bonds and introduce more defects during the reduction process [37]. The reduction evolution of GO by NaBH 4 was also recorded by FTIR spectroscopy analysis (Figure 1d). For the GO sample, the FTIR spectrum of GO shows a prominent peak at 3435 cm −1 , which resulted from the presence of the O-H stretching of intercalated water. The peaks at about 1735 cm −1 and 1628 cm −1 are assigned to carboxyl/carbonyl (C=O) and aromatic (C=C) stretching groups, respectively [38]. The peak that occurs at 1380 cm −1 is associated with the -CH 3 group, while the peak that occurs at 1104 cm −1 stems from -C-O-stretching. The deformation vibration of the -CH outside the aromatic ring is responsible for the low-intensity absorption bands that can be seen between 900 and 700 cm −1 . When GO is reduced by NaBH 4 at RT, it is clearly seen that the C=O and C-O species do not disappear, demonstrating that the reduction process is not completed. When the reduced temperature up to 105 • C, there are fewer peaks observed in RGO-105-50, suggesting an efficient removal of oxygen-containing functional groups.
peak that occurs at 1380 cm is associated with the −CH3 group, while the peak that occurs at 1104 cm −1 stems from −C−O− stretching. The deformation vibration of the −CH outside the aromatic ring is responsible for the low-intensity absorption bands that can be seen between 900 and 700 cm −1 . When GO is reduced by NaBH4 at RT, it is clearly seen that the C=O and C−O species do not disappear, demonstrating that the reduction process is not completed. When the reduced temperature up to 105 °C, there are fewer peaks observed in RGO-105-50, suggesting an efficient removal of oxygen-containing functional groups. The chemical states of different elements and the presence of functional groups of GO and RGO samples were probed by XPS. Carbon and oxygen are the most prominent elements in the survey scan spectra of all samples (Figure 1e), while there is a Na 1s peak in sample RGO-105-50, which is derived from the residual NaBH4. The high resolution C 1s region spectra of RGO-105-50 demonstrates four types of carbon bonds at binding energies of 284.5, 286.1, 286.7, and 288.4 eV, corresponding to C−C, C−O, C=O, and O−C=O, respectively [39]. There are the same types of carbon bonds in GO, RGO-RT-15, RGO-RT-50, and RGO-RT-100 ( Figure S3a-d). It can be found that the peak intensity of C−O bonds displays a remarkable decrease contrasted with that of C−C as the molar concentrations of NaBH4 increase. When the reduced temperature is 105 °C, C−O group is effectively removed even at the mid-range molar concentration of NaBH4. Table S1 lists the areas of the contributing peaks. During chemical reduction, it is of particular note that NaBH4 can effectively reduce C−O, while having a moderate capability in reduction of C=O, and is The chemical states of different elements and the presence of functional groups of GO and RGO samples were probed by XPS. Carbon and oxygen are the most prominent elements in the survey scan spectra of all samples (Figure 1e), while there is a Na 1s peak in sample RGO-105-50, which is derived from the residual NaBH 4 . The high resolution C 1s region spectra of RGO-105-50 demonstrates four types of carbon bonds at binding energies of 284.5, 286.1, 286.7, and 288.4 eV, corresponding to C-C, C-O, C=O, and O-C=O, respectively [39]. There are the same types of carbon bonds in GO, RGO-RT-15, RGO-RT-50, and RGO-RT-100 ( Figure S3a-d). It can be found that the peak intensity of C-O bonds displays a remarkable decrease contrasted with that of C-C as the molar concentrations of NaBH 4 increase. When the reduced temperature is 105 • C, C-O group is effectively removed even at the mid-range molar concentration of NaBH 4 . Table S1 lists the areas of the contributing peaks. During chemical reduction, it is of particular note that NaBH 4 can effectively reduce C-O, while having a moderate capability in reduction of C=O, and is basically invalid for the reduction of O-C=O groups. Figure 1g and Figure S3e When GO was reduced at room temperature, the C:O ratio decreased slightly as the NaBH 4 concentration increased. However, C:O atomic ratio of RGO-105-50 increased dramatically to 4.04, confirming reduction temperature is more conducive to the reduction effect [40]. The morphology and structure of the GO and RGO nanosheets were further investigated by TEM. Figure 2 depicts the TEM images of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50, all of these displaying half-transparent sheet-like structure. When lignite was pretreated by HTHP, the size of lignite-derived GO was found to be remarkably increased to 2.8 µm larger than that of GO (size: 40-200 nm) obtained by leaching low-grade coal with HNO 3 [26]. As with the enhanced reduction degree, the RGO show the less obvious folding features, affirming the removement of oxygen functional groups in agreement with the FTIR and XPS results. Using the HRTEM image (Figure 2d), the inter-planar distance is calculated to be between 0.34 and 0.36 nm, which is close to the value of 0.337 nm for graphite [35]. served in C 1s spectra. The C:O atomic ratios of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50, obtained by XPS analysis, are 1.95, 1.65, 1.65, 1.59, and 4.04. When GO was reduced at room temperature, the C:O ratio decreased slightly as the NaBH4 concentration increased. However, C:O atomic ratio of RGO-105-50 increased dramatically to 4.04, confirming reduction temperature is more conducive to the reduction effect [40].
The morphology and structure of the GO and RGO nanosheets were further investigated by TEM. Figure 2 depicts the TEM images of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50, all of these displaying half-transparent sheet-like structure. When lignite was pretreated by HTHP, the size of lignite-derived GO was found to be remarkably increased to 2.8 μm larger than that of GO (size: 40-200 nm) obtained by leaching low-grade coal with HNO3 [26]. As with the enhanced reduction degree, the RGO show the less obvious folding features, affirming the removement of oxygen functional groups in agreement with the FTIR and XPS results. Using the HRTEM image (Figure 2d), the inter-planar distance is calculated to be between 0.34 and 0.36 nm, which is close to the value of 0.337 nm for graphite [35].  The electrochemical energy storage performance of the RGO materials were firstly evaluated by CV and GCD measurements using the 3-electrode system (Figures S5a,b, S6a,b, S7a,b, and S8a,b). As we can see, the CV curves of the RGO materials possess the characteristic of quasirectangular nature. In addition, the GCD plots exhibit quasi-linear nature which is in the agreement with the CV results. Both the CV and GCD curves indicate the RGO-105-50 electrode has the highest specific capacitance. Specifically, the gravimetric capacitance of the RGO-105-50 is 30.1 F g −1 at 0.05 A g −1 , superior to RGO-RT-15 (1.1 F g −1 ), RGO-RT-50 (1.5 F g −1 ), and RGO-RT-100 (5.9 F g −1 ), as shown in Figures S5c, S6c, S7c, and S8c. The Nyquist plot of RGO-105-50 exhibits an even steeper slope of the linear plots in the low-frequency regions, suggesting the high reduction degree is beneficial to facilitate ion diffusion in energy storage process (Figures S5d, S6d, S7d, and S8d).  The electrochemical energy storage performance of the RGO materials were firs evaluated by CV and GCD measurements using the 3-electrode system ( Figures S5a  S6a,b, S7a,b, and S8a,b). As we can see, the CV curves of the RGO materials possess t characteristic of quasi-rectangular nature. In addition, the GCD plots exhibit quasi-lin nature which is in the agreement with the CV results. Both the CV and GCD curves in cate the RGO-105-50 electrode has the highest specific capacitance. Specifically, the gra imetric capacitance of the RGO-105-50 is 30.1 F g −1 at 0.05 A g −1 , superior to RGO-RT (1.1 F g −1 ), RGO-RT-50 (1.5 F g −1 ), and RGO-RT-100 (5.9 F g −1 ), as shown in Figures S5c, S  S7c, and S8c. The Nyquist plot of RGO-105-50 exhibits an even steeper slope of the line plots in the low-frequency regions, suggesting the high reduction degree is beneficial facilitate ion diffusion in energy storage process ( Figures S5d, S6d, S7d, and S8d).
The typical electrochemical properties of the lignite-derived GO and RGO materi were also measured in two-electrode cell configurations using PVA/H3PO4 electrolyte schematic diagram of flexible supercapacitor is depicted in Figure 4a, containing electro film, H3PO4/PVA gel electrolyte, and Au-coated PET substrate. The cyclic voltamme method was used to evaluate the electrochemical performance of the lignite-derived G and RGO supercapacitors at different scan rates ranging from 2 mV s −1 to 10 V s −1 (Figu S9a-d, S10a-d, S11a-d, S12a-d, and S13a-d). It can be seen that, with the exception of t CV curves of GO at low scan rates of less than 0.01 V s −1 , other CV curves exhibit nea ideal rectangular profile. The deviation from a rectangular shape is most likely caused the irreversible reaction of oxygen functional groups on graphene with residual water the electrolyte. As shown in Figures S9e, S10e, S11e, S12e, and S13e, the energy stora performance of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50 were furth measured by GCD. The curves of RGO-105-50 show more like symmetric linear sha than other lignite-derived RGO supercapacitors, indicating a good capacitive perf mance. Meanwhile, the electrochemical performance of the c-RGO-105-50 electrode w also conducted under the same test condition (Figure S14a-e). Figure 4b shows the C curves of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50 acquired at a sc rate of 0.05 V s −1 . Obviously, the current density of RGO-105-50 is higher than those of G and other RGO reduced with NaBH4 at RT, implying improved electrochemical perf mance due to efficient reduction degree. Besides, the discharge time of RGO-105-50 longer than those of other four kinds of electrodes (Figure 4c). The typical electrochemical properties of the lignite-derived GO and RGO materials were also measured in two-electrode cell configurations using PVA/H 3 PO 4 electrolyte. A schematic diagram of flexible supercapacitor is depicted in Figure 4a, containing electrode film, H 3 PO 4 /PVA gel electrolyte, and Au-coated PET substrate. The cyclic voltammetry method was used to evaluate the electrochemical performance of the lignite-derived GO and RGO supercapacitors at different scan rates ranging from 2 mV s −1 to 10 V s −1 (Figures S9a-d, S10a-d, S11a-d, S12a-d, and S13a-d). It can be seen that, with the exception of the CV curves of GO at low scan rates of less than 0.01 V s −1 , other CV curves exhibit nearly ideal rectangular profile. The deviation from a rectangular shape is most likely caused by the irreversible reaction of oxygen functional groups on graphene with residual water in the electrolyte. As shown in Figures S9e, S10e, S11e, S12e, and S13e, the energy storage performance of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50 were further measured by GCD. The curves of RGO-105-50 show more like symmetric linear shape than other lignite-derived RGO supercapacitors, indicating a good capacitive performance. Meanwhile, the electrochemical performance of the c-RGO-105-50 electrode was also conducted under the same test condition (Figure S14a-e). Figure 4b shows the CV curves of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50 acquired at a scan rate of 0.05 V s −1 . Obviously, the current density of RGO-105-50 is higher than those of GO and other RGO reduced with NaBH 4 at RT, implying improved electrochemical performance due to efficient reduction degree. Besides, the discharge time of RGO-105-50 is longer than those of other four kinds of electrodes (Figure 4c).
The specific capacitance values of the lignite-derived GO and RGO electrodes obtained from the CV curves at different scan rates are presented in Figure 4d. The RGO-105-50 electrode has the highest specific capacitance (30.6 F cm −3 ) at 2 mV s −1 , which is superior to GO (1.5 F cm −3 ), RGO-RT-15 (6.7 F cm −3 ), RGO-RT-50 (7.3 F cm −3 ), and RGO-RT-100 (16.5 F cm −3 ). Besides, the RGO-105-50 electrode maintains maximum capacitance in the entire measurement region, while the difference in specific capacitance of these supercapacitors becomes inconspicuous as the scan rate increasing to 10 V s −1 . Excitingly, the storage capacitance of RGO-105-50 is slightly better than that of c-RGO-105-50 (28.02 F cm −3 at 2 mV s −1 as shown in Figure S14f). Furthermore, the gravimetric capacitance of the RGO-105-50 is 51.5 F g −1 better than that of the GO and other RGO supercapacitors. The electrical resistivity of GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50 films were measured by a four-point probe method. The electrical resistivity of RGO-105-50 possesses the minimum value, which is one of the reasons that RGO-105-50 has the best energy storage performance. To further investigate the origin of the excellent energy storage performance, EIS is used to evaluate the electron/ion transport kinetics of electrodes [41]. Figure 5a depicts the obtained Nyquist plots of the GO and RGO materials at frequencies ranging from 10 5 to 10 −2 Hz. In the low-frequency zone, the straight-line part of the RGO-105-50 is more parallel to the imaginary axis than that of GO and other RGO materials, indicating that this electrode has a higher mobility of electron and ion. In addition, the lack of a semicircle at the high-frequency zone of the RGO-105-50 electrode suggests a low charge transfer resistance, which is beneficial for the rapid formation of an electrical double layer in the electrode. The equivalent series resistance (ESR), which includes the resistance of the electrolyte, intrinsic resistance of substrate, and contact resistance, can be obtained by the intersection of the curve with the real axis at high frequency. According to the inset of Figure 5a, the ESR values for GO, RGO-RT-15, RGO-RT-50, RGO-RT-100, and RGO-105-50 are 8.7 Ω, 5.2 Ω, 4.7 Ω, 4.4 Ω, and 2.7 Ω, respectively, indicating that the RGO-105-50 material has high electrode conductivity in agreement with the result of electrical resistivity tests (Figure 4f). It can be demonstrated that the RGO-105-50 electrode shows the advantage of ion diffusion during the electrochemical energy storage process, benefiting from the efficient reduction of C=O and C-O groups induced by NaBH 4 .
To further display the practical applications of the as-obtained lignite-derived RGO flexible supercapacitors, the Ragone plot of RGO-105-50 compared with previously reported energy storage devices, such as a lithium thin-film battery (4 V/500 mA h) and an aluminum electrolytic capacitor (3 V/300 mF) [42], is shown in Figure 5b. Note that the RGO-105-50 device could achieve a high energy density of 4.24 mW h cm −3 and a power density of 2.89 W cm −3 , surpassing commonly studied energy storage devices based on two-dimensional and carbon-based materials, such as reduced graphene film by CH 4 plasma reduction (1.6 W cm −3 and 0.15 mWh cm −3 ) [43], laser-scribed graphene (3.8 W cm −3 and 0.11 mWh cm −3 ) [44], nanocomposites MoS 2 @Ni(OH) 2 (11 W cm −3 and 5.2 mWh cm −3 ) [45], three-dimensional graphene and metal-organic framework composite (0.02 W cm −3 and 0.03 mWh cm −3 ) [46], CNT film (0.01 W cm −3 and 0.08 mWh cm −3 ) [47], and black phosphorus nanosheets (8.83 W cm −3 and 2.47 mWh cm −3 ) [29]. The comparison of energy storage performance of RGO-105-50 with other electrode materials is shown in Table S2. To further display the practical applications of the as-obtained lignite-derived RGO flexible supercapacitors, the Ragone plot of RGO-105-50 compared with previously reported energy storage devices, such as a lithium thin-film battery (4 V/500 mA h) and an aluminum electrolytic capacitor (3 V/300 mF) [42], is shown in Figure 5b. Note that the RGO-105-50 device could achieve a high energy density of 4.24 mW h cm −3 and a power density of 2.89 W cm −3 , surpassing commonly studied energy storage devices based on two-dimensional and carbon-based materials, such as reduced graphene film by CH4 plasma reduction (1. One of the most critical properties to be evaluated before putting flexible supercapacitors into practical application is their long-term cycling stability. In this work, the RGO-105-50 device was tested for stability at 0.5 V s −1 . The CV curve of RGO-105-50 after 20,000 long-term cycles (Figure 5c) possesses a larger integral area than its initial CV curve, confirming its remarkable cycling stability. As depicted in Figure 5d, the retention rate of RGO-105-50 shows slight growth over 20,000 cycles, finally reaching 112.3% of the initial capacitance, superior to that of c-RGO-105-50 fabricated by the same processing technic (101.8% after 20,000 long-term cycles as shown in Figure S14g). The enhancement of the capacity after long cycling is due to the adequate penetration of the electrolyte and the gradual activation of the surface. Furthermore, the CV curves of RGO-105-50 were measured under various bending angles in the range of 0-180 • (Figure 5e). Notably, RGO-105-50 flexible supercapacitor subjected to different bending angles show similar CV curves even at extreme degrees of bending (Figure 5f), and the capacitance remains 79.5% of the capacitance of flat state without destroying the structural integrity, illustrating the excellent mechanical flexibility and potential application in flexible energy storage devices.
To further explore the micro-structure and surface chemical compositions of RGO-105-50 after 20,000 cycles, SEM and XPS were adopted. After electrochemical measurements, RGO nanoflakes tend to stack on each other, displaying a smooth surface ( Figure S15a). From the XPS spectra of C 1s and O 1s in Figure S15b-c, the C-O bond is significantly enhanced in complex electrochemical energy storage processes. The change in relative peak strength of C-O bond is believed to offer a good pseudocapacitance effect, which is responsible for the improvement of the retention rate of RGO-105-5 [48]. In addition, Figure S16 depicts the obtained Nyquist plots of RGO-105-50 after 20,000 cycles, the EIS test of RGO-105-50 after 20,000 cycles shows that the ESR increases to 8.8 Ω and the low frequency straight line is closer to the real axis, indicating the slower ion diffusion process and lower ion transport.

Conclusions
In summary, we have demonstrated an optimizing route for the mass production of large-size GO derived from low-value lignite, that benefits from the highly ordered graphite material prepared by the HTHP method. Lignite-derived RGO can be obtained by the efficient reduction of NaBH 4 , in which the reduction temperature is more conducive to the degree of reduction compared with molar concentrations of reductive agent. The RGO-105-50 electrode prepared by vacuum filtration features porosity and high conductivity, which improves the contact area between electrode and electrolyte and facilitates electron/ion transport. As is expected, the RGO-105-50 electrode exhibits a high volumetric capacitance of 30.6 F cm −3 , an excellent volumetric energy density of 4.2 mW h cm −3 , a power density of 5.9 W cm −3 , and exceptional flexibility, i.e., 79.5% of the initial specific capacitance even under 180 • bending. Remarkably, the specific capacitance of the assembled RGO-105-50 flexible supercapacitors remains about 112.3% after 20,000 cycles, which is higher than that of RGO prepared by the chemical reduction of purchased GO (101.8% retention of the initial capacitance after 20,000 cycles). The current work not only demonstrates an effective route for the synthesis of the large-size RGO sheets in the manner of large quantity and low cost, but also shows that the lignite-derived RGO can be supposed as a promising candidate of high-performance electrode materials for flexible supercapacitors.

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