Review Article  Open Access
Ian T. Witting, Francesco Ricci, Thomas C. Chasapis, Geoffroy Hautier, G. Jeffrey Snyder, "The Thermoelectric Properties of Type Bismuth Telluride: Bismuth Selenide Alloys ", Research, vol. 2020, Article ID 4361703, 15 pages, 2020. https://doi.org/10.34133/2020/4361703
The Thermoelectric Properties of Type Bismuth Telluride: Bismuth Selenide Alloys
Abstract
Alloying bismuth telluride with antimony telluride and bismuth selenide for  and type materials, respectively, improves the thermoelectric quality factor for use in room temperature modules. As the electronic and thermal transports can vary substantially, the alloy composition is a key engineering parameter. The type Bi_{2}Te_{3}Se alloy lags its type counterpart in thermoelectric performance and does not lend itself as readily to simple transport modeling which complicates engineering. Combining literature data with recent results across the entire alloy composition range, the complex electronic structure dynamics and trends in lattice thermal conductivity are explored. Spinorbit interaction plays a critical role in determining the position and degeneracy of the various conduction band minima. This behavior is incorporated into a twoband effective mass model to estimate the transport parameters in each band. An alloy scattering model is utilized to demonstrate how phonon scattering behaves differently on either side of the intermediate ordered compound Bi_{2}Te_{2}Se due to chalcogen site occupancy preference. The parametrization of the electronic and thermal transports presented can be used in future optimization efforts.
1. Introduction
Bismuth telluride is the dominant thermoelectric material for applications near room temperature due to its inherently low lattice thermal conductivity and high electronic weighted mobility [1]. Its performance is ultimately limited by the deleterious effects of thermally generated minority carriers as a result of its small band gap (0.14 eV) [2]. This effect can be partially mitigated by doping more heavily than what would be ideal considering the transport of majority carriers alone. Alloying with Sb_{2}Te_{3} for type or Bi_{2}Se_{3} for type materials can also reduce bipolar effects by increasing the band gap, but the impact on thermal and electric transports must also be considered [3–6]. Thermal conductivity can be significantly reduced by the enhanced scattering of high frequency phonons by introducing mass contrast and bonding changes [7–9]. Charge carrier mobility can also be disrupted during alloying while changes in electronic structure could prove beneficial or detrimental [10, 11]. Engineering bismuth telluridebased thermoelectrics require simultaneous consideration of these details.
The type Bi_{2}Te_{3}Sb_{2}Te_{3} alloy has been well characterized and modeled. The thermal conductivity variation with composition is well described using a mass contrast alloy scattering model where composition varies on a single cation site [12]. The valence band structure for both binary compounds is similar with the primary difference being a switch in the order of the two valence band maxima [1]. These band edges cross in energy near Bi_{0.5}Sb_{1.5}Te_{3} resulting in an enhancement of the thermoelectric power factor [13–15]. This convergence in energy combined with a significant reduction in thermal conductivity leads to this being the optimum composition.
The type Bi_{2}Te_{3}Bi_{2}Se_{3} alloy is less straightforward. Mass contrast alloy scattering modeling does not match thermal conductivity values reported in the literature, and there is wide variability in qualitative trends reported [8, 9, 12, 16]. The electronic band structure of Bi_{2}Te_{3} is comprised of conduction and valence band extrema with high valley degeneracy while the structure of Bi_{2}Se_{3} is far simpler with a direct gap between singly degenerate extrema at the point [1, 17]. This dramatic difference in Fermi surface complexity along with a peak in the optical band gap near the intermediate composition Bi_{2}Te_{2}Se suggests that complex band dynamics are in play [4, 5]. There has recently been a great deal of interest in the electronic structure of alloys as these materials are also topological insulators [18–23]. Accordingly, these studies have primarily focused upon the surface states and the variation with composition of the energy level of the Dirac point relative to band edges. While the bulk structures of the endmembers Bi_{2}Te_{3} and Bi_{2}Se_{3} have been extensively studied, a comprehensive description of the dynamics in the alloy has not been reported.
In this report, we aim to explain long standing unresolved issues in the alloy system. Combining experimental results from the literature and new electronic structure calculations, we describe how the electronic structure changes with composition. Electronic structure calculations for this system are highly sensitive to the magnitude of the spinorbit interaction which has complicated resolving these issues. This is then used to parameterize the electronic transport in type in an effective mass model. In addition, we present an alloy scattering model of the lattice thermal conductivity which provides clarity regarding compositional trends by accounting for the occupancy preference of two inequivalent chalcogen sites.
2. Crystal Structure
Bismuth telluride, bismuth selenide, and all intermediate alloys have the tetradymite crystal structure in the symmetry group (Figure 1) [24, 25]. The material is comprised of repeating quintuple layers of where the number in parentheses designates two inequivalent chalcogen () sites. Bismuth atoms are octahedrally coordinated by chalcogen atoms, and the site atoms are octahedrally coordinated by Bi. The site atoms are covalently bonded with three Bi atoms and by weaker van der Waals bonds with three other atoms. The Bi bond length is close to the expected covalent bond length, while the Bi bond is longer and near the value expected for ionic bonding [25, 26]. The site is preferentially occupied by Se [25–28]. This site selectivity has been attributed to a reduction in bond angle strain and to the more electronegative Se occupying the site with the highest cation coordination [28, 29]. This is a key difference from the type Bi_{2}Te_{3}Sb_{2}Te_{3} alloy system where Bi and Sb occupy equivalent sites without preference.
The difference in bonding at each chalcogen site and the occupancy preference has significant consequences for the electrical and thermal properties of the system. The preference results in the formation of an ordered compound at Bi_{2}Te_{2}Se which acts as a demarcation point between which the chalcogen site is being changed during alloying. There is some disagreement in the literature as to whether the ordered compound is stable or if solid state phase separation occurs as observed in Bi_{2}Te_{2}S [29]. In Bi_{2}Te_{2}S, this separation occurs to relieve bond angle strains that would otherwise occur in the ordered Te^{(1)}–Bi–S^{(2)}–Bi–Te^{(1)} compound. Ordering of the Bi_{2}Te_{2}Se phase has been observed in numerous studies [25, 26, 30, 31], but others report that the phase will separate into two compounds with compositions near and if annealed at temperatures near 300°C for extended periods [32, 33]. If phase separation does occur, it may be kinetically inhibited at temperatures typically used for thermoelectric or topological insulator applications. This report assumes that this inhibition occurs in typical thermoelectric samples.
The layered structure of the tetradymites results in significant anisotropy in transport properties. In the type alloy system, the ratio of the electrical conductivity within the planes to that in the direction parallel to the axis ranges from 3 to 7, while the thermal conductivity anisotropy ratio is closer to 2 [34–36]. As the Seebeck coefficient is nearly isotropic when not demonstrating intrinsic conduction, the ratio is 23 necessitating the use of single or oriented polycrystals to maximize performance. This contrasts with the type Bi_{2}Te_{3}Sb_{2}Te_{3} alloy system where the electrical and thermal conductivity ratios are both near 23 resulting in a nearly isotropic [37]. This isotropy in the type system allows for significant improvements in to be made by nanostructuring for the reduction of lattice thermal conductivity. [38, 39] Nanostructuring can produce benefits in the type alloys as well, but they are smaller and require additional processing steps to preferentially orient the grains [40–44]. Figureofmerit values with nanotexturing have been reported between in type Bi_{0.5}Sb_{1.5}Te_{3}, but only in the range of 11.2 for alloys with . The experimental data presented in this study come from single or highly oriented large grain polycrystals, and properties are measured along the high direction in the planes.
3. Electronic Transport
The electronic transport properties vary significantly with alloy composition as shown in Figure 2. Pisarenko plot fits of Seebeck versus Hall carrier concentration using a single valley effective mass model assuming deformation potential scattering find that the Seebeck effective mass, , decreases from 1.06 m_{e} for Bi_{2}Te_{3} to 0.25 m_{e} for Bi_{2}Se_{3} [1, 8, 35, 45–50]. Similarly, the weighted mobility, μ_{w}, which sets the maximum achievable power factor, decreases monotonically from 590 to 170 cm^{2}V^{1} s^{1} going from Bi_{2}Te_{3} to Bi_{2}Se_{3} [16]. This is likely due to a loss of Fermi surface complexity by a decrease in valley degeneracy and/or by a decrease of the conduction anisotropy of each valley [51].
In the Bi_{2}Te_{3}Sb_{2}Te_{3} alloy system, a peak in the effective mass and an abrupt change in the band gap slope with composition is observed near Bi_{0.5}Sb_{1.5}Te_{3} [6, 52]. This has been attributed to the crossing in energy with composition of two valence bands [13]. In the type alloy system, there is a peak in the band gap with alloy composition very near Bi_{2}Te_{2}Se and also a peak in Seebeck effective mass of ~1.30 m_{e} near Bi_{2}Te_{2.5}Se_{0.5} (Figure 2) [1, 4, 5, 8, 53–55]. As the compositions for two peaks do not coincide, the transport behavior cannot be described by a simple crossing of two conduction bands and a more detailed investigation is required.
4. Electronic Structure
The electronic structure of the alloy system has been studied extensively due to the importance of these materials as thermoelectrics and topological insulators; however, the complete picture of dynamics of key band edges has not been laid out. This section summarizes the experimental and theoretical work for these alloys which is then used to inform our picture of how the band structure evolves with composition. It should be noted that both experimental and theoretical characterizations of the band structure are performed at low temperatures, and positions and shapes of band extrema may be shifted at temperatures where these materials operate as thermoelectrics [4, 56].
Experimental characterization of the Bi_{2}Te_{3} conduction band Fermi surface has found the conduction band minimum to be sixfold degenerate and nonparabolic [57–64]. Calculations have found this minimum to be within the bisectrix plane (Figure 3(a), highlighted) and slightly displaced from the line, a point we designate as [1]. Shubnikovde Haas (SdH) measurements also detected a second, heavier band edge only 25 meV above the CBM or at a carrier concentration of [64]. This pocket has been calculated to be twofold degenerate and lying on the line, a point designated here as . When both the and pockets are occupied, the conduction band Fermi surface forms a “double tripod” structure which has been observed by angleresolved photoemission spectroscopy (ARPES) [20]. A third pocket higher in energy than is predicted to be present in a position slightly displaced from the line making it sixfold degenerate. Experimental characterization of this edge has not been made; however, it cannot be ruled out that at operating temperatures and elevated doping, it potentially plays a role in transport.
ARPES and SdH measurements have found the CBM in Bi_{2}Se_{3} to be a single ellipsoidal valley centered at and elongated along the trigonal axis [22, 65, 66]. At a Fermi level 160 meV above the CBM or , a steep increase in with the Fermi level is observed which may indicate the filling of a second band edge having a larger densityofstates [65, 67]. Electronic structure calculations for Bi_{2}Se_{3} (Figures 3(c) and 4(a)) qualitatively match with experimental observations. Near the CBM edge, the Fermi surface is a simple ellipsoid centered at . Our calculations find a second band edge along the line (Figure 3(e)), similar to the CBM of Bi_{2}Te_{3} at which may explain the rise in densityofstates observed by SdH 160 meV above the CBM. To the best of the authors’ knowledge, this second edge has not been directly imaged by ARPES. Chen et al. used ARPES to observe the conduction band up to a Fermi level of 150 meV above the minimum at [22]. This is below where the second edge was potentially detected by SdH.
Köhler et al. performed Shubnikovde Haas measurements on alloy samples near the binary endpoints: and [68]. On the Bi_{2}Te_{3} side, the energy separation between the CBM () and the second band edge () linearly decreases to zero at . Within this composition range, the orientation and masses of the ellipsoidal CBM does not appear to change significantly. For the seleniderich alloys, the SdH measurements found a slight decrease of ~15 meV in the energy offset at between the CBM () and the proposed second conduction band (). The precise magnitude of the offset decrease could not be determined due to the error margins in the experiment. ARPES performed on an alloy crystal with a composition near Bi_{2}Te_{1.5}Se_{1.5} does not clearly show a conduction band pocket near ; however, there is a triangular Fermi surface centered at which could be a trigonally warped or pocket [69].
Band structure calculations on the ordered Bi_{2}Te_{2}Se compound disagree regarding the order and location of the conduction band extrema. Reports can be found with the CBM occurring along the line and a second edge along the [70], the same relative locations but in reversed order [71–73], or having only a single valley at the point [2]. These discrepancies arise due to the details of how these various calculations are performed and the sensitivity of results near the band gap to the effects of the spinorbit interaction (SOI) [2, 17, 18]. Further, the true location of band extrema can be missed if calculations were only performed along high symmetry directions. This is problematic as detailed studies have found the CBM at in Bi_{2}Te_{3} to be near but not on the line [1, 17, 74–80]. Evaluating the results of these calculations is complicated by a lack of experimental observations of the conduction band Fermi surface in the literature. ARPES has been performed on Bi_{2}Te_{2}Se; however, that study focused on the Dirac point in the surface states very near the VBM and did not capture the bulk conduction band [23].
In the absence of SOI, both the conduction band minimum (CBM) and valence band maximum (VBM) would be single valleys located at , making these materials poor thermoelectrics and having no topological insulating properties. The CBM would consist primarily of Bi 6p states while the VBM would be comprised of Se 4p and Te 5p states. Instead, the SOI shifts the energy levels of these bands such that the Bi states and chalcogen states reverse their ordering at [2, 19, 81, 82]. Where the bands overlap near , an anticrossing occurs due to their opposite parity thereby opening a band gap. Now at , the lowest energy conduction band states have chalcogen character and the valence band more Bi character. A schematic of the SOIinduced overlap effect is shown in Figure 3(b). In Bi_{2}Te_{3}, the overlap is significant enough to invert the curvature of the bands at and shift the band extremum to a lower symmetry and higher valley degeneracy points. While the overlap does occur in Bi_{2}Se_{3}, it is not enough to invert the curvature at and shift the extrema to different locations in the Brillouin zone [18, 19].
The magnitude of the SOI induced overlap increases monotonically from Bi_{2}Se_{3} to Bi_{2}Te_{3} when alloying [83]. Alloying between Bi_{2}Te_{3} and Bi_{2}Se_{3} can be qualitatively understood as modulating the magnitude of the SOIinduced band overlap. This is illustrated by electronic structure calculations in Figures 3(c)–3(f) for Bi_{2}Se_{3} with varying degrees of SOI. Doubling the SOI for Bi_{2}Se_{3} inverts the curvature at creating the doubly degenerate valley at z. Furthermore, the energy level of the pocket has moved above that of the sixfold degenerate pocket at . Note that this band structure closely resembles that of Bi_{2}Te_{3} where the SOI effect is greater (Figures 4(c) and 4(d)).
Knowing the effects of varying SOI allows for understanding of trends in the band gap with alloy composition. The opening of the band gap when substituting Te into Bi_{2}Se_{3} results from the CBM at Γ moving higher in energy and the VBM lower due to the increased effects of the SOIinduced anticrossing. This explains why the band gap increases instead of the decrease expected when substituting less electronegative anions [84]. The peak in band gap near Bi_{2}Te_{2}Se is not caused by the crossing of the CBM of Bi_{2}Te_{3} and the CBM of Bi_{2}Se_{3} as has been previously suggested [85]. The pocket is formed from the curvature inversion of , and the crossing of and was found at by SdH [68], coinciding with the observed peak in Seebeck effective mass (Figure 2). Any movement of the CB edge along also cannot explain the peak in band gap as it is calculated to be higher in energy than and for all compositions (Figure 3). Thus, the relative shifts of valence band edges must also be considered. In Bi_{2}Te_{3}, the VBM is located near the line with a second edge at , while in Bi_{2}Se_{3}, the VBM is found at . Note that the valence band at also inverts and shifts its extremum to in more Terich alloys. At some composition, the valley near the line and valence band edge must cross the and this may be the origin of the band gap peak. Recent interband absorption studies on this alloy system have found that at Bi_{2}Te_{2}Se, there is a change in the location of the lowest energy direct transition [55]. This may be a switch from the VB to CB transition in to the VB to CB for . As the and cross at in the conduction band, this would support a crossing between the valence band counterparts at .
5. Effective Mass Modeling of Conduction Band Electronic Transport
The shifts in band edge energies discussed in the previous section have been summarized schematically in Figure 5(a). With this picture established, an effective mass model can be fit that parameterizes the electronic transport in the and conduction band edges across the whole range of compositions. While charge transport in most materials is far more complex than that of nearly free electrons, an effective mass model that parameterizes the transport with empirical parameters equivalent to those of bands with a parabolic dispersion can make accurate predictions and provide useful physical insight [86]. The parameters of this model are determined empirically and should be interpreted as such.
The Seebeck coefficient, ; electrical conductivity, ; weighted mobility, ; and Hall coefficient, , for carriers in a single band where acoustic phonon scattering dominates can be described using
Here, is a Fermi integral of order and a function of the reduced chemical potential , is the average longitudinal elastic constant, is the valley degeneracy, is the acoustic deformation potential, is the conductivity effective mass, and is the Seebeck effective mass. When multiple band edges are contributing to transport, their combined contributions can be expressed using
A more complete effective mass model for this system would contain the variation with composition for each band of the Seebeck and conductivity masses, the deformation potentials, and the band degeneracies in addition to changes in lattice stiffness and alloy scattering energies. Unfortunately, most of these parameters have not been reported in the literature, and the large number of unknown variables leads to fit solutions which are not unique. Despite this difficulty, some conclusions may still be drawn regarding the band parameters.
Kohler reported that the lowest conduction band in Bi_{2}Te_{3} was sixfold degenerate with a densityofstates mass of 0.27 m_{e} and a second conduction band was present only 25 meV higher having a much higher mass estimated to be ~3 m_{e} [64]. Single band edge effective mass modeling of an extensive amount of literature data for type Bi_{2}Te_{3} found the relationship between conductivity and Seebeck could be well described by a weighted mobility of and the Hall and Seebeck by a Seebeck effective mass of [1]. If instead a twoconduction band model is fit to the same experimental data while using the masses and energy offsets reported by Kohler, the CBM at has a while the second minimum at has . The exact ratio of the calculated weighted mobilities will depend upon the true mass and offset of the pocket; however, a twoband system with a CBM of 0.27 m_{e} behaving like a single band with a mass of 1.06 m_{e} necessitates the second band having both a higher mass and weighted mobility. This result is contrary to expectation considering the pocket has higher degeneracy and lower effective mass and . This suggests that the acoustic deformation potential of the pocket must be significantly larger than that of the . When fitting a similar model for type Bi_{2}Te_{2.7}Se_{0.3}, Konstantinov et al. found that the deformation potential of the sixfold pocket must be 68 times larger than the second, higher mass band edge [87]. Furthermore, the higher mass and weighted mobility implies that the large majority of carriers and nearly half of the electrical conductivity in type Bi_{2}Te_{3} is attributable to the pocket. This is a key departure from the typical description of the conduction band as simply being sixfold degenerate.
In Bi_{2}Se_{3}, the lowest band is singly degenerate at with a band edge mass of 0.155 m_{e}; however, the nonparabolicity of the band makes this a poor descriptor of the densityofstates beyond carriers [67]. Instead, a slightly higher Seebeck mass would be more accurate for modeling at thermoelectric doping levels. Single band effective mass modeling finds that a mass of 0.25 m_{e} accurately predicts the Seebeck data for Hall carrier concentrations between and carriers. A potential second band edge was detected 160 meV above the CBM at a carrier concentration near carriers whose mass was only reported as “high” with no estimated value. The absence of any significant deviation with increased doping from the trend using the single band effective mass indicates that the CBM at dominates transport in type Bi_{2}Se_{3}. The weighted mobility of any second band must be significantly lower than that of the for a single band mass model to fit. Our electronic structure calculations predict this second band to be located at (Figure 3(e)). As is a sixfold degenerate point, the lower weighted mobility could result from some combination of a large conductivity mass or a large deformation potential.
We present in Figure 5(b) the result of a parametrization of the transport in this alloy system using a twoband effective mass model which follows the general trends in deformation potentials and masses previously discussed. The fitting parameters used for each band edge are the Seebeck effective mass, ; the deformation potential, ; and the ratio of the valley degeneracy to conductivity effective mass, . The valley degeneracy to conductivity mass ratio is not broken into its constituents as the value of valley degeneracy is not clearly defined for the as it transitions between in Bi_{2}Se_{3} and in Bi_{2}Te_{3}. Similarly, predicting the conductivity mass variations requires more significant knowledge than is available of pocket anisotropy changes with alloy composition. The values presented here are meant to serve as reasonable estimates based upon available data to be used for further experimental and theoretical verification.
Linear variations are assumed for fitting parameters between the binary compositions. Wherever possible, data from experimental results and band structure calculations are used to set parameters of the model. It must be noted that many of these experimental values were measured at temperatures lower than 300 K where the transport data was taken. Band offsets and masses could shift with temperature, and the degree to which this occurs is not clear. The values solved for in fitting are given in Table 1. These values were obtained by minimizing a sum of squared errors between the model and data points of Figure 2. The band gap and electronic transport in  and type bismuth telluride alloys show important differences in compositional trends due to the details of their complex band dynamics. In , a change in slope of band gap with composition coincides with a peak in Seebeck mass and weighted mobility when modeled using a single band [6, 13–15]. This correlation is explained by a crossing of two valence bands in energy with composition [6]. In , a peak in band gap does not coincide in composition with the peak in effective mass and no peak in weighted mobility is observed [4, 8, 16, 53, 54]. This suggests that the underlying causes differ for the mass and band gap peaks. The solid lines in this figure are simple guides to the eye, while the dashed lines represent single band evaluation of the simplified twoconduction band effective mass model discussed within the text and evaluated for a constant doping level of .

We attempted to include in our model the effects of alloy scattering of charge carriers as described by Harrison and Hauser [11]; however, best fit values obtained were very small (<1 meV) and did not improve the model fit. We therefore exclude these effects from the model. It is not clear why alloy scattering does not play a prominent role in transport as it does in other semiconductor alloys [88]. One possibility is that deformation potential scattering by phonons is strong enough to minimize the impact of other scattering mechanisms. This appears to be the case for the lack of significant ionized impurity scattering near room temperature in Bi_{2}Te_{3} due to the relatively soft lattice and large static dielectric permittivity [1]. The relaxation times for carrier scattering have and dependence for alloy and deformation potential scatter, respectively [11, 89]. The impact of alloy scattering therefore may be more important at temperatures below 300 K where the data used in this study was collected.
Shown in Figure 5(b) are the Seebeck masses and resulting weighted mobilities for the two conduction bands which produce the fits of Figure 2 (dashed curves). The weighted mobility of both pockets decreases with increasing Se content of the alloy, however, for different reasons. The pocket presumably remains sixfold degenerate, but its mass increases (decreasing ) leading to a reduction in mobility. The pocket decreases in from 2 in Bi_{2}Te_{3} to 1 in Bi_{2}Se_{3}; however, its mass decreases substantially such that the weighted mobility would still be expected to increase. This suggests that the deformation potential of the band increases from Bi_{2}Te_{3} to Bi_{2}Se_{3}. Our model fits this increase as from 5 to 16 eV. An improved fit could be obtained by not assuming a linear variation in parameters across the entire alloy composition; however, experimental and theoretical data do not provide sufficient guidance to avoid overfitting. Future studies could refine the understanding of this important thermoelectric and topological insulator material system.
6. Thermal Conductivity of Alloys
Much of thermoelectric material engineering pertains to reducing a material’s lattice thermal conductivity while still maintaining a high weighted mobility. This is possible because of the difference in order of magnitudes of the mean free paths of charge carriers (small) and phonons (small to large). Nevertheless, most methods of reducing thermal conductivity also reduce carrier mobility somewhat, and evaluation of the success in balancing electrical and thermal engineering can be performed using the quality factor [3]. Successful reduction of lattice thermal conductivity involves introducing effective scattering across the entire range of phonon frequencies responsible for heat transport [90]. Alloying reduces thermal conductivity by enhancing phonon scattering due to the influence of mass fluctuations and localized strains. With a relaxation time proportional to , it is most effective at scattering high frequency phonons [7]. Note that the relaxation time of other phonon scattering mechanisms common in thermoelectrics has different frequency dependencies and thus scatters different wavelengths of phonons: boundary scattering , dislocation core , and dislocation strain [91]. These contributions are highly dependent upon material processing. As we are concerned with the pure material properties of these alloys, the data considered within this section is from single or highly oriented large grain polycrystals and thus should be influenced by alloy scattering alone.
The effect on lattice thermal conductivity of alloying Bi_{2}Te_{3} and Bi_{2}Se_{3} is shown in Figure 6(a). There are large variations in the lattice thermal conductivity reported in the literature; however, most studies report a global minimum near Bi_{2}Te_{2.5}Se_{0.5}. To gain further insight, the alloy scattering model of Callaway and von Baeyer and Klemens was fit to the data of each study as shown in Figure 6(b) [92, 93]. Such models have been useful in explaining the thermal conductivity trends in lead chalcogenides [3], bismuth antimony telluride [12], SiGe [94], and half Heuslers [95].
The lattice thermal conductivity of the alloy (Equation (8)), , is expressed as a function of the of the “pure” compound which is the linear extrapolation between the two alloy endpoints, and of the disorder scattering parameter, , which is defined in Equation (9). The disorder scattering parameter is a function of the Debye temperature, ; the average atomic volume, ; the average speed of sound, ; and the scattering parameter, . This scattering parameter is found in multiple forms in the literature. Due to the occupancy preference of the chalcogen sites and the presence of an ordered compound at Bi_{2}Te_{2}Se, the Bi_{2}Te_{3}Se alloy system can be viewed as two separate alloy systems: Bi_{2}Te_{3}Bi_{2}Te_{2}Se and Bi_{2}Te_{2}Se and Bi_{2}Se_{3}. In order to account for alloying on two different chalcogen sublattices, we use the scattering parameter form of Yang et al. which was used to describe the effects of alloying on each site in ZrNiSn [95].
Here, the scattering parameter (Equation (10)) is broken into the summation of a mass fluctuation (Equation (11)), , and a strain field fluctuation parameter (Equation (12)), . The fluctuation parameters sum over all sublattices, three for tetradymites, and the occupancy of the th sublattice having a degeneracy of . Both are functions of the fractional occupancy, , by species ; the average atomic mass of the sublattice, ; and the average mass of the compound, . The mass fluctuation parameter varies with the difference in mass, , between the alloying species, and the strain field on the difference in atomic radius, , and a parameter . This parameter was introduced by Abeles who derived it being dependent upon the Grüneisen parameter, bulk modulus of the matrix and the sphere of the alloying atom, and Poisson’s ratio. Exact calculation of all these parameters is difficult. Thus, in practice, is a phenomenological, adjustable parameter determined when fitting experimental data and ranges between 0 (no strain scattering) and 100 [3, 95]. The parameter arises from a slight modification of the original scattering parameter derivation of Klemens which instead contains a term related to the average stiffness constant of nearest neighbor bonds [7]. While it may appear arbitrary in usage, the magnitude of can be interpreted as a general indicator of the importance of strain scattering on lattice thermal conductivity due to alloying. This suggests strain engineering by careful alloying could be an avenue for further thermal conductivity reduction, as has recently been reported for Eu and Mndoped PbTe [96].
While there is discrepancy in the reported endpoint lattice thermal conductivities of Bi_{2}Te_{3}, Bi_{2}Te_{2}Se, and Bi_{2}Se_{3}, the leastsquarederror fit parameters for each alloying range listed in Table 2 are similar. This is the case even for the data of Rosi et al. that found a local maximum in near Bi_{2}Te_{1.5}Se_{1.5} [9]. No other report to our knowledge has observed a similar maximum. The majority of Rosi et al.’s data still fall within the same bounds of the model as other studies which suggests the local maximum is an outlier. Using the average endpoint and fitting parameters for each study, a model was generated considering mass fluctuation and strain effects (Figure 6(b)). The key role of Se/Te site preference emerges from this analysis. Between where the composition changes only the occupancy of the (1) site, the variation in lattice thermal conductivity is primarily due to mass fluctuation effects (red dashed curve). This is not the case for , where some effect of strain and/or changes in bonding must be present (shaded region indicating reasonable bounds). The inability of mass fluctuation scattering alone to explain alloying effects contrasts with the Bi_{2}Te_{3}Sb_{2}Te_{3}type alloy where only mass fluctuation scattering is required [12, 97]. In the analogous compound Bi_{2}Te_{2}S, the more electronegative sulfur atom preferentially occupies the (2) site. At this composition, large bond angle strains prevent the formation of the ordered compound, and instead, some of the (2) site are occupied by Te [29]. Similar but smaller strains must occur when alloying Se on the (2) site as between Bi_{2}Te_{3} and Bi_{2}Te_{2}Se.

7. Conclusion
Transport modeling of the type alloy system proves more challenging than the type alloys and complicates engineering thermoelectric materials. Differences in the SOI between the telluride and selenide lead to changes in location of the conduction band minimum with composition as well as changes in valley degeneracy and effective mass of some band edges. We have summarized the current state of experimental and theoretical knowledge of these band dynamics and parametrized the transport using a simple effective mass model to supply estimates of their values. Reports on the lattice thermal conductivity in the type alloy have yielded unclear results regarding trends and location of minimum values. Using an alloy scattering model incorporating the ordered Bi_{2}Te_{2}Se compound provides clarity and identifies bond strain as a key factor between Bi_{2}Te_{3} and Bi_{2}Te_{2}Se. It is our hope that the analyses provided here will aid future efforts to characterize and engineer type alloys.
8. Experimental Methods
8.1. Electronic Structure Calculations
The Bi_{2}Te_{3} (), Bi_{2}Se_{3}, and Bi_{2}Te_{2}Se () compounds were taken in their trigonal crystal structure (space group R3m) with experimental cell parameters from the Materials Project as shown in Table 3 [98].

For the composition , a supercell was built starting from the primitive cell of composition . The supercell parameters were taken as a linear interpolation between the supercell parameters for and compositions. Then, four Se atoms were replaced with four Te atoms. Finally, a smaller primitive cell with only 10 atoms was found.
The electronic band structures were calculated with Density Functional Theory (DFT), using the Vienna Ab initio Simulation Package (VASP) [99, 100], using the PerdewBurkeErnzerhof (PBE) [101] generalized gradient approximation (GGA) functional and adopting the projector augmentedwave (PAW) [102, 103] approach. Spinorbit interaction (SOI) was included in all the calculations to obtain the complexity of the band structure. Also, the source code of VASP was edited in order to change the weight of SOI in the case of Bi_{2}Se_{3}. In order to get a band structure that fits better the experimental band structure, we set a SOC weight equal to 0.7, corresponding to 100% in Figure 3. Van der Waals interactions were not included.
The band structure on the bisectrix plane for the Bi_{2}Te_{3} and Bi_{2}Se_{3}, shown in Figure 4, was calculated with a nonselfconsistent field calculation on a grid of 10,404 points.
The band structure on the high symmetry path, shown in Figure 4, was calculated using 50 points in each segment. The information about the atomic character of the bands is shown using colors.
The band structure of the Bi_{2}Te_{2.5}Se_{0.5} compound was calculated using BandUP software [104, 105]. This software allows to unfold the bands of a supercell back to the standard high symmetry path of the primitive cell. It is worth mentioning that the unfolding procedure returns a weight of the projection of the eigenvalues to the specific path. Thus, the band structure in Figure 4(b) is plotted in gray scale representing the weight of the projection. The same technique was used to plot the bands on the bisectrix plane, but only the maximal weight of the projection was considered. This allowed a plot in two dimensions while using color to represent the band energy. Also, to limit the computational cost, a point mesh covering only the region occupied by the two pockets of interest was used.
All data analyses, such as band structure plotting and supercell generation, were carried out using the pymatgen python package [106].
8.2. Effective Mass Modeling
In a single band system with a parabolic energy dependence on , the empirically determined Seebeck mass and weighted mobility are independent of the doping level. In a system where multiple bands are contributing to transport, the Seebeck mass and weighted mobility determined using a single band model will vary with doping level. This complicates comparison as our goal in this report is to explain the variation of a single band effective mass with alloy composition. To allow for comparison, the dashed curves in Figure 2 were produced at a constant total carrier concentration of . Charge neutrality was used to determine the Fermi level in the multiband system at this doping level and subsequently calculate the transport using Equations (1) through (7). Once the total Seebeck, conductivity, and Hall values were calculated using the multiband model, single band mass and weighted mobility parameters could be then calculated. Note that while a constant dopant concentration was used to produce the dashed curves of Figure 2, the multiband model parameters were best fit considering the various measured Hall values from the cited works.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This research was carried out under a contract with the National Aeronautics and Space Administration and was supported by the NASA Science Missions Directorate’s Radioisotope Power Systems Technology Advancement Program. The authors also gratefully acknowledge thermoelectrics research at Northwestern University through the award 70NANB19H005 from U.S. Department of Commerce, National Institute of Standards and Technology as part of the Center for Hierarchical Materials Design (CHiMaD). FR acknowledges support from the LOCOTED (Low Cost ThermoElectric Devices) project funded by the Walloon Region (Programmes FEDER).
References
 I. T. Witting, T. C. Chasapis, F. Ricci et al., “The Thermoelectric Properties of Bismuth Telluride,” Advanced Electronic Materials, vol. 5, no. 6, article 1800904, 2019. View at: Publisher Site  Google Scholar
 H. Shi, D. Parker, M. H. Du, and D. J. Singh, “Connecting thermoelectric performance and topologicalinsulator Behavior:Bi_{2}Te_{3} and Bi_{2}Te_{2} Se from first principles,” Physical Review Applied, vol. 3, no. 1, pp. 1–10, 2015. View at: Publisher Site  Google Scholar
 H. Wang, A. D. Lalonde, Y. Pei, and G. J. Snyder, “The criteria for beneficial disorder in thermoelectric solid solutions,” Advanced Functional Materials, vol. 23, no. 12, pp. 1586–1596, 2013. View at: Publisher Site  Google Scholar
 D. L. Greenaway and G. Harbeke, “Band structure of bismuth telluride, bismuth selenide and their respective alloys,” Journal of Physics and Chemistry of Solids, vol. 26, no. 10, pp. 1585–1604, 1965. View at: Publisher Site  Google Scholar
 I. G. Austin and A. Sheard, “Some optical properties of Bi_{2}Te_{3}Bi_{2}Se_{3} alloys,” Journal of Electronics and Control, vol. 3, no. 2, pp. 236237, 1957. View at: Publisher Site  Google Scholar
 R. Sehr and L. R. Testardi, “The optical properties of ptype Bi_{2}Te_{3}Sb_{2}Te_{3} alloys between 215 microns,” Journal of Physics and Chemistry of Solids, vol. 23, no. 9, pp. 1219–1224, 1962. View at: Publisher Site  Google Scholar
 P. G. Klemens, “The scattering of lowfrequency lattice waves by static imperfections,” Proceedings of the Physical Society Section A, vol. 68, no. 12, pp. 1113–1128, 1955. View at: Publisher Site  Google Scholar
 U. Birkholz, “Untersuchung der intermetallischen verbindung Bi_{2}Te_{3} sowie der festen lösungen Bi_{2x}Sb_{x}Te_{3} und Bi_{2}Te_{3x}Se_{x} hinsichtlich ihrer eignung als material für halbleiterThermoelemente,” Zeitschrift für Naturforsch, vol. 13, no. 9, pp. 780–792, 1958. View at: Publisher Site  Google Scholar
 F. D. Rosi, B. Abeles, and R. V. Jensen, “Materials for thermoelectric refrigeration,” Journal of Physics and Chemistry of Solids, vol. 10, no. 23, pp. 191–200, 1959. View at: Publisher Site  Google Scholar
 Y. Pei, H. Wang, and G. J. Snyder, “Band engineering of thermoelectric materials,” Advanced Materials, vol. 24, no. 46, pp. 6125–6135, 2012. View at: Publisher Site  Google Scholar
 J. W. Harrison and J. R. Hauser, “Alloy scattering in ternary IIIV compounds,” Physical Review B, vol. 13, no. 12, pp. 5347–5350, 1976. View at: Publisher Site  Google Scholar
 H. J. Goldsmid, “Recent studies of bismuth telluride and its alloys,” Journal of Applied Physics, vol. 32, no. 10, pp. 2198–2202, 1961. View at: Publisher Site  Google Scholar
 H.S. Kim, N. A. Heinz, Z. M. Gibbs, Y. Tang, S. D. Kang, and G. J. Snyder, “High thermoelectric performance in (Bi_{0.25}Sb_{0.75})_{2}Te_{3} due to band convergence and improved by carrier concentration control,” Materials Today, vol. 20, no. 8, pp. 452–459, 2017. View at: Publisher Site  Google Scholar
 R. S. Erofeev and V. N. Ovechkina, “Special features of the energy spectrum of solid solutions of the Bi_{2}Te_{3}Sb_{2}Te_{3} system,” Izvestiya Akademii Nauk SSSR, Neorganicheskiye Materialy, vol. 17, no. 10, pp. 1780–1784, 1965. View at: Google Scholar
 V. S. Gaidukova, R. S. Erofeev, and V. N. Ovechkina, “Characteristics of the energy spectrum of solid solutions in the SbTeBiTe system,” Izvestiya Akademii Nauk SSSR, Neorganicheskiye Materialy, vol. 17, no. 2, pp. 244–247, 1981. View at: Google Scholar
 C. H. Champness, W. B. Muir, and P. T. Chiang, “Thermoelectric properties of ntype Bi2Te3–Bi2Se3ALLOYS,” Canadian Journal of Physics, vol. 45, no. 11, pp. 3611–3626, 1967. View at: Publisher Site  Google Scholar
 J. P. Heremans and B. Wiendlocha, “Tetradymites: Bi_{2}Te_{3}related materials,” in Materials Aspect of Thermoelectricity, C. Uher, Ed., pp. 39–94, CRC Press, 2016. View at: Publisher Site  Google Scholar
 J. P. Heremans, R. J. Cava, and N. Samarth, “Tetradymites as thermoelectrics and topological insulators,” Nature Reviews Materials, vol. 2, no. 10, article 17049, 2017. View at: Publisher Site  Google Scholar
 H. Zhang, C.X. Liu, X.L. Qi, X. Dai, Z. Fang, and S.C. Zhang, “Topological insulators in Bi_{2}Se_{3}, Bi_{2}Te_{3} and Sb_{2}Te_{3} with a single Dirac cone on the surface,” Nature Physics, vol. 5, no. 6, pp. 438–442, 2009. View at: Publisher Site  Google Scholar
 Y. L. Chen, J. G. Analytis, J.H. Chu et al., “Experimental realization of a threedimensional topological insulator, Bi_{2}Te_{3},” Science, vol. 325, no. 5937, pp. 178–181, 2009. View at: Publisher Site  Google Scholar
 D. O. Scanlon, P. D. C. King, R. P. Singh et al., “Controlling bulk conductivity in topological insulators: key role of antisite defects,” Advanced Materials, vol. 24, no. 16, pp. 2154–2158, 2012. View at: Publisher Site  Google Scholar
 Y. L. Chen, J.H. Chu, J. G. Analytis et al., “Massive Dirac fermion on the surface of a magnetically doped topological insulator,” Science, vol. 329, no. 5992, pp. 659–662, 2010. View at: Publisher Site  Google Scholar
 M. Neupane, S.Y. Xu, L. A. Wray et al., “Topological surface states and Dirac point tuning in ternary topological insulators,” Physical Review B, vol. 85, no. 23, article 235406, 2012. View at: Publisher Site  Google Scholar
 P. W. Lange, “Ein vergleich zwischen Bi_{2}Te_{3} und Bi_{2}Te_{2}S,” Naturwissenschaften, vol. 27, no. 8, pp. 133134, 1939. View at: Publisher Site  Google Scholar
 S. Nakajima, “The crystal structure of Bi_{2}Te_{3−x}Se_{x},” Journal of Physics and Chemistry of Solids, vol. 24, no. 3, pp. 479–485, 1963. View at: Publisher Site  Google Scholar
 J. A. Bland and S. J. Basinski, “The crystal structure of Bi_{2}Te_{2}Se,” Canadian Journal of Physics, vol. 39, no. 7, pp. 1040–1043, 1961. View at: Publisher Site  Google Scholar
 T. Zhu, L. Hu, X. Zhao, and J. He, “New insights into intrinsic point defects in V_{2}V_{I3} thermoelectric materials,” Advanced Science, vol. 3, no. 7, article 1600004, 2016. View at: Publisher Site  Google Scholar
 J. R. Drabble and C. H. L. Goodman, “Chemical bonding in bismuth telluride,” Journal of Physics and Chemistry of Solids, vol. 5, no. 1–2, pp. 142–144, 1958. View at: Publisher Site  Google Scholar
 L. Pauling, “The formula, structure, and chemical bonding of tetradymite, Bi_{14}Te_{13}Se_{8} and the phase Bi_{14}Te_{15}S_{6},” American Mineralogist, vol. 60, no. 11–12, pp. 994–997, 1975. View at: Google Scholar
 S. Misra and M. B. Bever, “On the solid solutions of bismuth telluride and bismuth selenide,” Journal of Physics and Chemistry of Solids, vol. 25, no. 11, pp. 1233–1241, 1964. View at: Publisher Site  Google Scholar
 S. Jia, H. Ji, E. ClimentPascual et al., “Lowcarrierconcentration crystals of the topological insulator Bi_{2}Te_{2}Se,” Physical Review B, vol. 84, no. 23, pp. 1–7, 2011. View at: Publisher Site  Google Scholar
 J. F. Dumas, G. Brun, B. Liautard, J. C. Tedenac, and M. Maurin, “New contribution in the study of the Bi_{2}Te_{3}Bi_{2}Se_{3} system,” Thermochimica Acta, vol. 122, no. 1, pp. 135–141, 1987. View at: Publisher Site  Google Scholar
 H. Ghoumari Bouanani, D. Eddike, B. Liautard, and G. Brun, “Solid state demixing in Bi_{2}Se_{3}Bi_{2}Te_{3} and Bi_{2}Se_{3}In_{2}Se_{3} phase diagrams,” Materials Research Bulletin, vol. 31, no. 2, pp. 177–187, 1996. View at: Publisher Site  Google Scholar
 M. Carle, P. Pierrat, C. LahalleGravier, S. Scherrer, and H. Scherrer, “Transport properties of n type Bi_{2}(Te_{1−x}Se_{x})_{3} single crystal solid solutions (x ⩽ 0.05),” Journal of Physics and Chemistry of Solids, vol. 56, no. 2, pp. 201–209, 1995. View at: Publisher Site  Google Scholar
 M. Stordeur, K. K. Ketavong, A. Priemuth, H. Sobotta, and V. Riede, “Optical and electrical investigations of ntype Bi_{2}Se_{3} single crystals,” Physica Status Solidi B: Basic Solid State Physics, vol. 169, no. 2, pp. 505–514, 1992. View at: Publisher Site  Google Scholar
 J.P. Fleurial, L. Gailliard, R. Triboulet, H. Scherrer, and S. Scherrer, “Thermal properties of high quality single crystals of bismuth telluride–Part I: experimental characterization,” Journal of Physics and Chemistry of Solids, vol. 49, no. 10, pp. 1237–1247, 1988. View at: Publisher Site  Google Scholar
 S. Scherrer and H. Scherrer, “Bismuth telluride, antimony telluride, and their solid solutions,” in CRC Handbook of Thermoelectrics, D. M. Rowe, Ed., CRC Press, 1995. View at: Publisher Site  Google Scholar
 S. I. Kim, K. H. Lee, H. A. Mun et al., “Thermoelectrics. Dense dislocation arrays embedded in grain boundaries for highperformance bulk thermoelectrics,” Science, vol. 348, no. 6230, pp. 109–114, 2015. View at: Publisher Site  Google Scholar
 B. Poudel, Q. Hao, Y. Ma et al., “Highthermoelectric performance of nanostructured bismuth antimony telluride bulk alloys,” Science, vol. 320, no. 5876, pp. 634–638, 2008. View at: Publisher Site  Google Scholar
 L. P. Hu, X. H. Liu, H. H. Xie, J. J. Shen, T. J. Zhu, and X. B. Zhao, “Improving thermoelectric properties of n type bismuthtelluridebased alloys by deformationinduced lattice defects and texture enhancement,” Acta Materialia, vol. 60, no. 11, pp. 4431–4437, 2012. View at: Publisher Site  Google Scholar
 X. Yan, B. Poudel, Y. Ma et al., “Experimental studies on anisotropic thermoelectric properties and structures of ntype Bi_{2}Te_{2.7}Se_{0.3},” Nano Letters, vol. 10, no. 9, pp. 3373–3378, 2010. View at: Publisher Site  Google Scholar
 R. J. Mehta, Y. Zhang, C. Karthik et al., “A new class of doped nanobulk highfigureofmerit thermoelectrics by scalable bottomup assembly,” Nature Materials, vol. 11, no. 3, pp. 233–240, 2012. View at: Publisher Site  Google Scholar
 W.S. Liu, Q. Zhang, Y. Lan et al., “Thermoelectric property studies on Cudoped ntype Cu_{x}Bi_{2}Te_{2.7}Se_{0.3} nanocomposites,” Advanced Energy Materials, vol. 1, no. 4, pp. 577–587, 2011. View at: Publisher Site  Google Scholar
 L. Hu, T. Zhu, X. Liu, and X. Zhao, “Point defect engineering of highperformance bismuthtelluridebased thermoelectric materials,” Advanced Functional Materials, vol. 24, no. 33, pp. 5211–5218, 2014. View at: Publisher Site  Google Scholar
 J. Navrátil, P. Lošťák, and J. Horák, “Transport coefficient of galliumdoped Bi_{2}Te_{3} single crystals,” Crystal Research and Technology, vol. 26, no. 6, pp. 675–681, 1991. View at: Publisher Site  Google Scholar
 P. Lošťák, L. Beneš, S. Civiš, and H. Süssmann, “Preparation and some physical properties of Bi_{2−x}In_{x}Se_{3} single crystals,” Journal of Materials Science, vol. 25, no. 1, pp. 277–282, 1990. View at: Publisher Site  Google Scholar
 V. A. Kulbachinskii, N. Miura, H. Nakagawa et al., “Conductionband structure of Bi_{2−x}Sb_{x}Se_{3} mixed crystals by Shubnikov–de Haas and cyclotron resonance measurements in high magnetic fields,” Physical Review B, vol. 59, no. 24, pp. 15733–15739, 1999. View at: Publisher Site  Google Scholar
 J. Navrátil, J. Horák, T. Plecháček et al., “Conduction band splitting and transport properties of Bi_{2}Se_{3},” Journal of Solid State Chemistry, vol. 177, no. 4–5, pp. 1704–1712, 2004. View at: Publisher Site  Google Scholar
 P. Janíček, Č. Drašar, L. Beneš, and P. Lošt’ák, “Thermoelectric properties of Tldoped Bi2Se3single crystals,” Crystal Research and Technology, vol. 44, no. 5, pp. 505–510, 2009. View at: Publisher Site  Google Scholar
 P. Lošt’ák, Č. Drašar, H. Süssmann, P. Reinshaus, R. Novotný, and L. Beneš, “Preparation and some physical properties of (Bi_{1−x}Sb_{x})_{2}Se_{3} single crystals,” Journal of Crystal Growth, vol. 179, no. 1–2, pp. 144–152, 1997. View at: Publisher Site  Google Scholar
 Z. M. Gibbs, F. Ricci, G. Li et al., “Effective mass and Fermi surface complexity factor from ab initio band structure calculations,” NPJ Computational Materials, vol. 3, no. 1, p. 8, 2017. View at: Publisher Site  Google Scholar
 M. Stordeur, M. Stölzer, H. Sobotta, and V. Riede, “Investigation of the valence band structure of thermoelectric (Bi_{1−x}Sb_{x})_{2}Te_{3} single crystals,” Physica Status Solidi B: Basic Solid State Physics, vol. 150, no. 1, pp. 165–176, 1988. View at: Publisher Site  Google Scholar
 G. N. Gordiakova, G. V. Kokosh, and S. S. Sinani, “Investigation of thermoelectric properties of solid solutions Bi_{2}Te_{3}Bi_{2}Se_{3},” Soviet PhysicsTechnical Physics, vol. 3, no. 1, pp. 1–14, 1958. View at: Google Scholar
 S. Wang, G. Tan, W. Xie et al., “Enhanced thermoelectric properties of Bi_{2}(Te_{1−x}Se_{x})_{3}based compounds as ntype legs for lowtemperature power generation,” Journal of Materials Chemistry, vol. 22, no. 39, p. 20943, 2012. View at: Publisher Site  Google Scholar
 A. Dubroka, O. Caha, M. Hronček et al., “Interband absorption edge in the topological insulators Bi_{2}(Te_{1−x}Se_{x})_{3},” Physical Review B, vol. 96, no. 23, article 235202, 2017. View at: Publisher Site  Google Scholar
 I. G. Austin, “The optical properties of bismuth telluride,” Proceedings of the Physical Society, vol. 72, no. 4, pp. 545–552, 1958. View at: Publisher Site  Google Scholar
 J. R. Drabble, “Galvanomagnetic effects in ptype bismuth telluride,” Proceedings of the Physical Society, vol. 72, no. 3, pp. 380–390, 1958. View at: Publisher Site  Google Scholar
 J. R. Drabble, R. D. Groves, and R. Wolfe, “Galvanomagnetic effects in ntype bismuth telluride,” Proceedings of the Physical Society, vol. 71, no. 3, pp. 430–443, 1958. View at: Publisher Site  Google Scholar
 L. R. Testardi, P. J. Stiles, and E. Burstein, “De Haasvan Alphen and high field galvanomagnetic studies of the Bi_{2}Te_{3} valence band structure,” Solid State Communications, vol. 1, no. 2, pp. 28–34, 1963. View at: Publisher Site  Google Scholar
 R. B. Mallinson, J. A. Rayne, and R. W. Ure Jr., “Concentration dependence of de Haasvan Alphen effect in ntype Bi_{2}Te_{3},” Physics Letters A, vol. 24, no. 13, pp. 713714, 1967. View at: Publisher Site  Google Scholar
 H. A. Ashworth, J. A. Rayne, and R. W. Ure, “Transport properties of Bi_{2} Te_{3},” Physical Review B, vol. 3, no. 8, pp. 2646–2661, 1971. View at: Publisher Site  Google Scholar
 R. B. Mallinson, J. A. Rayne, and R. W. Ure, “de Haasvan Alphen effect in ntype Bi_{2}Te_{3},” Physics Review, vol. 175, no. 3, pp. 1049–1056, 1968. View at: Publisher Site  Google Scholar
 H. Köhler, “Nonparabolicity of the highest valence band of Bi_{2}Te_{3} from Shubnikovde Haas effect,” Physica Status Solidi B: Basic Solid State Physics, vol. 74, no. 2, pp. 591–600, 1976. View at: Publisher Site  Google Scholar
 H. Köhler, “NonParabolic E(k) relation of the lowest conduction band in Bi_{2} Te_{3},” Physica Status Solidi B: Basic Solid State Physics, vol. 73, no. 1, pp. 95–104, 1976. View at: Publisher Site  Google Scholar
 H. Köhler and G. Landwehr, “Constant energy surfaces of ntype bismuth selenide from the Shubnikovde Haas effect,” Physica Status Solidi B: Basic Solid State Physics, vol. 45, no. 2, pp. K109–K112, 1971. View at: Publisher Site  Google Scholar
 K. Eto, Z. Ren, A. A. Taskin, K. Segawa, and Y. Ando, “Angulardependent oscillations of the magnetoresistance in Bi_{2}Se_{3} due to the threedimensional bulk Fermi surface,” Physical Review B, vol. 81, no. 19, article 195309, 2010. View at: Publisher Site  Google Scholar
 H. Köhler, “Conduction band parameters of _{Bi2See} from Shubnikovde Haas investigations,” Physica Status Solidi B: Basic Solid State Physics, vol. 58, no. 1, pp. 91–100, 1973. View at: Publisher Site  Google Scholar
 H. Köhler, W. Haigis, and A. von Middendorff, “Shubnikovde Haas investigations on ntype Bi_{2}(Se_{x}Te_{1−x})_{3},” Physica Status Solidi B: Basic Solid State Physics, vol. 78, no. 2, pp. 637–642, 1976. View at: Publisher Site  Google Scholar
 J.L. Mi, M. Bremholm, M. Bianchi et al., “Phase separation and bulk pn transition in single crystals of Bi_{2}Te_{2}Se topological insulator,” Advanced Materials, vol. 25, no. 6, pp. 889–893, 2013. View at: Publisher Site  Google Scholar
 P. Larson, V. A. Greanya, W. C. Tonjes, R. Liu, S. D. Mahanti, and C. G. Olson, “Electronic structure of Bi_{2}X_{3} (X = S, Se, T) compounds: comparison of theoretical calculations with photoemission studies,” Physical Review B, vol. 65, no. 8, article 085108, 2002. View at: Publisher Site  Google Scholar
 X.Q. Dai, B. Zhao, J.H. Zhao, Y.H. Li, Y.N. Tang, and N. Li, “Robust surface state of intrinsic topological insulator Bi_{2}Te_{2}Se thin films: a firstprinciples study,” Journal of Physics: Condensed Matter, vol. 24, no. 3, article 035502, 2012. View at: Publisher Site  Google Scholar
 L.L. Wang and D. D. Johnson, “Ternary tetradymite compounds as topological insulators,” Physical Review B, vol. 83, no. 24, article 241309, 2011. View at: Publisher Site  Google Scholar
 I. P. Rusinov, I. A. Nechaev, and E. V. Chulkov, “Manybody effects on the width of the band gap in Bi_{2}Te_{2}X (X = Te, Se, S) topological insulators,” JETP Letters, vol. 98, no. 7, pp. 397–402, 2013. View at: Publisher Site  Google Scholar
 S. J. Youn and A. J. Freeman, “Firstprinciples electronic structure and its relation to thermoelectric properties of Bi_{2}Te_{3},” Physical Review B, vol. 63, no. 8, article 085112, 2001. View at: Publisher Site  Google Scholar
 P. Larson, “Effect of p_{1/2} corrections in the electronic structure of Bi_{2}Te_{3} compounds,” Physical Review B, vol. 68, no. 15, article 155121, 2003. View at: Publisher Site  Google Scholar
 M. Kim, A. J. Freeman, and C. B. Geller, “Screened exchange LDA determination of the ground and excited state properties of thermoelectrics: Bi_{2}Te_{3},” Physical Review B, vol. 72, no. 3, pp. 2–5, 2005. View at: Publisher Site  Google Scholar
 G. Wang and T. Cagin, “Electronic structure of the thermoelectric materials Bi_{2}Te_{3} and Sb_{2}Te_{3} from firstprinciples calculations,” Physical Review B, vol. 76, no. 7, article 075201, 2007. View at: Publisher Site  Google Scholar
 J. M. Crowley, J. TahirKheli, and W. A. Goddard III, “Accurate Ab initio quantum mechanics simulations of Bi_{2}Se_{3} and Bi_{2}Te_{3} topological insulator surfaces,” Journal of Physical Chemistry Letters, vol. 6, no. 19, pp. 3792–3796, 2015. View at: Publisher Site  Google Scholar
 B. Y. Yavorsky, N. F. Hinsche, I. Mertig, and P. Zahn, “Electronic structure and transport anisotropy of Bi_{2}Te_{3} and Sb_{2},” Physical Review B, vol. 84, no. 16, article 165208, 2011. View at: Publisher Site  Google Scholar
 T. J. Scheidemantel, C. AmbroschDraxl, T. Thonhauser, J. V. Badding, and J. O. Sofo, “Transport coefficients from firstprinciples calculations,” Physical Review B, vol. 68, no. 12, article 125210, 2003. View at: Publisher Site  Google Scholar
 P. Larson, S. D. Mahanti, and M. G. Kanatzidis, “Electronic structure and transport _{ofBi2Te3andBaBiTe3},” Physical Review B, vol. 61, no. 12, pp. 8162–8171, 2000. View at: Publisher Site  Google Scholar
 G. A. Thomas, D. H. Rapkine, R. B. Van Dover et al., “Large electronicdensity increase on cooling a layered metal: Doped Bi_{2}Te_{3},” Physical Review B, vol. 46, no. 3, pp. 1553–1556, 1992. View at: Publisher Site  Google Scholar
 H. Lee and O. V. Yazyev, “Interplay between spin–orbit coupling and crystalfield effect in topological insulators,” Journal of Physics: Condensed Matter, vol. 27, no. 28, article 285801, 2015. View at: Publisher Site  Google Scholar
 W. G. Zeier, A. Zevalkink, Z. M. Gibbs, G. Hautier, M. G. Kanatzidis, and G. J. Snyder, “Thinking like a chemist: intuition in thermoelectric materials,” Angewandte Chemie International Edition, vol. 55, no. 24, pp. 6826–6841, 2016. View at: Publisher Site  Google Scholar
 G. R. Hyde, H. A. Beale, I. L. Spain, and J. A. Woollam, “Electronic properties of Bi_{2}Se_{3} crystals,” Journal of Physics and Chemistry of Solids, vol. 35, no. 12, pp. 1719–1728, 1974. View at: Publisher Site  Google Scholar
 A. F. May and G. J. Snyder, “Introduction to modeling thermoelectric transport at high temperatures,” in Materials, preparation, and characterization in thermoelectrics: Its Energy Harvesting, p. 11, CRC press, 2012. View at: Publisher Site  Google Scholar
 P. P. Konstantinov, L. V. Prokofieva, M. I. Fedorov et al., “Transport coefficients of nBi_{2}Te_{2.7}Se_{0.3} in a twoband model of the electron spectrum,” Semiconductors, vol. 39, no. 9, p. 1023, 2005. View at: Publisher Site  Google Scholar
 E. Bellotti, F. Bertazzi, and M. Goano, “Alloy scattering in AlGaN and InGaN: a numerical study,” Journal of Applied Physics, vol. 101, no. 12, article 123706, 2007. View at: Publisher Site  Google Scholar
 J. Bardeen and W. Shockley, “Deformation potentials and mobilities in nonpolar crystals,” Physics Review, vol. 80, no. 1, pp. 72–80, 1950. View at: Publisher Site  Google Scholar
 E. S. Toberer, A. Zevalkink, and G. J. Snyder, “Phonon engineering through crystal chemistry,” Journal of Materials Chemistry, vol. 21, no. 40, article 15843, 2011. View at: Publisher Site  Google Scholar
 H.S. Kim, S. D. Kang, Y. Tang, R. Hanus, and G. J. Snyder, “Dislocation strain as the mechanism of phonon scattering at grain boundaries,” Materials Horizons, vol. 3, no. 3, pp. 234–240, 2016. View at: Publisher Site  Google Scholar
 J. Callaway and H. C. von Baeyer, “Effect of point imperfections on lattice thermal conductivity,” Physics Review, vol. 120, no. 4, pp. 1149–1154, 1960. View at: Publisher Site  Google Scholar
 P. G. Klemens, “Thermal resistance due to point defects at high temperatures,” Physics Review, vol. 119, no. 2, pp. 507–509, 1960. View at: Publisher Site  Google Scholar
 E. F. Steigmeier and B. Abeles, “Scattering of phonons by electrons in germaniumsilicon alloys,” Physics Review, vol. 136, no. 4A, pp. A1149–A1155, 1964. View at: Publisher Site  Google Scholar
 J. Yang, G. P. Meisner, and L. Chen, “Strain field fluctuation effects on lattice thermal conductivity of ZrNiSnbased thermoelectric compounds,” Applied Physics Letters, vol. 85, no. 7, pp. 1140–1142, 2004. View at: Publisher Site  Google Scholar
 Y. Wu, P. Nan, Z. Chen et al., “Manipulation of band degeneracy and lattice strain for extraordinary PbTe thermoelectrics,” Research, vol. 2020, article 8151059, 12 pages, 2020. View at: Publisher Site  Google Scholar
 H. J. Goldsmid, Introduction to Thermoelectricity, Second., vol. 121, Springer Berlin Heidelberg, Berlin, Heidelberg, 2016. View at: Publisher Site
 A. Jain, S. P. Ong, G. Hautier et al., “Commentary: the Materials Project: a materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, article 011002, 2013. View at: Publisher Site  Google Scholar
 G. Kresse and J. Hafner, “Ab initiomolecular dynamics for liquid metals,” Physical Review B, vol. 47, no. 1, pp. 558–561, 1993. View at: Publisher Site  Google Scholar
 G. Kresse and J. Furthmüller, “Efficient iterative schemes forab initiototalenergy calculations using a planewave basis set,” Physical Review B, vol. 54, no. 16, pp. 11169–11186, 1996. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996. View at: Publisher Site  Google Scholar
 P. E. Blöchl, “Projector augmentedwave method,” Physical Review B, vol. 50, no. 24, pp. 17953–17979, 1994. View at: Publisher Site  Google Scholar
 G. Kresse and D. Joubert, “From ultrasoft pseudopotentials to the projector augmentedwave method,” Physical Review B, vol. 59, no. 3, pp. 1758–1775, 1999. View at: Publisher Site  Google Scholar
 P. V. C. Medeiros, S. Stafström, and J. Björk, “Effects of extrinsic and intrinsic perturbations on the electronic structure of graphene: retaining an effective primitive cell band structure by band unfolding,” Physical Review B, vol. 89, no. 4, article 041407, 2014. View at: Publisher Site  Google Scholar
 P. V. C. Medeiros, S. S. Tsirkin, S. Stafström, and J. Björk, “Unfolding spinor wave functions and expectation values of general operators: introducing the unfoldingdensity operator,” Physical Review B, vol. 91, no. 4, article 041116, 2015. View at: Publisher Site  Google Scholar
 S. P. Ong, W. D. Richards, A. Jain et al., “Python Materials Genomics (pymatgen): a robust, opensource python library for materials analysis,” Computational Materials Science, vol. 68, pp. 314–319, 2013. View at: Publisher Site  Google Scholar
 C. H. Champness, P. T. Chiang, and P. Parekh, “Thermoelectric properties of Bi_{2}Te_{3}–Sb_{2}Te_{3} alloys,” Canadian Journal of Physics, vol. 43, no. 4, pp. 653–669, 1965. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Ian T. Witting et al. Exclusive Licensee Science and Technology Review Publishing House. Distributed under a Creative Commons Attribution License (CC BY 4.0).