Hao Liu,
Xiao-Juan Liu,
Feng-Ying Dong and
Xin-Zhi Sun*
College of Chemistry and Pharmaceutical Sciences, Qingdao Agricultural University, Qingdao 266109, China. E-mail: xzsun@qau.edu.cn
First published on 25th September 2020
It is a great challenge to ideally integrate graphene with its unique two-dimensional (2D) and porous structure into the pseudocapacitive materials. In this paper, a simple technique, i.e. direct-laser-writing (DLW), was developed to fabricate microsupercapacitors (MSCs) with excellent electrochemical performance, marked as Ni–Co–S/laser induced graphene (LIG) that exhibit a high areal specific capacitance of 680 mF cm−2 at the current density of 1 mA cm−2. A symmetric MSC device was assembled using Ni–Co–S/LIG as a positive electrode and active carbon (AC) as the negative electrode, and exhibited a high areal energy density of 56.9 μW h cm−2 at the power density of 800 μW cm−2, and excellent cycling stability maintaining 89.6% of the areal specific capacitance after 8000 cycles. The synergistic effect of bimetallic Ni–Co–S and the LIG with the 2D structure results in the excellent electrochemical performance. This work demonstrates a method to integrate Ni–Co–S pseudocapacitive materials into porous graphene with a direct-laser-writing technique. The produced integrated materials possess high energy density that can be used in MSCs.
In recent years, graphene-based nanomaterials have attracted great attention for their unique two-dimensional (2D) and porous structure.13,14 However, commercialization of graphene-based nanomaterials is difficult due to its complicated, costly, and tedious process.15–17 For example, graphene films were produced by using laser scribing of hydrated graphene oxide (GO) film, which showed predominant electrochemical performance.18–21 But it is still difficult to synthesize GO into graphene film in large scale.22,23 In recent years, a simple, fast and direct method is developed for the preparation of 2D porous graphene film from polymers by using direct-laser-writing (DLW) technique. Since then, the produced so-called laser-induced-graphene (LIG)24 has been widely used in MSCs.25–29
The energy density of the supercapacitor device is determined by the specific capacitance and the voltage window (E = CV2/2).30 Therefore, in order to improve the energy storage, pseudocapacitive materials are loaded onto the graphene to offer pseudocapacitance from the redox reaction generated by metal ions. Recently, transition metal sulfide (TMS) electrode materials have attracted much attention due to their predominant pseudocapacity.31–33 For example, the fabrication of cobalt-nickel-sulfide as electrode materials has been applied in improving supercapacitor performance owing to fast redox reactions, higher conductivity, diverse valence state and availability.34 Pang and his group members state in detail the application of different sulfides in energy storage.35 Compared with nonmetallic sulfide, the combination of nickel and cobalt in bimetal sulfide gives a wide voltage window and improves the energy storage performance.36,37
In this study, we developed a simple and one-step method towards preparation of bimetal nickel-cobalt-sulfide composite from the precursor, and LIG from the carbon cloth by using direct-laser-writing (DLW) technique. The bimetal Ni–Co–S/LIG exhibits an outstanding areal specific capacitance of 680 mF cm−2 at the current density of 1 mA cm−2, which is superior to the substrate LIG electrode materials. In addition, the asymmetric MSCs were fabricated with Ni–Co–S/LIG as positive electrode and active carbon (AC) as negative electrode materials, respectively. Moreover, the asymmetric MSCs device can reach a maximum energy density of 56.9 μW h cm−2 at a power density of 800 μW cm−2, which exceeds most of the reported asymmetric MSCs containing nickel or cobalt or LIG materials.
X-ray photoelectron spectroscopy (XPS) was carried on a Thermo ESCALAB 250Xi device, which Al Kα X-ray beams was used as excitation source (hv = 1486.6 eV). All the spectra were adjusted with C1s peaks at 284.5 eV as reference.
The crystalline size (La) in “a” axis was obtained from the different ratio of the intensity of the G peak (IG) and D peak (ID) using Raman spectroscopic data. The equation following was used to calculate the value of La.24
(1) |
The specific areal (CA, in mF cm−2) capacitance of the negative electrode and device were obtained from the galvanostatic discharge graph at different densities according to the following formular.24
(2) |
Electrochemical impedance spectroscopy was performed using a sinusoidal signal with a frequency scope from 0.01 to 105 Hz with a magnitude of 10 mV. The solution of resistance (Rs) was got with the voltage drop (in the beginning of the discharge curve) using the following equation.
(3) |
The asymmetric MSCs device was built with Ni–Co–S/LIG as positive electrode and the commercial activate carbon (AC) as negative electrode, which were separated with the nonwovens separator (NKK-MPF 30AC-100, Japan) in 2 M KOH electrolyte.
To get the optimal electrochemical performance of asymmetric MSCs device, the charge balance of negative and positive should abide by the rule of following relationship q+ = q−. The charge storage in each electrode is related to the specific capacitance (C), the voltage window obtained from GCD curve (ΔV) and the mass of the active materials (m), showing in the equation below39
q = C × ΔV × m | (4) |
Hence, the mass of AC was calculated according to the above-mentioned formula.
The AC electrode was fabricated according to the procedure as followed. The negative materials were composed of AC, the binder (such as polyvinylidene fluoride, PVDF) and acetylene black, with the mass ratio of 8:1:1, respectively. Firstly, AC and acetylene black were mixed uniformly by grinding, and then a drop of the solvent (such as 1-methyl-2-pyrrolidone) was added while maintaining the mixing until the homogeneous slurry was formed. Secondly, the slurry was painted on a piece of nickel foam (the painting area is about 1 × 1 cm2), and the foam was pressed under 10 MPa for 1 min. Last, the nickel foam was dried at 80 °C for 12 h. The mass loading of the AC on nickel foam was 7.5–8.5 mg cm−2. Compared to the AC, the contribution of the nickel foam is ignorable under the same measurement conditions.
The specific areal energy density and specific areal power density of the asymmetric MSCs device were obtained according to two equations shown below24
(5) |
(6) |
The crystal structure and compositions of the as-prepared samples are characterized by X-ray diffraction (XRD), as shown in Fig. 2a. All diffraction peaks of Ni3S2/LIG are indexed to (110), (003), (202) and (113) plane reflections of rhombohedral phase Ni3S2 (JCPDS 044-1418). Except the peaks at 25.89° and 44.30° which are from the C·C substrate, other diffraction peaks of Ni–Co–S/LIG are indexed to cubic Co9S8 (JCPDS 003-0631) and rhombohedral phase Ni3S2 (JCPDS 044-1418). Obviously, the XRD patterns of Ni–Co–S/LIG differed from the XRD patterns of the carbon cloth and Ni–Co precursor as shown in Fig. S1.† The element compositions can also be confirmed by EDS analysis. As shown in Fig. 2c, Ni–Co–S/LIG consisted of five elements including C, Ni, Co, S and O. No diffraction peaks of oxide (MO, M = Ni or Co) can be observed in the XRD pattern of Ni–Co–S/LIG, which may indicate their low content or poor crystal structure. Ni3S2/LIG and Co9S8/LIG have been prepared by the same method under the same preparation condition with Ni–Co–S/LIG. As shown in Fig. S2a,† the XRD patterns of Ni–Co–S/LIG were the sum of that of Ni3S2/LIG and Co9S8/LIG. Fig. S2b† showed the XRD patterns with the different amount of Co mingling in the Ni–Co–S/LIG, indicating these four samples had the same composition.
Fig. 2 The XRD pattern of Ni–Co–S/LIG (a), Raman spectroscopy of C·C, C·C-starch and Ni–Co–S/LIG (b) and EDX spectrum of Ni–Co–S/LIG (c). |
The samples as-prepared were directly written using 10.6 μm CO2 laser scribing. The existence of graphene can be confirmed by the Raman spectrum measurement. In the grapheme structure, the D, G and 2D bands associated with the edge defect, the highly ordered graphite and characteristic of few layer graphene, respectively.41 As shown in Fig. 2b, three characteristic peaks of the D, G and 2D bands of LIG, starch-LIG and Ni–Co–S/LIG are appeared at 1350, 1580 and 2700 cm−1 approximately (Table S1†), which indicates the presence of graphene. Obviously, the 2D band of Ni–Co–S/LIG is more intense than that of the other two. The ratio of IG/ID expresses the amount of defects on the LIG and the degree of graphitization.42 As shown in Table S1,† the IG/ID ratio of Ni–Co–S/LIG (1.34) is higher than that of LIG (1.05) and starch-LIG (0.99), which indicates the existence of Ni–Co–S enhancing the degree of graphitization.43,44
To obtain the appropriate laser power, LIG was prepared under different power ranging from 2.4 to 7.2 W at the laser scan rate of 9 mm s−1. Raman spectroscopy is a powerful tool to obtain crystalline size (La) along “a” axis of graphitic materials by analyzing ratios of the integrated intensities of G and D peaks (IG/ID).45 Fig. S3a† shows typical Raman spectra of LIG prepared using laser powers from 2.4 W to 7.2 W. The relationship of the laser power and the IG/ID ratio is plotted in the Fig. S3c.† The La values are obtained using eqn (1), which is shown in Fig. S3d.† The IG/ID ratio and La value reach the highest when the laser power is 6.0 W, which are 1.75 and 34 nm, respectively. The higher temperature is beneficial to the surface temperature of the products46. However, too high temperature results in generating more D peaks that is not desired.24 Therefore, 6.0 W is chosen as the proper laser power based on its higher IG/ID ratio. Fig. S3b† shows the Raman spectra of Ni–Co–S/LIG with different Co content, and their peaks of D, G, 2D bands and the ratio of IG/ID are basically identical (Table S2†).
Typical morphologies of the samples are obtained by measuring FESEM and TEM. The mixture dipping in the carbon cloth contained Ni2+ and Co2+ that the concentration of Ni2+ fixed at 10 mmol L−1, while the percent content of Co2+ was changed from 2% to 15%. As shown in Fig. S4,† the degree of tightness of Ni–Co–S/LIG increased with smaller interstitials between the Ni–Co–S/LIG particles as the percent content of Co2+ increased from 2% to 15%. The typical SEM images of other samples are shown in Fig. 3. There are two differences between LIG and starch-LIG, as shown in Fig. 3a–d. Firstly, the surface of starch-LIG is the smoother than that of LIG. Secondly, the carbon clothes without starch have cracked more seriously than starch-LIG, and the diameter becomes smaller. A possible reason resulting in these two differences is that the carbon cloth was protected by the starch covering above the surface. This is consistent with the results of the Raman spectra shown in Fig. 2b. Furthermore, the surface of LIG and starch-LIG have marked big difference compared to the carbon cloth without laser scribbling (Fig. S5†).
Fig. 3 The SEM images of LIG (a and b), starch-LIG (c and d), Ni–S/LIG (e and f), Ni–Co–S/LIG (g and h). |
Metal ions mingling in the starch solutions form M-S/LIG (M = Ni or Ni–Co). Obviously, M-S/LIG was quite different than LIG or starch-LIG (Fig. 3). Plenty of particles appear on the surface of the carbon cloth, shown in Fig. 3f and h, which were confirmed as the mixture of Ni3S2/LIG and Co9S8/LIG through XRD characterization (Fig. 2a). The appearing of Co9S8 particles in Ni–Co–S/LIG leads to the surface of Ni–Co–S/LIG denser than that of Ni–S/LIG, since the surface of Co–S/LIG is tighter under the same concentration than Ni content in Ni–S/LIG (Fig. S6†).
Fig. 4 shows the typical TEM, HR-TEM and SAED pattern of Ni–Co–S/LIG. The different magnification of TEM images of Ni–Co–S/LIG are shown in Fig. 4a and b. The dark gray circular particles in the Fig. 4b are confirmed as the Ni–Co–S particles according to Fig. S7,† then the other light colour materials are the LIG. Majority of the of Ni–Co–S particles has a diameter of around 10 nm, as shown in Fig. 4b and e few other particles have the diameter of 25 nm approximately, indicating the Ni–Co–S particles were uniformly prepared. Fig. 4c is representative selected area electron diffraction (SAED) pattern of Ni–Co–S/LIG. The existence of concentric rings indicates the polycrystalline structure, which the rings can be indexed to the (110) and (300) diffractions of rhombohedral phase Ni3S2. The high-resolution TEM (HRTEM) images of the Ni–Co–S/LIG are shown in Fig. 4d and e. From Fig. 4d, the lattice fringe has 0.33 nm interplanar space, which is ascribed to the (220) plane of LIG corresponding to the light colour region LIG in Fig. 4b and the green colour of Fig. S7b.† As shown in Fig. 4e, the lattice fringe of 0.191 nm and 0.287 nm can be indexed to (511) diffraction of cubic Co9S8 and (110) diffraction of rhombohedral phase Ni3S2, respectively, corresponding to the dark gray region Ni–Co–S particles in Fig. 4b. The less intense XRD signal (Fig. 2a) and diffuse concentric rings (Fig. 4c) may be caused by the small particle size and poor polycrystalline structure. The TEM image of Ni–Co–S/LIG and the corresponding elemental mapping are shown in Fig. S8,† with all elements are well distributed in Ni–Co–S/LIG, indicating the high loading of Ni–Co–S particles in LIG. This conclusion can be drawn as well by the EDS spectrum (Fig. S9†).
The detailed information for the surface valence state of elements of as-prepared products was obtained using the X-ray photoelectron spectroscopy (XPS) and the results are shown in Fig. 5 and Fig. S10.† The Ni 2p spectra of XPS is well fitted with two spin–orbit splitting, which has doublet characteristics of Ni2+ and Ni3+. As shown in Fig. 5a, the binding energies of Ni 2p peaks at 856.98 eV and 875.11 eV correspond to Ni2+ located in the octahedral sites.47 The binding energies at 855.55 eV and 873.30 eV can be assigned to Ni3+, which located in tetrahedral sites.48 The intensity of the double peaks is stronger than that of Ni2+, indicating Ni3+ is the dominant state of Ni near the surface of the sample and comes from Ni2S3, which is consistent with the XRD analysis. The Co 2p spectrum of XPS is similar to Ni spectra shown in Fig. 5b. The binding energies of stronger peaks at 781.28 eV and 796.38 eV can be assigned to Co3+,47 corresponding to Co 2p3/2 and 2p1/2, respectively. The other spin–orbit doublets locating at 783.60 eV for Co 2p3/2 and 797.57 eV for Co 2p1/2 can be ascribed to Co2+.36 There are two shake-up satellites (denoted as Sat.) existing in both Ni and Co 2p spectra of XPS. The satellite peaks of Co 2p are relatively weak, which is ascribed to the dominance of Co3+.49 Based on the conclusion above, the compositions of as-prepared Ni–Co–S/LIG include Ni2+, Ni3+, Co2+ and Co3+, which can contribute as positive materials with high electrochemical property.50 The S 2p spectra shown in Fig. 5c possesses two peaks round at 162.68 eV for S 2p3/2 and 163.78 eV for S 2p1/2, which correspond to the S in low coordination and metal–sulfur bonds.48 The O 2p scan is shown in Fig. 5d, and two peaks observing at 532.15 eV and 533.22 eV correspond to O 2p3/2 and 2p1/2, respectively.
CoS + OH− ⇌ CoSOH + e− |
CoSOH + OH− ⇌ CoSO + H2O + e− |
NiS + OH− ⇌ NiSOH + e− |
The CV curves of Ni–Co–S/LIG have higher current and larger integral areas than that of Ni–S/LIG, illustrating that Ni–Co–S/LIG have higher capacitance. The excellent electrochemical property of Ni–Co–S/LIG may be attributed the addition of the Co element. The CV curves of a series of Ni–Co–S/LIG with different Co content are shown in Fig. S12,† with the scan rates changing from 10 to 100 mV s−1. All of the CV curves show typical couple of peaks and the shapes of the CV curves are maintained well along the scan rate rising to 100 mV s−1. In addition, the peak current rises/enhances gradually along with the increase of Co content or the scan rate. Obviously, the anodic peaks in CV curves (Fig. S12 a–d†) shift to higher voltage with the increase of Co content at the same scan rate, which consistent with the electrochemical behaviour of Co–S and Ni–S reported previously.52,53. Differently, the CV curves of substrate LIG and starch-LIG did not retain original shape with obvious distortion along with the scan rate from 10–100 mV s−1 (Fig. S13†) for the polarization.
The GCD curves of the as-prepared samples at the current density of 1 mA cm−2 are shown in Fig. 6b. Obviously, the discharge time of Ni–Co–S/LIG is much longer than that of the other three materials, in accordance with the result of CV curves in Fig. 6a. To further evaluate the charge transfer resistance and the diffusive resistance of electrolyte, the electrochemical impedance spectroscopy (EIS) measurement of as-fabricated samples were conducted from 0.1 MHz to 100 MHz, and the results are as shown in Fig. 6c. All the plots show a semicircle and a line appearing in the high and low frequency range, respectively. The diffuse resistance of electrolyte can be calculated from the slope of the curves in the low frequency range, which is defined as the Warburg impedance (W). The Ni–Co–S/LIG has small slopes of the straight lines, indicating the high ion diffusion resistance, which may be attributed to the addition of second cobalt ions. The diameter of semi-circle is related to the charge transfer resistance (Rct). The large the diameter, the greater the Rct. Similar with Warburg impedance, the Rct value of Ni–Co–S/LIG is larger than the others owing to Co ions. The intercept with the X-axis at high frequency is related with solution resistance (Rs). From Fig. 6c, it can be seen all the samples have small Rs. The galvanostatic charge/discharge curves of the Ni–Co–S/LIG were carried out at diverse current densities with the working potential range of 0–0.4 V (Fig. 6d). All the discharge curves at different current densities exhibit the distinct plateau region, which is the characteristic of the typical pseudocapacitive electrode in line with the result of CV curves (Fig. 6a).
The GCD profiles of Ni–Co–S/LIG with different Co content are shown in Fig. S14.† All discharge curves have plateau region, indicating Ni–Co–S/LIG possesses the pseudocapacitive character regardless of the Co content.
The areal specific capacitance of LIG, starch-LIG, Ni–S/LIG and Ni–Co–S/LIG (Co-10%) were calculated to be 4.2, 1.85, 75 and 680 mF cm−2 at a low current density of 1 mA cm−2, respectively, which was higher than cobalt-based composite.54 And the capacitance retention rate of Ni–S/LIG and Ni–Co–S/LIG were 73% and 90%, respectively, indicating the increase of capacitance retention with the addition of Co ion. The areal specific capacitance of LIG, starch-LIG, Ni–S/LIG and Ni–Co–S/LIG were calculated from the discharge time and plotted in Fig. 6e at different current density. Ni–Co–S/LIG shows a exceeding large areal specific capacitance of 680 mF cm−2, which is 162 times larger than the substrate LIG (4.2 mF cm−2) and about 368 times larger than that of starch-LIG (1.85 mF cm−2) at the 1 mA cm−2. The areal specific capacitance of starch-LIG is lower than LIG, attributing to the protection of starch when the carbon cloth is written using a CO2 laser. Also, Ni–Co–S/LIG has higher capacity than Ni–S/LIG, which may be the synergistic effect of Co ion. The value of areal specific capacitance decreases gradually along with the increasing current density from 1 to 10 mA cm−2. It is noteworthy that the capacitance of Ni–Co–S/LIG can reach a capacitive retention of 88.3% at the large current density of 10 mA cm−2, which is much higher than that of the substrate LIG (4.2 mF cm−2 at 1 mA cm−2), starch-LIG (1.85 mF cm−2 at 1 mA cm−2) and Ni–S/LIG (75 mF cm−2 at 1 mA cm−2). What is more, Ni–Co–S/LIG has excellent performance compared to a plenty of previously reported LIG-based electrode materials. To the best of our knowledge, such an excellent rate performance has superior to most pseudocapacitor-type electrode materials including carbon materials.50,52,55–57 The cycling stability of the Ni–Co–S/LIG was further investigated under the GCD process running for 10000 cycles at a current density of 5 mA cm−2, as shown in Fig. 6f. Remarkably, the areal specific capacitance retention can reach ca. 92.1% of the initial value after 10000 cycles, illustrating excellent long-life span cycling stability. And the 92.1% capacitance retention rate is superior to that of many electrode materials reported previously.52,57,58
The energy storage of the as-fabricated samples used as electrode materials in the practical application is evaluated to assemble Ni–Co–S/LIG as the positive electrode and active carbon as the negative electrode into the asymmetric MSCs device, as schematically illustrated in Fig. 7a. The mass ratio of Ni–Co–S/LIG and AC is balanced according to the eqn (4) before assembling the asymmetric MSCs device. Take full advantage of the complementary voltage windows ranges of Ni–Co–S/LIG and AC, as shown in Fig. 7b, the working potential of the asymmetric MSCs device can be enlarged to 1.6 V. To search the optimum operating potential window of asymmetric MSCs device, a series of CV measurements at the different potential windows in 2 M KOH electrolyte are conducted, shown in Fig. 7c. Finally 1.6 V is chosen as the maximum voltage window. When the voltage window exceeds 1.6 V, an obvious polarization is observed, which is attributed to the oxygen evolution reaction referring to the secondary reaction.59,60 As shown in Fig. 7d, all CV curves of the asymmetric MSCs device at various scan rates from 10 to 100 mV s−1 at 0–1.6 V are observed a couple of redox peaks appearing, which is ascribed to the reversible redox reaction accompanied with the electrolyte ions insertion/extraction. The shape of CV curve is almost maintained at the high scan rate of 100 mV s−1, which illustrates the electrode materials possess excellent rate property as a consequence of fast electronic and ions transportation in the electrode materials. The typical GCD profiles at various current densities ranging from 1 to 10 mA cm−2 are shown in Fig. 7e. The plateaus of all the GCD curves are the characteristic of pseudocapacitors, which is in accordance with the CV analysis. The areal specific capacitance is calculated from typical GCD curves at a series of current densities, which is plotted in Fig. 7f. The areal specific capacitance of Ni–Co–S/LIG (Co-10%)//AC asymmetric MSCs device is 160 mF cm−2 at the current density of 1 mA cm−2, which reaches 141 mF cm−2 with the 90.1% of capacitance retention at the current density of 10 mA cm−2, higher than that of Ni–Co–S/LIG (Co-0%)//AC asymmetric MSCs device (70.6%).
As shown in Fig. S15,† the minor change of IR drop at different charge/discharge current densities and the gentle slope can be attributed to small internal resistance and excellent electric and ionic conductivities of the bimetal transition metal sulfides deposited in Ni–Co–S/LIG. Fig. S16† shows that the two asymmetric MSCs devices have small Rs. It is a little difference that the Rct value of Ni–Co–S/LIG (Co-0%)//AC asymmetric MSCs is larger than that of Ni–Co–S/LIG (Co-10%)//AC asymmetric MSCs.
The cycling performance of the Ni–Co–S/LIG//AC asymmetric MSCs was investigated at a high current density of 5 mA cm−2, and the result was shown in Fig. 7g. It is observed that the long-term cycling stability of the Ni–Co–S/LIG//AC asymmetric MSCs device retains 89.6% of its initial areal specific capacitance after 8000 cycles, which indicates the as-fabricated asymmetric MSCs device possesses desirable cycling stability. Other than some collapse found after 8000 cycles, no obvious change was observed between the morphology of the Ni–Co–S/LIG before and after charging shown in Fig. S17.†
Fig. 7h shows a Ragone plot, illustrating the relationship between the energy density and the power density of Ni–Co–S/LIG//AC asymmetric MSCs device. It can be clearly seen that the areal energy density of the asymmetric MSCs device is 56.9 μW h cm−2 at the areal power density of 800 μW cm−2. Moreover, an areal energy density of 50.13μWh cm−2 can be obtained at an areal power density of 8000 μW cm−2. These results exceed a plenty of previously-reported supercapacitors in literatures (Table S3†), such as LIG-MnO2-MSCs (32.4μWh cm−2 at 2334 μW cm−2),61 LIG-PANI-MSCs (8.0μWh cm−2 at 649 μW cm−2),61 LIG-FeOOH//LIG-MnO2 (9.6 μW h cm−2 at 11853 μW cm−2),61 MnO2/CNT/nylon FSCs (2.6 μW h cm−2 at 66.9 μW cm−2)62, LIG (0.8 μW h cm−2 at 2290 μW cm−2),61 LIG-MoS2 (2.8 μW h cm−2 at 440 μW cm−2)25 just name a few. Hence, the Ni–Co–S/LIG//AC asymmetric MSCs device provides an effective strategy for building high-energy and high-power density asymmetric supercapacitors. To further investigate the practical application, one or two asymmetric MSCs devices can light 17 red and 8 yellow commercial (2.0 V) LEDs, as shown in Fig. S18,† indicating the as-prepared asymmetric MSCs device has high power and energy performance. Obviously, the brightness of three asymmetric MSCs devices lighting the 25 commercial LEDs is higher than that of two asymmetric MSCs devices, which illustrates three asymmetric MSCs devices can store more energy.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ra06376b |
This journal is © The Royal Society of Chemistry 2020 |