Combination of Carrier Concentration Regulation and High Band

Degeneracy for Enhanced Thermoelectric Performance of. Cu3SbSe4. Dan Zhang, Junyou Yang*, Qinghui Jiang*, Zhiwei Zhou, Xin Li, Jiwu Xin, Abdul Basit, ...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Combination of Carrier Concentration Regulation and High Band Degeneracy for Enhanced Thermoelectric Performance of Cu3SbSe4 Dan Zhang, Junyou Yang, Qinghui Jiang, Zhiwei Zhou, Xin Li, Jiwu Xin, Abdul Basit, Yangyang Ren, Xu He, Weijing Chu, and Jingdi Hou ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.7b08121 • Publication Date (Web): 09 Aug 2017 Downloaded from http://pubs.acs.org on August 10, 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 22

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

Combination of Carrier Concentration Regulation and High Band Degeneracy for Enhanced Thermoelectric Performance of Cu3SbSe4 Dan Zhang, Junyou Yang*, Qinghui Jiang*, Zhiwei Zhou, Xin Li, Jiwu Xin, Abdul Basit, Yangyang Ren, Xu He, Weijing Chu, Jingdi Hou

State Key Laboratory of Material Processing and Die & Mould Technology, Huazhong University of Science & Technology, Wuhan 430074, PR China.

ABSTRACT The effect of Al-, Ga- and In-doping on the thermoelectric (TE) properties of Cu3SbSe4 has been comparatively studied on the basis of theoretical prediction and experimental validation. It is found that tiny Al/Ga/In substitution leads to a great enhancement of electrical conductivity with high carrier concentration and also large Seebeck coefficient due to the preserved high band degeneracy and thereby a remarkably high power factor. Ultimately, coupled with the depressed lattice thermal conductivity, all three elements (Al/Ga/In) substituted samples have obtained a highly improved thermoelectric performance with respect to un-doped Cu3SbSe4. In compared with the samples at the same Al/In doping level, the slight Ga-doped sample presents better TE performance over the wide temperature range, and the Cu3Sb0.995Ga0.005Se4 sample presents a record high ZT value of 0.9 among single doped Cu3SbSe4 at 623K, which is about 80% higher than that of pristine Cu3SbSe4. This work offers an alternative approach to boot the TE properties of Cu3SbSe4 by selecting efficient dopant to weaken the coupling between electrical conductivity and Seebeck coefficient.

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 2 of 22

KEYWORDS: Cu3SbSe4; Al/Ga/In doping; carrier concentration; high band degeneracy; decoupling; thermoelectric properties

1. INTRODUCTION The R&D for alternative energy technologies have been continually promoted by the increasing global energy consumptions and resulted environmental issues.1,

2

Thermoelectric (TE) materials have drawn extensive attention as a key solution in recent years because they can realize the direct conversion between heat and electricity. The conversion efficiency of a TE device is decided by the material figure-of-merit ZT = S2σT/κ, where S, σ, κ and T are the Seebeck coefficient, electrical conductivity, thermal conductivity (Comprising the lattice contribution κL and the electronic part κe) and the absolute temperature, respectively. Most of the state-of-the-art TE materials that have high ZT value more or less include some toxic, rare, or expensive elements which potentially

hamper

their

commercialization

process.3-9

Therefore,

Cu-based

semiconductors have attracted great interests recently as a cheap and Pb-free robust TE candidate.10-12 Owing to advantages of large carrier effective mass (~1.1me) and small band-gap (~0.29eV),13 stannite semiconductor Cu3SbSe4 is stepping into sight as a cheap and Pb-free Cu-based promising TE material. Many efforts have been made to substitute Ge or Sn for Sb to improve its intrinsically poor electrical conductivity,13-17 however, less improvement has been obtained in TE performance mainly due to the simultaneously great diminution of Seebeck coefficient when doping with Ge or Sn, and the separation

ACS Paragon Plus Environment

Page 3 of 22

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

of the degenerate valence bands is the main cause for the rapid decrease of Seebeck coefficient in the heavy Ge/Sn-doped Cu3SbSe4 (Figure S1).18,

19

Therefore, it is a

challenge to optimize the interrelated electrical conductivity and Seebeck coefficient simultaneously for high performance Cu3SbSe4 TE material. Herein, we systematically explored the doping effect of M (M= Al, Ga, In) on TE properties of Cu3SbSe4 by means of theoretical prediction and experimental validation. After tiny efficient Al/Ga/In substitution, the electrical resistivity decreases obviously while the Seebeck coefficient stays high because almost no destruction to the high degeneracy of valence band,20 and the lattice thermal conductivity also decreases by the introduced point defects scattering to phonons. As a result, a greatly enhanced ZT value has been obtained over the wide temperature range for all the Cu3Sb1-xMxSe4 (M=Al, Ga, and In) samples, and the maximum ZT value approaches 0.9 at 623K in Cu3Sb0.995Ga0.005Se4, as well as we know, it is the highest value ever reported for single doped Cu3SbSe4-based TE materials at 623K.

2 EXPERIMENTAL SECTION 2.1 Synthesis. Cu3Sb1-xMxSe4 (x= 0, 0.005, 0.01, 0.015, 0.02; M= Al, Ga, In) samples were synthesized via high temperature melting, ball milling and hot-pressing method. Stoichiometric mixtures of elemental powders with high purity were firstly weighted and sealed in evacuated silica tubes, and heated at 1173K for 12 h, then cooled to 773 K in the furnace and quenched. The quenched ingots were directly annealed at 573K for 2 days to make them more homogeneous. After that, they were pulverized, sieved and ball

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

milled. In the end, the obtained powders were hot-pressed at 673 K for 1h under 120 MPa. 2.2 Characterization and Measurement. The phase structure was characterized by X-ray diffraction (XRD) on a step scan mode with Cu Kα radiation (PANalytical, Almelo, Netherlands). The microstructure of the samples was analyzed by field-emission scanning electron microscopy (FESEM, NanoSEM 450 ) equipped with energy-dispersive X-ray spectroscopy (EDS). Transmission electron microscopy (TEM) observation was performed on a JEM-2100 microscope operated with an acceleration voltage of 200 kV. The simultaneous measurement of Seebeck coefficient and electrical resistivity was carried out on a commercial measurement system (Namicro-III, Wuhan, China) from room temperature to 623K. The thermal conductivity (κ) was calculated as the product of material density determined by the Archimedes’ method, thermal diffusivity measured by a laser flash method (LFA 427, Netzsch Instruments, Germany) and specific heat derived by differential scanning calorimeter (Diamond DSC, PerkinElmer Instruments, USA). Carrier concentration and mobility of the sample were measured on a Hall effect measurement system (HMS 5500, Ecopia, South Korea) under constant magnetic fields of ± 0.5 T. 2.3 Theory Calculation. The first-principles calculations have been carried out on the Vienna ab initio simulation package (VASP) using the Perdew-Burke-Ernzerhof (PBE) 21

generalized gradient approximation (GGA) and projector augmented wave (PAW) 22

pseudopotentials. Hubbard-like term (+U) was also introduced in order to deal with the localization of d-electrons on Cu atoms.23-25 A supercell model of 2×2×1 tetragonal unit

ACS Paragon Plus Environment

Page 4 of 22

Page 5 of 22

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

cell was selected within our computation resources, a plane-wave cutoff energy of 400 eV, 6×6×6 Monkhorst-Pack k-point mesh and an energy conversion threshold of 1.0×10–5eV per atom were used for calculations.

3 RESULTS AND DISCUSSION 3.1 Electronic Structure Calculation and Analysis

Figure 1. Band structure of pristine (a) and doped Cu3SbSe4: (b) AlSb; (c) GaSb; (d) InSb.

As Al/Ga/In possesses one valence electron less than Ge/Sn, its substitution for Sb probably provide more p-type carriers in comparison to that of Ge/Sn. Hence, the electronic structure of Al-, Ga- and In-doped Cu3SbSe4 were calculated and the results were shown in Figure 1-2. As shown in Figure 1a, the undoped Cu3SbSe4 presents three-fold degenerate bands at Gamma point, which results in the intrinsically large Seebeck coefficient.24 It is shown that Al, Ga and In tend to occupy Sb site in Cu3SbSe4 due to the smallest calculated forming energy than those of other sites (Cu & Se).25, 26

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

The Fermi level moves down-ward into valence band on doping with Al/Ga/In (Figure1b-d), thus resulting in higher hole concentration and the enhanced electrical conductivity in Cu3SbSe4. However, the above degenerate valence bands diverge in heavy doped Cu3SbSe4 which will certainly lead to the diminution of Seebeck coefficient, and the similar split of degenerate bands has also been reported in other heavy doped Cu3SbSe4.17 Thus, selecting highly efficient dopants to greatly lower the substitution concentration should be a feasible attempt to improve the electrical conductivity while maintain the large Seebeck coefficient with little divergence of degenerate bands. Moreover, our theoretical calculation results (Figure 1-2 and Figure S1) also indicate that Al, Ga and In dopants are more efficient in contributing carriers in Cu3SbSe4 with respect to the widely used Ge/Sn elements, because they shift the Fermi level deeper into the valence band even at the same doping level. Therefore, Al/Ga/In substitution is implemented to enhance the electrical properties and thereby improve the TE properties of Cu3SbSe4 in the following sections.

Figure 2. The change of DOS after doping with Al/Ga/In/Ge/Sn.

ACS Paragon Plus Environment

Page 6 of 22

Page 7 of 22

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

3.2 Effect of Ga-doping content on the Electrical and Thermal Transport Properties of Cu3SbSe4

Figure 3. Powder XRD patterns of Cu3Sb1-xGaxSe4 samples (a), a representative Rietveld refinement result using X-ray diffraction data from matrix (b), the results of lattice parameters (c), and room temperature carrier mobility µ & carrier concentration nh as a function of the Ga-doping content (d).

Considering the smallest difference in electronegativity and moderate ionic radius difference between Ga (χ(Ga) = 1.81, r(Ga3+) = 0.47Ǻ) and Sb (χ(Sb) = 2.05, r(Sb5+)= 0.6Ǻ) in compared with Al (χ(Al) = 1.68, r(Al3+)= 0.39Ǻ) and In (χ(In) = 1.78, r(In3+)= 0.62Ǻ), Ga probably has higher solid solubility in Cu3SbSe4,27 thus Ga-doping is firstly studied. Figure 3a shows the powder XRD pattern of Ga-doped Cu3SbSe4 samples. One can observe that all the diffraction peaks match very well with Cu3SbSe4 (No. 01-085-0003) and no visible impurity phase appears within detection limit, revealing all the Ga-substituted samples are single phase Cu3SbSe4. The calculated lattice parameters from Rietveld refinement on Fullprof software28 (Figure 3c) decrease when Ga enters

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

into Sb site because the smaller ionic radius of Ga than Sb and the resulted shrinkage of unit cell. Noticeably, they level off when x(Ga) exceeds 0.015, indicating that the solution limit of Ga in Cu3SbSe4 is around 0.015. Moreover, Figure 3d shows that Ga can contribute approximately 1.5 holes per atom at low doping content, and the carrier concentration (nh) saturates when x exceeds the solution limit. At the same time, the carrier mobility (µ) decreases after Ga-doping due to the extra ionized impurity scattering.

Figure 4. FESEM images of x(Ga) = 0.01 (a) and the corresponding elemental EDS mapping (b-e), the enlarged image (f) of selected region in (a), and the TEM image (g) with an inserted high-resolution transmission electron microscopy (HRTEM) picture.

ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22

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

Figure 5. Temperature dependent thermoelectric properties for the Cu3Sb1-xGaxSe4 (x = 0, 0.005, 0.01, 0.015, 0.02) bulk samples: (a) electrical resistivity, (b) Seebeck coefficient, (c) power factor, (d-e) total and lattice thermal conductivities, and (f) ZT value.

The microstructure of the x(Ga) = 0.01 sample is analyzed by FESEM and TEM (Figure 4). The multi-scaled grains in bulk sample are beneficial to suppress lattice thermal conductivity by integrating phonon scattering in wide range (Figure 4a).5, 20 Elements of Cu, Sb, Ga and Se are found to have a homogeneous distribution in the EDS mapping images (Figure 4b-e), which implies Ga has doped into Cu3SbSe4 uniformly without

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

obvious second phase aggregate. The nano-scaled particles can be clearly observed in the enlarged SEM and TEM images (Figure 4f-g), which were confirmed as Cu3SbSe4 by HRTEM lattice image. Figure 5 displays the TE properties of Cu3Sb1-xGaxSe4 (x = 0, 0.005, 0.01, 0.015, 0.02) bulk samples. As shown in Figure 5a-b, both the electrical resistivity (ρ) and Seebeck coefficient (S) increase with the rising of temperature after Ga doping, indicating the Ga-doped samples have changed from the intrinsic semiconductor of pristine Cu3SbSe4 into a degenerate semiconductor and this consists well with the results of first-principles calculations (Figure 1c). The ρ and S are decreased as the fraction of Ga increases mainly due to the increased carrier concentration (Figure 3d), and the reduction in ρ for the doped samples is very remarkable as compared with un-doped Cu3SbSe4 due to the efficient p-type doping of Ga, for example, it decreases from 195.88µΩ.m (x = 0) to 58.54µΩ.m (x = 0.005) only with tiny amount of Ga substitution. As a consequence, the power factor has greatly improved in the wide temperature region (Figure 5c) for the Ga-doped samples and the maximum one is up to 1312.73µWm-1K-2, however, its enhancement starts to decline when x(Ga) is above 0.005 mainly due to the diminution of S from the split of the degenerate bands at higher doping concentration (Figure 1c). As described in Figure 5d, the thermal conductivity (κ) decreases slightly after Ga-doping, and the lattice thermal conductivity (κL) obtained by subtracting the carrier contribution κe from κ was also presented in Figure 5e, where κe is estimated using the Wiedemann-Franz law as the way in our previous work,12, 29 one can observe that κL decreases after doping with Ga due to the increased point-defects scattering of phonons,

ACS Paragon Plus Environment

Page 10 of 22

Page 11 of 22

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

and it agrees with the relationship of κL ∝ 1/T, which suggests that the dominant scattering in phonon transport is the acoustic phonon-phonon Umklapp scattering. Finally, a peak ZT value of 0.90 at 623K is achieved for the x(Ga) = 0.15 sample, which is about 80% higher than that of pristine Cu3SbSe4. 3.3 Distinct Impact of Al, Ga and In Substitution on the Thermoelectric Properties of Cu3SbSe4

Figure 6. Powder XRD spectrums of Cu3Sb1-xMxSe4 (M= Al, Ga and In) samples (a) and the results of lattice parameters from Rietveld refinement (b).

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

Figure 7. Temperature dependent thermoelectric properties for the Cu3Sb1-xMxSe4 (x = 0.005, 0.01;M = Al, Ga, In) bulk samples: (a) electrical resistivity, (b) Seebeck coefficient, (c) power factor (PF), (d) change of Seebeck coefficient with Hall carrier concentration at room temperature (i. e. Pisarenko curve), (e) total and lattice thermal conductivities, and (f) ZT value. The total thermal conductivities are displayed in the inserted figure of 7(e).

As the solid solubility of Ga is very low and high Ga doping density can quickly diminish the Seebeck coefficient in Cu3SbSe4, the substitution concentration of other same group elements (Al and In) is designed below 0.01. Figure 6a displays the powder XRD pattern of the bulk Al-, Ga- and In-doped Cu3SbSe4 samples. All appeared

ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22

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

diffraction peaks exhibit good match with the Cu3SbSe4 (No. 01-085-0003) and no visible impurity phase can be observed within detection limit, indicating all the bulk Al/Ga/In-doped Cu3SbSe4 samples are single phase. The refined lattice parameters are also shown in Figure 6b. The smallest change of lattice parameters in Al-doped samples shows the effective Al doped into the matrix lattice is lowest, thus, Al has lowest solid solubility in Cu3SbSe4 mainly due to the largest ionic radius and electronegativity difference between Sb; Element In should also has lower solid solubility than that of Ga as the increase in In-doping concentration from 0.005 to 0.01 makes tiny change in lattice parameters, which can be ascribed to the relatively larger electronegativity difference between In-Sb as comparison to that of Ga-Sb. Therefore, their solid solution limit in Cu3SbSe4 are in sequence of Ga > In > Al.

Figure 8. The ZT values at 623K for the present Cu3Sb0.995A0.005Se4 (A= Al, Ga, In) and other single doped Cu3SbSe4.

The TE properties of Al- Ga- and In-doped Cu3SbSe4 are compared and presented in Figure 7. One can see that all doped samples have become degenerate semiconductor

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

(Figure 7a-b) and agrees well with our theoretical calculations (Figure 1b-d). The values of ρ and S in doped samples are in order of Al > In > Ga, which consists with their corresponding solid solubility in Cu3SbSe4 and the resultant carrier concentration. The greatly decreased ρ and relatively large S for all doped samples enable the remarkable enhancement in power factor (Figure 7c). The S versus nh curve in Figure 7d clearly gives us an insight into the advantage of the highly efficient Al/Ga/In substitution compared with the widely used Ge/Sn doping: 0.5 percent of Al/Ga/In substitution can provide comparable or even higher carrier concentration than that of one percent Ge-doping and also has little deviation from the Pisarenko curve due to the nearly non-changed band structure.30,

31

The preserved high degeneracy for the x = 0.005

samples, which can be indirectly reflected from the Pisarenko curve, is beneficial to large Seebeck coefficient. At the same time, the high doping concentration, such as x(Ga)= 0.015/0.02; x(Ge)= 0.01/0.02 and x(Sn)= 0.025/0.05, make the S deteriorate rapidly with obvious deviation from the Pisarenko curve due to severe split of the degenerate bands at VBM (Figure 1 and S1).32 In addition, the calculated κL decreases after doping with Al/Ga/In because of the increased point-defects scattering to phonons and thus κ also has a slight decline. As a consequence, the ZT values of Al/Ga/In substituted samples obtain great improvement in wide temperature range with respect to un-doped matrix, and the maximum one increases from 0.83 for Cu3Sb0.995Al0.005Se4 to 0.88 for Cu3Sb0.995In0.005Se4 and then to 0.90 for Cu3Sb0.995Ga0.005Se4 sample. As compared in Figure 8, to the best of our knowledge, the peak ZT value of 0.9 is the highest value for single doped Cu3SbSe4 at 623K in literature13,14,16,17,33 and even higher

ACS Paragon Plus Environment

Page 14 of 22

Page 15 of 22

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

than the maximum ZT value at 650K of Ge-S co-doped Cu3SbSe4.15 Moreover, considering the high thermal conductivity in Al/Ga/In doped Cu3SbSe4, there is still plenty of room to further boot the TE performance of this promising Pb- & Te-free TE material by phonon engineering, such as nano-structuring,34 introduction of nano-sized second phase.35-37

CONCLUSION In summary, the detailed effect of Al-, Ga- and In- doping on the TE properties of Cu3SbSe4 has been comparatively studied by combining theory and experiment. We demonstrate in this intrinsically high band degeneracy Cu3SbSe4 compound that tiny Al/Ga/In substitution (x= 0.005) leads to an enhanced electrical conductivity with high carrier concentration, a large Seebeck coefficient with preserved high band degeneracy and thereby a remarkably high power factor.

As compared the samples at the same

Al/In doping level, slight Ga-doped sample has better TE performance in wide temperature range mainly due to the high electrical conductivity from the high solid solubility in Cu3SbSe4, and the achieved maximum ZT values of 0.9 at 623K is about 80% higher than that of pristine Cu3SbSe4, and is also the reported highest ZT value for single doped Cu3SbSe4 at the same temperature in literature.

ASSOCIATED CONTENT Supporting Information Band structure of the doped Cu3SbSe4: (a) GeSb, (b) SnSb. This material is available free of charge via the Internet at http://pubs.acs.org.

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

AUTHOR INFORMATION Corresponding Author *Email: [email protected] *Email: [email protected] Notes The authors declare no competing financial interest.

ACKNOWLEDGMENTS This work is co-financed by National Natural Science Foundation of China (Grant NO. 51572098 and 51632006), National Basic Research Program of China (Grant NO. 2013CB632500), Natural Science Foundation of Hubei province (Grant NO. 2015CFB432), Open Fund of State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology (NO. 2016-KF-5). The technical assistance from the analytical and testing center of HUST is likewise gratefully acknowledged.

REFERENCES [1] Bell, L. E. Cooling, Heating, Generating Power, and Recovering Waste Heat with Thermoelectric Systems. Science 2008, 321, 1457-1461. [2] DiSalvo, F. J. Thermoelectric Cooling and Power Generation. Science 1999, 285, 703-706. [3] Tan, G.; Shi, F.; Hao, S.; Zhao, L.; Uher, C.; Wolverton, C.; Dravid, V.; Kanatzidis, M. Non-equilibrium Processing Leads to Record High Thermoelectric Figure of Merit in PbTe-SrTe. Nat. Commun. 2016, 7, 12167.

ACS Paragon Plus Environment

Page 16 of 22

Page 17 of 22

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

[4] Zhao, L. D.; Lo, S. H.; Zhang, Y.; Sun, H.; Tan, G.; Uher, C.; Wolverton, C.; Dravid, V.; Kanatzidis, M. Ultralow Thermal Conductivity and High Thermoelectric Figure of Merit in SnSe Crystals. Nature 2014, 508, 373-377. [5] Rogl, G; Grytsiv, A.; Rogl, P.; Peranio, N.; Bauer, E.; Zehetbauer, M.; Eibl, O. N-Type Skutterudites (R,Ba,Yb)yCo4Sb12 (R= Sr, La, Mm, DD, SrMm, SrDD) Approaching ZT≈ 2.0. Acta Mater. 2014, 63, 30-43. [6] Li, W.; Zheng, L.; Ge, B.; Lin, S.; Zhang, X.; Chen, Z.; Chang, Y.; Pei, Y. Promoting SnTe as An Eco-Friendly Solution for p-PbTe Thermoelectric via Band Convergence and Interstitial Defects. Adv. Mater. 2017, 29, 1605887. [7] Luo, Y.; Yang, J.; Li, G.; Liu, M.; Xiao, Y.; Fu, L.; Li, W.; Zhu, P.; Peng, J.; Gao, S.; Zhang, J. Enhancement of the Thermoelectric Performance of Polycrystalline In4Se2.5 by Copper Intercalation and Bromine Substitution. Adv. Energy Mater. 2014, 4, 1300599. [8] Kim, S. I.; Lee, K. H.; Mun, H. A.; Kim, H. S.; Hwang, S. W.; Roh, J. W.; Yang, D. J.; Shin, W. H.; Li, X. S.; Lee, Y. H.; Snyder, G. J.; Kim, S. W. Dense Dislocation Arrays Embedded in Grain Boundaries for High-Performance Bulk Thermoelectric, Science 2015, 348, 109-114, [9] Luo, Y.; Yang, J.; Liu, M.; Xiao, Y.; Fu, L.; Li, W.; Zhang, D.; Zhang, M.; Cheng, Y. Multiple Heteroatom Induced Carrier Engineering and Hierarchical Nanostructures for High Thermoelectric Performance of Polycrystalline In4Se2.5. J. Mater. Chem. A 2015, 3, 1251-1257. [10] Qiu, P.; Shi, X.; Chen, L. Cu-based Thermoelectric Materials. Energy Storage

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

Mater. 2016, 3, 85-97. [11] Zhang, R.; Chen, K.; Du, B.; Reece, M. Screening for Cu-S Based Thermoelectric Materials Using Crystal Structure Features. J. Mater. Chem. A 2017, 5, 5013-5019. [12] Zhang, D.; Yang, J.; Jiang, Q.; Fu, L.; Xiao, Y.; Luo, Y.; Zhou, Z. Ternary CuSbSe2 Chalcostibite: Facile Synthesis, Electronic-Structure and Thermoelectric Performance Enhancement. J. Mater. Chem. A 2016, 4, 4188-4193. [13] Wei, T.; Wang, H.; Gibbs, Z. M.; Wu, C.; Snyder, G. J.; Li, J. Thermoelectric Properties of Sn-Doped P-Type Cu3SbSe4: A Compound with Large Effective Mass and Small Band Gap. J. Mater. Chem. A 2014, 2, 13527-13533. [14] Skoug, E.; Cain, J.; Majsztrik, P.; Kirkham, M.; Lara-Curzio, E.; Morelli, D. Doping Effects on the Thermoelectric Properties of Cu3SbSe4. Sci. Adv. Mater. 2011, 3, 602-606. [15] Skoug, E.; Cain, J.; Morelli, D. High Thermoelectric Figure of Merit in the Cu3SbSe4- Cu3SbS4 Solid Solution. Appl. Phys. Lett. 2011, 98, 261911. [16] Yang, C.; Huang, F.; Wu, L.; Xu, K. New Stannite-Like P-Type Thermoelectric Material Cu3SbSe4. J. Phys. D-Appl. Phys. 2011, 44, 295404. [17] Liu, Y.; García, G.; Ortega, S.; Cadavid, D.; Palacios, P.; Lu, J.; Ibáñez, M.; Xi, L.; Roo, J.; López, A.; Martí-Sánchez, S.; Cabezas, I.; Mata, M.; Luo, Z.; Dun, C.; Dobrozhan, O.; Carroll, D.; Zhang, W.; Martins, J.; Kovalenko, M.; Arbiol, J.; Noriega, G.; Song, J.; Wahnón, P.; Cabot, A. Solution-Based Synthesis and Processing of Sn-and Bi-Doped Cu3SbSe4 Nanocrystals, Nanomaterials and Ring-Shaped Thermoelectric Generators. J. Mater. Chem. A 2017, 5, 2592-2602.

ACS Paragon Plus Environment

Page 18 of 22

Page 19 of 22

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

[18] Pei, Y.; Wang, H.; Snyder, G. Band Engineering of Thermoelectric Materials. Adv. Mater. 2012, 24, 6125-6135. [19] 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, 3848-3853. [20] Jin, H.; Wiendlocha, B.; Heremans, J. P-type Doping of Elemental Bismuth with Indium, Gallium and Tin: A Novel Doping Mechanism in Solids. Energy Environ. Sci. 2015, 8, 2027-2040. [21] Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979. [22] Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. [23] Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558-561. [24] Do, D. T.; Mahanti, S. D. Bonds, Bands, and Band Gaps in Tetrahedrally Bonded Ternary Compounds: The Role of Group V Lone Pairs. J. Phys. Chem. Solids 2014, 75, 477-485. [25] Do, D. T.; Mahanti, S. D. Theoretical Study of Defects Cu3SbSe4: Search for Optimum Dopants for Enhancing Thermoelectric Properties. J. Alloys Compd. 2015, 625, 346-354. [26] Shen, J.; Zhang, X.; Chen, Z.; Lin, S.; Li, J.; Li, W.; Li, S.; Chen, Y.; Pei, Y. Substitutional Defects Enhancing Thermoelectric CuGaTe2. J. Mater. Chem. A 2017, 5,

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 20 of 22

5314-5320. [27] Shannon, R. D. Revised Effective Ionic Radii and Systematic Studies of Interatomic Distances in Halides and Chalcogenides. Acta Crystallogr. 1976, 32, 751-767. [28] Roisnel, T.; Rodríquez-Carvajal, J. WinPLOTR: A Windows Tool for Powder Diffraction Pattern Analysis. Mater. Sci. Forum 2001, 378, 118-123. [29] Zhang, D.; Yang, J.; Jiang, Q.; Zhou, Z.; Li, X.; Basit, A.; Ren, Y.; He, X. Multi-cations Compound Cu2CoSnS4: DFT Calculating, Band Engineering and Thermoelectric Performance Regulation. Nano Energy 2017, 36, 156-165. [30] Toberer, E. S.; May, A. F.; Melot, B. C.; Flage-Larsen, E.; Snyder, G. J. Electronic Structure and Transport in Thermoelectric Compounds AZn2Sb2 (A= Sr, Ca, Yb, Eu). Dalton Trans. 2010, 39, 1046-1054. [31] Fu, C.; Zhu, T.; Liu, Y.; Xie, H.; Zhao, X. Band Engineering of High Performance P-Type FeNbSb Based Half-Heusler Thermoelectric Materials for Figure of Merit zT> 1. Energy Environ. Sci. 2015, 8, 216-220. [32] 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, 1600007. [33] Wei T.; Li F.; Li J.; Enhanced Thermoelectric Performance of Nonstoichiometric Compounds Cu3- xSbSe4 by Cu Deficiencies. J. Electron. Mater. 2014, 43, 2229-2238. [34] Fu, L.; Yang, J.; Peng, J.; Jiang, Q.; Xiao, Y.; Luo, Y.; Zhang, D.; Zhou, Z.; Zhang,

ACS Paragon Plus Environment

Page 21 of 22

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

M.; Cheng, Y.; Cheng, F. Enhancement of Thermoelectric Properties of Yb-Filled Skutterudites by An Ni-Induced “Core-Shell”Structure. J. Mater. Chem. A 2015, 3, 1010-1016. [35] Ge, Z.; Qin, P.; He, D.; Chong, X.; Feng, D.; Ji, Y.; Feng, J.; He, J. Highly Enhanced Thermoelectric Properties of Bi/Bi2S3 Nanocomposites. ACS Appl. Mater. Interfaces. 2017, 9, 4828-4834. [36] Zhao, W.; Liu, Z.; Wei, P.; Zhang, Q.; Zhu, W.; Su, X.; Tang, X.; Yang, J.; Liu, Y.; Shi, J.; Chao, Y.; Lin, S.; Pei, Y. Magnetoelectric Interaction and Transport Behaviours in Magnetic Nanocomposite Thermoelectric Materials. Nat. Nanotechnol. 2016, 12, 55-60. [37] Yang, J.; Aizawa, T.; Yamamoto, A.; Ohta T. Effects of Interface Layer on Thermoelectric Properties of a Pn Junction Prepared via the BMA-HP Method. Mater. Sci. Eng. B, 2001, 85, 34-37.

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

Table of Contents Graphic 33x13mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 22 of 22