Received: 04 Jan 2019
Revised: 25 Jan 2019
Accepted: 27 Jan 2019
Published online: 27 Jan 2019
Mingxiang Hu,1 Qian Lv1 and Ruitao Lv1,2*
1 State Key Laboratory of New Ceramics and Fine Processing, School of Materials Science and Engineering, Tsinghua University, Beijing 100084, China
2 Key Laboratory of Advanced Materials (MOE), School of Materials Science and Engineering, Tsinghua University, Beijing 100084, China
*E-mail: lvruitao@tsinghua.edu.cn
Heteroatom doping is an efficient way to tailor the electrochemical properties of carbon-based energy storage materials. However, synthesizing nitrogen-doped carbon materials with a high doping level (e.g. over 10.0 at.% N content) and controllable doping contents is very challenging. Herein, a tablet-sintering method using two nitrogen sources (ammonia and urea) is proposed to realize the regulation of nitrogen contents in the carbon materials from 4.8 at.% to almost 12.0 at.%. Optimized sample with 10.1 at.% nitrogen content demonstrates excellent pseudocapacitive behaviors with 104.1 F g-1 and 176.3 F g-1 for the lithium- and sodium- ion capacitorsin organic electrolytes at 1 A g-1. Moreover, highly doped nitrogen could promote the cycling stability of carbon-based materials high voltage (up to 5.0 V) in the ionic liquid electrolyte. This may open a new avenue for the chemical doping of carbon-based materials and their applications in high-performance electrochemical devices.
2.2 Sample characterizations
The crystalline structures of obtained samples were characterized by Xray diffraction (XRD, Rigaku, D/Max-2500/PC) measurement using Cu-Kα radiation. Raman spectra (LabRAM,HR Evolution) were used to characterize the defects in carbon-based materials. The morphologies of as-synthesized samples were investigated with a field-emission scanning electron microscopy (FE-SEM, Hitachi S4800) and a transmission electron microscopy (TEM, JEOL 2010), both of which were equipped with energy dispersive spectrometers (EDS). The specific surface area (SSA) and pore structures were characterized by the nitrogen adsorption/desorption at 77 K (Belsorp-Max/Mini). The pore size distributions were calculated by the density functional theory (DFT) method. X-ray photoelectron spectroscopy (XPS) data were collected with an ESCLab 250X followed by calibrating the binding energy with C1s (284.8 eV).

Scheme 1. Schematic illustration of a bi-source doping synthesis of nitrogen-doped graphene oxide. Graphene oxide (GO) and urea powders are mixed and pressed into a tablet, then sintered in NH3 atmosphere. Here urea powders (solid) and ammonia (gas) are used as the “internal” and “external” nitrogen sources for the controllable doping of GO (solid), respectively.
2.3 Electrochemical tests
Electrochemical properties of as-obtained samples were tested in the stainless-steel 2032-coin cells. The electrodes with ~1.2 mg cm-2 loading mass were prepared by coating the slurries containing the active materials (90 wt.%) and polyvinylidene diflouride (PVDF) (10 wt.%) onto the surface of aluminum foils and dried at 120 ℃ in vacuum overnight. Half-cells were assembled in an argon-filled glove box with metal Na/Li foils as the counter electrodes. Electrolytes for LICs and NICs are LiPF6 (1 M) and NaClO4 (1 M) in the ethylene carbonate (EC) and dimethyl carbonate (DMC) (volume ratio of EC: DMC=1:1), respectively. The NICs were tested at the potential range of 1.5-4.2 V (vs. Na+/Na). Because of the lower standard potential of Li, LICs were tested at the potential ranging from 1.8 V to 4.5 V (vs. Li+/Li) on Land battery testing systems. Meanwhile, the ionic liquid electrolyte was obtained by dissolving lithium bis(trifluoromethylsulphonyl)imide (LiTFSI) in 1-ethyl-3-methylimidazolium-bis(trifluoromethylsulfonyl)imide) (EMITFSI). It was used to investigate the electrochemical properties under high-voltage (5 V). The electrochemical properties of as-obtained samples in aqueous solution (6 M KOH as electrolyte) were also investigated systematically. To elucidate the energy storage mechanism, the electrochemical impedance spectroscopy (EIS) with the frequency from 100 kHz to 10 mHz and cyclic voltammetry (CV) techniques at related voltage windows were conducted on an electrochemical workstation (EC-Lab, Bio-shop).

Fig. 2 (a) Inductively coupled plasma mass spectrometry (ICP-MS) analysis of C, N and O elements in NGO, G1U0 and G1U2 samples. (b) Nitrogen and oxygen elemental contents in different samples. (c) Schematic diagram of nitrogen doped graphene; Pyridinic N (N-6), pyrrolic N (N-5) and graphitic N (N Q). (d) High-resolution XPS spectra of N 1s in the G1U2 samples.
The nitrogen doping levels of G1U2, G1U0 and NGO by ICP-MS analysis are 10.4 at.%, 7.1 at.% and 4.6 at.% (Fig. 2 (a)), which are in agreement with the XPS analysis. As demonstrated in Fig. 2 (b), the nitrogen contents will gradually increase with the increased amounts of urea from 4.8 at.% to almost 12.0 at.%. Meanwhile, the surface oxygencontaining groups are also partly removed with the increased urea amounts as shown in Fig. 2 (b). By deconvoluting the C1s peaks (Figs. 3(a-c)), the contents of C=C species calculated based on the integral area ratios are 73.2 %, 67.1 % and 72.8 % in G1U2, G1U0 and NGO, respectively. The increased contents of C=C bonds will increase the conductivity of samples thus lead to better electrochemical properties.26, 32 High-resolution N1s spectra of G1U2, G1U0 and NGO are shown in Figs. 3 (d-f). Different nitrogen species are illustrated in Fig. 2 (c), which are mainly the pyridinic nitrogen (N-6), pyrrolic nitrogen (N-5) and graphitic nitrogen (N-Q). Bandosz et al39 reported that the pyridinic nitrogen and pyrrolic nitrogen will contribute to an enhanced capacitance due to their positive charges compared with graphitic nitrogen. The integral calculations of peak areas confirm that the total contents of pyridinic nitrogen and pyrrolic nitrogen in G1U2, G1U0 and NGO are 93.75 %, 84.52 % and 77.62 %, respectively. The ratios of pyridinic nitrogen and pyrrolic nitrogen will increase with the increase of the amount of solid nitrogen source (urea), which will be important for the enhancement of the electrochemical performance of as-obtained samples. As shown in Fig. 2(d), the ratios of pyridinic nitrogen (N-6), pyrrolic nitrogen (N-5) and graphitic nitrogen (N-Q) in G1U2 are 6.25 %, 44.49 % and 49.26 %, respectively. It is expected that the G1U2 could deliver a competitive performance as electrode in capacitors. As-synthesized sample was furtherly investigated as energy storage materials in LICs/NICs and aqueous supercapacitors. Firstly, the electrochemical properties of all these samples are investigated as the cathodes of NICs in the potential range from 1.5 V to 4.2 V at l A g-1. As demonstrated in Fig. S6, the capacity of NGO is 42 mAh g-1, much higher than that (17 mAh g-1) of G1U0. Tablet-sintering samples show that the capacities increase first and then decrease as shown in Fig. S6. As the cathode of NIC, the G1U2 sample demonstrates the best electrochemical capacity among all these doped GO samples with ~82 mAh g-1. Galvanostatic charge-discharge (GCD) curves of G1U2, G1U0 and NGO as the cathodes of LICs and NICs in organic electrolyte systems are shown in Fig. 4. In Fig. 4(a), G1U2 could deliver about 104.1 F g-1 specific capacitance at 1 A g-1 in the voltage window ranging from 1.8 V to 4.5 V in LICs. Meanwhile, both G1U0 and NGO only deliver 36.1 F g-1 capacitances in LICs, indicating the enhanced electrochemical capacitances by high content of nitrogen doping. As shown in Fig. 4(b), the cyclic voltammetry (CV) curves of G1U2, G1U0 and NGO electrodes in LICs display rectangular shapes, which are the typical EDLC behaviors. Besides, several humps exist in the CV curves (Fig. 4(b)), demonstrating the surface redox reactions of pseudocapacitances. The Nyquist plots (Fig. 4(c)) of half-cell LICs with open circuit states exhibit 45o slopes in the low frequency region, agreeing with the EDLC behaviors of carbon-based materials. The solution resistances (R ) of G1U2, G1U0 and NGO electrodes in LICs are 4.5 Ω, 4.7 Ω and 12.1 Ω, respectively, indicative of prior contacts between the electrolyte and electrodes. Different from the Ohmic contact resistances, the charge-transfer resistances (R) are usually the ct measurement of ion conductivity, whose values can be calculated by the diameters of the semicircles in Nyquist plots in high-frequency region. Clearly, R of G1U2 is about 10 Ω, much smaller than those of G1U0 ct (102 Ω) and NGO (53 Ω). Performances of as-obtained samples were also tested in the NIC systems using the same methods and the results are shown in Fig. 4(d). The GCD measurements of G1U2, G1U0 and NGO are carried out in the potentials ranging from 1.5 V to 4.2 V which is different from LICs due to the higher standard hydrogen electrode (SHE) potential of Li. It can becalculated that the capacitances of G1U2, G1U0 and NGO are 176.3 F g-1, 27.9 F g-1 and 77.8 F g-1, respectively. The capacitance of G1U2 is competitive compared with some recent works of carbon-based materials in capacitors.40-43

Fig. 3 High-resolution X-ray photoelectron spectroscopy (XPS) spectra of G1U2; G1U0 and NGO samples: (a-c) C1s; (d-f) N1s.
Different from the performance in LICs, the electrochemical performance of NGO is better than that of G1U0 in NICs, which may come from the inadequate contact between the surface of GO powders and ammonia molecules after a compaction process of G1U0. Moreover, this phenomenon could be explained from the EIS results in Fig. 4(f), which reveals the same trends as the Nyquist plots of LICs (Fig. 4(c)). Correspondingly, CV curves of G1U2, G1U0 and NGO in NICs shown in Fig. 4(e) display the same characteristics as the CV curves of LICs. The CV curve with rectangle shape indicates the existence of EDLCs while the humps overlaid the curves are from the pseudocapacitive reactions. Comparing the G1U2 electrochemical results of NICs and LICs in Fig. S7, it is found that the capacitance of NIC is about 60 % larger than that of LIC. Because of the higher standard hydrogen electrode (SHE) potential (Li+/Li: -3.0 V; Na+/Na: -2.7 V) of Li, the voltage potentials of LICs and NICs are set to be the same based on Li+/Li as shown in Figs. S7 (a, b). As EIS plots shown in Fig. S7(c), the charge transfer resistances of G1U2 in NIC and LIC are about 10, no obvious differences in both devices. After 3500 cycles and 2000 cycles in LICs and NICs, G1U2 could still deliver ~60 mAh g-1 and ~80 mAh g-1 capacities under 1 A g-1, respectively (Fig. S8). No obvious morphology changes can be found in the sample before and after the LIC and NIC tests, as shown in Fig. S9.

Fig. 4 Electrochemical properties of G1U2, G1U0 and NGO samples. (a) Comparison of Galvanostatic charge-discharge (GCD) curves as the cathodes of lithium ion capacitors (LICs) at 1 A g-1; (b) CV curves in LICs at as can rate of 50 mV/s; (c) Comparison of EIS curves of as-obtained samples; (d) Comparison of GCD curves as the cathodes of sodium ion capacitors (NICs) at 1 A g-1; (e) Cyclic voltammetry (CV) curves in NICs at the scan rate of 50 mV/s; (f) Comparison of Electrochemical impedance spectroscopy (EIS) curves of as-obtained samples in NICs.
To investigate the capacity retention and the time constants of G1U2 in both LIC and NIC, the capacitors were firstly charged under 1 A g-1 and then discharged under various current densities as exhibited in Fig. 5. The time constant (τ ) was calculated using based on the 0 following formula:

where R is the resistance of electrode, C is the capacitance, U is the voltage drop, and I is the current density. As demonstrated in Figs. 5(a) and (c), it can be calculated that the time constants of LIC and NIC are about 4.16 s and 10.58 s, respectively. In a stable charge and discharge process, the charge and discharge capacities are mutually influenced by each other. High current discharge will cause a decrease in charge capacity at low currents as demonstrated in Figs. 5(b) and (d). The rate capacitances of G1U2 working as cathode for NIC are compared with some state-of-the-art carbon-based NIC electrodes as exhibited in Fig. 5(e).42,44-47 It can be seen that G1U2 could not only deliver high capacitance in low current densities, but also possess high capacity retention under high current densities.

Fig. 5 (a), (c) GCD profiles of G1U2 in LIC and NIC with 1 A g-1 charge current densities and different discharge current densities; (b), (d) Summary of charge and discharge capacities of samples at various current densities and their related Coulombic efficiencies. (e) Comparison of rate capacitance of G1U2 in NIC with other carbon-based materials. NPHCS: nitrogen-doped microporous hard carbons;44 SCN A: peanut skin derived carbon nanosheets;45 NCCF-Acid: Nano cellular carbon foams tested in acid.47
The energy density (E) of the hybrid device can be evaluated based on the formula as below:

where U represents the working potential, I represents the working current, and m refers to the mass of active materials. Obviously, E is proportional to the working potential (U), indicating that the increase of U would sharply enhance the electrochemical properties of the obtained devices. Herein, optimized samples G1U2 were investigated in the ionic liquid in the voltage ranging from 2.5 V to 5.0 V as shown in Fig. 6. When the potential reaches 5.0 V, the anions intercalation would take place. When working at 1 A g-1, the Coulombic efficiency of G1U2 in ionic liquid is only about 87 %, which is mainly caused by the low Coulombic efficiency of anions intercalation (11.8 %) as marked in Fig. 6(a). Beltrop et al.48 proposed that the size of anions and the cut-off potential will obviously influence the Coulombic efficiency of batteries. Moreover, both the electrochemical exfoliation and poor wetting ability between electrolyte and electrode will lead to the ultralow Coulombic efficiency in this system.49, 50 G1U2 could deliver ~134.7 F g-1, 78.9 F g-1 and 47.8 F g-1capacities at 1 A g-1, 2 A g-1 and 4 A g-1 as demonstrated in Figs. 6(b) and (c), indicating the poor rate performance in ionic liquid. Meanwhile, the kinetics of G1U2 in ionic liquid were investigated using CV method under various scan rates as shown in Fig. 6(d). Obviously, the areas of CV curves are proportional to the scan rates as shown in Fig. 6(d). The peak currents in CV curves obey a power-law relationship with the scan rates as follows:

where i is the peak current, v is the scan rate, and b values could be used to evaluate the reaction process of electrode. When b=0.5, it refers to a typical Faradaic intercalation process controlled by semi-infinite linear diffusion; when b=1, it represents a capacitive charge storage without diffusion-control mechanism.11, 51 As shown in Fig. 6(e), when working at 3.5 V in ionic liquid, the b values of G1U2 in oxidization and reduction processes are 0.74 and 0.68, revealing the pseudocapacitive behaviors of G1U2 sample. The electrochemical properties of G1U2 in ionic liquid were furtherly investigated by EIS analysis as shown in Fig. 6(f). Nyquist plots of G1U2 demonstrate that the solution resistance (Rs) in this system is about 21 Ω, much bigger than that in LIC (4.5 Ω) and NIC (9.0 Ω) systems. In the inset of Fig.6(f), the semicircles in high frequency region display a flat shape which are caused by the rough surface of highly nitrogen-doped GO, consistent with the pseudocapacitance introduced surface reaction process.24, 52

Fig. 6 (a) GCD curve of G1U2 in lithium ionic liquid at 1 A g-1 with the working voltages between 2.5 V and 5.0 V; (b) GCD curves of G1U2 in lithium ionic liquid at 1, 2 and 4 A g-1; (c) Rate capacities of G1U2; (d) CV curves of G1U2 at various scan rates; (e) Determination of the b values for various at 3.5 V for G1U2 in the ionic liquid electrolyte; (f) EIS data of G1U2 in ionic liquid and its fitted curve.
Fig. 7 reveals the electrochemical properties of G1U2 in aqueous solution (6 M KOH) between -1.0 V and 0 V vs. Hg/HgO. CV curves under different scan rates (Fig. 7(a)) display typical EDLC behaviors with rectangular shapes. G1U2 could deliver about 303, 285, 280, 254 and 234 F g-1 capacitances at 1, 2, 4, 20 and 40 A g-1 current densities. Excellent cycling property of G1U2 can be also obtained as shown in Fig. 7(c). Compared with the initial GCD curve, there is only 12.5% decrease in capacitance after 3000 cycles in aqueous solution at a high current density of 4 A g-1. The capacitance under specific potential can be calculated from the slopes of the lines between the current density and the scan rates based on the CV curves in Fig. 7(a). As shown in Fig. 7(d), the capacitances at different potentials are almost the same at about 300 F g-1, corresponding with the results of GCD curves. Benefiting from the ultralow resistance in aqueous solution (Fig. S10), G1U2 exhibits excellent rate performances with the increased scan rates as shown in Fig. 7(e). Even under the scan rate of 400 mV/s, it could still deliver over 240 F g-1 capacitance. As shown in Fig. 7(f), compared with some recent works on carbon-based materials as electrodes for SCs,53-57 although some materials exhibit comparable performance to that of G1U2 at low current densities (e.g. 1, 2 A g-1), as-synthesized G1U2 demonstrates superior capacitance than the others at high current densities (e.g. 20, 40 A g-1), which indicates their excellent rate performance.

Fig. 7 (a) CV curves of G1U2 in aqueous solution under various scan rates; (b) Galvanostatic discharge/charge (GCD) profiles of G1U2 in aqueous solution at current densities from 1 to 40 A g-1; (c) GCD curves of G1U2 in the first and 3000th cycles under 4 A g-1 in aqueous solution; (d) The sweep rate dependence on currents of G1U2 to characterize charge storage mechanisms. (e) Rate capacitance of G1U2 at different scan rates in aqueous electrolyte; (f) Comparison of rate capacitance of G1U2 in supercapacitor (SC) with other carbon-based materials. NC700: Cross-linked, N-doped, pillared-layered porous carbons obtained at 700 oC;53 PZCC: NaH2PO4-ZnCl2-treated cellulose-derived carbon;57 PFC-700: Perilla frutescenleaf-derived carbon nanosheets;56 AC: activated carbon;55 UPCS C-700: Ultrathinporous carbon shell synthesized by carbonization at 700 oC.54
52. X. Wang, G. Sun, P. Routh, D. H. Kim, W. Huang and P. Chen, Chem. Soc. Rev., 2014, 43, 7067-7098.
53. W. Tian, H. Zhang, H. Sun, M. O. Tadé and S. Wang, Carbon, 2017, 118, 98-105.
54. W. Yang, W. Yang, F. Ding, L. Sang, Z. Ma and G. Shao, Carbon, 2017, 111, 419-427.
55. K. Fic, G. Lota, M. Meller and E. Frackowiak, Energy Environ. Sci., 2012, 5, 5842-5850.
56. B. Liu, Y. Liu, H. Chen, M. Yang and H. Li, J. Power Sources, 2017, 341, 309-317.
57. J. Yi, Y. Qing, C. Wu, Y. Zeng, Y. Wu, X. Lu and Y. Tong, J. Power Sources, 2017, 351, 130-137.