3.1 Characterization of Z-V5S8/CoS-CNFs
Figure 1 illustrates the synthesis process of free-standing Z-V5S8/CoS-CNFs anode materials. The sample was prepared through a pre-oxidation treatment (Fig. S1a) followed by a subsequent one-step sulfidation process (Fig. S1b). The morphology and crystal structure of Z-V5S8/CoS-CNFs were characterized using SEM and TEM. As depicted in Fig. 2a-2c, SEM analysis revealed the smooth surface of the CNFs. However, the surface exhibited roughness after pre-oxidation, likely attributed to the expulsion of small gas molecules from the material during the pre-oxidation process[24].
The structure of Z-V5S8/CoS-CNFs presents a one-dimensional homogeneous nanofiber structure with an average diameter of 200 nm (Fig. S2). This small size is pivotal in shortening the Na+ transport path during charge and discharge processes[25]. As the metal concentration increases, there is a slight trend of increasing fiber diameter. Additionally, there is a noticeable increase in the number of metal particles adhering to the fiber surface. Beyond a metal concentration of 0.2 mol L− 1, spherical and lamellar crystals become distinctly visible on the fiber surface (Fig. S3). These correspond to cobalt sulfides and vanadium sulfides, respectively, indicating that the carbon nanofiber fails to effectively bind the metal sulfide crystals that form during the sulfidation process.
Subsequently, the 0.2-V5S8/CoS-CNFs were further characterized by TEM. In Fig. 2d, it can be observed that V5S8 and CoS are uniformly dispersed within the CNFs, with the darker granular shapes representing the nanoparticles of the two metal sulfides, aligning with the original design. In Fig. 2e, the two distinctive lattice stripes of the substances correspond to the (420) crystal plane of V5S8 and the (002) crystal plane of CoS, with interplanar distances of 0.227 nm and 0.276 nm, respectively. The clearly defined interface between the two crystalline surfaces confirms the presence of a heterostructure. The crystal plane spacing periodogram (Fig. 2e) further supports the accurate determination of V5S8 and CoS lattice spacing. Interfaces and lattice mismatch (m) can be calculated using Eq. (1) [26]:
$$\text{m=}\frac{\left|{\text{d}}_{\text{1}}\text{sin}\text{θ-}{\text{d}}_{\text{2}}\right|}{\text{0.5×(}{\text{d}}_{\text{1}}\text{sin}\text{θ}\text{+}{\text{d}}_{\text{2}}\text{)}} \left(1\right)$$
Where m represents the lattice mismatch, d₁ and d₂ are the interplanar spacings associated with specific crystallographic planes, and θ denotes the angle of incidence, also known as the bragg angle.
For example, lattice fringes corresponding to the (420) plane of V5S8 and the (002) plane of CoS, observed at a 50 ° angle (as shown in Fig. 2e), exhibited a lattice mismatch (m) of 42.25%. Notably, significant lattice distortions at the interfaces of bimetallic sulfides can generate an abundance of crystal defects, thereby enhancing conductivity and providing more active sites for Na+ storage[27]. The selected area electron diffraction (SAED) plot in Fig. 2g reveals multiple concentric diffraction rings associated with V5S8 and CoS, thus confirming the formation of the V5S8/CoS heterostructure. The yellow dashed ring pattern corresponds to the (222) and (-224) crystal planes of V5S8, while the green dashed ring pattern corresponds to the (101) and (110) crystal planes of CoS (Fig. S4). Elemental mapping in Fig. 2h clearly demonstrates the uniform distribution of S, V, and Co within the 0.2-V5S8/CoS-CNFs, consistent with the observations in the TEM images.
To determine the sample composition, both CNFs (as shown in Fig. S5) and 0.2-V5S8/CoS-CNFs (as shown in Fig. 3a) underwent XRD analysis to investigate their physical phases. It is noteworthy that the diffraction peaks of 0.2-V5S8/CoS-CNFs exhibit a strong correspondence with the characteristic peaks of V5S8 (ICSD, Card No. 97-065-1379) and CoS (ICSD, Card No. 97-062-4857), thereby confirming the successful incorporation of transition metal sulfides (TMSs) into the CNFs. Notably, these diffraction peaks are characterized by their sharpness and high intensity, indicating a high degree of crystallinity.
In the Raman spectra, the characteristic D peak at 1351.792 cm− 1 and G peak at 1590.824 cm− 1 are evident. The D peak is associated with the disordered structure of carbon, vacancies, and other defects, while the G peak is related to ordered sp2 hybridized graphitic carbon[28]. Notably, the intensity ratio of the D and G peaks is often indicative of the concentration of defects within the carbon structure[29]. In the Raman spectra of CNFs, the intensity ratio of the D and G peaks was measured as 1.42. Upon the incorporation of metal sulfides, the intensity ratio for 0.2-V5S8/CoS-CNFs increased to 1.92 (as shown in Fig. 3b). This higher intensity ratio of the D and G peaks in 0.2-V5S8/CoS-CNFs indicates a higher concentration of defects, which is beneficial for improving the rate performance and storage capacity of the material. The mass fractions of V5S8 and CoS within the 0.2-V5S8/CoS-CNFs were determined to be 16.78% and 14.93%, respectively, based on the weight loss curves (Fig. S6).
3.2 Electrochemical Performance of Z-V5S8/CoS-CNFs in Sodium-ion Battery
To further elucidate the advantages of sodium storage, a 2032 button cell was assembled using Z-V5S8/CoS-CNFs as the anode and sodium metal as the counter electrode. The results of CV for the first three cycles of 0.2-V5S8/CoS-CNFs are presented in Fig. 4a. For comparison, CVs of Z-V5S8/CoS-CNFs (Z = 0, 0.1, 0.15, 0.3) are shown in Fig. S7. In the initial cathodic scan, two distinct reduction peaks at 1.107 V and 0.803 V are observed. The first reduction peak at 1.107 V is attributed to the insertion of sodium ions into V5S8 and CoS, while the subsequent peak at 0.803 V can be ascribed to the further insertion of Na+ ions, resulting in the formation of V, Co, and amorphous Na2S during the conversion reaction[30]. However, these two peaks disappear in subsequent scans, likely due to the formation of the solid electrolyte interphase (SEI) and irreversible embedding of Na+[31].
In the subsequent anodic scans, three oxidation peaks at 1.416 V, 1.648 V, and 1.924 V are observed, corresponding to the multi-step reversible reactions leading to the formation of V5S8 and CoS. Additionally, a peak at 0.051 V suggests the insertion of Na+ ions into amorphous carbon[32]. Notably, CV curves exhibit good overlap after the first cycle, indicating that the electrochemical reaction of the 0.2-V5S8/CoS-CNFs electrode is reversible, and the material exhibits excellent structural stability. The above reactions can be described by the following equations[30, 33]:
Discharge process:
V5S8+xNa++Xe-=NaxV5S8 (2)
NaxV5S8+(16-x)Na++(16-x)e-=8Na2S+5V (3)
CoS+xNa++xe-=NaxCoS (4)
NaxCoS+(2-x)Na++(2-x)e-=Na2S+Co (5)
Charge process:
8Na2S+5V= NaxV5S8+(16-x)Na++(16-x)e- (6)
NaxV5S8=V5S8+xNa++Xe- (7)
Na2S+Co= NaxCoS+(2-x)Na++(2-x)e- (8)
NaxCoS=CoS+xNa++xe- (9)
Fig. 4b and Fig. S8 illustrate the discharge-charge curves of 0.2-V5S8/CoS-CNFs and Z-V5S8/CoS-CNFs (Z=0, 0.1, 0.15, 0.3) at a current density of 0.2 A g-1. Notably, the initial discharge-charge capacity of 0.2-V5S8/CoS-CNFs is 398.64/320.93 mAh g-1, showcasing an impressive initial coulombic efficiency (ICE) of 80.5 % during the first discharge-charge cycle. Importantly, the ICE of 0.2-V5S8/CoS-CNFs surpasses that of previously reported electrospun sodium electrode materials based on TMSs, as detailed in Table S1. Furthermore, the discharge-charge curves for the subsequent four cycles exhibit excellent overlap, providing compelling evidence of the outstanding cycling performance exhibited by 0.2-V5S8/CoS-CNFs.
The cycling performance of the Z-V5S8/CoS-CNFs electrode material at 0.2 A g-1 is depicted in Fig. 4c. The 0.2-V5S8/CoS-CNFs samples exhibited an initial specific capacity of 398.64 mAh g-1 during the first discharge and maintained a stable capacity of 233 mAh g-1 after 250 cycles. In contrast, CNFs, 0.1-V5S8/CoS-CNFs, and 0.15-V5S8/CoS-CNFs samples retained capacities of 112.76 mAh g-1, 168.46 mAh g-1, and 187.45 mAh g-1, respectively, after 250 cycles. These results indicate lower sodium storage capacity in CNFs, 0.1-V5S8/CoS-CNFs, and 0.15-V5S8/CoS-CNFs due to their lower V5S8 and CoS content. Additionally, the specific capacity of the 0.3-V5S8/CoS-CNFs sample exhibited a continuous decay during cycling and was lower than that of the 0.2-V5S8/CoS-CNFs sample. This decay may be attributed to the larger-sized V5S8 and CoS agglomerates, which had limited capacity release capabilities[23]. These findings highlight that maintaining moderate amounts (less than 31.71 %) of V5S8 and CoS uniformly dispersed and efficiently confined by CNFs effectively suppresses volume expansion of the active substance and stabilizes the electrode structure[34]. Conversely, when the V5S8 and CoS content in the fibers increased to 50.76 %, significant precipitation and aggregation of V5S8 and CoS occurred (Fig. S2), impeding the CNFs' ability to facilitate volume changes.
The electrolyte gradually decomposes on the surfaces of V5S8 and CoS during cycling, leading to a rapid deterioration of sodium storage capacity[35]. A similar trend is observed in the discharge capacity statistics of Z- V5S8/CoS-CNFs (Fig. S9). Fig. 4d presents the rate performance, where the 0.2- V5S8/CoS-CNFs samples deliver a high specific capacity of 415 mAh g-1 in the first cycle. Subsequently, at current densities of 0.2, 0.5, 1, 2, and 5 A g-1, the specific capacities are 310, 268, 253, 229, and 201 mAh g-1, respectively. Upon resetting the current density to 0.2 A g-1, the capacity rapidly recovers to 264 mAh g-1, showcasing excellent rate performance. This can be attributed to the appropriate carbon content, which ensures a stable carbon skeleton, suppresses volume expansion, and maintains material stability during electrochemical cycling.
Fig. 4e illustrates the relationship between discharged capacitance and voltage at different current densities, where the voltage plateau remains consistent at all current densities. This indicates low electrochemical polarization in the samples, a favorable trait for their application in energy storage devices. This behavior results from the stable heterostructures of the active interface and minimal electrochemical polarization phenomenon[36]. Notably, the comparison of 0.2-V5S8/CoS-CNFs with previously reported anode materials (Fig. 4f) [37-43] underscores its exceptional rate performance.
As expected, 0.2-V5S8/CoS-CNFs exhibit remarkable long-term cycling stability (Fig. 4g) with a Coulombic efficiency of nearly 100 %. To further validate their enduring stability, SEM images of 0.2-V5S8/CoS-CNFs before and after cycling are provided in Fig. S10. Clearly, the electrode fiber structure remains well-preserved after cycling, with no noticeable changes in diameter or metal precipitation or agglomeration. This confirms the efficacy of the fiber structure in suppressing volume changes during cycling, aligning with the experimental design.
To gain further insights into the reaction kinetics of 0.2-V5S8/CoS-CNFs, we conducted CV experiments at various scanning rates, as illustrated in Fig. 5a. The CV curves exhibit consistent oxidation and reduction potentials, which align with the peak position at 0.1 mV s-1, indicating that there are no significant changes in electrochemical reaction kinetics across different scan rates. The peak currents (i) of the redox peaks exhibit variations with scan rate (v), and their relationship conforms to Eq. (10), where the value of b can provide insights into the electrochemical mechanism[44]. A value of b close to 0.5 signifies the dominance of diffusion-controlled behavior, while a value near 1.0 suggests the predominance of capacitive behavior. To determine the value of b, Eq. (10) is logarithmically transformed into Eq. (11). As depicted in Fig. 5b, the b-values for peaks 1 to 5 are 0.80, 0.86, 0.98, 0.65, and 0.82, respectively. All these b values fall within the range of 0.5-1.0, indicating the coexistence of both diffusion-controlled and capacitive behaviors during charging and discharging processes. This suggests that the electrode material exhibits pseudocapacitive properties, where both diffusion and capacitance mechanisms contribute significantly to the electrochemical processes[45].
In the provided equations Eq. (12), k1v represents the contribution of pseudocapacitance, while k2v1/2 represents the contribution of ion diffusion. Fig. 5c illustrates the calculated pseudocapacitance contributions at different scan rates. The observed increase in capacitance contribution with higher scan rates confirms the rapid transport of Na+ ions within the material. As depicted in Fig. 5d, at a scan rate of 1.0 mV s-1, the surface capacitive behavior of the material accounts for 89.84 % of the contribution, with diffusion contributing only 10.16 %. The pseudocapacitance contributions at other scan rates are shown in Fig. S11. The larger capacitance contribution is advantageous for facilitating the charging and discharging processes of the material at high rates[46]. A comparison reveals that 0.2-V5S8/CoS-CNFs exhibit higher pseudocapacitance contributions compared to CNFs at the same scan rate (Fig. S12). This indicates that the inclusion of V5S8 and CoS enhances the pseudocapacitance effect, consequently improving the kinetics of sodium storage.
The reaction kinetics were further characterized through EIS tests. As depicted in Fig. 5e, the Nyquist plot consists of semicircles and slanted lines. The semicircles in the high-frequency region are associated with system impedance (Rs) and charge transfer impedance (Rct), while the slanted straight lines in the low-frequency region indicate diffusion behavior, related to the diffusion impedance of Na+[47, 48]. It is evident that electrodes with higher metal concentrations exhibit smaller Rs and Rct. Table S1 provides specific values of impedance after fitting EIS theory. As the metal concentration increases, the sample demonstrates enhanced electrical conductivity, thereby accelerating the transport of electrons and ions, and consequently, the reaction kinetics during the electrochemical process. Additionally, the EIS analysis was further extended to estimate the diffusion coefficient of Na+ ions (DNa+), based on the Z’-ω(-1/2) (ω=2πf) curves in the low-frequency region. DNa+ is calculated using the following equation[49]:
Where, R, T, A, F, n, and C represent the gas constant, absolute temperature, electrode surface area, Faraday constant, molar electron transfer number, and Na+ molar concentration, respectively. The Warburg coefficient σ is determined from the slope of the straight line. There is an observed trend of decreasing sodium ion diffusion coefficient with increasing metal concentration, which could potentially contribute to the capacity decay observed in concentrations up to 0.3-V5S8/CoS-CNFs.
The diffusion kinetics of 0.2-V5S8/CoS-CNFs were further investigated using GITT. DNa+ coefficient can be calculated by the following equation[50]:
where τ represents the relaxation time, ∆Es stands for the change in steady-state potential after the current pulse, and ∆Et represents the potential change (V) under a stable electric energy pulse without iR drop (Fig. 6b). L corresponds to the length of the Na+ diffusion path, which is simplified to the thickness of the electrode. Fig. 6a illustrates the discharge-charge GITT curve of the 0.2-V5S8/CoS-CNFs electrode. By calculating the diffusion coefficient of sodium ions during discharge charging (Fig. 6d), it can be observed that an increase in the metal sulfide content leads to a decrease in the Na+ diffusion coefficient, consistent with the results of the EIS tests. The diffusion coefficients ranged between 10-8 and 10-7 cm2 s-1, which are significantly higher than those of other materials in the same field (Table. S3). The excellent kinetic properties of the material may be attributed to its unique structural design, which facilitates the migration of sodium ions.