Bi2PdO4: A Promising Thermoelectric Oxide with High Power Factor

Dec 23, 2016 - ABSTRACT: The search for new energy harvesting materials that directly convert (waste) heat into electricity has received increasing at...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/cm

Bi2PdO4: A Promising Thermoelectric Oxide with High Power Factor and Low Lattice Thermal Conductivity Jiangang He,† Shiqiang Hao,† Yi Xia,‡ S. Shahab Naghavi,† Vidvuds Ozoliņs,̌ ‡ and Chris Wolverton*,† †

Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, United States Department of Materials Science and Engineering, University of California, Los Angeles, California 90095, United States



S Supporting Information *

ABSTRACT: The search for new energy harvesting materials that directly convert (waste) heat into electricity has received increasing attention. Transition metal oxides are a promising class of thermoelectric (TE) materials that can operate at high temperature due to their chemical and thermal stability. However, the high lattice thermal conductivity, poor electrical conductivity, and low thermopower have significantly impeded their applications to date. Using first-principles calculations, we predict a known oxide Bi2PdO4 to be a highly efficient holedoped TE material with low lattice thermal conductivity and high power factor. These properties are due to (i) the strong anharmonicity stemming from Bi3+ 6s2 lone pair electrons (leading to low lattice thermal conductivity) and (ii) the flat-anddispersive valence band structure with high band degeneracy originating from the localized Pd2+ dz2 orbitals in the stacked square planar ligand field (leading to a large power factor). Our results highlight the possibility of oxides as potential TE materials and also afford a novel strategy of designing TE materials by synthesizing compounds which combine a lone pair active cation with a d8 cation in a stacked square planar ligand field.

T

clathrates,16 R-Heusler,17 and CsAg5Te318), anharmonicity (Cu12Sb4S1319 and SnSe4), and lone pair electrons.20,21 However, approaches for optimizing electrical performances are relatively limited due the complexity of the electronic structure. Recently, Bilc et al.22 demonstrated that the combination of flat (narrow carrier energy distribution, large effective mass) and dispersive (small effective mass) bands yields a high PF since S is proportional to the effective mass and σ is inversely proportional to the effective mass. Although the flat-and-dispersive band structure, which is energetically dispersive in some regions of the Brillouin zone and flat in others, is rare in materials, it has been observed in some full Heusler compounds.22 Oxides in general have high stability in air at high temperatures.23,24 Many oxides also are lightweight and environmentally friendly, which are desirable for commercial applications of energy harvesting devices. However, oxide thermoelectrics are often limited by high κL and low σ,25 leading to low ZT.26 Oxides with low κL have so far only been observed in complex crystal structures with a large number of atoms per primitive unit cell and low crystal symmetry, yielding a more tortuous path for phonon transport.6 A commonly used

hermoelectric materials, which directly interconvert heat and electricity, have important applications in the recovery of waste heat and refrigeration. However, although a long-term effort has been given to new TE material exploration and performance optimization,1 there are only very few materials with commercial applications. Commercial TE materials require high heat-electric conversion efficiency at a wide range of temperatures, stability (thermal, chemical, and mechanical) at high temperature and in air, and low cost. The figure of merit ZT of TE materials is defined as ZT = (S2σT)/ (κL + κe), where S, σ, κL, κe, and T are the Seebeck coefficient, electrical conductivity, lattice thermal conductivity, electron thermal conductivity, and temperature, respectively. Therefore, a high ZT material must possess the following combination of properties: high S and σ and low κL and κe. The power factor (PF, S2σ) of a compound mainly depends on its electronic structure: a large dispersive band (small effective mass) around the Fermi level gives rise to a large σ, while a sharp increase in the density of states (large effective mass) usually leads to large S.2 Owing to these requirements on material properties, TE materials with high PF and low κ (κL + κe) have, to date, mainly been discovered in the chalcogenide family, such as Bi2Te3, PbTe, SnSe, and AgSbTe2.3 For example, the ZT of the undoped SnSe bulk was reported to be 2.6 at 923 K,4 while Na doped SnSe has a high ZT in a broad temperature range.5 So far, many phonon-engineering strategies have been employed for reducing κL in semiconductors:6 nanostructuring (PbTe,7−11 half-Heusler12−14), rattling modes (skutterudites,15 © 2016 American Chemical Society

Special Issue: Computational Design of Functional Materials Received: October 4, 2016 Revised: December 22, 2016 Published: December 23, 2016 2529

DOI: 10.1021/acs.chemmater.6b04230 Chem. Mater. 2017, 29, 2529−2534

Article

Chemistry of Materials strategy of reducing κL in oxides is to employ the interfaces of layered oxides to scatter heat carrying phonons, such as NaCoO2,27 Ruddleson−Popper structures,28 and artificial superlattices.24 The main drawback of this method is the low σ through the interface due to the large inertial effective mass m*I .29 Although a high density of states effective mass (m*d ) is favorable for S,30 the band effective mass mb* is also high in these compounds due to the single-valley band structure and 31 low band degeneracy Nv since m*d = N2/3 v m* b . Large m* b leads τ to low electrical conductivity because σ ∝ m * (τ is the electron

octahedron connects with two neighboring BiO6 octahedra neighbors by edge-sharing along the c axis and six BiO6 octahedra by corner-sharing in the ab plane. Due to this connectivity and the strong distortion of the BiO6 octahedra, sizable tunnels are formed along the [001] and [110] directions. The large asymmetric distortion of the BiO6 octahedra also indicates the possibility of strong anharmonicity, which can significantly enhance heat carrying phonon scattering and reduce heat transport correspondingly20 (we will discuss κL later). The calculated phonon dispersion is shown in Figure 2. Three acoustic modes, which are largely responsible for heat

b

lifetime). On the basis of crystal field and crystal structure chemistry analysis, in this paper, we propose an alternative design strategy: we employ an oxide with asymmetrically distorted oxygen octahedra stemming from stereoactive lone pair electrons of the cation, Bi3+, and the electronic structure with a flat-and-dispersive band and high valley degeneracy originating from a d8 Pd2+ cation in a stacked square planar ligand field. This strategy is realized in a prototype compound, Bi2PdO4, which is predicted to be a promising p-type TE material with both low κL and high PF by means of firstprinciples calculations (see Supporting Information for computational details). Bi2PdO4 crystallizes in a tetragonal crystal structure (P4/ncc, No. 130), which is isostructural with Bi2CuO4. As shown in Figure 1, Pd2+ is coordinated in a square-planar geometry by

Figure 2. Phonon dispersion (left) and density of states (right) of Bi2PdO4. Inset is the first Brillouin zone.

transport, have zone-boundary frequencies around 30 cm−1 along the Γ−X (a, b axis) and Γ−M (diagonal of a and b axis) directions and 40 cm−1 along Γ−Z (c axis). The acoustic phonons in the ab plane (Γ−X and Γ−M) are softer than these along the c axis (Γ−Z). This softening is mainly due to the larger tunnel structure in the ab plane (i.e., tunnel along [001]) and the stereochemically active lone pair electron of Bi3+; see Figure 1. As shown in the phonon density of states (DOS) in Figure 2, these low frequency phonon modes have the largest contribution from the heavy Bi atoms. The three lowest optical modes at Γ are the polar modes A2u (38 cm−1) and Eu (39 cm−1). The lowest frequency modes at X, M, and Z points of the first Brillouin zone are X1 (27 cm−1), M1 (27 cm−1), and Z3Z4 (40 cm−1), respectively. All these low-frequency acoustic modes indicate small group velocities near the zone center (ΓX: 1520, 2400, and 2450 m/s; Γ-M: 840, 1480, 3280 m/s; Γ-Z: 1450, 1450, 4270 m/s) and low Debye temperatures, which facilitate a low κL (e.g., according to the Debye−Callaway model32,33). Accurate calculation of κL requires a full determination of anharmonic force constants. We use the recently developed first-principles compressive sensing lattice dynamics method34 and the first order perturbation theory with Boltzmann transport equation (BTE); see Supporting Information. The resulting lattice thermal conductivities, κL, are shown in Figure 3 as functions of temperature. The lattice thermal conductivities yy zz along the a (κxx L , or b, κL ) and c axes (κL ) are 2.2 and 2.8 W −1 −1 m K at 300 K, respectively, which are comparable with that of the high ZT chalcogenide material PbTe and is much lower

Figure 1. Crystal structure of Bi2PdO4 in top (a) and side view (b). The structure has tunnels along the c axis. (c) is the top view of the electron localization function (ELF) with isosurface level of 0.55. (d) is the distorted octahedra of BiO6 with ELF = 0.55. The stereoactive lone pair electrons of Bi3+ are clearly seen. The black lines indicate the unit cell.

four oxygen atoms, forming a PdO4 unit. The PdO4 units are stacked along the c axis, sharing two oxygens with the BiO6 octahedra that are asymmetrically distorted due to the stereoactive lone pair 6s2 electrons of Bi3+. The largest O−Bi−O angles in the BiO6 octahedra are 137° (in the ab plane) and 167° (along the c-axis). The corresponding Bi−O bond lengths are 2.765 and 2.376 Å, respectively. Each distorted BiO6 2530

DOI: 10.1021/acs.chemmater.6b04230 Chem. Mater. 2017, 29, 2529−2534

Article

Chemistry of Materials

as indicated by the large charge density lobe in the [001] tunnel in Figure 1d. We also computed the Born effective charges (Z*). The Z* diagonal elements for the Bi3+ cation (3.78; 3.78; 4.76) and oxygen anion (−1.84; −2.82; −2.28) are significantly larger than their normal charge states +3 and −2, respectively. This anomalous Z* indicates a strong dynamic charge transfer associated with the change of Bi−O bond lengths and the deviation of quadratic relation between energy and atom displacement38,39 and, therefore, large anharmonicity, which results in low κL. On the other hand, the avoided crossing between the polar optical mode Eu and the acoustic modes is clearly seen along the Γ−X direction but absent in the Γ−Z direction. This avoided crossing is found in clathrate materials with rattling modes and is indicative of strong coupling between optical and acoustic modes.6,40 The large gap at the avoided crossing point indicates a high coupling strength, which significantly increases the phonon scattering rates, reduces acoustic mode velocities, and again contributes to the low κL.6,41 In Figure 5, we show the phonon mode scattering rates of three phonon interaction in the absorption (Γ+: λ + λ′ → λ″) and emission (Γ−: λ → λ′ + λ″) processes. The peaks in the scattering rates plot show the scattering magnitude of first phonon mode (λ) induced by the second phonon mode (λ′). In the absorption process, a low-frequency phonon mode combines with other low-frequency modes, giving rise to a high-frequency optical mode. In the emission process, the phonon mode is only allowed to decompose into a lowerfrequency mode, thus restricting the second phonon mode (λ′) in the right lower triangle. Both processes satisfy energy and crystal-momentum conservation. In Figure 5a, the scattering peaks centered around (ωλ ≈ 45 cm−1, ωλ′ ≈ 45 cm−1) indicate that acoustic modes, which carry the most heat, are dominantly scattered by the low-frequency (∼40 cm−1) optical modes near the avoided crossing, consistent with high coupling strength indicated by the avoided crossing. Phonon modes with ωλ > 100 cm−1 show considerable scattering in both absorption and emission processes as can be seen in Figure 5a,b. Therefore, the low κL is a combined effect of strong anharmonic phonon− phonon interaction and small group velocity, stemming from the stereochemically active lone pair electrons of Bi3+ 6s2 and the tunnel structure. In the square planar ligand field, d orbitals split into four energy levels, dxz/dyz, dz2, dxy, and dx2−y2 from low to high energy. However, the presence of weak interactions with ligands along the c-axis, i.e., the elongated octahedral ligand field, increases the energy level of the dz2 orbital. The amount of energy level shift directly depends on the interaction strength (or the distance between the center metal and its ligands along the c-

Figure 3. Lattice thermal conductivity of Bi2PdO4 as a function of temperature.

than most TE oxides including the layered oxides PbPdO2 (9 W m−1 K−1) and NaxCoO2−δ (19 W m−1 K−1).27,35,36 To understand the low κL, in Figure 4 we plot the third-order force

Figure 4. Norm of the calculated third-order IFCs.

constants Φmnl =

∂ 3E ∂um∂un∂ul

(E and u are the total energy and

atom displacement, respectively) for different atom species m, n, and l. Since phonon scattering rates are roughly proportional to |Φ|2,37 a high value of Φmnl suggests large anharmonicity. Both ΦBi3 (

∂ 3E ) ∂uBi ∂uBi ∂uBi

and ΦO3 (

∂ 3E ) ∂uO∂uO∂uO

play the most

important role, consistent with the large anharmonicity stemming from the stereochemically active Bi 6s2 lone pair,

Figure 5. Phonon scattering rates of the three phonon process in (a) absorption (Γ+: λ + λ′ → λ″) and (b) emission (Γ−: λ → λ′ + λ″) processes. 2531

DOI: 10.1021/acs.chemmater.6b04230 Chem. Mater. 2017, 29, 2529−2534

Article

Chemistry of Materials

gaps of small band gap compounds,46,47 the band gap of Bi2PdO4 should be between them. As we discussed earlier, the top of valence band is composed mainly of the dz2 orbital states of the Pd2+ cation. As overlap of dz2 with the px and py orbitals of the in-plane oxygen ligands is negligible, a relatively flat band in the ab plane is expected. In fact, as seen in Figure 6, the valence band maximum along Γ− X−M−Γ is extremely flat, affording a high density of states (DOS) near the Fermi level. On the other hand, a proper overlap between Pd2+ dz2 orbitals along the c-axis leads to very dispersive bands along Γ−Z, X−R, and M−A (c-axis) directions, indicating small band effective masses and therefore high carrier mobility. To clearly see this effect, we applied uniaxial tensile strain along the c-axis to artificially increase the Pd−Pd distance from 3.083 to 3.43 Å. Figure S2 shows that the dz2 bandwidth and dz2/dxy band order are very sensitive to the Pd−Pd distance. Therefore, a proper Pd−Pd (square planar stacking) distance in Bi2PdO4 is very important to achieve such flat-and-dispersive band structure. As demonstrated by Bilc et al.,22 high PF along the dispersive band direction is expected once the dispersive and flat bands touch together, i.e., X and M in the case of Bi2PdO4; see Figure 6. Another advantage of the Bi2PdO4 band structure for a high PF is the high band degeneracy Nv, which is usually observed when multiple bands have similar energies at the same valley point (orbital degeneracy) or when the valley is located at a low symmetry point of the Brillouin zone of a high symmetry lattice (valley degeneracy). As we can see in Figure 6, the orbital degeneracy of the valence band along the X−M direction is twofold, and the valley degeneracies at M, X, and Y (between M and X) are 1, 2, and 4, respectively. Therefore, the band degeneracy at Y (Nv = 8) is comparable with the high Nv compounds usually only observed in cubic lattices, such as the half Heusler FeMSb (M = V and Nb) compound.49 In general, good TE materials with ZT > 1 almost always have Nv ≥ 6.50

axis). In the case of Bi2PdO4, the Pd−Pd distance between two stacked PdO4 square planar is 3.083 Å. Therefore, Pd2+ is effectively in an elongated octahedral ligand field, and the dz2 level switches order with dxy, becoming the top of the valence band. Since the ligand field splitting in 4d transition metals is large, the low spin state is always preferred. Consequently, the four low-energy levels are fully occupied by the eight electrons of the Pd2+ cation, with dz2 as the highest occupied orbital, resulting in an insulating ground state, as depicted in Figure 6.

Figure 6. DFT (PBE) band structure (left) and density of state (DOS) of Bi2PdO4. The inset is the decomposed charge density of the top of the valence band.

Therefore, Bi2PdO4 is a band semiconductor, as for the divalent Pd2+ cation in PbPdO2.36 The calculated DFT band gaps for PBE42 and HSE0643 exchange-correlation functionals are 0.15 and 1.41 eV, respectively, consistent with the semiconducting transport behavior observed in previous experimental measurements.44,45 Since the semilocal functionals such as PBE generally underestimate the band gap and HSE06 (25% Hartree−Fock exchange43) trends to overestimate the band

Figure 7. Seebeck coefficient S and |S| (a) and S2σ/τ (b) of Bi2PdO4 along the c-axis and SrTiO3 as a function of temperature and carrier concentration. Power factor (σS2) (c) and ZT (d) of Bi2PdO4 as function of chemical potential, compared with cubic SrTiO3. The absolute value of S is used for SrTiO3 as it is an n-type semiconductor. Electron relation time (τ) is assumed as 4.3 fs for SrTiO3 at 300 K. The experiment lattice thermal conductivity (κL = 10 W m−1 K−1 at 300 K, ref 48) of SrTiO3 is used. 2532

DOI: 10.1021/acs.chemmater.6b04230 Chem. Mater. 2017, 29, 2529−2534

Article

Chemistry of Materials

Office of Science, Basic Energy Sciences, under Grant DEFG02-07ER46433. This research was supported in part through the computational resources and staff contributions provided for the Quest high performance computing facility at Northwestern University which is jointly supported by the Office of the Provost, the Office for Research, and Northwestern University Information Technology.

Electron transport properties of Bi2PdO4 were calculated using Boltzmann transport theory within relaxation time approximation as implemented in BoltzTraP code.51 The evolutions of S and S2σ/τ with respect to temperature (T) and carrier concentration (n) are shown in Figure 7a,b, respectively. Hole-doped Bi2PdO4 has a slightly higher Seebeck coefficient along the c-axis (Szz) than |S| (absolute value of S) of the n-type oxide SrTiO3 in wide ranges of T and n. SrTiO3 is well-known as a high-Seebeck n-type oxide.52 Combined with the higher value of σzz (smaller m* along Γ−Z direction), Bi2PdO4 shows two to four times higher S2zzσzz/τ than SrTiO3 at the same T and n as shown in Figure 7b. Figure 7c,d shows the PF and ZT as a function of relaxation time at T = 300 K. Bi2PdO4 has a much higher PF and ZT than SrTiO3 with the same relaxation time τ = 4.3 fs (which is adopted from ref 53), due to the higher σzz and lower κL. We note that previous experimental work on PbPdO2 has shown that p-type transport can be achieved by substituting Pd2+ with Li+ because Pd2+ has similar ionic radius with Li+.36 Since the Pd2+ cation in Bi2PdO4 has a coordination environment (stacked square planar) very similar to that in PbPdO2, we surmise that p-type Bi2PdO4 is achieved by substituting Pd2+ with Li+ as well.36 In this paper, a promising p-type TE oxide with high power factor and low lattice thermal conductivity is demonstrated by using first-principles lattice dynamic and electron Boltzmann transport theory within the relaxation time approximation. The yy zz low lattice thermal conductivity (κxx L = κL = 2.2 and κL = 2.8 W −1 −1 m K ) and high power factor of Bi2PdO4 not only suggest the possibility of using air-stable oxides as efficient TE materials but also offer new insights for discovering and designing highZT TE materials by combing stereochemically active lone pair electrons and d8 cations in a stacked square planar ligand field. The other strategies for reducing lattice thermal conductivity, such as nanostructuring, can be used to further enhance ZT in this system as well.





(1) Snyder, G. J.; Toberer, E. S. Complex thermoelectric materials. Nat. Mater. 2008, 7, 105−114. (2) Mahan, G.; Sofo, J. The best thermoelectric. Proc. Natl. Acad. Sci. U. S. A. 1996, 93, 7436−7439. (3) Sootsman, J. R.; Chung, D. Y.; Kanatzidis, M. G. New and old concepts in thermoelectric materials. Angew. Chem., Int. Ed. 2009, 48, 8616−8639. (4) Zhao, L.-D.; Lo, S.-H.; Zhang, Y.; Sun, H.; Tan, G.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals. Nature 2014, 508, 373−377. (5) Zhao, L.-D.; Tan, G.; Hao, S.; He, J.; Pei, Y.; Chi, H.; Wang, H.; Gong, S.; Xu, H.; Dravid, V. P.; Uher, C.; Snyder, G. J.; Wolverton, C.; Kanatzidis, M. G. Ultrahigh power factor and thermoelectric performance in hole-doped single-crystal SnSe. Science 2016, 351, 141−144. (6) Toberer, E. S.; Zevalkink, A.; Snyder, G. J. Phonon engineering through crystal chemistry. J. Mater. Chem. 2011, 21, 15843−15852. (7) Wu, H. J.; Zhao, L.-D.; Zheng, F. S.; Wu, D.; Pei, Y. L.; Tong, X.; Kanatzidis, M. G.; He, J. Q. Broad temperature plateau for thermoelectric figure of merit ZT > 2 in phase-separated PbTe0.7S0.3. Nat. Commun. 2014, 5, 4515. (8) Gelbstein, Y.; Dashevsky, Z.; Dariel, M. In-doped Pb 0.5 Sn 0.5 Te p-type samples prepared by powder metallurgical processing for thermoelectric applications. Phys. B (Amsterdam, Neth.) 2007, 396, 16−21. (9) Gelbstein, Y.; Davidow, J. Highly efficient functional GexPb1−xTe based thermoelectric alloys. Phys. Chem. Chem. Phys. 2014, 16, 20120− 20126. (10) Gelbstein, Y. Phase morphology effects on the thermoelectric properties of Pb0.25Sn0.25Ge0.5Te. Acta Mater. 2013, 61, 1499−1507. (11) Dado, B.; Gelbstein, Y.; Mogilansky, D.; Ezersky, V.; Dariel, M. P. Structural Evolution Following Spinodal Decomposition of the Pseudoternary Compound (Pb0.3Sn0.1Ge0.6)Te. J. Electron. Mater. 2010, 39, 2165−2171. (12) Xie, W.; Weidenkaff, A.; Tang, X.; Zhang, Q.; Poon, J.; Tritt, T. M. Recent advances in nanostructured thermoelectric half-Heusler compounds. Nanomaterials 2012, 2, 379−412. (13) Kirievsky, K.; Shlimovich, M.; Fuks, D.; Gelbstein, Y. An ab initio study of the thermoelectric enhancement potential in nanograined TiNiSn. Phys. Chem. Chem. Phys. 2014, 16, 20023−20029. (14) Kirievsky, K.; Gelbstein, Y.; Fuks, D. Phase separation and antisite defects in the thermoelectric TiNiSn half-Heusler alloys. J. Solid State Chem. 2013, 203, 247−254. (15) Uher, C. In Semiconductors and Semimetals; Tritt, T. M., Ed.; Recent Trends in Thermoelectric Materials Research I; Elsevier: 2001; Vol. 69, pp 139−253. (16) Nolas, G. S.; Slack, G. A.; Schujman, S. B. In Semiconductors and Semimetals; Tritt, T. M., Ed.; Recent Trends in Thermoelectric Materials Research I; Elsevier, 2001; Vol. 69, pp 255−300. (17) He, J.; Amsler, M.; Xia, Y.; Naghavi, S. S.; Hegde, V. I.; Hao, S.; Goedecker, S.; Ozoliņ s,̌ V.; Wolverton, C. Ultralow Thermal Conductivity in Full Heusler Semiconductors. Phys. Rev. Lett. 2016, 117, 046602. (18) Lin, H.; Tan, G.; Shen, J.-N.; Hao, S.; Wu, L.-M.; Calta, N.; Malliakas, C.; Wang, S.; Uher, C.; Wolverton, C.; Kanatzidis, M. G. Concerted Rattling in CsAg5Te3 Leading to Ultralow Thermal Conductivity and High Thermoelectric Performance. Angew. Chem., Int. Ed. 2016, 55, 11431−11436.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.chemmater.6b04230. Computational details of the electrical transport, phonon dispersion, lattice thermal conductivity, HSE06 band structure, Fermi surface, and band structures under uniaxial strain (PDF)



REFERENCES

AUTHOR INFORMATION

Corresponding Author

*(C.W.) E-mail: [email protected]. ORCID

Jiangang He: 0000-0001-9643-3617 Chris Wolverton: 0000-0003-2248-474X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS J.H. acknowledges support via ONR STTR N00014-13-P-1056. S.H. acknowledges support by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Award Number DE-SC0014520. Y.X., S.S.N., V.O., and C.W. acknowledge support by the U.S. Department of Energy, 2533

DOI: 10.1021/acs.chemmater.6b04230 Chem. Mater. 2017, 29, 2529−2534

Article

Chemistry of Materials (19) Lu, X.; Morelli, D. T.; Xia, Y.; Zhou, F.; Ozoliņs,̌ V.; Chi, H.; Zhou, X.; Uher, C. High Performance Thermoelectricity in EarthAbundant Compounds Based on Natural Mineral Tetrahedrites. Adv. Energy Mater. 2013, 3, 342−348. (20) Skoug, E. J.; Morelli, D. T. Role of Lone-Pair Electrons in Producing Minimum Thermal Conductivity in Nitrogen-Group Chalcogenide Compounds. Phys. Rev. Lett. 2011, 107, 235901. (21) Nielsen, M. D.; Ozolins, V.; Heremans, J. P. Lone pair electrons minimize lattice thermal conductivity. Energy Environ. Sci. 2013, 6, 570−578. (22) Bilc, D. I.; Hautier, G.; Waroquiers, D.; Rignanese, G.-M.; Ghosez, P. Low-Dimensional Transport and Large Thermoelectric Power Factors in Bulk Semiconductors by Band Engineering of Highly Directional Electronic States. Phys. Rev. Lett. 2015, 114, 136601. (23) Walia, S.; Balendhran, S.; Nili, H.; Zhuiykov, S.; Rosengarten, G.; Wang, Q. H.; Bhaskaran, M.; Sriram, S.; Strano, M. S.; Kalantarzadeh, K. Transition metal oxides−Thermoelectric properties. Prog. Mater. Sci. 2013, 58, 1443−1489. (24) Koumoto, K.; Wang, Y.; Zhang, R.; Kosuga, A.; Funahashi, R. Oxide thermoelectric materials: a nanostructuring approach. Annu. Rev. Mater. Res. 2010, 40, 363−394. (25) 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, 4414−4426. (26) Gaultois, M. W.; Sparks, T. D.; Borg, C. K.; Seshadri, R.; Bonificio, W. D.; Clarke, D. R. Data-driven review of thermoelectric materials: performance and resource considerations. Chem. Mater. 2013, 25, 2911−2920. (27) Fujita, K.; Mochida, T.; Nakamura, K. High-temperature thermoelectric properties of NaxCoO2-δ single crystals. Jpn. J. Appl. Phys. 2001, 40, 4644. (28) Koumoto, K.; Terasaki, I.; Funahashi, R. Complex oxide materials for potential thermoelectric applications. MRS Bull. 2006, 31, 206−210. (29) Bilc, D. I.; Floare, C. G.; Zârbo, L. P.; Garabagiu, S.; Lemal, S.; Ghosez, P. First-principles Modelling of SrTiO3 based Oxides for Thermoelectric Applications. J. Phys. Chem. C 2016, 120, 25678− 25688. (30) Snyder, G. J.; Toberer, E. S. Complex thermoelectric materials. Nat. Mater. 2008, 7, 105−114. (31) Zeier, W. G.; Zevalkink, A.; Gibbs, Z. M.; Hautier, G.; Kanatzidis, M. G.; Snyder, G. J. Thinking Like a Chemist: Intuition in Thermoelectric Materials. Angew. Chem., Int. Ed. 2016, 55, 6826− 6841. (32) Asen-Palmer, M.; Bartkowski, K.; Gmelin, E.; Cardona, M.; Zhernov, A. P.; Inyushkin, A. V.; Taldenkov, A.; Ozhogin, V. I.; Itoh, K. M.; Haller, E. E. Thermal conductivity of germanium crystals with different isotopic compositions. Phys. Rev. B: Condens. Matter Mater. Phys. 1997, 56, 9431−9447. (33) Zhang, Y.; Skoug, E.; Cain, J.; Ozoliņs,̌ V.; Morelli, D.; Wolverton, C. First-principles description of anomalously low lattice thermal conductivity in thermoelectric Cu-Sb-Se ternary semiconductors. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 054306. (34) Zhou, F.; Nielson, W.; Xia, Y.; Ozoliņ s ,̌ V. Lattice Anharmonicity and Thermal Conductivity from Compressive Sensing of First-Principles Calculations. Phys. Rev. Lett. 2014, 113, 185501. (35) El-Sharkawy, A. A.; El-Azm, A. M. A.; Kenawy, M. I.; Hillal, A. S.; Abu-Basha, H. M. Thermophysical properties of polycrystalline PbS, PbSe, and PbTe in the temperature range 300−700 K. Int. J. Thermophys. 1983, 4, 261−269. (36) Lamontagne, L. K.; Laurita, G.; Gaultois, M. W.; Knight, M.; Ghadbeigi, L.; Sparks, T. D.; Gruner, M. E.; Pentcheva, R.; Brown, C. M.; Seshadri, R. High Thermopower with Metallic Conductivity in pType Li-Substituted PbPdO2. Chem. Mater. 2016, 28, 3367.

(37) 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. (38) Zhao, X.; Vanderbilt, D. Phonons and lattice dielectric properties of zirconia. Phys. Rev. B: Condens. Matter Mater. Phys. 2002, 65, 075105. (39) Waghmare, U. V.; Spaldin, N. A.; Kandpal, H. C.; Seshadri, R. First-principles indicators of metallicity and cation off-centricity in the IV-VI rocksalt chalcogenides of divalent Ge, Sn, and Pb. Phys. Rev. B: Condens. Matter Mater. Phys. 2003, 67, 125111. (40) Christensen, M.; Abrahamsen, A. B.; Christensen, N. B.; Juranyi, F.; Andersen, N. H.; Lefmann, K.; Andreasson, J.; Bahl, C. R.; Iversen, B. B. Avoided crossing of rattler modes in thermoelectric materials. Nat. Mater. 2008, 7, 811−815. (41) Li, W.; Carrete, J.; Madsen, G. K. H.; Mingo, N. Influence of the optical-acoustic phonon hybridization on phonon scattering and thermal conductivity. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 205203. (42) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (43) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2003, 118, 8207− 8215. (44) Bettahar, N.; Conflant, P.; Abraham, F. Effect of copper substitution on the electrical transport properties of (Bi,Pb)2MO4 (M = Pd, Pt) linear chain compounds. J. Alloys Compd. 1992, 188, 211− 214. (45) Chater, E.; Sellami, M.; Bettahar, N. Impact of palladium substitution for copper in Bi2CuO4 synthesized via sol-gel method. Eur. Phys. J.: Appl. Phys. 2011, 53, 11502. (46) Marques, M. A. L.; Vidal, J.; Oliveira, M. J. T.; Reining, L.; Botti, S. Density-based mixing parameter for hybrid functionals. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 83, 035119. (47) He, J.; Franchini, C. Screened hybrid functional applied to 3d0→ 3d8 transition-metal perovskites LaMO3 (M = Sc̆Cu): Influence of the exchange mixing parameter on the structural, electronic, and magnetic properties. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86, 235117. (48) Popuri, S. R.; Scott, A.; Downie, R.; Hall, M.; Suard, E.; Decourt, R.; Pollet, M.; Bos, J.-W. Glass-like thermal conductivity in SrTiO 3 thermoelectrics induced by A-site vacancies. RSC Adv. 2014, 4, 33720−33723. (49) Fu, C.; Zhu, T.; Pei, Y.; Xie, H.; Wang, H.; Snyder, G. J.; Liu, Y.; Liu, Y.; Zhao, X. High Band Degeneracy Contributes to High Thermoelectric Performance in p-Type Half-Heusler Compounds. Adv. Energy Mater. 2014, 4, 1400600. (50) Tang, Y.; Gibbs, Z. M.; Agapito, L. A.; Li, G.; Kim, H.-S.; Nardelli, M. B.; Curtarolo, S.; Snyder, G. J. Convergence of multivalley bands as the electronic origin of high thermoelectric performance in CoSb3 skutterudites. Nat. Mater. 2015, 14, 1223. (51) Madsen, G. K.; Singh, D. J. BoltzTraP. A code for calculating band-structure dependent quantities. Comput. Phys. Commun. 2006, 175, 67. (52) Mehdizadeh Dehkordi, A.; Bhattacharya, S.; Darroudi, T.; Graff, J. W.; Schwingenschlögl, U.; Alshareef, H. N.; Tritt, T. M. Large Thermoelectric Power Factor in Pr-Doped SrTiO3−δ Ceramics via Grain-Boundary-Induced Mobility Enhancement. Chem. Mater. 2014, 26, 2478−2485. (53) García-Fernández, P.; Verissimo-Alves, M.; Bilc, D. I.; Ghosez, P.; Junquera, J. First-principles modeling of the thermoelectric properties of SrTiO3/SrRuO3 superlattices. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86, 085305.

2534

DOI: 10.1021/acs.chemmater.6b04230 Chem. Mater. 2017, 29, 2529−2534