Scattering Mechanisms and Suppression of Bipolar Diffusion Effect in Bi2Te2.85Se0.15Ix Compounds

We investigated the anisotropic thermoelectric properties of the Bi2Te2.85Se0.15Ix (x = 0.0, 0.1, 0.3, 0.5 mol.%) compounds, synthesized by ball-milling and hot-press sintering. The electrical conductivities of the Bi2Te2.85Se0.15Ix were significantly improved by the increase of carrier concentration. The dominant electronic scattering mechanism was changed from the mixed (T ≤ 400 K) and ionization scattering (T ≥ 420 K) for pristine compound (x = 0.0) to the acoustic phonon scattering by the iodine doping. The Hall mobility was also enhanced with the increasing carrier concentration. The enhancement of Hall mobility was caused by the increase of the mean free path of the carrier from 10.8 to 17.7 nm by iodine doping, which was attributed to the reduction of point defects without the meaningful change of bandgap energy. From the electron diffraction patterns, a lattice distortion was observed in the iodine doped compounds. The modulation vector due to lattice distortion increased with increasing iodine concentration, indicating the shorter range lattice distortion in real space for the higher iodine concentration. The bipolar thermal conductivity was suppressed, and the effective masses were increased by iodine doping. It suggests that the iodine doping minimizes the ionization scattering giving rise to the suppression of the bipolar diffusion effect, due to the prohibition of the BiTe1 antisite defect, and induces the lattice distortion which decreases lattice thermal conductivity, resulting in the enhancement of thermoelectric performance.

Thermoelectric properties and X-ray diffraction (XRD) patterns of the hot-pressed samples were performed with the parallel (Pa) and perpendicular (Pe) directions to the press direction. The XRD patterns were obtained using the Cu-κα radiation (D8 advance, Bruker, Germany). The temperature-dependent electrical conductivity (σ) and Seebeck coefficient ( S ) were simultaneously measured under the helium atmosphere by a four-point probe method using a thermoelectric properties measurement system (ZEM-3, ULVAC-RIKO, Yokohama, Kanagawa, Japan).
The Hall carrier concentrations (n H ) and the Hall coefficient (R H ) were obtained by the relations of n H = −1/(R H e), where R H = ρ xy /H is the Hall coefficient, e is the elementary charge, ρ xy is the Hall resistivity, H is the applied magnetic field. The Hall resistivity (ρ xy ) was measured by the five-probe contact method under sweeping magnetic fields from −5 to 5 T using the physical property measurement system (PPMS Dynacool 14 T, Quantum Design, San Diego, CA, USA). The Hall mobility was obtained from the relation µ H = 1/n H eρ.
The total thermal conductivity (κ) was obtained by the relation of κ = ρ s λC p , where ρ s , λ, and C p are sample density, thermal diffusivity, and specific heat, respectively. The thermal diffusivity (λ) was measured by a laser flash method (LFA-457, NETZSCH, Germany). The sample density was obtained using the sample mass and dimensions. The specific heat (C p ) was estimated from high temperature fitting from the measurement data by the physical property measurement system (PPMS Dynacool 14 T, Quantum Design, San Diego, CA, USA). Figure 1a presents the XRD patterns of the hot-pressed Bi 2 Te 2.85 Se 0.15 I x (x = 0.0, 0.1, 0.3, 0.5 mol.%) compounds for two different planes of parallel (Pa, solid line) and perpendicular (Pe, dotted line) directions to the press direction, as indicated in the right inset. All the peaks were normalized with the highest (015) peaks. The XRD peaks were indexed with the Bi 2 Te 3 structure (Rhombohedral, space group No.166) with a small change of lattice parameter. specific heat (Cp) was estimated from high temperature fitting from the measurement data by the physical property measurement system (PPMS Dynacool 14 T, Quantum Design, San Diego, CA, USA). Figure 1a presents the XRD patterns of the hot-pressed Bi2Te2.85Se0.15Ix (x = 0.0, 0.1, 0.3, 0.5 mol.%) compounds for two different planes of parallel (Pa, solid line) and perpendicular (Pe, dotted line) directions to the press direction, as indicated in the right inset. All the peaks were normalized with the highest (015) peaks. The XRD peaks were indexed with the Bi2Te3 structure (Rhombohedral, space group No.166) with a small change of lattice parameter. As presented in the left-upper inset, the XRD patterns along the parallel direction showed the enhanced (00l) peaks, indicating that the bismuth telluride's quintuple layers were stacked along the parallel direction to the press direction, which was similar to those of the hot-deformed Bi2−xSbxTe2.3Se0.7 [13] and the well-aligned Bi2Te3−xSex nanocrystals [11]. The preferred orientation in polycrystal is frequently found in layer structured materials [23,28]. The anisotropic texture implies the anisotropic thermoelectric properties.

Results and Discussion
The c-axis lattice parameters decreased with increasing iodine doping concentration, while the a-axis lattice parameter was not sensitive to the iodine concentration, as shown in Figure 1b, which was due to the smaller atomic radius of I (0.133 nm) than Te (0.137 nm). The decrease of lattice parameters was found in the Bi2Te3−xIx compounds [29]. The formation energy calculation shows that the iodine can be substituted at the Te(2) sites [19]. Therefore, the reduced lattice parameters of the Bi2Te2.85Se0.15Ix (x = 0.1, 0.3, and 0.5 %) compounds can be understood by the iodine substitution at the Te sites. Figure 2a,b presents the temperature-dependent electrical conductivities σ(T) of the Bi2Te2.85Se0.15Ix (x = 0, 0.1, 0.3, 0.5 mol.%) compounds. The anisotropic measurement of σ(T) for x = 0.0 had two distinctive regions, as shown in Figure 2a; one was temperatureinsensitive (T ≤ 400 K) and the other was a region linearly increasing with temperature (T ≥ 420 K). The different temperature-dependency of σ(T) can be understood by the distinctively different scattering mechanism. The electrical conductivity is written as σ = neμ, where mobility is implicitly represented by the effective mass m* and average relaxation time <τ> by μ = e<τ>/m*. The relaxation time is a function of electronic energy E and temperature T by the relation; = ( / ) , where r is the scattering factor such that r = −1/2 for acoustic phonon scattering, r = 0 for a neutral impurity scattering, r = 1/2 for an optical phonon scattering, r = 1 for a mixed scattering, and r = 3/2 for an ionized As presented in the left-upper inset, the XRD patterns along the parallel direction showed the enhanced (00l) peaks, indicating that the bismuth telluride's quintuple layers were stacked along the parallel direction to the press direction, which was similar to those of the hot-deformed Bi 2−x Sb x Te 2.3 Se 0.7 [13] and the well-aligned Bi 2 Te 3−x Se x nanocrystals [11]. The preferred orientation in polycrystal is frequently found in layer structured materials [23,28]. The anisotropic texture implies the anisotropic thermoelectric properties.
The c-axis lattice parameters decreased with increasing iodine doping concentration, while the a-axis lattice parameter was not sensitive to the iodine concentration, as shown in Figure 1b, which was due to the smaller atomic radius of I (0.133 nm) than Te (0.137 nm). The decrease of lattice parameters was found in the Bi 2 Te 3−x I x compounds [29]. The formation energy calculation shows that the iodine can be substituted at the Te(2) sites [19]. Therefore, the reduced lattice parameters of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.1, 0.3, and 0.5%) compounds can be understood by the iodine substitution at the Te sites. Figure 2a,b presents the temperature-dependent electrical conductivities σ(T) of the Bi 2 Te 2.85 Se 0.15 I x (x = 0, 0.1, 0.3, 0.5 mol.%) compounds. The anisotropic measurement of σ(T) for x = 0.0 had two distinctive regions, as shown in Figure 2a; one was temperatureinsensitive (T ≤ 400 K) and the other was a region linearly increasing with temperature (T ≥ 420 K). The different temperature-dependency of σ(T) can be understood by the distinctively different scattering mechanism. The electrical conductivity is written as σ = neµ, where mobility is implicitly represented by the effective mass m* and average relaxation time <τ> by µ = e<τ>/m*. The relaxation time is a function of electronic energy E and temperature T by the relation; τ = τ 0 (E/k B T) r , where r is the scattering factor such that r = −1/2 for acoustic phonon scattering, r = 0 for a neutral impurity scattering, r = 1/2 for an optical phonon scattering, r = 1 for a mixed scattering, and r = 3/2 for an ionized impurity scattering [30][31][32][33]. The σ(T) of the pristine Bi 2 Te 2.85 Se 0.15 compound is attributed to the mixed scattering at a lower temperature than 400 K and the ionized impurity scattering (σ ∝ T 1.5 ) at high temperatures (T ≥ 420 K). The temperature-exponent behavior of Hall mobility may also follow the one of electrical conductivity because the carrier concentration is not sensitive to temperature in a degenerated semiconductor. So, it is reasonable to interpret the scattering exponent behavior from the temperature-dependent electrical conductivity. electrical conductivity.
The σ(T) of Bi2Te2.85Se0.15Ix (0.1 ≤ x ≤ 0.5 mol.%) exhibited a metallic or highly degenerated semiconducting behavior, as shown in Figure 2b. The electrical conductivity of the x = 0.5 mol.% sample (1717 S/cm at 323 K) was comparable with the one of the hotpressed Bi2Te2.85Se0.15I0.005 compound (about 1600 S/cm at 323 K) [17]. The electrical conductivities of the Pe-direction were higher than the Pa-direction values due to the preferred orientation of the layer structure of the bismuth telluride. From the fitting of the scattering exponent, the main scattering mechanism of the Bi2Te2.85Se0.15Ix compounds could be regarded as the acoustic phonon scattering [30,34]   The room-temperature electrical conductivities (left axis) were increased with increasing the iodine doping concentration for parallel (open black square) and perpendicular (closed black square) direction for the press direction, as shown in Figure  2c. The enhancement of electrical conductivity by iodine doping was mainly affected by The σ(T) of Bi 2 Te 2.85 Se 0.15 I x (0.1 ≤ x ≤ 0.5 mol.%) exhibited a metallic or highly degenerated semiconducting behavior, as shown in Figure 2b. The electrical conductivity of the x = 0.5 mol.% sample (1717 S/cm at 323 K) was comparable with the one of the hot-pressed Bi 2 Te 2.85 Se 0.15 I 0.005 compound (about 1600 S/cm at 323 K) [17]. The electrical conductivities of the Pe-direction were higher than the Pa-direction values due to the preferred orientation of the layer structure of the bismuth telluride. From the fitting of the scattering exponent, the main scattering mechanism of the Bi 2 Te 2.85 Se 0.15 I x compounds could be regarded as the acoustic phonon scattering [30,34] because the scattering exponents on the I-doped compounds become close to the σ ∝ T −1.5 , where r = −1.25 (x = 0.1 mol.%), −1.31 (x = 0.3 mol.%), and −1.36 (x = 0.5 mol.%).
The room-temperature electrical conductivities (left axis) were increased with increasing the iodine doping concentration for parallel (open black square) and perpendicular (closed black square) direction for the press direction, as shown in Figure 2c. The enhancement of electrical conductivity by iodine doping was mainly affected by the enhanced Hall carrier concentrations (right axis, red closed circle). From the electrical conductivity and Hall carrier concentration, the Hall mobilities µ H = 1/(neρ) were obtained. The Hall where m* is the effective mass of the carrier, r H is the Hall factor, η = E F /k B T is the reduced Fermi energy, and F n (η) is the n-th order Fermi integral given by [31,35] The effective masses of the carrier can be obtained from the Hall carrier concentration, while the reduced Fermi energy is obtained from the Seebeck coefficient as the following equation [32,35].
In the case of the acoustic phonon dominant scattering (r = −1/2), the Hall mobilities decreased with increasing Hall carrier concentration, as presented in Figure 2d. When we assume the mixed scattering (r = 1), the Hall mobility significantly increased with increasing Hall carrier concentration for perpendicular (solid black line) and parallel (black dotted line) directions, while the Hall mobilities in acoustic phonon scattering decreased with increasing Hall concentration, as shown in Figure 2d. The Hall mobilities were enhanced with increasing the I-doping concentration. The Hall mobility versus Hall carrier concentration showed that the acoustic phonon scattering was a dominant scattering mechanism. Figure 3a,b show the temperature-dependent negative Seebeck coefficients S(T) for the pristine (a) and I-doped Bi 2 Te 2.85 Se 0.15 I x (x = 0.1, 0.3, 0.5 mol.%) compounds (b). The negative Seebeck coefficients indicated the electronic transport of the carrier. The Seebeck coefficients of the pristine sample (x = 0.0) for the perpendicular (−216 µV/K) and parallel (−215 µV/K) directions were similar to the hot-pressed Bi 2 Te 2.85 Se 0.15 compounds (about −220 µV/K) at room temperature [36]. The S(T) of the parent compound was well described by the mixed scattering exponent (r = 1, blue dashed line) below 425 K, while the ionization scattering (r = 3/2) was dominant at a high-temperature range (T ≥ 425 K). This scattering mechanism was also consistent with the S(T) along the Pa-direction, as presented in the inset of Figure 3a.
The temperature-dependent behavior of the Seebeck coefficient on the I-doped compounds (0.1 ≤ x ≤ 0.5 mol.%) showed a different behavior with those of the pristine compound, as shown in Figure 3b. The symbols were experimental S(T) data for Pe-(closed symbols) and Pa-(open symbols) directions. The solid and dashed lines were a theoretical fitting with the experiment in terms of acoustic phonon scattering (r = −1/2). Therefore, the dominant scattering mechanism was changed from mixed (T ≤ 425 K) and ionization scattering (T ≥ 425 K) in the pristine compound to the acoustic phonon scattering over a wide temperature range in the iodine doped compounds. temperature, we used the corresponding scattering exponent within a single parabolic band assumption. The effective masses of the carriers in the I-doped samples were significantly increased (m* = 0.73 ~ 1.10, r = −1/2) as compared with the pristine sample (m* = 0.31, r = 1). The carrier effective mass could be enhanced by increasing the carrier scattering with the nanoparticles [37], resonant level formation [38], Rashba band splitting [19,39], nonparabolic band [40], the increase of band valley degeneracies [41,42] and the charge density wave (CDW) formation [43], etc.
where EF is the Fermi energy. The carrier mean free paths λ of the pristine sample Bi2Te2.85Se0.15 are about 10.8 nm (Pe-direction) and 8.0 nm Pa-direction, which are comparable with the zone melted n-type Bi2Te2.79Se0.21 (~ 7 nm) [45] and the hot-pressed CuI-doped Bi2Te2.7Se0.3 (10.82 nm) [10]. The mean free paths of the I-doped samples increased with increasing I-doping concentration up to 17.7 nm (Pe) and 13.0 nm (Pa). The larger λ in perpendicular direction than those in the parallel direction was due to the anisotropic electronic transport along with the in-plane electronic conductivity (perpendicular to the press direction).  Figure 3c shows the negative Seebeck coefficient −S as a function of Hall carrier concentration n H at room temperature. Because the mixed (r = 1) and acoustic scatterings (r = −1/2) are dominant in parent and iodine-doped compounds, respectively, at room temperature, we used the corresponding scattering exponent within a single parabolic band assumption. The effective masses of the carriers in the I-doped samples were significantly increased (m* = 0.73~1.10, r = −1/2) as compared with the pristine sample (m* = 0.31, r = 1). The carrier effective mass could be enhanced by increasing the carrier scattering with the nanoparticles [37], resonant level formation [38], Rashba band splitting [19,39], nonparabolic band [40], the increase of band valley degeneracies [41,42] and the charge density wave (CDW) formation [43], etc.
From the carrier mobility µ and the carrier effective mass m*, the mean free paths λ of the carrier can be calculated by the following equation, as shown in Figure 3d (left axis, black squares) [44,45].
where E F is the Fermi energy. The carrier mean free paths λ of the pristine sample Bi 2 Te 2.85 Se 0.15 are about 10.8 nm (Pe-direction) and 8.0 nm Pa-direction, which are comparable with the zone melted n-type Bi 2 Te 2.79 Se 0.21 (~7 nm) [45] and the hot-pressed CuI-doped Bi 2 Te 2.7 Se 0.3 (10.82 nm) [10]. The mean free paths of the I-doped samples increased with increasing I-doping concentration up to 17.7 nm (Pe) and 13.0 nm (Pa). The larger λ in perpendicular direction than those in the parallel direction was due to the anisotropic electronic transport along with the in-plane electronic conductivity (perpendicular to the press direction). The enhanced Hall mobilities of Bi 2 Te 2.85 Se 0.15 I x were affected by the increase of carrier mean free path rather than the Fermi energy E F (right axis, red circles), as shown in Figure 3d. Because the E F subsequently decreased with the substitution of iodine concentration for x = 0.1 mol.% with the small increase of E F for high doping concentration. The carrier mean free path (λ) can be increased by the preferred crystallographic orientation (texturing effect) in bismuth tellurides [10] and by the reduction of defects [46]. Since the texturing effect with increasing I-doping concentration was not significantly observed in the XRD patterns, the increase of carrier mean free path could be attributed to the reduction of defects. The Bi 2 (Te, Se) 3 compounds included the intrinsic point defects such as the Se vacancy. The Se-deficiency generated the Bi Te antisite defect at the Te(Se) sites [46]. The iodine atoms substituted at the Te(Se) sites, resulting in the inhibition of the vacancy and antisite defects [19]. Therefore, the enhanced mean free path of the Bi 2 Te 2.85 Se 0.15 I x compounds was attributed to the reduction of point defects.
The temperature-dependent total thermal conductivities κ total (T) of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.0, 0.1, 0.3, 0.5 mol.%) compounds are shown in Figure 4a. The κ total (T) of the pristine sample (x = 0.0) was increased with increasing temperature (black squares), which is typical behavior of bipolar diffusion effect [47]. The κ total (T) of the I-doped samples (x = 0.1, 0.3, and 0.5 mol.%) decreased with increasing temperature up to the 420 K and then increased again at high temperatures, implying the suppression of the bipolar diffusion effect by iodine doping. The total thermal conductivities of the Pe-direction were higher than those of the Pa-direction by the layer structure of the bismuth telluride [20,23]. The enhanced Hall mobilities of Bi2Te2.85Se0.15Ix were affected by the increase of carrier mean free path rather than the Fermi energy EF (right axis, red circles), as shown in Figure  3d. Because the EF subsequently decreased with the substitution of iodine concentration for x = 0.1 mol.% with the small increase of EF for high doping concentration. The carrier mean free path (λ) can be increased by the preferred crystallographic orientation (texturing effect) in bismuth tellurides [10] and by the reduction of defects [46]. Since the texturing effect with increasing I-doping concentration was not significantly observed in the XRD patterns, the increase of carrier mean free path could be attributed to the reduction of defects. The Bi2(Te, Se)3 compounds included the intrinsic point defects such as the Se vacancy. The Se-deficiency generated the BiTe antisite defect at the Te(Se) sites [46]. The iodine atoms substituted at the Te(Se) sites, resulting in the inhibition of the vacancy and antisite defects [19]. Therefore, the enhanced mean free path of the Bi2Te2.85Se0.15Ix compounds was attributed to the reduction of point defects.
The temperature-dependent total thermal conductivities κtotal (T) of the Bi2Te2.85Se0.15Ix (x = 0.0, 0.1, 0.3, 0.5 mol.%) compounds are shown in Figure 4a. The κtotal(T) of the pristine sample (x = 0.0) was increased with increasing temperature (black squares), which is typical behavior of bipolar diffusion effect [47]. The κtotal(T) of the I-doped samples (x = 0.1, 0.3, and 0.5 mol.%) decreased with increasing temperature up to the 420 K and then increased again at high temperatures, implying the suppression of the bipolar diffusion effect by iodine doping. The total thermal conductivities of the Pe-direction were higher than those of the Pa-direction by the layer structure of the bismuth telluride [20,23].  However, in contrast with that, the Pe-direction electrical conductivities and Seebeck coefficients of the x = 0.5 mol.% sample were comparable with the data of the Bi 2 Te 2.85 Se 0.15 I 0.005 compound [17], the thermal conductivities of the I-doped sample (0.5 mol.%) [17] were closer to the Pa-direction thermal conductivities rather than Pedirection values of the x = 0.5 mol.% sample, which was not mention the measurement direction of the thermoelectric properties. The measurement direction matching of the power factor and thermal conductivity are necessary for the practical and correct combination of the thermoelectric properties [6]. Therefore, we believe that the thermal conductivities on this I-doped Bi 2 Te 2.85 Se 0.15 are more reliable.
The electronic thermal conductivity can be calculated by the Wiedemann-Franz law (κ el = L 0 σT), where L 0 , σ, and T are the Lorenz number, electrical conductivity, and absolute temperature, respectively. In the case of simple metals, the Lorenz number is L 0 = (π 2 /3)·(k B /e) 2 = 2.45 × 10 −8 W·Ω·K −2 . Since the Lorenz number is incorrect in a correlated metal or a degenerated semiconductor, the Lorenz number is calculated by using the following equation [31]. Figure 4b presents the calculated temperature-dependent Lorenz numbers of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.0, 0.1, 0.3, and 0.5 mol.%) compounds. When we calculated the Lorenz number for the parent compound (x = 0.0) by mixed scattering (r = 1, blue star) and ionization scattering (r = 3/2, red star), the Lorenz numbers were about 2.4~2.5 × 10 −8 W·Ω·K −2 . On the other hand, for acoustic phonon scattering (r = − 1/2), the Lorenz numbers were in the range of (1.6~2.1) × 10 −8 W·Ω·K −2 . When we subtracted the electronic thermal conductivity κ el from total thermal conductivity κ total , the combined contribution of lattice κ L and bipolar thermal conductivities κ b were obtained, as shown in Figure 4c. All the samples have a bipolar diffusion effect, which is increased with increasing temperature at high temperatures (T ≥ 425 K). The bipolar diffusion effect was significantly suppressed in the I-doped Bi 2 Te 2.85 Se 0.15 I x compounds (x = 0.1, 0.3, and 0.5 mol.%) even a small amount of the iodine doping. It is generally known that bismuth tellurides have a bipolar effect due to the small band gap. [48,49] The co-excitation of electron and hole in the bipolar effect decreases the Seebeck coefficient, resulting in the deterioration of thermoelectric performance. The suppression of the bipolar diffusion effect has been found in the I-doped [25,49], Cu-doped [8], Ca-doped [48], CuI-doped bismuth tellurides [19], Bi 2 Te 3−x Se x [4,25] as well as the case of an increase of energy band gap [4,25].
The pristine sample (x = 0.0) showed a semiconducting behavior in electrical resistivity above 400 K. From the Arrhenius equation σ = σ 0 e −E a /k B T , where σ 0 is a preexponential factor, E a is the activation energy, and k B is Boltzmann constant [50], the activation energy E a and band gap energy E g are obtained by E a~0 .076 eV and E g~0 .152 eV (not shown), which is estimated by the relation of the band gap energy E a~Eg /2. In addition, the thermal band gap energy of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.1, 0.3, and 0.5 mol.%) compounds is obtained by about 0.14 eV from the E g = 2eS max T max relation [25], where S max is the maximum Seebeck coefficient at the temperature (T max ). The estimated band gap energies of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.0, 0.1, 0.3, and 0.5 mol.%) compounds are comparable with the optical band gap (E g = 0.14~0.18 eV) and thermal band gap (E g~0 .154 eV) of the Idoped Bi 2 (Te 1−x Se x ) [25] and the optical bandgap (0.14 eV) of Bi 2 Te 3−x I x . [49] The obtained band gap energies showed that the band gap energy was not changed significantly by iodine doping.
The suppression of the bipolar diffusion effect by iodine doping can be described by the inhibition of the vacancy and antisite defects due to stable substitution at the Te-site. The formation energy calculation results showed that the Bi Te1 anti-site defects were the most energetically favorable defect in Bi-rich composition. In this case, the minority hole carriers could be generated by the Bi Te1 antisite defects in the n-type bismuth telluride [51]. The pristine bismuth telluride can be regarded as the Bi-rich phase due to the Te evaporation during synthesis, which is well known as the origin of the native defects [4,36]. For example, the Se-deficient compound Bi 2 Te 2.3 Se 0.7−x clearly showed that the bipolar effects became significant with increasing Se deficiencies [12]. The iodine substitution caused the suppression of the bipolar effect in the I-doped sample (x = 0.1, 0.3, and 0.5 mol.%) at the Te(2) sites [19]. Therefore, the reduction of point defects in the Bi 2 Te 2.85 Se 0.15 I x is supported by the enhanced mean free path of the carriers.
It is noteworthy that the lattice and bipolar thermal conductivities κ L + κ b of the Bi 2 Te 2.85 Se 0.15 I x compounds (x = 0.0, 0.1, 0.3, and 0.5 mol.%) clearly showed that the iodine doping decreased the lattice and bipolar thermal conductivities significantly, as shown in Figure 4c,d. The lattice thermal conductivity of the x = 0.1 mol.% showed the unusual anisotropic lattice thermal conductivity, which was lower lattice thermal conductivity along the Pe-direction (in-plane dominant) than that of the Pa-direction (out-of-plane dominant) values above 473 K.
The decreased lattice and bipolar thermal conductivity and the unusual anisotropic lattice thermal conductivity can be explained by the lattice distortion (or CDW). Previously, we found the lattice distortion in the hot-deformed CuI-doped Bi 2 Te 2.7 Se 0.3 [20] and CDWlike behavior in the CuI-doped Bi 2 Te 2.1 Se 0.9 [27]. The lattice distortion is effective to decrease lattice thermal conductivity because the lattice thermal conductivity of the in-plane direction (strong covalent bonding layer) can be smaller than the out-of-plane direction lattice thermal conductivity (weak van der Waals bonding) like a charge density wave (CDW) effect [28,52,53].
The lattice distortion of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.1, 0.3, and 0.5 mol.%) compounds were found in the high-resolution transmission electron microscopy (HR-TEM) and electron diffraction (ED) patterns, as shown in Figure 5. The ED pattern of the parent compound showed clear spots along the [111] axis, as presented in the inset of Figure 5a. On the other hand, there were lattice distortion peaks for the iodine doped compounds. The modulation vectors were indicated from the (000) point, marked by the red arrow. The modulation vectors increased with increasing iodine doping concentration, as shown in Figure 5b-d, indicating the shorter range lattice distortion in real space for the higher I-doping concentration. The periodic lattice modulation was also found in the CuI-doped Bi 2 Te 2.7 Se 0.3 [20] and the CuI-doped Bi 2 Te 2.1 Se 0.9 [27]. The lattice distortion by the iodine doping significantly decreased the lattice thermal conductivity.
Temperature-dependent power factors S 2 σ(T) of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.0, 0.1, 0.3, 0.5 mol.%) compounds are presented in Figure 6a. The room-temperature power factor of the pristine sample was about 1.4 mW·m −1 ·K −2 for the Pe-direction (1.0 mW·m −1 ·K −2 for the Pa-direction). The power factors of the I-doped Bi 2 Te 2.85 Se 0.15 I x (x = 0.1, 0.3 and 0.5 mol.%) compounds were increased owing to the enhancement of electrical conductivity as compared with the pristine sample. Because of the anisotropic electrical conductivity, the power factor of the Pe-direction was higher than those of the Pa-direction. Figure 6b presents that the power factor reached 2.7 mW·m −1 ·K −2 at 300 K for x = 0.1 mol.% for the Pe-direction. The power factor of the pristine sample (1.13 mW·m −1 ·K −2 at 323 K, x = 0.0) was much lower than the one of Bi 2 Te 2.85 Se 0.15 (~2.10 mW·m −1 ·K −2 at 323 K) [17]. compounds were increased owing to the enhancement of electrical conductivity as compared with the pristine sample. Because of the anisotropic electrical conductivity, the power factor of the Pe-direction was higher than those of the Pa-direction. Figure 6b presents that the power factor reached 2.7 mW·m −1 ·K −2 at 300 K for x = 0.1 mol.% for the Pe-direction. The power factor of the pristine sample (1.13 mW·m −1 ·K −2 at 323 K, x = 0.0) was much lower than the one of Bi2Te2.85Se0.15 (~2.10 mW·m −1 ·K −2 at 323 K) [17].    3, and 0.5 mol.%) were higher than those of the pristine sample (x = 0.0) by the enhanced power factor and the decreased lattice and bipolar thermal conductivity. The zT values of the Pe-direction were higher than the Pa-direction from the anisotropic electrical conductivity. The zT value of the x = 0.1 mol.% was reached 0.72 at 423 K for the Pe-direction. Therefore, the Iodine doping suppresses the bipolar diffusion effect and increases the electrical conductivity, resulting in the enhancement of the power factor.
A recent study suggested how to find the maximum power factor without experimental optimization by defining the electronic quality factor. The electronic quality factor is defined by [54] where S r is the scaled Seebeck coefficient defined by |S|e/k B . Figure 7 shows the power factor S 2 σ with absolute Seebeck coefficient of the

Conclusions
In summary, we investigated the anisotropic thermoelectric properties of the Bi2Te2.85Se0.15Ix (x = 0.0, 0.1, 0.3, and 0.5 mol.%) compounds, synthesized by ball-milling and hot press sintering. The anisotropic thermoelectric properties were attributed to the structurally preferred orientation of the bismuth telluride along the (00l) direction with

Conclusions
In summary, we investigated the anisotropic thermoelectric properties of the Bi 2 Te 2.85 Se 0.15 I x (x = 0.0, 0.1, 0.3, and 0.5 mol.%) compounds, synthesized by ball-milling and hot press sintering. The anisotropic thermoelectric properties were attributed to the structurally preferred orientation of the bismuth telluride along the (00l) direction with the press direction. The electrical conductivities of the I-doped samples were mainly enhanced by the increased carrier concentration. The Hall mobility was enhanced with increasing carrier concentration due to the increase of the mean free path from 10.8 to 17.7 nm by iodine doping. The enhanced mean free path was caused by the reduction of point defects. We found that the carrier scattering mechanism was changed from the mixed (T ≤ 400 K) and ionization scattering (T ≥ 420 K) for pristine compound (x = 0.0) to the acoustic phonon scattering for iodine doped compounds (x = 0.1~0.5 mol.%). The electron diffraction results clearly showed that the I-doping generated the periodic lattice distortion in the bismuth telluride. For increasing I-doping concentration, the modulation vector in momentum space was increased, in other words, the shorter range of lattice distortion for the higher iodine doping concentration. The periodic lattice distortion could increase the carrier effective mass as the CDW effect. The ionization scattering of the carrier should be minimized to suppress the bipolar diffusion effect in thermal conductivity due to the prohibition of vacancy and antisite defects due to stable substitution at the Te-site, while the lattice distortion decreases the lattice thermal conductivity by iodine doping. Because of the enhanced power factor and the decreased lattice and bipolar thermal conductivity by the iodine doping, we observed the enhancement of zT values from 0.31 (x = 0.0) to 0.56 (x = 0.1 mol.%) at 300 K (0.72 at 423 K for Pe-direction). This research strongly suggests that the minimization of ionization scattering suppresses the bipolar diffusion effect by hindering vacancy and defect states, and lattice distortion decreases lattice thermal conductivity, resulting in the enhancement of thermoelectric performance.