High-Throughput Screening for Advanced Thermoelectric Materials

6 days ago - ACS Editors' Choice: Unprecedented Ultralow Detection Limit of Amines — and More! This week: Unprecedented ultralow detection limit of ...
0 downloads 0 Views 2MB Size
Forum Article www.acsami.org

Cite This: ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

High-Throughput Screening for Advanced Thermoelectric Materials: Diamond-Like ABX2 Compounds Ruoxi Li,† Xin Li,† Lili Xi,*,† Jiong Yang,*,† David J. Singh,‡ and Wenqing Zhang*,#,⊥ †

Materials Genome Institute, Shanghai University, 99 Shangda Road, Shanghai 200444, China Department of Physics and Astronomy, University of Missouri, Columbia, Missouri 65211, United States # Department of Physics and Shenzhen Institute for Quantum Science & Technology, Southern University of Science and Technology, No. 1088 Xueyuan Road, Shenzhen, Guangdong 518055, China ⊥ Center for Quantum Computing, Peng Cheng Laboratory, Shenzhen 518000, China Downloaded via UNIV OF LOUISIANA AT LAFAYETTE on April 26, 2019 at 20:32:02 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.



S Supporting Information *

ABSTRACT: High-throughput (HTP) calculations are a highly promising direction for the discovery of novel functional materials. Here we use an HTP framework to investigate the electronic structures and p-type thermoelectric properties of the ABX2 compounds with diamond-like structures. We show application of HTP both to identify compounds and also to identify underlying trends. A total of 65 entries out of 84 908 in the Materials Informatics Platform are selected for this study. The electronic structures and chemical-bonding analyses reveal that there exists a general conductive network consisting of the anion X sublattice, which dominates the electrical transport properties of the compounds. Electrical and thermal transport properties of the 41 pnictide and chalcogenide compounds with sufficient band gaps are studied. Pnictide compounds have relatively smaller Seebeck coefficients than the chalcogenide compounds. This is due to the smaller effective masses around the valence band maxima. The electrical conductivities and power factors, however, are better in pnictide compounds. This is because pnictide compounds have high electronic group velocities and electronic relaxation times. Combined with the predictions of lattice thermal conductivities based on the Slack model, 12 novel p-type and n-type ABX2 materials with high ZT values are predicted. KEYWORDS: high-throughput, conductive network, thermoelectric, diamond-like structures, chalcopyrite clathrates, half-Heuslers, and diamond-like compounds,9−22 among others, emerged as high-performance thermoelectrics. Some of these recently discovered compounds show exceptionally high TE efficiency, leading to the promise of applications in energy technologies23 as well as suggesting that there are other novel high-performance TE materials waiting to be discovered. In recent years, some progress has been made in the study of TE materials by high-throughput (HTP) screening. Using AFLOWLIB,24 Carrete et al. obtained the relationship between power factor and effective mass, through the calculations of more than 2500 compounds.25 On the basis of the Materials Project (MP),26 Zhu et al. screened more than 9000 materials and predicted TmAgTe2 and other XYZ2 compounds with high TE properties,27 but the related experimental ZT values were low. Subsequently, Aydemir et al. studied YCuTe2 through isoelectronic element substitution, and its performance was improved.28 Chen et al. recently studied the classification of TE

1. INTRODUCTION Thermoelectric (TE) materials, as a new type of energy conversion material, are receiving widespread attention because of the direct conversion between electricity and heat.1−3 The conversion efficiency of TE devices closely relates to the dimensionless TE figure of merit, ZT = S2σT/κ. Here S, σ, κ, and T are the Seebeck coefficient, electrical conductivity, thermal conductivity, and the absolute temperature, respectively. S2σ is the power factor (PF). The total thermal conductivity κ has both lattice and electronic components, κL and κe, respectively. The parameters that determine the ZT value are interrelated. For example, the increase of S causes the decrease of σ. These countercorrelations among σ, S, and κ greatly complicate the optimizing of ZT values.4−6 The electrical transport properties in TE materials are determined by the electronic states within a small energy range, of a few kBT, around the band edge (i.e., valence band maximum (VBM)) or conduction band minimum (CBM). The properties depend strongly on doping level, which must be optimized. In addition to the traditional TE materials such as Bi2Te3, SiGe, and PbTe,1,7,8 compounds with complex structures and multiple components, such as skutterudites, © XXXX American Chemical Society

Special Issue: Materials Discovery and Design Received: January 19, 2019 Accepted: April 15, 2019

A

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces materials through HTP computation and machine learning.29 Recently, Xi et al.20 predicted novel high-performance TE materials using an HTP workflow based on the Materials Informatics Platform (MIP)30 applied to a class of compounds. This was followed by experimental verification in chalcogenides with diamond-like structure. A number of the recently discovered TE materials contain “conductive networks”, which are composed frameworks involving a portion of the atoms that dominate the electrical transport properties. Other atoms, because their electronic states are not close to the band edges, only serve as carrier reservoirs. This is related to the Zintl concept but is more general. The existence of conductive networks in a compound facilitates dopability, while retaining mobility, and therefore enables the optimization of carrier concentrations using chemical strategies. Diamond-like materials are a new series of promising TE materials. They are characterized by high Seebeck coefficients and low thermal conductivities.12−21 Examples include AgGaTe2, with a reported ZT of 0.77 at 850 K,17 CuGaTe2, with a κl value as low as 0.89 W/m·K and a ZT of 1.4,18 and CuInTe2, with a ZT of 1.18.19 First-principles studies show that the VBM is characterized by an anion-dominated conductive network at the VBM in chalcogenide diamond-like compounds.20−22 S-, Se-, and Te-based chalcogenide diamond-like compounds have been widely reported. Compounds based on other anions, such as P-, As-, and Sb-based compounds, have been less studied as TE materials,31−33 although many of these have been studied as photovoltaic materials, with findings of reasonable semiconducting properties, such as charge carrier transport.34−36 Here we use the HTP calculations, recently done for chalcogenides,20 to study the TE properties of ternary ABX2 (A = IA, IIA, IIIA, VA, IB, and IIB; B = IA, IIIA, IVA, VA, and IIB, VB, VIII; X = IIIA, VA, VIA, and VIIA) compounds with diamond-like structures.

respectively, Edef is the deformation potential of the band edge, and G is the Young’s modulus. This method was applied in our recent HTP work of Xi et al.,20 where a good balance between efficiency and accuracy was obtained. As discussed in detail below, we found conductive networks at the VBM formed from the anion sublattices. The deformation potential Edef in this work was adopted from the calculated values of six compounds, CuInS2, CuInSe2, CuInTe2, MgSiP2, MgSiAs2, and MgSiSb2, for p-type properties and of the 41 investigated compounds for n-type properties (Tables S1 and S2).20 The Young’s moduli G, which enters the HTP calculation of the relaxation time (eq 1) for the investigated compounds (Table S4) was calculated by the Voigt-ReussHill (VRH) theory.51 Lattice thermal conductivities were calculated through the Slack model52 κL = A ·

M̅ θD3δ γ 2n2/3T

(2)

where M̅ , θD, δ, γ, and n are the average mass per atom in the crystal, the Debye temperature, the volume per atom, the Grüneisen parameter, and the number of atoms in the primitive unit cell, respectively. A is a constant determined by γ. The values are shown in Table S4. All the necessary parameters for κL values can be obtained from elastic constant calculations.53,54 The Slack model has been found to be suitable for diamond-like compounds.52,55 The electronic thermal conductivity is calculated by Boltzmann transport theory, as is the electrical conductivity and Seebeck coefficient.56

3. RESULTS AND DISCUSSION The ternary compounds ABX2 (A = IA, IIA, IIIA, VA, IB, and IIB; B = IA, IIIA, IVA, VA, IIB, VB, and VIII; X = IIIA, VA, VIA, and VIIA) with diamond-like structures can be derived from the zinc blende structure. Figure 1a shows the crystal structures of

2. COMPUTATIONAL METHODS The underlying density functional theory (DFT) calculations related to structure and elastic constants were carried out using the projector augmented wave (PAW) method, as implemented in the Vienna Ab initio Simulation Package (VASP).37,38 The Perdew−Burke−Ernzerhof (PBE) generalized gradient approximation (GGA) was used as the exchange-correlation functional.39A plane-wave energy cutoff of 400 eV and an energy convergence criterion of 10−4 eV for self-consistency were adopted. All the atomic positions were fully relaxed until the calculated Hellmann−Feynman force on each atom was less than 10−4 eV/Å. The electronic structures and transport properties were calculated by using the modified Becke-Johnson (MBJ) meta-GGA potential, and the Hubbard-U correction with Dudarev’s approach was used. This approach has been shown to be effective in other diamondlike TE materials.40−44 Chemical-bonding information was obtained using the band-resolved projected crystal orbital Hamilton populations (COHPs) as implemented in the Local Orbital Basis Suite Towards Electronic-Structure Reconstruction (LOBSTER) package.43−47 Electrical transport properties were calculated based on the Boltzmann transport theory.42 The relaxation time was treated by the deformation potential method and the constant electron−phonon coupling approximation, as implemented in Transoptic.20,48−50

τn−, k1 =

2 2πkBTEdef V ℏG

∑ δ(εn, k − εn′ , k′) n ′ , k′

(1) Figure 1. (a) Crystal structures of zinc blende and ternary compounds ABX2 with diamond-like structures. (b) The high-performance thermoelectric material screening workflow.

2 2πkBTEdef

where Cn,k = is the parameter relating to the electron−phonon V ℏG coupling matrix; n and k are the band indices and wavevector, B

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

Figure 2. (a) Band structures and density of states of ZnSnSb2 and CuInTe2. (b) Maximum p-type electronic fitness function (EFF) with respect to the carrier concentration for the maximum value. (c) The band-resolved projected crystal orbital Hamilton populations (COHPs) for Zn−Sb, Sb−Sb, and Sb onsite of ZnSnSb2 as well as Cu−Te, Te−Te,and Te onsite of CuInTe2.

we split the cation site. Differences in the atomic radii and charge states of the constituent elements lead to finite but generally small lattice distortions in the ABX2 structure. These can be described by a lattice distortion parameter η = c/2a, where a and c are the tetragonal lattice parameters.

zinc blende and ABX2. The zinc blende structure belongs to space group 216 (F4̅3m) and can be split to two distinct facecentered cubic (FCC) sublattices. Similarly, the ABX2 diamondlike compounds are based on a further splitting of one site in the zinc blende structure leading to space group of 122 (I4̅2d). Here C

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

EFFs are for ZnGeP2 and ZnGeAs2, which are also predicted to be very high-power-factor compounds. In order to further rationalize the chemical interaction between the constituent elements, we calculated the bandresolved projected COHPs47 for the ABX2 compounds with band gaps between 0.1 and 2.5 eV. Figure 2c shows the bandresolved projected COHPs for Zn−Sb, Sb−Sb, and Sb onsite (atomic-like) of ZnSnSb2 as well as Cu−Te, Te−Te,and Te onsite of CuInTe2. The results for the interatomic interactions in the ABX2 compounds investigated are shown in Figure S3. Positive and negative values correspond to the antibonding state and bonding states, respectively, and zero represents nonbonding. The VBM for ZnSnSb2 is mainly composed of the Sb onsite states, whereas the other interatomic interactions are nonbonding states. In the representative Cu containing compound, such as CuInTe2, the Cu−Te antibonding states are non-negligible. Following analysis of the band-resolved projected COHPs for the 41 ABX2 compounds, we conclude that for pnictides, the VBM is mainly composed of X onsite states or nonbonding states. For Cu-based chalcogenides, the major interactions are composed of Cu−X antibonding states and X onsite states. Following the analysis and procedure of ref 20, we find that the VBMs of the 41 ABX2 compounds investigated here all show a conductive network dominated by the anion sublattice. A conductive network indicates that in one compound, certain atoms or atomic groups dominate the electrical transport properties, whereas other atoms are out of the conductive network and serve as carrier reservoirs. The carrier concentration can thus be optimized by doping on the out-of-network atoms, which minimize the alteration on the electronic structure of the base compound. As shown in ref 20, nonbonding state compounds have slightly more delocalized wave functions and thus more beneficial conductive networks than those of antibonding states. Therefore, we can conclude that pnictides may have better p-type electrical transport properties than chalcogenides. The anion sublattice conductive network is reflected in the electrical transport properties. Figure 3a,b shows the Pisarenko relations for the p-type P-, As-, and Sb-based as well as S-, Se-, and Te-based ABX2 compounds at 700 K. The Pisarenko relations can be roughly categorized according to the different anions, indicating again the importance of the anion sublattice conductive networks. The calculated Seebeck effective masses from the Pisarenko analysis are shown in Figure 3a,b and are 2.93, 2.10, and 1.52 me for P-, As-, and Sb-based compounds and 4.59, 3.97, and 2.82 me for S-, Se-, and Te-based compounds. The relative trend of effective masses from 3p to 5p anions21 holds in the pnictide compounds. Figure 3c,d shows the theoretical electrical conductivities as a function of carrier concentrations at 700 K. As a test, the experimental electrical conductivity of Cu2SnS3 is labeled in Figure 3d.60 It can be seen that the experimental value is close to the theoretical value, which supports the accuracy of the methodology. The electrical conductivities are classified by anion sublattices; the electrical conductivities with the same anion sublattice have some deviations. This is understandable, because the electrical conductivities are more sensitive to the subtle variations of the band shapes. It is this property that is measured by the EFF and that allows the effective mass for conductivity to differ from the Pisarenko effective mass and decouple the usual contrary relationship of S and σ on mass. At the same carrier concentration range, P-, As-, and Sb-based compounds have higher electrical conductivities than those of

Figure 1b shows the HTP workflow in this study. We applied three search criteria on the MIP: (1) No. 122 space group; (2) eight atoms per unit cell; (3) atomic ratio of 1:1:2. These criteria selected 65 entries out of 84 908 in the MIP. The first step following this selection is band gap screening. Becuase TE materials are mostly small band gap semiconductors, compounds with the calculated band gap in the range of 0.1−2.5 eV are selected. A total of 41 compounds out of 65 meet this criterion. The band-resolved projected COHP is used to obtain the chemical-bonding characteristics. The p-type and n-type TE properties at 700 K are then obtained. The electronic fitness function (EFF), which is a measure of the extent to which a band structure can resolve the counter-correlation of S and σ (e.g., through complex nonparabolic bands) is calculated.50 Finally, high-performance TE compounds are predicted. Figure 2a shows the band structures and density of states (DOS) of two compounds, ZnSnSb2 and CuInTe2. Both of the VBMs are triply degenerate. These bands are composed of Sb p states in ZnSnSb2 and Cu d and Te p states in CuInTe2. The electronic structures of all the ABX2 compounds are shown in Figure S1. More than 90% of compounds have triple degenerate bands at the VBM. For Cu-based compounds, the VBMs are mainly composed of Cu d and X p states, whereas for the other compounds, the VBM has X p character. Although the band characters at the VBM of the two compounds in Figure 2a are different, their band shapes are similar. This implies that pnictide ABX2 compounds may have a similar conductive network as their chalcogenide counterparts. Furthermore, the VBM of ZnSnSb2 is more dispersive than that of CuInTe2, implying a lower band edge DOS. The pnictides in general show a smaller DOS near the VBM than the chalcogenides do, as seen in Figure S1. This is important for the electrical transport properties, as shown below. On the basis of Ioffe’s theory, the optimal band gaps for TE applications are around 10 kBTop, where Top is the operating temperature.57 At 700 K, the optimal band gap is around 0.7 eV. Considering the uncertainty of the band gap calculation in density functional theory, here we use the band gap range of 0.1∼2.5 eV for the further TE properties. We calculated the phonon spectra of these 41 compounds as shown in Figure S2. All the phonon spectra of these compounds have no imaginary frequencies, so these compounds are considered stable. Figure 2b shows the relationship between the maximum ptype electronic fitness function (EFF)58,59 at 700 K with respect to the carrier concentrations, categorized by the anion sublattice. The EFF identifies the complexity of electronic structures and was recently proposed by Singh and co-workers.58,59 Generally speaking, the higher the value, the better the electrical transport properties of the compound. The EFFs for the chalcogenides are within the range of 3∼7*10−20 W5/3 m s−1/3 K−2 (the EFF values can be seen in Table S3). Interestingly, these values are comparable to those of the sulfide compounds in ref 58, though the latter are not limited to the diamond-like structures. The pnictides, on the other hand, have distinctly larger EFFs than those of chalcogenides by up to 100%. This indicates that from a pure band structure point of view, the pnictides are generally superior to the chalcogenides. The three compounds with the highest EFFs are BeSiP2, BeSiAs2, and BeSiSb2. Interestingly, these are also the compounds with the highest predicted power factors, showing the importance of band structure to the TE properties. Unfortunately, these compounds may be unsuitable due to the presence of Be, which is both light (unfavorable for low κL, especially combined with Si) and toxic. The next highest D

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

Figure 4. (a) Calculated and experimental temperature dependent lattice thermal conductivities for CuInTe2 and AgInSe2. (b) Relationship between calculated lattice thermal conductivities (700 K) and mean atomic weights with different anion sublattice. (c,d) Predicted maximum ZT values for 41 p-type compounds at 700 K; the compounds with ZT over 0.6 for p-type and n-type are labeled.

Figure 4b shows the relationship between the κLs at 700 K and the mean atomic weights of 41 ABX2 compounds. Generally speaking, heavy constituents lead to low κLs. It is notable that κLs for the pnictide compounds are usually higher than those for the chalcogenide compounds. This is due to the generally higher elastic moduli for the pnictides (Table S4). Figure 4c,d shows the predicted maximum ZT values by using the calculated power factors and thermal conductivities. Here, total thermal conductivities, both lattice thermal conductivities and electron thermal conductivities, are considered for ZT value prediction. Figure 4c shows the predicted maximum ZT values as a relationship of corresponding carrier concentrations at 700 K for 41 p-type ABX2 compounds. Despite the high PFs, pnictides are comparable with the chalcogenide counterparts in ZTmax because of the cancellation from enhanced κes and κLs. In Figure 4c, 12 compounds with ZT > 0.6 are labeled. CuInTe2 and AgGaTe2 have been experimentally confirmed to be highperformance p-type TE materials, and their experiment ZT values at 700 K are 0.59 and 0.5, respectively, which are consistent with our theoretical results (0.66 and 0.67).17,19 Nomura et al. recently reported that ZnSnSb2 is a promising TE material.32 The other compounds, predicted TE materials CdSnSb2, AgInS2, AgGaSe2, ZnCdTe2, and LiInTe2, have been little studied in this context experimentally, but based on the results above, they are promising TE materials. By the same method, the predicted maximum ZT values of n-type ABX2 compounds at 700 K are shown in Figure 4d, and 12 compounds with ZT > 0.6 are labeled. Among them, ZnSnSb2, AgInS2, AgGaSe2, AgInSe2, and LiInTe2 have high ZT values for both ntype and p-type, which may result in good TE application prospects.

Figure 3. (a,b) Pisarenko plots for P-, As-, and Sb-based as well as S-, Se-, and Te-based compounds. (c,d) Calculated electrical conductivities for P-, As-, and Sb-based as well as S-, Se-, ans Te-based compounds at 700 K. Experimental value of Cu2SnS3 at 700 K is marked in pink.60 (e) The correlations between ΔCF and η. (f) The correlations between power factors and ΔCF.

S-, Se-, and Te-based compounds. The lighter VBMs of pnictides cause larger group velocities (lower transport effective mass) and less electron−phonon scattering phase space, both of which are beneficial to the electrical conductivities and PFs. As shown in Figure 3f, pnictides usually have higher PFs up to 100 μW/cm·K2, much larger than those of chalcogenides (up to 30 μW/cm·K2). The electronic structures and electrical transport properties indicate that there exists a universal conductive network dominated by the anion sublattice in these ABX2 compounds. For the compounds with the same anion sublattice, there are also minor differences in their electronic structures (e.g., the crystal field splitting energy ΔCF (= E(Γ5v) − E(Γ4v), MgSiP2 in Figure S1) at the VBM.16 As shown in Figure S1, the three bands at the VBM are split in some ABX2 compounds (i.e., ΔCF ≠ 0). Figure 3e shows a good linear relationship between ΔCF and the distortion parameter η. When η is close to unity, ΔCF approaches zero, which agrees well with the “the unity-η rule”.16 Figure 3f shows the relationship between calculated PF (with the same carrier concentrations corresponding to the maximum ZT values) at 700 K and ΔCF for the 41 ABX2 compounds. In the chalcogenides, PF is almost unchanged with ΔCF, whereas in the pnictides, there is a rough positive correlation between PF and ΔCF. As mentioned, the lattice thermal conductivities are obtained using the Slack model.58 Figure 4a shows the calculated as well as experimental lattice thermal conductivities at 300−900 K for CuInTe2 and AgInSe2.61,62 The experimental values agree well with theoretical values especially at high temperatures, which indicates that the Slack model is suitable for thermal conductivity prediction in these diamond-like compounds.

4. CONCLUSION An HTP workflow for high-TE-performance screening has been performed for ABX2 compounds with diamond-like structures based on the MIP, and results have been analyzed to extract trends and physical understanding. We find that there is a universal conductive network dominated by the anion sublattice E

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

edges the support program of the Shanghai Subject Chief Scientist (16XD1401100), the support from Guangdong I n n o v a t i o n R e s e a r c h T ea m P r o j e c t (G r a n t N o . 2017ZT07C062), and the Shenzhen Pengcheng-Scholarship program. Work at the University of Missouri is supported by the Department of Energy, Basic Energy Sciences, grant number DE-SC0019114.

in both pnictides and chalcogenides. The VBMs of these compounds are triply degenerate (neglecting spin-orbit) or slightly split bands and are mainly composed of antibonding or nonbonding states. The lattice distortion has a minor influence on the band structures for the same anion sublattice. The pnictides have relatively larger PFs than the chalcogenides because of the lighter bands near the VBM. This also leads to lower optimal carrier concentrations. The lattice thermal conductivities are in general negatively correlated with the mean atomic weights of the compounds. In addition, the pnictides have somewhat higher κLs than the chalcogenides, partly compensating for the higher power factors. Finally, 12 promising high-performance TE materials of both n-type and ptype are predicted. Some of these 12 (CuInTe2, AgGaTe2, and ZnSnSb2) have been found experimentally to be good TE materials, supporting the approach used. Some of the predicted compounds, ZnSnSb2, AgInS2, AgGaSe2, AgInSe2, and LiInTe2, have high ZT values for both n-type and p-type, which may have good TE application prospects. The HTP workflow in this study shows the effectiveness in screening novel functional materials.





(1) Goldsmid, H. J. Thermoelectric Refrigeration, Electrical Engineering, Springer: Boston, MA, 2013. (2) Sales, B. C. Smaller is Cooler. Science 2002, 295 (5558), 1248− 1249. (3) Rowe, D. M., Ed. Thermoelectrics Handbook: Macro to Nano; CRC Press, 2005. (4) Liu, W.; Yan, X.; Chen, G.; Ren, Z. Recent Advances in Thermoelectric Nanocomposites. Nano Energy 2012, 1 (1), 42−56. (5) 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 Theory−experiment Perspective. npj Comput. Mater. 2016, 2, 15015. (6) Zhu, T.; Liu, Y.; Fu, C.; Heremans, J. P.; Snyder, J. G.; Zhao, X. Compromise and Synergy in High-efficiency Thermoelectric materials. Adv. Mater. 2017, 29 (14), 1605884. (7) Ni, W. X.; Chen, W. M.; Buyanova, I. A.; Henry, A.; Hansson, G. V.; Monemar, B. Some Critical Issues on Growth of High Quality Si and SiGe Films Using a Solid-source Molecular Beam Epitaxy System. J. Cryst. Growth 1995, 157 (1−4), 242−247. (8) Heremans, J. P.; Jovovic, V.; Toberer, E. S.; Saramat, A.; Kurosaki, K.; Charoenphakdee, A.; Yamanaka, S.; Snyder, G. J. Enhancement of Thermoelectric Efficiency in PbTe by Distortion of the Electronic Density of States. Science 2008, 321 (5888), 554−557. (9) Shi, X.; Yang, J.; Salvador, J. R.; Chi, M.; Cho, J. Y.; Wang, H.; Bai, S.; Yang, J.; Zhang, W.; Chen, L. Multiple-filled Skutterudites: High Thermoelectric Figure of Merit Through Separately Optimizing Electrical and Thermal Transports. J. Am. Chem. Soc. 2011, 133 (20), 7837−7846. (10) Nolas, G. S., Ed. The Physics and Chemistry of Inorganic Clathrates; Springer Series in Materials Science; Springer, 2014; Vol. 199. (11) Yang, J.; Li, H. M.; Wu, T.; Zhang, W. Q.; Chen, L. D.; Yang, J. H. Evaluation of Half-Heusler Compounds as Thermoelectric Materials Based on the Calculated Electrical Transport Properties. Adv. Funct. Mater. 2008, 18 (19), 2880−2888. (12) Skoug, E. J.; Cain, J. D.; Morelli, D. T. High Thermoelectric Figure of Merit in the Cu3SbSe4-Cu3SbS4 Solid Solution. Appl. Phys. Lett. 2011, 98 (26), 261911. (13) Cheng, N.; Liu, R.; Bai, S.; Shi, X.; Chen, L. Enhanced Thermoelectric Performance in Cd doped CuInTe2 Compounds. J. Appl. Phys. 2014, 115 (16), 163705. (14) Zeier, W. G.; Zhu, H.; Gibbs, Z. M.; Ceder, G.; Tremel, W.; Snyder, G. J. Band Convergence in the Non-cubic Chalcopyrite Compounds Cu2MGeSe4. J. Mater. Chem. C 2014, 2 (47), 10189− 10194. (15) Qin, Y.; Qiu, P.; Liu, R.; Li, Y.; Hao, F.; Zhang, T.; Ren, D.; Shi, X.; Chen, L. Optimized Thermoelectric Properties in Pseudocubic Diamond-like CuGaTe2 Compounds. J. Mater. Chem. A 2016, 4 (4), 1277−1289. (16) Zhang, J.; Liu, R.; Cheng, N.; Zhang, Y.; Yang, J.; Uher, C.; Shi, X.; Chen, L.; Zhang, W. High-performance Pseudocubic Thermoelectric Materials from Non-Cubic Chalcopyrite Compounds. Adv. Mater. 2014, 26 (23), 3848−3853. (17) Yusufu, A.; Kurosaki, K.; Kosuga, A.; Sugahara, T.; Ohishi, Y.; Muta, H.; Yamanaka, S. Thermoelectric Properties of Ag1−xGaTe2 with Chalcopyrite Structure. Appl. Phys. Lett. 2011, 99 (6), No. 061902. (18) Plirdpring, T.; Kurosaki, K.; Kosuga, A.; Day, T.; Firdosy, S.; Ravi, V.; Snyder, G. J.; Harnwunggmoung, A.; Sugahara, T.; Ohishi, Y.; et al. Chalcopyrite CuGaTe2: a High-efficiency Bulk Thermoelectric Material. Adv. Mater. 2012, 24 (27), 3622−3626.

ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsami.9b01196.



REFERENCES

The deformation potential (Edef) values for p-type and ntype ABX2 compounds (X = S, Se, Te, P, As, Sb) are shown in Tables S1 and S2. The distortion parameters η, crystal field splitting energies ΔCF, band gaps Eg, maximum p-type electronic fitness functions, and the corresponding carrier concentrations p for 41 compounds are listed in Table S3, respectively. Lattice thermal conductivities κL, constants A, mean atomic weights, Young’s moduli G, predicted maximum ZT values and the corresponding carrier concentrations, Seebeck coefficients, and power factors PFs of the 41 investigated compounds at 700 K for both n-type and p-type are shown in Table S4. Figure S1 shows band structures and density of states (DOS) at 700 K for 65 ABX2 compounds. Figure S2 shows the phonon spectra of the 41 compounds. Figure S3 shows band-resolved projected crystal orbital Hamilton populations (COHPs) (A−A, A−X, B−X, X onsite) for 41 compounds (PDF)

AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. (L.X.) *E-mail: [email protected]. (J.Y.) *E-mail: [email protected]. (W.Z.) ORCID

Jiong Yang: 0000-0002-5862-5981 David J. Singh: 0000-0001-7750-1485 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the National Key Research and Development Program of China (Nos. 2017YFB0701600 and 2018YFB0703600), the Natural Science Foundation of China (Grant Nos. 51572167, 51632005, 11574333, 11674211, and 51761135127), and the 111 Project D16002. W.Q.Z. acknowlF

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

(37) Kresse, G.; Joubert, D. From ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59 (3), 1758−1775. (38) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54 (16), 11169. (39) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77 (18), 3865− 3868. (40) Singh, D. J. Electronic Structure Calculations with the TranBlaha Modified Becke-Johnson Density Functional. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 82 (20), 205102. (41) Wang, S.; Wang, Z.; Setyawan, W.; Mingo, N.; Curtarolo, S. Assessing the Thermoelectric Properties of Sintered Compounds via High-Throughput Ab-Initio Calculations. Phys. Rev. X 2011, 1 (2), 021012. (42) Zhang, Y.; Zhang, J.; Gao, W.; Abtew, T. A.; Wang, Y.; Zhang, P.; Zhang, W. Near-edge Band Structures and Band Gaps of Cu-based Semiconductors Predicted by The Modified Becke-Johnson Potential plus an on-site Coulomb U. J. Chem. Phys. 2013, 139 (18), 184706. (43) Maintz, S.; Deringer, V. L.; Tchougreeff, A. L.; Dronskowski, R. Analytic Projection From Plane-Wave and PAW Wavefunctions and Application to Chemical-Bonding Analysis in Solids. J. Comput. Chem. 2013, 34 (29), 2557−2567. (44) Dronskowski, R.; Blöchl, P. E. Crystal Orbital Hamilton Populations (COHP): Energy-resolved Visualization of Chemical Bonding in Solids Based on Density-Functional Calculations. J. Phys. Chem. 1993, 97 (33), 8617−8624. (45) Deringer, V. L.; Tchougréeff, A. L.; Dronskowski, R. Crystal Orbital Hamilton Population (COHP) Analysis as Projected From Plane-wave Basis sets. J. Phys. Chem. A 2011, 115 (21), 5461−5466. (46) Maintz, S.; Deringer, V. L.; Tchougréeff, A. L.; Dronskowski, R. LOBSTER: A Tool to Extract Chemical Bonding from Plane-Wave Based DFT. J. Comput. Chem. 2016, 37 (11), 1030−1035. (47) Sun, X.; Li, X.; Yang, J.; Xi, J.; Nelson, R.; Ertural, C.; Dronskowski, R.; Liu, W.; Snyder, G. J.; Singh, D. J.; Zhang, W. Achieving Band Convergence by Tuning the Bonding Ionicity in n-type Mg3Sb2. J. Comput. Chem. 2019. (48) Transoptic. http://www.mgi.shu.edu.cn/Portals/675/transoptic. zip. (49) Yang, J.; Xi, L.; Zhang, W.; Chen, L. D.; Yang, J. Electrical Transport Properties of Filled CoSb3 Skutterudites: A Theoretical Study. J. Electron. Mater. 2009, 38 (7), 1397−1401. (50) Chen, J.; Wang, D.; Shuai, Z. First-Principles Predictions of Thermoelectric Figure of Merit for Organic Materials: Deformation Potential Approximation. J. Chem. Theory Comput. 2012, 8 (9), 3338− 3347. (51) Toonder, J. M. J.; van Dommelen, J. A. W.; Baaijens, F. P. T. The Relation Between Single Crystal Elasticity and the Effective Elastic Behaviour of Polycrystalline Materials: Theory, Measurement and Computation. Modell. Simul. Mater. Sci. Eng. 1999, 7 (6), 909−928. (52) Slack, G. A. Nonmetallic Crystals With High Thermal Conductivity. J. Phys. Chem. Solids 1973, 34 (2), 321−335. (53) Sun, L.; Liu, B.; Wang, J.; Li, Z.; Wang, J. Theoretical Study on the Relationship Between Crystal Chemistry and Properties of Quaternary Y−Si−O−N Oxynitrides. J. Am. Ceram. Soc. 2016, 99 (7), 2442−2450. (54) Liu, B.; Cooper, V. R.; Wang, H.; Xiang, Q.; Li, J.; Gao, Y.; Yang, J.; Zhou, Y.; Liu, B. Discovery of ABO3 Peroskites as Thermal Barrier Coatings Through High-throghput First Principles Calculations. Mater. Res. Lett. 2019, 7 (4), 145−151. (55) Jia, T.; Chen, G.; Zhang, Y. Lattice Thermal Conductivity Evaluated Using Elastic Properties. Phys. Rev. B: Condens. Matter Mater. Phys. 2017, 95 (15), 155206. (56) Madsen, G. K.; Singh, D. J. BoltzTraP. A Code for Calculating Band-structure Dependent Quantities. Comput. Phys. Commun. 2006, 175 (1), 67−71. (57) Ioffe, A. F. Semiconductor Thermoelements and Thermoelectric Cooling, Physics Today: London, UK, 1957.

(19) Liu, R.; Xi, L.; Liu, H.; Shi, X.; Zhang, W.; Chen, L. Ternary compound CuInTe2: a promising Thermoelectric material with Diamond-Like Structure. Chem. Commun. 2012, 48 (32), 3818−3820. (20) Xi, L.; Pan, S.; Li, X.; Xu, Y.; Ni, J.; Sun, X.; Yang, J.; Luo, J.; Xi, J.; Zhu, W.; Li, X.; Jiang, D.; Dronskowski, R.; Shi, X.; Snyder, G. J.; Zhang, W. Discovery of High-Performance Thermoelectric Chalcogenides through Reliable High-Throughput Material Screening. J. Am. Chem. Soc. 2018, 140 (34), 10785−10793. (21) Shi, X.; Xi, L.; Fan, J.; Zhang, W.; Chen, L. Cu−Se Bond Network and Thermoelectric Compounds with Complex Diamondlike Structure. Chem. Mater. 2010, 22 (22), 6029−6031. (22) Xi, L.; Zhang, Y. B.; Shi, X. Y.; Yang, J.; Shi, X.; Chen, L. D.; Zhang, W.; Yang, J.; Singh, D. J. Chemical Bonding, Conductive Network, and Thermoelectric Performance of the Ternary Semiconductors Cu2SnX3 (X = Se, S) from First Principles. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86 (15), 155201. (23) Zhu, H.; He, R.; Mao, J.; Zhu, Q.; Li, C.; Sun, J.; Ren, W.; Wang, Y.; Liu, Z.; Tang, Z.; et al. Discovery of ZrCoBi Based half Heuslers with High Thermoelectric Conversion Efficiency. Nat. Commun. 2018, 9 (1), 2497. (24) Curtarolo, S.; Setyawan, W.; Hart, G. L. W.; Jahnatek, M.; Chepulskii, R. V.; Taylor, R. H.; Wang, S.; Xue, J.; Yang, K.; Levy, O.; Mehl, M. J.; Stokes, H. T.; Demchenko, D. O.; Morgan, D. AFLOW: An Automatic Framework for High-Throughput Materials Discovery. Comput. Mater. Sci. 2012, 58, 218−226. (25) Carrete, J.; Li, W.; Mingo, N.; Wang, S.; Curtarolo, S. Finding Unprecedentedly Low-Thermal-Conductivity Half-Heusler Semiconductors via High-Throughput Materials Modeling. Phys. Rev. X 2014, 4 (1), 011019. (26) The Materials Project. https://www.materialsproject.org/. (27) Zhu, H.; Hautier, G.; Aydemir, U.; Gibbs, Z. M.; Li, G.; Bajaj, S.; Poehls, J.-H.; Broberg, D.; Chen, W.; Jain, A.; White, M. A.; Asta, M.; Snyder, G. J.; Persson, K.; Ceder, G. Computational and Experimental Investigation of TmAgTe2 and XYZ2 Compounds, a New Group of Thermoelectric Materials Identified by First-principles High-throughput Screening. J. Mater. Chem. C 2015, 3 (40), 10554−10565. (28) Aydemir, U.; Poehls, J.-H.; Zhu, H.; Hautier, G.; Bajaj, S.; Gibbs, Z. M.; Chen, W.; Li, G.; Ohno, S.; Broberg, D.; Kang, S. D.; Asta, M.; Ceder, G.; White, M. A.; Persson, K.; Jain, A.; Snyder, G. J. YCuTe2: A Member of a New Class of Thermoelectric Materials with CuTe4-based Layered Structure. J. Mater. Chem. A 2016, 4 (7), 2461−2472. (29) Chen, W.; Pohls, J.-H.; Hautier, G.; Broberg, D.; Bajaj, S.; Aydemir, U.; Gibbs, Z. M.; Zhu, H.; Asta, M.; Snyder, G. J.; Meredig, B.; White, M. A.; Persson, K.; Jain, A. Understanding Thermoelectric properties from High-throughput Calculations: Trends, Insights, and Comparisons with Experiment. J. Mater. Chem. C 2016, 4 (20), 4414− 4426. (30) Materials Informatics Platform (MIP). http://mip.shu.edu.cn. (31) Kocak, B.; Ciftci, Y.; Surucu, G. Structural and Thermoelectronic Properties of Chalcopyrite MgSiX2 (X= P, As, Sb). J. Electron. Mater. 2017, 46 (1), 247−264. (32) Nomura, A.; Choi, S.; Ishimaru, M.; Kosuga, A.; Chasapis, T.; Ohno, S.; Snyder, G. T.; Ohishi, Y.; Muta, H.; Yamanaka, S.; Kurosaki, K. Chalcopyrite ZnSnSb2: A Promising Thermoelectric Material. ACS Appl. Mater. Interfaces 2018, 10 (50), 43682−43690. (33) Jaffe, J.; Zunger, A. Electronic Structure of the Ternary Pnictide Semiconductors ZnSiP2, ZnGeP2, ZnSnP2, ZnSiAs2, and MgSiP2. Phys. Rev. B: Condens. Matter Mater. Phys. 1984, 30 (2), 741. (34) Xiao, J.; Zhu, S.; Zhao, B.; Chen, B.; Liu, H.; He, Z. Computational Assessment of Promising Mid-infrared Nonlinear Optical Materials Mg−IV−V2 (IV= Si, Ge, Sn; V= P, As): A Firstprinciples Study. Mater. Res. Express 2018, 5 (3), No. 035907. (35) Ouahrani, T.; Ö ztekin Ç iftci, Y.; Mebrouki, M. Dynamical and Anisotropic Behavior of the MSiP2 (M= Be, Mg, Cd, Zn and Hg) compounds. J. Alloys Compd. 2014, 610, 372−381. (36) Fahad, S.; Murtaza, G.; Ouahrani, T.; Khenata, R.; Yousaf, M.; Omran, S. B.; Mohammad, S. Structural, Elastic, Electronic, Bonding, and Optical Properties of BeAZ2 (A= Si, Ge, Sn; Z= P, As) Chalcopyrites. J. Alloys Compd. 2015, 646, 211−222. G

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces (58) Xing, G.; Sun, J.; Li, Y.; Fan, X.; Zheng, W.; Singh, D. J. Electronic Fitness Function for Screening Semiconductors as Thermoelectric Materials. Phys. Rev. Mater. 2017, 1 (6), No. 065405. (59) He, X.; Singh, D. J.; Boon-on, P.; Lee, M.-W.; Zhang, L. Dielectric Behavior as a Screen in Rational Searches for Electronic Materials: Metal Pnictide Sulfosalts. J. Am. Chem. Soc. 2018, 140 (51), 18058− 18065. (60) Tan, Q.; Sun, W.; Li, Z.; Li, J.-F. Enhanced Thermoelectric Properties of Earth-abundant Cu2SnS3 via In Doping Effect. J. Alloys Compd. 2016, 672, 558−563. (61) Qiu, P.; Qin, Y.; Zhang, Q.; Li, R.; Yang, J.; Song, Q.; Tang, Y.; Bai, S.; Shi, X.; Chen, L. Intrinsically High Thermoelectric Performance in AgInSe2 n-Type Diamond-Like Compounds. Adv. Sci. 2018, 5 (3), 1700727. (62) Luo, Y.; Yang, J.; Jiang, Q.; Li, W.; Zhang, D.; Zhou, Z.; Cheng, Y.; Ren, Y.; He, X. Progressive Regulation of Electrical and Thermal Transport Properties to High-Performance CuInTe2 Thermoelectric Materials. Adv. Energy. Mater. 2016, 6 (12), 1600007.

H

DOI: 10.1021/acsami.9b01196 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX