3.1 Control variable method
3.1.1 Effects of the molar ratio of KMnO4 and MnSO4·H2O
The effects of the molar ratio of KMnO4 and MnSO4·H2O in the growth solution on the morphology, structure and electrochemical properties of MnO2@Ppy@NPCNFs were investigated. In the hydrothermal synthesis of MnO2, the reaction temperature was maintained at 160 ℃ and the reaction time was kept at 6 h. The molar ratios of KMnO4 and MnSO4·H2O were set to 1:3, 2:3, 3:3, 5:3 and 9:3, respectively.
The crystal types of MnO2 prepared using different molar ratios of KMnO4 and MnSO4·H2O were characterized by X-ray diffraction, and the XRD spectra of MnO2@Ppy@NPCNFs were shown in Fig. 2a. It could be seen that when the molar ratios were 1:3 and 2:3, the products had obvious diffraction peaks, indicating that they had good crystallinity and almost no amorphous structure. Compared with the standard card JCPDS NO.24-0735, the products belonged to the tetragonal system of β-MnO2 single-chain structure [31-32]. When the molor ratio was 3:3, there were characteristic diffraction peaks at 2θ of 12.2°, 17.9°, 28.7°, 37.3°, 41.8°, 49.6°, 59.7 ° and 69.5° in the products, which were consistent with the diffraction peaks of the standard card JCPDS NO.44-0141. Therefore, the products generated at 3:3 belonged to the tetragonal α -MnO2 system, and their structure were I4/m (87) spatial point group . When the molar ratio increased to 5:3, the products had characteristic peaks at 12.2°, 24.9°, 36.9° and 66.4°. Compared with the standard card JCPDS NO.80-1098, these four diffraction peaks corresponded to (001), (002), (201) and (021) crystal planes of δ-MnO2, respectively, illustrating the formation of δ-MnO2. The δ-MnO2 was a rhomboid crystal system with a layered structure, which belonged to potassium-containing hydrated MnO2 [33]. When the molar ratio reached 9:3, the products were still δ-MnO2, but the crystallinity of the products was lower than that of the products at 5:3.
SEM images of MnO2@Ppy@NPCNFs obtained from the growth solution with different molar ratios of KMnO4 and MnSO4·H2O were shown in Fig. 2b-f. According to XRD patterns of MnO2@Ppy@NPCNFs (Fig. 2a), Fig. 2b and c indicated that the products at 1:3 and 2:3 molar ratios were β-MnO2 rods, and they are similar in morphology and size, with lengths of 16-30 μm and widths of 2-3 μm. However, the surface of the products at 2:3 was smoother than that at 1:3. When the molar ratio increased to 3:3, the products were α-MnO2, which presented a needle-like structure with lengths of 200-900 nm and diameters of 25-40 nm, as displayed in Fig. 2d. Furthermore, as shown in Fig. 2e, the products at 5:3 belonging to δ-MnO2 appeared many flower-like nanospheres composed of a large number of nano-flakes interleaved and stacked with each other. The thickness of nano-flakes was about 10 nm, the lamellar spacing was about 4 nm, and the diameter of flower-like nanospheres was 170-240 nm. When the molar ratio reached 9:3, the products were still flower-like nanospheres, but compared with the products at 5:3, the folds and lamellar structure on the surface of nanospheres were reduced, the lamellar spacing increased to about 20 nm, and the diameter of flower-like nanospheres increased to 200-400 nm, as illustrated in Fig. 2f.
Accordingly, the above XRD and SEM characterization indicated that the molar ratio of KMnO4 and MnSO4·H2O in the growth solution had a great influence on the crystal type of MnO2 generated on Ppy@NPCNFs. When the content of oxidizer KMnO4 was higher than MnSO4·H2O, δ-MnO2 with the lamellar structure was produced. With the increased of MnSO4·H2O addition, the crystal type of MnO2 generated gradually changed from δ-MnO2 to α-MnO2, illustrating that the coexistence of δ-MnO2 required sufficient K+ to maintain the lamellae structure. With the decrease of K+, the lamellar structure of δ -MnO2 began to collapse and form α-MnO2 with [2X2] tunnel structure. Then, after the further decrease of K+, α-MnO2 transformed into β-MnO2 with [1X1] tunnel structure with higher crystallinity. Therefore, K+ could be regarded as a template for the lamellar structure and tunnel structure of MnO2 formed in the reaction process [34], and the final products were controlled by changing the K+ concentration in the reaction process [35, 36].
The electrochemical performance of MnO2@Ppy@NPCNFs prepared using different molar ratios of KMnO4 and MnSO4·H2O were displayed in Fig. 3. Fig. 3a showed the CV curves of the above products at the scanning speed of 10 mV s-1, which were approximately symmetrical rectangles, suggesting their good electrochemical reversibility and electrical double-layer capacitance (EDLC) characteristics [37-39]. Furthermore, these CV curves exhibited a broad weak redox peak, indicating that they had the characteristics of Faraday pseudocapacitance [40]. Moreover, the enclosed area of CV at 5:3 was the largest, suggesting that it had the highest SC. Fig. 3b displayed the GCD curves of each product at a constant current of 1 A g-1, which had good symmetry and presented a good linear relationship, further indicating their excellent cycling performance and electrochemical reversibility [41, 42]. According to the GCD curves of MnO2@Ppy@NPCNFs obtained from different molar ratios of KMnO4 and MnSO4·H2O at 1 A g-1 , their SC calculated was 562.12 F g-1 (1:3), 569.33 F g-1 (2:3), 574.35 F g-1 (3:3), 588.76 F g-1 (5:3) and 580.98 F g-1 (9:3), respectively.
In conclusion, the electrochemical performance of the obtained MnO2@Ppy@NPCNFs at 5:3 was the best. This was because at this time, the MnO2 loaded on Ppy@NPCNFs was flower-like nanosphere composed of δ-MnO2 with the best lamellar structure, and its surface had more folds, which was beneficial to increase the specific surface area, the contact area between electrode material and electrolyte, and the energy storage site of electrochemical reaction, thus improving its electrochemical performance. In addition, the interlayer spacing of the product δ-MnO2 at 5:3 was about 4 nm (Fig.2e), while the hydration radius of K+ was 0.33 nm, which was much smaller than the interlayer spacing of the product δ-MnO2. Therefore, the free interpenetration of K+ between the layers of δ-MnO2 could be carried out [43]. The special interlayer structure of δ-MnO2 can also form EDLC inside the δ-MnO2, which greatly increased the EDLC of MnO2@Ppy@NPCNFs, thus enhancing the overall electrochemical performance of MnO2@Ppy@NPCNFs. In contrast, α-MnO2 and β-MnO2 with tunnel structure were not conducive to the penetration of K+ in the electrolyte, resulting in relatively smaller capacitance [44].
3.1.2 Effects of the reaction temperature
According to the optimal molar ratio of KMnO4 and MnSO4·H2O determined by the above experiments, the effects of the reaction temperature on the MnO2 crystal form, morphology and electrochemical properties of MnO2@Ppy@NPCNFs were discussed. In this section, the molar ratio of KMnO4 and MnSO4·H2O in the growth solution was 5:3, and the reaction time was kept at 6 h. The reaction temperatures were set at 120 ℃, 140 ℃, 160 ℃, 180 ℃ and 200 ℃, respectively.
Fig.4a showed the XRD patterns of MnO2@Ppy@NPCNFs prepared at different reaction temperatures. It could be found when the reaction temperature was between 120-160 ℃, there were characteristic peaks at 12.2°, 24.9°, 36.9° and 66.4°, which corresponded to (001), (002), (201) and (021) crystal planes of δ-MnO2, respectively. But when the temperature was 120 ℃, the intensity of the diffraction peaks was weak, which was because the low reaction temperature led to the low nucleation activity of MnO2, resulting in incomplete crystal nucleus growth of MnO2. As the temperature increased, the characteristic diffraction peak intensity of δ-MnO2 increased gradually, indicating that its crystallinity increased. When the temperature rose to 180-200 ℃, the characteristic peaks of products illustrated that they were β-MnO2 with higher crystallinity, in which the diffraction peak intensity of the product fabricated at 200 ℃was stronger, indicating that β-MnO2 with a good crystal structure was formed.
Fig. 4b-f showed the SEM images of MnO2@Ppy@NPCNFs obtained at different reaction temperatures. It was found that the products prepared at 120-160 ℃ (Fig.4b-d) all consisted of a large number of flower-like nanospheres, which were δ-MnO2 according to the XRD results. And with the increase of temperature the diameter and lamellar spacing of nanospheres decreased. As shown in Fig.4b-d, the diameter of nanospheres was reduced from 280-450 nm (120 ℃) to 190-270 nm (140 ℃) and then to 170-240 nm (160 ℃) , and their corresponding lamellar spacings decreased from about 100 nm to about 20 nm and then to about 4 nm. When the temperature reached to 180 ℃ (Fig.4e), the products were transformed from δ-MnO2 to β-MnO2 by XRD analysis, the morphology of the products changed from nanospheres to nanorods with lengths of 0.7-1.3 μm and widths of about 300 nm, which had unstable structure and uneven surface. This might be due to the collapse of nanosphere lamellar structure caused by the higher reaction temperature. With the temperature further increased to 200 ℃ (Fig.4f), the products changed from unstable nanorods with rough surface to quite complete and stable microrods with smooth surface, and their lengths and widths were 15-30 μm and 2-3 μm, respectively. It might be because the shorter nanorods were dissolved due to the high reaction temperature, which made the micron structure uniform and stable.
Electrochemical performance
Fig. 5a showed the CV curves of products prepared at different reaction temperatures in a three-electrode system at a scanning rate of 10 mV s-1, which were similar to the CV curves in the above section, indicating that they had obvious EDLC characteristics and good electrochemical reversibility. And the MnO2@Ppy@NPCNFs obtained at 160 ℃ had the largest enclosed area of CV curve, suggesting its largest SC. Fig. 5b exhibited the GCD curves of these products at 1 A g-1, which had good linear symmetry in the potential range of -1.0-0 V (vs.SCE), further indicating their good coulombic efficiency and electrochemical reversibility. Based on these GCD curves, the SC of the products was 556.43 F g-1 (120 ℃), 578.90 F g-1 (140 ℃), 588.76 F g-1 (160 ℃), 539.12 F g-1 (180 ℃) and 511.32 F g-1 (200 ℃), respectively.
In conclusion, when the reaction temperature was 160 ℃, the electrochemical performance of MnO2@Ppy@NPCNFs was the best. This was because based on the XRD and SEM results, the flower-like δ-MnO2 nanospheres loaded on Ppy@NPCNFs at 160 ℃ had the smallest size, resulting in the smallest lamellar spacing and the largest number of nano-flakes, which provided the largest specific surface area and the largest number of active sites to promote their electrochemical performance.
3.1.3 Effects of the reaction time
Base on the previous experiment results, the effects of the reaction time on the MnO2 crystal form, morphology and electrochemical performances of MnO2@Ppy@NPCNFs were explored. Therefore, the molar ratio of KMnO4 and MnSO4·H2O was fixed at 5:3, the reaction temperature remained 160 ℃, and the reaction time was 2 h, 4 h, 6 h, 8 h and 10 h, respectively.
Combined with XRD and SEM analysis, as shown in Fig. 6, when the reaction time was 2-6 h (Fig. 6b-d), the products were all flower-like δ-MnO2 nanospheres , and the diameter and lamellar spacing of nanospheres decreased as the reaction time increased. When the reaction time was 8-10 h (Fig. 6e and f), the MnO2 crystal form of the products changed from flower-like δ-MnO2 nanospheres to needle-like α-MnO2 nanowires, and the diameter and length of nanowires decreased with the increase of reaction time.
The above experimental results were due to the spontaneous reaction of KMnO4 and MnSO4·H2O at a certain temperature and pressure during the whole reaction process, and the MnOx structural unit was first formed [45]. When the reaction time was long enough, nanospheres with large lamellar structure were gradually formed due to the concentration reaction, and their lamellar spacing was large. As the reaction time increased, the nanosheets gradually gathered and the lamellar spacing decreased, resulting in the formation of flower-like nanospheres. With the further increase of reaction time, the flower-like nanosphere size continued to decrease. Subsequently, in order to maintain the stable structure of the products and reduce the surface free energy, the lamellar structure of the flower-like nanospheres gradually collapsed and shrank, and then the nanowire structure was formed.
Electrochemical performance
Fig. 7 showed the CV (a) and GCD (b) curves of the products obtained after different reaction times. It could be found these products had good EDLC characteristics, electrochemical reversibility and coulombic efficiency, and MnO2@Ppy@NPCNFs prepared after 6 h had the largest enclosed area of CV curve, indicating its largest SC. In addition, according to their GCD curves, the SC of MnO2@Ppy@NPCNFs was 574.23 F g-1 (2 h), 581.98 F g-1 (4 h), 588.76 F g-1 (6 h), 558.32 F g-1 (8 h) and 543.35 F g-1 (10 h), respectively. Therefore, the electrochemical performance of MnO2@Ppy@NPCNFs obtained after 6 h was the best because of their flower-like δ-MnO2 nanospheres with the smallest size. In conclusion, when the molar ratio of KMnO4 and MnSO4·H2O was 5:3, the reaction temperature was 160 ℃ and the reaction time was 6 h, the electrochemical performance of MnO2@Ppy@NPCNFs obtained was the best and its SC was the largest .
3.2 Response surface methodology and experimental verification
Based on the above results of control variable method, the actual values of factors and their corresponding response values were shown in Table S2. The synthesis process of MnO2 was designed by Design-Expert software to study the effects of reaction temperature, molar ratio and reaction time on the electrochemical performance of MnO2@Ppy@NPCNFs, and a statistical model was established according to Eq. (3). Subsequently, the established model was applied to predict the optimal reaction parameters, and the experimental verification was carried out, thus obtaining the MnO2@Ppy@NPCNFs with the best electrochemical performance.
3.2.1 Statistical analysis
Analysis of variance was performed on the experimental data in Table S2, and the analysis results were shown in Table S3. According to the analysis results and Eq. (3), a quadratic regression model of the SC (Y) against the reaction temperature (A), the molar ratio of reactants (B) and the reaction time (C) could be written as follows:
Y = 563.83-17.40A+ 15.33B-6.58C+4.92AB-4.87AC-2.19BC-14.23A2-24.44B2-17.24C2 (4)
It could be seen from Table S3 that the F-value of this model was 82.82 and its P-value<0.0001, indicating an extremely significant level. P-value of the lack of fit greater than 0.05 showed the difference of the term was not significant, indicating that the model had a high fitting degree with the actual experimental values and was highly reliable. Moreover, the R2 value of 0.9907 and the Adj-R2 value of 0.9787 further exhibited this conclusion. In addition, it was found from Fig. S1 that there was a linear correlation between the predicted values obtained by the model and the experiment data, and the fitting degree between them was high.
3.2.2 Response surface analysis
The regression function can usually be represented by contour and 3D response surface plots, which is mainly used to characterize the interaction type between the two variables and the relationship between the response of each variable and experimental level. The elliptical contour shows that the interaction between variables is obvious, while the circular contour is the opposite [46,47]. In this experiment, the influences of three reaction parameters on the SC of MnO2@Ppy@NPCNFs were shown in Fig. S2. It could be found that all 3D response surface plots indicated that with the increase of each single reaction parameter, the SC of MnO2@Ppy@NPCNFs increased first and then decreased. In addition, the elliptical contour plots in Fig. S2 (a and b) illustrated that the interaction between the reaction temperature and molar ratio was significant when the reaction time remained constant, and there was a significant interaction between the reaction time and reaction temperature when the molar ratio was kept constant. Fig. S2c showed when the reactant temperature was kept constant, the interaction between the molar ratio and reaction time was not significant due to the circular contour plot. The results of response surface analysis were consistent with the regression equation results in Table S2, and there was an extreme value of SC in Fig. S2, indicated that the parameter range of the experiment was reasonable and the designed model had strong reliability.
3.2.3 Prediction of optimal reaction parameters
According to the established regression model (Eq. (4)), the optimal reaction parameters were obtained: the molar ratio of KMnO4 to MnSO·H2O was 1.93, the reaction temperature was 140.64 ℃, and the reaction time was 6.26 h. It was predicted that the SC of MnO2@Ppy@NPCNFs obtained under these optimal reaction parameters was 593.809 F g-1.
3.2.4 Experimental verification
Three groups of repeated experiments were carried out to verify the optimal reaction parameters predicted by the model. The GCD curves at 1 A g-1 of MnO2@Ppy@NPCNFs (1.93/140.64℃/6.26h) obtained by three verification experiments were shown in Fig. 8a. Based on these GCD curves, the SCs of the three experimental products calculated were 597.21 F g-1, 590.45 F g-1, and 596.65 F g-1, respectively, and the average SC was 595.77 F g-1. The error between the average experimental value and the theoretical value was very small, which indicated that the optimization results were reliable, and the optimal reaction parameters obtained had practical application value.
The MnO2@Ppy@NPCNFs (1.93/140.6℃/6.26h) was characterized by XRD. Compared with the sample with the largest SC obtained in the previous experiments (the molar ratio of 5:3 (1.67), the reaction temperature of 160 ℃, the reaction time of 6 h), it had a similar diffraction pattern. As shown in Fig. 8b, their characteristic diffraction peaks corresponded to (001), (002), (201) and (021) crystal planes, indicating that the outmost layer of MnO2@Ppy@NPCNFs (1.93/140.64℃/6.26h) was also δ-MnO2. According to its SEM images, the MnO2@Ppy@NPCNFs (1.93/140.64/6.26h) showed an obvious characteristic morphology of δ-MnO2, which were flower-like nanospheres with diameters of 70-180 nm and lamellar spacing of about 2.5 nm (Fig. 9a. And the product had an obvious multicore-shell ternary structure, as indicated in Fig. 9b, in which NPCNFs with the fiber diameter of about 280 nm served as the core layer, Ppy was used as the intermediate layer with a thickness of about 50 nm, and the loaded δ-MnO2 was the outermost shell layer with a thickness of about 100 nm. The innermost NPCNFs provided porous structure and conductive skeleton, and the intermediate layer Ppy enhanced the conductivity of the material. The outermost MnO2 not only improved the disadvantages of low life and unstable performance of the middle layer Ppy, but also provided more active sites. Therefore, the final product MnO2@Ppy@NPCNFs would have excellent electrochemical performances under the synergistic action of the three internal materials. As shown in (Fig. 9c), the lattice fringes with the spacing of 0.35 nm and 0.71 nm could be found in the HR-TEM image of MnO2@Ppy@NPCNFs, which were indexed as the (002) and (001) planes of δ-MnO2, respectively. The elemental distribution of MnO2@Ppy@NPCNFs was further verified by the elemental mapping image (Fig. 9d, where O was distributed in the whole fiber, C and N were evenly distributed in the inner layer of the fiber, while Mn was only distributed in the outermost layer of the fiber, indicating the multilayer structure of the fiber and the outermost layer was MnO2.
Fig. 10 showed the comparison of electrochemical properties of the three electrode materials, including MnO2@Ppy@NPCNFs (1.93/140.64℃/6.26h), MnO2@Ppy@NPCNFs (1.67/160.00℃/6.00h) and Ppy@NPCNFs). According to the CV curves (Fig. 10a) and GCD curves (Fig. 10b), it could be concluded that the three products all had good EDLC characteristics and electrochemical reversibility, and their CV curves all displayed a broad and weak redox peak at 10 mV s-1, especially MnO2@Ppy@NPCNFs, indicating that they had the characteristics of Faraday pseudocapacitance. In addition, according to the GCD curves, the SC of MnO2@Ppy@NPCNFs (1.93/140.64℃/6.26h) (596.65 F g-1) was 74.2% higher than that of Ppy@NPCNFs (342.13 F g-1). It was mainly due to the excellent capacitive characteristics of a large number of flower-like nanospherical δ-MnO2 loaded on the Ppy@NPCNFs surface, which mdae the electrode material have a large number of active sites and larger specific surface area, thus promoting the adsorption and desorption of ions in the electrolyte on the material surface.
Meanwhile, the electrochemical impedance spectra of the three products (Fig. 10c) showed that the semicircle diameter of MnO2@Ppy@NPCNFs (1.93/140.64℃/6.26h) in the high-frequency region was the smallest, and its slope of the inclined line in the low-frequency region was the largest, illustrating that it has the smallest ion transfer resistance and charge diffusion resistance, which made the electrode material had better capacitance performance [48,49]. The SC attenuation and coulomb efficiency variation of MnO2@Ppy@NPCNFs (1.93/140.64℃/6.26h) after 1000 cycles of GCD at 1 A g-1 (Fig. 10d) displayed that its SC retention rate reached 96.2% and its coulomb efficiency was stable at 91±3%, which were higher than those of Ppy@NPCNFs (95.4% and 86±2%) [18]. The main reason was that the metal oxide MnO2 made up for the structural instability of conductive polymer Ppy in a long time cycle. All of the above indicated that MnO2@Ppy@NPCNFs had excellent electrochemical performances, and its electrochemical performances were superior to most of the MnO2-loaded carbon materials previously reported (Table S4). According, it was an ideal choice for supercapacitor electrode materials.
The capacitance contribution of materials mainly includes two parts: EDLC (Cd) and pseudocapacitance (Cp), in which Cd is attributed to the contribution of specific surface area and inherent structural defects with micropores, while Cp is derived from the contribution of active sites. The total capacitance contribution of materials will be the sum of Cd and Cp. Trassatti method can be performed to investigate the capacitive storage behavior of electrode materials , and the total voltammetric charge (q) measured can be expressed as follows [50]:
where k1 and k2 are constants, v is a scan rate, q∞ is the storage charge at v→∞, which is equivalent to Cd, and q0 is the storage charge at v = 0.
The variation curve of specific capacitance for MnO2@Ppy@NPCNFs was determined through GCD curves at different current densities (Fig. S3), and its specific capacitance of MnO2@Ppy@NPCNFs (595.53 F g− 1 at 1 A g-1) was much higher than other samples [18], as shown in Fig. 11d. Then, according to the CV curves of the four samples at different scan rates (Fig. S4), their specific capacitances were determined by formulas (5) and (6), respectively, as shown in Fig. 11a-c. It could be seen that the MnO2@Ppy@NPCNFs electrode displayed higher Cd and Cp values than other samples, especially Cp value, illustrating that the large number of active sites provided by MnO2 on its outmost layer could significantly enhance the charge storage efficiency through the pseudocapacitance mechanism.
Consequently, a symmetric supercapacitor was assembled using MnO2@Ppy@NPCNFs as electrodes for evaluating its practicability. The CV and GCD curves of the symmetric supercapacitor measured in a two-electrode configuration were exhibited in Fig. 12a and b, respectively. Its energy density and power density (Fig. 12c) were calculated from the GCD curve according to Eqs. (1 and 2), and illustrated an energy density of 9.36 Wh kg−1 at a power density of 1000 W kg−1, which were significantly higher than those of the Ppy@NPCNFs-based symmetric supercapacitor (3.10 Wh kg−1 at 50 W kg−1) [18]. In addition, the retention rate of specific capacitance could reach 93.75% after 1000 cycles at 10 A g−1 (Fig.12d). The electrochemical performance comparison of MnO2@Ppy@NPCNFs-based symmetric supercapacitor with the MnO2-loaded supercapacitors previously reported was shown in Table S5.