Sb Doping of Metallic CuCr2S4 as a Route to Highly Improved

Feb 17, 2017 - Somnath Acharya , Sharmistha Anwar , Takao Mori , Ajay Soni ... Kan Chen , Cono Di Paola , Baoli Du , Ruizhi Zhang , Savio Laricchia , ...
0 downloads 0 Views 6MB Size
Article pubs.acs.org/cm

Sb Doping of Metallic CuCr2S4 as a Route to Highly Improved Thermoelectric Properties Atta Ullah Khan,† Rabih Al Rahal Al Orabi,‡,§ Amir Pakdel,† Jean-Baptiste Vaney,† Bruno Fontaine,∥ Régis Gautier,∥ Jean-François Halet,∥ Seiji Mitani,⊥,# and Takao Mori*,†,# †

National Institute for Materials Science (NIMS), MANA, 1-1-1 Namiki, Tsukuba 305-0044, Japan Department of Environmental Science and Engineering, Ewha Womans University, Seoul 120-750, Korea § Department of Physics and Science of Advanced Materials, Central Michigan University, Mt. Pleasant, Michigan 48858, United States ∥ Institut des Sciences Chimiques de Rennes, UMR 6226 CNRSUniversité de Rennes 1Ecole Nationale Supérieure de Chimie de Rennes, CS 50837, 35708 Rennes, France ⊥ National Institute for Materials Science (NIMS), 1-2-1 Sengen, Tsukuba 305-0047, Japan # Graduate School of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan ‡

S Supporting Information *

ABSTRACT: We report for the first time the thermoelectric properties of CuCr2−xSbxS4 (0.22 ≤ x ≤ 0.5). Although CuCr2S4 has been reported to be a metallic compound, addition of Sb shifts the material toward the semiconductor side. This is confirmed by band structure calculations of CuCr2−xSbxS4 (x = 0, 0.25, 0.5) models. Increasing Sb content enhances the power factor. However, beyond x = 0.3, further Sb addition lowers the electrical conductivity and power factor. A very interesting point is the simultaneous increase of the Seebeck coefficient as well as the electrical conductivity with increasing temperature, which acts like a variable range hopping (VRH) compounds but possesses much better properties than those having VRH. Samples were annealed for 48 h prior to thermoelectric properties measurements to have a reliable dimensionless figure of merit (ZT). An attractive ZT of 0.43 is obtained at ∼650 °C. The attractive thermoelectric properties we discovered by driving a metal compound into a semiconductor make this compound an interesting thermoelectric material especially because of the cheap constituent elements compared to those of typical state-of-the-art thermoelectric materials. Furthermore, this material is stable up to 650 °C at least, a relatively high temperature for sulfides. Additionally, we discovered a miscibility gap in this solid solution close to an Sb content of 0.15; although a detailed study dedicated entirely to this miscibility gap would be required, it will encourage the researchers to further explore this system. efforts are ongoing to find new and more efficient thermoelectric materials.4−7 Efficiency of thermoelectric materials is measured by the dimensionless figure of merit (ZT). It is calculated by the formula ZT = S2σT/κ, where S stands for the Seebeck coefficient, σ represents the electrical conductivity, and κ is the

1. INTRODUCTION With increasing consumption of fossil fuels and a relentless depletion of these resources, it has become very important to convert them more efficiently into usable forms of energy. About two-thirds of the energy of fossil fuels is wasted mainly in the form of waste heat.1 It is possible to convert this waste heat into useful energy by employing thermoelectric materials.2,3 This is why research focused on thermoelectric materials has increased by many folds in recent times and why intense © 2017 American Chemical Society

Received: December 19, 2016 Revised: February 16, 2017 Published: February 17, 2017 2988

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Article

Chemistry of Materials thermal conductivity. S2σ as a whole is termed as the power factor. As obvious from the formula, the power factor should be as high as possible while still keeping the thermal conductivity as low as possible. Different thermoelectric materials are required for different temperature ranges and applications. Most middle or low temperature range materials are chalcogenides or antimonides.5−7 Out of chalcogenides, sulfur is the cheapest and lightest element, which makes sulfides attractive compounds for bulk thermoelectric materials. Moreover, sulfur is one of the byproducts of oil exploration, and its production is far more than its consumption. Companies are piling up sulfur in the forms of pyramids,8 posing a huge potential threat to the environment. Therefore, it is very important to utilize sulfur to make some valuable products such as thermoelectric materials. Indeed, several sulfide compounds have been discovered with good thermoelectric properties.9−13 However, naturally occurring minerals, if suitable for thermoelectric applications, are the best option.11,13,14 In addition to a newly discovered sulfospinel cuprokalininite (CuCr2S4) found from metamorphic rocks in the South Baikal region of Russia,15 another mineral named florensovite with the formula Cu(Cr1.5Sb0.5)S4 is already reported to be found from the Baikal region as well.16 Sbdoped CuCr2S4 has been investigated regarding its interesting magnetic properties: CuCr 1.5 Sb 0.5S 4 is reported to be antiferromagnetic,17 while CuCr2S4 is described as ferromagnetic.18 Moreover, giant magnetoresistance has been reported for CuCr1.6Sb0.4S4.19 However, to the best of our knowledge, thermoelectric properties were not investigated for Sb-doped CuCr2S4. On the basis of these previous reports,17−19 we assume that replacing some Cr by Sb can cause a significant influence on thermoelectric properties. Although CuCr2S4 is reported to be metallic,20−22 a detailed literature survey did not show any report related to the effect of Sb addition on the thermoelectric properties except that CuCr1.5Sb0.5S4 is a semiconductor.17 In this study, for the first time, we systematically and comprehensively investigate the band structure and high temperature thermoelectric properties of CuCr2−xSbxS4 using experimental and theoretical tools. Because undoped CuCr2S4 is a metallic compound with the Fermi level below the top of the valence bands, we show that Sb doping provides extra electrons to fill some of the holes in the valence bands, leading to an enhanced Seebeck coefficient and reduction of the magnetic moment. We also show that the substitution of Sb with Cr atoms in CuCr2S4 leads to a high power factor for x = 0.3. This yields to a high figure of merit ZT equal to 0.43 at T ∼ 650 °C, which is one of the highest values reported for sulfurbased thermoelectric materials.

vacuum. Highly dense (∼97% of the theoretical density) disk-shaped pellets with dimensions of 10 mm diameter and 2 mm thickness were obtained. The obtained pellets were sealed in quartz tube and annealed at 650 °C for 48 h. Powder X-ray Diffraction (PXRD) and High Resolution Transmission Electron Microscopy (HRTEM). Samples were characterized by PXRD with Cu Kα (λ = 1.54056 Å) radiation operated at room temperature using Ultima III (Rigaku corporation, Tokyo, Japan) instrument. Lattice constants, volumes, and theoretical densities were calculated by the Rietveld refinement using the Fullprof Suite software. HRTEM was carried out using a JEM-2100F microscope (JEOL Ltd., Japan) equipped with an energy dispersive spectrometer (EDS). Electronic Transport Measurements. The sintered pellets were cut into ∼1.8 × 2.0 × 8 mm bars for simultaneous measurement of the electrical resistivity and the Seebeck coefficient. Measurements were performed under He atmosphere from room temperature to 873 K using a ULVAC-Riko-ZEM-5 instrument. Thermal Transport Measurements. The sintered pellets were polished into a cylinder shape ∼10 and ∼1.7 mm thick for thermal diffusivity measurements. The thermal diffusivity coefficient (D) was measured with an ULVAC TC-7000 instrument. The heat capacity (Cp) was measured using differential scanning calorimetry (DSC, Netzsch AST 449). The total thermal conductivity, k, was calculated using the formula k = DdCp, where d represents the density of the sample. The electronic thermal conductivity is proportional to the electrical conductivity according to the Wiedemann−Franz law (ke = σLT). The Lorenz number, L, was estimated as a function of temperature from the experimental Seebeck coefficients using the single parabolic band model and assuming that acoustic phonon scattering limits the mobility, as explained in previous literature.23 The electronic thermal conductivity was subtracted from the total thermal conductivity to obtain lattice thermal conductivity. Error bars were determined from measurements of different samples of the same composition. Hall Effect Measurements. The Hall coefficient RH was measured at 300 K on the sample CuCr1.5Sb0.5S4 with the AC transport option of the Quantum Design Physical Properties Measurement System (PPMS). A five-wire configuration was used by attaching Au wires with Ag epoxy onto the sample. The magnetic field was swept between −3 and +3 T. The hole concentration p and Hall mobility μH were estimated from p = (1/e)RH and μH = RH/ρ. The Hall effect could not be measured on the other samples due to parasitic effects owing to the magnetic nature of these latter compounds. First-Principles Calculations. Spin-polarized density functional theory (DFT) geometry optimizations of CuCr2−xSbxS4 (x = 0, 0.25, and 0.5) were carried out with the CASTEP8.0 code24 using the GGA in the parametrization of PBE functional.25 Cell parameters and atomic positions were both relaxed. Because of the presence of localized d electrons, an additional Hubbard-like term was introduced for Cr and Cu. The simplified Dudarev approach26 was used with Ueff = U − J = 3.5 eV for Cu 3d orbitals and 2.5 eV on the Cr 3d orbitals.27 All ultrasoft pseudopotentials were generated using the OTF_ultrasoftpseudopotential generator included in the program. The cutoff energy for plane-waves was set to 500 eV. The electronic wave function was sampled with 30 k-points in the first Brillouin zone using the Monkhorst−Pack method.28 For the electronic band structures, we used the full-potential linearized augmented plane wave (FLAPW) approach, as implemented in the WIEN2K code.29 The PBE+U approach was used with the Ueff values for 3d orbitals of Cu and Cr previously mentioned. A plane wave cutoff corresponding to RMTKmax = 7 was used. The radial wave functions inside the nonoverlapping 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 a Brillouin zone integration mesh of 500 k-points.

2. EXPERIMENTAL PROCEDURES Sample Preparation. CuCr2−xSbxS4 (0 ≤ x ≤ 0.5) samples were prepared by sealing high purity powders of Cr, Sb, Cu, and S in quartz tubes under vacuum. Tubes were heated from room temperature to 750 °C in 8 h, kept at 750 °C for 15 h, and then slowly cooled. The materials were ground in SiC mortars, compressed in a cold isostatic press (CIP) machine, sealed again in quartz under vacuum, heated to 650 °C, and kept there for 15 h to remove any stress induced during pressurized sintering. The obtained samples of CuCr2−xSbxS4 were ground into fine powders using a mortar and pestle to reduce the grains. The powders were then densified by a spark plasma sintering (SPS, Dr Sinter 1080, SPS-Syntex Inc., Tokyo, Japan) method at 720 °C for 5 min in a 10 mm diameter graphite die under an axial pressure of 80 MPa in

3. RESULTS AND DISCUSSION Powder X-ray diffraction measurements followed by the Rietveld refinement confirmed the crystallinity, single phase 2989

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Article

Chemistry of Materials

Figure 1. Rietveld refinement of the XRD pattern of CuCr1.5Sb0.5S4 sample, indicating a single phase material and confirming MgAl2O4-type structure. Reasonably low R factors indicate a good fit of the observed and calculated patterns.

Figure 2. Lattice parameters vs Sb content. Slight deviation from Vegard’s law is observed. Lattice parameters of a sintered pellet having Sb content of 0.17 after annealing showed a clear split (green triangles), indicating the presence of a miscibility gap at 600 °C.

formation, and replacement of Cr by Sb, forming a CuCr2−xSbxS4 solid solution. Figure 1 shows the Rietveld refinement of the XRD pattern of CuCr1.5Sb0.5S4, indicating a single-phase sample. Lattice parameters obtained (a = 1.0001 ± 0.0004 nm) were in good agreement with those reported in the literature (1.0018 nm).30 Figure 2 shows the graph of lattice parameters vs Sb content. Owing to the size of Sb much larger than that of Cr, their covalent radii are equal to 1.40 and 1.22 Å, respectively, it is obvious that the replacement of Cr by Sb will expand the lattice, resulting in an increase in lattice constants. Peaks in the XRD pattern shift to the lower 2θ, in accordance

with the lattice parameters expansion. In general, a slight positive deviation from Vegard’s law is observed. Lattice parameters of a sintered pellet with composition CuCr2−xSbxS4 (x = 0.17) after annealing show a clear split (for XRD pattern, see Supporting Information, Figure S2), indicating the presence of a miscibility gap at 600 °C. This miscibility gap is observed for the first time, as no systematic study is reported on this solid solution. High crystallinity as well as the presence of a miscibility gap were also confirmed by HRTEM. Images obtained by HRTEM, performed on a sample with an Sb content x = 0.3 (Figures 3a− 2990

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Article

Chemistry of Materials

Figure 3. HRTEM images of x = 0.3 (a−c), indicating a highly crystalline material and presenting the grain boundaries (a and b). Inset in panel c shows the interplanar distance taken for 10 planes, and an average value produced the interplanar distance for this particular plane. STEM helped us confirm the presence of a miscibility gap in the x = 0.17 sample by EDX analysis of points 1 and 2 (d). EDX spectra of points 1 (e) and 2 (f) clearly showed a variation in Cr and Sb content ratio.

metallic.22 Therefore, we focused mainly on the Sb-rich side of the solid solution. It has been previously reported that Sb possesses a valence state similar to that of Cr (+3) in CuCr2−xSbxS4.32,33 This disagrees with the charge of +3.53 computed for Sb in CuCr1.5Sb0.5S4, which suggests a +5 valence state for the pnictogen. Indeed, Warczewski et al. previously experimentally proposed a +5 state.17 The replacement of Cr by Sb hardly affects the local crystal structure. Moreover, shorter than expected Sb−S distances (Supporting Information, Table S1) indicate a rather strong bond.32 This could be a possible reason for a highly stable Sb-rich side of the solid solution under vacuum. The crystal structure of CuCr1.5Sb0.5S4 is shown in Figure 4. This structure consists of tetrahedra and octahedra made of sulfur. One eighth of the total tetrahedra are filled by Cu, and half of the total octahedra are filled by a Cr/Sb mixture.

c), indicated a highly crystalline, defect-free single-phase material. No intergranular phase was observed. However, scanning transmission electron microscopy (STEM), performed on the sample with Sb = 0.17, helped us to confirm the presence of a miscibility gap in this sample, which was already observed in the X-ray powder diffraction patterns. Energy dispersive X-ray (EDX) analysis (Figures 3e and f) of the two grains unambiguously showed a variation in Cr and Sb content ratio. We decided to focus on the thermoelectric properties of the Sb-rich side of the miscibility gap for two reasons: (a) the decomposition temperature of CuCr2S4 under vacuum is reported to be 445−550 °C.31 TC-7000 runs under continuous pumping and caused the samples with Sb content 0 and 0.075 to decompose partially. (b) CuCr2S4 is already reported to be 2991

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Article

Chemistry of Materials

as in the present case may have its role in low lattice thermal conductivity as well. HRTEM images (Figure 3a−c) confirmed that the low thermal conductivity observed in the present study was not a result of processing defects; instead, it is more likely an intrinsic property. Considering the electronic properties, addition of Sb is expected to lower the electrical conductivity compared to that of the metallic CuCr2S4 with a corresponding rise in the Seebeck coefficient. Variation of magnetic properties that suggest ferromagnetic (ferrimagnetic) to antiferromagnetic transition with increasing Sb content was also observed (Supporting Information, Figure S1). Figure 5 shows the measured thermoelectric properties of CuCr2−xSbxS4 (0.22 ≤ x ≤ 0.5). Replacement of Cr with Sb decreases the electrical conductivity and makes it a semiconductor compound. From Figure 5a, it is obvious that the electrical conductivity is inversely proportional to the Sb content at low temperatures and decreases with increasing Sb. Regarding the temperature dependence, generally, the electrical conductivity increases as the temperature increases. CuCr2S4 is metallic at room temperature with only a small Seebeck coefficient of 16 μV/K.22 Although, the mineral florensovite CuCr1.5Sb0.5S4 is reported to be a semiconductor, no information regarding the Seebeck coefficient is given yet.17 In current study, we observed that the addition of Sb improves the Seebeck coefficient significantly (Figure 5b). We investigated the thermoelectric properties of this composition and discovered a large Seebeck coefficient of >200 μV/K. Positive

Figure 4. Representation of the crystal structure of CuCr1.5Sb0.5S4. Pink tetrahedra are filled by Cu, while Cr/Sb occupies the center of the octahedra.

The rest of the tetra- and octahedra are empty. Spitzer34 suggested that the spinel structure, having primarily octahedral coordination and partial filling of the polyhedra, tends to have lower lattice thermal conductivity. Additionally, a large unit cell

Figure 5. Temperature dependence of the thermoelectric properties of CuCr2−xSbxS4 (0.22 ≤ x ≤ 0.5). (a) Electrical conductivity (σ) seems to drop rapidly after a certain Sb content. (b) The Seebeck coefficient is positive in the whole temperature range, pointing toward p-type material with holes as dominant charge carriers. (c) Total thermal conductivity (κ) with similar values at low temperatures but differentiating considerably at high temperatures. (d) Lattice thermal conductivity (κlat). (e) Power factor, which increases initially with increasing Sb content, reaches a maximum, and then decreases. (f) Figure of merit (ZT) with a highest value of 0.43 for Sb = 0.3. 2992

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Article

Chemistry of Materials

increase in both conductivity and Seebeck coefficient with temperature, the degenerate semiconductor nature may be playing a role, although systems such as Bi2Te3 often show metallic temperature dependence of the conductivity. The relatively low values of mobility indicate the presence of disorder, which is consistent with the crystal structure features mentioned above and the low thermal conductivity. This may also be affecting the temperature dependence of the conductivity, although it is not likely to be a conventional VRH system. The thermoelectric characteristics are overall attractive, and further investigations on the transport of other systems with high carrier concentration and some degree of disorder are merited. Computing the experimentally measured value would give a high effective mass of 5.9m0 with m0 being the electron mass.23 This value appears unexpectedly high, and a breakdown of the above formulation may be occurring, but an enhanced effective mass could nevertheless justify the combination of the high carrier concentration and the high thermopower. In the case of chalcopyrite, magnetic interaction was proposed as the origin of the heavy effective mass and large Seebeck coefficient (large power factor).11,13 Because the present system is also magnetic, a similar situation may be applicable; however, this case indicates that the ordering itself of the magnetic moments may not be critical because the ordering temperature is significantly lower. The explicit mechanism which can be related to, for example, existence of sizable magnetic moments in the system, spin fluctuations, and magnetic interactions cannot be rigorously proved yet. Probing the detailed linkage between the thermoelectric and magnetic properties remains an intriguing research topic to be pursued further. In any case, a heavy effective mass is indicated, and this is also consistent with the DFT calculations described next. To rationalize these physical measurements, DFT calculations were carried out. As shown by the X-ray diffraction experiments, substitution of Cr by Sb in the compound does not modify the description of the crystal structure in the cubic space group Fd3̅m. However, to model the substitution of Cr by Sb atoms in the compound, the crystal structure of CuCr2−xSbxS4 must be described in a lower symmetrical space group. Among the two rates of substitution of Cr with Sb that have been considered (0.25 and 0.5), several distributions can be envisioned. To easily compare the electronic structure of CuCr2−xSbxS4 (x = 0, 0.25, and 0.5), the monoclinic space group P2/m was considered for all models. As Gupta et al. did to study the electronic structure of CoxCu1−xCr2S4 nanocrystals,27 GGA+U calculations were carried out. Table 2 displays the optimized unit cell parameters and pertinent distances. As expected from the larger size of Sb compared to that of Cr and as observed experimentally, substitution of Cr with Sb increases the volume of the unit cell. The DFT optimized unit cell parameters are slightly larger than those resulting from the X-

values suggest that it is a p-type material. However, it possesses a very low electrical conductivity compared to that of the parent compound CuCr2S4 (1.11 × 105 S/m).22 It indicates that the Seebeck coefficient is quite sensitive to Sb content. Like the electrical conductivity, the Seebeck coefficient shows similar temperature dependence and increases in all cases with rising temperature. The only exception is the sample with the highest Sb content x = 0.5, in which it saturates around 350 °C, stays more or less unchanged up to 500 °C, and then starts decreasing. According to the Boltzmann equation, electrical conductivity and the Seebeck coefficient have an inverse relationship, so an improvement in one will result in a degradation of the other. However, the ideal case would be an enhancement in both properties simultaneously, and among typically studied thermoelectric material, there are not many cases where it is observed that both properties increase.35−38 Some notable examples are the VRH systems such as the higher borides, where both properties increase with increasing temperature.37−41 However, these are materials that are typically only thermoelectrically competitive at high temperatures. The current compound under study also shows this kind of behavior, especially in the case of an Sb content of 0.3 where both electrical conductivity and Seebeck coefficient are increasing almost linearly with increasing temperature. However, in this case, considering the relatively high electrical conductivity, it is unlikely to be a VRH system. Furthermore, the overall thermoelectric properties of the present system appear to be superior compared to those VRH systems, making this compound very attractive for further studies. As discussed above, addition of Sb lowers the electrical conductivity but enhances the Seebeck coefficient, so it is necessary to find a compromise between electrical conductivity and the Seebeck coefficient. A better term to represent is the power factor (S2σ), which should be optimized. Figure 5e indicates that the addition of Sb increases the power factor to a certain value and then starts decreasing. The optimized value of Sb from the graph is found to be x = 0.30. A fairly high power factor of ∼1 mW/mK2 is observed at 643 °C. It must be noted that the properties were measured on annealed samples (annealed at 650 °C after sintering) to remove any stress induced due to pressurized sintering. This kind of stress sometimes results in special properties which dissipate after a few runs. Therefore, our reported properties are supposed to stay more or less the same over multiple runs. The Hall effect was measured for the sample CuCr1.5Sb0.5S4, as interfering magnetic effects in the other compounds would not allow a reliable measurement. The measured Hall coefficient and calculations of the hole carrier concentration and hole mobility values are shown in Table 1. This sample shows a very high carrier concentration at 300 K of 2.41 × 1021 cm−3, a high value consistent with a metallic or degenerate semiconducting nature of the compound properties. The hole mobility of 4.06 cm2/V·s is rather low compared to those of other state of art thermoelectric compounds that can exceed 100 cm2/V·s. Regarding the trend discussed above and the

Table 2. Optimized Cell Parameters and Pertinent Distances (Å) CuCr2−xSbxS4 (x = 0, 0.25, and 0.5) Models

Table 1. Electronic Transport Parameters Obtained by Hall Effect Measurement of CuCr1.5Sb0.5S4 measured Hall coefficient (cm3/C) CuCr1.5Sb0.5S4

0.0025

carrier concentration p (cm−3) 2.41 × 10

21

Hall mobility (cm2/V·s) 4.06 2993

x

0

0.25

0.5

a b=c Cu−S Cr−S Sb−S

9.995 7.051 2.296 2.387

10.012 7.103 2.277−2.347 2.349−2.466 2.545−2.552

10.171 7.185 2.325−2.345 2.364−2.467 2.570−2.597

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Article

Chemistry of Materials ray analysis; this overestimation falls into the range of 1−3% that is often computed with the PBE functional. Substitutions of Cr with Sb lower the symmetry of the crystal structure. In the case of CuCr1.75Sb0.25S4, some Cu−S and Cr−S contacts are slightly shorter than those in the Sb-free phase, whereas some others reach 2.347 and 2.466 Å, respectively. A similar effect is computed for CuCr1.5Sb0.5S4. These distortions in the tetrahedral and octahedral environments of Cu and Cr, respectively, should decrease the thermal conductivity of the materials. Spin-polarized total and atom-projected density of states (DOS) of CuCr2−xSbxS4 (x = 0, 0.25, and 0.5) are sketched in Figure 6. Only ferromagnetic states were computed. DOS of Sb-free phase is similar to the one computed by Gupta et al.27 CuCr2S4 is computed to be metallic and magnetic; the cell moment is about 5.1 μB per formula unit. As shown in Figure 6, substituting Cr with Sb decreases the shift between spin-up and spin-down components of the DOS. This leads to the reduction of the magnetic moment (down to 4 μB for CuCr1.5Sb0.5S4) and qualitatively agrees with the magnetic data (Supporting Information, Figure S1). According to the DOS sketched in Figure 6, a semimetallic behavior is expected for these phases at very low temperature. Moreover, CuCr2−xSbxS4 (x = 0.25 and 0.5) should exhibit electrical conductivity lower than that of the Sb-free phase. Even if CuCr2−xSbxS4 (x = 0.25 and 0.5) phases exhibit semiconducting properties, their resistivity is very low, close to the one measured for semimetals. Moreover, one must keep in mind that our calculations are performed at zero temperature. Some charge-ordering or other charge reorganization among the Cu and Cr ions may occur at higher temperatures, leading to the creation of a semiconducting gap. As the content of Sb increases in the structure, electron filling of the valence band increases. This explains why the Seebeck coefficient increases with the Sb content at low temperature. Heremans et al. pointed out that large Seebeck coefficients can be obtained when the Fermi level is positioned where the DOS has a sharp slope.42 The comparison of the DOS sketched in Figure 6 shows that the DOS is higher at the Fermi level for the CuCr1.75Sb0.25S4 model. Moreover, among the three models, the DOS has the sharpest slope at the Fermi level for this CuCr1.75Sb0.25S4 model. The band structure of this model (Supporting Information, Figure S3), shows the presence of flat bands in the vicinity of the Fermi level. Although exact values cannot be determined from our calculations, these features are consistent with the large effective mass indicated from the transport measurements. Another important factor is the thermal conductivity. Apparently, it has a mixed trend and does not change linearly with Sb content (Figure 5c). Lattice thermal conductivity (Figure 5d) overall also shows a mixed trend. Although, while it is hard to explain this trend, there is a clue in bond lengths (Supporting Information, Table S1). (Cr/Sb)−S distances remain the same for Sb content of 0.22 and 0.3, resulting in similar thermal conductivity, but increase considerably in the case of Sb = 0.5. This increase in bond length seems to relax the structure and results in an increase in lattice thermal conductivity. On the other hand, Cu−S bond distances have a steady growth with increasing Sb content and do not seem to affect the thermal conductivity significantly. Overall, as mentioned above, the spinel structure with primarily octahedral coordination and partial filling of the polyhedra have a trend of lower lattice thermal conductivity34 with further disorder in our case coming from Cr/Sb mixing. Moreover, our effort to create vacancies on the sulfur site through sulfur deficient synthesis

Figure 6. Spin-polarized DFT DOS of CuCr2−xSbxS4 (x = 0, 0.25, and 0.5). No band gap at the Fermi level is observed, indicating that all the compositions are metallic at 0 K.

failed and ended up in multiphase samples. The lattice thermal conductivity shows similar temperature dependence to glassy/ disordered systems,43,44 namely little dependence or an increase with temperature at relatively high temperatures. In some systems, lone pairs have been proposed as a powerful mechanism to engender very low thermal conductivity.45−47 However, the thermal conductivity increasing largely for the heavy Sb doping rules out any dominant effect of the lone pairs in this case. Furthermore, the indication of the +5 valence state for Sb obtained from theoretical calculation rules out the existence of lone pairs. The minimum lattice thermal 2994

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Chemistry of Materials



ACKNOWLEDGMENTS Support from CREST, JST is acknowledged. We also thank the CNRS PICS Synthesis and Characterization of New NonOxides Functional Materials (Grant 2012 N°159769) for support of this work.

conductivity (κmin) uses the shortest scattering distance (l) of λ/2, where λ is the phonon wavelength. κmin can be approximated at high temperatures by the following equation:48 1 ⎛⎜ π ⎞⎟ 2⎝6⎠

1/3

κ min =

⎛ V ⎞−2/3 KB⎜ ⎟ (νs) ⎝n⎠



where KB is Boltzmann constant, V is the unit cell volume, n is the number of atoms per unit cell, and νs (2632 ms−1) is the average speed of sound.49 Using this formula, κmin is estimated at ∼0.2 W/mK. For all samples, the κlat values are higher than the κmin value, indicating that further reductions to κlat may be possible. Such a lowering might be achieved, for example, by reducing the grain size to enhance the boundary scattering of heat-carrying phonons.50,51 On the basis of all of this data, the dimensionless figure of merit of this series was calculated. It is obvious that the ZT increases with increasing Sb content up to a certain optimum amount of Sb (x = 0.30) and then starts decreasing again (Figure 5f). The highest value for ZT was found to be 0.43 at ∼643 °C, which is the best value among the many other related compounds of this series.22 This value is much higher than that of the other sulfur-based chalcopyrites, e.g., CuFeS2.11,13 It is slightly higher than TiS2, comparable to Cu0.1TiS2,12 but lower than the tetrahedrites.14

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.chemmater.6b05344. Measured magnetic properties and related text, X-ray diffractograms showing the miscibility gap, spin-up band structure calculations, and a table containing bond lengths (PDF)



REFERENCES

(1) He, J.; Kanatzidis, M. G.; Dravid, V. P. High performance bulk thermoelectrics via a panoscopic approach. Mater. Today 2013, 16, 166−176. (2) Koumoto, K.; Mori, T. Thermoelectric Nanomaterials; Springer: Heidelberg, 2013. (3) Slack, G. A. CRC Handbook of Thermoelectrics; Rowe, D. M., Ed.; CRC: New York, 1995; Chapter 34. (4) Mori, T.; Nolas, G. S. Contemporary innovations for thermoelectrics research and development. Scr. Mater. 2016, 111, 1−2. (5) Minnich, A. J.; Dresselhaus, M. S.; Ren, Z. F.; Chen, G. Bulk nanostructured thermoelectric materials: current research and future prospects. Energy Environ. Sci. 2009, 2, 466−479. (6) Li, J.-F.; Liu, W.-S.; Zhao, L.-D.; Zhou, M. High-performance nanostructured thermoelectric materials. NPG Asia Mater. 2010, 2, 152−158. (7) Gaultois, M. W.; Sparks, T. D.; Borg, C. K. H.; Seshadri, R.; Bonificio, W. D.; Clarke, D. R. Data-driven review of thermoelectric materials: performance and resource considerations. Chem. Mater. 2013, 25, 2911−2920. (8) Southern Fried Science. Alberta, Canada is the proud owner of the largest man-made pyramid on the planet. http://www. southernfriedscience.com/?p=13874 (accessed May 15, 2016). (9) Wan, C.; Gu, X.; Dang, F.; Itoh, T.; Wang, Y.; Sasaki, H.; Kondo, M.; Koga, K.; Yabuki, K.; Snyder, G. J.; Yang, R.; Koumoto, K. Flexible n-type thermoelectric materials by organic intercalation of layered transition metal dichalcogenide TiS2. Nat. Mater. 2015, 14, 622−627. (10) Guilmeau, E.; Bréard, Y.; Maignan, A. Transport and thermoelectric properties in copper intercalated TiS2 chalcogenide. Appl. Phys. Lett. 2011, 99, 052107. (11) Tsujii, N.; Mori, T. High thermoelectric power factor in a carrier-doped magnetic semiconductor CuFeS2. Appl. Phys. Express 2013, 6, 043001. (12) Maignan, A.; Guilmeau, E.; Gascoin, F.; Bréard, Y.; Hardy, V. Revisiting some chalcogenides for thermoelectricity. Sci. Technol. Adv. Mater. 2012, 13, 053003. (13) Ang, R.; Khan, A. U.; Tsujii, N.; Takai, K.; Nakamura, R.; Mori, T. Thermoelectricity generation and electron−magnon scattering in a natural chalcopyrite mineral from a deep-sea hydrothermal vent. Angew. Chem., Int. Ed. 2015, 54, 12909−12913. (14) Lu, X.; Morelli, D. T.; Xia, Y.; Zhou, F.; Ozolins, V.; Chi, H.; Zhou, X.; Uher, C. High performance thermoelectricity in earthabundant compounds based on natural mineral tetrahedrites. Adv. Energy Mater. 2013, 3, 342−348. (15) Reznitsky, L. Z.; Sklyarov, E. V.; Ushchapovskaya, Z. F.; Suvorova, L. F.; Polekhovsky, Y. S.; Dzerzanovsky, P.; Barash, I. G. Cuprokalininite, CuCr2S4, a new sulfospinel from metamorphic rocks of the sludyanka complex, South Baikal region. Geol. Ore Deposits 2011, 53, 758−766. (16) Reznitsky, L. Z.; Sklyarov, E. V.; Piskunova, L. F.; Ushchapovskaya, Z. F. Florensovite Cu(Cr1.5Sb 0.5)S4, a New Sulfospinel from the Baikal Region. Zap. Vsesoyuzn. Mineral. O-va 1989, 118, 57−65. (17) Warczewski, J.; Kusz, J.; Filimonov, D. S.; Kessler, Y. A.; Koroleva, L. I.; Mikheev, M. G.; Odintsov, A. G.; Aminov, T. G.; Busheva, E. V.; Shabunina, G. G. New antiferromagnetic semiconductor CuCr1.5Sb0.5S4. J. Magn. Magn. Mater. 1997, 175, 299−303. (18) Lotgering, F. K. Ferromagnetism in spinels: CuCr2S4 and CuCr2Se4. Solid State Commun. 1964, 2, 55−56. (19) Koroleva, L. I.; Demin, R. V.; Warczewski, J.; Krok-Kowalski, J.; Mydlarz, T.; Gilewski, A.; Pacyna, A. Giant Magnetoresistance in Spinel CuCr1.6Sb0.4S4. Phys. Stat. Sol. (a) 2002, 189, 853−857.

4. CONCLUSIONS The CuCr2−xSbxS4 (0 ≤ x ≤ 0.5) series was prepared and characterized. Lattice parameters vs Sb content are reported. A miscibility gap close to the Sb content of 0.15 is observed. Replacement of Cr by Sb lowers the thermal conductivity. We determined that the addition of Sb helps in enhancing the Seebeck coefficient. However, addition of Sb degrades the electrical conductivity. Overall, it helps in improving the power factor, and a high power factor at an Sb content of x = 0.3 is observed, resulting in an attractive ZT of 0.43. A MgAl2O4-type structure is quite common, and many different compounds comprised of this structure type have been reported. This report may further trigger intense research on the thermoelectric properties of this series of compounds because Sbdoped CuCr2S4 may be a serious candidate for consideration in medium temperature power generation where toxic and/or expensive element-free materials are required.



Article

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Takao Mori: 0000-0003-2682-1846 Notes

The authors declare no competing financial interest. 2995

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996

Article

Chemistry of Materials (20) Lotgering, F. K.; vanStapele, R. P.; vanderSteen, G. H. A. M.; Van Wieringen, J. S. Magnetic properties, conductivity and ionic ordering in Fe1−xCuxCr2S4. J. Phys. Chem. Solids 1969, 30, 799−804. (21) Koroleva, L. I. Antiferron states of charge carriers in the ferromagnetic semiconductors CuCr2S4‑xSex (0.5 ≤ x ≤ 1.5) with Curie points above room temperature. Zh. Eksp. Teor. Fiz. 1994, 106, 280−296. (22) Snyder, G. J.; Caillat, T.; Fleurial, J.-P. Thermoelectric Properties of selenide spinels. MRS Proceedings 2000, 626, Z3.3. (23) 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: Condens. Matter Mater. Phys. 2009, 80, 125205. (24) 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. Kristallogr. - Cryst. Mater. 2005, 220, 567−570. (25) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (26) Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.; Sutton, A. P. Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study. Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 57, 1505−1509. (27) Ramasamy, K.; Sims, H.; Gupta, R. K.; Kumar, D.; Butler, W. H.; Gupta, A. CoxCu1‑xCr2S4 nanocrystals: Synthesis, magnetism, and band structure calculations. Chem. Mater. 2013, 25, 4003−4009. (28) Monkhorst, H. J.; Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B: Condens. Matter Mater. Phys. 1976, 13, 5188−5192. (29) Blaha, P.; Schwarz, K.; Madsen, G. K. H.; Kvasnicka, D.; Luitz, J. WIEN2K, An Augmented Plane Wave + Local Orbitals Program for Calculating Crystal Properties; Techn. Universität Wien: Austria, 2001. (30) Kesler, Y. A.; Koroleva, L. I.; Mikheeva, M. G.; Odintsov, A. G.; Filimonov, D. S. A new antiferromagnetic semiconductor CuCr1.5Sb0.5S4. Inorg. Mater. 1993, 29, 135−137. (31) von Philipsborn, H.; Treitinger, L. Part B: Spinels, Fe Oxides, and Fe-Me-O Compounds. Landolt-Börnstein-Group III Condensed Matter 1980, 12b, 409. (32) Zajdel, P.; Kisiel, A.; Warczewski, J.; Konior, J.; Koroleva, L. I.; Krok-Kowalski, J.; Gusin, P.; Burattini, E.; Cinque, G.; Grilli, A.; Demin, R. V. The influence of the concentration of Sb ions onto the local crystal and electronic structures of CuCr2−xSbxS4 (x = 0.3, 0.4, 0.5) studied by XANES and EXAFS measurements and LAPW numerical calculations. J. Alloys Compd. 2005, 401, 145−149. (33) Krok-Kowalski, J.; Warczewski, J.; Koroleva, L. I.; Krajewski, K.; Gusin, P.; Duda, H.; Zajdel, P.; Pacyna, A.; Mydlarz, T.; Matyjasik, S.; Demin, R. V. On the influence of Sb concentration on the magnetization and magnetoresistivity in the spinel compounds CuCr2−xSbxS4 (where x = 0.3, 0.4, 0.5). J. Alloys Compd. 2004, 377, 53−58. (34) Spitzer, D. P. Lattice thermal conductivity of semiconductors: a chemical bond approach. J. Phys. Chem. Solids 1970, 31, 19−40. (35) Wang, P. L.; Kolodiazhnyi, T.; Yao, J.; Mozharivskyj, Y. Decoupling the electrical conductivity and Seebeck coefficient in the RE2SbO2 compounds through local structural perturbations. J. Am. Chem. Soc. 2012, 134, 1426−1429. (36) Neophytou, N.; Zianni, X.; Kosina, H.; Frabboni, S.; Lorenzi, B.; Narducci, D. Simultaneous increase in electrical conductivity and Seebeck coefficient in highly boron-doped nanocrystalline Si. Nanotechnology 2013, 24, 205402. (37) Mori, T.; Nishimura, T. Thermoelectric properties of homologous p- and n-type boron-rich borides. J. Solid State Chem. 2006, 179, 2908−2915. (38) Mori, T. High temperature thermoelectric properties of B12 icosahedral cluster-containing rare earth boride crystals. J. Appl. Phys. 2005, 97, 093703. (39) Hossain, M. A.; Tanaka, I.; Tanaka, T.; Khan, A. U.; Mori, T. Crystal growth and anisotropy of high temperature thermoelectric properties of yttrium borosilicide single crystals. J. Solid State Chem. 2016, 233, 1−7.

(40) Golikova, O. A.; Amandzhanov, N.; Kazanin, M. M.; Klimashin, G. M.; Kutasov, V. V. Electrical Activity of Impurities in βRhombohedral Boron. Phys. Stat. Sol. (a) 1990, 121, 579−586. (41) Mori, T.; Tanaka, T. Effect of transition metal doping and carbon doping on thermoelectric properties of YB66 single crystals. J. Solid State Chem. 2006, 179, 2889−2894. (42) Heremans, J. P.; Wiendlocha, B.; Chamoire, A. M. Resonant levels in bulk thermoelectric semiconductors. Energy Environ. Sci. 2012, 5, 5510−5530. (43) Slack, G. A. The thermal conductivity of nonmetallic solids. Solid State Physics; Seitz, F., Turnbull, D., Ehrenreich, H., Eds.; Academic Press: New York, 1979; Vol. 34. (44) Mori, T.; Martin, J.; Nolas, G. Thermal Conductivity of YbB44Si2. J. Appl. Phys. 2007, 102, 073510. (45) Skoug, E. J.; Morelli, D. T. Role of Lone-Pair Electrons in Producing Minimum Thermal Conductivity in Nitrogen-Group Chalcogenide Compounds. Phys. Rev. Lett. 2011, 107, 235901. (46) Nielsen, M. D.; Ozolins, V.; Heremans, J. P. Lone pair electrons minimize lattice thermal conductivity. Energy Environ. Sci. 2013, 6, 570−578. (47) Du, B.; Chen, K.; Yan, H.; Reece, M. J. Efficacy of lone-pair electrons to engender ultralow thermal conductivity. Scr. Mater. 2016, 111, 49−53. (48) Cahill, D. G.; Watson, S. K.; Pohl, R. O. Lower Limit to the Thermal Conductivity of Disordered Crystals. Phys. Rev. B: Condens. Matter Mater. Phys. 1992, 46, 6131−6140. (49) Kanomata, T.; Ido, H. Exchange Strictions of Chalcogen Spinels CuCr2X4 (X = S, Se and Te). J. Phys. Soc. Jpn. 1974, 36, 1322−1324. (50) Al Rahal Al Orabi, R.; Mecholsky, N. A.; Hwang, J.; Kim, W.; Rhyee, J.-S.; Wee, D.; Fornari, M. Band Degeneracy, Low Thermal Conductivity, and High Thermoelectric Figure of Merit in SnTe− CaTe Alloys. Chem. Mater. 2016, 28, 376−384. (51) Al Rahal Al Orabi, R.; Hwang; Lin, C.-C.; Gautier, R.; Fontaine, B.; Kim, W. J.; 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.



NOTE ADDED AFTER ASAP PUBLICATION This paper was published ASAP on February 28, 2017, with an inclusion of an equation in the Results and Discussion Section that should have been deleted. The corrected version was reposted on March 7, 2017.

2996

DOI: 10.1021/acs.chemmater.6b05344 Chem. Mater. 2017, 29, 2988−2996