Electronic Band Structure Engineering and Enhanced Thermoelectric

Jan 8, 2018 - Heating and cooling cycles were carried out to assess the thermal stability during the measurement and gave repeatable data. ..... This ...
1 downloads 17 Views 2MB Size
Subscriber access provided by READING UNIV

Article

Electronic band structure engineering and enhanced thermoelectric transport properties in Pb-doped BiCuOS oxysulfide Jean-Baptiste Labégorre, Rabih Al Rahal Al Orabi, Agathe Virfeu, Jacinthe Gamon, Philippe Barboux, Lauriane Pautrot-d'Alençon, Thierry Le Mercier, David Berthebaud, Antoine Maignan, and Emmanuel Guilmeau Chem. Mater., Just Accepted Manuscript • DOI: 10.1021/acs.chemmater.7b04989 • Publication Date (Web): 08 Jan 2018 Downloaded from http://pubs.acs.org on January 8, 2018

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.

Chemistry of Materials 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 33 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

Chemistry of Materials

Electronic band structure engineering and enhanced thermoelectric transport properties in Pb-doped BiCuOS oxysulfide

Jean-Baptiste Labégorre,1 Rabih Al Rahal Al Orabi,2,3* Agathe Virfeu,2 Jacinthe Gamon,4,5 Philippe Barboux,5 Lauriane Pautrot-d’Alençon,4 Thierry Le Mercier,4 David Berthebaud,1 Antoine Maignan,1 Emmanuel Guilmeau1,*

1

Laboratoire CRISMAT, UMR-CNRS 6508, ENSICAEN, UNICAEN, Normandie

Université, 6 Boulevard du Maréchal Juin, 14050 Caen Cedex 04, France 2

Solvay, Design and Development of Functional Materials Department, Axel’One, 87 avenue

des Frères Perret, 69192 Saint Fons, Cedex, France 3

Department of Physics, Central Michigan University, Mt. Pleasant, Michigan 48879, United

States 4

Solvay, Centre de Recherches d’Aubervilliers, 52 rue de la Haie-Coq, 93308 Aubervilliers

Cedex, France 5

Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie de

Paris, , 11 rue Pierre et Marie Curie, 75005 Paris, France *Email : [email protected], Tel : +33 231 45 13 67 *Email : [email protected], Tel : +33 428 27 10 80

ACS Paragon Plus Environment

1

Chemistry of Materials 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 33

Abstract In this paper, Bi1-xPbxCuOS samples (0 ≤ x ≤ 0.05) have been synthesized with a simple and scalable ball-milling process, followed by a reactive Spark Plasma Sintering. Our results highlight that, Pb for Bi substitution increases the charge carriers concentration by more than two orders of magnitude from 1.4 × 1017 cm-3 to 2.6 × 1019 cm-3 for x = 0 and x = 0.05 respectively. As a result, the electrical resistivity is divided by more than 50 at room temperature and the Seebeck coefficient drops from 707 µV K-1 to 265 µV K-1 where our experimental results are supported with Density Functional Theory (DFT) calculations. Electronic structure calculations show that, just below the top of the valence band, several other bands are present and may contribute to the transport properties with appropriate tuning of the heavy-light valence band and the position of the Fermi level. Pb doping increases the number of holes pockets and several band degeneracies appear around the Fermi level, leading to a drastic enhancement of the power factor up to 0.2 mW m-1 K-2 at 700 K. This is 5 times higher than the value of the pristine compound. The intrinsically low thermal conductivity of 0.7 W m-1 K-1 at 700 K is interpreted on the basis of vibrational properties calculations within the Density Functional Perturbation Theory (DFPT) approach. It indicates that soft acoustic modes along Γ-Z direction suggest weak interatomic bonding between the layers and possible strong anharmonicity. The power factor being enhanced with a minimal impact on the thermal conductivity, the figure of merit ZT reaches 0.2 at 700 K for x = 0.05. To the best of our knowledge, it is considered the best reported value among the family of oxysulfides

ACS Paragon Plus Environment

2

Page 3 of 33 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

Chemistry of Materials

Introduction Thermoelectric devices are promising candidates to reduce the CO2 emissions and fight against the energy crisis due to their ability to achieve direct and reversible energy conversion between heat and electricity. The conversion efficiency of a thermoelectric material is commonly described by the dimensionless figure of merit ZT (ZT = S2T/ρκ), where T is the absolute temperature (K), S is the Seebeck coefficient (µV K-1), ρ is the electrical resistivity (Ω m) and κ is the thermal conductivity (W m-1 K-1).1 However, current thermoelectric modules are confined to niche markets due to elements scarcity, high costs, toxicity and low thermal stability under air. During the last two decades, efforts have been devoted to develop thermoelectric oxides that are capable of addressing the aforementioned issues.2–5 A significant step has been especially crossed with the discovery of relatively low thermal conductivity and high electrical performances (p-type), that exists inside the natural superlattices of NaxCoO26 and Ca3Co4O9.7,8 Nevertheless, these oxides are still not considered to have good performance compared to other thermoelectric materials because of their relatively high electrical resistivity and large thermal conductivity. In 2010, another layered material BiCuOSe attracted tremendous attention due to its promising ZT in the medium temperature range (0.32 at 673 K).9 This oxychalcogenide superlattice crystallizes in the tetragonal P4/nmm space group and is built of two different layers that are stacked along the c axis alternately (Figure 1).10 The covalent anti-fluorite-like layer (Cu2Se2)2-, build up from edge-sharing CuSe4 tetrahedra, is the conductive pathway for charge carriers transport. The ionic (Bi2O2)2+ block is insulating and acts as a charge reservoir. The succession of these two weakly bonded layers behave like phonon scattering interfaces and give it a very low lattice thermal conductivity (~ 0.5 W m-1 K-1 at 773 K).11 On the other hand, the high electrical resistivity problem should be solved in order to enhance the merit figure of this material further. This has been achieved by increasing the carrier concentration with monovalent or

ACS Paragon Plus Environment

3

Chemistry of Materials 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 4 of 33

divalent cations (Na+, K+, Mg2+, Ca2+, Sr2+, Ba2+ or Pb2+)9,12–18 substitution inside the insulating (Bi2O2)2+ layer. Furthermore, other strategies have been used, such as the creation of copper vacancies in the conductive (Cu2Se2)2- layer,

19

as well as limiting the decrease of

the charge carriers mobility via modulation doping.20 Electronic band structure engineering through the isoelectronic substitution of Se by Te has also been developed successfully.21 These studies led to a substantial improvement of the thermoelectric performances to record ZT values of 1.4 (at 923 K) and 1.5 (at 873 K) for doped and textured samples22 and for dually doped samples respectively.23 Thus, this material became a viable alternative to the use of PbTe and nearly competes with the high performance of multiple-filled skutterudites (ZT = 1.7 at 850 K).24 Despite the achievements reached in the BiCuOSe system in the past years, there is a lack of studies on the isostructural sulfur counterpart due to its lower thermoelectric performances, i.e. a ZT ∼0.07 at 650 K.25 The low value of ZT is mainly attributed to (1) a lower carrier concentration (1017 cm-3 against ∼1019 cm-3 in BiCuOSe), (2) a wider band gap (1.1 eV

against 0.75 eV in BiCuOSe) which is much higher than the optimal value predicted by the Goldsmid-Sharp formula Eg = 2eSmaxTmax of 0.35 eV for thermoelectric materials,26 and (3) a higher thermal conductivity in BiCuOS than in BiCuOSe (S is much lighter than Se). Nevertheless, sulfur being significantly less toxic and more abundant than selenium, there would be undeniable benefits in replacing Se by S, while maintaining relatively thermoelectric performances. Pb doping has been well investigated in p-type BiCuOSe for its role in the significant improvement of the power factor compared to those of pristine BiCuOSe. However, to the best of our knowledge, electrical transport properties were solely investigated for Pb-doped BiCuOS.27–29 The similarities of the band and crystal structure between BiCuOS and BiCuOSe gave a hint that Pb-doping could modify the electronic properties and enhance the ACS Paragon Plus Environment

4

Page 5 of 33 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

Chemistry of Materials

power factor. In this paper, we present the effect of Pb doping on the band structure and electronic properties of BiCuOS, using both theoretical and experimental studies for the first time. We prove that Pb doping in BiCuOS improves the electronic properties by activating several hole-pockets with high effective masses, resulting in an enhanced power factor. This yields to a high figure of merit ZT of 0.2 at 700 K for a Pb concentration of 5%. To the best of our knowledge, it is the best ever achieved figure of merit value for an oxysulfide compound.

Fig. 1 Crystal structure of 2×2×2 supercells of BiCuOS (left) and Bi0.937Pb0.063CuOS (right) used for the computations.

Experimental details Bi1-xPbxCuOS (x = 0, 0.025, 0.0375, 0.05) powders were synthesized from stoichiometric mixtures of Bi2O3 (> 99.5 %, Alfa Aesar), Bi (> 99.99 %, Alfa-Aesar), PbO (> 99 %, Merck), Cu (> 99 %, < 625 mesh, Alfa Aesar) and S (> 99 %, < 325 mesh, Alfa Aesar). The precursors were weighed and 5.5 g of each composition were set inside a 45 mL tungsten carbide jar with 7 tungsten carbide balls (Ø 10 mm). The jar was closed in argon-filled glove box to ensure an inert atmosphere during the process. The mechanosynthesis was held in a planetary

ACS Paragon Plus Environment

5

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

Page 6 of 33

ball mill (Fritsch Pulverisette 7 Premium line) during 6 cycles of 30 min at 600 rpm. Finally, the resulting powder were ground, placed in 10 mm graphite dies and densified by Spark Plasma Sintering (SPS) under a uniaxial pressure of 64 MPa under vacuum at 748 K for 25 min (heating and cooling rates of 100 K/min). The relative densities of the samples were respectively 90.6 % (x = 0), 94.5 % (x =0.025), 94.9 % (x = 0.0375) and 94.1 % (x = 0.05) of the theoretical values. Powder X-Ray diffraction (PXRD) patterns were collected at room temperature using a PANalytical X-Pert Pro diffractometer equipped with a Cu Kα X-ray tube. The scanning angle 2θ was varied between 20° and 100°. PXRD patterns were refined by the Rietveld method with FullProf and WinPLOTR software packages.30,31 Details of the refinements are listed in Table S1. Microstructures and grain sizes of sintered samples were examined by scanning electron microscopy (SEM) on a ZEISS Supra 55. The electrical resistivity (ρ) and the Seebeck coefficient (S) were simultaneously measured on bar-shaped samples with typical size of 3.0 × 3.0 × 8.0 mm3 using the standard four-probe method (ULVAC-RIKO ZEM-3). The measurements were performed from room temperature to 673 K under partial helium pressure. Heating and cooling cycles were carried out to assess the thermal stability during the measurement and gave repeatable data. The thermal conductivity (κ) was calculated from the product of the geometrical density, heat capacity (Dulong-Petit) and thermal diffusivity obtained by laser flash diffusivity method on a Netzsch LFA 457 (with 6.0 × 6.0 × 1.0 mm3 samples) measured between 300 K and 673 K. Finally, Hall effect measurements by Van der Pauw method were carried out in a PPMS (Physical Properties Measurement System) system from Quantum Design to determine the carrier concentrations (nH = 1/eRH where e is the electronic charge) and mobilities at room temperature. Samples with dimensions of 3.0 × 3.0

× 0.3 mm3 were used with indium as contacts. The Hall coefficients were calculated from the linear fits of the Hall resistivities versus magnetic fields between -7 T and 7 T. It is

ACS Paragon Plus Environment

6

Page 7 of 33 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

Chemistry of Materials

noteworthy that the electrical and thermal properties were measured along the direction perpendicular to the compression applied during the sintering process. The estimated measurements uncertainties are 6% for the Seebeck coefficient, 8% for the electrical resistivity, 11% for the thermal conductivity and 16% for the final figure of merit, ZT.32

First-Principle Calculations Density Functional Theory (DFT) geometry optimizations were carried out with the CASTEP8.0 code33 using a set of ultrasoft pseudopotentials34,35 with the PBEsol exchange– correlation functional.36 Cell parameters and atomic positions were both relaxed. During the geometry optimizations, a convergence threshold of 0.02 eV Å−1 was used for the residual forces and 0.1 kbar for the pressure. The cut-off energy for plane-waves was set to 700 eV. For the electronic band structures we used the full-potential linearized augmented plane wave (FLAPW) approach, as implemented in the WIEN2K code.37 A plane-wave cutoff corresponding to RMTKmax= 7 was used in all calculations. The radial wave functions inside the non-overlapping muffin-tin spheres were expanded up to lmax=12. The charge density was Fourier expanded up to Gmax=16 Å-1. Total energy convergence was achieved with respect to the Brillouin zone (BZ) integration mesh with 500k-points. Because GGA exchangecorrelation functionals are known to underestimate experimental band gaps, we used the modified Becke-Johnson (mBJ) functional which leads to excellent agreement with the experimental values for the energy separation between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO).38 Since the Bi atoms in oxychalcogenide compounds have fairly high atomic number and mass, the relativistic effect can not be neglected, and we have included the effect of spin-orbit coupling to elucidate realistic electronic structure. We used 5000 k-points in the BZ to compute the band derivatives and the density of states. The doping effect was examined by constructing a ACS Paragon Plus Environment

7

Chemistry of Materials 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 8 of 33

2×2×2 (64 atoms) supercell in which one of the Bi atoms was replaced by a Pb atom and the corresponding Pb doping concentration is 6.3 at.%. The band structure of both undoped and Bi0.937Pb0.063CuOS was calculated using the 2×2×2 supercell under the same calculation conditions. Vibrational properties have been obtained by using Density Functional Perturbation Theory (DFPT)39 as implemented in the Quantum ESPRESSO (QE) package.40 The calculation has been performed by using a set of ultrasoft pseudopotentials35 with the PBEsol exchangecorrelation functional.36 The kinetic energy cutoff and the charge density cutoff are set to be 40 Ry and 480 Ry, respectively. An 8 x 8 x 4 k-point mesh is used to sample the Brillouin zone. The Grüneisen’ parameter can be defined as  = −

    

and

characterizes the relationship between phonon frequency and volume change. In the quasiharmonic DFT phonon calculations, the system volume is isotropically expanded by +4% from the DFT relaxed volume.

Results and discussion Crystal structure Figure 2a displays the PXRD patterns of Bi1-xPbxCuOS (0 ≤ x ≤ 0.05) sintered samples. All the major Bragg peaks can be indexed in the ZrSiCuAs structure with the tetragonal P4/nmm space group.41 The samples with x ≤ 0.025 are single phase, whereas the main peak (012) of metallic bismuth (Bi) is discernible at ~ 27.1° for the two highest lead contents, indicating < 1 wt.% (x = 0.0375) and ~ 1 wt.% (x = 0.05) of bismuth metal in the samples, respectively. The increase of the bismuth content as secondary phase probably originates from the decrease of Bi vacancies formation energy when Pb fraction rises, as discussed in Bi1-xPbxCuOTe samples

ACS Paragon Plus Environment

8

Page 9 of 33 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

Chemistry of Materials

by Tae-Ho An et al.42 It is noteworthy that minor amounts of unreacted Bi and Bi2O3 are found after the ball-milling process, as highlighted in Figure S1. During SPS, the assynthesized powders react and high purity bulk polycristalline samples are formed when x ≤ 0.025. This short process enables a reduction of the synthesis time as compared to sealed tube (> 10 h)27 or hydrothermal (> 55 h)43 synthesis.

Fig. 2 (a) PXRD patterns of the Bi1-xPbxCuOS (0 ≤ x ≤ 0.05) samples processed by SPS. The green stars indicate the main diffraction peak (012) of metallic bismuth (R-3m) observed in

ACS Paragon Plus Environment

9

Chemistry of Materials 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 10 of 33

the samples x = 0.0375 and x = 0.05. Rietveld refinements of PXRD patterns of (b) x = 0 and (c) x = 0.025 sintered samples.

The refined unit cell parameters of BiCuOS, a = 3.870(1) Å and c = 8.561(1) Å are in good agreement with the values reported for single crystals, i.e. a = 3.8705(5) Å and c = 8.561(1) Å.10 Moreover, the Rietveld refinements of PXRD patterns (See for Figures 2b and 2c, for x = 0 and x = 0.025 respectively; refinement parameters are reported in Table S1) evidenced a copper site occupancy ranging between 0.95 and 0.99 for all the samples. This is considered consistent with the ability of MCuOCh (M = La, Bi ; Ch = S, Se, Te) system to easily form copper vacancies, leading to non-intentional p-type doping in pristine samples.27,29,44,45 As shown in Figure 3, the cell parameters a and c, which are determined from Rietveld refinements, increase monotonically with the Pb fraction from x = 0 to 0.0375. The linear increase of both lattice parameters follows Vegard’s law and is consistent with the substitution of lead on bismuth site because Pb2+ (129 pm) has a larger ionic radius than Bi3+ (117 pm).46 However, c axis expands faster than a axis when lead content is raised. This behavior reveals the weakness of Coulombic attraction between the layers due to substitution. Indeed, the holes transferred to the conductive layer decrease the difference of electronic charge between the two blocks [(Bi1-xPbx)2O2](2-2x)+ and [Cu2S2](2-2x)-. Thus the Coulombic attraction decreases with x which leads to a larger increase of c axis when the lead amount rises.47 The decrease of the lattice constants between x = 0.0375 and x = 0.05 could be due to Bi vacancies as evidenced by the increase of the Bi metal secondary phase. Such decrease of lattice parameters is due to Bi deficiencies which is already noted for BiCuOSe samples.48 It seems that the intrinsic nature of BiCuOS to form Cu and Bi vacancies hinders the synthesis of highly pure samples when the Pb fraction is high. This low doping level is in sharp contrast

ACS Paragon Plus Environment

10

Page 11 of 33 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

Chemistry of Materials

with the selenide-derivative, firstly reported with a Pb2+ solubility lower than 7 at.% before being synthesized until 20 at.%.18,49 Finally, a typical SEM micrograph of the fractured cross-section of the densified Bi1xPbxCuOS

sample is displayed in Figure 4. It reveals a fine and dense microstructure with an

average grain size between 100-200 nm. Such small grain size is mostly explained by the combination of mechanosynthesis and SPS processes. Indeed, as showed in Figure S1, the wide width of the diffraction peaks suggest that the crystallites are in a nanometric scale after the mechanosynthesis process. The subsequent grain growth during SPS is also limited due to the short time process (25 min) and the low sintering temperature to prevent from sulfur volatilization.

Fig. 3 Lattice parameters of the Bi1-xPbxCuOS (0 ≤ x ≤ 0.05) series after sintering, obtained by Rietveld refinements of the PXRD patterns.

ACS Paragon Plus Environment

11

Chemistry of Materials 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 12 of 33

Fig. 4 Typical SEM micrograph of Bi0.95Pb0.05CuOS sintered sample. The microstructure is representative of all the samples.

Thermoelectric properties The temperature dependences of the electrical resistivity (ρ) in the Bi1-xPbxCuOS series are illustrated in Figure 5a. The electrical resistivity of all the samples decreases upon heating, indicating a semiconducting behavior. Pristine BiCuOS exhibits high resistivity over the whole temperature range, from 5.3×104 mΩ cm at 300 K to 5.7×102 mΩ cm at 700 K. The electrical resistivity is decreased gradually by the Pb for Bi substitution, and reached minimum values of 1.0×103 mΩ cm at 300 K and 5.6×101 mΩ cm at 700 K for Bi0.95Pb0.05CuOS sample. Aiming to prove the impact of lead doping over the carrier concentration, Hall effect measurements were performed. The Hall coefficients obtained by the Van der Pauw method are all positive at 300 K, demonstrating that holes are the dominant carriers in the samples. The room temperature Hall carrier concentration (nH) of BiCuOS is approximately 1.4×1017 cm-3 (Figure 6). As discussed earlier in the structural part, the copper vacancies are considered to be the source of hole carriers in the undoped sample.27,29,44,45 The carrier concentration increases by more than two orders of magnitude and reaches 2.6×1019 cm-3 for x = 0.05. The simultaneous increase of the lattice parameters and of the charge ACS Paragon Plus Environment

12

Page 13 of 33 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

Chemistry of Materials

carriers concentration demonstrates that Pb2+ was successfully incorporated into the BiCuOS lattice. Thus, assuming that Pb substitution adopts the one carrier model, the decrease of the electrical resistivity originates from the holes generated according to the equation: 

×     + ℎ

The charge carriers mobility (µH) depicted in Figure 6 unveil a very low value, close to 1 cm2 V-1 s-1 in the undoped sample. This one is decreased further to 0.3 cm2 V-1 s-1 for x = 0.05. The drop of the mobility can be attributed to the increase of the ionized impurity scattering, to the rise of the charge carriers concentration and hole effective mass, as discussed below.50 By comparison, the mobility of the undoped selenide analogue is close to ~20 cm2 V-1 s-1.47 This significant contrast originates from a lower hybridization of Cu 3d and S 3p orbitals due to more localized 3p orbitals compared to Se 4p. Furthermore, the small grain size produced by the mechanical alloying process increases the grain boundaries density and the number of charge barriers. This is confirmed by mobilities as low as 2.8 cm2 V-1 s-1 in ball-milled BiCuOSe samples with 1-2 µm grains.11

ACS Paragon Plus Environment

13

Chemistry of Materials 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 14 of 33

Fig. 5 Temperature dependences of (a) electrical resistivity (b) Seebeck coefficient (c) power factor (d) thermal conductivity; blue dashed line is the calculated lattice thermal conductivity (for x = 0) and (e) ZT in the Bi1-xPbxCuOS (0 ≤ x ≤ 0.05) series.

ACS Paragon Plus Environment

14

Page 15 of 33 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

Chemistry of Materials

Fig. 6 Temperature dependences of the charge carriers concentration (nH) and mobility (µH) in the Bi1-xPbxCuOS (0 ≤ x ≤ 0.05) samples at 300 K.

As shown in Figure 5b, the temperature dependence of the Seebeck coefficient is positive for all samples over all the temperature range and corroborates p-type transport behavior. The pristine sample shows a large Seebeck coefficient close to 707 µV K-1 at 300 K. This value decreases to 545 µV K-1 when the temperature is raised to 700 K as expected for insulating or semiconductor materials. The Seebeck values of all the Pb substituted samples adopt a degenerated semiconductor behavior and increases with the temperature. The value reaches 339 µV K-1 at 700 K for the highest lead content. The magnitude of the Seebeck coefficient decreases with the increasing of Pb fraction and drops to 265 µV K-1 at 300 K for x = 0.05, which is consistent with the increase in carrier concentration. It must be noted that the Seebeck coefficient and the carrier concentration tend to saturate above x = 0.0375, when the secondary phase (Bi) appears. Nonetheless, the carrier concentration being governed by different structural features, such as cationic vacancies (Bi, Cu), Pb/Bi substitution, sulfur vacancies, it makes difficult to establish clearly the role of Bi metal on the TE properties.

ACS Paragon Plus Environment

15

Chemistry of Materials 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 16 of 33

Using our experimental Seebeck coefficient at 300 K, the chemical potential ( ! ) can be

estimated using equation 1 where kB is the Boltzmann constant, " is the acoustic phonon scattering (" = 0) and Fj(!) are the Fermi integrals given in equation 2. The hole effective mass can then be determined from equation 3 using the measured carrier concentration (Np) : $=

%& ) + *+, * - ( − -. 1 ' , + *+* - 5

+0 - = 1

6

Ϛ0 3Ϛ 2 , + 'Ϛ4- )D C

=> :) 8 = < B )%& ; ?@+,⁄) - ∗

3

As shown in Table 1, the large Seebeck coefficient is related to the significant effective mass of BiCuOS (1.71 me), which is higher than the value reported in BiCuOSe sample (1.1 me).13 Furthermore m* progressively increases from 1.71 me to 2.1 me (me being the free electronic mass) with increasing Pb content from 0 to 0.05, which is probably resulting from the increasing contribution of several hole pockets in the multi-valley valence band of BiCuOS.

Table 1 Calculated effective mass from experimental Seebeck coefficient values for the different compounds of Bi1-xPbxCuOS (0 ≤ x ≤ 0.05) series. Compounds m*/me

x=0 1.71

x = 0.025 1.81

x = 0.0375 1.85

x = 0.05 2.2

As shown in Figure 5c, the maximum power factor (PF = S2/ρ) is obtained for the Bi0.95Pb0.05CuOS sample at 700 K (0.2 mW m-1 K-2) and is 5 times higher than the undoped sample value (4.0 × 10-2 mW m-1 K-2). However, the limited carrier concentration and the low ACS Paragon Plus Environment

16

Page 17 of 33 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

Chemistry of Materials

electronic mobility restrain this enhancement and better values could probably be reached with higher effective doping. Aiming to clarify the improved electronic properties, we performed DFT calculations to determine the electronic structure of BiCuOS and 6.3 at.% Pb-doped BiCuOS (Figure 1). As mentioned in the introduction, the band structure of pristine BiCuOS exhibits features that do not favor high electronic properties for thermoelectric applications. This limit is related to the large band gap of the pristine BiCuOS (1.1 eV), which is much higher than the optimal values predicted by Goldsmid-Sharp formula Eg = 2eSmaxTmax of 0.35 eV for thermoelectric materials.26 As shown in Figure 7, the band structure calculations performed for pristine and 6.3 at.% doped BiCuOS confirm the indirect band-gaps of the compounds with valence band maximum (VBM) on the Γ-M line and conduction band minimum (CBM) located at the Γ point of the Brillouin zone. The computed band gap of BiCuOS compound is determined to be 1.08 eV which is in very good agreement with experimental measurements, estimated between 1.07 and 1.1 eV.27,51 Replacing Bi by Pb decreases the overall ionic character of the system with the effect of reducing the band gap to 0.55 eV for the 6.3 at.% doped sample. This value is high enough in order to prevent the bipolar conduction and maintain a high Seebeck coefficient. Our calculations show that, just below the top of the valence band, several other bands are present and may contribute to the transport properties with appropriate tuning of the heavy-light valence band and the position of the Fermi level. Since the substitution of Bi by Pb decreases the number of the electrons, the Fermi level moves downward as the content of Pb increases. As shown in the Figure 7b, Pb doping increases the number of holes pockets and several band degeneracies appear around the Fermi level. In order to clarify this concept, we calculated the Fermi surface for 6.3 at.% Pb doped BiCuOS (Figure 7c). It can be viewed that the number of hole pockets is enlarged, and the increase of

ACS Paragon Plus Environment

17

Chemistry of Materials 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 18 of 33

Nv implies an enhancement of the power factor. Also, similar band modifications were found in Sn-/Pb-doped BiCuOSe and Ga-doped SnTe.52–54

Fig. 7 Electronic band structure of BiCuOS (a) and Bi0.937Pb0.063CuOS (b). In Bi0.937Pb0.063CuOS, the hole pockets formed on R-A line and A-Z line are represented (c). The projected density of states (PDOS) of BiCuOS and Bi0.937Pb0.063CuOS are represented in Figure 8. From the computed density of states, it can be inferred that the top of the valence band mainly results from the hybridization between Cu 3d and S 3p orbitals while the minimum of the conduction band is dominated by Bi 6p states. Compared to BiCuOSe, the hybridization is higher between Cu 3d and Se 4p and explains the decrease of carrier mobility in the sulfur derivative. It is clear that, Pb-derived states seem to contribute substantially to the energy region near the top of the valence band (Figure 8b). The Pb contribution is dominated by the 6s components into the states close to the valence band (Figure S2). Such delocalized 6s orbitals from the lone pair electrons in Pb are considered as an important

ACS Paragon Plus Environment

18

Page 19 of 33 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

Chemistry of Materials

criteria to enhance the electronic transport as proved in the variation of the carrier concentration and mobility with Pb concentration. Similar results have been found in PbBiCuOSe doped sample.52

Fig. 8 Projected density of states of BiCuOS (a) and Bi0.937Pb0.063CuOS (b). Figure 5d illustrates the temperature dependences of the thermal conductivity in the Bi1xPbxCuOS

(0 ≤ x ≤ 0.05) series. The thermal conductivity of all the samples is close to 1.1 W

m-1 K-1 at room temperature and falls down to 0.67–0.71 W m-1 K-1 at 700 K. The intrinsically low thermal conductivity is characteristic of the BiCuOCh (Ch = S, Se, Te) family and represents the main advantage of these compounds for thermoelectric applications. This particularity originates from the structure itself which involves weak chemical bonding between the oxide and chalcogenide layers (i.e. low Young’s modulus), high anharmonic interactions of phonons (i.e. large Grüneisen parameter), a low vibration frequency provided by Cu atoms and strong phonons scattering at layers interfaces.55 Furthermore, point defects induced by copper/bismuth vacancies may play a role in the reduction of the lattice thermal conductivity. The total thermal conductivity (κ) stems from the contributions of the lattice thermal conductivity (κL) and the electronic thermal conductivity (κe) according the relation κ = κL + κe. The electrical thermal conductivity is proportional to the electrical conductivity

ACS Paragon Plus Environment

19

Chemistry of Materials 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 33

according to the Wiedmann-Franz relation κe = LTσ. Here, the Lorenz number L, was estimated as function of temperature from the experimental Seebeck coefficients using single parabolic band model (L and κe are respectively displayed in Figures S3 and S4).56 The calculated Lorenz numbers for all samples are in the range of 1.45 × 10-8 W Ω K-2 and 1.55

× 10-8 W Ω K-2, which are lower than the metallic limit of 2.45 × 10-8 W Ω K-2.57

The lattice thermal conductivity curves show no obvious trend related to the lead percent doping, indicating a weak effect of the induced point defects over the phonon scattering. This can be explained by the small mass fluctuation between Pb (207.20 g mol-1) and Bi (208.98 g mol-1) and the moderate atomic size mismatch of Pb2+ (129 pm) and Bi3+ (117 pm). Due to the low carrier concentration (2.6×1019 cm-3), the electronic heat transport contribution is limited below 2.7 % of the total thermal conductivity at 700 K even for the highest lead concentration. Thus the lattice thermal conductivity is the dominant contribution to the total thermal conductivity. In order to understand the low thermal conductivity of BiCuOS, we calculated the Grüneisen parameters and phonons dispersion of BiCuOS compound. As displayed in Figure 9a, the three acoustic modes have low boundary frequency around 65 cm-1 along Γ-X (010 direction) and Γ-M (110 direction), and 50 cm-1 along Γ-Z (001 direction). It was shown that low acoustic modes frequencies are beneficial for thermoelectric materials and lead to very low thermal conductivity.58,59 All these low frequencies acoustics modes indicate small group velocities around center (Γ-X, Γ-M and Γ-Z) and low Debye temperatures (Table 2), which lead to low lattice thermal conductivities. It can be seen that the acoustic modes are soft along Γ-Z direction are suggesting weak interatomic bonding between the layers and possible strong anharmonicity. It is known that the Grüneisen parameters are useful to estimate the lattice thermal anharmonicity and to interpret the physical nature of the lattice thermal conductivity.

ACS Paragon Plus Environment

20

Page 21 of 33 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

Chemistry of Materials

Thus, we calculated the Grüneisen parameters (γ) of the acoustic modes in function of the Brillouin zone (Figure 9b), and found a high value γ ~ 8 along Γ-Z direction. This indicates a very strong anharmonicity between the layer. Similar results have been found in BiCuOSe.60

Fig. 9 (a) Phonon dispersion of BiCuOS compound. The black, red and blue lines highlight transverse (TA/TA’) and longitudinal (LA) acoustic modes, respectively. The green lines highlight the optical modes. (b) Acoustic phonon dispersion of transverse (TA/TA’) and longitudinal (LA) modes (straight lines), associated Grüneisen parameters (dots) in function of Brillouin zone of BiCuOS compound.

Using the phonon dispersion, we can calculate the longitudinal (LA) and transverse (TA/TA’) velocities (the slope of the acoustic phonon dispersion around the Γ point) and their Debye temperature (FG ) calculated using Debye model given by : FG = H I/K  L6N O P

QD R

where H are the sound velocity associated to the longitudinal and transverse modes, and V is ACS Paragon Plus Environment

21

Chemistry of Materials 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 22 of 33

the average volume per atom. The averages longitudinal (LA) and transverse (TA/TA’) velocities, Debye temperature and Grüneisen parameter (γ) of each acoustic dispersion are presented in Table 2. The calculated Grüneisen parameters are in very good agreement with previous reported values in the literature,61,62 and close to those average Grüneisen parameters of SnSe.58 It is noteworthy that recently, Vaqueiro et al. found that, the low frequencies phonon modes and large Grüneisen parameters have largest contribution from the Cu atoms.55 On the basis of the Debye-Callaway model,63 we calculated the lattice thermal conductivity (Figure 5d). The calculated values of 1.1–0.65 W m-1 K-1 in the temperature range of 300– 700 K are in very good agreement with the experimental values, ranging between 1.1 and 0.7 W m-1 K-1. Table 2 Average transverse (TA/TA’) and longitudinal (LA) Grüneisen parameters (γ), Debye temperatures (θ) and phonon velocities of BiCuOS (ν).

BiCuOS

γTA

γTA’

γLA

θTA (K)

θTA’ (K)

θLA (K)

3.627

3.1674

4.2

95

95

206

νTA (m/s) 1534

νTA’ (m/s) 1534

νLA (m/s) 3490

Finally, we computed the theoretical minimum thermal conductivity κmin using the shortest scattering distance within the model proposed by Cahill et al.64 and obtained a value of ∼0.23 W m-1 K-1. For all samples, the κL values remain significantly higher than the κmin value, indicating that further reduction may be possible, and might for example be achieved by enhancing structural disordering. The combination of electrical and thermal transport properties allows to estimate the figure of merit (ZT) as a function of temperature in the Bi1-xPbxCuOS series (Figure 5e). The pristine sample exhibits a value close to 0.04 at 700 K, which is lower than the 0.07 (at 650 K)

ACS Paragon Plus Environment

22

Page 23 of 33 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

Chemistry of Materials

recently reported by Zhu et al.25 However, the higher carrier concentration (4.0 × 1018 cm-3) found in their samples indicates a higher concentration of cationic vacancies with nonintentional p-type doping. With the rising of the Pb fraction, the charge carriers concentration gradually increases and a maximum ZT of 0.2 at 700 K is obtained for x = 0.05. To our knowledge, only few thermoelectric oxysulfides were reported with rather low ZT values. For instance, Bi4O4S3 and LaPbBiS3O reach respectively 0.03 and 0.0023 at 300 K,65,66 whereas LaOBiS2 peaks around 0.07 at 743 K.67 The ZT value reported for Bi0.95Pb0.05CuOS in the present study is significantly higher and thereby, becomes the best ever achieved for an oxysulfide compound. Conclusion The thermoelectric properties of the Bi1-xPbxCuOS series (0 ≤ x ≤ 0.05) have been studied both experimentally and theoretically. The Pb for Bi substitution increases the charge carriers concentration and enables to increase substantially the electrical conductivity while moderately decreasing the Seebeck coefficient. This leads to a maximum power factor for x = 0.05 sample. The intrinsically low thermal conductivity provided by weak interatomic bonding between the layers and possible strong anharmonicity, is imperceptibly increased by the cationic substitution. Finally, the combination of low thermal conductivity with a power factor of 0.2 mW m-1 K-2 at 700 K increased the figure of merit ZT up to 0.2 in 5at.% Pbdoped BiCuOS sample. This value is 5 times higher with respect to BiCuOS. To the best of our knowledge, it is considered as the best ever reported value among the family of oxysulfides. Despite the obtained low ZT value, the approaches developed to optimize the electrical and thermal properties of BiCuOS appears as an interesting playground, which can be extended to other oxychalcogenides and materials in the future.

ACS Paragon Plus Environment

23

Chemistry of Materials 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 24 of 33

References (1)

Rowe, D. M. Thermoelectrics Handbook: Macro to Nano, ed. D. M.; CRC Press: Boca Raton, FL, 2006.

(2)

Ohtaki, M.; Araki, K.; Yamamoto, K. High Thermoelectric Performance of Dually Doped ZnO Ceramics. J. Electron. Mater. 2009, 38, 1234–1238.

(3)

Bérardan, D.; Guilmeau, E.; Maignan, A.; Raveau, B. In2O3:Ge, a Promising N-Type Thermoelectric Oxide Composite. Solid State Commun. 2008, 146, 97–101.

(4)

Lu, Z.; Zhang, H.; Lei, W.; Sinclair, D. C.; Reaney, I. M. High-Figure-of-Merit Thermoelectric La-Doped A-Site-De Ficient SrTiO3 Ceramics. Chem. Mater. 2016, 28, 925–935.

(5)

Bocher, L.; Aguirre, M. H.; Logvinovich, D.; Shkabko, A.; Robert, R.; Trottmann, M.; Weidenkaff, A. CaMn1-xNbxO3 (x < 0.08) Perovskite-Type Phases As Promising New High-Temperature N -Type Thermoelectric Materials. Inorg. Chem. 2008, 47 (18), 8077–8085.

(6)

Terasaki, I.; Sasago, Y.; Uchinokura, K. Large Thermoelectric Power in NaCo2O4 Single Crystals. Phys. Rev. B - Condens. Matter Mater. Phys. 1997, 56 (20), R12685– R12687.

(7)

Masset, A. C.; Michel, C.; Maignan, A.; Hervieu, M.; Toulemonde, O.; Studer, F.; Raveau, B.; Hejtmanek, J. Misfit-Layered Cobaltite with an Anisotropic Giant Magnetoresistance: Ca3Co4O9. Phys. Rev. B 2000, 62 (1), 166–175.

(8)

Li, S.; Funahashi, R.; Matsubara, I.; Ueno, K.; Yamada, H. High Temperature Thermoelectric Properties of Oxide CaCo9O28. J. Mater. Chem. 1999, 9, 1659–1660.

ACS Paragon Plus Environment

24

Page 25 of 33 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

Chemistry of Materials

(9)

Zhao, L. D.; Berardan, D.; Pei, Y. L.; Byl, C.; Pinsard-Gaudart, L.; Dragoe, N. Bi1−xSrxCuSeO Oxyselenides as Promising Thermoelectric Materials. Appl. Phys. Lett. 2010, 97, 92118-1-092118–3.

(10)

Kusainova, A. M.; Berdonosov, P. S.; Akselrud, L. G.; Kholodkosvkaya, L. N.; Dolgikh, V. A.; Popovkin, B. A. New Layered Compounds with the Genereal Composition (MO)(CuSe), Where M = Bi, Nd, Gd, Dy, and BiOCuS : Syntheses and Crystal Structure. J. Solid State Chem. 1994, 112, 189–191.

(11)

Li, F.; Li, J.-F.; Zhao, L.-D.; Xiang, K.; Liu, Y.; Zhang, B.-P.; Lin, Y.-H.; Nan, C.-W.; Zhu, H.-M. Polycrystalline BiCuSeO Oxide as a Potential Thermoelectric Material. Energy Environ. Sci. 2012, 5, 7188–7195.

(12)

Li, J.; Sui, J.; Pei, Y.; Meng, X.; Berardan, D.; Dragoe, N.; Cai, W.; Zhao, L.-D. The Roles of Na Doping in BiCuSeO Oxyselenides as a Thermoelectric Material. J. Mater. Chem. A 2014, 2, 4903–4906.

(13)

Sun Lee, D.; An, T.-H.; Jeong, M.; Choi, H. S.; Lim, Y. S.; Seo, W.-S.; Park, C.-H.; Park, C.; Park, H.-H. Density of State Effective Mass and Related Charge Transport Properties in K-Doped BiCuOSe. Appl. Phys. Lett. 2013, 103, 232110.

(14)

Li, J.; Sui, J.; Barreteau, C.; Berardan, D.; Dragoe, N.; Cai, W.; Pei, Y.; Zhao, L.-D. Thermoelectric Properties of Mg Doped P -Type BiCuSeO Oxyselenides. J. Alloys Compd. 2013, 551, 649–653.

(15)

Pei, Y.-L.; He, J.; Li, J.-F.; Li, F.; Liu, Q.; Pan, W.; Barréteau, C.; Bérardan, D.; Dragoe, N.; Zhao, L.-D. High Thermoelectric Performance of Oxyselenides : Intrinsically Low Thermal Conductivity of Ca-Doped BiCuSeO. NPG Asia Mater. 2013, 5, e47.

ACS Paragon Plus Environment

25

Chemistry of Materials 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

(16)

Page 26 of 33

Feng, D.; Zheng, F.; Wu, D.; Wu, M.; Li, W.; Huang, L.; Zhao, L.; He, J. Investigation into the Extremely Low Thermal Conductivity in Ba Heavily Doped BiCuSeO. Nano Energy 2016, 27, 167–174.

(17)

Li, J.; Sui, J.; Pei, Y.; Barreteau, C.; Berardan, D.; Dragoe, N.; Cai, W.; He, J.; Zhao, L.-D. A High Thermoelectric Figure of Merit ZT > 1 in Ba Heavily Doped BiCuSeO Oxyselenides. Energy Environ. Sci. 2012, 5, 8543–8547.

(18)

Luu, S. D. N.; Vaqueiro, P. Synthesis, Structural Characterisation and Thermoelectric Properties of Bi1-xPbxOCuSe. J. Mater. Chem. A 2013, 1, 12270–12275.

(19)

Liu, Y.; Zhao, L.; Liu, Y.; Lan, J.; Xu, W.; Li, F.; Zhang, B.; Berardan, D. Remarkable Enhancement in Thermoelectric Performance of BiCuSeO by Cu Deficiencies. J. Am. Chem. Soc. 2011, 133, 20112–20115.

(20)

Pei, Y.-L.; Wu, H.; Wu, D.; Zheng, F.; He, J. High Thermoelectric Performance Realized in a BiCuSeO System by Improving Carrier Mobility through 3D Modulation Doping. J. Am. Chem. Soc. 2014, 136 (39), 13902–13908.

(21)

Liu, Y.; Lan, J.; Xu, W.; Liu, Y.; Pei, Y.-L.; Cheng, B.; Liu, D.-B.; Hin, Y.-H.; Zhao, L.-D. Enhanced Thermoelectric Performance of a BiCuSeO System via Band Gap Tuning. Chem. Commun. 2013, 49, 8075–8077.

(22)

Sui, J.; Li, J.; He, J.; Pei, Y.-L.; Berardan, D.; Wu, H.; Dragoe, N.; Cai, W.; Zhao, L.D. Texturation Boosts the Thermoelectric Performance of BiCuSeO Oxyselenides. Energy Environ. Sci. 2013, 6 (10), 2916.

(23)

Liu, Y.; Zhao, L.-D.; Zhu, Y.; Liu, Y.; Li, F.; Yu, M.; Liu, D.; Xu, W.; Lin, Y.-H.; Nan, C.-W. Synergistically Optimizing Electrical and Thermal Transport Properties of BiCuSeO via a Dual-Doping Approach. Adv. Energy Mater. 2016, 6, 1502423.

ACS Paragon Plus Environment

26

Page 27 of 33 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

Chemistry of Materials

(24)

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, 7837–7846.

(25)

Zhu, H.; Su, T.; Li, H.; Pu, C.; Zhou, D.; Zhu, P.; Wang, X. High Pressure Synthesis, Structure and Thermoelectric Properties of BiCuChO (Ch = S, Se, Te). J. Eur. Ceram. Soc. 2017, 37, 1541–1546.

(26)

Goldsmid, H. J.; Sharp, J. W. Estimation of the Thermal Band Gap of a Semiconductor from Seebeck Measurements. J. Electron. Mater. 1999, 28 (7), 869–872.

(27)

Hiramatsu, H.; Yanagi, H.; Kamiya, T.; Ueda, K.; Hirano, M.; Hosono, H. Crystal Structures, Optoelectronic Properties, and Electronic Structures of Layered Oxychalcogenides MCuOCh ( M = Bi,La ; Ch = S, Se, Te): Effects of Electronic Configurations of M3+ Ions. Chem. Mater. 2008, 20, 326–334.

(28)

Pal, A.; Kishan, H.; Awana, V. P. S. Synthesis and Structural Details of BiOCu1−xS : Possible New Entrant in a Series of Exotic Superconductors ? J. Supercond. Nov. Magn. 2010, 23, 301–304.

(29)

Berthebaud, D.; Guilmeau, E.; Lebedev, O. I.; Maignan, A.; Gamon, J.; Barboux, P. The BiCu1−xOS Oxysulfide: Copper Deficiency and Electronic Properties. J. Solid State Chem. 2016, 237, 292–299.

(30)

Rodríguez-Carvajal, J. Recent Advances in Magnetic Structure Determination by Neutron Powder Diffraction. Phys. B Condens. Matter 1993, 192 (1–2), 55–69.

(31)

Roisnel, T.; Rodríguez-Carvajal, J. WinPLOTR: A Windows Tool for Powder Diffraction Pattern Analysis. Mater. Sci. Forum 2001, 378–381, 118–123.

ACS Paragon Plus Environment

27

Chemistry of Materials 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

(32)

Page 28 of 33

Alleno, E.; Bérardan, D.; Byl, C.; Candolfi, C.; Daou, R.; Decourt, R.; Guilmeau, E.; Hébert, S.; Hejtmanek, J.; Lenoir, B.; Masschelein, P.; Ohorodnichuk, V.; Pollet, M.; Populoh, S.; Ravot, D.; Rouleau, O.; Soulier, M. A Round Robin Test of the Uncertainty on the Measurement of the Thermoelectric Dimensionless Figure of Merit of Co0.97Ni0.03Sb3. Rev. Sci. Instrum. 2015, 86, 11301.

(33)

Clark, S. J.; Segall, M. D.; Pickard, C. J.; Hasnip, P. J.; Probert, M. I. J.; Refson, K.; Payne, M. C. First Principles Methods Using CASTEP. Z. Krist. 2005, 220, 567–570.

(34)

Vanderbilt, D. Soft Self-Consistent Pseudopotentials in a Generalized Eigenvalue Formalism. Phys. Rev. B 1990, 41 (11), 7892–7895.

(35)

Garrity, K. F.; Bennett, J. W.; Rabe, K. M.; Vanderbilt, D. Pseudopotentials for HighThroughput DFT Calculations. Comput. Mater. Sci. 2014, 81, 446–452.

(36)

Perdew, J. P.; Ruzsinszky, A.; Csonka, G. I.; Vydrov, O. A.; Scuseria, G. E.; Constantin, L. A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. Phys. Rev. Lett. 2008, 100, 136406.

(37)

Blaha, P.; Schwarz, K.; Madsen, G.; Kvasnicka, D.; Luitz, J. WIEN2k, an Augmented Plane Wave plus Local Orbitals Program for Calculating Crystal Properties, Techn. Uni.; Schwarz, K., Ed.; 2001.

(38)

Tran, F.; Blaha, P. Accurate Band Gaps of Semiconductors and Insulators with a Semilocal Exchange-Correlation Potential. Phys. Rev. Lett. 2009, 102 (22), 226401.

(39)

Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and Related Crystal Properties from Density-Functional Perturbation Theory. Rev. Mod. Phys. 2001, 73, 515–562.

ACS Paragon Plus Environment

28

Page 29 of 33 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

Chemistry of Materials

(40)

Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I.; Dal Corso, A.; de Gironcoli, S. QUANTUM ESPRESSO : A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys. Condens. Matter 2009, 21, 395502.

(41)

Johnson, V.; Jeitschko, W. ZrCuSiAs : A “ Filled ” PbFCl Type. J. Solid State Chem. 1974, 11, 161–166.

(42)

An, T.-H.; Lim, Y. S.; Choi, H.-S.; Seo, W.-S.; Park, C.-H.; Kim, G.-R.; Park, C.; Lee, C. H.; Shim, J. H. Point Defect-Assisted Doping Mechanism and Related Thermoelectric Transport Properties in Pb-Doped BiCuOTe. J. Mater. Chem. A 2014, 2, 19759–19764.

(43)

Sheets, W. C.; Stampler, E. S.; Kabbour, H.; Bertoni, M. I.; Cario, L.; Mason, T. O.; Marks, T. J.; Poeppelmeier, K. R. Facile Synthesis of BiCuOS by Hydrothermal Methods. Inorg. Chem. 2007, 46 (25), 10741–10748.

(44)

Ueda, K.; Hosono, H. Crystal Structure of LaCuOS1-xSex Oxychalcogenides. Thin Solid Films 2002, 411, 115–118.

(45)

Hiramatsu, H.; Kamiya, T.; Tohei, T.; Ikenaga, E.; Mizoguchi, T.; Ikuhara, Y.; Kobayashi, K.; Hosono, H. Origins of Hole Doping and Relevant Optoelectronic Properties of Wide Gap P-Type Semiconductor , LaCuOSe. J. Am. Chem. Soc. 2010, 132 (42), 15060–15067.

(46)

Shannon, R. D. Revised Effective Ionic Radii and Systematic Studies of Interatomie Distances in Halides and Chaleogenides. Acta Crystallogr. 1976, A32, 751–767.

(47)

Barreteau, C.; Bérardan, D.; Amzallag, E.; Zhao, L.; Dragoe, N. Structural and Electronic Transport Properties in Sr-Doped BiCuSeO. Chem. Mater. 2012, 24, 3168–

ACS Paragon Plus Environment

29

Chemistry of Materials 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 30 of 33

3178. (48)

Dong, S.-T.; Lv, Y.-Y.; Zhang, B.-B.; Zhang, F.; Yao, S.; Chen, Y. B.; Zhou, J.; Zhang, S.-T.; Gu, Z.-B.; Chen, Y.-F. Strong Correlation of the Growth Mode and Electrical Properties of BiCuSeO Single Crystals with Growth Temperature. Cryst. Eng. Comm 2015, 17 (32), 6136–6141.

(49)

Pan, L.; Bérardan, D.; Zhao, L.; Barreteau, C.; Dragoe, N. Influence of Pb Doping on the Electrical Transport Properties of BiCuSeO. Appl. Phys. Lett. 2013, 102, 23902.

(50)

Debye, P. P.; Conwell, E. M. Electrical Properties of N-Type Germanium. Phys. Rev. 1953, 93 (4), 693–706.

(51)

Richard, A. P.; Russell, J. A.; Zakutayev, A.; Zakharov, L. N.; Keszler, D. A.; Tate, J. Synthesis, Structure, and Optical Properties of BiCuOCh (Ch=S, Se, and Te). J. Solid State Chem. 2012, 187, 15–19.

(52)

Lan, J.-L.; Liu, Y.-C.; Zhan, B.; Lin, Y.-H.; Zhang, B.; Yuan, X.; Zhang, W.; Xu, W.; Nan, C.-W. Enhanced Thermoelectric Properties of Pb-Doped BiCuSeO Ceramics. Adv. Mater. 2013, 25, 5086–5090.

(53)

Yang, Y.; Liu, X.; Liang, X. Thermoelectric Properties of Bi1-xSnxCuSeO Solid Solutions. Dalt. Trans. 2017, 46, 2510–2515.

(54)

Al Rahal Al Orabi, R.; Hwang, J.; Lin, C. C.; Gautier, R.; Fontaine, B.; Kim, W.; Rhyee, J. S.; Wee, D.; Fornari, M. Ultralow Lattice Thermal Conductivity and Enhanced Thermoelectric Performance in SnTe:Ga Materials. Chem. Mater. 2017, 29, 612–620.

(55)

Vaqueiro, P.; Al Rahal Al Orabi, R.; Luu, S. D. N.; Guélou, G.; Powell, A. V.; Smith,

ACS Paragon Plus Environment

30

Page 31 of 33 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

Chemistry of Materials

R. I.; Song, J.-P.; Wee, D.; Fornari, M. The Role of Copper in the Thermal Conductivity of Thermoelectric Oxychalcogenides: Do Lone Pairs Matter? Phys. Chem. Chem. Phys. 2015, 17, 31735–31740. (56)

May, A. F.; Toberer, E. S.; Saramat, A.; Snyder, G. J. Characterization and Analysis of Thermoelectric Transport in N-Type Ba8Ga16-xGe30+x. Phys. Rev. B 2009, 80, 125205.

(57)

Kim, H. S.; Gibbs, Z. M.; Tang, Y.; Wang, H.; Snyder, G. J. Characterization of Lorenz Number with Seebeck Coefficient Measurement. APL Mater. 2015, 3, 41506.

(58)

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 (7496), 373–377.

(59)

Al Rahal Al Orabi, R.; Orisakwe, E.; Wee, D.; Fontaine, B.; Gautier, R.; Halet, J.-F.; Fornari, M. Prediction of High Thermoelectric Potential in AMN2 Layered Nitrides: Electronic Structure, Phonons, and Anharmonic Effects. J. Mater. Chem. A 2015, 3 (18), 9945–9954.

(60)

Saha, S. K. Exploring the Origin of Ultralow Thermal Conductivity in Layered BiOCuSe. Phys. Rev. B 2015, 92, 041202(R).

(61)

Pan, L.; Xia, Q.; Ye, S.; Ding, N.; Liu, Z. First Principles Study of Electronic Structure, Chemical Bonding and Elastic Properties of BiOCuS. Trans. Nonferrous Met. Soc. China 2012, 22, 1197–1202.

(62)

Liu, G.; Sun, H.; Zhou, J.; Li, Q.; Wan, X. G. Thermal Properties of Layered Oxychalcogenides BiCuOCh (Ch = S, Se, and Te): A First-Principles Calculation. J. Appl. Phys. 2016, 119, 185109.

ACS Paragon Plus Environment

31

Chemistry of Materials 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

(63)

Page 32 of 33

Callaway, J. Model for Lattice Thermal Conductivity at Low Temperatures. Phys. Rev. 113AD, 4, 1046–1051.

(64)

Cahill, D. G.; Watson, S. K.; Pohl, R. O. Lower Limit to the Thermal Conductivity of Disordered Crystals. Phys. Rev. B 1992, 46, 6131–6140.

(65)

Tan, S. G.; Li, L. J.; Liu, Y.; Tong, P.; Zhao, B. C.; Lu, W. J.; Sun, Y. P. Superconducting and Thermoelectric Properties of New Layered Superconductor Bi4O4S3. Phys. C Supercond. its Appl. 2012, 483, 94–96.

(66)

Mizuguchi, Y.; Nishida, A.; Omachi, A.; Miura, O. Thermoelectric Properties of New Bi-Chalcogenide Layered Compounds. Cogent Phys. 2016, 3, 1156281.

(67)

Mizuguchi, Y.; Omachi, A.; Goto, Y.; Kamihara, Y.; Matoba, M.; Hiroi, T.; Kajitani, J.; Miura, O. Enhancement of Thermoelectric Properties by Se Substitution in Layered Bismuth- Chalcogenide LaOBiS2-xSex. J. Appl. Phys. 2014, 116, 163915.

ACS Paragon Plus Environment

32

Page 33 of 33 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

Chemistry of Materials

Table of Content Thermoelectric properties in the Bi1-xPbxCuOS series have been studied both experimentally and theoretically. A record ZT of 0.2 is obtained in the family of oxysulfides.

ACS Paragon Plus Environment

33