crystals-13-00453

crystals-13-00453

Crystals 2023, 13(3), 453; doi:10.3390/cryst13030453

Article

Size Effect on the Thermal Conductivity of a Type-I Clathrate

Monika Lužnik 1,*, Günther Lientschnig ¹, Mathieu Taupin ¹, Andreas Steiger-Thirsfeld ², Andrey Prokofiev ¹ and Silke Paschen 1,*Orcid

¹

Institute of Solid State Physics, TU Wien, Wiedner Hauptstraße 8-10, 1040 Wien, Austria

²

USTEM, TU Wien, Wiedner Hauptstraße 8-10, 1040 Wien, Austria

*

Correspondence: budnowski@ifp.tuwien.ac.at (M.L.); paschen@ifp.tuwien.ac.at (S.P.)

Academic Editor: Konstanze R. Hahn

Received: 7 February 2023 / Revised: 27 February 2023 / Accepted: 28 February 2023 / Published: 5 March 2023

Abstract

:

Clathrates are a materials class with an extremely low phonon thermal conductivity, which is a key ingredient for a high thermoelectric conversion efficiency. Here, we present a study on the type-I clathrate La 1.2 Ba 6.8 Au 5.8 Si 38.8 □ 1.4 directed at lowering the phonon thermal conductivity even further by forming mesoscopic wires out of it. Our hypothesis is that the interaction of the low-energy rattling modes of the guest atoms (La and Ba) with the acoustic modes, which originate mainly from the type-I clathrate framework (formed by Au and Si atoms, with some vacancies □), cuts off their dispersion and thereby tilts the balance of phonons relevant for thermal transport to long-wavelength ones. Thus, size effects are expected to set in at relatively long length scales. The structuring was carried out using a top-down approach, where the wires, ranging from 1260 nm to 630 nm in diameter, were cut from a piece of single crystal using a focused ion beam technique. Measurements of the thermal conductivity were performed with a self-heating 3 ω technique down to 80 K. Indeed, they reveal a reduction of the room-temperature phonon thermal conductivity by a sizable fraction of ∼40 % for our thinnest wire, thereby confirming our hypothesis.

Keywords:

clathrates; thermal conductivity; thermoelectrics; size effects; mesowires

1. Introduction

In view of global warming and the recent energy crisis, sustainable alternatives to fossil fuels are more important than ever [1,2]. Whereas at present thermoelectric generators (TEGs) play only a minor role in the wide group of renewable energy sources, their ability to turn a temperature difference into electricity makes them uniquely suited to be used in combination with other generators to recover energy that is otherwise lost as waste heat [3]. To boost the energy output, however, it is necessary to maximize their thermoelectric conversion efficiency. This in turn requires the used materials to have a high thermoelectric figure of merit Z T .

A material’s (dimensionless) thermoelectric figure of merit is given by Z T = S 2 σ κ T , with the absolute temperature T, the Seebeck coefficient S, the electrical conductivity σ and the thermal conductivity κ . S, σ and κ are interrelated via their dependencies on the charge carrier concentration, which are such that maximizing Z T amounts to a compromise [4]. However, the fact that the thermal conductivity is determined by both an electronic and a phonon contribution opens up the possibility to reduce κ without significantly degrading the electronic properties. This is the basis for the so-called phonon glass–electron crystal concept, which has triggered the research on several materials classes. Clathrates, for instance, form cage-like structures that encapsulate heavy guest atoms (Figure 1a). The rattling of these guest atoms disrupts the heat-carrying phonons, thus leading to a remarkably low phonon thermal conductivity, while leaving the electrical conductivity essentially unaffected [5,6,7,8,9].

Using inelastic neutron and X-ray scattering, an avoided crossing between the acoustic phonon modes and the flat rattling modes was identified [10,11,12], as illustrated in Figure 1b. The strong interaction of these modes, together with the associated thermodynamic signatures, the particular temperature dependence of the thermal conductivity and the observed universal scaling of the room-temperature phonon thermal conductivity with the product of sound velocity and Einstein temperature of the lowest-lying rattling mode, were taken as evidence for a phonon Kondo effect as the underlying mechanism [13]. The suppression of high-frequency acoustic phonons not only causes remarkably low phonon thermal conductivities, but also means that low-frequency (long-wavelength) phonons have a greater share in the heat transport. In fact, calculations suggest that in certain clathrates more than 90% of the heat is carried by phonons with a wavelength longer than 10 nm [14]. Therefore, even mesoscopic structures (about 100 nm to a few μ m) made from those clathrates are expected to lead to a sizable reduction of the thermal conductivity.

We set out to test this conjecture by measuring individual type-I clathrate mesowires, fabricated from a clathrate single crystal using a focused ion beam (FIB). We study the electrical resistivity and thermal conductivity at temperatures between 300 K and 80 K, where the lattice thermal conductivity of bulk type-I clathrates was shown to be dominated by intrinsic phonon–phonon Umklapp scattering; notably, for two well-controlled sample series—one with a varying vacancy content and one with a varying charge carrier concentration—the curves collapsed above about 50 K, demonstrating that scattering of phonons from vacancies and electrons does not play a significant role at higher temperatures [13]. To avoid uncertainties in the thermal conductivity measurement, which are common in this temperature range, we implement a 3 ω technique in which a suspended mesowire acts as both the heater and the thermometer [15]. Our results confirm that, even with relatively thick wires (diameters < 600 nm), a sizable reduction of the phonon thermal conductivity is achieved.

A number of previous studies explored size effects on the thermoelectric properties of clathrates in polycrystalline bulk samples, using melt spinning [16,17,18], ball milling and (high-pressure torsion) hot pressing [18,19,20], high-pressure and high-temperature (HPHT) sintering [21,22], as well as approaches that introduced a secondary phase to stabilize small grains [23]. Whereas some of the results were promising, firm conclusions on the relation between the grain size and the lattice thermal conductivity remained elusive. This calls for a proof-of-concept study on single wires, as performed in the present work.

Figure 1. (a) Crystal structure [24] of a type-I clathrate, with the guest atoms (green) sitting inside cages formed by essentially covalently bonded host atoms (blue and orange). In the case of La 1.2 Ba 6.8 Au 5.8 Si 38.8 □ 1.4 (La-BAS), the green spheres denote Ba and La, the blue ones Si and the orange ones Au. The vacancies □ are distributed between Si and Au; (b) sketch of the phonon dispersion of type-I clathrates, as proposed in [13]. An acoustic phonon mode (the dashed line shows the dispersion assumed in the Debye model) hybridizes with a flat rattling mode (blue), resulting in a severe flattening of the former at large wave vectors q (red). This shifts the characteristic energy scale from the Debye energy ℏ ω D to the Einstein energy ℏ ω E of the lowest-lying rattling mode.

Figure 1. (a) Crystal structure [24] of a type-I clathrate, with the guest atoms (green) sitting inside cages formed by essentially covalently bonded host atoms (blue and orange). In the case of La 1.2 Ba 6.8 Au 5.8 Si 38.8 □ 1.4 (La-BAS), the green spheres denote Ba and La, the blue ones Si and the orange ones Au. The vacancies □ are distributed between Si and Au; (b) sketch of the phonon dispersion of type-I clathrates, as proposed in [13]. An acoustic phonon mode (the dashed line shows the dispersion assumed in the Debye model) hybridizes with a flat rattling mode (blue), resulting in a severe flattening of the former at large wave vectors q (red). This shifts the characteristic energy scale from the Debye energy ℏ ω D to the Einstein energy ℏ ω E of the lowest-lying rattling mode.

2. Materials and Methods

2.1. Samples

For our study, we chose La 1.2 Ba 6.8 Au 5.8 Si 38.8 □ 1.4 (La-BAS) single crystals, which were grown in-house using a floating zone technique. They were studied using scanning electron microscopy (SEM) with energy-dispersive X-ray spectroscopy (EDX) and wavelength-dispersive X-ray spectroscopy (WDX) analyses, transmission electron microscopy (TEM) as well as powder X-ray diffraction (XRD). This thorough analysis confirmed the phase purity of the grown crystals down to the nanometer scale [25]. As a result of the decreased charge carrier concentration, the substitution with La enhances the thermopower of the Ba 8 Au x Si 46 − x series substantially [25]. However, while the phonon thermal conductivity κ ph ∼ 1.8 Wm − 1 K − 1 at 300 K is within the range of the Ba 8 Au x Si 46 − x series [26], it is rather high compared to other clathrates, leaving room for improvement. We thus found this well-characterized compound to be ideally suited for our mesostructuring study.

A piece of crystal was investigated using Laue X-ray diffraction. We confirmed that it was indeed a single crystal (Figure 2a) and oriented it (Figure 2b). A crystal oriented along the 1–10 direction was then used to fabricate several mesowires with diameters of approximately 630 nm, 920 nm and 1260 nm.

2.2. Nanofabrication Process

The mesowires were fabricated with a top-down approach, in which the Ga + focused ion beam (FIB) of a FEI Quanta 200 3D Dual Beam system is used to form mesowires out of bulk single crystals. To protect the sample surface, first a platinum metal layer is deposited over the area of interest. Two trenches, one on each side of the metal line, are then cut with the FIB, so that a thin platelet, a so-called lamella, is obtained. In order to remove the lamella from the surrounding bulk crystal, it is tilted and cut free except for a small area on the side. It is then welded to a micro-manipulator needle using ion beam induced tungsten deposition. Finally, the remaining bridge between the lamella and the bulk sample is cut, and the lamella is transferred to and mounted on a transmission electron microscope (TEM) grid (Figure 3b), where the FIB is used to shape thin wires out of it (Figure 3c). The amorphous surface layer that forms due to bombardment with the high-energy Ga + ions is removed by an Ar + soft-milling procedure once the mesowires have almost the desired thickness. Next, the individual wires are again welded to the micro-manipulator tip and cut off from the rest of the material. The samples are then transported to a measuring platform, which was prepared beforehand using optical lithography. It consists of Ti/Au stripes (5 nm/45 nm thick, respectively) on an etched Si/SiO 2 substrate and can hold up to four wires. The wires are contacted there, using a combination of ion beam induced and electron beam induced tungsten deposition (Figure 3d).

2.3. Measurement Method

The mesowires were measured at temperatures between 300 K and 80 K in an LN 2 -flow cryostat, which is equipped with two radiation shields inside a vacuum chamber. Using a 4-point configuration, a 3 ω method was deployed to measure the sample’s thermal conductivity. In this technique, which is commonly used for microwires [27], nanotubes [28] and nanowires [29,30], the suspended sample acts as both heater and thermometer. This makes it less sensitive to radiation losses than other measurement methods, which is a big advantage for high-temperature measurements. Figure 4b shows a schematic of the set-up: a lock-in amplifier (Signal Recovery 7265) provides a sinusoidal voltage signal, which is turned into an AC current using a current source (Anmesys AMS220). The current is fed through the sample, and the voltage is picked up on the inner contacts. The V ( I ) dependence of the first harmonic signal provides information on the electrical resistance, whereas the thermal conductivity can be obtained from the third harmonic using the equation [15]

V 3 ω ≈ 4 I 3 R R ′ L π 4 A κ 1 + ( 2 ω γ ) 2 ,

(1)

where I is the root mean square of the applied AC current, ω is its angular frequency, R and R ′ are the sample’s resistance and its temperature derivative, γ is the thermal time constant and L and A are the length between the voltage contacts and the cross-section of the sample, respectively. This equation not only gives direct information on the thermal conductivity κ , but also via the thermal time constant γ , on the specific heat C P .

The measurement setup and program were tested using a 25 μ m thick Pt wire. The frequency dependent third harmonic at 300 K is shown in Figure 5a. The fitted curve yields κ = 71.9 Wm − 1 K − 1 , which is in excellent agreement with the literature value for Pt, κ = 71.6 Wm − 1 K − 1 [31]. The difference between the obtained specific heat ( C P = 137 Jkg − 1 K − 1 ) and the literature value ( C P = 133 Jkg − 1 K − 1 [31]) is slightly larger, but still only about 3%. With these results, we deemed the setup suitable for the mesowire measurements.

Compared to the microwire measurements, the characteristic frequency dependence of the mesowires is shifted to higher frequencies. As it can no longer be fully detected by the lock-in amplifier, we decided to use a different approach for determining the thermal conductivity: in the low-frequency limit, Equation (1) takes the form

V 3 ω ≈ 4 I 3 R R ′ L π 4 A κ ,

(2)

i.e., the third harmonic becomes nearly frequency independent. Although the information on the specific heat is lost, the thermal conductivity can be readily obtained by varying the excitation current at a fixed frequency. Figure 5b shows that the third harmonic signal of a 920 nm mesowire is linear in I 3 for the whole temperature range. The slopes of the linear fits are used to calculate κ .

3. Results

Figure 6a shows the temperature-dependent thermal conductivity of La-BAS mesowires with three diameters of roughly 630 nm, 920 nm and 1260 nm. The results are the average of measurements of two different mesowires, the dispersion being represented by the shaded area. The left axis shows the phonon thermal conductivity, which was obtained by subtracting the electron thermal conductivity from the measured total thermal conductivity. It is rather flat over the measured temperature range, with a small decrease at higher temperatures. The room temperature phonon thermal conductivity κ ph ranges from 1.29 Wm − 1 K − 1 for the thickest wires to 1.08 Wm − 1 K − 1 for the wires with 630 nm diameter. Compared to the bulk phonon thermal conductivity [25], this corresponds to a reduction of roughly 40% for the thinnest wires.

The electron thermal conductivity is plotted on the right axis of Figure 6a. It was estimated using the Wiedemann–Franz law κ el = L T / ρ , with the Lorenz number L = 2.44 × 10 − 8 V 2 K − 2 and ρ as shown in Figure 6b. With the highest obtained value being κ el ∼ 0.25 Wm − 1 K − 1 at room temperature, it is apparent that the thermal conductivity is governed by the phonon contribution.

As expected for a metal, the electrical resistivity of the La-BAS mesowires decreases with decreasing temperature. Comparing Figure 6a,b, one can see that, while the thermal conductivity decreases when reducing the diameter, the electrical resistivity somewhat increases. In terms of the overall thermoelectric conversion efficiency, the question of whether the positive effects outweigh the negative ones will thus have to be addressed.

4. Discussion

The phonon thermal conductivity of clathrates can be fitted using a modified Callaway model, as proposed in [13]. It uses the formalism of the well established Callaway model [32], but accounts for the suppression of high-frequency phonons by substituting the Debye temperature with the Einstein temperature of the lowest-lying rattling mode,

κ ph = k B 4 T 3 2 π 2 v s ℏ 3 ∫ 0 θ E / T τ x 4 e x ( e x − 1 ) 2 d x with x ≡ ℏ ω k B T ,

(3)

where τ is the relaxation rate obtained from three main contributions and Matthiessen’s rule

τ − 1 = τ D − 1 + τ U − 1 + τ B − 1 with τ D − 1 = A ω 4 , τ U − 1 = B ω 2 T , τ B − 1 = v s / λ .

(4)

The contributions τ D − 1 , τ U − 1 and τ B − 1 are the defect scattering rate, the high-temperature Umklapp scattering rate according to Klemens [33] and the boundary scattering rate, respectively.

Using Equations (3) and (4), and the values θ E = 82.3 K and v s = 3330 m/s [34], a good fit to the bulk phonon thermal conductivity of La-BAS is obtained, as shown in Figure 7a. From the fit, one obtains the defect scattering constant A = 5.74 × 10 3 K − 4 s − 1 , the Umklapp scattering constant B = 38.8 × 10 3 K − 3 s − 1 and the boundary scattering length λ = 4.11 × 10 − 5 m.

We then simulate curves assuming that the FIB processing changed the boundary scattering length to the diameters of the mesowires (1260 nm, 920 nm and 630 nm) but left the other scattering processes unaffected. These curves match the measured data quite well at high temperatures. Deviations are observed below about 270 K. They are likely due to additional defects introduced during the mesostructuring process. As shown in [13], defect scattering in type-I clathrates dominates at low temperatures but can be discarded as the dominating effect at room temperature.

For comparison, we also fit the bulk data with the original Callaway model with θ D = 348.3 K [34] (Figure 7b). This fit, though it naturally yields different fitting parameters, is of similar quality as the one using the modified Callaway model. However, the simulated curves, obtained by again changing the boundary scattering length to match the mesowire diameters, differ much less from the bulk curve and fail to capture the high-temperature mesowire data. This confirms the validity of the assumption that the (Kondo-like) interaction between the acoustic modes and the rattling modes effectively cuts off the acoustic phonon spectrum and thus shifts the weight of heat-carrying phonons to lower wavelength ones.

Focusing on the high-temperature region, where effects of point defect scattering were shown to be unimportant [13], we visualize the difference between these two models by plotting the phonon thermal conductivity vs. the boundary scattering length. Figure 8 shows the respective simulations at 300 K, together with the data points for bulk and mesowires. It is clear that only the modified Callaway model, where the Einstein temperature takes the role of the Debye temperature, captures our data. The fact that the boundary-scattering effect is governed by the Einstein and not by the Debye temperature makes type-I clathrates uniquely suited for size effect approaches. Moving on, it will be of great interest to test whether the thermal conductivity of clathrate mesowires with a lower Einstein temperature does in fact show an even stronger size dependence.

Finally, we estimate the effect of the mesostructure on the thermoelectric efficiency of La-BAS. Using the compound’s Seebeck coefficient measured in [25], we calculate the figure of merit, as shown in Figure 9. While Z T is almost identical for the mesowires with diameters of 1260 nm and 920 nm, a notable enhancement is observed for the thinnest mesowires (of 630 nm diameter), caused by the strong reduction of phonon thermal conductivity at this diameter.

For materials with more conventional phonon spectra, where the Debye temperature governs the phonon thermal conductivity, size effects are observed for much shorter dimensions. For instance, in Bi 2 Te 3 wires, a reduction of the room temperature phonon thermal conductivity of only ∼25% is observed for a diameter of 300 nm [35], and, in Bi 2 Se 3 , a reduction of ∼29% was found for a diameter of 200 nm [36]. The room temperature total thermal conductivities of 182 nm PbTe wires [37] and 50–100 nm PbSe wires [38] show a reduction of ∼50%, but, in the latter work, the decrease is partly attributed to the rough surface. While the total thermal conductivity of InAs nanowires of 40 nm diameter shows an 80% decrease, the decrease expected for 400 nm wires amounts to only ∼26% [39]. Most studies on Si have focused on extremely thin wires [40]. In a study comparing somewhat larger structures [41], the change in the total thermal conductivity between films with characteristic dimensions of 250 nm and 150 nm amounts to less than 10%, whereas, for nanowires of 50 nm and 20 nm, it differs by more than 45%, showing that substantial size effects in Si require nanostructures as opposed to mesostructures. As our work shows, this is in contrast to type-I clathrates, where sizable effects can already be achieved with mesostructures. Our proof-of-concept study opens the door for new approaches to enhance Z T in type-I clathrates, using mesostructures as opposed to nanostructures that are much more difficult to stabilize.

Author Contributions

Investigation, formal analysis and visualization, M.L.; investigation and software, G.L.; resources, M.T. (measuring platforms), A.S.-T. (mesowires), A.P. (bulk sample); writing—original draft preparation, M.L.; writing—review and editing, G.L. and S.P.; conceptualization, supervision, project administration and funding acquisition, S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Austrian Science Fund (FWF), Grant Nos. 29279, W1243 and I5868–FOR5249.

Data Availability Statement

Original data are available from the corresponding authors upon reasonable request.

Acknowledgments

Open Access Funding by the Austrian Science Fund (FWF).

Conflicts of Interest

The authors declare no conflict of interest.

References

IEA. World Energy Outlook 2022; International Energy Agency: Paris, France, 2022. [Google Scholar]

Dong, F.; Li, Y.; Gao, Y.; Zhu, J.; Qin, C.; Zhang, X. Energy transition and carbon neutrality: Exploring the nonlinear impact of renewable energy development on carbon emission efficiency in developed countries. Resour. Conserv. Recycl. 2022, 177, 106002. [Google Scholar] [CrossRef]

Pourkiaei, S.M.; Ahmadi, M.H.; Sadeghzadeh, M.; Moosavi, S.; Pourfayaz, F.; Chen, L.; Pour Yazdi, M.A.; Kumar, R. Thermoelectric cooler and thermoelectric generator devices: A review of present and potential applications, modeling and materials. Energy 2019, 186, 115849. [Google Scholar] [CrossRef]

Snyder, G.J.; Toberer, E.S. Complex thermoelectric materials. Nat. Mater. 2008, 7, 105–114. [Google Scholar] [CrossRef] [PubMed]

Nolas, G.S.; Slack, G.A.; Schujman, S.B. Semiconductor clathrates: A phonon glass electron crystal material with potential for thermoelectric applications. Semiconduct. Semimet. 2001, 69, 255. [Google Scholar] [CrossRef]

Paschen, S.; Carrillo-Cabrera, W.; Bentien, A.; Tran, V.H.; Baenitz, M.; Grin, Y.; Steglich, F. Structural, transport, magnetic, and thermal properties of Eu₈Ga₁₆Ge₃₀. Phys. Rev. B 2001, 64, 214404. [Google Scholar] [CrossRef]

Paschen, S.; Pacheco, V.; Bentien, A.; Sanchez, A.; Carrillo-Cabrera, W.; Baenitz, M.; Iversen, B.B.; Grin, Y.; Steglich, F. Are type-I clathrates Zintl phases and ‘phonon glasses and electron single crystals’? Phys. B Condens. Matter 2003, 328, 39. [Google Scholar] [CrossRef]

Bentien, A.; Christensen, M.; Bryan, J.D.; Sanchez, A.; Paschen, S.; Steglich, F.; Stucky, G.D.; Iversen, B.B. Thermal conductivity of thermoelectric clathrates. Phys. Rev. B 2004, 69, 045107. [Google Scholar] [CrossRef]

Takabatake, T.; Suekuni, K.; Nakayama, T.; Kaneshita, E. Phonon-glass electron-crystal thermoelectric clathrates: Experiments and theory. Rev. Mod. Phys. 2014, 86, 669. [Google Scholar] [CrossRef]

Christensen, M.; Abrahamsen, A.B.; Christensen, N.B.; Juranyi, F.; Andersen, N.H.; Lefmann, K.; Andreasson, J.; Bahl, C.R.H.; Iversen, B.B. Avoided crossing of rattler modes in thermoelectric materials. Nat. Mater. 2008, 7, 811–815. [Google Scholar] [CrossRef]

Euchner, H.; Pailhès, S.; Nguyen, L.T.K.; Assmus, W.; Ritter, F.; Haghighirad, A.; Grin, Y.; Paschen, S.; de Boissieu, M. Phononic filter effect of rattling phonons in the thermoelectric clathrate Ba₈Ge40+xNi6−x. Phys. Rev. B 2012, 86, 224303. [Google Scholar] [CrossRef]

Pailhès, S.; Euchner, H.; Giordano, V.M.; Debord, R.; Assy, A.; Gomès, S.; Bosak, A.; Machon, D.; Paschen, S.; de Boissieu, M. Localization of propagative phonons in a perfectly crystalline solid. Phys. Rev. Lett. 2014, 113, 025506. [Google Scholar] [CrossRef]

Ikeda, M.S.; Euchner, H.; Yan, X.; Tomeš, P.; Prokofiev, A.; Prochaska, L.; Lientschnig, G.; Svagera, R.; Hartmann, S.; Gati, E.; et al. Kondo-like phonon scattering in thermoelectric clathrates. Nat. Commun. 2019, 10, 887. [Google Scholar] [CrossRef]

Kaur, P.; Madsen, G.K.H.; Bera, C. Thermoelectric figure of merit and thermal conductivity of type-I clathrate alloy nanowires. MRS Commun. 2019, 9, 370–374. [Google Scholar] [CrossRef]

Lu, L.; Yi, W.; Zhang, D.L. 3ω method for specific heat and thermal conductivity measurements. Rev. Sci. Instrum. 2001, 72, 2996–3003. [Google Scholar] [CrossRef]

Yan, Y.; Tang, X.; Li, P.; Zhang, Q. Microstructure and thermoelectric transport properties of type-I clathrates Ba₈Sb₂Ga₁₄Ge₃₀ prepared by ultrarapid solidification process. J. Electron. Mater. 2009, 38, 1278–1281. [Google Scholar] [CrossRef]

Tomeš, P.; Yan, X.; Kastner, R.; Svagera, R.; Waas, M.; Eilertsen, J.; Weidenkaff, A.; Paschen, S. Thermoelectric properties of meltspun Ba₈Cu₅(Si,Ge,Sn)₄₁ clathrates. J. Alloys Compd. 2016, 654, 300–307. [Google Scholar] [CrossRef]

Christian, R.; Ikeda, M.; Lientschnig, G.; Prochaska, L.; Prokofiev, A.; Tomeš, P.; Yan, X.; Zolriasatein, A.; Bernardi, J.; Schachinger, T.; et al. Nanostructured clathrates and clathrate-based nanocomposites. Phys. Status Solidi A 2016, 213, 784–801. [Google Scholar] [CrossRef]

Zolriasatein, A.; Yan, X.; Bauer, E.; Rogl, P.; Shokuhfar, A.; Paschen, S. Influence of PCA on thermoelectric properties and hardness of nanostructured Ba–Cu–Si clathrates. Mater. Des. 2015, 87, 883–890. [Google Scholar] [CrossRef]

Yan, X.; Falmbigl, M.; Rogl, G.; Grytsiv, A.; Prokofiev, A.; Bauer, E.; Rogl, P.; Zehetbauer, M.; Paschen, S. High-pressure torsion to improve thermoelectric efficiency of clathrates? J. Electron. Mater. 2013, 42, 1330–1334. [Google Scholar] [CrossRef]

Sun, B.; Jia, X.; Zhao, J.; Li, Y.; Liu, H.; Ma, H. Effects of pressure on the microstructure and simultaneous optimization of the electrical and thermal transport properties of Yb0.5Ba7.5Ga₁₆Ge₃₀. Inorg. Chem. 2018, 57, 3323–3328. [Google Scholar] [CrossRef]

Sun, B.; Zhao, J.; Li, Y.; Cao, L.; Yang, Y.; Fan, X.; Liu, X.; Wang, C.; Huang, X.; Wang, X.; et al. Effect of HPHT processing on structural and thermoelectric properties of low-cost type-I clathrate Ba₈Cu₆Si₄₀. J. Phys. Chem. C 2020, 124, 9082–9088. [Google Scholar] [CrossRef]

Yan, X.; Bauer, E.; Rogl, P.; Bernardi, J.; Prokofiev, A.; Paschen, S. Thermoelectric properties and stability of nanocomposites type-I clathrate Ba-Cu-Si with SiC. Z. Anorg. Allg. Chem. 2020, 646, 1055–1061. [Google Scholar] [CrossRef]

Momma, K.; Izumi, F. VESTA3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]

Prokofiev, A.; Sidorenko, A.; Hradil, K.; Ikeda, M.; Svagera, R.; Waas, M.; Winkler, H.; Neumaier, K.; Paschen, S. Thermopower enhancement by encapsulating cerium in clathrate cages. Nat. Mater. 2013, 12, 1096–1101. [Google Scholar] [CrossRef]

Candolfi, C.; Aydemir, U.; Baitinger, M.; Oeschler, N.; Steglich, F.; Grin, Y. High temperature thermoelectric properties of the type-I clathrate Ba₈AuxSi46−x. J. Appl. Phys. 2012, 111, 043706. [Google Scholar] [CrossRef]

Bhatta, R.P.; Annamalai, S.; Mohr, R.K.; Brandys, M.; Pegg, I.L.; Dutta, B. High temperature thermal conductivity of platinum microwire by 3ω method. Rev. Sci. Instrum. 2010, 81, 114904. [Google Scholar] [CrossRef]

Choi, T.Y.; Poulikakos, D.; Tharian, J.; Sennhauser, U. Measurement of thermal conductivity of individual multiwalled carbon nanotubes by the 3ω method. Appl. Phys. Lett. 2005, 87, 013108. [Google Scholar] [CrossRef]

Chien, C.H.; Lee, P.C.; Tsai, W.H.; Lin, C.H.; Lee, C.H.; Chen, Y.Y. In-situ observation of size and irradiation effects on thermoelectric properties of Bi-Sb-Te nanowire in FIB trimming. Sci. Rep. 2016, 6, 23672. [Google Scholar] [CrossRef]

Li, L.; Jin, C.; Xu, S.; Yang, J.; Du, H.; Li, G. Thermal conductivity of a single Bi0.5Sb1.5Te₃ single-crystalline nanowire. Nanotechnology 2014, 25, 415704. [Google Scholar] [CrossRef]

Properties of solids. In CRC Handbook of Chemistry and Physics; David, R.L. (Ed.) CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]

Callaway, J. Model for lattice thermal conductivity at low temperatures. Phys. Rev. 1959, 113, 1046–1051. [Google Scholar] [CrossRef]

Klemens, P.G. Thermal conductivity and lattice vibrational modes. In Solid State Phys.; Seitz, F., Turnbull, D., Eds.; Academic Press: Cambridge, MA, USA, 1958; Volume 7, pp. 1–98. [Google Scholar] [CrossRef]

Ikeda, M.S. Mechanism of the Thermal Conductivity of Type-I Clathrates. Ph.D. Thesis, TU Wien, Vienna, Austria, 2015. [Google Scholar]

Muñoz Rojo, M.; Abad, B.; Manzano, C.V.; Torres, P.; Cartoixà, X.; Alvarez, F.X.; Martín Gonzalez, M. Thermal conductivity of Bi₂Te₃ nanowires: How size affects phonon scattering. Nanoscale 2017, 9, 6741–6747. [Google Scholar] [CrossRef] [PubMed]

Dedi; Lee, P.-C.; Wei, P.-C.; Chen, Y.-Y. Thermoelectric characteristics of a single-crystalline topological insulator Bi₂Se₃ nanowire. Nanomaterials 2021, 11, 819. [Google Scholar] [CrossRef] [PubMed]

Roh, J.W.; Jang, S.Y.; Kang, J.; Lee, S.; Noh, J.-S.; Kim, W.; Park, J.; Lee, W. Size-dependent thermal conductivity of individual single-crystalline PbTe nanowires. Appl. Phys. Lett. 2010, 96, 103101. [Google Scholar] [CrossRef]

Liang, W.; Rabin, O.; Hochbaum, A.I.; Fardy, M.; Zhang, M.; Yang, P. Thermoelectric properties of p-type PbSe nanowires. Nano Res. 2009, 2, 394–399. [Google Scholar] [CrossRef]

Swinkels, M.Y.; van Delft, M.R.; Oliveira, D.S.; Cavalli, A.; Zardo, I.; van der Heijden, R.W.; Bakkers, E.P.A.M. Diameter dependence of the thermal conductivity of InAs nanowires. Nanotechnology 2015, 26, 385401. [Google Scholar] [CrossRef]

Boukai, A.I.; Bunimovich, Y.; Tahir-Kheli, J.; Yu, J.K.; Goddard Iii, W.A.; Heath, J.R. Silicon nanowires as efficient thermoelectric materials. Nature 2008, 451, 168–171. [Google Scholar] [CrossRef] [PubMed]

Ju, Y.S. Phonon heat transport in silicon nanostructures. Appl. Phys. Lett. 2005, 87, 153106. [Google Scholar] [CrossRef]

Figure 2. (a) Laue X-ray diffraction pattern of the used La-BAS single crystal; (b) oriented single crystal.

Figure 2. (a) Laue X-ray diffraction pattern of the used La-BAS single crystal; (b) oriented single crystal.

Figure 3. Scanning electron microscopy images of (a) the bulk single crystal with cut out lamellas; (b) a lamella; (c) mesowires cut from the lamella; and (d) a single mesowire on the measuring platform.

Figure 3. Scanning electron microscopy images of (a) the bulk single crystal with cut out lamellas; (b) a lamella; (c) mesowires cut from the lamella; and (d) a single mesowire on the measuring platform.

Figure 4. (a) measurement setup; (b) schematic of the measurement circuit, adapted from [15].

Figure 4. (a) measurement setup; (b) schematic of the measurement circuit, adapted from [15].

Figure 5. (a) Frequency dependence of the third harmonic voltage signal of a Pt test wire at 300 K, measured with 3.5 mA; (b) current dependence of the third harmonic of a 920 nm La-BAS mesowire at different temperatures. Prior to measuring the V ( I ) dependencies, frequency sweeps were performed to select a measuring frequency at which the 3 ω signal is nearly frequency independent and thus Equation (2) is valid. For better visibility, not all measured temperatures are included in this plot.

Figure 5. (a) Frequency dependence of the third harmonic voltage signal of a Pt test wire at 300 K, measured with 3.5 mA; (b) current dependence of the third harmonic of a 920 nm La-BAS mesowire at different temperatures. Prior to measuring the V ( I ) dependencies, frequency sweeps were performed to select a measuring frequency at which the 3 ω signal is nearly frequency independent and thus Equation (2) is valid. For better visibility, not all measured temperatures are included in this plot.

Figure 6. Temperature-dependent (a) phonon thermal conductivity (full symbols, left axis) and electron thermal conductivity (open symbols, right axis), as estimated from the Wiedemann–Franz law and (b) electrical resistivity.

Figure 6. Temperature-dependent (a) phonon thermal conductivity (full symbols, left axis) and electron thermal conductivity (open symbols, right axis), as estimated from the Wiedemann–Franz law and (b) electrical resistivity.

Figure 7. Temperature-dependent phonon thermal conductivity of a bulk La-BAS single crystal and mesowires fabricated from the same batch (symbols), the best fit to the bulk single crystal (full blue line) and simulated curves (dashed lines) for which the boundary scattering length is fixed to the mesowire diameters and all other parameters are kept unchanged, using (a) the modified Callaway model with θ E = 82.3 K and (b) the original Callaway model with θ D = 348.3 K. Data for the bulk single crystal are taken from [34]. Only the modified Callaway model can capture the overall high-temperature behavior. Deviations at low temperatures are attributed to small changes in the defect scattering rate due to the FIB processing.

Figure 7. Temperature-dependent phonon thermal conductivity of a bulk La-BAS single crystal and mesowires fabricated from the same batch (symbols), the best fit to the bulk single crystal (full blue line) and simulated curves (dashed lines) for which the boundary scattering length is fixed to the mesowire diameters and all other parameters are kept unchanged, using (a) the modified Callaway model with θ E = 82.3 K and (b) the original Callaway model with θ D = 348.3 K. Data for the bulk single crystal are taken from [34]. Only the modified Callaway model can capture the overall high-temperature behavior. Deviations at low temperatures are attributed to small changes in the defect scattering rate due to the FIB processing.

Figure 8. Phonon thermal conductivity as a function of boundary scattering length at 300 K, simulated using the modified (black) and the original Callaway model (red). Squares denote mesowire data, the circle the bulk data point taken from [34].

Figure 8. Phonon thermal conductivity as a function of boundary scattering length at 300 K, simulated using the modified (black) and the original Callaway model (red). Squares denote mesowire data, the circle the bulk data point taken from [34].

Figure 9. Z T as a function of T. The behavior of the mesowires with a 630 nm diameter at high temperatures points towards an improvement as a result of mesostructuring.

Figure 9. Z T as a function of T. The behavior of the mesowires with a 630 nm diameter at high temperatures points towards an improvement as a result of mesostructuring.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).