Interaction of Skutterudites with Contact Materials: A Metallurgical Analysis

More than hundred diffusion couples between p- and n-type skutterudites and various materials were prepared and interaction zones were investigated after heat treatment at 600 °C for 1100 h. The constitution of reaction/diffusion zones was discussed in terms of: (a) atom site preference in the skutterudite lattice, (b) phase equilibria in multicomponent systems and (c) particularities of the crystal structure of intermediate phases. It could be shown that phase composition and thermo-mechanical properties of bonding can be engineered by chemical substitution. The results obtained allowed the determination of several necessary criteria for the development of chemically and mechanically stable diffusion barriers/couples for skutterudite based thermoelectric (TE) modules.


Introduction
Thermoelectric (TE) energy conversion enables a direct transformation of heat into electricity. This technology has a great potential in various fields, like in waste heat recovery or in increasing the efficiency of industrial processes to not only reduce primary energy but also to minimize CO 2 emissions. Besides tellurides, half-Heuslers, silicides, Zintl phases, oxides and more, skutterudites are most promising high efficient TE materials for middletemperature applications up to 600-700°C. To our knowledge, the highest dimensionless figure of merit (ZT) values for these bulk materials are ZT = 1.4 for p-type DD y Fe 3 CoSb 1.2 [1] (DD = didymium, a natural mixture of Nd and Pr) and ZT = 1.9 for n-type (Sr,Ba,Yb) y Co 4 Sb 12 . [2] Skutterudites can be produced from relatively cheap and abundant starting materials, and commercial scale production for p-DD y Fe 3 CoSb 12 (ZT 800K = 1.2), p-MM y Fe 3-CoSb 12 (ZT 800K = 1.0), and n-(MM,Sm) y Co 4 Sb 12 (ZT 800K = 1.3) was established by Treibacher Industry AG (TIAG, Austria). With such high ZT values, and proven long term thermal stability [3] a thermoelectric power generation already attracts interest for practical recuperation of waste heat into electricity, and requests the development of interface materials for durable TE modules.
A thermoelectric module is a complex set of different materials, which needs to be engineered in a highly sophisticated manner in order to provide a long term flawless and constant performance with respect to thermoelectric efficiency, which in turn directly depends on a chemically and metallurgically unaltered materials' constitution. The materials have to withstand thermal stresses, created due to rapid temperature cycling, and temperature gradients across the module, therefore the mechanical properties of TE materials are equally important, [4,5] especially the coefficient of thermal expansion, as for skutterudites it is known that it differs for the p-and n-type alloys. [6] Various pure metals and binary intermetallics (Fig. 1), ternary (Fe-Ni-Cr [7] and Co-Fe-Ni [8] ), silicides Co-Si, [9,10] antimonides Ni-Sb, [11] Bi 2 Te 3 -based TE materials [12] and nitrides TaN [13,14] were already used as contact materials for skutterudite-based TE modules.
These numerous investigations are mainly focused on the production and test of the TE durability of the TE modules. Various skutterudite compositions, contact materials and conditions have been tested, but systematic investigations of chemical aspects of the interaction of different materials with skutterudites are missing. We believe that fundamental information is requested for a successful development of durable skutterudite-based TE modules, as well as for a better understanding of the chemistry of these materials.
The current work is devoted to systematic investigations of interactions of skutterudites with different materials at the same conditions. Several materials were tested in view of two aspects: (a) ability of the elements to enter the skutterudite lattice (formation of skutterudites-structural chemical aspect) and (b) interaction of the contact materials with skutterudites and formation of intermediate phases (thermodynamic and phase equilibria aspects).
In this article we will report results of the interaction of the unfilled skutterudite CoSb 3 with various pure metals. Then, the research was extended to single-filled n-Ba y Co 4 Sb 12 and n-Sr y Co 4 Sb 12 , p-Ce y Fe 3 CoSb 12 and to multi-filled p-MM y Fe 3 CoSb 12 (MM = mischmetal, a natural mixture of La, Ce, Pr and Nd). Finally, we report the test results of the interaction of binary intermetallics with filled skutterudites, and compare results with pure elements.
In a first step, the experiments on the interaction of pure elements and nitride TiN with the unfilled skutterudite CoSb 3 were performed by means of powder metallurgical techniques. Powders were sieved to grain sizes of 20-50 lm, and mixed with CoSb 3 in a molar ratio 1:2 and then they were cold-pressed and heat-treated in evacuated quartz ampoules at 600°C for 14 days. The product of interaction was investigated by means of x-ray powder diffraction (XPD), collected from a HUBER-Guinier image plate with monochromatic CuK a1 -radiation (k = 0.154056 nm).
Metals and intermetallics for the investigation of the diffusion couples were arc-melted on a water-cooled copper hearth in Ti-gettered argon from elemental ingots with a minimal purity better than 99 mass%. The diameter of the so obtained pseudo-spherical specimens ([ = 1.9-2.1 mm) was measured and these balls were embedded into the skutterudite matrixes via hot-pressing (56 MPa, 600°C, 2 h in 5 N Ar employing a uniaxial hot press   [36] B (Mo-Ti), [29] C (Ag-Pd), [34] D (Mo-Ti), [31,35] E (Al-Ti), [7,[37][38][39] F (Fe-Ni), [40] G (Al-Ni) [41] W200/250-2200-200-KS). The tablets with a diameter of 10 mm were vacuum-sealed in quartz tubes and heattreated at 600°C for 1100 h. The annealed cylinders were ground to the appearance of the embedded balls, and then polished using standard procedures. Microstructure and chemical compositions were analysed by scanning electron microscopy/backscattered electrons (SEM/BSE) and electron probe microanalysis (EPMA) via an INCA Penta FETx3-Zeiss SUPRA TM 55VP equipment with an EDX detector.

Formation and General Consideration on Chemical Activity of the Skutterudites
Skutterudites crystallise in the structure type of CoAs 3 or LaFe 4 P 12 (space group Im 3), as a filled variant. The chemical formula of these compounds can be described as E y T 4 X 12 , where T is a transition element of the VIII group in 8c (, , ); the electropositive element (E) occupies the icosahedral hole in 2a (0, 0, 0) which can be filled by a number of electropositive elements such as Ca, Sr, Ba, Tl, lanthanides and actinides, or under special conditions by Y, Hf, Pb and Sn (high pressure or multilayer precursor methods) and by K and Na (in nanowire state); X is a pnicogen or one of the IVA-VIA elements in the 24 g site (0, y, z). Some elements like Sn were found to be able to occupy the 24 g site (Co 4 Sb 12-y Sn y , CoSn 6 Te 6 , La 1-x Co 4 Sb 12-y Sn y , Yb 1-x Co 4 Sb 12-y Sn y , Tl 1-x Co 4-Sb 12-y Sn y ) as well as the 2a site (Sn y Ni 4 Sb 12-x Sn x ). The icosahedral cage in the skutterudite structure may be filled also by electronegative fillers such as I and Br.
The ability of these elements to enter different positions in the crystal lattice of the skutterudite structure is shown in Fig. 2(a), and has to be used as guide for the selection of the materials being in contact with skutterudite in the TE module. Taking into account that all highlighted elements may substitute atoms in the skutterudite lattice, one may suggest exchange reactions between materials that contain these elements and skutterudites. Respectively TE properties in such diffusion zones will be affected. Our analysis on alloying skutterudites with different elements show that the strongest effect on the transport properties is observed in case of a substitution in the 24 g site. Therefore, elements, such as Ga, In, Ge, Sn, Pb, Se and Te (site preference for 24 g) as constituent for the contact materials should be used only in form of highly stable compounds that may warrant a low interaction with skutterudites. Metals of the VIII group (substitution of the 8c site) can be used when the diffusivity of these elements in the skutterudite lattice will be suppressed by additional intermediate phase(s) formed during the interaction. Elements in the empty cells (Fig. 2a) are not known to enter the skutterudite lattice and therefore, they are the most promising components for contact materials. Considering this, besides d-elements, we also investigated the interaction with such electrically and thermally conductive materials as graphite and titanium nitride.

Interaction with Unfilled CoSb 3
Evaluation of the reaction products between pure elements and CoSb 3 by means of x-ray powder diffraction revealed an interaction in almost all cases: Ti, Zr, Hf, V-medium/ strong interaction, Nb and Ta-low interaction, with the exception of W, Re, graphite and TiN-no interaction. These observations (Fig. 2b) are in good correlation with the site preference in the skutterudite lattice (Fig. 2a), and with the constitution of the binary systems of these elements with antimony. The elements may be divided into five groups. The first group consists of noble metals (Mn, Ru, Pt, Ag and Cu) that show a very strong interaction. Note that these elements are main constituents of high temperature brazing alloys and sometimes they are used for contacting skutterudite based TE modules. Such a strong interaction might be explained via the existence of a low temperature eutectic in the binary systems with antimony. Most likely a similar behaviour is valid for Au, Zn and Cd. The second group that shows a strong interaction consists mainly of 8th group's elements, which have the ability to occupy the 8c site in the skutterudite structure. Medium interaction was observed for the majority of metals from groups IV-VI that do not tend to replace atoms in the crystal lattice of the skutterudite and these elements form rather stable antimonides. No interaction was observed in case of M = W, Re and C. This correlates with the absence of intermediate compounds in the binary M-Sb systems and ternary M-Co-Sb systems. Despite titanium interacting rather well with CoSb 3 , no interaction was observed for TiN, indicating that this compound coexists in a thermodynamic equilibrium with the skutterudite.

Interaction of Filled Skutterudites with Pure Elements
In a next step, our investigations were extended to the interaction of filled p-type skutterudites (Ce y Fe 3 CoSb 12 and MM y Fe 3 CoSb 12 ) and n-type skutterudites (Ba y Co 4 Sb 12 and Sr y Co 4 Sb 12 ) by means of diffusion couples. Phase constituents, their compositions and thickness of the reaction zone were determined by SEM ( Fig. 3 and 4) and EPMA ( Table 1). The thickness of the diffusion layer was measured as the distance from the initial ball to the phase boundary of the skutterudite matrix, however, it has to be noted that, in cases of interactions with elements that prefer to occupy the 8c crystallographic site (Fig. 2a), the composition of the skutterudite matrix close to this phase boundary is significantly enriched by these elements. The thickness of this diffusion zone was not established due to the absence of a phase boundary. It is only possible to estimate an extension of this zone, which appears in SEM images as skutterudite matrix with increased density. The thickness of the diffusion layer was normalized by using the initial diameter of the balls and assuming their spherical shape. Thickness of the diffusion layer and phases formed during interaction and compositions of intermediate phases are listed in Table 1. As arbitrary characteristic for the quality of the diffusion layers we used the term ''goodness'' that includes thickness, porosity, formation of cracks and quality of the bonding to metal and skutterudite matrix. Solubility of the third component in the parent phases, which formed in the diffusion layer, was also accounted for in order to estimate the ability of the metals to diffuse into these phases. We have to note that the conditions used for the diffusion experiment (the spheroid metal specimens were embedded into the skutterudite matrix) produced much more thermo-mechanical stress (crack formation) in comparison to a planar contact between electrode and TE materials in the module. Titanium A half-Heusler (HH) phase TiCoSb forms in contact with n-type skutterudite Ba 0.3 Co 4 Sb 12 (Fig. 3a).
Furthermore, we observe a very extensive zone of BaSb 3 , Ti 11-x Sb 8-y and very thin lavers, which, considering their brightness in backscatter electrons on SEM images and the constitution of the binary Ti-Sb phase diagram, [42] might be attributed to Ti 2 Sb and Ti 5 Sb 3 . Note, that no binary TiSb and TiSb 2 were determined in the diffusion zone being in good correlation with the constitution of the ternary system Ti-Co-Sb (Fig. 5a). The absence of TiSb most likely has its origin in particularities of phase equilibria in the quaternary system with Ba-Co-Ti-Sb, suggesting that equilibria with BaSb 3 are more thermodynamically stable. In contrast to n-type, a thick layer of TiSb 2 was detected in the diffusion zone, formed with p-MM y Fe 3 CoSb 12 (Fig. 4). In addition, we observed the formation of Ti 11-x Sb 8-y and Ti 3 Sb, but not of TiSb.
Vanadium Almost all binary V-Sb phases [47] formed in the diffusion zone. In addition, we detected the formation of an unreported compound, Ba 5 V 12 Sb 19?x for which the crystal structure was resolved to be isotypic with Ba 5-Ti 12 Sb 19?x . [48] The thicknesses of the diffusion layers were about 40 lm; they were very dense and showed a good bonding to skutterudite and parent metals. A significant solubility of Fe and Co was observed only for binary V 3 Sb 2 and ''V 5 Sb 4 ''. Despite a big similarity in the phase equilibria in ternary systems (Ti,V)-Co-Sb ( Fig. 5a and b), the constitution of the reaction zone differs significantly. Ternaries Ba 5 M 12 Sb 19?x exist also with Ti and Nb, [48] but we observed this structure only in the reaction couple with vanadium. In the other investigated systems binary antimonides and/or HH-phases were observed (Table 1) indicating that these phase equilibria are more energetically favourable.
Zirconium Both samples prepared with zirconium were lost during grinding. Probably, this happened due to almost no interaction (low bonding) with p-and n-type materials.
Hafnium, Niobium and Tantalum Only binary (Hf,Nb,Ta)Sb 2 were detected in the reaction zone for both types of materials. These diffusion layers were rather dense, but brittle. The n-type material showed stronger interaction with Hf (higher thickness of the diffusion layer). A very similar constitution of the diffusion zone was observed for Nb and Ta, but the thickness of the MSb 2 layer was significantly lower ( Table 1). The presence of MSb 2 in the diffusion zone is in accordance with the ternary systems (Hf,Nb)-Co-Sb ( Fig. 4c and d, no phase diagram is known for Ta-Co-Sb) and shows that the equilibrium MCoSb-MSb 2 separates CoSb 3 from the equilibria with other phases. This phase equilibrium is not present in the ternary systems (Ti,V)-Co-Sb ( Fig. 4a and b) indicating that the binary phases (Zr,Nb,Ta)Sb 2 are efficient diffusion barriers. Furthermore, the reduced thickness of the interaction zone for the samples with niobium and tantalum indicated that the low symmetry monoclinic structure of binary (Nb,Ta)Sb 2 (OsGe 2 -structure type) is more resistant against diffusion of foreign atoms than both, orthorhombic (Zr,Hf)Sb 2 (TiAs 2 -structure type) and tetragonal (Ti,V)Sb 2 (CuAl 2 -structure type). A small solubility of the third element was measured in all these phases, which could indicate a suppressed ability for foreign elements to diffuse through these crystal lattices. In addition to these two aspects (phase equilibria and crystal structure), there is also a certain correlation between thickness of diffusion layers MSb 2 and melting points of these compounds: VSb 2 (869°C), (Nb,Ta)Sb 2 (1089°C),  TiSb 2 (1051°C), ZrSb 2 (1363°C) and HfSb 2 (980°C). [42,47] Chromium Diffusion zones with a thickness of about 70 lm consisted of CrSb. No Fe or Co was detected in the diffusion layer. The diffusion layer was very brittle.
Molybdenum Brittle diffusion layers of Mo 3 Sb 7 with relatively good bonding ability were observed. No other binary molybdenum antimonides are known.
Tungsten No interaction was observed. Due to a significant mismatch of thermal expansion coefficients big cracks formed in the skutterudite matrix.
Iron, Cobalt and Nickel These elements showed a very strong interaction with p-and n-skutterudites. Diffusion layers with a thickness in the range of 100-200 lm were observed for Fe and Co whilst the interaction zone expanded up to 500 lm for Ni. Despite these extensive diffusion layers, we want to emphasize the good quality of the contact between the diffusion layer and the skutterudites. The composition of the skutterudites near the boundary of the reaction zone shifted due to the Fe/Co/Ni substitution in the 8c site with a respective change of the filling level in the 2a site. This exchange results in the formation of an additional quantity of the skutterudite that diminishes the porosity of the initial skutterudite matrix.
Ruthenium A strong interaction with the formation of RuSb 2 and Fe/Co/Ru substitution was found. For the n-type material we observe a complete substitution with the formation of BaRu 4 Sb 12 . The diffusion layer was rather thin (below 100 lm) and well bonded to (Ru) but the contact to the skutterudite matrix was poor.
Palladium and Platinum Both noble metals showed a very strong interaction with skutterudites. Palladium balls were completely resolved in the skutterudite matrix. The interaction of the platinum with skutterudites was similar to  [43] (b) M = V, [44] (c) M = Zr [45] and (d) M = Nb [46] that observed for nickel, and resulted in a very thick diffusion zone with numerous phases (Table 1).

Interaction of Filled Skutterudites with Binary Compositions
Thermal expansion coefficients for the metals, that show a low interaction with skutterudites, are significantly lower (a = 5 7 9 9 10 -6 K -1 ) than typical for p-type skutterudites (a = 11 7 12 9 10 -6 K -1 ) and are closer to that of n-type (a = 8 7 9 9 10 -6 K -1 ). [6] Only TiN has a thermal expansion coefficient (a = 9.4 9 10 -6 K -1 ) that fits rather well with both types of the TE materials, and this suggests necessarily development of individual diffusion barriers for different types of TE-legs. A significant mismatch in the thermal expansion coefficients of the materials is a main problem for the thermomechanical stability of the bonding, and it could be that materials which are effective protective barriers will not be suitable for practical application. Taking into account that thermal expansion might be adjusted by chemical substitution, our next experiment concerned binary compositions formed by elements that have low thermal expansion coefficients and did not enter the skutterudite lattice (Ti, V, Nb, Ta, Cr, Mo, W) with those that form skutterudites, but have a high bonding compatibility and high thermal expansion coefficients (Fe, Co and Ni; a & 13 9 10 -6 K -1 ).
The thickness of the diffusion layer and the bonding ability to skutterudites were compared with those observed for pure metals. Binary eutectics were used as compositions for this investigation, and in the cases of absence of eutectics (V-Fe and Fe-Cr), we used solid solutions. Structures of the diffusion zones with p-and n-skutterudites are shown in Figures S1-S3 (Supplementary  Information); the thickness of the diffusion zones is listed in Table 2. The quality of the interface between the diffusion zones and skutterudites/alloys was improved in almost all cases compared to those with pure elements (Fig. 3 and 4 and Table 1), and this effect we attribute to a better matching of the thermal expansion coefficients of the bonded materials. Effect of the alloying on the thickness of the interaction layer is rather different for p-and n-type materials. Dependences (Fig. 6c and d and 7) are not smooth, because the solid solutions exists only in the binary systems Fe-V and Fe-Cr, whilst the other investigated systems contain various phases with different interaction ability. Therefore, lines connecting experimental points in Fig. 6 and 7 are shown only as guides for the eyes.
The strongest ability to suppress interaction and diffusion was observed for binary alloys with V, Nb and Ta. A beneficial effect was observed for the n-type materials, for which in most cases the thickness of the diffusion zone was lower than for pure metals. A strong suppression of the interaction with n-type materials was also observed in case of Co 72 Mo 28 , Ni 44 Cr 56 and Co 59 Cr 41 . Note, that the effect of W-alloying on the suppression of the interaction is not as strong as one may suggest, considering the absence of an interaction between skutterudites and pure tungsten. The same concerns exist for the interaction of Nb depleted compositions with p-type materials (Fig. 6d). Whilst alloying by this element has generally a very beneficial effect, these two phase compositions interact with p-skutterudite even stronger than pure elements. This deteriorative effect cannot be explained only by particularities of the crystal structure of the second phase in these samples. Niobium-based binary Laves phases (with MgNi 2 and MgZn 2 type) coexist in equilibria with (Co) and (Fe) whilst NbNi 3 with TiCu 3 structure forms in the Nb-Ni system. The same structure types form also in binary systems with Ta and other elements, but their effect on the interaction with skutterudites differs (Fig. 7). The constitution of the reaction zone (Figs S2-S3, Supplementary information) was in many cases much more complicated than in case of the interaction with pure elements ( Fig. 3 and 4, and Table 2), and several unknown phases were detected and will be the subject of forthcoming investigations. Complete phase characterization will be published after a complete structural identification, however it became evident that the formation of particular phases in the diffusion zone were governed by the phase equilibria whilst the extensions of the diffusion zones depended on their crystal structure and the possibility to dissolve foreign atoms. The chemical stability of the contact (absence of interaction) was warranted under the condition that phases were in a thermodynamic equilibrium, and small extensions of their homogeneity regions of their phases indicated that they could act as efficient diffusion barriers.
These considerations are also helpful for the development of segmented TE legs. As example, the ternary systems (Ti,V,Zr)-Co-Sb contain two ''thermoelectric'' phases: CoSb 3 and HH-(Ti,V,Zr)CoSb. According to the constitution of these ternary systems, [43][44][45] these phases coexist in equilibrium (Fig. 5), providing a thermodynamically stable contact. In addition, HH phases are in equilibrium with metals or with M-Co intermetallics, which in the case of a good thermal and electrical conductivity may serve as chemically stable electrode materials, as the following examples show: the formation of HH was observed in the diffusion couple between n-type materials: (a) TiCoSb and Ti (in contact with skutterudite) and (b) V(CoFe)Sb and Fe 70 V 30 (in contact with the intermetallic phases). The interaction of Fe 70 V 30 with p-type material revealed the formation of HH-VFeSb (in contact with the intermetallic phase). Furthermore, the reaction of Fe 84 Ti 16 with both skutterudites results in the formation of a new off-stoichiometric compound, TiFe 1.33 Sb, [49] whilst off-stoichiometric HH-Nb 1-x CoSb forms during the interaction of Fe-Nb alloys with n-type skutterudites and it also coexists in equilibrium with CoSb 3 (Fig. 5d) according to the calculated phase diagram Nb-Co-Sb at 0 K. [46] The same considerations are valid for other promising candidates for diffusion barrier-antimony rich binary MSb 2 , which form equilibria with CoSb 3 (Fig. 5) and show only small solubility for the elements that are the main constituents of p-and n-skutterudites (Co, Fe and Ni). It is important to note that equilibria, existing in this ternary system, are present also in higher order systems, and therefore this methodology may be used for the determination of a proper diffusion barrier in any Sb-based TE material.
Besides the chemical compatibility of the contacts and particularities of the crystal structure of the contacted materials, one has to consider the thermo-mechanical stability of the bonding (mechanical properties, thermal expansion coefficients, etc.), and the possibility to engineer these properties via substitution without changing the phase equilibria.

Conclusions
Investigations of the interaction of antimony based skutterudites with pure metals and various binary compositions revealed two main chemical criteria for the development of a chemically stable and efficient diffusion barrier (or electrode) for the operation in contact with skutterudites at high temperatures. The material (a) has to be in thermodynamic equilibrium with the skutterudite, and (b) show a low solubility of elements that are constituents of electrode materials and skutterudites. The investigation revealed several materials, such as: pure metals, Nb, Ta, W and Re, high antimonides of transitional elements of the IV and V groups and nitrides (TiN). Based on the experimental results obtained and the analysis of literature data on the atom site preference in the skutterudites, five groups of elements were defined according to their ability to interact with skutterudites. It was shown that the quality of the interfaces between phases in the diffusion zone may be adjusted by chemical substitution, and this provides a better thermo-mechanical stability of contacts in TE modules. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.