Approximation of Composition and Temperature Dependent Heat Conductivity and Optimization of Thermoelectric Energy Conversion in Silicon–Germanium Alloys

We analyze the efficiency as thermoelectric energy converter of a silicon–germanium alloy with composition and temperature dependent heat conductivity. The dependency on composition is determined by a non-linear regression method (NLRM), while the dependency on temperature is approximated by a first-order expansion in the neighborhood of three reference temperatures. The differences with respect to the case of thermal conductivity depending on composition only are pointed out. The efficiency of the system is analyzed under the assumption that the optimal energy conversion corresponds to the minimum rate of energy dissipated. The values of composition and temperature which minimize such a rate are calculated as well.


State of the Art
Silicon-germanium alloys of the type Si c Ge 1−c , with composition c ∈ [0, 1] changing inside the system, are extensively used in modern technology [1] as, for instance, in design of thermoelectric energy generators [2]. The thermoelectric efficiency of energy generators η el = P el Q tot , with P el the electric power output, andQ tot thermal power entering in the system, is an increasing function the so called figure-of-merit Z = 2 σ e λ , where is the Seebeck coefficient, σ e the electrical conductivity, and λ the thermal conductivity [3,4]. It is easily verified the the physical dimension of Z is the inverse of temperature, namely K −1 . Hence, since it is sometime convenient to work with non-dimensional quantities, some authors denote as figure-of-merit the quantity ZT, with T the temperature. Critical analysis of assigning the figure-of-merit is given in [5,6]. If one remains in the frame of linear thermodynamics, for a system whose two sides are constantly maintained at different temperatures T h (the hottest temperature) and T c (the coldest one), it can be proved that the maximum efficiency is [4] η max = η C 1 − 1/ξ 1 + 1/ξ (1) wherein η C = 1 − T c / T h is the classical Carnot efficiency, and ξ ≡ √ ZT +1. Thus, the higher ξ, (and, hence, the smaller λ), the higher η max . However, any strategy for the reduction in λ should take into account the results of phonon hydrodynamics, wherein it is assumed that the transport of heat is due to the phonons [7]. Phonons are quasiparticles, which in a solid crystal form a rarefied gas, whose kinetic equation can be obtained similarly The frequencies ν N and ν R of normal and resistive processes, respectively, determine the characteristic relaxation times τ N = 1/ν N and τ R = 1/ν R . Purely diffusive heat transport takes over when there are many more R processes than N processes, i.e., when ν R tends to infinity and τ R = 1/ν R tends to zero. If, instead, there are only few R processes and many more N processes, then τ N = 1/ν N tends to zero, and a wavelike energy transport (second sound propagation) may occur. The total relaxation time τ can be calculated by according to the Mathiessen rule 1/τ = 1/τ R + 1/τ N , while the thermal conductivity is given by λ = c v τv 2 /3, where is the mass density, c v the specific heat per unit mass at constant volume, and v is the average of the phonons' speed. Thus, a reduction in λ would produce a reduction of τ, i.e., an increment of the phonon scattering [7,8], with a consequent increment of dissipation. Thus, numerator and denominator in the expression of η el cannot be controlled independently, and one should look for reductions of λ which produce moderate increment of phonon scattering.
Materials with composition-dependent thermal conductivity are often used to enhance Z [2,[9][10][11] since, by grading appropriately the stochiometry, λ can be reduced, achieving so a consequent increment of Z [4,7]. We also observe that a further improvement of Z could be obtained by increasing the product σ e as function of c. At the moment we are not aware of the data concerning the dependency of and σ e on c, but such a task could be considered in future studies. Herein, we limit ourselves to take into account the dependency of and σ e on temperature reported in [12], and use three different values of those quantities for three different operational temperatures. Different physical quantities as, for instance, the electronic part of the thermal conductivity, could be considered as well. The role of such material parameter has been considered in [13,14]. Indeed, the problem under consideration depends on several parameters of different nature. Herein, in order to obtain applicable results, we limited ourselves to consider a simple but meaningful case.
The thermoelectric efficiency of graded systems has been investigated in Refs. [15][16][17][18][19], where the dependency of the performance of a thermoelectric energy generator as function of the two independent parameters c and x, with c as the composition of the system, and x as the square root of the mean value of the temperature gradient applied to its boundaries, has been investigated extensively. We used the square root of the mean value of the temperature gradient and not the temperature gradient itself, in order to obtain a more manageable expression for the rate of energy dissipated. The first independent parameter, i.e., the composition cannot be tuned externally, since c is fixed after manufacturing the system. The second parameter, instead related to the applied temperature gradient, can be tuned externally. We have pointed out that the analysis of those systems yields new information on how manufacturing homogeneous thermoelectric generators made by silicon-germanium alloys, by determining the composition and the temperature gradient which optimize their efficiency. It is worth observing that our aim in [15][16][17][18][19] was not to obtain explicit values of the efficiency for given independent parameters, because such values depend on several physical quantities which can be changed arbitrarily, such as, for instance, the geometry of the system and the external electric field. This fact motivated our choice of determining the values of the independent physical parameters which minimize the energy dissipated along the thermoelectric process, whatever is the value of such energy. One could wonder if the minimum of energy dissipated corresponds to the optimal efficiency of the thermoelectric process. Our answer is affirmative, and is based on the following considerations. In the following, we disregard all the losses introduced in the production ofQ tot and in the management of the generated difference of electrical potential, and we focus only on the thermodynamic process inside a thermoelectric wire of length L. It consists in the generation of an electric potential after that an amount of heat per unit timeQ tot entered the system. Such a heat produces dissipation by Joule effect which, in any point z of the system and at any time t, is represented by the rate of energy dissipated E (c(z), x(z, t)) ≡Ē (z, t). Then, in a thermoelectric process of duration t 0 , the total energy dissipated is given by so that, being the integrand a positive quantity, the right-hand side of Equation (2) is minimum if, and only if, the integrand function is minimum in any point of the domain of integration [19]. We underline again that the previous analysis regards only the process of thermoelectric energy conversion inside the wire. A more complete analysis should include the dissipation due to the production ofQ tot , and that due to the transport and management of the obtained difference of electric potential, as well as the details of the processing parameters. However, such an analysis is outside the scopes of the present research, and is more pertinent to the field of engineering.
In order to make our investigation applicable to real cases, in [15][16][17][18][19] we used the experimental data on thermal conductivity of silicon-germanium alloys as function of the composition, at T = 300 K, T = 400 K, and T = 500 K [20][21][22]. From the disposition of the experimental points in the plane (c, λ(c)) it is evident that the thermal conductivity is very steep in the two regions close to c = 0, and c = 1. Indeed, around c = 0, due to the small value of c, the alloy can be considered as doped Ge, while, close c = 1, due to the high value of c, the alloy can be considered as doped Si. Thus, since doped crystalline semiconductors have reduced thermal conductivity with respect to the alloy, such a rapid decreasing is expected. Between those two steep parts of the curve, the experimental values of λ have small variation, so that λ presents a kind of wide minimum between the two dilute zones c ∈ [0, 0.1] and c ∈ [0.9, 1]. From the mathematical point of view, such a behavior is well represented by the sum of two exponential functions with A , B , D , E , F and G as unknown (non independent) parameters, to be determined by NLRM, under the constraints λ(0) = λ Ge and λ(1) = λ Si . Thus, as first step, we have looked for a fitting curve represented by the Equation (3), and applied the following iterative procedure [23]: • Starting from the disposition of the experimental data in the plane (c, λ(c)), we have assigned an initial estimated value of each parameter entering Equation (3); • We have generated the curve defined by the initial values of the parameters; • We have calculated the sum of the squares (i.e., the sum of the squares of the vertical distances of the experimental points from the curve); • We have adapted the parameters in such a way that the curve was as close as possible to the experimental points; • We have stopped the calculations when we have observed a negligible difference of results in successive iterations.
In this way, we have obtained the best-fit curve of the data as Some more details of NLRM applied to obtain a fit of the experimental data are given in Ref. [16]. Therein, it is illustrated the method for minimizing the mean distance between the experimental points and the fitting ones. Of course, such a distance cannot be reduced to zero, and this error affects the coefficients M, N, P, Q, and, as a consequence, the function λ given in Equations (4) and (5). We expect that this error will influence also the optimal values of c, T, and λ calculated in Section 3.
In [15][16][17][18][19], we have studied the effects of the action of an electric field E on a graded Si c Ge 1−c wire of length L, crossed by an electric current i. The right-hand side (z = L), was supposed to be at the temperature T h , while the left-hand side, (z = 0), was supposed to be at the temperature T c . We assumed that i was flowing uniformly from left to right, and that the heat rateQ tot was entering uniformly into the hot side of the system, giving rise to a heat flux q.
In our investigations, we applied the basic equations of thermoelectricity, whose physical meaning is discussed, for instance, in [4,24]. One of the key quantities in thermoelectricity is the form of the heat flux. Classically, the Fourier-like constitutive equation where Π is the Peltier coefficient, and Πi is the additional heat flux due to the circulation of the electric current i [4,7,25,26] is postulated. However, heat transport theory is currently broadening its field of applicability since, owing to the miniaturization, new phenomenologies, beyond the classical Fourier theory of heat conduction, have been discovered [27,28].
Those new phenomena depend on the relationship between the mean free path of the heat carriers , and the characteristic dimension of the conductor L, expressed by the Knudsen number Kn = / L. Fourier's law is valid when / L 1, namely, for L. However, Kn can increase for a reduction in L, as in miniaturization technologies. Thus, when the mean free path of the heat carriers is comparable to the characteristic dimension of the conductor, i.e., Kn 1, more complicated transport laws for the heat flux are necessary [19,27,28]. In the present investigation, the constitutive equation for q has been supposed to be where l denotes a characteristic-length vector, and b, (< 1), is a dimensionless physical parameter entering the effective thermal conductivity λ e f f ≡ λ(1 − b) [25,26]. In stationary situations, Equation (7) arises in thermomass (TM) theory [29,30] of heat transport. In a non-stationary case, the thermomass heat transport equation reads wherein τ tm is a relaxation time [25,26,29,30]. In thermomass description, the heat flux is generated by a gas of heat carriers, characterized by an effective mass density and flowing through the medium under the action of a thermomass-pressure gradient. This gas is made by massive quasi-particles of heat carriers, named thermons, which are nothing but the vibrations of the molecules generated by heating the conductor, with null rest-mass and dynamic mass which may be calculated from the Einstein's mass-energy duality. In gases and liquids, the thermons are supposed to be attached to the molecules or atoms of the medium. In solids, the thermomass gas coincides with the phonon gas for crystals, attached on the electron gas for pure metals, or just both of them for systems in which the heat carriers are phonons and electrons. The physical parameters entering Equation (8) are [19] τ tm = λ 2γρc 2 v T with the dimensionless parameter γ being the Grüneisen constant, standing for a dimensionless number which is called thermal Mach number of the drift velocity relative to the thermal-wave speed in the heat-carrier collection, and a characteristic-length vector. In fact, the physical dimensions of |l| are meters, as it can be directly inferred by the dimensional analysis of Equation (7). It characterizes the strength of the non-Fourier effects introduced by Equation (7) and, for conceivable values of q, attains values which are always much smaller than those of the mean-free path of the thermons [25,26]. Equation (7) is non-linear in the heat flux and also accounts for first-order non-local effects through the term ∇q · l. In [15][16][17][18][19], we applied Equation (7), since it contains the meaningful concept of effective thermal conductivity, accounting for the experimental evidence that the thermal conductivity is not independent of the heat flux. The constitutive equation for the current density is with E as the electric field, herein regarded as an external force applied to the system [4]. Finally, Equations (7) and (9) must be coupled with the energy-rate equation [4] ρ ∂u ∂t where u is the specific internal energy, and the quantity E · i is the rate of energy production due to the circulation of electric current. According to second law of thermodynamics, for such a system, the energy dissipated along the process, locally can be written as where σ s denotes the local entropy production. Our analysis has been carried out under the hypotheses that λ(c, T) λ(c, T h ), that both q and E depend only on the position on the longitudinal axis z, and that q and i are parallel. Then, by using Equation (11), together with the constitutive equations for the heat flux and for the electric current, after some manipulations the local rate of energy dissipated can be written as follows [16]: where E and l denote the mean values on the interval [0, L] of |E| and |l|, respectively, is the square root of the mean value of the temperature gradient. By Equation (12), we infer that the effective thermal conductivity λ(c)(1 − b) influences, in a meaningful way, the rate of energy dissipated. It is worth noticing that in deriving the expression (12) of the energy dissipated we did not use the second Kelvin relation Π = T because when the heat flux is given by Equation (7), such a relation could no longer be valid, as proved in [31].
Then, under the assumption that the optimal thermoelectric energy conversion corresponds to the minimum of the function E (c, x), we have determined the couples (c opt , x opt ) which minimize E (c, x) in different situations. Moreover, in correspondence of each minimum, we have obtained the value λ opt of the thermal conductivity. For the sake of illustration, in Table 1 are summarized the results obtained in Ref. [18]. Table 1. Results for (λ opt , c opt ) in Ref. [18], for L = 3 mm. To our best knowledge, there are not similar investigations in literature, so that we cannot compare the present results with pre-existing ones. Their experimental confirmation could follow by their possible application to design and manufacturing of thermoelectric energy converters.

The Present Research
Another parameter which can be externally controlled is the operational temperature, namely, the temperature of the conductor during the process. Thus, for the same system considered in Ref. [18], it seems important to obtain the couples (c opt , T opt ) which minimize the energy dissipated, once x is fixed. To achieve that task, we need to recalculate the energy dissipated as function of c and T, at constant x. A direct inspection of Equation (12) suggests that this can be obtained if, and only if, we are able to express λ as function of c and T. On the other hand, since at different temperatures correspond a different phonon scattering, we expect a different behavior of λ with respect to the case of constant temperature. Such a hypothesis seems to be confirmed by Table 2, wherein it is shown the thermal conductivity of a pure Si and pure Ge wire of length L = 3 mm at the temperatures T = 300 K, T = 400 K, and T = 500 K [18]. Table 2. Thermal conductivity of pure Si and pure Ge, for a wire of length L = 3 mm (Ref. [18]). We note a marked difference of the values of λ for the different temperatures. Thus, a dependency of λ on temperature must be taken into account. Since we have no experimental data of λ as function of temperature for Si − Ge wires of length L = 3 mm, as a first step we look for an approximation of it, in the neighborhood of the reference temperatures T = 300 K, T = 400 K, and T = 500 K. In more detail, we look for λ in the form  (2), with a suitable correction factor to be determined in order to keep the points of minimum of the temperature close to the reference temperature. For ϕ(T) and γ(T), instead, we use the zeroth order approximation given by Equation (4). We guess that, since the functions m(T), n(T), p(T), and q(T) enter the exponents appearing in the expression of λ(c, T), even a small variation of them will produce a sensible variation in λ, so that we expect a remarkable difference with respect to the isothermal case. This fact is in accordance with our hypothesis that variations in temperature influence remarkably the phonon scattering and, as a consequence, the thermal conductivity. Thus, around T = 300 K, for instance, we write n(T) = h 3 n 0 + h 4 n 1 (T −300) (15) p(T) = h 5 p 0 + h 6 p 1 (T −300) (16) q(T) = h 7 q 0 + h 8 q 1 (T −300) (17) where h i , i = 1, . . . 8, are the correction factors we are looking for, and m i , n i , p i , and q i , i = 0, 1, are the coefficients of the expansion. The quantities h i , i = 1, . . . 8, do not derive by a theoretical calculation, but are chosen empirically in order to obtain values of λ which remain limited for T ranging in the neighborhood of the reference temperatures. In future studies, we plan to investigate the possible existence of optimal values of them. Notice that we do not impose the conditions h 1 m 0 = M, h 3 n 0 = N, h 5 p 0 = P, and h 7 q 0 = Q, which would guarantee that for T = 300 K the fit in Equation (2) is recovered. We are aware that such a strategy could appear unusual at a first look. Our choice is motivated by the experimental evidence that, in materials in which λ depends on composition and temperature, there is a different rate of phonon scattering with respect to materials in which λ is insensitive to the temperature variation, so that we foresee different values of λ for the same composition. Our hypothesis is confirmed by numerical investigations, because the conditions h 1 m 0 = M, h 3 n 0 = N, h 5 p 0 = P, and h 7 q 0 = Q do not produce either an acceptable fit or a minimum value of the heat conductivity. Anyway, we underline again that here we do not obtain any fit of λ, because we do not have data on its dependency on T for Si − Ge wires of length L = 3 mm. Herein, we obtain only a first-order approximation of λ as function of T, and investigate how such temperature dependency influences the thermoelectric efficiency. Such a situation will be analyzed in more detail in Section 3. Then, we do not impose any limitation to Equations (14)- (17), and look for correction factors which yield a physically acceptable approximation of λ. In this way, we obtain a dependency of λ on c which is qualitatively similar to that obtained in [18], but having different values. In future studies, we aim at exploring how much the numerical error influences such differences, and if we can reduce them by improving our approximation of λ. The procedure illustrated above is also applied for the other two temperatures under consideration, by using different correction factors. In this way, we obtain an approximation of the function λ(c, T) which allows to express the rate of energy dissipated E as function of c and T. Then, we determine the conditions ensuring the optimal efficiency of the thermoelectric energy conversion by calculating the points of minimum of E (c, T). The plots of λ(c, T) for different temperatures are obtained as well.
The paper runs as follows.
In Section 2, we apply the approximation procedure described above to obtain the thermal conductivity of a wire of length L = 3 mm as function of c and T.
In Section 3, we study the first and second derivatives of λ(c, T) to calculate its minima around the temperatures under consideration. We prove that for each temperature there exists one, and only one, couple (c opt , T opt ) which minimizes E (c, T). Moreover, we discuss the results in view of our approximation of λ and explain how they can be used in manufacturing thermoelectric energy converters.

Approximation of Thermal Conductivity
Following the way paved in Ref. [18], herein we study the effects of the action of an electric field E on a graded Si c Ge 1−c wire of length L, crossed by an electric current i. The right-hand side (z = L), is supposed to be at the temperature T h , while the left-hand side, (z = 0), is supposed to be at the temperature T c . Moreover, i is flowing uniformly from left to right, and an amountQ tot of heat per unit time is entering uniformly into the hot side of the element, giving rise to a heat flux q.
In order to obtain an approximation of the heat conductivity, we use the values of the material constants M, N, P, and Q determined in [18]. Such values are quoted in Table 3. Table 3. The material parameters M, N, P, and Q in Equation (2) for a Si c Ge 1−c wire of length L = 3 mm (see Ref. [18]). Then, up to the second-order approximation of M(T), N(T), P(T), and Q(T), we can write

Temperature
and, consequently, To proceed further, we use the values of pure Si and pure Ge [18], which are quoted in Table 2.

Results
For the rigid conductor considered here (L = 3 mm), in Ref. [18] we have calculated the minima of the rate on energy dissipated at a fixed operational temperature, by regarding such a rate as a function of the composition c and of x (see Equation (12)). In the present paper, instead, we carry on our analysis by supposing that x is fixed, while the rate of energy dissipated depends on composition c and temperature T, through the heat conductivity λ(c, T). To calculate the couples (c opt , T opt ) which minimize E (c, T) around each temperature, we first observe that Equation (12) can now be rewritten as where and are constant. Then, the points of minimum of E (c, T) coincide with those of λ(c, T). Before calculating these minima some preliminary considerations are in order.
First of all, we observe that by Figures 1-3 it is evident that the heat conductivity, as function of composition, is more steep in the intervals 0 ≤ c ≤ 0.1 and 0.9 ≤ c ≤ 1. In fact, since the phonon scattering increases when a small amount of impurities of (Si in almost pure Ge, or Ge in almost pure Si), a subsequent reduction in thermal conductivity takes over for values of c close to the boundaries of the interval [0, 1]. If, instead, the quantities of Ge and Si are almost comparable, no further appreciable reduction in λ takes over, so that λ is almost constant in the interval 0.1 ≤ c ≤ 0.9. Thus, we expect that c opt lies in this interval. Moreover, since we used a first-order series expansion of the material functions M(T), N(T), P(T), and Q(T) around one of the reference temperatures T = 300 K, T = 400 K, and T = 500 K, we look for T opt in the neighborhood of such temperatures. The considerations above are confirmed by the results in Table 4, wherein the values of c opt , T opt , and λ opt are reported for the three temperatures. In order to give a first approximated estimation of the energy dissipated along a thermoelectric process, in Table 5 we show the calculated values of E opt , corresponding to the three reference temperatures, for the system considered in [12]. To obtain such values, we supposed T c = T opt , and T h = T opt +0.1 K, where T opt is the optimal temperature shown in Table 4, corresponding to the reference temperatures T = 300 K, T = 400 K, and T = 500 K, respectively. Such a small difference of temperature could appear unusual, but one must take into account that at very small scale even small difference of temperature produce very high gradients. Thus, our choice allows to obtain realistic values of E opt . For the sake of calculation, we have approximated the optimal values of temperature by omitting the decimal digits. For and σ e , we used the experimental data in Ref. [12]. As we observed above, in the genuinely non-linear regime the SKR is no longer true, in general [31]. Thus, in our computation, we used three different values of Π close but not equal to T h , namely Π(325 K) = 2.8 × 10 −2 V, Π(425 K) = 4 × 10 −2 V, and Π(527 K) = 6 × 10 −2 V. For E we have chosen the value of 1 V/m, and for l the value of 30 nm. Such a value derives by the fact that in thermomass theory the value |l| is ever smaller of the mean free path of phonons in the system considered [25,26], and that such a value has been estimated to be around 40 nm for Si-Ge alloys at room temperature [32]. For the material parameter b < 1 entering the effective thermal conductivity, we have used the empirical value 1 − b = 10 −2 . Table 5. Values of E opt calculated by using Equation (12) and the experimental data for the system considered in Ref. [12].

Temperature (K)
E By the results above, we infer a substantial coincidence of the values of c opt in the neighborhood of the different temperatures. Furthermore, the difference T opt − T is almost the same for the three reference temperatures. We note that by the plots of λ emerges that it is maximum for pure crystals, i.e., in the presence of an ordered crystal lattice that is capable to reduce the phonon scattering. In the two dilute zones c ∈ [0, 0.1] and c ∈ [0.9, 1], which correspond to doped Ge and doped Si, respectively, λ starts to decrease (in [0.9, 1] λ decreases going from Si to Ge), because the doping increases the phonon scattering. For c 0.5, the ordered structure is lost, and the phonon scattering is maximum, due to the disordered spatial disposition of the molecules of the alloy. Thus, it is plausible to suppose that the higher the alloy disorder, the lower the thermal conductivity, and this explains why the minima of λ are all very close to c = 0.5.
We note a discrepancy of the value of λ opt around T = 300 K with respect to the values around T = 400 K and T = 500 K. This is due to the fact that, as it can be seen by Table 1, the values of λ Si and λ Ge at T = 300 K are higher with respect to the values at T = 400 K and T = 500 K. This influences the corresponding values of the functions ϕ and γ and, as a consequence, the values of λ, as it is evident in Figures 1-3. Finally, this difference in the values of λ, produces also a difference in the values of the rate of energy dissipated, as shown in Table 5.
For the sake of comparison, in Figures 4-6, we show the curves λ = λ(c), as obtained in Ref. [18], (blue curve) and λ = λ(c, 300 K), as obtained in the present investigation (red curve). Doing that, one should take into account that the red curves do not follow by a fit of λ, but they follow by a first-order approximation of λ as function of T. Perhaps, such an approximation can be improved and the differences between the two curves can be reduced. We aim at doing that in future investigations. We can observe the two following marked differences: • The values of λ = λ(c, 300 K) are higher of the values of λ = λ(c), even in the dilute intervals 0 ≤ c ≤ 0.1 and 0.9 ≤ c ≤ 1; • In two narrow intervals close to c = 0 and c = 1, the curve λ = λ(c, 300 K) is less steep of the curve λ = λ(c). Thus, in the composition and temperature dependent representation of λ, the effect of phonon scattering in the decrement of λ is still present but is less dramatic, even for doped Si and doped Ge. We will test further such a result in next investigations.
The analysis developed in the present research yields new information on how manufacturing homogeneous thermoelectric generators made by Silicon-Germanium alloys, by determining the composition and the operational temperature which optimize their efficiency. In fact, we may conclude that: • Around T = 300 K, the optimal thermoelectric energy conversion is achieved at T = 325.268 K in a homogeneous silicon-germanium wire Si c Ge 1−c , with c = 0.492488; • Around T = 400 K, the optimal thermoelectric energy conversion is achieved at T = 425.473 K in a homogeneous silicon-germanium wire Si c Ge 1−c , with c = 0.479984; • Around T = 500 K, the optimal thermoelectric energy conversion is achieved at T = 527.669 K in a homogeneous silicon-germanium wire Si c Ge 1−c , with c = 0.49022. Our next investigation will be devoted to the following tasks: • To improve the approximation developed here, perhaps using a second-order series expansion; • To look for experimental data on the dependency of Seebeck coefficient and thermal conductivity on the composition, in order to obtain a fit also for those quantities; • To look for more data concerning the dependency of the material functions considered here on temperature, and try to get a fit as well.
In further researches, we also aim at investigating the same type of problem considered here, if the constitutive equation for the heat flux changes in q = −λ∇ T +µq · ∇q + 2 (∆q + ∇∇ · q) + Πi where µ is a material function and is the mean free path of the phonons, which in solid crystals are responsible of the heat transport. The previous constitutive relation follows by a non-linear generalization of the celebrated Guyer-Krumhansl equation [33], which, from the macroscopic point of view, can be obtained within the framework of Extended Irreversible Thermodynamics [7,34,35]. In addition to non-linear effects, represented by the non-linear dependency of λ on composition and temperature and by the term µq · ∇q, Equation (50) takes also into account non-local interactions of phonons through the additional term 2 (∆q + 2∇∇ · q). In general, such a term is negligible for bulk systems, as that analyzed in the present paper, but becomes important in nanosystems, where the characteristic length of the system is comparable to the mean free path of the heat carriers [7,27,28]. Thus, at nanometric scale, it would be interesting to analyze the problem considered here by using Equation (50) instead of Equation (8), in order to investigate the influence of non-local effects on the efficiency of thermoelectric energy conversion.
Author Contributions: Conceptualization, V.A.C. and P.R.; formal analysis, V.A.C. and P.R.; investigation, V.A.C. and P.R.; writing-original draft, V.A.C. and P.R.; writing-review and editing, V.A.C. and P.R. All authors have read and agreed to the published version of the manuscript.