Thermoelectric Response in Single Quintuple Layer Bi2Te3 - ACS

Oct 5, 2016 - Because Bi2Te3 belongs to the most important thermoelectric materials, the successful exfoliation of a single quintuple layer has opened...
0 downloads 8 Views 2MB Size
Thermoelectric Response in Single Quintuple Layer Bi2Te3 S. Sharma and U. Schwingenschlögl* Physical Science and Engineering Division (PSE), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia ABSTRACT: Because Bi2Te3 belongs to the most important thermoelectric materials, the successful exfoliation of a single quintuple layer has opened access to an interesting two-dimensional material. For this reason, we study the thermoelectric properties of single quintuple layer Bi2Te3 by considering both the electron and phonon transport. On the basis of first-principles density functional theory, the electronic and phononic contributions are calculated by solving Boltzmann transport equations. The dependence of the lattice thermal conductivity on the phonon mean free path is evaluated along with the contributions of the acoustic and optical branches. We find that the thermoelectric response is significantly better for p- than for n-doping. By optimizing the carrier concentration, at 300 K, a ZT value of 0.77 is achieved, which increases to 2.42 at 700 K.

W

the c-axis, referred to as quintuple layers (QLs). Weak van der Waals (vdW) interaction between the QLs makes it possible to obtain single-QL samples. Mechanically exfoliated single-QL samples and ultrathin films show enhanced thermoelectric response.36,37 Density functional theory has been used to describe the behavior of single-QL Bi2Te3,38 and the thermal conductivities of perfect and nanoporous few-QL films have been studied by molecular dynamics simulations,39 finding that nanoporousity reduces the thermal conductivity. The thermoelectric parameters of single-QL Bi2Te3 have been calculated by the Landauer approach using ab initio band structure data.40 Employing the experimental value of κlat (1.5 W K−1 m−1), the authors have predicted an increase of ZT at room temperature by a factor of 10 as compared to bulk Bi2Te3. The thermoelectric behavior of single-QL Bi2Te3 also has been investigated by means of Boltzmann equations for the electronic and phononic transport, based on molecular dynamics simulations,41 predicting ZT = 0.6 at room temperature. Recently, the same approach has been used for studying the electronic transport in few-QL Bi2Te3 nanofilms,42 employing κlat values from ref 39. Single-QL Bi2Te3 (ZT ≈ 1) turns out to be superior to thicker films. All previous theoretical studies have predicted κlat using molecular dynamics simulations, though the accuracy depends critically on the accuracy of the interatomic potentials.30,39 The problem can be avoided in density functional theory, as applied in the present work in combination with Boltzmann transport theory. We obtain κlat by solving the Boltzmann transport

hile the consumption of energy is increasing, the amount of available natural resources is depleting quickly. Search for clean and renewable energy resources thus is becoming more and more critical. One very promising approach is the recovery of waste heat with the help of thermoelectric materials.1−6 The conversion efficiency of such materials depends on the dimensionless figure of merit, ZT = S2σT/κ, where S, σ, κ, and T are the Seebeck coefficient, electrical conductivity, thermal conductivity (sum of electronic (κele) and lattice (κlat) contributions), and temperature, respectively. Achieving a high ZT value requires maximizing the power factor S2σ and minimizing κ. However, the interdependence between these transport coefficients makes optimization a challenging task. Hicks and Dresselhaus have theoretically demonstrated that the reduction of the dimensionality of a material can enhance ZT strongly.7,8 In recent years, rapid progress in nanotechnology has provided a large variety of low-dimensional materials, such as graphene,9 phosphorene,10,11 silicene,12 stanene, 13 transition metal dichalcogenides, 14,15 and MXenes,16,17 which are currently attracting growing interest from the thermoelectric point of view. For example, monolayer transition metal dichalcogenides have been predicted to improve ZT strongly as compared to the corresponding bulk materials,14,15 and a superior thermoelectric response has been reported for SnSe nanosheets.18 The thermoelectric properties of Bi2Te3 and its alloys have been studied extensively both experimentally19−28 and theoretically.7,29−34 The material has a layered rhombohedral structure with lattice parameters a = 0.438 nm and c = 3.045 nm (space group R3̅m),35 consisting of layers of five atomic planes with the sequence Te(1)−Bi−Te(2)−Bi−Te(3) along © XXXX American Chemical Society

Received: July 20, 2016 Accepted: October 3, 2016

875

DOI: 10.1021/acsenergylett.6b00289 ACS Energy Lett. 2016, 1, 875−879

Letter

http://pubs.acs.org/journal/aelccp

Letter

ACS Energy Letters equation for phonons (without any assumption about the phonon lifetimes) iteratively by means of the ShengBTE code.43 This approach is parameter-free and requires only crystal structure information. It has been used successfully for predicting κlat of GeSb2Te4 (experimental value of 1.3 W K−1 m−1, theoretical value of 1.5 W K−1 m−1).44 For describing the electronic transport, we employ also Boltzmann transport theory. Different density functionals are compared, in particular, the Perdew−Burke−Ernzerhof (PBE) functional45 and the computationally expensive hybrid Heyd−Scuseria− Ernzerhof (HSE06) functional.46 We use density functional theory as implemented in the Vienna ab initio simulation package (projector augment wave approach).47 A vacuum slab of more than 15 Å thickness is used to build a two-dimensional model and to prevent artificial interaction because of the periodic boundary conditions. In order to account for the vdW interaction, the optB86b-vdW functional is employed for structure optimization48 because previous results on bulk Bi2Te3 indicate best performance under the vdW functionals with respect to the description of the lattice parameters.34 The energy cutoff of the plane wave expansion is set to 450 eV, and the Brillouin zone is sampled on a Monkhorst−Pack 11 × 11 × 1 k-mesh. Both the lattice constants and atomic positions are relaxed until the forces on the atoms have declined to 0.01 eV/Å and the total energy changed to 0.001 meV. The transport coefficients are obtained by Boltzmann transport theory within the constant relaxation time approximation. For this purpose, the electronic band structure is calculated at the PBE level on a 45 × 45 × 1 k-mesh (following ref 42, while we find that a 21 × 21 × 1 k-mesh is already well converged) and at the HSE06 level on a 21 × 21 × 1 k-mesh, taking into account the spin−orbit coupling. The effect of doping is studied within the rigid band approach as implemented in the BoltzTrap code.49 Moreover, the phonon dispersion relation is calculated using the phonopy code.50 For evaluating κlat, the required harmonic interatomic force constants are obtained by density functional perturbation theory51 and the third-order interatomic force constants by a finite difference method, using a 4 × 4 × 1 supercell and a 3 × 3 × 1 k-mesh. For the third-order force constants, the interaction up to fourth nearest neighbors is considered. Using both the harmonic and third-order force constants as input, the Boltzmann transport equation for phonons is solved selfconsistently by means of the ShengBTE code43 in order to determine the phonon relaxation times and κlat. A 95 × 95 × 1 q-mesh is employed to obtain the phonon lifetimes. Long-range electrostatic interactions are included by evaluating the dielectric tensor and Born effective charges using density functional perturbation theory. The fully optimized structure of single-QL Bi2Te3 is shown in Figure 1. The obtained lattice constants are a = b = 4.38 Å, with a QL thickness of 7.60 Å, in agreement with previous calculations.38,41 The present work is the first to fully optimize the structure of single-QL Bi2Te3 using a vdW functional. The modifications with respect to the PBE results (a = b = 4.43 Å, with a QL thickness of 7.56 Å) amount to about 1%. Because a very accurate description of the electronic band structure is required to predict transport coefficients, we compare results obtained by the PBE and HSE06 functionals. Figure 2 shows for the PBE functional a band gap of 310 meV, while the HSE06 functional yields a value of 570 meV.

Figure 1. Side (left) and top (right) views of single-QL Bi2Te3. Red and blue spheres represent Bi and Te atoms, respectively.

Figure 2. Electronic band structures (hexagonal Brillouin zone) and densities of states at the PBE and HSE06 levels.

For temperatures of 300, 500, and 700 K, Figure 3 shows the variations of S, σ, S2σ, and ZT with the electron and hole concentrations. Because σ is obtained with respect to the electronic relaxation time τ, τ in principle can be determined by fitting to the experimental value of σ at a particular carrier concentration. However, due to a lack of experimental data for single-QL Bi2Te3, earlier studies have used τ of bulk Bi2Te3.41,42 In order to enable comparison to these studies, we employ the same values of 2.13 × 10−14 s (300 K), 1.23 × 10−14 s (500 K), and 1.12 × 10−14 s (700 K). Furthermore, κele is obtained from σ using the Wiedemann−Franz relation, κele = LσT (L being the Lorenz number), which is a good approximation for Bi-based bulk and low-dimensional materials,22,26,52,53 while ref 30 has raised doubts about the accuracy of the Boltztrap code in the metallic limit. It is clear from a comparison of the PBE and HSE06 results in Figure 3 that the accurate description of the band structure by the HSE06 functionals has important effects on the predicted transport coefficients. According to Figure 2, the band forming the valence band edge is much flatter than that forming the conduction band edge, that is, holes have larger effective masses than electrons. Hence, p-doping leads to higher S values than n-doping. At high temperature and low carrier concentration, the minority carriers lower S and enhance σ due to the small band gap.54 For the HSE06 functional at 300 K, a 876

DOI: 10.1021/acsenergylett.6b00289 ACS Energy Lett. 2016, 1, 875−879

Letter

ACS Energy Letters

Figure 3. Transport coefficients as functions of the carrier concentration for p-doping (left) and n-doping (right) at the PBE and HSE06 levels.

Figure 4. (a) Harmonic phonon dispersion relation in the hexagonal Brillouin zone and (b) lattice thermal conductivity as a function of the temperature along with data from previous molecular dynamics studies.39,41

Figure 5. (a) Lattice thermal conductivity as a function of the temperature and (b) cumulative lattice thermal conductivity as a function of the mean free path.

877

DOI: 10.1021/acsenergylett.6b00289 ACS Energy Lett. 2016, 1, 875−879

ACS Energy Letters



ZT value of 0.77 is obtained for a hole concentration of 4.0 × 1020 cm−3, while at 700 K, the maximal ZT value of 2.42 appears for a hole concentration of 3.2 × 1020 cm−3. Employing κlat of ref 41 would result in ZT values of 0.66 at 300 K and 1.96 at 700 K, whereas κlat of ref 39 would yield ZT values of 1.06 at 300 K and 3.01 at 700 K. The phonon dispersion of single-QL Bi2Te3 is shown in Figure 4a, agreeing well with previous results.41 Absence of negative frequencies confirms structural stability. The group velocities of the acoustic branches near the Γ-point turn out to be 1192, 1859, and 2920 m/s, while those of the three lowest optical branches are less than 5 m/s. In Figure 4b, the obtained values of κlat are compared to previous results based on molecular dynamics simulations. These results show a large variation, which is probably due to the employed interatomic potentials. In any case, the description of anharmonic force constants by empirical interatomic potentials is less accurate than our ab initio treatment. The fact that three low-frequency optical branches appear in close proximity to the acoustic branches (see Figure 4) makes significant contributions to the heat conduction possible. However, Figure 5a shows that the contributions of the acoustic branches clearly dominate, as a consequence of the small group velocities of the optical branches. In order to understand the effect of the sample size on the thermal transport, we analyze the cumulative lattice thermal conductivity (κlat c ) as a function of the phonon mean free path for temperatures of 300 and 700 K (see Figure 5b). The values are obtained by cumulating the thermal conductivities of all phonons up to a given mean free path. At 300 K, for example, it turns out that half of the contributions to κlat are from phonons with mean free paths below 320 nm, indicating that appropriate nanostructuring can be used to lower κlat. We note that the maximal electron/hole mean free path in the studied temperature range is smaller than 13 nm, which shows that the electron transport is not affected by the nanostructuring. The thermoelectric performance of single-QL Bi2Te3 has been investigated by means of the Boltzmann transport equations for the electrons and phonons. The electronic band structure, phonon dispersion relation, and interatomic force constants have been calculated using density functional theory. Results have been obtained and compared for the PBE and HSE06 functionals. The improved description of the band structure by the HSE06 hybrid functional is found to be of key importance for an accurate prediction of the thermoelectric properties, which demonstrates that electronic exchange effects play an unexpectedly large role in Bi2Te3. The lattice contribution to the thermal conductivity is dominated by the acoustic phonon branches despite the fact that low-frequency optical phonon branches appear. Beyond, the calculated dependence of the lattice thermal conductivity on the phonon mean free path suggests that the thermoelectric performance of single-QL Bi2Te3 can be enhanced efficiently by nanostructuring. Predictions of the lattice thermal conductivity of single-QL Bi2Te3 by means of molecular dynamics simulations, as reported in the literature, show large variations (for example, from 1.63 to 2.97 W K−1 m−1 at 300 K), which are likely due to issues with the employed interatomic potentials. The method used in the present work, on the other hand, does not suffer from this limitation.

Letter

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The research reported in this publication was supported by funding from King Abdullah University of Science and Technology (KAUST).



REFERENCES

(1) Rowe, D. M. CRC Handbook of Thermoelectrics: Macro to Nano. CRC: Boca Raton, FL, 2005. (2) Dresselhaus, M. S.; Chen, G.; Tang, M. Y.; Yang, R.; Lee, H.; Wang, D.; Ren, Z.; Fleurial, J.-P.; Gogna, P. New Directions for LowDimensional Thermoelectric Materials. Adv. Mater. 2007, 19, 1043− 1053. (3) Snyder, G. J.; Toberer, E. S. Complex Thermoelectric Materials. Nat. Mater. 2008, 7, 105−114. (4) Goldsmid, H. J. Introduction to Thermoelectricity; Springer Series in Materials Science; Springer: Berlin, Germany, 2009. (5) Alam, H.; Ramakrishna, S. A Review on the Enhancement of Figure of Merit from Bulk to Nano-Thermoelectric Materials. Nano Energy 2013, 2, 190−212. (6) Yang, J.; Xi, L.; Qiu, W.; Wu, L.; Shi, X.; Chen, L.; Yang, J.; Zhang, W.; Uher, C.; Singh, D. J. On the Tuning of Electrical and Thermal Transport in Thermoelectrics: An Integrated TheoryExperiment Perspective. npj Comput. Mater. 2016, 2, 15015. (7) Hicks, L. D.; Dresselhaus, M. S. Effect of Quantum-Well Structures on the Thermoelectric Figure of Merit. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 12727−12731. (8) Hicks, L. D.; Dresselhaus, M. S. Thermoelectric Figure of Merit of a One-Dimensional Conductor. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 16631−16634. (9) Balandin, A. A. Thermal Properties of Graphene and Nanostructured Carbon Materials. Nat. Mater. 2011, 10, 569−581. (10) Lv, H. Y.; Lu, W. J.; Shao, D. F.; Sun, Y. P. Enhanced Thermoelectric Performance of Phosphorene by Strain-Induced Band Convergence. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 90, 085433. (11) Fei, R.; Faghaninia, A.; Soklaski, R.; Yan, J.-A.; Lo, C.; Yang, L. Enhanced Thermoelectric Efficiency Via Orthogonal Electrical and Thermal Conductances in Phosphorene. Nano Lett. 2014, 14, 6393− 6399. (12) Yang, K.; Cahangirov, S.; Cantarero, A.; Rubio, A.; D’Agosta, R. Thermoelectric Properties of Atomically Thin Silicene and Germanene Nanostructures. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 125403. (13) Peng, B.; Zhang, H.; Shao, H.; Xu, Y.; Zhang, X.; Zhu, H. Low Lattice Thermal Conductivity of Stanene. Sci. Rep. 2016, 6, 20225. (14) Huang, W.; Da, H.; Liang, G. Thermoelectric performance of MX2 (M = Mo,W; X = S, Se) monolayers. J. Appl. Phys. 2013, 113, 104304. (15) Wickramaratne, D.; Zahid, F.; Lake, R. K. Electronic and Thermoelectric Properties of Few-Layer Transition Metal Dichalcogenides. J. Chem. Phys. 2014, 140, 124710. (16) Naguib, M.; Mochalin, V. N.; Barsoum, M. W.; Gogotsi, Y. 25th Anniversary Article: MXenes: A New Family of Two-Dimensional Materials. Adv. Mater. 2014, 26, 992−1005. (17) Khazaei, M.; Arai, M.; Sasaki, T.; Estili, M.; Sakka, Y. TwoDimensional Molybdenum Carbides: Potential Thermoelectric Materials of the MXene Family. Phys. Chem. Chem. Phys. 2014, 16, 7841−7849. (18) Wang, F. Q.; Zhang, S.; Yu, J.; Wang, Q. Thermoelectric Properties of Single-Layered SnSe Sheet. Nanoscale 2015, 7, 15962− 15970. 878

DOI: 10.1021/acsenergylett.6b00289 ACS Energy Lett. 2016, 1, 875−879

Letter

ACS Energy Letters (19) Wright, D. A. Thermoelectric Properties of Bismuth Telluride and Its Alloys. Nature 1958, 181, 834. (20) Goldsmid, H. J. Recent Studies of Bismuth Telluride and Its Alloys. J. Appl. Phys. 1961, 32, 2198−2202. (21) Yim, W. M.; Fitzke, E. V.; Rosi, F. D. Thermoelectric Properties of Bi2Te3-Sb2Te3-Sb2Se3 Pseudo-Ternary Alloys in the Temperature Range 77 to 300 K. J. Mater. Sci. 1966, 1, 52−65. (22) Venkatasubramanian, R.; Siivola, E.; Colpitts, T.; O’Quinn, B. Thin-Film Thermoelectric Devices with High Room-Temperature Figures of Merit. Nature 2001, 413, 597−602. (23) da Silva, L. W.; Kaviany, M.; Uher, C. Thermoelectric Performance of Films in the Bismuth-Tellurium and AntimonyTellurium Systems. J. Appl. Phys. 2005, 97, 114903. (24) Poudel, B.; Hao, Q.; Ma, Y.; Lan, Y. C.; Minnich, A.; Yu, B.; Yan, X.; Wang, D. Z.; Muto, A.; Vashaee, D.; et al. HighThermoelectric Performance of Nanostructured Bismuth Antimony Telluride Bulk Alloys. Science 2008, 320, 634−638. (25) Soni, A.; Yanyuan, Z.; Ligen, Y.; Aik, M. K. K.; Dresselhaus, M. S.; Xiong, Q. Enhanced Thermoelectric Properties of Solution Grown Bi2Te3−xSex Nanoplatelet Composites. Nano Lett. 2012, 12, 1203− 1209. (26) Saleemi, M.; Toprak, M. S.; Li, Shanghua.; Johnsson, M.; Muhammed, M. Synthesis, Processing, and Thermoelectric Properties of Bulk Nanostructured Bismuth Telluride (Bi2Te3). J. Mater. Chem. 2012, 22, 725−730. (27) Goldsmid, H. J. Bismuth Telluride and Its Alloys as Materials for Thermoelectric Generation. Materials 2014, 7, 2577−2592. (28) Yang, L.; Chen, Z.-G.; Hong, M.; Han, G.; Zou, J. Enhanced Thermoelectric Performance of Nanostructured Bi2Te3 through Significant Phonon Scattering. ACS Appl. Mater. Interfaces 2015, 7, 23694−23699. (29) Mishra, S. K.; Satpathy, S.; Jepsen, O. Electronic Structure and Thermoelectric Properties of Bismuth Telluride and Bismuth Selenide. J. Phys.: Condens. Matter 1997, 9, 461−470. (30) Huang, B.-L.; Kaviany, M. Ab initio and Molecular Dynamics Predictions for Electron and Phonon Transport in Bismuth Telluride. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 125209. (31) Park, M. S.; Song, J.-H.; Medvedeva, J. E.; Kim, M.; Kim, I. G.; Freeman, A. J. Electronic Structure and Volume Effect on Thermoelectric Transport in p-Type Bi and Sb Tellurides. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 155211. (32) Luo, X.; Sullivan, M. B.; Quek, S. Y. First-Principles Investigations of the Atomic, Electronic, and Thermoelectric Properties of Equilibrium and Strained Bi2Se3 and Bi2Te3 Including Van Der Waals Interactions. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86, 184111. (33) Ibarra-Hernández, W.; Verstraete, M. J.; Raty, J.-Y. Effect of Hydrostatic Pressure on the Thermoelectric Properties of Bi2Te3. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 90, 245204. (34) Cheng, L.; Liu, H. J.; Zhang, J.; Wei, J.; Liang, J. H.; Shi, J.; Tang, X. F. Effects of Van Der Waals Interactions and Quasiparticle Corrections on the Electronic and Transport Properties of Bi2Te3. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 90, 085118. (35) Feutelais, Y.; Legendre, B.; Rodier, N.; Agafonov, V. A Study of the Phases in the Bismuth-Tellurium System. Mater. Res. Bull. 1993, 28, 591−596. (36) Teweldebrhan, D.; Goyal, V.; Balandin, A. A. Exfoliation and Characterization of Bismuth Telluride Atomic Quintuples and QuasiTwo-Dimensional Crystals. Nano Lett. 2010, 10, 1209−1218. (37) Goyal, V.; Teweldebrhan, D.; Balandin, A. A. MechanicallyExfoliated Stacks of Thin Films of Bi2Te3 Topological Insulators with Enhanced Thermoelectric Performance. Appl. Phys. Lett. 2010, 97, 133117. (38) Li, X.; Ren, H.; Luo, Y. Electronic Structure of Bismuth Telluride Quasi-Two-Dimensional Crystal: A First Principles Study. Appl. Phys. Lett. 2011, 98, 083113. (39) Qiu, B.; Ruan, X. Thermal Conductivity Prediction and Analysis of Few-Quintuple Bi2Te3 Thin Films: A Molecular Dynamics Study. Appl. Phys. Lett. 2010, 97, 183107.

(40) Zahid, F.; Lake, R. Thermoelectric Properties of Bi2Te3 Atomic Quintuple Thin Films. Appl. Phys. Lett. 2010, 97, 212102. (41) Zhang, J.; Liu, H. J.; Cheng, L.; Wei, J.; Shi, J.; Tang, X. F.; Uher, C. Enhanced Thermoelectric Performance of a Quintuple Layer of Bi2Te3. J. Appl. Phys. 2014, 116, 023706. (42) Zhou, G.; Wang, D. Few-Quintuple Bi2Te3 Nanofilms as Potential Thermoelectric Materials. Sci. Rep. 2015, 5, 8099. (43) Li, W.; Carrete, J.; Katcho, N. A.; Mingo, N. ShengBTE: A Solver of the Boltzmann Transport Equation for Phonons. Comput. Phys. Commun. 2014, 185, 1747−1758. (44) Ibarra-Hernádez, W. Ab-Initio Study of Thermoelectricity of Layered Tellurium Compounds. Ph.D. Thesis, Universitè de Liége, 2015. (45) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (46) Heyd, J.; Scuseria, G. E. Efficient Hybrid Density Functional Calculations in Solids: Assessment of the Heyd-Scuseria-Ernzerhof Screened Coulomb Hybrid Functional. J. Chem. Phys. 2004, 121, 1187−1192. (47) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758−1775. (48) Klimeš, J.; Bowler, D. R.; Michaelides, A. Van Der Waals Density Functionals Applied to Solids. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 83, 195131. (49) Madsen, G. K. H.; Singh, D. J. BoltzTraP. A Code for Calculating Band-Structure Dependent Quantities. Comput. Phys. Commun. 2006, 175, 67−71. (50) Togo, A.; Oba, F.; Tanaka, I. First-Principles Calculations of the Ferroelastic Transition Between Rutile-Type and CaCl2-type SiO2 at High Pressures. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 78, 134106. (51) Gonze, X.; Lee, C. Dynamical Matrices, Born Effective Charges, Dielectric Permittivity Tensors, and Interatomic Force Constants from Density-Functional Perturbation Theory. Phys. Rev. B: Condens. Matter Mater. Phys. 1997, 55, 10355−10368. (52) Cui, J.; Xiu, W.; Xue, H. High Thermoelectric Properties of ptype Pseudobinary (Cu4Te3)x(Bi0.5Sb1.5Te3)1−x Alloys Prepared by Spark Plasma Sintering. J. Appl. Phys. 2007, 101, 123713. (53) Kim, K. T.; Choi, S. Y.; Shin, E. H.; Moon, K. S.; Koo, H. Y.; Lee, G.-G.; Ha, G. H. The Influence of CNTs on the Thermoelectric Properties of a CNT/Bi2Te3 Composite. Carbon 2013, 52, 541−549. (54) May, A. F.; McGuire, M. A.; Singh, D. J.; Ma, J.; Delaire, O.; Huq, A.; Cai, W.; Wang, H. Thermoelectric Transport Properties of CaMg2Bi2, EuMg2Bi2, and YbMg2Bi2. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 035202.

879

DOI: 10.1021/acsenergylett.6b00289 ACS Energy Lett. 2016, 1, 875−879