As seen in Fig. 1a, the crystal structure of the pristine Cu3SbSe4 is the tetragonal structure (I-42m space group) [38]. The X-ray diffraction (XRD) patterns of the synthesized Cu3Sb1 − xBixSe4 (x = 0–0.04) compounds are shown in Fig. 1b. The major diffraction peaks match well with the standard peaks of Cu3SbSe4 (JCPDS No.85 − 0003) phase, which belongs to the zinc blende structure with tetragonal system. As the concentration of Bi element increases to 0.03, some weak diffraction peaks at the 2θ angle of ~ 24o, 31o and 65o are observed, indicating the Bi and Cu2 − xSe phases [52]. Previous work also demonstrated that Bi phase was precipitated from the matrix as Bi concentration exceeds 0.02 in Cu3SbSe4 [45]. It can be deduced that the addition of Bi is propitious to the precipitation of Cu2 − xSe phase, which is beneficial to decrease in lattice thermal conductivity.
The morphologies and compositions of Cu3SbSe4-based samples were examined by SEM and EDS analyses. As depicted in Fig. 2, SEM images of the fracture surfaces (Fig. 2a-d) indicate that Cu3SbSe4 materials are isotropic and some impurity phases are embedded into the Cu3SbSe4 matrix. EDS result demonstrates that the impurity phase contains of Cu and Se elements (Fig. 2e) and the atomic ratio (%) of copper to selenium is 61.31: 32.56, implying that the impurity phase is copper selenide of Cu2 − xSe, which is consistent with XRD results. The back-scattered electron (BSE) image of Cu3Sb0.985Bi0.015Se4 sample clearly present the micro or even nanoscale impurity phases (Fig. S1), which are useful to reduce κlat by scattering long-wave phonons [46]. The Cu2 − xSe with extremely low κlat can decrease the κtot of Cu3SbSe4 to some degree [52, 53].
XPS is used to check the valence state of elements and chemical composition of the Cu3SbSe4-based samples [54]. The experimental data of XPS spectra peaks of Cu3Sb0.985Bi0.015Se4 samples are presented in Fig. 3. XPS spectra confirm the presence of Cu, Sb, Se and Bi elements in the Cu3SbSe4-based compounds. As shown in Fig. 3a, Cu 2p core electron peak of specimen is split into two peaks, and the XPS spectra of Cu 2p3/2 (931.2 eV) and the Cu 2p1/2 (952.5 eV) peaks can confirming indicate the existence of Cu1+ in the sample [55]. The Sb 3d spectra can be deconvoluted into two peaks of Sb 3d 5/2 (528.4 eV) and Sb 3d 3/2 (537.7 eV) (Fig. 2b), which can be indexed to Sb5+ [56]. The Se 3d spectra can be also split into two peaks at 53.1 eV (3d5/2) and 54.0 eV (3d3/2) (Fig. 2c), which are matched with Se2− of oxidation state [57]. The Bi 4f spectra do not display distinguished changes. Consequently, the components of as-sintered Bi-doped samples are speculated to be consisted of Cu3SbSe4, Cu2 − xSe, which is in agreement with XRD and SEM results.
The presence of secondary phase is clearly visualized by employing a high-resolution TEM. The low-magnification images presented in Fig. 4a demonstrate that some impurity phases are distributed in the matrix. As shown in Fig. 4b, a crossed fringe, with interplanar spacing values of 0.5164 nm, which can be indexed to (101) planes of Cu3SbSe4 [40]. The 0.2986 nm, 0.1041 nm and 0.1157 nm of spacing are in response to (200), (511) and (422) plane of Cu2Se, respectively [53]. It can be seen that some amorphous phases exist in the matrix (yellow A presents amorphous phase and C presents crystallized phase). In addition, the interfaces (blue dot lines) among Cu3SbSe4, Cu2Se and amorphous phases can be clearly observed. All of which are good for reducing lattice thermal conductivity. The SAED pattern, corresponding to (105), (220) planes of Cu3SbSe4, and (111), (220) planes of Cu2Se, which may be further supposed the as-prepared sample is consisted of Cu3SbSe4 and Cu2 − xSe [52–54].
To validate our experimental design and investigate the change of electronic properties, the density functional theory (DFT) calculations were performed [13]. To ensure the optimal doping model is consistent with the experiment and guarantee the rationality for theory calculation, the Cu48Sb16Se64 and Cu48Sb15Bi1Se64 structures were constructed. The crystal structure of Bi-doped Cu48Sb16Se64 (including the top and side views) is shown in Fig. S2. Bi-doped Cu48Sb16Se64. Bi can form four 2.74Å Bi-Se bonds with the surrounding Se atoms, which is 0.2 Å longer than the original Sb-Se bond, mainly caused by the larger atomic radius of the Bi atom. The DOS of pristine and Bi doped Cu48Sb16Se64 is presented in Fig. 5a. After Bi doping, the Fermi level changes from 3.99 eV to 3.77 eV, thus improving electrical transport characteristics. The band gap (Eg) of Cu48Sb16Se64 is 0.25 eV, which is basically consistent with previous report [58], and increase to 0.27 eV for Cu48Sb15Bi1Se64. More importantly, the DOS near the Fermi level (0 eV) are increased, which is composed of the hybridized Cu 3d, Sb 5s, Se 4p and Bi 6s orbitals (Fig. 5b), leading to large Seebeck coefficients
The temperature-dependent electrical transport properties of Cu3SbSe4-based samples are depicted in Fig. 6. As seen in Fig. 6a, the ρ of pristine Cu3SbSe4 sample decreases with elevating temperature, indicating that it is non-degenerate semiconductor. However, Bi-doped samples become partially degenerate semiconductors, especially the samples with low Bi content (x < 0.015). Clearly, the ρ of all Bi-doped Cu3SbSe4 compounds is lower than that of pristine Cu3SbSe4 due to the increased carrier concentration (Tab. S1). In comparison with Sb5+ (0.62 Å), the bigger radius of Bi5+ (0.74 Å) can generate Cu and/or Sb vacancies [45]. Meanwhile, the ρ of the samples increases with increasing Bi content (x < 0.3), however, the ρ sharply decreases as x = 0.03 and 0.04 because of the relatively high carrier mobility and carrier concentration. In addition, the ρ shows the temperature-independent behavior at high temperature because the intrinsic excitation gradually becomes dominant at 573–673 K, which is similar to the transport behavior of non-degenerate semiconductor.
Figure 6b demonstrates the temperature-dependent S of Cu3SbSe4-based samples. It can be seen that they are p-type conductivity derived from the positive values of S, and the charge carriers in these compounds are holes. After Bi doping, the S first decreases and then rises with increasing temperature, achieving a minimum S, which should be attributed to the α-β phase transition of Cu2 − xSe phase [52, 53]. Although the increased hole concentration (Tab. S1) deteriorates S, Bi doping can enhance DOS near Fermi level (Ef) and enlarge Eg according to the DFT calculations, resulting in the relatively large S [45].
Figure 6c presents the temperature-dependent PF of Cu3SbSe4-based samples. It can be seen that the PF values of pristine Cu3SbSe4 sample increase with elevating temperature. Besides, PF values of Bi-doped compounds decrease slightly and then increases with rising temperature, and are higher than that of pristine Cu3SbSe4 at high temperature, which is attributed to the enhancement of electrical conductivity after Bi-doping. However, the excessive Bi content (x ≥ 0.02) does not further enhance PF values due to the deteriorated Seebeck coefficient. The maximum PF value of ~ 1000 µWcm− 1K− 2 is achieved at 673 K for Cu3Sb0.985Bi0.015Se4 compound.
(c) power factor of Cu3SbSe4-based samples.
The temperature-dependent thermal transport properties of the Cu3SbSe4-based samples are displayed in Fig. 7. The D decreases with the increase of temperature and the room-temperature D drops from ~ 1.7 mm2s− 1 for pristine Cu3SbSe4 to ~ 1.45 mm2s− 1 for Cu3Sb0.97Bi0.03Se4 sample after Bi doping (Fig. 7a). As presented in Fig. 7b, the κtot reduces with rising temperature within the range of measurement due to the enhanced lattice vibration scattering at high temperature [11, 59]. Besides, the κtot of all the Cu3SbSe4-based samples decreases after Bi doping, drops from ~ 3.11 Wm− 1K− 1 at 300 K and ~ 1.03 Wm− 1K− 1 at 673 K for pristine Cu3SbSe4 to ~ 2.81 Wm− 1K− 1 at 300 K and ~ 0.69 Wm− 1K− 1 at 673 K for Cu3Sb0.985Bi0.015Se4 sample. Generally, the κlat can be calculated by subtracting the electronic part κele from the κtot. According to the Wideman–Franz law Eq. (1):
where, L is the Lorenz number and can be calculated from Eq. (2) [60, 61]:
The calculated L values of the Cu3SbSe4-based samples are presented in Fig. S3 and L ranges from 1.51 to 2.13 W Ω K− 2. The increased κele of Bi-doped compounds is shown in Fig. 7c, which can be closely related to electrical conductivity derived from the rising carrier concentration and high carrier mobility. As seen in Fig. 7d, the prominent scattering mechanism of Cu3SbSe4-based samples is the phonon-phonon Umklapp scattering, which can be inferred from the relationship of κlat ∝ T− 1. As expected, the Bi-doping can decrease κlat in the measured temperature region. On one hand, the isovalent substitutions on Sb can reduce κlat because of stronger mass and strain field fluctuations, which is similar to the case that InSb point defects in Cu3SbSe4 can scatter high frequency phonons [39]. On the other hand, the heavy element Bi has much larger ionic radius in comparison with Sb, leading to the reduced κlat. The heat-carrying acoustic phonons can be blocked by softening the chemical bonding interaction in crystal due to the expansion of unit cell [46]. Furthermore, the impurity phase Cu2Se and amorphous phase with extremely low κlat can also decrease κlat of Bi-doped Cu3SbSe4 [52, 53]. Consequently, the reduction of κtot should be attributted to the depressed κlat caused by the strengthened phonon scattering of multiscale defects involving point defects, nanoprecipitates, amorphous phases and grain boundaries [62, 63]. The lowest κlat ~ 0.53 Wm− 1K− 1 and κtot ~ 0.62 Wm− 1K− 1 is achieved at 673 K for the Cu3Sb0.98Bi0.02Se4 sample.
The temperature-dependent ZT of Cu3SbSe4-based samples with various Bi concentration is presented in Fig. 8a. Benefiting from the synergistically optimized electrical and thermal transport properties, the ZT values obtains an enhancement at high temperature for Bi-doped compounds, which can reach the maximum ZT value of ~ 0.95 at 673 K for the Cu3Sb0.985Bi0.015Se4 sample. To further analyze our experimental results, the previous ZTmax values of other single-doped Cu3SbSe4 samples are listed in Fig. 8b [13, 38–45]. Accordingly, isovalent heavy element Bi doping in Cu3SbSe4 compounds can simultaneous optimize electrical and thermal parameters, which is highly effective in improving TE properties.