Next Article in Journal
Magnetoelectric Coupling Effects in Tb-Doped BiFeO3 Nanoparticles
Next Article in Special Issue
Short-Period Skyrmion Crystals in Itinerant Body-Centered Tetragonal Magnets
Previous Article in Journal
Controlling the Magnetic Properties of La0.9A0.1Mn0.9Cr0.1O3 (A: Li, K, Na) Powders and Ceramics by Alkali Ions Doping
Previous Article in Special Issue
Interacting with Futuristic Topological Quantum Materials: A Potential Candidate for Spintronics Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Conductivity for p–(Bi, Sb)2Te3 Films of Topological Insulators

Ioffe Institute, Russian Academy of Sciences, 26 Politekhnicheskaya, St. Petersburg 194021, Russia
*
Author to whom correspondence should be addressed.
Magnetochemistry 2023, 9(6), 141; https://doi.org/10.3390/magnetochemistry9060141
Submission received: 20 January 2023 / Revised: 14 May 2023 / Accepted: 22 May 2023 / Published: 25 May 2023
(This article belongs to the Special Issue Spin and Charge Transport in Novel Quantum and Topological Materials)

Abstract

:
In this study, we investigated the temperature dependencies of the total, crystal lattice, and electronic thermal conductivities in films of topological insulators p–Bi0.5Sb1.5Te3 and p–Bi2Te3 formed by discrete and thermal evaporation methods. The largest decrease in the lattice thermal conductivity because of the scattering of long-wavelength phonons on the grain interfaces was observed in the films of the solid-solution p–Bi0.5Sb1.5Te3 deposited by discrete evaporation on the amorphous substrates of polyimide without thermal treatment. It was shown that in the p–Bi0.5Sb1.5Te3 films with low thermal conductivity, the energy dependence of the relaxation time is enhanced, which is specific to the topological insulators. The electronic thermal conductivity was determined by taking into account the effective scattering parameter in the relaxation time approximation versus energy in the Lorentz number calculations. A correlation was established between the thermal conductivity and the peculiarities of the morphology of the interlayer surface (0001) in the studied films. Additionally, the total κ and the lattice κL thermal conductivities decrease, while the number of grains and the roughness of the surface (0001) increase in unannealed films compared to annealed ones. It was demonstrated that increasing the thermoelectric figure of merit ZT in the p–Bi0.5Sb1.5Te3 films formed by discrete evaporation on a polyimide substrate is determined by an increase in the effective scattering parameter in topological insulators due to enhancement in the energy dependence of the relaxation time.

1. Introduction

Currently, high-performance thermoelectrics based on bismuth and antimony chalcogenides with a tetradymite structure [1,2] are attracting much attention as promising 3D topological insulators (TIs). Topological surface states (TSS) in TIs arise due to the inversion of the energy gap edges caused by strong spin–orbit interaction [3,4,5]. In addition, the bulk becomes an insulator, and the electrons on the surface acquire unusual spin momentum locked with linear dispersion and spin polarization specific to the Dirac fermions [6,7]. The linear dispersion and tight interaction between spin and momentum, owing to time-reversal symmetry, prevent fermions from backscattering on non-magnetic defects. This promotes an enhancement in the mobility of charge carriers [8]. Such topological phenomena expand the application possibilities of TIs in various fields of physics [9,10,11,12], including thermoelectricity [4,5,8,13].
An enhancement in thermoelectric performance in TI films based on Bi2Te3 is associated with an increase in the energy dependence of the spectral distribution of the mean free paths of both phonons and electrons [8,13]. Nontrivial energy band structure and large spin–orbit coupling of Dirac fermions in topological insulators account for the enhancement in the energy dependence of the spectral distribution of the mean free path of the electrons lF(E) [4,5]. The estimations of lF(E) in TIs have shown that the electrons possess a wider spectrum than phonons. The effect of electron energy filtering in TIs [14,15] provides an increase in the Seebeck coefficient due to the scattering of charge carriers on the grain interfaces in the films [4,13].
The strong dependence on energy of the mean free path of phonons in the energy range, which is significantly narrower than the electron energy range, accounts for the intensive phonon scattering on an interface between two grains in a polycrystalline material. Furthermore, the enhancement in the phonon scattering results in a decrease in the crystal lattice thermal conductivity in chalcogenide films [8,13,16,17]. In ref. [18], the effect of phonon scattering on grain interfaces on the reduction in the lattice thermal conductivity were considered in solid solutions based on bismuth telluride. Depending on features of the phonon spectrum, the acoustic phonons with low frequency and long wavelength have larger mean free path values than high-frequency ones [18].
The largest heat transfer in the films is determined by long-wavelength phonons that mainly affect the decrease in the lattice thermal conductivity κL due to enhanced scattering. However, at low temperatures up to the Debye temperature (TD = 145 K), the decrease in κL is determined mainly by phonons scattering on acceptor antisite defects and on impurity atoms in solid solutions [19,20]. Furthermore, an additional decrease in κL occurs in the layered films of Bi2Te3 and its solid solutions when phonons are scattered on van der Waals interfaces in the superlattices with a period of about 2 nm, consisting of two inverted quintiples located between the Te(1) layers [21]. With increasing temperature, the effect of scattering on the κL by point defects is reduced while the scattering of phonons on the grain interfaces becomes dominant. However, near room temperature, the contribution of phonon–phonon scattering becomes apparent [16,22,23,24].
A specific feature of topological thermoelectrics is the residual bulk electrical conductivity related to the presence of bulk defects [25,26]. The reduction in bulk conductivity occurs as a result of the optimization of the thermoelectric composition owing to mutual compensation of contributions of acceptor and donor intrinsic defects, providing the increase in TSS contribution to the total conductivity. However, in Bi2Te3-based films, the residual bulk conductivity cannot be completely eliminated.
The systematization of resistivity values in bismuth and antimony chalcogenide films allows one to distinguish, depending on the residual resistivity, two main directions of using TIs [3,5,27,28,29]. Materials with low values of residual resistivity are of high interest for property studies and applications in thermoelectricity [5,28,29]; similarly, materials with high values are significant in the research [27] and development of topological devices [3], as well as for enabling studies of relativistic phenomena, discussed in [30]. A quantitative assessment of these two types of TIs is presented in [27,28,29].
The listed properties of TIs determine a decrease in thermal conductivity and an increase in the thermoelectric power coefficient (Seebeck coefficient). The properties also provide an enhancement in the thermoelectric performance of chalcogenide films compared to bulk thermoelectrics despite a slight reduction in electrical conductivity. In this work, the effect of scattering mechanisms on total κ, lattice κL, and electronic κe thermal conductivities was investigated. This effect is associated with the peculiarities of the energy dependence of the relaxation time in the submicron p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films obtained by discrete and thermal evaporation methods. The selection of formation technique and film composition allows for the optimization of the scattering mechanisms of phonons and electrons, which promote a decrease in the lattice thermal conductivity κL due to the effect of the interlayer surface morphology, grain interfaces, and van der Waals superlattice.

2. Film Deposition Technique and Structure

The polycrystalline films of solid solutions p–(Bi,Sb)2Te3 and p–Bi2Te3 were obtained by the methods of discrete and thermal evaporation in an isothermal chamber which provides a vacuum of 1 × 10−6 torr and homogeneous temperature distribution along the substrate plane. Fresh muscovite mica cleavage planes with a thickness of 3–10 µm and polyimide films with a thickness of 6–20 µm were used as substrates.
To obtain the films by the thermal evaporation method, the starting material was loaded into a quartz crucible heated by a molybdenum coil. The films deposition rate was 10–15 Å/s. During the film formation by discrete evaporation, the initial material in the form of powder with a grain size of about 10 µm was passed by small portions into heated quartz crucible, where it instantly evaporated.
The variation in temperature of the crucible and substrate material of mica (muscovite) and polyimide showed that the optimal temperatures to obtain films with the specified composition are 800–850 °C for the evaporator and 250–300 °C for the substrate. Table 1 represents the technological parameters of the formed p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films for thermal conductivity studies.
To describe the structure of Bi2Te3 and its solid solutions, a primitive rhombohedral or hexagonal unit cell is used. According to X-ray powder diffraction data, the hexagonal unit cell of the Bi2Te3 has a and c parameters of 4.3805 Å and 30.487 Å, respectively [31]. The Te atoms possess either six Bi neighboring atoms or three Bi and Te atoms, allowing us to distinguish the sequence of simple layers and more complex formations consisting of five layers, called quintuples. The considered hexagonal unit cell contains three such quintuples. The Te and Bi atomic layers in the quintuple alternate in the sequence of (-Te(1)-Bi-Te(2)-Bi-Te(1)-), while the chemical bonds in the layers are mainly covalent with some degree of ionicity. It is known that in the p–Bi0.5Sb1.5Te3 solid solutions, the Sb atoms substitute Bi. In these materials, intrinsic antisite point defects on the sites of tellurium BiTe and impurity defects caused by Sb→Bi atom substitutions in solid solutions are revealed. The quintuple boundaries are the interlayer van der Waals surfaces or cleavage planes (0001). The quintuples are bonded by weak van der Waals forces, which provide a slight exfoliation of the crystal along the (0001) planes perpendicular to the crystallographic c axis. The layered structure of the Bi2Te3 crystals and solid solutions determines the significant anisotropy of the transport properties.
In the considered materials, the interlayer surface (0001) possesses the minimum value of free energy [32,33]. It facilitates the formation of stable crystal seeds of bismuth and antimony chalcogenides, which are formed with an orientation predominantly along the c axis perpendicular to the substrate plane, even on substrates with a large mismatch of the crystal lattice parameters. In the process of heat treatment at T = 390 °C in the Ar atmosphere, an intense selective evaporation of the Te atoms from the grain interfaces occurs, in addition to the secondary recrystallization. The depletion of the grain interfaces by tellurium affects the thermoelectric properties and in turn leads to an increase in the charge carrier concentration in the p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films [16]. It should be noted that studies on the morphology of the interlayer surface (0001) of p–Bi2Te3 films using conductive atomic force microscopy [34] showed an increase in local electrical conductivity in the region of the grain interfaces.

3. Atomic Force Microscopy Study of Interlayer Surface (0001) in the Films

The investigation of the interlayer surface (0001) morphology in the p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films was carried out using atomic force microscopy (AFM) in semicontact mode.
Figure 1, Figure 2 and Figure 3 present typical two-dimensional (2D) and three-dimensional (3D) morphology images, and the profiles and the histograms of nanofragment distribution on the surface (0001) depending on the height for all studied samples. As an example, the p–Bi0.5Sb1.5Te3 films were analyzed (samples 1 and 4, Table 1). In these films, the relief of the surface (0001) is composed of separate nanofragments, islands, terraces consisting of coalescent islands, and rows containing dislocations (Figure 1). The observed relief is formed by the diffusion processes and elastic stresses, which result in the deformation of the interlayer surface during the film formation. Separate nanofragments arising on the (0001) surface are scattering centers for phonons, and they lead to a decrease in the thermal conductivity of films based on chalcogenides of bismuth and antimony.
The average heights of nanofragments on the surface (0001) (Figure 1a) are about 36 nm and have similar values for profiles obtained on various surface regions (Figure 2), indicating the homogeneity of the sample surface relief. Further information about the surface relief (0001) was found from the analysis of histograms (Figure 3), which were determined from the morphology images of the p–Bi0.5Sb1.5Te3 films (samples 1 and 4, Table 1, Figure 1).
One can estimate the influence of annealing from the analysis of the distribution of nanofragments on the interlayer surface (0001) depending on their height, considering the p–Bi0.5Sb1.5Te3 films obtained by discrete evaporation as an example (Figure 3). The number of nanofragments with sizes ranging from 16 to 28 nm (Figure 3, curves 1–3) was maximal for p–Bi0.5Sb1.5Te3 film subjected to annealing (Table 1, sample 4). For the unannealed film, the dimensions of most nanofragments were from 31 to 45 nm (Table 1, sample 1, Figure 3, curves 4–6).
The average heights of nanofragments Ra and the root mean square of height deviations of nanofragments Rq (i.e., roughness) vary from 5.5 to 7 nm in the annealed p–Bi0.5Sb1.5Te3 film, while in the unannealed film, Ra and Rq increase from 10 to 15.85 nm. Thus, the influence of annealing results in both a decrease in the number of nanofragments of maximum sizes and a decrease in Ra and Rq values compared with the unannealed p–Bi0.5Sb1.5Te3 film.
Fourier images of the morphology of the film surface (0001) for annealed p–Bi0.5Sb1.5Te3 deposited by the discrete evaporation on a mica substrate were obtained using the Fast Fourier Transform (FFT). These images represent the intensities distribution of two-dimensional reciprocal space being centered at the Γ point of the Brillouin zone (Figure 4a,b). The spectral components of intensities near the Brillouin zone center are associated with the interference of quasiparticle excitations of surface electrons on defects [35,36,37]. In these images, the spectral components are compressed in the vicinity of point Γ (Figure 4). The components of higher orders on Fourier images of the interlayer surface in solid solutions based on Bi2Te3 were found using scanning tunneling microscopy with atomic resolution [38].
The grain parameters were determined from the analysis of the surface morphology images of the p–Bi0.5Sb1.5Te3 films obtained by discrete evaporation (Table 2). It was shown that the film deposited on a mica substrate that was annealed possesses a larger average area of grains than the unannealed ones deposited on a polyimide substrate. In addition, annealing nearly halves the number of grains from 69 to 37 in the studied films, accompanied by an increase in average grain areas, which is in good agreement with [16,17,39].

4. Thermal Conductivity

The thermoelectric properties of the formed films were measured using the Physical Property Measurement System (PPMS) controlled by the Quantum Design Thermal Transport Option (TTO) software. The thermal conductivity of a material was measured by a TTO system as the temperature drops along the sample as a known quantity of heat that passes through the sample. The Seebeck coefficient was obtained by TTO as an electrical voltage drop caused by the temperature gradient across the studied sample. The electrical resistivity measurements were carried out by the PPMS system using the standard four-probe method. The PPMS allows for simultaneous measurements of thermoelectric properties over a wide temperature range from 2 to 390 K with a temperature change rate of ±0.5 K/min. The typical measurement accuracy of the Seebeck coefficient and thermal conductivity is ±5%, and the typical precision of electrical resistivity is 0.01% for 1 Ω and 200 µA.
The total thermal conductivity of the p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films can be represented as: κ = κL + κe, where κL and κe are the crystal lattice and the electronic thermal conductivities, respectively. According to the Wiedemann–Franz law, the relationship between the electronic thermal conductivity κe and electrical conductivity σ is given by κe = L(r, η) σT, where L(r, η) is the Lorentz number:
L ( r , η ) = k e 2 ( r + 7 / 2 ) F r + 5 / 2 ( η ) ( r + 3 / 2 ) F r + 1 / 2 ( η ) ( r + 5 / 2 ) 2 F r + 3 / 2 2 ( η ) ( r + 3 / 2 ) 2 F r + 1 / 2 2 ( η )
where k is the Boltzmann’s constant, e is the charge of an electron, r is the scattering parameter, and η is the reduced Fermi level. Fr + n(η) is the Fermi function at n = 0.5, 1.5, and 2.5 in the form:
F s ( η ) = 0 x s e x η + 1 d x
where s = r + n.
The Lorentz number (1) was calculated for a single parabolic band with arbitrary degeneracy in the isotropic relaxation time approximation as a power law:
τ = τ0Er
The constant τ0 does not depend on energy E, and r is the current value of the effective scattering parameter reff [40]. The parameter reff and the reduced Fermi level η were determined by the least squares method from the experimental values of the thermoelectric power coefficient α(r, η) for the studied films and the degeneracy parameter βd(r, η) in accordance with (2–5). The latter was calculated within the framework of the many-valley model of the energy spectrum from the ratios of isotropic factors in the expressions for electrical conductivity, Hall conductivity, and magnetoconductivity as discussed in [40]. The expressions for the thermoelectric power coefficient α(r, η) and the degeneracy parameter βd(r, η) in the isotropic relaxation time approximation (3):
α = k e ( r + 5 / 2 ) F r + 3 / 2 ( η ) ( r + 3 / 2 ) F r + 1 / 2 ( η ) η
β d r , η = ( 2 r + 3 / 2 ) 2 F 2 r + 1 / 2 2 ( η ) ( r + 3 / 2 ) ( 3 r + 3 / 2 ) F r + 1 / 2 ( η ) F 3 r + 1 / 2 ( η )
The dependences Fr+n(η), reff(η), and α(η) (as shown in Figure 5, Figure 6 and Figure 7, curve 7) illustrate the variation ranges of the reduced Fermi level η and the reff parameter corresponding to the experimental values of the thermopower coefficient α in the temperature range of 40–300 K. The obtained reff parameter (Figure 6) differs from the value of r = −0.5, which is specific to the acoustic scattering mechanism. However, in bulk thermoelectrics, |reff| is typically less than in the films [40].
The increase in |reff| observed in these films determines the enhancement in the relaxation time energy dependence (3) in the films of 3D-TIs [15,41,42]. The strong dependence of relaxation time on energy in TIs is related to the properties of the surface states of the Dirac fermions, which are protected by time-reversal symmetry from backscattering on non-magnetic defects. It promotes the creation of channels that are perfect for rapid, low-dissipative electrical transport and improvement in thermoelectric properties [15]. The mechanism of Dirac fermions scattering and the dependence of τ(E) in TIs have been studied in [41,42] within the framework of non-equilibrium Boltzmann transport theory. Estimations of τ(E) were carried out for cases of long-range Coulomb scattering (LCS) and short-range delta scattering (SDS). It was found that the dependence in the limit gives τ(E) ∝ E for LCS and τ(E) ∝ E−1 for SDS.
In the Dirac semimetal system, such as graphene with Fermi energy (EF) located near the neutral Dirac point, the LCS dominates [41,42,43], and the relaxation time reaches the limit of τ(E) ∝ E. For 3D-TIs based on bismuth and antimony chalcogenides with low thermal conductivity, the limiting dependence of τ(E) ∝ E−1 corresponds to the SDS on neutral impurity or phonon scattering [41,42].
The reff values closest to the limiting τ(E) ∝ E−1 were obtained in an unannealed film of p–Bi0.5Sb1.5Te3 formed by discrete evaporation on a polyimide substrate (Figure 6). Here, the reff varies from −0.925 to −0.77 depending on the reduced Fermi level η (Figure 6, curve 1). According to scanning tunneling spectroscopy data [38], in the films of similar composition with an energy gap of 240 meV, the Fermi energy relative to the conduction band edge is EF = 105 meV, and the Dirac point at ED = −130 meV is located close to the valence band edge [15,41,42]. In this case, the dependence of τ(E) becomes close to the limiting τ(E) ∝ E−1 [15,43]. The more noticeable difference in reff from −1 in the p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films, as shown in Figure 6, is assumed to be caused by the potential fluctuations arising from bulk carriers [15,43].
The Lorentz number L(r, η) (1) employed in the calculation of κe was obtained using the data shown in Figure 5, Figure 6 and Figure 7. The values of the Lorentz number in Figure 7 are consistent with [44] and are significantly lower than L0 = π2k2/3e2 = 2.44 10−8 WΩK−2 for complete degeneracy since the studied films partially degenerate at thermoelectric power α = 242–200 μV K−1 (see Table 1). According to estimations the use of L in calculations for the case of complete degeneracy can lead to errors of up to 40% [45].
Temperature dependences in the total thermal conductivity κ in the films of p–Bi0.5Sb1.5Te3 and p–Bi2Te3 were measured on samples with a thickness of about 1 μm. The effects associated with the surface states of the Dirac fermions can be observed not only in ultra-thin samples, owing to large quantum phase coherence length lφ related to the inelastic electron scattering processes. The value of lφ is usually significantly higher than the electron mean free path lF, which allows the presence of the topological surface states to be revealed by transport property studies in both the nanometer-thick films [7,46,47,48] and also in submicron-thick films [15]. The obtained reff(η) (Figure 6) demonstrates an enhancement in the energy dependence of the relaxation time (3) that confirms the possibility of investigation of TSS in submicron-thick films.
The temperature dependences of the total κ, lattice κL and electron κe thermal conductivities of the p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films shown in Figure 8, Figure 9 and Figure 10 demonstrate that these values depend on the composition, film deposition method, and subsequent thermal treatment. It is known that in TIs, a tight coupling between spin and momentum prevents the backscattering of fermions on nonmagnetic impurities and defects. However, the residual bulk conductivity associated with the presence of bulk intrinsic defects [25,26] remains, which determines the scattering processes in the films under study.
The typical decrease in thermal conductivity in the films with an increase in temperature is consistent with the Peierls theory of three-phonon scattering processes, according to which κ ∝ T−1 at T ≥ TD [49]. The weakening of the thermal conductivity temperature dependencies compared to T−1 is determined by additional scattering processes, which are influenced by the film formation technique, film composition, and thickness.
It should be noted that calculations of the electronic thermal conductivity κe were carried out, taking into account the parameter reff(η) (Figure 6) [40]. A significant decrease in κ and κL, accompanied by the weakening of temperature dependence of κ(T) and κL(T), was observed in the p–Bi0.5Sb1.5Te3 solid solution (Figure 8 and Figure 9, samples 1 and 2), deposited by discrete evaporation on a polyimide substrate. However, the largest decrease in κ and κL in the whole studied temperature range was obtained in the films that were not subjected to heat treatment (Figure 8 and Figure 9, sample 1).
The values of κ and κL were increased and the dependence of κ(T) and κL(T) was enhanced in the low-temperature range in the p–Bi0.5Sb1.5Te3 films, deposited with the discrete evaporation method on the mica substrates followed by annealing in the Ar atmosphere (Figure 8 and Figure 9, sample 3). A further increase in κ and κL and rise in κ(T) and κL(T) were observed in the unannealed p–Bi0.5Sb1.5Te3 films, formed by thermal evaporation on the polyimide substrate (Figure 8 and Figure 9, sample 4). In the unannealed p–Bi2Te3 films deposited on the polyimide substrate obtained by discrete evaporation, the values of κ and κL were significantly higher (Figure 9 and Figure 10, sample 5) than in the p–Bi0.5Sb1.5Te3 solid solution (Figure 8 and Figure 9, sample 1). The largest values of κ and κL in the p–Bi2Te3 films were obtained when thermal evaporation was used (Figure 8 and Figure 9, sample 6). The total κ and lattice κL thermal conductivity values (Figure 8 and Figure 9) were in good agreement with the data for the Bi2Te3 films obtained by mechanical exfoliation [39], for the films deposited onto the polyimide substrates by thermal evaporation [50], and for nanocrystalline composite p–Bi0.52Sb1.48Te3 solid solutions [17].
The considered dependences of κ(T) and κL(T) show that the optimization of technological parameters affects the intensity of the scattering of long-wavelength phonons on the grain interfaces in the films. Thus, the optimization allows for a significant reduction in the values of κ and κL at temperatures above the Debye temperature of TD in the unannealed p–Bi0.5Sb1.5Te3 films, obtained by discrete evaporation on the polyimide substrate [16,22,23,24]. The decrease in κ and κL and the weakening of the κ(T) and κL(T) at temperatures below 100 K is determined by the scattering of phonons on intrinsic antisite defects of BiTe and impurity defects that originated during the formation of the p–Bi0.5Sb1.5Te3 films [19,20]. Additionally, in these films with low values of κ and κL, the thermoelectric power coefficient α increases (Figure 8 and Figure 9, samples 1 and 4, Table 1) due to the effect of energy filtering [14,15].
As shown in Figure 1, Figure 3, Figure 8 and Figure 9, a correlation was established between the thermal conductivity in the p–Bi0.5Sb1.5Te3 films obtained by discrete evaporation and the images of surface morphology (0001). In the unannealed p–Bi0.5Sb1.5Te3 films (Figure 8 and Figure 9, curve 1), the average heights Ra, the root mean square deviations of the heights Rq of nanofragments, and the area and number of grains on the (0001) surface morphology images are increased. Meanwhile, the values of the thermal conductivities of κ and κL are reduced compared to the annealed films (Figure 8 and Figure 9, curves 2 and 4).
Moreover, studies on the total thermal conductivity κ versus thickness were carried out on the p–Bi2Te3 and p–Bi0.5Sb1.5Te3 films obtained by discrete evaporation. It could be considered that at the film thickness of 0.5 µm a tendency towards a decrease in thermal conductivity is revealed in some samples of both the p–Bi2Te3 and the p–Bi0.5Sb1.5Te3 films (Figure 10). A significant decrease in the thermal conductivity κ to 13 and 25 10−3 W cm−1 K−1 was obtained in [21] at 100 K in ultra-thin p–Bi2Te3 films with a thickness of 3.03 and 9.09 nm, respectively.
The dependence of thermal conductivity on the thickness of p–Bi0.5Sb1.5Te3 films (Figure 10) associated with the influence of Dirac fermions surface states is consistent with micro-Raman spectroscopy studies [51]. The appearance of inactive longitudinal phonons A1u2 in the Raman spectra of p–Bi0.5Sb1.5Te3 films becomes noticeable at a film thickness of about 500 nm. The intensity of the A1u2 mode and the intensity ratios of the inactive longitudinal A1u2 and active transverse Eg2 optical phonon modes are found to increase with decreasing thickness [51]. Such a dependence of the intensities of longitudinal and transverse modes on the film thickness is determined by the distortion of the inversion symmetry of the crystal, and characterizes the effects associated with electronic surface states of Dirac fermions [52,53]. An inactive A1u2 mode of weak intensity was also observed in micro-Raman spectra of Bi2Te3 measured at pressures of about 3 GP [54], at which a topological surface phase transition was observed [55].
The properties of the topological surface states of Dirac fermions in micron-thick Sn−(Bi,Sb)2(Te,S)3 films with a tetradymite structure were investigated in thicker films [15,56]. The estimations of the ratios of surface conductivity Gs to the total conductivity Gtot in the films showed that Gs/Gtot = 80% for film thicknesses of 3 µm at temperatures up to 210 K [15] and for 1 µm thick films at T up to 125 K [56], which confirms the domination of the surface transport properties.
In addition to the dependences of the total thermal conductivity κ on thickness in p–Bi2Te3 and p–Bi0.5Sb1.5Te3, similar dependences of resistivity ρ were considered for Bi2Te3 films grown by the molecular-beam epitaxy (MBE) method [25]. Table 3 shows that the values of ρ decreased with increasing thickness.
In the submicron p–Bi2Te3 films studied in this work, the resistivity was significantly lower, and the electrical conductivity was, accordingly, higher [25] (Table 3). This was achieved by optimizing the film formation technological parameters, which resulted in a decrease in bulk point defects accompanied by an increase in electrical conductivity due to an enhancement in the contribution of surface electrical conductivity. Thus, such thickness-dependent analysis of ρ should reasonably be made for films obtained under similar growth conditions.
Temperature dependences of the electronic thermal conductivity κe are determined by the corresponding electrical conductivity σ(T) (Figure 11, inset), which includes the contribution of both the metallic surface and the bulk electrical conductivity of TIs films [25,26,27]. The variation in the technological parameters of the film formation allows us to reduce the bulk conductivity by optimizing the scattering of electrons on intrinsic antisite defects of BiTe and impurity substitution defects of Sb → Bi atoms in the p–Bi0.5Sb1.5Te3 solid solution. In the low-temperature range, an enhancement in κe with an increase in the temperature from 40 to 80 K was observed (Figure 11, curves 1–3, inset), and a significant decrease in the angular coefficients of dlnσ/dlnT was obtained compared to the interval of 100–300 K in the p–Bi0.5Sb1.5Te3 films. The revealed feature of κe(T) at low temperatures is explained by the enhancement in the energy dependence of relaxation time τ due to the rise in |reff|, which was used in the Lorentz number (1) calculations.
The increase in |reff| with decreasing Fermi level η (Figure 6) is associated with the strengthening of electron scattering on intrinsic and impurity substitution defects at low temperatures [19,20]. Intensive scattering on defects with a temperature decrease from 80 to 40 K determines the reduction in electrical conductivity σ and corresponding decrease in electronic thermal conductivity κe. The “hump” on the κe(T) dependence observed in the temperature range of 80–110 K can be explained by potential fluctuations due to the influence of bulk carriers. These fluctuations result in the deviation of the energy dependence of relaxation time from the limiting τ(E) ∝ E−1 [15,41,42], and therefore to the decrease in the |reff| parameter with rise in η (Figure 6).
The temperature dependences σ(T) (Figure 11, inset) were analyzed using the ρ = 1/σ, obtained as resistivity ratios between values of ρ measured at room (300 K) and low (20 K) temperatures in accordance with [29,57]. The observed decrease in the slopes σ(T) at low temperatures down to 80 K (Figure 11, inset) turned out to be insufficient to change the σ(T) of the films from a metallic conductivity type to a semiconductor one, in contrast to [27], where the ratios of the ρ(10 K)/ρ(300 K) were considered.
Our studies have shown that the ratios of ρ(300 K)/ρ(20 K) (Figure 11) as well as the total κ and the lattice κL thermal conductivities (Figure 8 and Figure 9) obtained in the p–Bi0.5Sb1.5Te3 are lower than in the p–Bi2Te3 regardless of the film formation technique used. A decrease in the ratios ρ(300 K)/ρ(20 K) in the p–Bi0.5Sb1.5Te3 films obtained by discrete evaporation is accompanied by a decrease in κ and κL (Figure 8 and Figure 9, samples 1, 2, and 4).
The carried-out estimations of the thermoelectric figure of merit ZT, where Z = α2σ/κ, showed that in the p–Bi0.5Sb1.5Te3 films formed by discrete evaporation, the highest ZT values near room temperature were 1, 1.12, and 1.08 in samples 1, 2, and 4 (Table 1), respectively. In the p–Bi0.5Sb1.5Te3 films (samples 1 and 2) deposited on a polyimide substrate, the ratios of ρ(300 K)/ρ(20 K) were lower than in sample 4 obtained on a mica substrate, which possessed high electrical conductivity (Figure 11, curve 4, inset). In the p–Bi2Te3 films (samples 5 and 6, Table 1) with high ratios ρ(300 K)/ρ(20 K) of about 10, the ZT values decreased from 0.87 upon discrete evaporation to 0.65 upon thermal evaporation. In the bulk samples of the p–Bi0.5Sb1.5Te3 solid solution (Figure 11, curves 7 and 8, inset), the ratio is about 8. Thus, the resistivity ratios are ρ(300 K)/ρ(20 K) and can be used to analyze the thermoelectric properties of topological insulator films based on p–Bi2Te3.
It is well known from recent studies that the values of the thermoelectric figure of merit in chalcogenide films of TIs can vary over a wide range depending on the formation technique and composition. Therefore, near room temperature, the values of ZT were obtained as 0.9 in the p–Bi2Te3 [50] and the p–Bi0.5Sb1.5Te3 [57] films, 0.37 in p–Bi0.5Sb1.5Te3 films [58], and 1.44 in the p–type Bi2Te3/Bi0.5Sb1.5Te3 superlattice films with a low interfacial resistance [59].
More detailed information about the properties of the surface states of Dirac fermions was obtained from studies on quantum oscillations of the magnetoresistance in strong magnetic fields [60,61,62], angle-resolved photoemission spectroscopy [63,64], and differential tunneling conductance spectra using scanning tunneling spectroscopy [35,38,65,66]. The existence of the surface states of Dirac fermions in the p–Bi0.5Sb1.5Te3 and the p–Bi2Te3 submicron films was confirmed from studies of differential tunneling conductance spectra by scanning tunneling spectroscopy at room temperature [38]. In the studied p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films, the values of the thermoelectric power coefficient are α = 178 μV K−1 and 272 μV K−1, and the surface concentration of fermions ns is 4 1012 cm−2 and 1.2 1013 cm−2, respectively. The quantum Hall effect, the parameters of the surface states of Dirac fermions, and the metallic conductivity on the surface of the p–Bi0.5Sb1.5Te3 and the p–Bi2Te3 films, obtained from studies of galvanomagnetic properties in strong magnetic fields, confirm the topological properties of submicron-thick films [15,56] and the prospects for practical application.
Currently, high-pressure methods are employed for the optimization of thermoelectric properties. In the developed p–Bi0.5Sb1.5Te3 films, considered in this study, topological phase transitions (TPTs) were found at pressures of P = 3–4 GPa as a result of pressure-dependent thermoelectric property investigations [34]. Near the TPTs at room temperature, the power factor increases by more than two times, while the thermal conductivity rises slightly [67], which makes the considered films promising as the p-branch for the development of thermoelectric modules [68]. The module consists of a specially designed high-pressure cell compressed between two opposing anvils.
Multicomponent solid solutions Bi2(Te, Se, S)3 were used as the n-branch of thermoelements. The performance parameters of the high-pressure module can be optimized by changing the pressure applied to its elements [68]. Additionally, the thermoelectric parameters can be enhanced due to the appearance of additional transverse and longitudinal magneto-thermoelectric effects in an external magnetic field [69].

5. Conclusions

Investigations of the thermal conductivity of the p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films of TIs formed by discrete and thermal evaporation techniques were carried out. The study revealed that the total κ, lattice κL, and electronic κe thermal conductivities depend not only on the composition and optimization of the film formation parameters but also on the influence of the energy dependence of the relaxation time. Additionally, the effective scattering parameter of reff was taken into account while calculating the κe values using the Lorentz number L(r, η). The unannealed films of the p–Bi0.5Sb1.5Te3 solid solution formed by discrete evaporation on polyimide substrates exhibited a significant decrease in κ, κL, and κe values.
The reduction in the crystal lattice thermal conductivity κL in TI films is associated with the enhancement in the energy dependence of the mean free path of phonons, which leads to the intensive scattering of long-wavelength phonons on the grain interfaces. At temperatures near and above TD up to room temperature, the contribution of phonon–phonon scattering increases. In the range of low temperatures, the main reason for the decrease in κL is the scattering of short-wavelength phonons on intrinsic point defects and impurities. Furthermore, an additional decrease in κL also occurs due to the scattering of phonons on the interfaces of the Te(1) layers of the van der Waals gap.
The reduction in κe with temperature is explained by a decrease in the electrical conductivity σ in the films caused by electron scattering on the intrinsic antisite BiTe and impurity substitution defects of Sb → Bi. The peculiarity of the dependence κe(T) in the temperature range of 80–110 K in the form of a hump is explained by the influence of potential fluctuations in bulk carriers. These fluctuations lead to a decrease in the effective scattering parameter |reff| with increasing Fermi level due to the deviation of the energy dependence of relaxation time from the limiting τ(E) ∝ E−1.
The study also investigated the effect of heat treatment on the relief properties of the interlayer surface (0001) studied by AFM and on the variation in the thermal conductivity in the p–Bi0.5Sb1.5Te3 films. The results show that the total κ and the crystal lattice κL thermal conductivities are reduced, but the number of grains on the surface (0001) and the nanofragment parameters (the average heights Ra and the root mean square deviations in the heights Rq) are increased in the unannealed films compared to the annealed ones.
Finally, the low thermal conductivity and high electrical conductivity obtained in annealed p–Bi0.5Sb1.5Te3 films formed by discrete evaporation on a polyimide substrate provide an increase in the thermoelectric figure of merit. The observed increase in ZT up to 1.12 near room temperature is associated with an enhancement in the energy dependence of the relaxation time in TIs caused by an increase in the effective scattering parameter reff.

Author Contributions

Conceptualization, L.N.L.; Data curation, L.N.L., O.A.U. and Y.A.B.; Formal analysis, L.N.L. and O.A.U.; Investigation, Y.A.B. and I.V.M.; Methodology, V.A.D., I.V.M. and V.N.P.; Resources, L.N.L., Y.A.B., V.A.D. and I.V.M.; Software, O.A.U.; Writing—original draft, L.N.L. and O.A.U.; Writing—review & editing, L.N.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Russian Foundation for Basic Research grant number 20-08-00464.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nolas, G.S.; Sharp, J.; Goldsmid, H.J. Thermoelectrics: Basic Principles and New Materials Developments; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2001. [Google Scholar]
  2. Rowe, D.M. Modules, Systems, and Applications in Thermoelectrics, 1st ed.; CRC Press: Boca Raton, FL, USA; London, UK; New York, NY, USA, 2012. [Google Scholar]
  3. Gilbert, M.J. Topological electronics. Commun. Phys. 2021, 4, 70. [Google Scholar] [CrossRef]
  4. Heremans, J.P.; Cava, R.J.; Samarth, N. Tetradymites as thermoelectrics and topological insulators. Nat. Rev. Mater. 2017, 2, 17049. [Google Scholar] [CrossRef]
  5. Ngabonziza, P. Quantum transport and potential of topological states for thermoelectricity in Bi2Te3 thin films. Nanotechnology 2022, 33, 192001. [Google Scholar] [CrossRef] [PubMed]
  6. Hasan, M.Z.; Kane, C.L. Colloquium: Topological insulators. Rev. Mod. Phys. 2010, 82, 3045–3067. [Google Scholar] [CrossRef]
  7. Chen, Y.L.; Analytis, J.G.; Chu, J.-H.; Liu, Z.K.; Mo, S.-K.; Qi, X.L.; Zhang, H.J.; Lu, D.H.; Dai, X.; Fang, Z.; et al. Experimental Realization of a Three-Dimensional Topological Insulator, Bi2 Te3. Science 2009, 325, 178–181. [Google Scholar] [CrossRef] [PubMed]
  8. Xu, N.; Xu, Y.; Zhu, J. Topological insulators for thermoelectrics. npj Quantum Mater. 2017, 2, 51. [Google Scholar] [CrossRef]
  9. Lee, J.; Koo, J.; Jhon, Y.M.; Lee, J.H. A femtosecond pulse erbium fiber laser incorporating a saturable absorber based on bulk-structured Bi2Te3 topological insulator. Opt. Express 2014, 22, 6165–6173. [Google Scholar] [CrossRef]
  10. Liu, H.; Zhu, X.; Sun, X.; Zhu, C.; Huang, W.; Zhang, X.; Zheng, B.; Zou, Z.; Luo, Z.; Wang, X.; et al. Self-Powered Broad-band Photodetectors Based on Vertically Stacked WSe2/Bi2Te3 p–n Heterojunctions. ACS Nano 2019, 13, 13573–13580. [Google Scholar] [CrossRef]
  11. Sun, R.; Yang, S.; Yang, X.; Vetter, E.; Sun, D.; Li, N.; Su, L.; Li, Y.; Gong, Z.-Z.; Xie, Z.-K.; et al. Large Tunable Spin-to-Charge Conversion Induced by Hybrid Rashba and Dirac Surface States in Topological Insulator Heterostructures. Nano Lett. 2019, 19, 4420–4426. [Google Scholar] [CrossRef]
  12. Jiang, Z.; Chang, C.-Z.; Masir, M.R.; Tang, C.; Xu, Y.; Moodera, J.S.; MacDonald, A.H.; Shi, J. Enhanced spin Seebeck effect signal due to spin-momentum locked topological surface states. Nat. Commun. 2016, 7, 11458. [Google Scholar] [CrossRef]
  13. Liu, T.-H.; Zhou, J.; Li, M.; Ding, Z.; Song, Q.; Liao, B.; Fu, L.; Chen, G. Electron mean-free-path filtering in Dirac material for improved thermoelectric performance. Proc. Natl. Acad. Sci. USA 2018, 115, 879–884. [Google Scholar] [CrossRef] [PubMed]
  14. Zebarjadi, M.; Esfarjani, K.; Dresselhaus, M.S.; Ren, Z.F.; Chen, G. Perspectives on thermoelectrics: From fundamentals to device applications. Energy Environ. Sci. 2012, 5, 5147–5162. [Google Scholar] [CrossRef]
  15. Matsushita, S.Y.; Ichimura, K.; Huynh, K.K.; Tanigaki, K. Large thermopower in topological surface state of Sn-BSTS topological insulators: Thermoelectrics and energy-dependent relaxation times. Phys. Rev. Mater. 2021, 5, 014205. [Google Scholar] [CrossRef]
  16. Puneet, P.; Podila, R.; Karakaya, M.; Zhu, S.; He, J.; Tritt, T.M.; Dresselhaus, M.S.; Rao, A.M. Preferential Scattering by Interfacial Charged Defects for Enhanced Thermoelectric Performance in Few-layered n-type Bi2Te3. Sci. Rep. 2013, 3, 3212. [Google Scholar] [CrossRef] [PubMed]
  17. Szczech, J.R.; Higgins, J.M.; Jin, S. Enhancement of the thermoelectric properties in nanoscale and nanostructured materials. J. Mater. Chem. 2011, 21, 4037–4055. [Google Scholar] [CrossRef]
  18. Zheng, J.; Kodera, Y.; Xu, X.; Shin, S.; Chung, K.M.; Imai, T.; Ihnfeldt, R.V.; Garay, J.E.; Chen, R. Suppressing thermal conductivity of nano-grained thermoelectric material using acoustically hard nanoparticles. J. Appl. Phys. 2021, 130, 235106. [Google Scholar] [CrossRef]
  19. Nam, H.; Xu, Y.; Miotkowski, I.; Tian, J.; Chen, Y.P.; Liu, C.; Hasan, M.Z.; Zhu, W.; Fiete, G.A.; Shih, C.-K. Microscopic investigation of Bi2−xSbxTe3−ySey systems: On the origin of a robust intrinsic topological insulator. J. Phys. Chem. Solids 2019, 128, 251–257. [Google Scholar] [CrossRef]
  20. Bessas, D.; Sergueev, I.; Wille, H.-C.; Perßon, J.; Ebling, D.; Hermann, R.P. Lattice dynamics in Bi2Te3 and Sb2Te3: Te and Sb density of phonon states. Phys. Rev. B 2012, 86, 224301. [Google Scholar] [CrossRef]
  21. Park, K.H.; Mohamed, M.; Aksamija, Z.; Ravaioli, U. Phonon scattering due to van der Waals forces in the lattice thermal conductivity of Bi2Te3 thin films. J. Appl. Phys. 2015, 117, 015103. [Google Scholar] [CrossRef]
  22. Poudel, B.; Hao, Q.; Ma, Y.; Lan, Y.; Minnich, A.; Yu, B.; Yan, X.; Wang, D.; Muto, A.; Vashaee, D.; et al. High-Thermoelectric Performance of Nanostructured Bismuth Antimony Telluride Bulk Alloys. Science 2008, 320, 634–638. [Google Scholar] [CrossRef]
  23. Zhu, T.; Hu, L.; Zhao, X.; He, J. New Insights into Intrinsic Point Defects in V2VI3 Thermoelectric Materials. Adv. Sci. 2016, 3, 1600004. [Google Scholar] [CrossRef] [PubMed]
  24. Wei, P.; Yang, J.; Guo, L.; Wang, S.; Wu, L.; Xu, X.; Zhao, W.; Zhang, Q.; Zhang, W.; Dresselhaus, M.S.; et al. Minimum Thermal Conductivity in Weak Topological Insulators with Bismuth-Based Stack Structure. Adv. Funct. Mater. 2016, 26, 5360–5367. [Google Scholar] [CrossRef]
  25. Zhang, J.; Chang, C.-Z.; Zhang, Z.; Wen, J.; Feng, X.; Li, K.; Liu, M.; He, K.; Wang, L.; Chen, X.; et al. Band structure engineering in (Bi1−xSbx)2Te3 ternary topological insulators. Nat. Commun. 2011, 2, 574. [Google Scholar] [CrossRef] [PubMed]
  26. Knispel, T.; Jolie, W.; Borgwardt, N.; Lux, J.; Wang, Z.; Ando, Y.; Rosch, A.; Michely, T.; Grüninger, M. Charge puddles in the bulk and on the surface of the topological insulator BiSbTeSe2 studied by scanning tunneling microscopy and optical spectroscopy. Phys. Rev. B 2017, 96, 195135. [Google Scholar] [CrossRef]
  27. Ren, Z.; Taskin, A.A.; Sasaki, S.; Segawa, K.; Ando, Y. Optimizing Bi2−xSbxTe3−ySeysolid solutions to approach the intrinsic topological insulator regime. Phys. Rev. B 2011, 84, 165311. [Google Scholar] [CrossRef]
  28. Ramachandran, B.; Wu, K.K.; Kuo, Y.K.; Guo, L.S.; Wang, L.M. Compositional effects on the low-temperature transport properties of non-stoichiometric Bi2TexSey-based crystals. J. Phys. D Appl. Phys. 2015, 50, 025302. [Google Scholar] [CrossRef]
  29. Chiatti, O.; Riha, C.; Lawrenz, D.; Busch, M.; Dusari, S.; Sánchez-Barriga, J.; Mogilatenko, A.; Yashina, L.V.; Valencia, S.; Ünal, A.A.; et al. 2D layered transport properties from topological insulator Bi2Se3 single crystals and micro flakes. Sci. Rep. 2016, 6, 27483. [Google Scholar] [CrossRef]
  30. Yar, A.; Naeem, M.; Khan, S.U.; Sabeeh, K. Hybridization effects on wave packet dynamics in topological insulator thin films. J. Phys. Condens. Matter 2017, 29, 465002. [Google Scholar] [CrossRef]
  31. Francombe, M.H. Structure-cell data and expansion coefficients of bismuth telluride. Br. J. Appl. Phys. 1958, 9, 415–417. [Google Scholar] [CrossRef]
  32. Medlin, D.L.; Ramasse, Q.M.; Spataru, C.D.; Yang, N.Y.C. Structure of the (0001) basal twin boundary in Bi2Te3. J. Appl. Phys. 2010, 108, 043517. [Google Scholar] [CrossRef]
  33. Virk, N.; Yazyev, O.V. Dirac fermions at high-index surfaces of bismuth chalcogenide topological insulator nanostructures. Sci. Rep. 2016, 6, 20220. [Google Scholar] [CrossRef] [PubMed]
  34. Lukyanova, L.N.; Usov, O.A.; Volkov, M.P.; Makarenko, I.V. Topological Thermoelectric Materials Based on Bismuth Telluride. Nanobiotechnol. Rep. 2021, 16, 282–293. [Google Scholar] [CrossRef]
  35. He, X.; Li, H.; Chen, L.; Wu, K. Substitution-induced spin-splitted surface states in topological insulator (Bi1−xSbx)2Te3. Sci. Rep. 2015, 5, 8830. [Google Scholar] [CrossRef] [PubMed]
  36. Ko, W.; Jeon, I.; Kim, H.W.; Kwon, H.; Kahng, S.-J.; Park, J.; Kim, J.S.; Hwang, S.W.; Suh, H. Atomic and electronic structure of an alloyed topological insulator, Bi1.5Sb0.5Te1.7Se1.3. Sci. Rep. 2013, 3, 2656. [Google Scholar] [CrossRef] [PubMed]
  37. Sessi, P.; Otrokov, M.M.; Bathon, T.; Vergniory, M.G.; Tsirkin, S.S.; Kokh, K.A.; Tereshchenko, O.E.; Chulkov, E.V.; Bode, M. Visualizing spin-dependent bulk scattering and breakdown of the linear dispersion relation in Bi2Te3. Phys. Rev. B 2013, 88, 161407. [Google Scholar] [CrossRef]
  38. Lukyanova, L.N.; Makarenko, I.V.; Usov, O.A.; Dementev, P.A. Scanning tunneling spectroscopy of the surface states of Dirac fermions in thermoelectrics based on bismuth telluride. Semicond. Sci. Technol. 2018, 33, 055001. [Google Scholar] [CrossRef]
  39. Goyal, V.; Teweldebrhan, D.; Balandin, A.A. Mechanically-exfoliated stacks of thin films of Bi2Te3 topological insulators with enhanced thermoelectric performance. Appl. Phys. Lett. 2010, 97, 133117. [Google Scholar] [CrossRef]
  40. Lukyanova, L.N.; Boikov, Y.A.; Danilov, V.A.; Usov, O.A.; Volkov, M.P.; Kutasov, V.A. Thermoelectric and galvanomagnetic properties of bismuth chalcogenide nanostructured heteroepitaxial films. Semicond. Sci. Technol. 2015, 30, 015011. [Google Scholar] [CrossRef]
  41. Nomura, K.; MacDonald, A.H. Quantum Transport of Massless Dirac Fermions. Phys. Rev. Lett. 2007, 98, 076602. [Google Scholar] [CrossRef]
  42. Nomura, K.; Koshino, M.; Ryu, S. Topological Delocalization of Two-Dimensional Massless Dirac Fermions. Phys. Rev. Lett. 2007, 99, 146806. [Google Scholar] [CrossRef]
  43. Chiba, T.; Takahashi, S. Transport properties on an ionically disordered surface of topological insulators: Toward high-performance thermoelectrics. J. Appl. Phys. 2019, 126, 245704. [Google Scholar] [CrossRef]
  44. Kim, H.-S.; Gibbs, Z.M.; Tang, Y.; Wang, H.; Snyder, G.J. Characterization of Lorenz number with Seebeck coefficient measurement. APL Mater. 2015, 3, 041506. [Google Scholar] [CrossRef]
  45. Toberer, E.S.; Baranowski, L.L.; Dames, C. Advances in Thermal Conductivity. Annu. Rev. Mater. Res. 2012, 42, 179–209. [Google Scholar] [CrossRef]
  46. Cao, H.; Tian, J.; Miotkowski, I.; Shen, T.; Hu, J.; Qiao, S.; Chen, Y.P. Quantized Hall Effect and Shubnikov–de Haas Oscillations in Highly Doped Bi2Se3: Evidence for Layered Transport of Bulk Carriers. Phys. Rev. Lett. 2012, 108, 216803. [Google Scholar] [CrossRef] [PubMed]
  47. Shahil, K.M.F.; Hossain, M.Z.; Teweldebrhan, D.; Balandin, A.A. Crystal symmetry breaking in few-quintuple Bi2Te3 films: Applications in nanometrology of topological insulators. Appl. Phys. Lett. 2010, 96, 153103. [Google Scholar] [CrossRef]
  48. Chen, Y.L.; Chu, J.-H.; Analytis, J.G.; Liu, Z.K.; Igarashi, K.; Kuo, H.-H.; Qi, X.L.; Mo, S.K.; Moore, R.G.; Lu, D.H.; et al. Massive Dirac Fermion on the Surface of a Magnetically Doped Topological Insulator. Science 2010, 329, 659–662. [Google Scholar] [CrossRef]
  49. Klemens, P.G. Thermal conductivity and lattice vibrational modes. Solid State Phys. 1958, 7, 1–98. [Google Scholar] [CrossRef]
  50. Goncalves, L.; Couto, C.; Alpuim, P.; Rolo, A.; Völklein, F.; Correia, J. Optimization of thermoelectric properties on Bi2Te3 thin films deposited by thermal co-evaporation. Thin Solid Films 2010, 518, 2816–2821. [Google Scholar] [CrossRef]
  51. Luk’yanova, L.N.; Bibik, A.Y.; Aseev, V.A.; Usov, O.A.; Makarenko, I.V.; Petrov, V.N.; Nikonorov, N.V.; Kutasov, V.A. Surface morphology and Raman spectroscopy of thin layers of antimony and bismuth chalcogenides. Phys. Solid State 2016, 58, 1440–1447. [Google Scholar] [CrossRef]
  52. Teweldebrhan, D.; Goyal, V.; Balandin, A.A. Exfoliation and Characterization of Bismuth Telluride Atomic Quintuples and Quasi-Two-Dimensional Crystals. Nano Lett. 2010, 10, 1209–1218. [Google Scholar] [CrossRef]
  53. Plucinski, L.; Herdt, A.; Fahrendorf, S.; Bihlmayer, G.; Mussler, G.; Döring, S.; Kampmeier, J.; Matthes, F.; Bürgler, D.E.; Grützmacher, D.; et al. Electronic structure, surface morphology, and topologically protected surface states of Sb2Te3 thin films grown on Si(111). J. Appl. Phys. 2013, 113, 053706. [Google Scholar] [CrossRef]
  54. Ovsyannikov, S.V.; Morozova, N.V.; Korobeinikov, I.V.; Lukyanova, L.N.; Manakov, A.Y.; Likhacheva, A.Y.; Ancharov, A.I.; Vokhmyanin, A.P.; Berger, I.F.; Usov, O.A.; et al. Enhanced power factor and high-pressure effects in (Bi,Sb)2(Te,Se)3 thermoelectrics. Appl. Phys. Lett. 2015, 106, 143901. [Google Scholar] [CrossRef]
  55. Korobeinikov, I.V.; Morozova, N.V.; Lukyanova, L.N.; Usov, O.A.; Ovsyannikov, S.V. On the Power Factor of Bismuth-Telluride-Based Alloys near Topological Phase Transitions at High Pressures. Semiconductors 2019, 53, 732–736. [Google Scholar] [CrossRef]
  56. Xu, Y.; Miotkowski, I.; Liu, C.; Tian, J.; Nam, H.; Alidoust, N.; Hu, J.; Shih, C.-K.; Hasan, M.Z.; Chen, Y.P. Observation of topological surface state quantum Hall effect in an intrinsic three-dimensional topological insulator. Nat. Phys. 2014, 10, 956–963. [Google Scholar] [CrossRef]
  57. Parashchuk, T.; Kostyuk, O.; Nykyruy, L.; Dashevsky, Z. High thermoelectric performance of p-type Bi0.5Sb1.5Te3 films on flexible substrate. Mater. Chem. Phys. 2020, 253, 123427. [Google Scholar] [CrossRef]
  58. Liao, C.-N.; Su, X.-W.; Liou, K.-M.; Chu, H.-S. Electrical and thermal transport properties of electrically stressed Bi–Sb–Te nanocrystalline thin films. Thin Solid Films 2011, 519, 4394–4399. [Google Scholar] [CrossRef]
  59. Park, N.-W.; Lee, W.-Y.; Yoon, Y.-S.; Kim, G.-S.; Yoon, Y.-G.; Lee, S.-K. Achieving Out-of-Plane Thermoelectric Figure of Merit ZT = 1.44 in a p-Type Bi2Te3/Bi0.5Sb1.5Te3 Superlattice Film with Low Interfacial Resistance. ACS Appl. Mater. Interfaces 2019, 11, 38247–38254. [Google Scholar] [CrossRef]
  60. Liu, H.; Liu, S.; Yi, Y.; He, H.; Wang, J. Shubnikov–de Haas oscillations in n and p type Bi2Se3 flakes. 2D Mater. 2015, 2, 045002. [Google Scholar] [CrossRef]
  61. Bao, L.; He, L.; Meyer, N.; Kou, X.; Zhang, P.; Chen, Z.-G.; Fedorov, A.V.; Zou, J.; Riedemann, T.M.; Lograsso, T.A.; et al. Weak Anti-localization and Quantum Oscillations of Surface States in Topological Insulator Bi2Se2Te. Sci. Rep. 2012, 2, 726. [Google Scholar] [CrossRef]
  62. Tu, N.H.; Tanabe, Y.; Satake, Y.; Huynh, K.K.; Le, P.H.; Matsushita, S.Y.; Tanigaki, K. Large-area and transferred high-quality three-dimensional topological insulator Bi2−xSbxTe3−ySey ultrathin film by catalyst-free physical vapor deposition. Nano Lett. 2017, 17, 2354–2360. [Google Scholar] [CrossRef]
  63. Sumida, K.; Ishida, Y.; Zhu, S.; Ye, M.; Pertsova, A.; Triola, C.; Kokh, K.A.; Tereshchenko, O.E.; Balatsky, A.V.; Shin, S.; et al. Prolonged duration of nonequilibrated Dirac fermions in neutral topological insulators. Sci. Rep. 2017, 7, 14080. [Google Scholar] [CrossRef] [PubMed]
  64. Alpichshev, Z.; Analytis, J.G.; Chu, J.-H.; Fisher, I.R.; Chen, Y.L.; Shen, Z.X.; Fang, A.; Kapitulnik, A. STM Imaging of Electronic Waves on the Surface of Bi2Te3: Topologically Protected Surface States and Hexagonal Warping Effects. Phys. Rev. Lett. 2010, 104, 016401. [Google Scholar] [CrossRef] [PubMed]
  65. Beidenkopf, H.; Roushan, P.; Seo, J.; Gorman, L.; Drozdov, I.; Hor, Y.S.; Cava, R.J.; Yazdani, A. Spatial fluctuations of helical Dirac fermions on the surface of topological insulators. Nat. Phys. 2011, 7, 939–943. [Google Scholar] [CrossRef]
  66. Ko, W.; Park, J.; Jeon, I.; Kim, H.W.; Kwon, H.; Oh, Y.; Kim, J.S.; Suh, H.; Hwang, S.W.; Chung, C. Local potential fluctuation of topological surface states in Bi1.5Sb0.5Te1.7Se1.3 observed by Landau level spectroscopy. Appl. Phys. Lett. 2016, 108, 083109. [Google Scholar] [CrossRef]
  67. Jacobsen, M.; Sinogeikin, S.; Kumar, R.; Cornelius, A. High pressure transport characteristics of Bi2Te3, Sb2Te3, and BiSbTe3. J. Phys. Chem. Solids 2012, 73, 1154–1158. [Google Scholar] [CrossRef]
  68. Korobeinikov, I.V.; Morozova, N.V.; Lukyanova, L.N.; Usov, O.A.; Kulbachinskii, V.A.; Shchennikov, V.V.; Ovsyannikov, S.V. Stress-controlled thermoelectric module for energy harvesting and its application for the significant enhancement of the power factor of Bi2Te3-based thermoelectrics. J. Phys. D Appl. Phys. 2018, 51, 025501. [Google Scholar] [CrossRef]
  69. Morozova, N.V.; Korobeinikov, I.V.; Ovsyannikov, S.V. Strategies and challenges of high-pressure methods applied to thermoelectric materials. J. Appl. Phys. 2019, 125, 220901. [Google Scholar] [CrossRef]
Figure 1. The (a,b) 2D and (c,d) 3D images of surface morphology (0001) of the unannealed p–Bi0.5Sb1.5Te3 film obtained by discrete evaporation on a polyimide substrate (a,c) and of the annealed one deposited on a mica substrate (b,d).
Figure 1. The (a,b) 2D and (c,d) 3D images of surface morphology (0001) of the unannealed p–Bi0.5Sb1.5Te3 film obtained by discrete evaporation on a polyimide substrate (a,c) and of the annealed one deposited on a mica substrate (b,d).
Magnetochemistry 09 00141 g001
Figure 2. The profiles (1–4) of the surface morphology (0001) images of the p–Bi0.5Sb1.5Te3 films (sample 1, Table 1) obtained along arbitrary and mutually perpendicular directions in Figure 1a: 1, 2 and 3, 4. Average height of nanofragments on the surface (0001) in nm: 1–36.35; 2–36.24; 3–35.54; 4–36.49.
Figure 2. The profiles (1–4) of the surface morphology (0001) images of the p–Bi0.5Sb1.5Te3 films (sample 1, Table 1) obtained along arbitrary and mutually perpendicular directions in Figure 1a: 1, 2 and 3, 4. Average height of nanofragments on the surface (0001) in nm: 1–36.35; 2–36.24; 3–35.54; 4–36.49.
Magnetochemistry 09 00141 g002
Figure 3. Height distribution of the nanofragments on the interlayer surface (0001) in the p–Bi0.5Sb1.5Te3 films depending on the heights (1–3) in the annealed film obtained by discrete evaporation on a mica substrate and (4–6) the unannealed film deposited on a polyimide substrate.
Figure 3. Height distribution of the nanofragments on the interlayer surface (0001) in the p–Bi0.5Sb1.5Te3 films depending on the heights (1–3) in the annealed film obtained by discrete evaporation on a mica substrate and (4–6) the unannealed film deposited on a polyimide substrate.
Magnetochemistry 09 00141 g003
Figure 4. Fourier images of the surface morphology (0001) of the annealed p–Bi0.5Sb1.5Te3 film obtained by the discrete evaporation method on a mica substrate. (a) Two-dimensional and (b) three-dimensional Fourier images, respectively.
Figure 4. Fourier images of the surface morphology (0001) of the annealed p–Bi0.5Sb1.5Te3 film obtained by the discrete evaporation method on a mica substrate. (a) Two-dimensional and (b) three-dimensional Fourier images, respectively.
Magnetochemistry 09 00141 g004
Figure 5. Fermi functions Fr+n(η) for the p–Bi0.5Sb1.5Te3 (1−4) and p–Bi2Te3 (5, 6) films. Sample numbers in figure and subsequent figures correspond to Table 1.
Figure 5. Fermi functions Fr+n(η) for the p–Bi0.5Sb1.5Te3 (1−4) and p–Bi2Te3 (5, 6) films. Sample numbers in figure and subsequent figures correspond to Table 1.
Magnetochemistry 09 00141 g005
Figure 6. The dependences of the effective scattering parameter of the charge carriers (reff) on the reduced Fermi level η in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Figure 6. The dependences of the effective scattering parameter of the charge carriers (reff) on the reduced Fermi level η in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Magnetochemistry 09 00141 g006
Figure 7. The dependences in the thermoelectric power coefficient α and the Lorentz number L on reduced Fermi level η in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Figure 7. The dependences in the thermoelectric power coefficient α and the Lorentz number L on reduced Fermi level η in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Magnetochemistry 09 00141 g007
Figure 8. The temperature dependences of the total thermal conductivity κ in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Figure 8. The temperature dependences of the total thermal conductivity κ in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Magnetochemistry 09 00141 g008
Figure 9. Temperature dependencies in the lattice thermal conductivity κL in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Figure 9. Temperature dependencies in the lattice thermal conductivity κL in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films.
Magnetochemistry 09 00141 g009
Figure 10. Dependences of the total thermal conductivity κ on thickness in p–Bi2Te3 (1) and p–Bi0.5Sb1.5Te3 (2) films obtained by discrete evaporation at T = 77 K.
Figure 10. Dependences of the total thermal conductivity κ on thickness in p–Bi2Te3 (1) and p–Bi0.5Sb1.5Te3 (2) films obtained by discrete evaporation at T = 77 K.
Magnetochemistry 09 00141 g010
Figure 11. The temperature dependencies of the electronic thermal conductivity κe and the electrical conductivity σ (inset) in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films. The designations of the samples (1–6) correspond to those in the Table 1, and 7 and 8 are the bulk samples of the p–Bi0.5Sb1.5Te3 solid solution. In the films of p–Bi0.5Sb1.5Te3 (1–3, inset), the angular coefficients of dlnσ/dlnT = 0.1 in the temperature interval T = (40–80) K and 1–1.2 in the interval T = (100–300) K. Resistivity ratios ρ(300 K)/ρ(20 K) were obtained as ρ = 1/σ from the inset: 1–5, 2–6.1, 3–5.3, 4–7.9, 5–10.2, 6–10.2.
Figure 11. The temperature dependencies of the electronic thermal conductivity κe and the electrical conductivity σ (inset) in the p–Bi0.5Sb1.5Te3 (1–4) and p–Bi2Te3 (5, 6) films. The designations of the samples (1–6) correspond to those in the Table 1, and 7 and 8 are the bulk samples of the p–Bi0.5Sb1.5Te3 solid solution. In the films of p–Bi0.5Sb1.5Te3 (1–3, inset), the angular coefficients of dlnσ/dlnT = 0.1 in the temperature interval T = (40–80) K and 1–1.2 in the interval T = (100–300) K. Resistivity ratios ρ(300 K)/ρ(20 K) were obtained as ρ = 1/σ from the inset: 1–5, 2–6.1, 3–5.3, 4–7.9, 5–10.2, 6–10.2.
Magnetochemistry 09 00141 g011
Table 1. The p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films formation technique parameters.
Table 1. The p–Bi0.5Sb1.5Te3 and p–Bi2Te3 films formation technique parameters.
No.Formation Technique Substrate Heat Treatment Thermoelectric Power Coefficient α, μV K−1
p–Bi0.5Sb1.5Te3
1discrete evaporation polyimideunannealed242
2discrete evaporation polyimideannealed215
3thermal evaporation polyimideunannealed 200
4discrete evaporation muscoviteannealed223
p–Bi2Te3
5discrete evaporation polyimideunannealed 234
6thermal evaporation muscoviteunannealed 203
Table 2. The average grain area <S> and areas of the grains of various sizes S1–S4 determined from the images of surface morphology (0001) in the unannealed (sample 1) and the annealed (sample 4) p–Bi0.5Sb1.5Te3 films.
Table 2. The average grain area <S> and areas of the grains of various sizes S1–S4 determined from the images of surface morphology (0001) in the unannealed (sample 1) and the annealed (sample 4) p–Bi0.5Sb1.5Te3 films.
No.
in Table 1
<S>,
µm2
S1,
µm2, %
S2
µm2, %
S3
µm2, %
S4
µm2, %
10.1350.001
45%
0.002–0.008
43%
0.01–0.05
10%
2.6–8.7
2%
41.6250.002–0.077
72%
0.15–0.92
18%
(9.5–41.25)
10%
Table 3. Dependence of resistivity ρ and electrical conductivity σ on Bi2Te3 film thickness.
Table 3. Dependence of resistivity ρ and electrical conductivity σ on Bi2Te3 film thickness.
SampleT, Kt ρ, µΩ mσ, Ω−1cm−1
[25]10030 nm5.331877
[25]10016 nm3.113217
[25]2030 nm4.442252
[25]2016 nm1.7765630
51001 µm2.334275
61001 µm2.24532
5201 µm0.9810230
6201 µm1.128930
Samples 5 and 6 of the data in Table 3 are taken from Figure 11, inset.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lukyanova, L.N.; Boikov, Y.A.; Usov, O.A.; Danilov, V.A.; Makarenko, I.V.; Petrov, V.N. Thermal Conductivity for p–(Bi, Sb)2Te3 Films of Topological Insulators. Magnetochemistry 2023, 9, 141. https://doi.org/10.3390/magnetochemistry9060141

AMA Style

Lukyanova LN, Boikov YA, Usov OA, Danilov VA, Makarenko IV, Petrov VN. Thermal Conductivity for p–(Bi, Sb)2Te3 Films of Topological Insulators. Magnetochemistry. 2023; 9(6):141. https://doi.org/10.3390/magnetochemistry9060141

Chicago/Turabian Style

Lukyanova, Lidia N., Yuri A. Boikov, Oleg A. Usov, Viacheslav A. Danilov, Igor V. Makarenko, and Vasilii N. Petrov. 2023. "Thermal Conductivity for p–(Bi, Sb)2Te3 Films of Topological Insulators" Magnetochemistry 9, no. 6: 141. https://doi.org/10.3390/magnetochemistry9060141

APA Style

Lukyanova, L. N., Boikov, Y. A., Usov, O. A., Danilov, V. A., Makarenko, I. V., & Petrov, V. N. (2023). Thermal Conductivity for p–(Bi, Sb)2Te3 Films of Topological Insulators. Magnetochemistry, 9(6), 141. https://doi.org/10.3390/magnetochemistry9060141

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop