Open Access Article
Nurfarhana Ahmad Musri*a,
Yoganash Putthisigamanya,
Puvaneswaran Chelvanathana,
Nadhrah Md Yatimb,
Farah Liana Mohd Redzuanc and
Ubaidah Syafiq
*a
aSolar Energy Research Institute (SERI), Universiti Kebangsaan Malaysia, Bangi, Selangor, Malaysia. E-mail: ubaidahsyafiq@ukm.edu.my
bFaculty of Science and Technology, Universiti Sains Islam Malaysia, Bandar Baru Nilai, Negeri Sembilan, Malaysia
cMalaysia-Japan International Institute of Technology, Universiti Teknologi Malaysia, Kuala Lumpur, Malaysia
First published on 7th October 2025
Thin film thermoelectric generators (TEGs) offer a compact and scalable solution for harvesting waste heat in microelectronic and flexible systems. However, their performance is highly influenced by material properties, fabrication parameters, and device geometry. This work presents a systematic optimization of Bi2Te3 and Sb2Te3 thin films for a high-performance TEG. The thin films were deposited via radio frequency (RF) magnetron sputtering under varying power conditions and characterized using XRD, FESEM, EDX, Hall measurements, and Seebeck coefficient analysis. The optimal sputtering powers were identified as 75 W for Bi2Te3 and 30 W for Sb2Te3, yielding high crystallinity and balanced electrical conductivity with favourable Seebeck values. A thin film TEG consisting of six n-type/p-type leg pairs was fabricated using these conditions and subjected to post-deposition annealing to improve performance. The TEG annealed at 200 °C demonstrated a peak power output of 0.84 μW at ∼120 °C ΔT, indicating enhanced crystallinity and reduced internal resistance. Dimensional optimization further revealed that wider and shorter TE legs significantly improve the output by minimizing internal resistance. The results highlight the importance of integrating material, process, and device-level optimization for the development of efficient and scalable thin film thermoelectric systems.
Bismuth telluride (Bi2Te3) and antimony telluride (Sb2Te3) are among the most extensively studied thermoelectric materials for near-room-temperature applications, owing to their high Seebeck coefficients, moderate electrical conductivity, and inherently low thermal conductivity.9,10 The efficiency of a thermoelectric material is governed by its figure of merit (ZT), which is strongly influenced by electrical conductivity (σ), Seebeck coefficient (S), and thermal conductivity (κ), defined as ZT = (S2σT)/κ, where S2σ can be defined as power factor (PF). Numerous studies have focused on enhancing the thermoelectric performance to widen its potential and application. Recent efforts have shifted toward thin-film configurations, where these properties can be finely tuned through nanoscale control, improved carrier mobility, and phonon scattering mechanisms. Compared to bulk counterparts, thin-film TEGs offer several advantages, including enhanced power density, mechanical flexibility, and better integration with microelectronic and wearable systems.11,12
However, achieving optimal performance in thin-film TEGs requires a multifaceted approach. Key factors such as sputtering parameters, post-deposition treatments, and geometric design must be precisely engineered to maximize electrical output while minimizing internal resistance and thermal losses. Recent advances in Bi2Te3 and Sb2Te3 thin films have demonstrated that careful control of film microstructure, stoichiometry and strain can markedly improve thermoelectric performance. For example, Zheng et al.13 demonstrated a textured, flexible Bi2Te3 thin-film design that achieved record room-temperature ZT values (∼1.2) while retaining mechanical robustness, underscoring the importance of crystallographic texture and film processing. In Sb2Te3-based systems, strategies such as lattice-strain engineering and heterostructuring have been shown to enhance the Seebeck coefficient and power factor by modifying carrier scattering and phonon transport.14 Other than that, Haidar et al.15 investigated the effects of sputtering temperature, co-sputtering techniques, and annealing treatments on the structural and thermoelectric properties of these materials for micro-scale TEGs. While their work provides valuable insights into the material-level optimization for vertical μTEG architectures, it does not address planar thin-film device configurations. Nevertheless, collectively, these material-level developments motivate our integrated experimental study to optimize Bi2Te3/Sb2Te3 thin film TEG performance.
Prior studies have shown that device geometry can strongly influence output power and power density. Kalil et al.16 conducted a detailed performance comparison of bulk TEGs by analyzing different leg geometries under constant leg volume, demonstrating that variations in leg shape and aspect ratio significantly affect power output and thermal management. However, their study focuses exclusively on bulk TEGs, where heat flow and mechanical constraints differ considerably from those in thin-film systems. Several studies combine numerical optimization with experiments to identify optimal leg geometries for thin film or wearable TEGs. For instance, Kuang et al.17 investigated annular thin film TEG architectures and highlighted the combined role of material properties and geometry on device performance, reporting that annealing improves the film Seebeck coefficient and conductivity, while geometry sets the achievable output power density. Mathematical and numerical studies by Newbrook et al.18 further showed that fill factor and interconnect conductivity critically determine thin film TEG performance. Building on this line of work, this study extends the understanding of thin film TEG design by focusing on planar devices, fabricated using optimized sputtering and annealing conditions, and introduces dimensional optimization of the thermoelectric legs.
In this work, we present a comprehensive investigation on the optimization of the deposition parameters for the fabrication of Bi2Te3/Sb2Te3 thin film TEGs. This study explores the influence of sputtering power on film crystallinity and thermoelectric properties, the impact of post-deposition annealing on device performance, and the role of leg geometry in enhancing power output. By systematically examining each of these factors, this research aims to establish a set of process guidelines for the fabrication of high-performance thin film thermoelectric devices, contributing towards the development of a scalable and cost-effective energy harvesting solution.
X-ray diffraction (XRD, X’pert-pro PANalytical X-ray diffractometer) using Cu Kα radiation (λ = 1.5406 Å) was used to investigate the phase structure of all the as-deposited thin films. Rietveld refinement of the experimental XRD data was done to calculate lattice parameters a and c. The field-emission scanning electron microscopy (FESEM/JEOL, JSM7800F, operating at 3.0 kV) topographic and cross-sectional images were used to evaluate the surface microstructure and film thickness. Energy dispersive X-ray spectroscopy (EDX) attached to the FESEM was used to confirm the film composition. Moreover, the Hall effect measurement system (HMS ECOPIA 3000) with a fixed magnetic field of 0.57 T was employed to test the carrier concentration and electrical conductivity. The temperature-dependent in-plane S was measured using a Linseis LSR-3, with rectangular samples (about 2 cm × 1.25 cm) mounted for measurement. The absolute S was determined in a two-contact configuration relative to a Pt standard, applying a temperature gradient of approximately 25 °C.
![]() | ||
| Fig. 2 (a) The fabricated Bi2Te3/Sb2Te3 thin film TEG. (b) Thermal image of a Bi2Te3/Sb2Te3 thin film TEG in operation. (c) Picture of the custom IVP measurement setup. | ||
m (166) with lattice parameters a = 4.3852 Å and c = 30.4830 Å. The observed reflections illustrate preferential growth along the c-axis with a sharp dominant peak growing anisotropically along the (0 1 5) direction, indicating a highly crystalline structure with a few small peaks corresponding to the (0 0 6), (1 0 10) and (1 1 0) directions. The intensity of these peaks systematically increases with RF power, indicating enhanced crystallinity as the RF power increases from 75 W to 135 W. This trend suggests that a higher RF power provides sufficient energy to atoms, promoting surface mobility and the formation of larger, well-ordered grains. This trend can also be seen more clearly in the Sb2Te3 spectra.
![]() | ||
| Fig. 3 XRD patterns of (a) Bi2Te3 films deposited at 75, 90, 105, 120, and 135 W and (b) Sb2Te3 films deposited at 10, 20, 30, 40, and 50 W. | ||
For Sb2Te3 samples (Fig. 3(b)), the peaks aligned well with the respective rhombohedral Sb2Te3 phase (PDF #01-071-0393) in the space group R
m (166) with lattice parameters a = 4.2640 Å and c = 30.4580 Å. The absence of extraneous peaks in the spectra confirmed that the films are phase-pure Sb2Te3 with no detectable secondary phases or impurities. The peaks are relatively small and broad, corresponding to the (0 1 5) direction especially at a low RF power due to the low thickness of the film, as proven by the FESEM cross-section image to be discussed. The intensity of these peaks also increases with RF power from 10 W to 50 W. At lower powers, the peaks are broader and less intense, reflecting poorer crystallinity and smaller crystallite sizes. In contrast, the peaks at 50 W are the sharpest with a notable increment at (1 1 0), highlighting optimal conditions for achieving high-quality crystalline films. Overall, the systematic improvement in peak sharpness and intensity with increasing power highlights the critical role of power optimization in determining the structural properties of the films.
Table 1 presents the crystallite sizes of Bi2Te3 and Sb2Te3 thin films sputtered at different RF powers. A clear trend emerges where the crystallite size increases from 8.52 nm (sample A1) to 9.59 nm (sample E1) as the sputtering power is increased. This behaviour can be attributed to the higher energy imparted to the sputtered atoms at elevated RF powers, which enhances atom mobility on the substrate surface.22 As a result, the atoms have more opportunity to diffuse, leading to larger grains during film growth. Although the overall differences in crystallite size are within a narrow range, even slight increases in grain size can influence the electrical and thermal properties of Bi2Te3. Consequently, selecting an appropriate sputtering power becomes critical for optimizing film quality and, by extension, the thermoelectric performance of the deposited layers. A similar trend emerges in the Sb2Te3 films: the crystallite size increases from 0.91 nm (sample A2) at 10 W to 3.60 nm (sample E2) at 50 W. Compared to Bi2Te3, Sb2Te3 films exhibit notably smaller crystallite dimensions. This discrepancy can be attributed to several factors, including the lower overall power range used for Sb2Te3 deposition, differences in material composition, and variations in growth kinetics. These results underscore the importance of tailoring sputtering parameters to each specific material system to achieve the desired microstructure and, consequently, optimal thermoelectric performance.
| Material | Sample | RF power, W | Crystallite size, nm | Thickness, μm |
|---|---|---|---|---|
| Bi2Te3 | A1 | 75 | 8.52 | 1.429 |
| B1 | 90 | 8.57 | 1.742 | |
| C1 | 105 | 8.88 | 1.909 | |
| D1 | 120 | 9.46 | 1.945 | |
| E1 | 135 | 9.59 | 2.077 | |
| Sb2Te3 | A2 | 10 | 0.91 | 0.139 |
| B2 | 20 | 1.28 | 0.365 | |
| C2 | 30 | 2.46 | 0.458 | |
| D2 | 40 | 3.08 | 0.640 | |
| E2 | 50 | 3.60 | 0.938 | |
The FESEM micrograph shows that the as-deposited Bi2Te3 and Sb2Te3 films cover the entirety of the substrate in a continuous and homogenous manner. As can be seen in Fig. 4(a), Bi2Te3 films comprise hexagonal grains, confirming the crystal structure of Bi2Te3. The grains are densely packed with varying sizes, indicating a polycrystalline nature. The surface of the Sb2Te3 film (Fig. 4(b)) exhibits more fine, rounded grains and is consistent in size compared to that of the Bi2Te3 film, indicating a nearly amorphous phase resulting in a low and broad peak in XRD. The FESEM surface images presented correspond to Bi2Te3 deposited at 75 W for 60 minutes and Sb2Te3 deposited at 30 W for 120 minutes. Since the surface morphology of the other films exhibited similar structural characteristics, their images are not included in the main text. Likewise, the FESEM cross-section images used for thickness estimation are omitted, as the measured values are summarized in Table 1. The complete set of FESEM images for all samples is available in the supplementary file for reference. From the cross-sectional views shown in Fig. 4(b) and (e), the interface between the film and the substrate looks distinct, suggesting good adhesion properties. The thickness of the Bi2Te3 and Sb2Te3 films shows a clear correlation with RF power, where a higher RF power results in increased film thickness. For Bi2Te3, the thickness ranges from 1.429 μm at 75 W to 2.077 μm at 135 W, indicating a steady increase with increasing power. Similarly, for Sb2Te3, the thickness increases from 0.139 μm at 10 W to 0.938 μm at 50 W. The observed trends are consistent with typical RF sputtering behaviour, where increased power enhances the kinetic energy of the sputtered species, leading to higher deposition rates.23,24 However, while thicker films may improve mechanical stability, they can also influence electrical and thermal transport properties, which must be carefully optimized for thermoelectric applications.
The EDX spectra of both Bi2Te3 and Sb2Te3 thin films provide crucial insights into the elemental composition and purity of the deposited materials. In the case of Bi2Te3, the spectrum shows prominent peaks corresponding to bismuth (Bi) and tellurium (Te), with their respective weight percentages of 51.0% and 42.8%. Additionally, small traces of carbon (4.2%) and oxygen (2.0%) are observed, likely due to surface contamination or oxidation. The composition of the films was calculated from the EDX spectra, shown in Fig. 4(c), giving an atomic percent of 41.9% Bi: 58.1% Te. The near-stoichiometric Bi–Te ratio suggests a well-formed Bi2Te3 structure, which is desirable for thermoelectric applications. Similarly, the Sb2Te3 spectrum exhibits dominant peaks of antimony (Sb) and tellurium (Te), with weight percentages of 38.1% and 61.9%, respectively. The composition of the films was calculated from the EDX spectra, shown in Fig. 4(f), giving an atomic percent of 40.0% Sb:60.0% Te. The estimated values show that the film gave a consistent Sb:Te ratio of 2
:
3, indicating that stoichiometric Sb2Te3 was successfully fabricated. The EDX spectra presented here are representative of the Bi2Te3 and Sb2Te3 thin films, as the atomic compositions for the other samples exhibit similar elemental distributions. Due to this consistency, the EDX results for the remaining samples are not included in the main text but will be provided in the supplementary file for reference.
![]() | ||
| Fig. 5 Conductivity and carrier concentration of (a) n-type Bi2Te3 (shown as negative values to indicate electron carriers) and (b) p-type Sb2Te3 thin films. | ||
The carrier concentration, represented by the blue squares, exhibits an increasingly negative trend with a higher RF power. The carrier concentration measured using Hall effect measurements is derived from the Hall coefficient (RH). The sign and magnitude of the Hall coefficient provide crucial information about the type and density of charge carriers present in the material.27 Specifically, a positive Hall coefficient indicates that the majority charge carriers are holes, identifying the material as p-type, whereas a negative Hall coefficient signifies that majority charge carriers are electrons, indicating n-type conductivity. Hence, the negative magnitude in the figure indicates the n-type conductivity of Bi2Te3, and the carrier concentration value is actually increasing at a higher RF power. At a lower RF power (75 W, sample A1), the carrier concentration is approximately 5.71 × 1020 cm−3, and it changes to around 2.03 × 1021 cm−3 at 135 W (sample E1).
The same trends were observed for Sb2Te3 thin films (Fig. 5(b)). Both σ and carrier concentration increase consistently with a higher RF power. Specifically, sample A2 shows the lowest σ and carrier concentration, while sample E2 demonstrates the highest values for both parameters. The magnitude of the carrier concentration shows positive values, indicating p-type conductivity for all samples. Overall, the results confirm that increasing RF power enhances the structural and electrical properties of Sb2Te3 thin films by improving crystallinity and increasing thickness, leading to higher σ and carrier concentration.
The S plot of Bi2Te3 thin films (Fig. 6(a)) deposited at various RF powers demonstrates the dependence of thermoelectric properties on deposition parameters. The negative values of absolute S across all samples confirmed that the films exhibit n-type conductivities, where electrons serve as the dominant charge carriers, consistent with Hall effect measurements. The S value increases with increasing temperature for all samples. This trend indicates a favourable thermoelectric response with temperature, which is associated with grain growth at elevated temperatures,28 leading to reduced carrier scattering.29However, it should be noted that a comprehensive evaluation of thermoelectric performance requires consideration of electrical and thermal conductivities in addition to S.
At near room temperature, sample A1 shows the highest absolute S, approximately −60 μV K−1, indicating superior thermoelectric performance compared to other samples. The higher magnitude of the S suggests a more significant contribution of charge carrier diffusion to the thermoelectric voltage, which is beneficial for energy conversion efficiency. The higher S of sample A1 can be correlated with its relatively lower carrier concentration, as shown in Fig. 5(a). According to the Pisarenko relation, the S is inversely related to carrier concentration, whereby a reduction in carrier density enhances the energy dependence of carrier transport and thus increases the Seebeck value.30 A lower RF power during deposition resulted in a finer grain structure and a moderate carrier concentration, balancing σ and thermoelectric performance. In contrast, samples deposited at higher RF powers, such as sample E1, show reduced S, due to increased carrier concentrations that lower the thermoelectric voltage.
Fig. 6(b) demonstrates the S plot of Sb2Te3 thin films, deposited using RF powers ranging from 10 W (sample A2) to 50 W (sample E2). The positive values of S for all samples confirm the p-type conductivity of the thin films, where holes act as the primary charge carriers. This observation aligns with the intrinsic behaviour of Sb2Te3, a well-known p-type thermoelectric material. Sample A2 exhibits an unusually high S near room temperature, exceeding the expected range, due to instrument limits. The conductivity of sample A2 is too low, falling below the resistance limit of the Seebeck instrument. As a result, the measured value lies outside the expected range, making the said measurement unreliable. However, as the temperature increases, the S stabilizes and shows improved consistency, suggesting that the carrier transport mechanisms become more favourable at higher temperatures. This trend indicates a strong dependence of thermoelectric performance on the deposition conditions, with a low RF power significantly affecting carrier concentration or scattering. The S of the Sb2Te3 thin films exhibits a distinct trend where it decreases with increasing temperature up to a certain point and then starts to increase again at higher temperatures. These results suggest a relatively stable thermoelectric performance for Sb2Te3, although the slight decrease indicates some temperature dependence, attributed to the amorphous–crystalline transition. Similar behaviour has also been reported by Fang et al.30 and V. Damodara et al.31,33
The σ and S of the thin films were determined to calculate the power factor (PF) of each film. PF becomes a decisive parameter for the selection of thin film materials due to its configuration. A low thermal conductivity is ideal for high efficiency; however, the substrate typically dominates thermal transport in thin film configurations because of the significant difference in thickness.32 Table 2 shows the calculated power factors of each sample at ambient temperature.
| Material | Sample | Conductivity, σ (Ω cm)−1 | Absolute Seebeck coefficient, S (μV K−1) | Power factor, PF (μW K−2 cm−1) |
|---|---|---|---|---|
| Bi2Te3 | A1 | 108.96 | 72.84 | 0.5781 |
| B1 | 150.6 | 62.02 | 0.5793 | |
| C1 | 239.46 | 69.43 | 1.1543 | |
| D1 | 295.92 | 64.03 | 1.2132 | |
| E1 | 360.46 | 65.33 | 1.5384 | |
| Sb2Te3 | A2 | 0.01067 | (OOR) | |
| B2 | 0.32594 | 423.43 | 0.1584 | |
| C2 | 1.0696 | 336.11 | 0.1208 | |
| D2 | 1.9206 | 301.44 | 0.1745 | |
| E2 | 3.055 | 241.56 | 0.1783 | |
The S and σ values greatly influenced the PF of thermoelectric materials, calculated by PF = S2σ. However, these two properties exhibit an inverse relationship due to their mutual dependence on carrier concentration, where an increase of one promotes a decrease of the other.33 The trend in the data reflects the well-known trade-off between these two properties.
For Bi2Te3, both σ and S exhibit variations across different samples, influencing the power factor. The trend shows that as conductivity increases (from 108.96 (Ω cm)−1 in sample A1 to 360.46 (Ω cm)−1 in sample E1), S generally decreases (from 72.84 μV K−1 to 65.33 μV K−1). However, this reduction in S is not steep enough to significantly lower the power factor. As a result, the power factor increases progressively from sample A1 (0.5781 μW K2 cm−1) to sample E1 (1.5384 μW K2 cm−1). This suggests that increasing σ plays a dominant role in improving the thermoelectric performance of Bi2Te3.
For Sb2Te3, the opposite trend is observed compared to Bi2Te3. The S is significantly higher in samples with lower σ, particularly in sample A2, but σ is very low (0.01067(Ω cm)−1), leading to a very small power factor (0.0163 μW K2 cm−1). As σ increases from sample A2 to sample E2, S decreases sharply, from 1236.84 μV K−1 to 241.56 μV K−1. While σ increases significantly (from 0.01067 to 3.055 (Ω cm)−1), the strong reduction in S limits the power factor improvement. Consequently, the power factor increases but at a much lower rate than in Bi2Te3, reaching a maximum of 0.1783 μW K2 cm−1 in sample E2. These results highlight that optimizing thermoelectric performance depends on finding a balance between σ and S.
The results obtained above guide the selection of the deposition parameters for device fabrication. While certain parameters yield higher PF values, other factors such as stability, reproducibility, and temperature-dependent behavior must also be considered. Achieving a balance between these aspects ensures that the fabricated device performs efficiently across a range of operating conditions. For the fabrication of the TEG, sample A (Bi2Te3) was selected due to its higher S near room temperature, which is critical for applications operating across ambient-to-moderate temperature gradients. This characteristic ensures sustained thermoelectric performance and contributes to a higher thermoelectric figure of merit (ZT). For the Sb2Te3 layer, sample C was chosen owing to its stable and consistent thermoelectric behaviour across the measured temperature range, offering a reliable balance between efficiency and operational stability. Although some samples exhibited higher PF, the selection was primarily guided by S trends, given that Hall measurements were limited to room temperature. The chosen samples demonstrated an optimal combination of initial S values, thermal stability, and structural integrity, ensuring enhanced efficiency, reliability, and practical applicability of the optimized TEG device.
![]() | ||
| Fig. 8 (a) IV and (b) IP plots of the optimized Bi2Te3/Sb2Te3 thin film TEG, annealed at 100 and 200 °C. | ||
Fig. 8(b) exhibits the IP plot, illustrating the power output as a function of current for the optimized Bi2Te3/Sb2Te3 thin-film TEG annealed at different temperatures. The curves show distinct power–current characteristics, indicating the influence of annealing temperature on the TEG's performance. Among the three samples, the TEG annealed at 200 °C (A200) exhibits the highest power output, peaking at 0.84 μW at around 0.012 mA. This indicates that higher annealing temperatures enhance the σ and overall thermoelectric efficiency of the device. This improvement is due to enhanced crystallinity, reduced grain boundary scattering, and optimized carrier concentration, which collectively contribute to better charge transport properties, confirming that the reduced internal resistance contributes to improved electrical performance. In contrast, the TEG annealed at 100 °C (A100) demonstrates a slight reduction in power output compared to the optimized sample, 0.25 μW at around 0.003 mA, aligning with the IV results where increased resistance may have caused power dissipation. This behavior can be correlated with the S trends shown in Fig. 6, where Bi2Te3 (sample A1) shows only a modest increase in Seebeck value when the temperature increases to 100 °C, while Sb2Te3 (sample C2) experiences a more pronounced reduction. Since the overall device voltage and power are determined by the combined contributions of both p-type and n-type legs, the significant degradation in the Sb2Te3 S offsets the minor improvement in Bi2Te3, resulting in a net decrease in thermoelectric performance at 100 °C. Thus, annealing at 200 °C appears to optimize the material properties, leading to enhanced thermoelectric efficiency.
Fig. 9(b) shows the IP curve of Sb2Te3-only thin film TEGs with different dimensions. The wider leg (4.5 mm) produces the highest power output, peaking at 0.19 μW, due to its lower internal resistance, which allows for a higher current flow. On the other hand, narrower legs (1.5 mm) exhibit a significantly lower power output compared to the standard dimension (0.11 μW), peaking at 0.07 μW due to their higher resistance, which restricts current flow and reduces the overall efficiency of the device. For longer legs, the maximum power is reduced, 0.10 μW for 25 mm and 0.05 μW for 35 mm compared to the standard dimensions, as their increased resistance limits current flow despite maintaining a good temperature gradient. From the IVP plot, it is evident that wider legs provide better σ, while longer or narrower legs increase resistance, limiting current output. The results from the dimensional optimization provided valuable and clear insights into the impact of dimension control on the performance of the TEGs, providing a crucial foundation for future research.
| This journal is © The Royal Society of Chemistry 2025 |