Silicon As an Unexpected n-Type Dopant in BiCuSeO Thermoelectrics

Aug 3, 2017 - Page 1. Silicon As an Unexpected n‑Type Dopant in BiCuSeO Thermoelectrics. Jiahong Shen. †,‡ and Yue Chen*,†. †. Department of...
0 downloads 12 Views 2MB Size
Subscriber access provided by BOSTON UNIV

Letter

Silicon as an unexpected n-type dopant in BiCuSeO thermoelectrics Jiahong Shen, and Yue Chen ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.7b06872 • Publication Date (Web): 03 Aug 2017 Downloaded from http://pubs.acs.org on August 6, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Applied Materials & Interfaces is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Silicon as an unexpected n-type dopant in BiCuSeO thermoelectrics Jiahong Shen1,2 and Yue Chen1,* 1. Department of Mechanical Engineering, The University of Hong Kong, Pokfulam Road, Hong Kong SAR, China 2. Department of Materials Science, Fudan University, 220 Handan Road, Shanghai 200433, China

Supporting Information Placeholder ABSTRACT: As a promising thermoelectric material,

BiCuSeO is of great interest for energy conversion. A higher figure of merit in n-type BiCuSeO than that in the p-type was predicted from theory, suggesting a need of in-depth investigations on the doping effects. In this work, the influences of group IV elements (Si, Ge, Sn, and Pb) on the electronic structures of BiCuSeO are studied from first principles. Despite the similar electronegativities of the group IV elements, Si is found to be an n-type dopant, being distinctly different from Ge, Sn, and Pb, which exhibit typical p-type behaviors. Detailed analysis on the doping effects is performed based on a recently developed band unfolding technique. Furthermore, Si-doped BiCuSeO is shown to have a higher power factor than p-type BiCuSeO from the Boltzmann transport theory.

sistance. Nonetheless, the ZT values of many widely studied TE oxides, including SrTiO36, In2O37, and Ca3Co4O98 remain relatively low.

KEYWORDS: thermoelectrics, defect formation energy, electronic structure, charge density, Boltzmann transport theory

Thermoelectric (TE) materials play an important role in energy conservation. Their applications including refrigeration and thermal power generation1 provide a prospective solution to environmental issues related to energy. The TE performance is determined by the dimensionless figure of merit (ZT), defined as ZT = (S2σT)/κ, where S, σ, κ, and T are Seebeck coefficient, electrical conductivity, thermal conductivity, and absolute temperature, respectively2. Several classes of bulk materials have been demonstrated to exhibit high ZT values, including PbTe3, BiSbTe4, and Bi2Te35; however, the applications of these materials are restrained due to the limitations in chemical and thermal stabilities at high temperatures. Thermoelectric oxides display high thermal and chemical stabilities, and good oxidation re-

Figure 1. (a) Primitive unit cell of BiCuSeO with a layered crystal structure. (b) Theoretical power factor of pristine BiCuSeO normalized by the relaxation time as a function of carrier concentration. (c) Formation energies of different doped BiCuSeO systems calculated from density functional theory.

Recently, the layered quaternary oxyselenide BiCuSeO with its naturally layered structure (see Fig. 1a), was found to be a good TE material in high temperature range due to its intrinsically low lattice thermal conductivity (~0.45 Wm-1K-1 at 923 K9). Different approaches have been conducted to further improve the energy conversion efficiency of BiCuSeO, including the introduction of Cu vacancies10, carrier concentration modification11 and band gap tuning12. A high ZT value of 1.4 was reported by Sui et al13 in Ba-doped BiCuSeO at 923 K. Most studies availbale in the literatures focus on p-type

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

BiCuSeO; e.g., the doping of alkali metals (Na11 and K14), alkaline earth metals (Mg15, Ca9, Sr16, and Ba13), Pb17, La12, and Ni18 on the Bi sublattice site, the doping of Zn19, Cd19, and Ag20 on the Cu sublattice site, and the doping of Te21 on the Se sublattice site. N-type BiCuSeO is relatively less studied, but Yang et al22 predicted that its ZT value can be higher than p-type BiCuSeO by 32% because of the increase of power factor. BiCuSeO adopts a crystal structure with the tetragonal P4/nmm space group. The fluorite-like (Bi2O2)2+ layers and the (Cu2Se2)2- layers are alternatively stacked along the c axis of the unit cell23. As the (001) plane of BiCuSeO has the lowest surface energy, c-axis-oriented growth of BiCuSeO has been achieved in experiments17. We have calculated the electrical transport properties along the caxis from Boltzmann transport theory. It is seen from Fig. 1b that the power factor of n-type BiCuSeO is much higher than that of p-type BiCuSeO. Here, we investigate the doping effects of group IV elements (Si, Ge, Sn, and Pb) on the electrical transport properties and the electronic structures of BiCuSeO. It was found in experiments24 that Pb exhibits a +2 charge state and Bi shows a +3 charge state in BiCuSeO; thus, the substitution of Pb for the Bi atom leads to a p-type doping behavior. Given that the electronegativities of group IV elements are similar (1.916, 1.994, 1.824, and 1.854 in Pauling units25 for Si, Ge, Sn, and Pb, respectively), it is expected that they all have p-type doping behaviors when they occupy the same sublattice site in BiCuSeO. Interestingly, Si is found to be a promising ntype dopant whereas other elements in this group lead to a p-type semiconductor as expected. Effective band structures and charge density distributions are calculated to provide insight into the different doping effects. Electrical transport properties of the doped systems are also calculated and compared to available experimental measurements. The lattice constants of BiCuSeO obtained from our calculations are a = 3.955 Å and c = 9.083 Å, and the experimental values are a = 3.928 Å and c = 8.929 Å26. It is well-known that the generalized gradient approximation (GGA) usually overestimates the lattice constants of solids27; thus, the small differences between our theoretical lattice constants and the experimental values are reasonable. To investigate the doping effects, one atom in a 3×3×1 supercell is substituted, corresponding to a doping concentration of 5.56%. The formation energy is calculated according to the following equation28 to determine the preferable substitutional lattice site:               ∆  ∆

where R represents the doping elements (Si, Ge, Sn, or Pb) and X represents the different sublattice sites (Bi, Cu, Se, or O).  and   are the total energies of the R-doped supercell and the pristine supercell, respectively.  and  are the energies of one R atom and one X

Page 2 of 12

atom in their respective bulk phases.  represents the charge states.  refers to the energy at the top of the valence band of the pristine system and ∆ is the Fermi energy relative to  . ∆ is the correction term to aligh the reference potentials of pristine and doped BiCuSeO. The results corresponding to neutral charge states (q=0) are shown in Fig. 1c; it is seen that all doping elements tend to occupy the Bi site because of the lowest formation energies. With a similar atomic radius to Bi, Pb has a relatively low formation energy comparing to Sn and Ge. Interestingly, Si-doped system has the lowest formation energy among the four elements. The valence states of Bi, Cu, Se and O atoms in BiCuSeO are +3, +1, -2 and -2, respectively. The common valence state of both Si and Ge is +4, while the common valence states of both Sn and Pb are +4 and +2. Thus, we have considered the charge states of +1, 0 and -1 when the dopants occupy the Bi site and the charge states of +3, +1 and 0 when the dopants occupy the Cu site. The Se and O sites are not considered for charged dopants because these sites are unlikely to be occupied due to their negative charges, which is also evidenced in the neutral defect formation energy calculations. It is seen from Fig. 2 that the defect formation energies corresponding to the Bi site are always lower than those corresponding to the Cu site for all doping systems; therefore, the Bi site is more likely to be occupied by the dopants regardless of the charge states and chemical potentials. In the following discussion, we focus on the systems in which the doping atoms occupy their most preferable lattice site.

Figure 2. Formation energies of (a) Si-doped, (b) Ge-doped, (c) Sn-doped and (d) Pb-doped BiCuSeO corresponding to different charge states of dopants. The solid lines represent the occupations of the Bi site, and the dashed lines represent the occupations of the Cu site.

Electronic structures are calculated to analyze the doping effects of group IV elements. As can be seen in Fig. 3a, Ge-, Sn-, and Pb-doped BiCuSeO have a similar density of states (DOS) with Fermi level moving slightly into the valence band, indicating that these elements are p-type dopants. This agrees with the expectation that Ge and Sn have similar doping effects with Pb because they

ACS Paragon Plus Environment

Page 3 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

have comparable electronegativities and they all tend to occupy the Bi sublattice site. We choose Pb-doped BiCuSeO as the representative p-type system since it has already been synthesized in experiments24. Most interestingly, the DOS of Si-doped BiCuSeO is distinctly different from other systems (see Fig. 3b); Fermi level moves into the conduction band, leading to an n-type semiconductor. The n-type behavior of Si-doped BiCuSeO indicates that Si may exhibit a +4 charge state, in comparison with the Pb-doped system where Pb exhibits a +2 charge state24. Because of the similar electronegativities of the group IV elements, the distinctly different doping effects of Si are unexpected and require further investigations. To unveil the detailed modifications to the band dispersions, we have calculated the unfolded effective band structures of Pb and Si-doped BiCuSeO, as shown in Fig. 3c and 3d. The conduction band minimum (CBM) of both systems is between the Γ and Ζ points, while the valence band maximum (VBM) is in the Μ-Γ and Z-R directions. In Pb-doped system, the energy of VBM exceeds Fermi level, corresponding to the p-type behavior observed in DOS. The effective band structures of Sn- and Ge-doped BiCuSeO are similar with that of the Pb-doped system (see Supporting Information). From the effective band structure of the Si-doped system, it is seen that Fermi level is higher than the CBM, which is in consistent with the n-type behavior observed in DOS. In addition, the weights of electronic bands near Fermi level are significantly modified due to the introduction of Si impurity in BiCuSeO, indicating changes in the electrical transport properties.

Bader charge analysis29 is performed to study the charge states of different atoms. As can be seen from Table 1, Ge, Sn, and Pb donate fewer electrons than Bi, rationalizing the p-type behavior of the doped systems where holes dominate the electrical transport. On the other hand, Si donates more electrons than Bi, indicating that electrons are the majority carriers and the system becomes an n-type semiconductor. The Bader charges of Bi, Cu, Se, and O atoms are less affected by the doping of Ge, Sn, or Pb comparing to Si; particularly, Si induces a large change in the Bader charge of O, which suggests that Si may have a strong interaction with the O atom.

Figure 3. Density of states of (a) Ge-, Sn-, and Pb-doped BiCuSeO and (b) Si-doped BiCuSeO. Effective band structures of (c) Pb-doped and (d) Si-doped BiCuSeO. Fermi level is located at 0 eV

Figure 4. Charge density distributions in Pb-doped BiCuSeO; (a) top view, (b) side view, and (c) the (001) plane crossing the oxygen atomic layer. Charge density distributions in Si-doped BiCuSeO; (d) top view, (e) side view, and (f) the (001) plane crossing the oxygen atomic layer. The isosurface level is 0.12 e/Å3.

Table 1. Bader charge of different dopants and of the Bi, Cu, Se, and O atoms nearest to the doping atoms.

pristine Si-doped Ge-doped Sn-doped Pb-doped

Bi 1.59 1.54 1.60 1.60 1.60

Cu 0.30 0.31 0.32 0.31 0.31

Se -0.79 -0.85 -0.74 -0.76 -0.78

O -1.10 -1.31 -1.11 -1.12 -1.10

R(dopants) 2.80 1.19 1.27 1.15

The total charge density distributions of Pb- and Sidoped BiCuSeO are shown in Fig. 4. These two systems have similar charge density distributions around the Bi, Cu, and Se atoms, in agreement with the Bader charge analysis. Nonetheless, the charge densities around the O atoms are significantly different in Pb- and Si-doped BiCuSeO; a stronger chemical bond is found to exist between the Si and O atoms while a weaker interaction is observed between the Pb and O atoms. The Si and O interaction is so strong that the Si atom is attracted to the oxygen layer, rationalizing the lowest formation energy of Si-doped BiCuSeO. It can be clearly seen from Fig. 4f that these strong Si-O bonds result in a transfer of electrons from Si to O, corresponding to the increase of Bader charge of the O atom.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

For 5.56% Pb-doped BiCuSeO, we obtain a carrier concentration of 7.82  10 cm-3, which agrees reasonably with the experimental result24 of 8.39  10 cm-3 for a 5% doping concentration. Based on the relation "  #$ %& '%&/) , we know that relaxation time " decreases as temperature T increases and " is proportional to '%&/) , in which ' represents the carrier concentration30. The experimental electrical conductivity24 of 5% Pb doping concentration is 1.11  10+ Ωm%& at 700 K. Comparing this value to the electrical transport properties calculated from Boltzmann transport theory, we obtain "  7.392  10%. $ %& '%&/) ; inserting this value to the above equation, we determine the value of C, which is 7.392  10%. Kscm%& . Based on the expression of relaxation time, it is seen from Fig. 5a that our theoretical electrical conductivity and Seebeck coefficient are in good agreements with experiments.

Page 4 of 12

trations. The predicted increase of power factor is promising for achieving high performance n-type BiCuSeO. In summary, the doping effects of group IV elements (Si, Ge, Sn, and Pb) on the electronic structures and electrical transport properties of BiCuSeO are investigated from first-principles calculations. The site preferences of dopants are determined using the defect formation energies and it is found that all of these dopants tend to occupy the Bi sublattice site regardless of the charge states and chemical potentials. The doping effects of Si are found to be distinctly different from other elements from group IV; Si is found to be an n-type dopant in BiCuSeO while other elements (Ge, Sn, and Pb) are p-type dopants. A recently developed band unfolding technique is utilized to obtain the effective band structures of the doped systems; significant band modifications are observed in Si-doped BiCuSeO near Fermi level. From Bader charge analysis and charge density distributions, we rationalize the different doping effects of Si. Despite the similar electronegativities of the group IV elements, they exhibit different charge states in BiCuSeO. Electrical transport properties are calculated using the Boltzmann transport theory and Si-doped BiCuSeO is found to be a promising n-type thermoelectric material with an enhanced power factor. ASSOCIATED CONTENT

Figure 5. (a) Seebeck coefficient and electrical conductivity of Pb-doped polycrystalline BiCuSeO at 700 K together with the experimental results24 (b) Absolute values of the Seebeck coefficients as functions of carrier concentration at 700 K (dashed lines) and 900 K (solid lines). (c) Electrical conductivities normalized by the relaxation time at different carrier concentrations. (d) Power factors normalized by the relaxation time.

The electrical transport properties along the c-axis are shown in Fig. 5b-5d. It is seen that the absolute value of Seebeck coefficient of Si-doped n-type BiCuSeO is lower than those of the pristine and Pb-doped p-type BiCuSeO. On the other hand, the electrical conductivity normalized by relaxation time of Si-doped BiCuSeO is much higher than those of the Pb-doped and pristine BiCuSeO (see Fig. 5c). As can be seen from Fig. 3c and 3d, the effective mass of the top of the valence band is larger than that of the bottom of the conduction band; therefore, the carrier mobility in Si-doped BiCuSeO is higher. Because conductivity is proportional to carrier mobility, the Si-doped BiCuSeO has a higher conductivity. The power factors of these systems are compared in Fig. 5d; it is seen that the Si-doped BiCuSeO has a larger power factor over a broad range of carrier concen-

Supporting Information. Computational details, formation energies, DOS and band structures from full relaxations (Fig. S1); effective band structures and charge density distributions of pristine, Sn-, and Ge-doped BiCuSeO (Fig. S2); effective band structures of Si-doped BiCuSeO at different doping concentrations (Fig. S3). The Supporting Information is available free of charge on the ACS Publications website.

AUTHOR INFORMATION Corresponding Author

[email protected] Funding Sources

Early Career Scheme of RGC

ACKNOWLEDGMENT We are grateful for the financial support from RGC under project numbers 27202516 and 17200017, and the research computing facilities offered by ITS, HKU.

REFERENCES (1) Bell, L. E. Cooling, Heating, Generating Power, and Recovering Waste Heat with Thermoelectric Systems. Science 2008, 321, 1457-1461. (2) Snyder, G. J.; Toberer, E. S. Complex Thermoelectric Materials. Nat. Mater. 2008, 7, 105-114. (3) Zhang, Q.; Cao, F.; Liu, W.; Lukas, K.; Yu, B.; Chen, S.; Opeil, C.; Broido, D.; Chen, G.; Ren, Z. Heavy Doping and Band Engineering by

ACS Paragon Plus Environment

Page 5 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Potassium to Improve the Thermoelectric Figure of Merit in p-Type PbTe, PbSe, and PbTe1-ySey. J. Am. Chem. Soc. 2012, 134, 10031-10038. (4) Poudel, B.; Hao, Q.; Ma, Y.; Lan, Y.; Minnich, A.; Yu, B.; Yan, X.; Wang, D.; Muto, A.; Vashaee, D.; Chen, X.; Liu, J.; Dresselhaus, M. S.; Chen, G.; Ren, Z. High-thermoelectric Performance of Nanostructured Bismuth Antimony Telluride Bulk Alloys. Science 2008, 320, 634-638. (5) Saleemi, M.; Toprak, M. S.; Li, S.; Johnsson, M.; Muhammed, M. Synthesis, Processing, and Thermoelectric Properties of Bulk Nanostructured Bismuth Telluride (Bi2Te3). J. Mater. Chem. 2012, 22, 725-730. (6) Ohta, H.; Kim, S.; Mune, Y.; Mizoguchi, T.; Nomura, K.; Ohta, S.; Nomura, T.; Nakanishi, Y.; Ikuhara, Y.; Hirano, M.; Hosono, H.; Koumoto, K. Giant Thermoelectric Seebeck Coefficient of Twodimensional Electron Gas in SrTiO3. Nat. Mater. 2007, 6, 129-134. (7) Ohta, H.; Seo, W. S.; Koumoto, K. Thermoelectric Properties of Homologous Compounds in the ZnO-In2O3 System. J. Am. Ceram. Soc. 1996, 79, 2193-2196. (8) Masuda, Y.; Nagahama, D.; Itahara, H.; Tani, T.; Seo, W. S.; Koumoto, K. Thermoelectric Performance of Bi- and Na-substituted Ca3Co4O9 Improved Through Ceramic Texturing. J. Mater. Chem. 2003, 13, 1094-1099. (9) Pei, Y.-L.; He, J.; Li, J.-F.; Li, F.; Liu, Q.; Pan, W.; Barreteau, C.; Berardan, D.; Dragoe, N.; Zhao, L.-D. High Thermoelectric Performance of Oxyselenides: Intrinsically Low Thermal Conductivity of Ca-doped BiCuSeO. Npg Asia Mater. 2013, 5, e47. (10) Liu, Y.; Zhao, L.-D.; Liu, Y.; Lan, J.; Xu, W.; Li, F.; Zhang, B.-P.; Berardan, D.; Dragoe, N.; Lin, Y.-H.; Nan, C.-W.; Li, J.-F.; Zhu, H. Remarkable Enhancement in Thermoelectric Performance of BiCuSeO by Cu Deficiencies. J. Am. Chem. Soc. 2012, 134, 3312-3312. (11) Li, J.; Sui, J.; Pei, Y.; Meng, X.; Berardan, D.; Dragoe, N.; Cai, W.; Zhao, L.-D. The Roles of Na Doping in BiCuSeO Oxyselenides As a Thermoelectric Material. J. Mater. Chem. A 2014, 2, 4903-4906. (12) Liu, Y.; Ding, J.; Xu, B.; Lan, J.; Zheng, Y.; Zhan, B.; Zhang, B.; Lin, Y.; Nan, C. Enhanced Thermoelectric Performance of La-doped BiCuSeO by Tuning Band Structure. Appl. Phys. Lett. 2015, 106, 233903. (13) Sui, J.; Li, J.; He, J.; Pei, Y.-L.; Berardan, D.; Wu, H.; Dragoe, N.; Cai, W.; Zhao, L.-D. Texturation Boosts the Thermoelectric Performance of BiCuSeO Oxyselenides. Energy Environ. Sci. 2013, 6, 2916-2920. (14) Lee, D. S.; An, T.-H.; Jeong, M.; Choi, H.-S.; Lim, Y. S.; Seo, W.-S.; Park, C.-H.; Park, C.; Park, H.-H. Density of State Effective Mass and Related Charge Transport Properties in K-doped BiCuOSe. Appl. Phys. Lett. 2013, 103, 232110. (15) Li, J.; Sui, J.; Barreteau, C.; Berardan, D.; Dragoe, N.; Cai, W.; Pei, Y.; Zhao, L.-D. Thermoelectric Properties of Mg Doped p-Type BiCuSeO Oxyselenides. J. Alloys Compd. 2013, 551, 649-653. (16) Barreteau, C.; Berardan, D.; Amzallag, E.; Zhao, L.; Dragoe, N. Structural and Electronic Transport Properties in Sr-Doped BiCuSeO. Chem. Mater. 2012, 24, 3168-3178.

(17) Wu, X.; Wang, J.-L.; Zhang, H.; Wang, S.; Zhai, S.; Li, Y.; Elhadj, D.; Fu, G. Epitaxial Growth and Thermoelectric Properties of c-Axis Oriented Bi1-xPbxCuSeO Single Crystalline Thin Films. Crystengcomm 2015, 17, 8697-8702. (18) Liu, Y.; Liu, J.; Zhang, B.; Lin, Y. Thermoelectric Properties of Ni Doped p-Type BiCuSeO Oxyselenides. In High-Performance Ceramics Viii, Pan, W.; Gong, J., Eds. 2014; Vol. 602-603, pp 906-909. (19) Luu, S. D. N.; Vaqueiro, P. Thermoelectric Properties of BiOCu1xMxSe (M = Cd and Zn). Semicond. Sci. Technol. 2014, 29, 064002. (20) Tan, S. G.; Lei, H.; Shao, D. F.; Lv, H. Y.; Lu, W. J.; Huang, Y. N.; Liu, Y.; Yuan, B.; Zu, L.; Kan, X. C.; Song, W. H.; Sun, Y. P. Enhanced Low Temperature Thermoelectric Performance of Ag-doped BiCuSeO. Appl. Phys. Lett. 2014, 105, 082109. (21) Barreteau, C.; Berardan, D.; Zhao, L.; Dragoe, N. Influence of Te Substitution on the Structural and Electronic Properties of Thermoelectric BiCuSeO. J. Mater. Chem. A 2013, 1, 2921-2926. (22) Yang, J.; Yang, G.; Zhang, G.; Wang, Y. X. Low Effective Mass Leading to An Improved ZT Value by 32% for n-Type BiCuSeO: A Firstprinciples Study. J. Mater. Chem. A 2014, 2, 13923-13931. (23) Hiramatsu, H.; Kamioka, H.; Ueda, K.; Hirano, M.; Hosono, H. Electrical and Photonic Functions Originating from Low-dimensional Structures in Wide-gap Semiconductors LnCuOCh (Ln = lanthanide, Ch = chalcogen): A Review. J. Ceram. Soc. Jpn. 2005, 113, 10-16. (24) Pan, L.; Berardan, D.; Zhao, L.; Barreteau, C.; Dragoe, N. Influence of Pb Doping on the Electrical Transport Properties of BiCuSeO. Appl. Phys. Lett. 2013, 102, 023902. (25) Mann, J. B.; Meek, T. L.; Allen, L. C. Configuration Energies of the Main Group Elements. J. Am. Chem. Soc. 2000, 122, 2780-2783. (26) Hiramatsu, H.; Yanagi, H.; Kamiya, T.; Ueda, K.; Hirano, M.; Hosono, H. Crystal Structures, Optoelectronic Properties, and Electronic Structures of Layered Oxychalcogenides MCuOCh (M = Bi, La; Ch = S, Se, Te): Effects of Electronic Configurations of M3+ Ions. Chem. Mater. 2008, 20, 326-334. (27) Haas, P.; Tran, F.; Blaha, P. Calculation of the Lattice Constant of Solids with Semilocal Functionals. Phys. Rev. B 2009, 79, 085104. (28) Li, X. P.; Xia, C. X.; Song, X. H.; Du, J.; Xiong, W. Q. N- and P-type Dopants in the InSe Monolayer Via Substitutional Doping. J. Mater. Sci. 2017, 52, 7207-7214. (29) Henkelman, G.; Arnaldsson, A.; Jonsson, H. A Fast and Robust Algorithm for Bader Decomposition of Charge Density. Comput. Mater. Sci. 2006, 36, 354-360. (30) Ong, K. P.; Singh, D. J.; Wu, P. Analysis of the Thermoelectric Properties of N-type ZnO. Phys. Rev. B 2011, 83, 5, 115110.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 12

TOC Graphic

6

ACS Paragon Plus Environment

Page 7 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(a) Primitive unit cell of BiCuSeO with a layered crystal structure. (b) Theoretical power factor of pristine BiCuSeO normalized by the relaxation time as a function of carrier concentration. (c) Formation energies of different doped BiCuSeO systems calculated from density functional theory. 79x78mm (600 x 600 DPI)

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Formation energies of (a) Si-doped, (b) Ge-doped, (c) Sn-doped and (d) Pb-doped BiCuSeO corresponding to different charge states of dopants. The solid lines represent the occupations of the Bi site, and the dashed lines represent the occupations of the Cu site. 469x303mm (72 x 72 DPI)

ACS Paragon Plus Environment

Page 8 of 12

Page 9 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Density of states of (a) Ge-, Sn-, and Pb-doped BiCuSeO and (b) Si-doped BiCuSeO. Effective band structures of (c) Pb-doped and (d) Si-doped BiCuSeO. Fermi level is located at 0 eV. 74x76mm (600 x 600 DPI)

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Charge density distributions in Pb-doped BiCuSeO; (a) top view, (b) side view, and (c) the (001) plane crossing the oxygen atomic layer. Charge density distributions in Si-doped BiCuSeO; (d) top view, (e) side view, and (f) the (001) plane crossing the oxygen atomic layer. The isosurface level is 0.12 e/Å^3. 59x40mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 10 of 12

Page 11 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(a) Seebeck coefficient and electrical conductivity of Pb-doped polycrystalline BiCuSeO at 700 K together with the experimental results21 (b) Absolute values of the Seebeck coefficients as functions of carrier concentration at 700 K (dashed lines) and 900 K (solid lines). (c) Electrical conductivities normalized by the relaxation time at different carrier concentrations. (d) Power factors normalized by the relaxation time. 61x47mm (600 x 600 DPI)

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC graphic 35x17mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 12 of 12