RSC Advances Check for updates Cite this: RSC Adv., 2019, 9, 21451 Received 29th May 2019 Accepted 1st July 2019 DOI: 10.1039/c9ra04056k rsc.li/rsc-advances Study of thermal stability of p-type skutterudites DD0.7Fe3CoSbi2 by Knudsen effusion mass spectrometry František Zelenka,ab Pavel Brož,©*ab Jan Vřeštál,ab Jiří Buršík,©c Gerda Rogl,©de Andriy Grytsiv©de and Peter Rogl©de The temperature and phase stability of p-type skutterudites, DDo.7Fe3CoSbi2, manufactured via various preparation techniques, all exhibiting a high ZT-level, have been studied by means of thermal analysis and Knudsen effusion mass spectrometry. The results from phase transformation measurements and characteristics of the evaporation of antimony, as the volatile element, supported by microstructure observations and by diffusion profiles are summarized and discussed in view of a full understanding of the degradation processes and knowledge of the long term operation stability of the bulk and nano-structured thermoelectrics studied. It was found out that the antimony evaporation is a complex diffusion kinetic process resulting in a stable Sb level dependent on the preparation route. The studied p-type skutterudites, DD0.7Fe3CoSbi2, have proven their long term stability in thermoelectric devices at a maximum operation temperature of 600 °C. Complementary data on the structural, physical and mechanical properties of the materials are presented as well. 1 Introduction Thermoelectric materials have gained increased interest for their ability to motionlessly exploit waste heat (renewable energy) and convert it to electricity. The research in the field of new prospective bulk and nano-structured thermoelectrics is, therefore, of great importance for materials engineering in the automotive and aerospace industry, electronics and other areas. Improvement of thermoelectric materials and production processes concerns their thermoelectric properties, their thermodynamic and phase equilibria and also their long term thermal stability, which strongly depends on the content of volatile constituents as the current advanced materials contain elements like Cd, Zn, Sb, Sr and others. These elements can evaporate under operational conditions, leading to structure changes and to a negative influence of the thermoelectric properties. In this work thermal and phase stability of bulk and nanostructured p-type DD0.7Fe3CoSb12 skutterudites prepared by various procedures including ball milling, hot and cold "Masaryk University, Faculty of Science, Department of Chemistry, Kotlářská 2, 61137, Brno, Czech Republic. E-mail: broz@chemi.muni.cz bMasaryk University, Central European Institute of Technology, CEITEC, Kamenice 753/5, 62500, Brno, Czech Republic institute of Physics of Materials, Czech Academy of Sciences, Žižkova 22, 61662, Brno, Czech Republic dInstitute of Materials Chemistry, University of Vienna, Wahringerstrafie 42, A-1090 Vienna, Austria "Christian Doppler Laboratory for Thermoelectricity, Vienna, Austria pressing and high pressure torsion were studied by a combination of differential thermal analysis (DTA) and Knudsen effusion mass spectrometry (KEMS). The measurements of the evaporation characteristics of antimony by the KEMS method, which comprises the main part of the work, were complemented by measurement of the diffusion profiles. The data were subsequently used for evaluation of the diffusion kinetics of the Sb evaporation and for determination of the long term thermal stability of the materials. 2 Skutterudites Skutterudites present us with tantalizing potential applications in thermoelectrics because of their unique structure.15 Skutterudites are compounds with a body centered cubic crystal structure, space group Im3,6 and the simplest compounds i.e. unfilled skutterudites, have the formula TX3 (T is a transition element of the VIII group, e.g. Fe, Co, Ni... and X is a pnictogen atom, e.g. Sb). Electropositive elements E (E is a rare earth, an actinoid, an alkaline-earth, or an alkaline metal) and electronegative elements such as iodine and others7 may fill the icosahedral voids, formed by tilted octahedra, in the so called filled skutterudites with the chemical formula E^X^, where y is the filling fraction. Due to the unusually low energy vibrational modes of the filler atoms, called 'rattlers', in the oversized cages, the physics of filled skutterudites is essentially governed by the interplay of filler ions and their host structure.814 The efficiency of advanced thermoelectric technology is determined by the dimensionless figure of merit, ZT = S2T/pX This journal is © The Royal Society of Chemistry 2019 RSC Adv., 2019, 9, 21451-21459 | 21451 RSC Advances Paper where S is the Seebeck coefficient, which is proportional to a voltage A.V generated by a thermal gradient AT (AV/AT), T is the absolute temperature, p is the electrical resistivity and A the thermal conductivity, the latter being the sum of the electronic conductivity, Ae, and the lattice thermal conductivity, Aph (A = Ae + Aph). Reducing Aph is an important issue to achieve a high ZT. In this work, didymium, DD, consisting of Pr (4.76 mass%) and Nd (95.24 mass%), was used as natural double filler. Pr and Nd have a relative large atomic mass in respect to a small ionic radius providing conditions for rattling of these atoms in the icosahedral voids thereby reducing significantly the thermal conductivity via enhanced phonon scattering. The introduction of structural disorder via T-metal substitution (Fe/Co) further reduces Aph (via so-called alloy scattering of phonons) and consequently enhances ZT.15'*" Ball milling (BM) or high energy ball milling (HBM) is employed to produce small particle sizes in order to additionally increase the thermoelectric performance via nano-structuring.29 The "top down methods" like mechanical milling (BM, HBM) followed by hot pressing (HP) are reasonable ways to prepare materials with a crystallite size of the order of several 100 nm and crystal imperfections (e.g. dislocation densities of about 1012 irT2). With these methods nano-structured multifilled p- and n-type skutterudites have reached remarkable high ZT values (ZT > 1.4). Furthermore Severe Plastic Deformation (SPD) with High Pressure Torsion (HPT), as one of its major representatives, can achieve materials with ultra fine grains in the sub-micrometer or even nanometer range. With this method an extremely large plastic strain under a significantly high hydrostatic pressure is imposed on a sample without any essential change of the sample's geometry, which leads to (i) ultra fine grained structures and (ii) an enhanced formation of lattice defects -vacancies and particularly of dislocations (for details see ref. 28, 30-32). After HPT the lattice defects and grain boundaries reduce the lattice thermal conductivity substantially via phonon scattering and thus were found to enhance ZT at least by 20%.33 3 Experimental techniques 3.1. Preparation of p-type skutterudites TfDyFe3CoSb12 p-type skutterudites with the nominal composition DD07Fe3-CoSb12 were prepared from a master alloy Fe3CoSb12 via an optimized melting reaction technique from stoichiometric amounts of high purity elements (Fe, 99.5%, wire, Co, 99.9%, powder < 150 micron, Sb, 99.95%, crystals) by mixing, sealing into evacuated quartz tubes, melting at 950 °C followed by air quenching (details are given in ref. 22, 24, 25 and 29). Then the stoichiometric amount of DD (from Treibacher Industries AG, 99.9%) was added. The samples were then sealed under vacuum into quartz tubes, quickly heated to 600 °C, then slowly heated to 720 °C (2 °C min-1) and after a fast rise of the temperature (5 °C min-1), melted at 950 °C and air quenched. The reguli were ground in the glove box to obtain particles < 50 (xm, followed by BM (for details of ball-milling conditions see ref. 28) in Ar-filled tungsten carbide vessels (Fritsch planetary mill Pul-verisette 4) and hot-pressed, HP, (FCT uniaxial hot-press HP W 200/250 2200-200-KS, Ar, 700 °C, 56 MPa, 30 min). After 21452 | RSC Adv., 2019, 9, 21451-21459 determination of the composition by EPMA this sample is labeled as DD0.7Fe31Co0.gSb12-BM-HP. Additionally HPT was performed with an equipment from W. Klement, Austria: the sample was processed with 4 GPa and one revolution at 400 °C (further on referred to as HPT). For all details of HPT processing of skutterudites the reader is referred to ref. 29-31. Accordingly this sample is labeled as DD0.7Fe3.1Co0.gSb12-BM-HP-HPT. Furthermore two p-type skutterudite samples were prepared from commercially produced powders from Treibacher Industries AG (TIAG). For one sample (DD07Fe3CoSb12-HP) the powder was hot pressed under the same conditions as described above without BM (labelled as DD0.7Fe3CoSb12-HP (TIAG)). For the other sample (labelled as DD0.7Fe3CoSb12-CP-HPT (TIAG)) the original powder was cold pressed (CP) followed by HPT processing at 400 °C under nitrogen (to avoid oxidation) applying 4 GPa and 1 revolution (for details see ref. 33). A surplus of 3 mass% Sb was added in order to compensate losses in the preparation. In the hot pressing step most of this antimony was squeezed out again and only minor amounts remained at the grain boundaries. 3.2. Study of thermal and phase stability The study of thermal and phase stability of the materials prepared was performed by means of DTA and KEMS on a Netzsch STA 409 CD/3/403/5/G apparatus, a specially-adapted type of the commercial STA 409 CD-QMS 403/5 Skimmer Coupling Instrument.34'35 3.2.1 DTA measurements. Three heating and cooling cycles at a heating and cooling rate of 5 °C min-1 were performed during the DTA to exactly specify phase transformations in the material. The alloy DD0.7Fe31Co0.gSb12-BM-HP (prepared by common techniques and available in relatively sufficient amount from all the alloys) was chosen for the measurement. In order to prevent Sb evaporation and its potential negative influence on the measured data, the samples were sealed in evacuated quartz ampoules having the same shape like the original alumina crucibles, as discussed in ref. 35. Small pieces with the total weight of approximately 0.5 g were used for every measurement. Argon (purity of 99.9999) at the flow rate of 70 ml min-1 was used as carrier gas. 3.2.2 KEMS measurements. Evaporation characteristics of Sb via KEMS measurements for the alloy DD0.7Fe31Co0.gSb12-BM-HP were investigated in the temperature interval of 530-625 °C; at first for shorter times to gain preliminary information on the stability of the material. Detailed measurements of all alloys at temperatures 530, 560 and 590 °C were consequently performed until limiting values of ion current intensities and hence limiting values of vapour pressure of Sb were achieved, in order to exactly specify the long time thermal stability. The samples DD0.7Fe3.1Co0.9Sb12-BM-HP and DD0.7Fe3CoSb12-HP (TIAG) were studied at the first two temperatures only due to their limited amount. Samples in form of tablets having approximately 400 mg for the preliminary measurement and 100-250 mg for the long time measurements were placed in the Knudsen cell and measured. Measurement of pure Sb was performed after the measurement of each sample and the data This journal is © The Royal Society of Chemistry 2019 Paper RSC Advances were used for evaluation of the evaporation processes of Sb from the investigated material. Pieces of Sb having approximately 400 mg were used for the measurement. Monitoring of monomer ions of the Sb isotopes generated at ionization energy of 35 eV was preferred due to their largest response from all measured oligomers, as found during the instrument calibration by measurement of the heats of sublimation of selected elements34 and as also used in our recent study of the thermal stability of primary CoSb3 skutterudite.36 3.3. Measurement of diffusion profiles Depth profiles of element distribution were measured on samples after thermal exposure at 530, 560 and 590 °C using a Tescan LYRA 3XMU SEMxFIB scanning electron microscope (SEM) working at acceleration voltage of 20 kV and equipped with an X-Max 80 Oxford Instruments detector for energy dispersive X-ray analysis (EDX). Bulk samples were mechanically polished to remove several tens of micrometres of the surface layer and final stages of polishing were done very gently using a Buehler VibroMet vibratory polisher to produce a sample edge as sharp as possible for near-surface profile measurements. Regarding the sample edge curvature and the interaction volume of EDX analyses, the depth of the region influenced could be measured with an uncertainty of about 1 \im. 3.4. Determination of the structure; physical and mechanical properties measurements Details of measurements (X-ray intensity data, density, crystallite size, electrical resistivity, Seebeck coefficient, thermal conductivity, hardness) are summarized in our previous articles on the subject.2831 4 Results and discussion 4.1. Physical, structural and mechanical properties Whereas the crystallite size for DD0.7Fe3.1Co0.gSb12-BM-HP (153 nm) is reduced to 54 nm after HPT, the dislocation density is enhanced for the HPT processed sample DDo.7Fe31Coo.gSb12-BM-HP-HPT from 2.8 x 1013 irT2 to 2.6 x 1014 irT2. These changes are reflected in the physical properties (see below). The relative density (rfrei = 98.9%) was higher for DDo.7Fe31Co0.9-Sb12-BM-HP than for DDo.7Fe3.1COo.9Sb12-BM-HP-HPT (_drel = 97.5%), because of the introduction of fine micro cracks during HPT. For DDo.7Fe31COo.9Sb12-BM-HP a ZT = 1.2 at 725 K was reported, which is one of the highest values for a p-type skutterudite. The HPT processed sample, due to smaller grains and higher dislocation density, DDo^Fe^Coo.gSb^-BM-HP-HPT had a higher electrical resistivity, which was overcompensated by a much lower thermal conductivity, especially in the temperature range of about 500-600 K, when the electrical resistivity decreases from a maximum value. This behavior together with a slightly enhanced Seebeck coefficient raises ZT to ZT > 1.4 at 725 K (for details see ref. 28, 30 and 31). The lattice parameters of DD0.7Fe3 iCOo.gSb^-BM-HP and DDo.yFea.iCOo.gSb^-BM-HP-HPT were a = 0.91123(2) nm and a = 0.91128(3) nm, but after the measurement with the This journal is © The Royal Society of Chemistry 2019 Knudsen cell the lattice parameter of DDo.7Fe31Coo.gSb12 is slightly smaller: a = 0.91118(2) nm. The Hall-Petch-type relation applies for HPT processed skutterudites: the smaller the crystallite size, the harder the alloy. Regardless of the fact that the higher the density, the harder the material, DD0.7Fe31-Co0.9Sb12-BM-HP has HV 0.1 ~ 460, whereas DDo.7Fe31COo.g-Sb12-BM-HP-HPT reveals HV 0.1 > 500 in the center and HV 0.1 > 600 in the rim area, the part of the sample exposed to the highest shear strain during HPT processing (for details see ref. 28, 30 and 31). DD0.7Fe3CoSb12-HP and DD0.7Fe3CoSb12-CP-HPT follow the trend of DDo.7Fe3.1COo.9Sb12-BM-HP and DDo.7Fe3.1Co0.9Sb12-BM-HP-HPT with a smaller grain size and a higher dislocation density, a slightly bigger lattice parameter and a lower relative density for the HPT processed sample (see Table 1; data adopted from ref. 28, 30 and 31). For convenience to the reader Fig. 1 summarizes the ZT values as a function of temperature and mechanical properties (hardness) vs. density of the materials for which the Sb-evaporation has been studied in this work. 4.2. Differential thermal analysis Information on the phase stability of DD0.7Fe3CoSb12 and respective phase transitions is crucial in the process of preparation and investigation of the thermal stability of the material. DTA yields valuable information for the selection of a maximum temperature up to which the thermal stability can be investigated. The DTA measurement was done on the DD0.7Fe31-Co0.gSb12-BM-HP alloy as mentioned in Section 3. This alloy prepared by common procedures can be taken as the representative for all the investigated alloys. The first measurement was realized up to 1000 °C in order to gain a complete series of phase transformation temperatures (see Fig. 2a). The second measurement realized up to 700 °C was aimed at a precise identification of the ternary transformation close to the binary eutectic isotherm (liquid <-> CoSb3 + Sb, see also ref. 36) at about 620 °C, which is a limiting temperature factor for the thermal stability of the CoSb3 skutterudite based thermoelectric materials (see ref. 36). The respective DTA record is shown in Fig. 2b. Samples, before and after the measurements, are shown in Fig. 3. Whilst the complete series of phase transformations led to sample decomposition, it is not the case of the eutectic reaction only. The temperatures of the phase transitions, as identified for the DD0.7Fe3CoSb12 skutterudite, are somewhat different in comparison with those in the binary Co-Sb36'37 and especially in the ternary Fe-Co-Sb38 system. As for the Fe-Co-Sb Table 1 Lattice parameter, a, in nm; crystallite size, cs, in nm; dislocation density, did, in 1014 irT2; relative density, drel, in %; ZT. Data adopted from ref. 28, 30 and 31 Compound a cs did dre\ ZT DD0.7Fe3.1Co0.9Sb12-BM-HP 0.91123(2) 152 0.28 98.7 1.2 DDo.yFes.iCOo.gSb^-BM-HP-HPT 0.91128(3) 53 2.6 98.1 1.4 DD07Fe3CoSb12-HP (TIAG) 0.91086(4) 76 0.9 99.7 1.1 DDo.7Fe3CoSb12-CP-HPT (TIAG) 0.91126(3) 44 3.6 98.0 1.3 RSC Adv., 2019, 9, 21451-21459 | 21453 RSC Advances Paper N a) quartz ampoule b) quartz ampoule Fig. 1 Summary of figure of merit data ZT vs. temperature as well as Vickers hardness vs. density for the compounds with various treatments studied in this work. Fig. 2 DTA measurement of the alloy DD07Fe3 jCoo^Sb^-BM-HP to (a) 1000 °C and (b) 700 °C at a heating and cooling rate of 5 °C min-1 and Ar flow rate of 70 ml min-1 First heating (green line), second heating (blue line) and third heating (brown line) (for interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article). system, one can notice that the Co1_xFexSb3 skutterudite is only stable for Fe concentrations lower than x ~ 0.26 at 580 °C,38 whereas the DD-flller stabilizes the DD0 7Fe3CoSb12 skutterudite (and even DDFe4Sb12). The phase stability of the DD07Fe3-CoSb12 skutterudite is mainly governed by Nd, which before after Fig. 3 Quartz ampoule with sample DD0.7Fe3 iCoo.9Sbi2-BM-HP before and after DTA measurement to (a) 700 and (b) 1000 °C. constitutes 95.24 mass% of DD. The stabilization effect of Nd for a ternary skutterudite NdFe4Sb12 is already obvious in the Nd-Fe-Sb system.39 41 However, it should be noted that although single crystal data have shown various degrees of defects in the Nd sublattice (Nd0 832Fe4Sb12,42 Nd0.73Fe4Sb12 (ref. 43)), phase diagram studies usually display full stoichiometry for Nd1Fe4Sb12.39'41 Data from the thermal analysis gained in this work can be used for thermodynamic and phase equilibria studies of quaternary Nd-Fe-Co-Sb or quinary Pr-Nd-Fe-Co-Sb systems, however, the aim of this paper is not a study of phase transformations and phase equilibria in these higher order systems. According to this information and on the basis of our DTA data, the compound DD07Fe3CoSb12 should decompose between 720 to 780 °C. This is also consistent with (a) phase diagram calculations by Wang et al.39 indeed revealing a phase decomposition of NdFe4Sb12 at 777 °C, (b) experimental data of Zeng et al.41 and Evers et al." who observed the skutterudite at 500 °C and 600 °C, respectively, as well as (c) experimental data of Nasir et al. ,40 which document the nonexistence of Nd1_xFe4Sb12 at 800 °C. The heat effects at high temperatures visible on heating and cooling curves of the DTA record (Fig. 2a) include the decomposition of the DD07Fe3-CoSb12 skutterudite and phase transformations with other intermetallic phases, originating during cooling, and with liquid phase (see the signals between 800 and 900 °C on the cooling curve). 4.3. Knudsen effusion mass spectrometry The preliminary effusion measurement, performed on the alloy DD0 7Fe3 iCo0 9Sb12-BM-HP in a wide temperature range as a sequence of vapour pressure measurements at various increasing temperatures, showed different sublimation behavior below and above the eutectic temperature of the system, visible as a sudden increase of the ion current intensity above the phase transition temperature about 610 °C on Fig. 4. In comparison with the vapour pressure measurement of pure Sb as the reference material, the results of which are given in ref. 34, performed after the measurement of skutterudite, no obvious linear temperature dependence of the logarithm of the partial pressure of Sb was found. Besides, the measurement of the skutterudite revealed an exponential decrease of the partial pressure of Sb at constant temperatures (see Fig. 4), and as 21454 I RSC Adv., 2019, 9, 21451-21459 This journal is © The Royal Society of Chemistry 2019 Paper RSC Advances 4.7E-11 4.2E-11 3.7E-11 < 3.2E-11 t 2.7E-11 g 2.2E-11 1.7E-11 1.2E-11 7.0E-12 -Mass 121 -Mass 123 -Temperat ure \ i. \ j Es S n r w r 1 1 0 >, L t v. 1C i V * -J k. 640 620 600 580 560 540 520 500 50 150 250 Time [min] 350 450 Fig. 4 Vapour pressure measurement of Sb at increasing temperature steps (upper line) in the temperature interval of 530-625 °C for the sample DDo.7Fe31Coo.9Sb12-BM-HP. Ion current intensities of monomers 121Sb and 123Sb (in the order of decreasing intensity) acquired at an electron energy of 35 eV (background values create the minima). Multiple abrupt drops of ion current correspond to closing the shutter between the Knudsen cell and the ion source, which was done to check the stability of the measured signals and to reduce potential contamination of the ion source by Sb. expected, falling finally down to limiting values of the partial pressure in infinite time, which indicates a complex diffusion kinetic process and the structure stability of the skutterudite. The exponential decrease of the partial pressure of Sb shows an easy and fast initial release of free or weakly bound Sb from the material, most likely the small amount of surplus free Sb appearing during the preparation of the material, as mentioned in Section 3.1 sample preparation. The limiting values can be assigned to the partial pressure of Sb from the skutterudite DD0.7Fe3CoSb12. For a theoretical description of the exponential behavior, the diffusion equation (eqn (1)) was used y = a\é l(-a2(x-a3))+a4 (1) where y and x introduce the ion current intensity and time and ax-aA are fitting parameters. The equation yielded similar limiting values of the vapour pressure of Sb of about 0.5 Pa. These values were taken as preliminary information for subsequent measurements at selected constant temperatures. The subsequent long time measurements led to experimentally determined limiting values of the vapour pressure of Sb (see Fig. 5a-c as an example for alloy DD0.7Fe3.1Co0.gSb12-BM-HP-HPT and Table 2), which were used for the calculation of mass depletion of samples by Sb evaporation. For this purpose the Hertz-Knudsen equation (eqn (2)) was applied, d«s PsbAoC at ^/2nRTMsb' (2) where nsb is the number of moles of Sb, t is the time in seconds, psb is the partial pressure of antimony, A0 is the area of the effusion orifice and C is its Clausing factor, R is the molar gas constant, T is the temperature in Kelvin and Msb is the molar mass of antimony. 1.1E-12 y 9.0E-13 ■M ratu c 0J 1_ 7.0E-13 CoSb3 + Sb. No linear temperature dependence of This journal is © The Royal Society of Chemistry 2019 RSC Adv., 2019, 9, 21451-21459 | 21457 RSC Advances Paper logarithm of partial pressure of Sb was found by KEMS as it was in case of pure Sb, measured as the reference material after every sample measurement. The partial pressure is exponentially reaching a limiting value, which can be assigned to the partial pressure of the skutterudite structure. The Sb evaporation is a complex diffusion kinetic process resulting in a stable state defined by the limiting value, which is dissimilar for the same skutterudite prepared by different procedures or by their combinations. Nevertheless, the differences were found to be small. KEMS is a suitable method for in situ monitoring of volatile components and shortens the time for reaching a stable limiting state in comparison with annealing performed at common operation pressure (1 bar) requiring very long time periods and not allowing direct measurement of the volatile species. Concentration of all elements in the studied samples after KEMS measurement is almost constant and X-ray diffraction pattern before and after the KEMS measurement is unchanged, confirming the structure stability of the skutterudite structure. The studied p-type skutterudites, DD0 7Fe3CoSb12, create confidence in their long term use in thermoelectric devices at a maximum operation temperature of 600 °C. Conflicts of interest There are no conflicts to declare. Acknowledgements The research reported herein was financially supported by the Czech Science Foundation under the project GA 17-12844S as well as by the project CEITEC 2020 (LQ1601) from the Ministry of Education, Youth and Sports of the Czech Republic under the National Sustainability Programme II. The support of OEAD for the joint bilateral project 2019-CZ/01 is gratefully acknowledged. The authors thank Treibacher Industries AG, TIAG, Austria (Dr Markus Hochenhofer) for supplying the skutterudite powder. References 1 D. M. Rowe, Handbook of Thermoelectrics, CRC Press, Boca Raton, FL, 2006. 2 J. S. Sakamoto, H. Schock, T. Caillat, J. P. Fleurial, R. Maloney, M. Lyle, T. Ruckle, E. Timm and L. Zhang, Skutterudite-based thermoelectric technology for waste heat recovery: progress towards a 1 kW generator, Set Adv. Mater., 2011, 3, 621-632. 3 D. M. Rowe, Modules, Systems, and Applications in Thermoelectrics, CRC Press, London, 2012. 4 X. Zhou, G. Wang, L. Guo, H. Chi, G. Wang, Q. Zhang, Ch. Chen, T. Thompson, J. Sakamoto, V. P. Dravid, G. Cao and C. Uher, Hierarchically structured Ti02 for Ba-filled skutterudite with enhanced thermoelectric performance, /. Mater. Chem. A, 2014, 2, 20629-20635. 5 X. Meng, Y. Liu, B. Cui, D. Qin, J. Cao, W. Liu, Z. Liu, W. Cai and J. Sui, High thermoelectric performance of single phase p-type cerium-filled skutterudites by dislocation engineering,/. Mater. Chem. A, 2018, 6, 20128-20137. 6 W. Jeitschko and D. B. Braun, LaFe4P12 filled with CoAs3-type structure and isotypic Ianthanoid transition metal polyphosphides, Acta Crystallogr., Sect. B: Struct. Crystallogr. Cryst. Chem., 1977, 33, 3401-3406. 7 B. Duan, J. Yang, J. R. Salvador, Y. He, B. Zhao, S. Wang, P. Wei, F. S. Ohuchi, W. Zhang, R. P. Hermann, O. Gourdon, S. X. Mao, Y. Cheng, Ch. Wang, J. Liu, P. Zhai, X. Tang, Q. Zhang and J. Yang, Electronegative guests in CoSb3, Energy Environ. Set, 2016, 9, 2090-2098. 8 L. Hicks and M. Dresselhaus, Thermoelectric figure of merit of a one-dimensional conductor, Phys. Rev. B: Condens. Matter Mater. Phys., 1993, 47(24), 16631-16634. 9 G. P. Meisner, D. T. Morelli, S. Hu, J. Yang and C. Uher, Structure and lattice thermal conductivity of fractionally filled skutterudites: Solid solutions of fully filled and unfilled end members, Phys. Rev. Lett, 1998, 80, 3551-3554. 10 H. Anno, K. Matsubara, Y. Notohara, T. Sakakibara and H. Tashiro, Effects of doping on the transport properties of CoSb3,/. Appl. Phys., 1999, 86, 3780-3786. 11 R. Venkatasubramanian, E. Siivola, T. Colpitts and B. O'Quinn, Thin-film thermoelectric devices with high room-temperature figures of merit, Nature, 2001, 413, 597-602. 12 L. D. Chen, T. Kawahara, X. F. Tang, T. Goto, T. Hirai, J. S. Dyck, W. Chen and C. Uher, Anomalous barium filling fraction and n-type thermoelectric performance of BayCo4Sb12,/. Appl. Phys., 2001, 90(4), 1864-1868. 13 G. Chen and A. Shakouri, Heat transfer in nanostructures for solid-state energy conversion, /. Heat Transfer, 2002, 124, 242-252. 14 G. S. Nolas, G. Fowler and J. Yang, Assessing the role of filler atoms on the thermal conductivity of filled skutterudites, /. Appl. Phys., 2006, 100(4), 43705. 15 B. C. Sales, D. Mandrus and R. K. Williams, Filled skutterudite antimonides: A new class of thermoelectric materials, Science, 1996, 22, 1325-1327. 16 B. C. Sales, D. Mandrus, B. C. Chakoumakos, V. Keppens and J. R. Thompson, Filled skutterudite antimonides: Electron crystals and phonon glasses, Phys. Rev. B: Condens. Matter Mater. Phys., 1997, 56, 15081. 17 S. Nolas, J. L. Cohn and G. A. Slack, Effect of partial void filling on the lattice thermal conductivity of skutterudites, Phys. Rev. B: Condens. Matter Mater. Phys., 1998, 58(1), 164-170. 18 V. L. Kuznetsov, L. A. Kutznetsova and D. M. Rowe, Effect of partial void filling on the transport properties of NdxCo4Sb12 skutterudites,/. Phys.: Condens. Matter, 2003, 15, 5035-5048. 19 X. Shi, H. Kong, C. P. Li, C. Uher, J. Yang, J. R. Salvador, H. Wang, L. Chen and W. Zhang, Low thermal conductivity and high thermoelectric figure of merit in n-type BaxYbyCo4Sb12 double-filled skutterudites, Appl. Phys. Lett., 2008, 92, 182101. 20 Y. Z. Pei, J. Yang, L. D. Chen, W. Zhang, J. R. Salvador and J. Yang, Improving thermoelectric performance of caged 21458 | RSC Adv., 2019, 9, 21451-21459 This journal is © The Royal Society of Chemistry 2019 Paper RSC Advances compounds through light-element filling, Appl. Phys. Lett., 2009, 95, 042101. 21 L. Zhang, A. Grytsiv, M. Kerber, P. RogI, E. Bauer and M. Zehetbauer, Thermoelectric performance of mischmetal skutterudites MmyFe4_xCoxSb12 at elevated temperatures, /. Alloys Compd., 2010, 490, 19-25. 22 G. RogI, A. Grytsiv, E. Bauer, P. RogI and M. Zehetbauer, Thermoelectric properties of novel skutterudites with didymium: DDy(Fe1_xCox)4Sb12 and DDy(Fe1_xNix)4Sb12, Intermetallics, 2010, 18, 57-60. 23 P. F. Qiu, J. Yang, R. H. Liu, X. Shi, X. Y. Huang, G. J. Snyder, W. Zhang and L. D. Chen, High-temperature electrical and thermal transport properties of fully filled skutterudites RFe4Sb12 (R: Ca, Sr, Ba, La, Ce, Pr, Nd, Eu, and Yb),J. Appl. Phys., 2011, 109, 063713. 24 G. RogI, A. Grytsiv, P. RogI, E. Bauer, M. B. Kerber, M. Zehetbauer and S. Puchegger, Multifilled Nanocrystalline p-type didymium - skutterudites with ZT>1, Intermetallics, 2010, 18, 2435-2444. 25 G. RogI, A. Grytsiv, P. RogI, E. Bauer and M. Zehetbauer, Enhanced thermoelectric figure of merit in p-type DD^Fe^ xCox)4Sb12, Solid State Phenom., 2011, 170, 240-243. 26 C. Zhou, L. Zhang and J. Sakamoto, Thermoelectric Properties of p-type skutterudite nanocomposites, in Nanoscale Thermoelectrics, lecture Notes in Nanoscale Science and Technology, ed. X. Wang and Z. M. Wang, Springer Int. Publishing, Switzerland, 2014, pp. 271-302. 27 G. RogI, A. Grytsiv, P. RogI, E. Bauer and M. Zehetbauer, A new generation of p-type didymium skutterudites with high ZT, Intermetallics, 2011, 19, 546-555. 28 G. RogI, A. Grytsiv, P. RogI, E. Bauer, M. Hochenhofer, R. Anbalagan, R. C. Mallik and E. Schafler, Nanostructuring of p- and n-type skutterudites reaching figures of merit approximately 1.3 and 1.6 respectively, Acta Mater., 2014, 76, 434-448. 29 L. Zhang, A. Grytsiv, M. Kerber, P. RogI, E. Bauer, M. J. Zehetbauer, J. Wosik and G. E. Nauer, MmFe4Sb12-and CoSb3-based nano-skutterudites prepared by ball milling: Kinetics of formation and transport properties, /. Alloys Compd., 2009, 481, 106-115. 30 G. RogI, P. RogI, E. Bauer and M. Zehetbauer, Severe Plastic Deformation, a tool to enhance thermoelectric performance, in Thermoelectric Nanomaterials, ed. K. Kuomoto and M. Mori, Springer-Verlag, Berlin Heidelberg, 2013, pp. 193-254. 31 G. RogI, A. Grytsiv, P. RogI, E. Royanian, E. Bauer, J. Horky, D. Setman, E. Schafler and M. Zehetbauer, Dependence of thermoelectric behaviour on severe plastic deformation parameters: a case study on p-type skutterudite DD0 60Fe3CoSb12, Acta Mater., 2013, 61, 6778-6789. 32 R. C. Mallik, R. Anbalagan, G. RogI, E. Royanian, P. Heinrich, E. Bauer, P. RogI and S. Suwas, Thermoelectric properties of Fe0.2Co3.8Sb12_xTexskutterudites,ActaMater., 2013, 61, 6698-6711. 33 G. RogI, A. Grytsiv, R. Anbalagan, J. Bursik, M. Kerber, E. Schafler, M. Zehetbauer, E. Bauer and P. RogI, Direct SPD-processing to achieve high-ZT skutterudites, Acta Mater., 2018, 159, 352-363. 34 P. Broz and F. Zelenka, Specially-adapted type of Netzsch STA Instrument as a tool for Knudsen effusion mass spectrometry, Int. J. Mass Spectrom., 2015, 383, 13-22. 35 P. Brož, F. Zelenka, J. Sopoušek, M. Hejdukova, T. Kuběnová and O. Zobač, Thermal analysis and Knudsen effusion mass spectrometry combined in a specially-adapted commercial skimmer coupled instrument (Netzsch), CALPHAD: Comput. Coupling Phase Diagrams Thermochem., 2019, 65, 86-92. 36 P. Brož, F. Zelenka, Z. Kohoutek, J. Vřeštál, V. Vykoukal, J. Buršík, A. Zemanová, G. RogI and P. Rogl, Study of thermal stability of CoSb3 skutterudite by Knudsen effusion mass spectrometry, CALPHAD: Comput. Coupling Phase Diagrams Thermochem., 2019, 65, 1-7. 37 Y. Zhang, Ch. Li, Z. Du and T. Geng, The thermodynamic assessment of the Co-Sb system, CALPHAD: Comput. Coupling Phase Diagrams Thermochem., 2008, 32, 56-63. 38 L. D. Dudkin and N. K. Abrikosov, Alloying of the semiconductive compound cobalt antimonide CoSb3 (in Russian), Fiz. Tverd. Tela, 1959, 1(1), 142-151. 39 W. Wang, T. Liu, H. Song, S. Xue, W. Zhang, H. Zhang, L. Zhang, Y. Jia, N. Wang, X. Cheng and Ch. Zeng, A Thermodynamic Modeling of the Fe-Nd-Sb System, Mater. Trans., 2016, 57(2), 103-111. 40 N. Nasir, A. Grytsiv, P. RogI, D. Kaczorowski and H. S. Effenberger, The system Nd-Fe-Sb: Phase equilibria, crystal structures and physical properties, Intermetallics, 2010, 18, 2361-2376. 41 L. Zeng, P. Qin, L. Nong, J. Zhang and J. Liao, The 773K isothermal section of the Nd-Fe-Sb ternary system, /. Alloys Compd., 2007, 437, 84-86. 42 C. B. H. Evers, W. Jeitschko, L. Boonk, D. J. Braun, T. Ebel and U. D. Scholz, Rare earth and uranium transition metal pnictides with LaFe4P12 structure,/. Alloys Compd., 1995, 224, 184-189. 43 E. Bauer, St. Berger, A. Galatanu, Ch. Paul, M. Delia Mea, H. Michor, G. Hilscher, A. Grytsiv, P. RogI, D. Kaczorowski, L. Keller, T. Hermannsdörfer and P. Fischer, Magnetic behaviour of PrFe4Sb12 and NdFe4Sb12 skutterudites, Phys. Rev. B: Condens. Matter Mater. Phys., 2002, 312-313, 840-842. 44 G. RogI, A. Grytsiv, E. Bauer and P. RogI, Thermoelectric Sb-based Skutterudites for medium temperatures - recent advances, in Advanced Thermoelectrics: Materials, Contacts, Devices and Systems, ed. Z. Ren, Y. Lan and Q. Zhang, CRC Press, Boca Raton, FL, 2018, pp. 193-230. This journal is © The Royal Society of Chemistry 2019 RSC Adv., 2019, 9, 21451-21459 | 21459