The Influence of Reduced Graphene Oxide on the Texture and Chemistry of N,S-Doped Porous Carbon. Implications for Electrocatalytic and Energy Storage Applications

In this work, we study the influence of reduced graphene oxide (rGO) on the morphology and chemistry of highly porous N,S-doped carbon cryogels. Simultaneously, we propose an easily upscalable route to prepare such carbons by adding graphene oxide (GO) in as-received suspended form to the aqueous solution of the ι-carrageenan and urea precursors. First, 1.25–5 wt% GO was incorporated into the dual-doped polymer matrix. The CO2, CO, and H2O emitted during the thermal treatments resulted in the multifaceted modification of the textural and chemical properties of the porous carbon. This facilitated the formation of micropores through self-activation and resulted in a substantial increase in the apparent surface area (up to 1780 m2/g) and pore volume (up to 1.72 cm3/g). However, adding 5 wt% GO led to overactivation. The incorporated rGO has an ordering effect on the carbon matrix. The evolving oxidative species influence the surface chemistry in a complex way, but sufficient N and S atoms (ca. 4 and >1 at%, respectively) were preserved in addition to the large number of developing defects. Despite the complexity of the textural and chemical changes, rGO increased the electrical conductivity monotonically. In alkaline oxygen reduction reaction (ORR) tests, the sample with 1.25 wt% GO exhibited a 4e− mechanism and reasonable stability, but a higher rGO content gradually compromised the performance of the electrodes. The sample containing 5 wt% GO was the most sensitive under oxidative conditions, but after stabilization it exhibited the highest gravimetric capacitance. In Li-ion battery tests, the coulombic efficiency of all the samples was consistently above 98%, indicating the high potential of these carbons for efficient Li-ion insertion and reinsertion during the charge–discharge process, thereby providing a promising alternative for graphite-based anodes. The cell from the 1.25 wt% GO sample showed an initial discharge capacity of 313 mAh/g, 95.1% capacity retention, and 99.3% coulombic efficiency after 50 charge–discharge cycles.


Introduction
The fascinating microstructural and surface properties of biomass-based carbon materials justify the extensive interest in their wide range of application in various fuel cells [1][2][3], supercapacitors [4,5], and alkaline, zinc, or ammonium ion batteries [6][7][8][9][10][11][12][13][14]. Compared ter the access of the electrolyte delivering dissolved oxygen into the pore system [48]. The thiophenic sulfur compounds appear to be particularly active in enhancing the physisorption of the oxygen from electrolytes in small carbon pores [48]. The S atom incorporated in aromatic rings causes a slightly positive charge on the neighboring C atoms [48,49] and brings additional hydrophobicity to the surface, thus promoting the adsorption of molecular oxygen [42].
DFT calculations revealed that dual S and N doping resulted in a large number of active sites through the redistribution of spin and charge densities, thus synergistically improving the ORR performance [49]. The atomic level distribution of (hetero)atoms in the porous carbon matrix significantly improved the rate of the ORR, which was otherwise too slow for practical applications [50,51]. They also affected the performance of supercapacitors [52]. They reduce or may even prevent the need for Pt loading, which is an economic obstacle in commercializing fuel cells. N, P, and S atoms in biomass-based carbon electrodes act as extra active Li-ion storage sites [53,54]. The kinetics of the electrochemical reaction in heteroatom doped carbons are much faster, resulting in a higher capacity and better rate performance. The C-S defect sites induce torsion in the graphitic layer structure, and the extended interlayer spacing facilitates the insertion of lithium ions during the charge-discharge cycling of the cells [55].
Nanocarbons (carbon nanotubes, carbide, and carbonitride graphene) themselves may form an electrically conductive interconnected aerogel matrix through van der Waals interactions [56,57]. These particles can be also used as templating agents [58]. Graphene, due to its outstanding properties (high electrical and thermal conductivity and good mechanical strength) plays a particularly important role in the development of a new generation of batteries. Their surface area, mechanical flexibility, and broad electrochemical window [59][60][61] make graphene-based electrodes promising for fuel cells [62], batteries [63][64][65], supercapacitors [66][67][68][69], etc. However, in polar media, graphene tends to aggregate [70]. Its hydrophilic derivative, graphene oxide (GO), owing to its oxygen-containing functional groups, is easier to disperse in the precursors. The thermal/annealing treatment during the synthesis converts the GO to reduced GO (rGO), and thus, the desired structure and physicochemical properties of the electrocatalytically more attracting graphene are at least partly restored [71]: Incorporation of 0.23-0.46 wt% well-dispersed GO into the resorcinol-formaldehyde (RF) polymer xerogel precursor resulted in a carbon xerogel of high porosity and of excellent electrical conductivity: when used as an electrode in aqueous supercapacitors at high current density, the capacitance and the power were enhanced by 25 and 100%, respectively, compared to the undoped carbon xerogel [72]. It was also concluded that the electrical conductivity seemed to play a more important role than the specific surface area in the electrode performance of mesoporous carbon xerogels. GO loadings exceeding 2 wt% resulted in lower energy and power densities in spite of the increased electric conductivity [73]. In Li-ion batteries, graphene-doped carbon/carbon aerogel electrodes can accommodate lithium more easily than the common graphite anode [59], as in addition to the intercalation mechanism, they exhibit fast lithium adsorption [31,[74][75][76], defect trapping [75], build-up of faradaic capacitance [77,78], etc.
In spite of the numerous studies on the performance of electrodes fabricated from carbons doped either with two heteroatoms or with carbon nanoparticles, their simultaneous application is rare [79,80], and none of the studies addressed the influence of the nanoparticle concentration.
Earlier we observed that while incorporation of N or reduced GO alone in a resorcinolformaldehyde-based carbon aerogel did not alter the 2e − the mechanism of the ORR, when both of them were present, the dominant pathway changed from the slow 2e − to the more efficient 4e − route [81]. Recently, we found that a carbon electrode fabricated from an ι-carrageenan-urea biopolymer hydrogel precursor displayed almost the same oxygen reduction reaction (ORR) performance as a commercial 20 wt % Pt/C electrode. The electrode was found to be stable and exhibited a 4-electron behavior in alkaline ORR [82].
Based on these findings, in this work, we focus on the influence of rGO on the morphology and chemistry of highly porous carbon materials obtained from the same biopolymer cryogel. Simultaneously, we propose an easily upscalable route to prepare rGO doped N,S containing biopolymer based carbons. In reference works, the incorporated carbon nanoparticle content was studied in a wide range [72,73,79,80]. Previously, we found that in responsive 3D poly(N-isopropylacrylamide) networks the percolation threshold was around 5 wt% [83]. Therefore, in this work the GO content was varied between 1.25 and 5 wt%. The texture of the carbon cryogels was characterized by electron microscopy imaging, N 2 and CO 2 adsorption measurements, Raman spectroscopy, and powder X-ray diffraction (XRD). X-ray photoelectron spectroscopy (XPS), SEM-supported electron dispersive spectroscopy (EDS), ultimate elemental analysis, and FTIR were used to study the surface and bulk chemistry. The samples were probed as electrodes (i) in a fuel cell under condition, similar to alkaline anion exchange membrane fuel cells (AAEMFCs) and (ii) in Li-ion batteries.

Synthesis
rGO containing N,S double-doped carbon aerogels were synthesized, using ι-carrageen an (commercial grade, Sigma Aldrich, Budapest, Hungary), urea pearls (98%, Sigma Aldrich, Budapest, Hungary), and aqueous GO suspension with a method adapted from Li et al. [17]. A detailed scheme of the synthesis procedure is presented (Figure 1).
Earlier we observed that while incorporation of N or reduced GO alone in a resorcinol-formaldehyde-based carbon aerogel did not alter the 2e − the mechanism of the ORR, when both of them were present, the dominant pathway changed from the slow 2e − to the more efficient 4e − route [81]. Recently, we found that a carbon electrode fabricated from an ι-carrageenan-urea biopolymer hydrogel precursor displayed almost the same oxygen reduction reaction (ORR) performance as a commercial 20 wt % Pt/C electrode. The electrode was found to be stable and exhibited a 4-electron behavior in alkaline ORR [82].
Based on these findings, in this work, we focus on the influence of rGO on the morphology and chemistry of highly porous carbon materials obtained from the same biopolymer cryogel. Simultaneously, we propose an easily upscalable route to prepare rGO doped N,S containing biopolymer based carbons. In reference works, the incorporated carbon nanoparticle content was studied in a wide range [72,73,79,80]. Previously, we found that in responsive 3D poly(N-isopropylacrylamide) networks the percolation threshold was around 5 wt% [83]. Therefore, in this work the GO content was varied between 1.25 and 5 wt%. The texture of the carbon cryogels was characterized by electron microscopy imaging, N2 and CO2 adsorption measurements, Raman spectroscopy, and powder X-ray diffraction (XRD). X-ray photoelectron spectroscopy (XPS), SEM-supported electron dispersive spectroscopy (EDS), ultimate elemental analysis, and FTIR were used to study the surface and bulk chemistry. The samples were probed as electrodes (i) in a fuel cell under condition, similar to alkaline anion exchange membrane fuel cells (AAEMFCs) and (ii) in Li-ion batteries.

Synthesis
rGO containing N,S double-doped carbon aerogels were synthesized, using ι-carrageenan (commercial grade, Sigma Aldrich, Budapest, Hungary), urea pearls (98%, Sigma Aldrich, Budapest, Hungary), and aqueous GO suspension with a method adapted from Li et al. [17]. A detailed scheme of the synthesis procedure is presented (Figure 1). Figure 1. Scheme of the synthesis procedure. The GO is added in as-received suspended form, and the biopolymer is formed in aqueous medium, making the synthesis route green and easily upscalable.

Characterization Methods
Scanning electron micrographs were taken by a Zeiss Sigma 300 FESEM field emission scanning electron microscope (Carl Zeiss QEC GmbH, Oberkochen, Germany). Lowtemperature (−196.15 • C) nitrogen adsorption measurements were performed after 24 h degassing at 110 • C on a NOVA 2000e (Quantachrome, Boynton Beach, FL, USA) automatic volumetric instrument. The apparent surface area S BET was determined using the Brunauer-Emmett-Teller (BET) model [86]. The pore volume V 0.98 was estimated from the amount of vapor adsorbed at p/p 0 = 0.98, assuming that the adsorbed gas fills the pores as liquid. The Dubinin-Radushkevich (DR) plot [87] was used to calculate the micropore volume V micro . The pore size distribution was computed by the Quenched Solid Density Functional Theory (QSDFT) for slit/cylindrical pore geometry [88]. Carbon dioxide adsorption was measured at 0 • C up to atmospheric pressure with an AUTOSORB-1 (Quantachrome, Boynton Beach, FL, USA) computer-controlled analyzer. Evaluation of the primary adsorption data was performed with the Quantachrome AsiQwin software (version 3.0). Raman spectra were obtained using a LabRAM (Horiba Jobin Yvon, Kyoto, Japan) instrument. The laser source was a λ = 532 nm Nd-YAG (laser power at the focal point 15 mW). A 0.6 OD filter was used to reduce the power of the beam. Parameter optimization and data analysis were performed by LabSpec 5 software. Powder X-ray diffractograms (XRD) were obtained in the range 2θ = 10 • -130 • with an X'Pert Pro MPD (PANalytical Bv., Almelo, The Netherlands) X-ray diffractometer using an X'celerator type detector and monochromatic Cu Kα radiation with a Ni filter foil (λ = 1.5406 Å).
The ultimate elemental analysis was performed in triplicate on an Elementar Vario Micro (CHNS) instrument. The lack of metal traces was confirmed with SEM/EDS (JEOL JSM 6380LA, Tokyo, Japan). X-ray photoelectron spectra were recorded on a Kratos XSAM 800 spectrometer operating in fixed analyzer transmission mode, using Mg K α1,2 (1253.6 eV) excitation. The pressure of the analysis chamber was lower than 1 × 10 −7 Pa. Survey spectra were recorded in the 150-1300 eV range in 0.5 eV steps. The photoelectron lines of C1s, O1s, N1s, and S2p were measured in 0.1 eV steps with 1 s dwell time. The spectra were referenced to the energy of the C1s line of the sp 2 type graphitic carbon, set at 284.3 ± 0.1 eV binding energy (BE). Peak decomposition was performed after Shirleytype background removal using a Gaussian-Lorentzian peak shape with a 70:30 ratio. Details of the applied fitting procedure are described elsewhere [89]. Quantitative analysis, based on integrated peak intensity, was performed by the XPS MultiQuant program [90], applying the conventional infinitely thick layer model using the experimentally determined photoionization cross-section data of Evans et al. [91] and the asymmetry parameters of Reilman et al. [92]. Attenuated total reflectance Fourier transform infrared spectroscopy (FTIR-ATR) were recorded on powdered carbons in the 4000-400 cm −1 wavenumber range at a resolution of 4 cm −1 by 32 scans using a Tensor 27 (Bruker Optik GmbH, Leipzig, Germany) spectrophotometer equipped with a Platinum ATR unit A225. The crystal was made of diamond having a refractive index of 2.4. For the background signal, the measured medium was air. Since absorption by the powders was very strong, a moderate polynomial baseline correction and smoothing were required.

Electrochemical Performance Tests
The electrical conductivity was estimated by a laboratory made instrument. The powdered samples were gradually compressed (0.5-5 MPa) in a rigid polytetrafluoroethylene (PTFE) tube (1 cm 2 internal cross section) with two copper bolts (one on each side of the Nanomaterials 2023, 13, 2364 6 of 22 tube). The conductivity was estimated from the pressure dependence of the resistance [93]. The electrocatalytic ORR tests were performed using a glassy carbon (GC) rotating disc electrode (RDE, Pine Research Instrumentation, Durham, NC, USA). The ink for the working electrodes was prepared by dispersing 2 mg powdered annealed carbon (CA, CAGO50, CAGO100, or CAGO200) in a mixture of 1.6 mL MilliQ water, 0.4 mL isopropyl alcohol, and 8 µL 5% Nafion ® solution. After 30 min sonication, the ink was pipetted onto the dry-mirror-polished GC and dried at room temperature. The loading varied between 50 and 400 µg/cm 2 . The measurements were implemented in 0.1 M KOH electrolyte using three-electrode systems with a hydrogen electrode as a reference and Pt wire as a counter electrode in a three-compartment PTFE cell [82]. All potentials are given vs. the reversible hydrogen electrode (RHE).
The working electrodes for the Li-ion storage studies were prepared by mixing the carbon samples and polyvinylidene fluoride binder (PVDF 99.9%, Solvay, Budapest, Hungary) in a weight ratio of 95:5 in a ball mill for 1 h. After adding an adequate amount of N-methyl-2-pyrrolidone (NMP 99.9%, Sigma-Aldrich, Budapest, Hungary), the slurry was cast on to a copper foil using an automatic film applicator (BYK Gardner GmbH, Geretsried, Germany). The coated film was dried overnight at 70 • C in a vacuum oven. Then, 14 mm diameter circular discs were cut from the dried film using a manual disc cutter (Berg &Schmidt GmbH, Hamburg, Germany). CR2032 coin-type cells were assembled in an argon-filled glove box using the as-prepared coated disc and Li metal disc (99.9%, Sigma-Aldrich, Budapest, Hungary) as working and counter electrodes, respectively. Whatman glass fiber served as the separator, and 1.0 M lithium hexafluorophosphate (LiPF 6, 99.95%, Sigma-Aldrich, Budapest, Hungary) in ethylene carbonate (EC, 98%, Sigma-Aldrich, Budapest, Hungary) and diethyl carbonate (DEC, 99%, Sigma-Aldrich, Budapest, Hungary) in 1:1 v% served as the electrolyte. The cells were assembled using a manual crimping machine (MTI MSK 110, Xiang'an, China). The electrochemical performance of the half-cell was performed using an electrochemical workstation (Biologic VMP 300, Seyssinet-Pariset, France) at room temperature. The galvanostatic charge-discharge cycling of the half-cell was conducted in a potential window of 0-3 V (vs. Li + /Li).

Development of the Morphology during the Annealing
The FESEM images demonstrate the complexity of the pore structure of the annealed carbon gels ( Figure 2). The early inclusion of GO already in the precursor suspension influences the morphology developing during gelation and cryogenic drying and is reflected also in the final annealed state (Figure 2c). It is apparent that the reduced GO formed during the heat treatments was evenly distributed in the carbon matrix. The typical nanopores within the CA matrix were 2-3 nm ( Figure S1). The yield of the synthesis steps is reported in Table S1. The low temperature (−196 °C) N2 adsorption isotherms of the annealed carbon samples having various GO content are shown in Figure 3a. All the isotherms are of composite Type IV + Type II with an H4 hysteresis loop displaying a sharp step at p/p0 = 0.45 [94]. The shape and type of the hysteresis loops suggest an interconnected pore network of micro-, meso-, and macropores. Numerical data deduced from the measured isotherms are presented in Table 1. Since the macropores are not totally filled with condensed nitrogen, the liquid equivalent volume V0.98 was determined at p/p0 = 0.98. The low temperature (−196 • C) N 2 adsorption isotherms of the annealed carbon samples having various GO content are shown in Figure 3a. All the isotherms are of composite Type IV + Type II with an H4 hysteresis loop displaying a sharp step at p/p 0 = 0.45 [94]. The shape and type of the hysteresis loops suggest an interconnected pore network of micro-, meso-, and macropores. Numerical data deduced from the measured isotherms are presented in Table 1. Since the macropores are not totally filled with condensed nitrogen, the liquid equivalent volume V 0.98 was determined at p/p 0 = 0.98. The low temperature (−196 °C) N2 adsorption isotherms of the annealed carbon samples having various GO content are shown in Figure 3a. All the isotherms are of composite Type IV + Type II with an H4 hysteresis loop displaying a sharp step at p/p0 = 0.45 [94]. The shape and type of the hysteresis loops suggest an interconnected pore network of micro-, meso-, and macropores. Numerical data deduced from the measured isotherms are presented in Table 1. Since the macropores are not totally filled with condensed nitrogen, the liquid equivalent volume V0.98 was determined at p/p0 = 0.98.    The most spectacular effect of the GO was in the micropore region of the isotherm, which ran almost parallel in the mesoporous region. Incorporation of GO increased not only the carbon content of the polymer matrix but also the volatile and gaseous products that were released during the thermal treatments (see also Equation (1)). The CO, CO 2 , and water act as self-activating agents during their path through the material and alter both its texture and chemistry. The 1.25-2.5 wt% GO samples (CAGO50 and CAGO100) increased the pore volumes and the apparent surface area almost proportionally, while addition of 5 wt% GO (CAGO200) resulted in a drop in both the surface area and pore volume. In this case, the amount of evolving oxidative reactants was too high and overactivated the sample [70]. Additionally, the ultra-micropore region was revealed by the CO 2 adsorption measurements performed at 0 • C. The combined pore size distribution curves in Figure 3b indicate the presence of ultra-micropores as well as significant micro-and narrow mesoporosity in the range of 1.7-6 nm. The increasing amount of oxidative volatile gases evolving during the thermal decomposition of the GO may also be responsible for the decreasing trend in ultra-microporosity. The Raman spectra in Figure S2 exhibit the iconic D (~1350 cm −1 ; defects, edges, and disordered carbon sites) and G (~1580 cm −1 ; vibration of sp 2 -hybridized graphitic carbon) band regions typical of carbon materials [95]. The first and second order regions of the spectra were deconvoluted for further analysis ( Figure S2b, Tables S2 and S3 [96][97][98][99]). The D band (around 1350 cm −1 ) indicates the aromatic clusters with edges/borders of mainly amorphous carbon (sp 3 -bonded) structures [100]. The blue shift of the G band may be caused by internal stress developing when the incorporated GO and/or newly developing clusters collide during the annealing treatment at 1000 • C [100,101]. The D band, associated with disordered graphenic lattices (around 1620 cm −1 ) appeared only in the GO-free CA sample. On the other hand, the second-order region was detected only in the rGO-containing samples. The trend in the I D /I G ratio along with the relative sharpening of the G band and the reduced intensity D2033 in the GO-added samples confirm an ordering effect of the forming rGO in the carbon matrix [9].

Surface Chemistry of the Samples
The effect of the GO on the surface composition was studied by XPS ( Table 2). The successful removal of the metal impurities was confirmed by SEM-supported ED spectroscopy ( Figure S3). As expected, the O content of the incorporated GO affected not only the porous texture but also the surface chemistry of the annealed samples. The complexity of the thermal degradation and high temperature chemical reactions are responsible for the change in the overall C or O content and the trend of the heteroatom/carbon ratios. The composite photoelectron lines (C1s, O1s, N1s, and S2p) were decomposed into different chemical states. The self-oxidation altered the surface composition as well as the chemical environment. Due to the complexity of the samples, however, the chemical states could not be resolved exactly. The suggested states referred to a broad range of chemical bonds. The shape of the C1 (sp 2 C) component was asymmetric, as in pure and low heteroatom-containing graphite and graphene. The tail of the asymmetric peak overlapped with the peaks of the C-O(N) components; using this line-shape also gave a satisfactory oxygen + nitrogen balance (the ratio of the measured and calculated concentrations) for all samples. Figure 4 presents a typical decomposition scheme of the C1s, O1s, N1s, and S2p spectral regions, while Tables 3 and 4 show the quantitative results of the decomposition, together with the binding energy ranges of the components and their assignations to various chemical states.
gether with the binding energy ranges of the components and their assignations to various chemical states.
The decomposition of the C1s signal revealed three different carbon species. The most abundant was the sp 2 form (C1). The presence of the π-π * shake-up satellite also confirmed the assignation of this component. The concentration of C1 was slightly higher in the rGO-containing samples. The other components could be assigned to carbon atoms connected to one (C2) or two (C3) heteroatoms.  Table 3. Decomposition of the C1s and O1s regions of the photoelectron spectra: binding energy ranges, chemical state assignations, and surface compositions (atomic %).

C1s
O1s   The decomposition of the C1s signal revealed three different carbon species. The most abundant was the sp 2 form (C1). The presence of the π-π * shake-up satellite also confirmed the assignation of this component. The concentration of C1 was slightly higher in the rGO-containing samples. The other components could be assigned to carbon atoms connected to one (C2) or two (C3) heteroatoms.
The O1 and O2 oxygen components were found in all samples; they could be assigned to the S-O and to various C-O bonds, respectively. The O3 component, present in the high rGO-containing samples, could be connected to either highly oxidized carbon (anhydride or carbonate). Owing to their low concentration, the corresponding carbon species could not be resolved in the spectra. Nitrogen was also present in three chemical states. N1 was the simple C-N bond, while N2 was bonded to carbon with further heteroatoms. The quantities of these nitrogens were commensurable. The low-concentration high-binding energy N3 component was a kind of quaternary ammonium. Due to the synthesis procedure, it is assumed that the nitrogen atoms were homogeneously distributed in the matrix and also decorated the walls of the micropores and ultra-micropores [102]. The sulfur content, present as C-S sulfide (S1) and in oxidized forms (S2; sulfone, sulfate), was the least affected by the treatments, but it should be noted that all the sulfur in the CAGO50 sample was sulfidic, while in the others, the contribution of the S2 forms reached 20-30%.
The FTIR measurements give complementary information about the chemical structures of the samples with diverse GO content ( Figure S4, Table 4). From the rich spectra, three main regions were considered to monitor the chemical changes in the bulk [103,104]. All spectra displayed a strong adsorption band at 1554 cm −1 , assigned to the -C=C-skeletal vibration of aromatic regions. The 1749 cm −1 band of the GO-loaded samples implied the presence of carbonyl groups (-C=O), stemming from aldehydes, ketones, esters, and carboxylic acids. The band at 1372 cm −1 , which was absent in the starting material, was ascribed to the -OH in-plane bending of phenols. If we compare the amount of the latter two signals to the corresponding carbon skeleton peak at 1554 cm −1 , we can conclude that the ratio of oxygen-containing groups increases with the increasing GO content.
In summary, addition of GO to the ι-carrageenan-urea carbon xerogel resulted in multifaceted modification of the textural and chemical properties of the samples. The thermal decomposition of GO enhanced self-activation and facilitated the formation of micropores which, for the CAGO50 and CAGO100 samples, also resulted in a considerable increase in the apparent surface area. With 5 wt% GO, too much oxygen was probably introduced into the matrix which destroyed the finer porosity. The annealing at 1000 • C removed part of the heteroatoms but left behind a high number of defects. It is expected that these new active sites (N-S-C) "constructed" from edged thiophene S, graphitic N, and pentagon defects, as well as the 3D porous architecture, result in an enhanced electrocatalytic activity [17].

Effect of GO on ORR Performance
Like other groups, we also found that the incorporation of rGO improved the electric conductivity of the cryogels [73,105]. In our samples, the improvement was proportional to the amount of GO added, in spite of the lower apparent surface area of CAGO200 ( Figure S5).
The ORR performance of the samples was investigated under conditions similar to the alkaline anion exchange membrane fuel cells (AAEMFCs). The linear sweep voltammograms (LSVs) of the four samples are plotted in Figure 5. The high potential regions revealed not only the slight differences between the various samples but also their alteration during the consecutive cycles. CAGO200 was found to be particularly unstable. The addition of GO gradually shifted the onset potential. CAGO50 exhibited the highest half-wave potential, which gradually decreased back to the value observed with the GO-free CA carbon. In Table 5, we compare the onset and half-wave potentials (E 1/2 ) of these curves with the results of other groups. mograms (LSVs) of the four samples are plotted in Figure 5. The high potential regions revealed not only the slight differences between the various samples but also their alteration during the consecutive cycles. CAGO200 was found to be particularly unstable. The addition of GO gradually shifted the onset potential. CAGO50 exhibited the highest halfwave potential, which gradually decreased back to the value observed with the GO-free CA carbon. In Table 5, we compare the onset and half-wave potentials (E1/2) of these curves with the results of other groups.  The Koutecky-Levich (KL) equation [114] was used to estimate the correlation between the current density j and the rotation rate ω where j k is the kinetic current density, j lim is the limiting diffusion current density, n is the number of electrons transferred in ORR per oxygen molecule, F is the Faraday constant, D is the diffusion coefficient of oxygen in the electrolyte, ν is the kinematic viscosity of the electrolyte, and C is the concentration of oxygen in the electrolyte. The KL plots of CAGO50 ( Figure 6) displayed a slope closest to the theoretical curve corresponding to the 4e − mechanism, i.e., the incorporation of 1.25% GO improved the performance, whereas further GO had the opposite impact. At higher GO loadings, the lines became steeper, nearer to the theoretical 2e − slope. Based on these findings from our samples, CAGO50 seems to have the optimum pore texture and chemistry for ORR under the investigated conditions. Although the electrical conductivity increased linearly with the increasing GO loading, no similar change in ORR was apparent. This agrees well with the work by Ramos-Fernández et al., who also observed that increasing GO content improved the electric conductivity of carbon xerogels, but enhanced porosity may spoil the electrochemical performance by disrupting the rGO network throughout the carbon matrix [73].
nearer to the theoretical 2e slope. Based on these findings from our samples, CAGO50 seems to have the optimum pore texture and chemistry for ORR under the investigated conditions. Although the electrical conductivity increased linearly with the increasing GO loading, no similar change in ORR was apparent. This agrees well with the work by Ramos-Fernández et al., who also observed that increasing GO content improved the electric conductivity of carbon xerogels, but enhanced porosity may spoil the electrochemical performance by disrupting the rGO network throughout the carbon matrix [73]. The (in)stability of the electrodes was further explored with CV responses measured after a series of treatments simulating application conditions. Based on the KL evaluation the GO-free CA and the CAGO50 and CAGO200, showing the lowest and highest slopes were selected for stability studies. The fresh electrodes made from the carbons, after running a CV cycle (curve 1) were exposed to further 50 cycles in argon-treated oxygen-free KOH electrolyte and then to another 50 cycles in oxygen-saturated KOH. Finally, an ORR test was run (sweep rate: 5 mV/s, rotation rates: 400, 625, 900, 1250 rpm). The voltammograms, curves 2, 3, and 4, obtained after each set of treatments are compared in Figure 7. (The effect of the electrode loading is discussed in the Supplementary Material) The slightly distorted rectangular shapes may stem from the pseudo-faradaic reactions of the corresponding functional groups [115]. The width of the hysteresis is not correlated with the BET surface area (Table 1), but the presence of the GO apparently reduces the gravimetric capacitance. The coincidence of curves 1 and 2 in CA carbon (Figure 7) revealed that this sample was stable under oxygen-free conditions. After an additional 50 cycles in oxygen-saturated electrolyte (curve 3), the widening of the CV curve implies that the electrocatalytically active surface area increased in the presence of oxygen. The consecutive four ORR cycles further modified the electrode surface (curve 4). The (in)stability of the electrodes was further explored with CV responses measured after a series of treatments simulating application conditions. Based on the KL evaluation the GO-free CA and the CAGO50 and CAGO200, showing the lowest and highest slopes were selected for stability studies. The fresh electrodes made from the carbons, after running a CV cycle (curve 1) were exposed to further 50 cycles in argon-treated oxygen-free KOH electrolyte and then to another 50 cycles in oxygen-saturated KOH. Finally, an ORR test was run (sweep rate: 5 mV/s, rotation rates: 400, 625, 900, 1250 rpm). The voltammograms, curves 2, 3, and 4, obtained after each set of treatments are compared in Figure 7. (The effect of the electrode loading is discussed in the Supplementary Material) The slightly distorted rectangular shapes may stem from the pseudo-faradaic reactions of the corresponding functional groups [115]. The width of the hysteresis is not correlated with the BET surface area (Table 1), but the presence of the GO apparently reduces the gravimetric capacitance. The coincidence of curves 1 and 2 in CA carbon (Figure 7) revealed that this sample was stable under oxygen-free conditions. After an additional 50 cycles in oxygen-saturated electrolyte (curve 3), the widening of the CV curve implies that the electrocatalytically active surface area increased in the presence of oxygen. The consecutive four ORR cycles further modified the electrode surface (curve 4).
The performance of the CAGO50 carbon was very similar, but it seemed somehow even less stable under oxygenic conditions. The CAGO200 electrode was unstable already in oxygen-free conditions (curves 1 and 2). The change in the surface was even more pronounced after being cycled in the oxygen-saturated electrolyte, which apparently stabilized the surface: no further changes were seen after these ORR cycles (curve 4). These observations indicate that the electron transfer mechanism deduced from the KL plot can only be considered as an approximation. The kinetics of the reduction reaction in the presence of the heteroatoms is indeed more sophisticated and further influenced by the rGO. The chemical species locally evolving in the ORR may initiate various processes including restructuring, activation, degradation, etc., of the electrode surface. The plots in Figure 7 clearly show that even if the electrode was stabilized under oxygen-free conditions, it might later undergo destabilization processes during the ORR. Nevertheless, the widening of the CV curves marked an increase in the electrocatalytically active surface area. The deviation of the corresponding curves confirmed that the addition of GO to the precursor affected the stability of the electrodes. As expected from the KL plot, the CAGO200 sample was more sensitive to the oxidative conditions, but after stabilization this sample exhibited the highest capacitance. The trend of the gravimetric capacitance in this figure was different from the trend expected from the apparent surface area (1070, 1479, and 933 m 2 /g for CA, CAGO50, and CAGO200, respectively), implying that the accessibility of the nitrogen molecules is different from that of the species involved in the electrocatalytic processes [116]. (3) after an additional 50 cycles in oxygen-saturated electrolytes; (4) after an additional ORR test with rotating electrode (sweep rate: 5 mV/s, four rotation rates). The column graph compares the gravimetric capacitances of the carbons corresponding to CV 1 (black), 2 (red), 3 (green), and 4 (blue), respectively.
The performance of the CAGO50 carbon was very similar, but it seemed somehow even less stable under oxygenic conditions. The CAGO200 electrode was unstable already in oxygen-free conditions (curves 1 and 2). The change in the surface was even more pronounced after being cycled in the oxygen-saturated electrolyte, which apparently stabilized the surface: no further changes were seen after these ORR cycles (curve 4). These observations indicate that the electron transfer mechanism deduced from the KL plot can only be considered as an approximation. The kinetics of the reduction reaction in the presence of the heteroatoms is indeed more sophisticated and further influenced by the rGO. The chemical species locally evolving in the ORR may initiate various processes including restructuring, activation, degradation, etc., of the electrode surface. The plots in Figure 7 clearly show that even if the electrode was stabilized under oxygen-free conditions, it might later undergo destabilization processes during the ORR. Nevertheless, the widening of the CV curves marked an increase in the electrocatalytically active surface area. The deviation of the corresponding curves confirmed that the addition of GO to the precursor affected the stability of the electrodes. As expected from the KL plot, the CAGO200 sample was more sensitive to the oxidative conditions, but after stabilization this sample exhibited the highest capacitance. The trend of the gravimetric capacitance in this figure was different from the trend expected from the apparent surface area (1070, 1479, and 933 m 2 /g (3) after an additional 50 cycles in oxygen-saturated electrolytes; (4) after an additional ORR test with rotating electrode (sweep rate: 5 mV/s, four rotation rates). The column graph compares the gravimetric capacitances of the carbons corresponding to CV 1 (black), 2 (red), 3 (green), and 4 (blue), respectively.

Effect of the GO on Li-Ion Battery Application
In order to reveal the effect of GO on the performance of the carbon cryogel as anode material in Li-ion batteries, galvanostatic charge-discharge cycling tests were performed. In the half-cell setup, the CA, CAGO50, CAGO100, or CAGO200 sample was the working electrode, and the Li metal disc served as the counter and reference electrode, i.e., no further conductive agent was used. In the galvanostatic test, a constant current I is applied to the cell and the potential is evaluated as a function of time t. The specific capacity C s during complete charge-discharge is given by The nature of the charge-discharge plots largely depends on the surface chemistry (microstructure, porosity, heteroatom content, etc.) and the adsorption capacity of the carbon host. Figure 8 shows the cell voltage vs. the specific capacity charge-discharge cycling plots of the samples under investigation. While highly ordered graphitized carbon anodes show a steep charge-discharge curve [117], sloping shapes similar to Figure 8 are distinctive for amorphous biomass-derived carbon anodes in Li-ion battery applica-tions [33,118,119]. The broad peaks in the XRD patterns ( Figure S7) reveal that the carbon cryogels studied in this work were typically amorphous; therefore, charge-discharge curves with this shape were expected, i.e., the CA, CAGO50, CAGO100, and CAGO200 have the potential to store the Li-ions gradually and consistently over their operating potential. A comparable charge-discharge plateau was measured also on RF polymer-based carbon xerogels [120]. The latter microporous carbon xerogels were successfully tested as anode material in lithium-ion cells exhibiting a specific capacity 288 mAh/g and retaining 97% of the initial capacity after 100 charge-discharge cycles.
show a steep charge-discharge curve [117], sloping shapes similar to Figure 8 are disti tive for amorphous biomass-derived carbon anodes in Li-ion battery applicatio [33,118,119]. The broad peaks in the XRD patterns ( Figure S7) reveal that the carbon cr gels studied in this work were typically amorphous; therefore, charge-discharge cur with this shape were expected, i.e., the CA, CAGO50, CAGO100, and CAGO200 have potential to store the Li-ions gradually and consistently over their operating potential comparable charge-discharge plateau was measured also on RF polymer-based carb xerogels [120]. The latter microporous carbon xerogels were successfully tested as ano material in lithium-ion cells exhibiting a specific capacity 288 mAh/g and retaining 9 of the initial capacity after 100 charge-discharge cycles.
The cyclability of the carbons is presented here as the total charge-discharge capac (or capacity retention) as a function of the number of cycles (Figure 9). The undoped c bon, CA showed an initial discharge capacity of 332 mAh/g at a constant current dens 100 mA/g. The cell showed good capacity retention and high coulombic efficiency (92. and 99.7%, respectively) even after 50 charge-discharge cycles. The cyclability of the carbons is presented here as the total charge-discharge capacity (or capacity retention) as a function of the number of cycles (Figure 9). The undoped carbon, CA showed an initial discharge capacity of 332 mAh/g at a constant current density 100 mA/g. The cell showed good capacity retention and high coulombic efficiency (92.9% and 99.7%, respectively) even after 50 charge-discharge cycles.
For comparison, 218 mAh/g at a current density 0.372 A/g was reported by Wang et al. [6]. Among the GO doped carbons, CAGO50 exhibited the best performance in terms of the cycling stability. The cell showed an initial discharge capacity of 313 mAh/g and a highly reversible capacity of 297 mAh/g after 50 cycles, corresponding to a 95.1% capacity retention. This good stability agreed well with the impedance analysis results: this carbon had the lowest impedance. The enhanced performance of CAGO50 compared to the pristine CA can be explained in terms of the surface area and porosity [121]. The higher surface area of CAGO50 improved the electrolyte accessibility and ion diffusion in the pores of the electrode, which in turn resulted in a higher coulombic efficiency and cell stability [30,121]. The charge-discharge plateau recorded in the first and 50th cycle was almost identical, confirming the resilient stability of the electrode [122]. Following this, the CAGO200 showed an initial discharge capacity of 240 mAh/g and a capacity retention of 93.6% after 50 cycles (Figure 9).

Figure 8.
Galvanostatic charge-discharge profiles of CA, CAGO50, CAGO100, and CAGO200 at a constant current density of 100 mA/g, cycle 1 black, cycle 10 red, cycle 20 olive, cycle 30 blue, cycle 40 dark yellow, and cycle 50 cyan. Figure 9. Long-term cycling performance of CA, CAGO50, CAGO100, and CAGO200 at a constant current density of 100 mA/g. Capacity retention black, Coulombic efficiency red.
For comparison, 218 mAh/g at a current density 0.372 A/g was reported by Wang et al. [6]. Among the GO doped carbons, CAGO50 exhibited the best performance in terms of the cycling stability. The cell showed an initial discharge capacity of 313 mAh/g and a highly reversible capacity of 297 mAh/g after 50 cycles, corresponding to a 95.1% capacity retention. This good stability agreed well with the impedance analysis results: this carbon The overall electrochemical performance of CAGO200 in the Li-ion battery tests was better than the undoped CA and slightly lower than the CAGO50. The CAGO100 presented the poorest electrochemical performance in spite of having the largest surface area. This apparent surface area, however, was also combined with a larger pore volume in the wider pore size region, which seemed to decrease the initial coulombic efficiency and reduce the reversible capacity of the battery [30]. The initial discharge capacity of CAGO100 was 255 mAh/g, but the cell maintained its stability poorly. The low retention capacity (67.2% after 50 cycles) might be due to a kinetic limitation in the porous electrode [121]. The compatibility of the pore and the Li + ion sizes is also responsible for the charge capacity. CAGO100 and CAGO200 showed the highest and lowest pore volume, respectively, both in the micro-and mesopore regions. On the one hand, wider pores may not be completely filled by the Li + ions, and on the other, low porosity limits the entrance of Li + containing electrolyte into the pores. CAGO50, which seems to have the optimum pore size compatible with Li + and is capable of effectively accommodating the Li + ions, exhibited the highest overall capacity among the GO-doped samples.
As shown in Figure 9, the coulombic efficiency of all the samples was consistently maintained above 98.5% indicating the efficient Li-ion insertion and reinsertion potential of these carbons during the charge-discharge process [122]. CA and CAGO100 displayed lower coulombic efficiency than CAGO50 and CAGO200 (>99.3%) indicating that the addition of GO was clearly beneficial for more stable electrochemical performance in Li-ion batteries, but the GO doping level must be carefully tuned.
It may be concluded that the carbons investigated here can be utilized as potential anodes in Li-ion batteries. According to our results, the pore structure and the chemical composition corresponding to CAGO50 proved to be the best to achieve an acceptably high reversible discharge capacity and electrode stability under prolonged cycling. The electrochemical performance of these carbon electrodes is comparable to the previously reported results on carbon electrode anodes in Li-ion batteries [123][124][125][126][127].

Conclusions
Addition of GO to the precursor biopolymer successfully increased the surface area and the porosity of the carbon material up to 2.5 wt% GO (up to S BET 1780 m 2 /g and V 0.98 1.7 cm 3 /g), but the high amount of gases released when 5 wt% GO was added destroyed the finer pores. The oxidative species evolving from the GO during the heat treatments influenced the surface chemistry in a complex way, but sufficient O, N, and S heteroatoms were preserved in various chemical forms. The highest concentration of sp 2 carbons was detected on the surface of sample CAGO50 (1.25 wt% added GO). Although the volatile thermal degradation products of the GO affected the porous texture and the surface chemistry in a subtle way, the residual rGO resulted in a monotonically increasing electrical conductivity, in spite of the notable drop in the surface area. The rGO affected the electron transfer mechanism and compromised the stability of the electrodes in ORR tests. The CAGO50 sample displayed the advantageous 4emechanism. The CAGO200 (5 wt% GO) sample was the most sensitive under oxidative conditions, but after conditioning, it exhibited the highest capacitance in ORR. In the Li-ion battery tests, the coulombic efficiency of all the samples was consistently above 98%. The cell made from CAGO50 showed an initial discharge capacity of 313 mAh/g, 95.1% capacity retention, and 99.3% coulombic efficiency after 50 charge-discharge cycles. Despite having the largest surface area sample, the CAGO100 performed more poorly than the CAGO50 in both application tests. This is a clear indication that concerted fine-tuning of the pore morphology and surface chemistry is required for both advanced electrocatalytic and electrochemical performance. The asprepared GO-doped carbon samples with their ultra-micropores and high surface area performed well and could be promising electrode materials for both applications.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/nano13162364/s1, Figure S1: SEM image of the CA matrix showing nanopores of width 2-3 nm; Figure S2: Raman spectra of the annealed carbon samples (a); deconvoluted Raman spectrum of CAGO50 (b); Figure S3: SEM/EDS spectra of the samples (at magnification ×100); Figure S4: FTIR spectra of the carbon samples; Figure S5: Influence of added GO on the specific conductivity and apparent surface area; Figure S6 Table S1: Yield of the synthesis steps (%); Table S2: Raman shifts from the deconvoluted spectra ( Figure S2b) [cm −1 ]; Table S3: Selected data deduced from the deconvoluted Raman spectra ( Figure S2b); Table S4: Ratio of the various carbon species based on FTIR spectra ( Figure S4). Ref. [128] is cited in the supplementary material in Table  S2.

Data Availability Statement:
The data presented in this study are available on request.