Electronic Transport in Alloys with Phase Separation (Composites) ()
1. Introduction
The performance of a thermoelectric material for cooling of power generation applications, or more generally, for energy conversion, are directly related to the dimesionless figure of merit defined by
ZT=S2Tρκ, (1)
where
ρ=1/σ . S,
κ ,
σ and
ρ are the Seebeck (or thermopower) coefficient, thermal conductivity, electrical conductivity, and electrical resitivty, respectively. T is the absolute temperature. The larger is ZT, the larger the thermoelectric performance of the material is. For many decades,
ZT≈1 was a practical upper limit realized in real materials. A further increase of ZT was limited by the fact that the transport coefficients occurring in Equation (1) are generally not independent from each other. In spite of these restrictions, during the last decade, progress was achieved by so-called nanostructured materials such as superlattices, quantum dots, nanocomposites, and nanowires.
In order to achieve further progress in the field of nanostructured materials with improved ZT, (1) a fundamental understanding of the carrier transport in these complex materials is necessary, and (2) effective experimental methods for designing, producing and measuring new material compositions with nanocomposite-structures are to be applied.
Regarding point (1), today it is generally accepted that application of the established classical theories for calculation of the electronic transport as the Boltzmann transport equation (BTE) can no longer be applied for the nanocomposite materials, because many of the characteristic lengths in the nanocomposite materials are smaller than the electron de Broglie wavelength (see, e.g., [1] , Figure 6 therein). Instead, modelling transport in nanocomposites is generally assumed to require more powerful tools such as non-equilibrium Green’s functions, which explicitly take account of electron wave effects [1] .
Regarding point (2), the combinatorial development of materials seems to be a proper method for experimental studies. From this method (applied to a certain material system) a large range of different compositions can be realized concurrently on one large substrat, for instance by deposition of thin films by simultaneous co-sputtering from two or three targets on a large substrat. Subsequently the transport coefficients can be measured with high lateral resolution applying measuring equipment such as the Potential-Seebeck Microprobe, the 3-Omega-method and the 4-point method for S,
κ , and
σ , respectively. The advantage of these methods is the fact that quick and precise measurements of the transport coefficients at identical positions over the complete substrate is possible. However, there is also a disadvantage of this method. Because in a continuous sample with a lateral concentration gradient the electrochemical potential
˜μ is the same in the whole sample, both the topological structure (atomic configuration) and electronic structure may be completelly different from those ones, if the samples with a certain composition are produced individually.
In the following we shall challenge these two points of view. We shall show that (1), the classical theories can, after all, be applied successfully for nanocomposites, particularly with respect to the electronic transport and its relation with the atomic structure and (2), the combinatorial development of materials can provide unfeasible results if applied noncritically.
On the way to this awareness, studies of amorphous alloys, as for instance a-Cr1−xSix alloys, have played an important role. Additional microstructure analyses of a series of amorphous transition-metal-alloys [2] [3] [4] [5] [6] , have shown that, for large concentration ranges, they are real nanocomposites composed of different (amorphous) phases, where the phase grains are found to be of the order of 1 - 2 nm [5] [6] . In correspondence to these analytical results, the a-Cr1−xSix thin films can be described as nanocomposites consisting of the phases
a-Cr1−xASixA (≡phase A) and
a-Cr1−xBSixB (≡phase B) where
xA≈0.25 and
xB≈1 .
In Sections 2.1-2.4 formulas for
σ ,
κ , S and R (the Hall coefficient) in composites will be derived and, in Section 2.5, compared with other published formulas. In Section 2.6 percolation elements will be additionally introduced in the formulas for two-phase composites. In Section 3.1 the classical formulas for the transport coefficients of the phases will be summarized applicable to large phase grains, which in Section 3.2 are extended to nanocomposites. In Sections 4-6 it will be shown that under certain conditions a discontinuity in the concentration dependence of the thermopower can occur, that the classical thermopower formula is to be supplemented by an additional term to be complete and that a noncritical application of the method of the combinatorial development of materials can provide unfeasible results.
The puzzling phenomenons: why there are simple metals with positive thermopower? and what is the reason for the phenomenon of the “Giant Hall effect”? and what is the reason for the fact that amorphous composites can exist at all? have been considered in the Section 5 and Section 3.2.1.
2. Electronic Transport in Composites
Let us consider a two phase-composite consisting of the phases
i=A,B in a symmetrical fashion regarding the average geometric form of the phase grains and without preferred orientations. Let us assume that each phase can be characterized by a set of transport coefficients
σi ,
κe,i ,
Si and
Ri for the phase i, which are the electrical conductivity, electronic contribution to the thermal conductivity, Seebeck coefficient and the Hall coefficient, respectively, in the phase i. The corresponding transport coefficients of the composite,
σ ,
κe , S and R are to be calculated, if the
σi ,
κe,i ,
Si and
Ri are known. The discussion will be restricted to small temperature gradients, small and constant electric and magnetic fields,
E and
H , respectively.
2.1. Electrical Conductivity
Applying effective medium theory (EMT), let us derive the relation between
σi and
σ , the electrical conductivities in the phases i and the composite, respectively. The strategy underlying the EMT is the following: a single phase grain of the phase i is considered to be completely embedded in an effective medium consisting of the two phases randomly arranged and characterized by the total transport coefficients. At the boundary face between this single phase grain and the surrounding effective medium continuity of the current densities and potentials and their gradients are to be saved.
The local electric current density
Ji can be written as
Ji=σiEi=−σigradφi. (2)
In analogy to Equation (2) we write for the electric current density J in the specimen
J=σE=−σgradφ. (3)
φi and
φ are the local and total electrostatic potential and
Ei=−gradφi and
E=−gradφ are the local and total electric field, respectively.
Now we demand that the total current density is equal to the average of the local current density,
J=〈Ji〉, (4)
and the same for the potential
gradφ=〈gradφi〉. (5)
Let us assume a spherical inclusion of the phase i with the radius
r0 , embedded in an uniform medium with the average conductivity
σ and that for this enclosed phase i the local transport equation Equation (2) holds. With Equations (3)-(5) we get
J=−σ〈gradφi〉=−〈σigradφi〉. (6)
The local potential
φi obeys the Laplace equation
Δφi=0, (7)
where the boundary conditions
φ=φi, (8)
Jn=Jni (9)
at
r=r0 are to be fulfilled, which is equivalent to
σgradrφ=σigradrφi (10)
at
r=r0 .
gradr is the gradient into the radial direction.
Jni and
Jn are the normal components of the current density in the sphere i and the surrounding effective medium, respectively. Equation (7) has the solution
φi=φ0+raicosϑ (11)
within the sphere i, and
φ=φ0+(br+cr30/r2)cosϑ (12)
within the effective medium.
ϑ is the angle between the direction of
E and the position vector
r with
|r|=r .
ai,b,c and
φ0 are constants. With
gradrφi=aicosϑ (13)
following from Equation (11), Equation (6) can be written as
σ〈ai〉=〈aiσi〉. (14)
With the boundary condition Equation (8) it follows that
ai=b+c, (15)
and with Equations (10)-(12) for
r=r0 ,
σ(b−2c)=σiai. (16)
Equations (15), (16) resolved for
ai providing
ai=3bσσi+2σ (17)
and introduced in Equation (14) provide
σ〈1σi+2σ〉=〈σiσi+2σ〉. (18)
Replacing the averages
〈…〉 by
Σiυi… we get
∑i υiσi−σσi+2σ=0, (19)
where
υi is the volume fraction of the phase i.
Equation (19) is the EMT formula for
σ . It holds for composites with more than two phases as well. The first authors who derived Equation (19) independently from each other, were Odelevskii [7] and Landauer [8] .
2.2. Thermal Conductivity
The schema shown in Section 2.1 can be applied to all the other transport parameters provided that a current density can be defined which is a function of a gradient of a potential, analoguously to Equation (2). For the thermal conductivity the corresponding relations read
JQ=−κegradT (20)
and
JQi=−κe,igradTi, (21)
where
JQ and
JQi are the total and local electronic thermal current density, respectively. Applying the same formalism as described in Section 2.1, however where the potentials
φ and
φi are replaced by the potentials T and
Ti , respectively, it follows that
∑i υiκe,i−κeκe,i+2κe=0. (22)
For the total thermal conductivity composed of both an electronic contribution
κe and an nonelectronic one
κne , it follows an analogous formula
∑i υiκi−κκi+2κ=0, (23)
where
κi is given by
κi=κe,i+κne,i (24)
provided that interactions between the various modes of heat transport can be neglected.
κne,i is the non-electronic contribution to
κi .
2.3. Thermopower
The following derivation will be based on J and
JS , the electric and entropy-flux density, respectively. The corresponding local flux densitities,
Ji and
JS,i , can be written as [9] [10]
Ji=σi[grad(˜μi/|e|)−SigradT], (25)
JS,i=σiSigrad(˜μi/|e|)−(σiS2i+κe,i/T)gradT, (26)
where
|e| is the elementary charge. (
−|e| is the charge of the electron.)
˜μi is the electrochemical potential in the phase i.
According to the strategy underlying the EMT, we demand continuity of the entropy-flux density and the electrochemical potential and their gradients, at the boundary face between a single phase grain and its surrounding (effective medium), where additionally
J=Ji=0 is to be fulfilled. Setting
Ji=0 in Equation (25), and inserting into Equation (26), one obtains for the local entropy-flux density,
JS,i=−κe,iSiTgrad(˜μi/|e|). (27)
In analogy to Equation (27) we write for the total entropy-flux density
JS in the specimen
JS=−κeSTgrad(˜μ/|e|), (28)
where
˜μ is the electrochemical potential of the composite.
The Equations (27) and (28) have the same structure as Equations (2) and (3); that is why, we can apply the same formalism as described in Section 2.1, however with other starting equations,
JS=〈JS,i〉, (29)
grad˜μ=〈grad˜μi〉, (30)
and we obtain finally
∑i υiκe,i/Si−κe/Sκe,i/Si+2κe/S=0. (31)
Considering the Wiedemann-Franz rule, Equation (31) can be transformed to
∑i υiσi/Si−σ/Sσi/Si+2σ/S≈0. (32)
Equation (32) was derived on basis J &
JQ , Equation (31) on basis J &
JS (Sonntag, [11] and [12] , respectively). For the one-band models, it makes little difference whether one chooses to interpret thermoelectric phenomenons on the basis of J &
JQ or J &
JS . However, for two-band or multiband models characteristic for semiconductors, the use of J &
JQ leads to complications, which are not to be expected, if J &
JS is used as a basis (Harman and Honig [9] , p. 28).
2.4. Hall Coefficient
Let us consider a non-magnetic two-phase composite. Under the same conditions as assumed in Section 2.1 the local electric current density in a single grain of the phase i (i = A or B) can be written as
Ji=↔σiEi, (33)
where
Ei and
↔σi are the electric field and the magnetoconductivity tensor [13] in this grain. For the electric current density outside of this grain we write analogously
J=↔σE, (34)
where
E and
↔σ are the electric field and the magnetoconductivity tensor outside of this grain (effective medium). For the determination of the coefficients in
↔σi we start with the equation for
Ji under the influence of an electric and magnetic field, [9] [14] [15]
Ji=e2iK11,iEi+e3imiK12,i(Ei×B)+e4im2iK13,iB(Ei⋅B). (35)
Krs,i are the transport integrals defined by
Krs,i=−43mi∫Erτsi1+μ2iB2∂f0i(E,T)∂ENi(E)dE (36)
with the Fermi-Dirac distribution function,
f0i(E,T)=11+eE−˜μikBT. (37)
Ni(E) ,
mi and
τi are the density of states, the effective mass and the relaxation time, respectively, of the carriers in the phase i. E and
kB are the energy and the Boltzmann constant, respectively.
ei=−|e| and
+|e| for electrons and holes, respectively. The third summand in Equation (35) disappears only if
Ei (or
Ei ) is always perpendicular to
B . In a composite, however,
B and
Ei (or
Ei ), are generally not perpendicular to each other because of the spherical boundary between a phase grain and its surrounding. Without loss of generality, the external fields applied to the sample,
Eext and
B , have the directions of the X and Z axes, respectively. Then Equation (33) and Equation (35) lead to
↔σi=σi1+ξ2i(1ξi0−ξi10001+νiξ2i), (38)
where
ξi≡μiB=σiRiB . Analogously we write for
↔σ ,
↔σ=σ1+ξ2(1ξ0−ξ10001+νξ2) (39)
with
ξ≡μB=σRB .
νi=1+cosαi and
ν=1+cosα , where
αi and
α are the angle between
Ei and
B , respectively between
E and
B .
μ and
μi are the average Hall mobility in the composite and the local Hall mobility in the phase i, respectively.
At the interface between a single phase grain and its surrounding continuity of the normal components of the current density and the tangential components of the potential gradient are to be fulfilled. For the limiting case
B=0 , this demand is fulfilled by
f(σ,σi)≡σAσB+σσA(3υA−1)+σσB(3υB−1)−2σ2=0 (40)
following from the EMT formula for
σ , Equation (19).
For the case
B≠0 , the tensor properties of
↔σi and
↔σ , Equation (38) and Equation (39), are to be taken into account. Equation (40) expressed in tensor form reads
↔σA↔σB+↔σ↔σA(3υA−1)+↔σ↔σB(3υB−1)−2↔σ↔σ=0, (41)
where the identities
↔σA↔σB=↔σB↔σA and
↔σ↔σi=↔σi↔σ have been used. Equation (41) determines the coefficients of Equation (39) as a function of the coefficients of Equation (38). Inserting Equation (38) and Equation (39) into Equation (41) and comparing coefficients for the tensor elements, we get
ξ=σAσB(ξA+ξB)+σσAξA(3υA−1)+σσBξB(3υB−1)4σ2−σσA(3υA−1)−σσB(3υB−1), (42)
following from the tensor elements
σxy or
σyx , where quadratic and higher powers of
ξ ,
ξi are neglected, i.e., Equation (42) and the following Equations (43), (44) are low-field approximations. Within this approximation the parameters
νi and
ν do not have an influence on the result. From the tensor elements
σxx ,
σyy , or
σzz , Equation (40) follows.
Substituting
ξ and
ξi in Equation (42) by R and
Ri and considering Equation (40) we get the R formula for two-phase composites:
R=σ2ARA[σB+σ(3υA−1)]+σ2BRB[σA+σ(3υB−1)]σ(σAσB+2σ2) (43)
The same formalism can also be applied to composites with more than two phases leading to relatively complex formulas for R. A self-contained and more manageable description of these R formulas is given by
(Rσ2∂∂σ+∑i=A,B,⋯Riσ2i∂∂σi)f(σ,σi)=0 (44)
with
f(σ,σi)=(∏i=A,B,⋯(σi+2σ))(∑i=A,B,⋯υiσi−σσi+2σ). (45)
Equations (43)-(45) were firstly published by Sonntag [16] .
2.5. Comparison with Other Transport Formulas for Composites
2.5.1. Thermopower
As we are interested in a direct comparison between the different thermopower formulas for composites or heterogeneous materials, the transport coefficients in the phase,
κi ,
Si , and
σi , are set to be constant, although this is not realized in real composites caused by the condition of a common electrochemical potential as well as because of the specific features at
υi<1/3 as discussed in Section 3.2.2. Additionally,
κi=κe,i is set.
The thermopower formula derived by Airapetiants ( [17] , Equation (12) therein), called
SAir , is
SAir=∑iSiσigi1−2σ∑igi (46)
with
gi=3υiκ(2σ+σi)(2κ+κi). (47)
In Figure 1,
SAir vs
υB is shown for a composite consisting of two phases, A and B, with the coefficients,
SA=−18.5 μV/K ,
SB=−370 μV/K ,
σA=7450 Ω−1⋅cm−1 ,
σB=7 Ω−1⋅cm−1 ,
κA=0.022 W⋅K−1⋅cm−1 ,
κB=0.070 W⋅K−1⋅cm−1 which were applied by Airapetiants [17] , already, where
κ and
σ are calculated by Equation (23) and Equation (19), respectively. For comparison the thermopower is calculated with our formula, Equation (31) and Equation (23) for
κ applying the same numbers for
κi and
Si [
σi does not occur in Equation (31)]. The result S is drawn in Figure 1 as well.
SAir(υB) , Equation (46), and our
S(υB) , Equation (31), agree very well in the concentration range
υB<0.62 , however, with increasing
υB beyond
υB=0.62 there is increasing difference between
SAir(υB) and
S(υB) .
Figure 1. Comparison of
SAir ,
SHer , and S vs.
υB calculated by Equation (46), Equation (56), and Equation (31), respectively, where
SA=−18.5 μV/K ,
σA=7450 Ω−1⋅cm−1 ,
κA=0.022 W⋅K−1⋅cm−1 ,
SB=−370 μV/K ,
σB=7 Ω−1⋅cm−1 ,
κB=0.070 W⋅K−1⋅cm−1 [the numbers for
Si ,
σi , and
κi are taken from [17] (Figure 2 therein); the quotients of them,
SB/SA ,
σB/σA , and
κB/κA correspond nearly to those applied in [webman] (Figure 2 therein).
The thermopower formula derived by Webman et al. ([webman], Equations (2.16), (2.17) therein),
SWeb=6κ∑iυiSiDi1−3∑iυiκiDi (48)
1We have replaced the averages
〈…〉 in the original equations by Webman et al. [webman] and Herring [herring] by
∑i υi… , Equations (48) and (55), respectively.
with1
Di=σi(κi+2κ)(σi+2σ), (49)
provides the same concentration dependence as those of Airapetiants [17] ,
SAir(υB) , shown in Figure 1, if
κ and
σ in Equations (48), (49) are again calculated by Equation (23) and Equation (19), respectively. This result seems to be surprizing, because at the end of the derivation of theirs formulas Webman et al. have written: “Our effective-medium result … differs from that previously derived by Airapetiants [17] , as the latter treatment did not involve a completely self-consistent averaging procedure.”
Also the thermopower formula derived by Halpern [19] (Equation (10) therein),
SHalp=SAσAκB−SBσBκAσAκB−σBκA+σAσBκ(SB−SA)σ(σAκB−σBκA), (50)
provides the same concentration dependence as those by Airapetiants [17] ,
SAir ; the same refers also to the thermopower formula derived by Balagurov ( [20] , Equation (7) therein),
SBal=SA+σAσB(SA−SB)σAκB−σBκA(κAσA−κσ) (51)
as well as to those by Bergman and Levy ( [21] , Equation (50) and Equation (9) therein),
SBerg=SB+(SA−SB)γ/σ−γB/σBγA/σA−γB/σB, (52)
with
γ=κ+TσS2, (53)
γi=κi+TσiS2i (54)
if
κ and
σ in Equations (50)-(53) are again calculated by Equation (23) and Equation (19), and
γi and
γ in Equation (52) are replaced by
κi and
κ , respectively, i.e., the second term in Equation (53) and Equation (54) is neglected [22] [23] (quantity of higher order).
Additionally, in Figure 1 it is drawn the concentration dependence of
SHer calculated by the formula by Herring ( [24] , second equation of Equation (38) therein]),
SHer=〈S〉−〈(S−〈S〉)(13κ〈κ〉+23ρ〈ρ〉)〉 (55)
with
ρ=σ−1 , derived for a randomly inhomogeneous medium. If
〈S〉 ,
〈κ〉 and
〈ρ〉 are interpreted as
∑i υiSi ,
∑i υiκi and
∑i υiρi , respectively, then Equation (55) leads to
SHer=∑i υiSi−∑j υj(Sj−∑i υiSi)(13κ∑iυiκi+23ρ∑iυiρi) (56)
with
ρi=σ−1i . It is noteworthy that Bergman and Levy [21] have concluded that for a two-component isotropic composite their result, Equation (52), agree in first order of
σA−σB and
γA−γB with that by Herring [24] , Equation (55), but disagree in the next order.
Fishchuk [25] received for the conducting range of a two-component system in the form of a semiconductor with random insulating inclusions [25] ,
S=SA , in correspondence with
SAir and S in Figure 1.
Summarizing, for the example considered in Figure 1, the surprising result follows: Although the formulas for
SAir,SWeb,SHalp,SBal , and
SBerg look very differently, all provide the same concentration dependence of the thermopower in the complete concentration range
0≤υB≤1 , whereas there are considerable differences to both
SHer(υB) and
S(υB) , Equation (31).
Which are the reasons for the differences between the thermopower formulas considered? All the thermopower formulas cited contain
Si ,
κi and
σi , whereas our Equation (31) contains only
Si and
κe,i (or
Si and
σi , Wiedemann-Franz rule is used, Equation (32)). Also the thermopower formula derived by Xia and Zeng [26] for polycrystals, contains
Si ,
κi and
σi .
The reason for the differences is the neglection of the quadratic term of the thermopower-coefficient in the heat current density (“quantity of higher order”) before
Ji=0 is set. Neglection of this term, “
σiS2i ”, leads to an additional term in the resulting
JS,i , Equation (27), after setting
Ji=0 ; this additional term depends on
σ as well. Such a term does not occur, if all the terms are maintained, before
Ji=0 is set. This statement is independent of the question whether the EMT formula is derived on the basis of J &
JQ or J &
JS .
The starting equations applied by Webman et al. [18] (Equations (2.2b) and (2.2a) therein) read
J=−σgradφ+PgradT, (57)
JQ=−κgradT+PTgradφ, (58)
where
P=Sσ . Setting
J=0 , it follows that
gradT=gradφS, (59)
and replacing
gradT in Equation (58) by Equation (59) we get
JQ=(−κS+σST)gradφ, (60)
Equation (60) contains
σ ,
κ and S. This does not occur if the complete formula,
JQ=−(κ+S2σT)gradT+PTgradφ, (61)
is applied as the starting formula instead of Equation (58). In Equation (28) only
κe and S occur.
Also the other authors have applied analogous formulas: Bergman and Levy ( [21] , Equation (50) and Equation (9) therein); Xia and Zeng [26] (Equation (5) therein). Halpern has also neglected the quadratic term in S (first line in [19] on page L219), although he comes to the reasonable result that for (two-phase) metal-insulator composites,
S≈SA , where
SA is the thermopower of the metallic phase. This result agrees with the result of Sonntag for metal-insulator composites ( [11] , Equation (59) therein).
The basis of all the derivations done in [11] [18] [19] [21] [26] was J &
JQ . As pointed out by Harman and Honig ( [9] , p. 28), for desribing thermoelectric phenomenons,
JS is more universal as
JQ . That is why, in [12] J &
JS has been applied as a basis. There are differences between Equation (31) derived in [12] on basis J &
JS and the formula derived on basis J &
JQ (Equation (30) in [11] ). While in the last one the chemical potentials of the phases,
˜μ0i , occur, in the former one only the common electrochemical potential,
˜μ , occurs mediated through the
Si described by Equation (73).
2.5.2. Hall Coefficient
The EMT formula for the Hall coefficient of two-phase composites derived by Cohen and Jortner [27] is
∑i=A,Bυiσ2iRi−σ2R(σi+2σ)2=0. (62)
As will be argued in the following, Equation (62) seems to be a good approximation for two-phase composites if
σA≈σB , but not if
σA≫σB , as typical for metal-insulator composites.
For three examples of two-phase composites, in Figure 2(b), Figure 2(d), and Figure 2(f), the concentration dependence of R related to its value at
υA=1 is shown, calculated by Equation (43), and compared with Equation (62), denoted as “C&J”. In Figure 2(a), Figure 2(c), and Figure 2(e), the corresponding concentration dependence of the Hall mobility
μ (
=σR ) related to its value at
υA=1 is shown, where
σ is calculated by Equation (19). There are two essential differences between the the two solutions Equation (62) and Equation (43):
(1) The most striking difference appears in Figure 2(a) and Figure 2(c): The “C&J” curves decrease dramatically with increasing
υA and pass through a pronounced minimum at
υA=1/3 , although
μA=μB and
μA>μB , respectively. In contrast, the
μ curves calculated by Equation (43) agree with the expectation: Figure 2(a):
μ agrees with
μA for all
υA ; Figure 2(c) and Figure 2(e):
μ increases and decreases with increasing
υA , respectively.
A possible interpretation for such dramatic decrease of
μ at
υA=1/3 (“C&J” curves) could be additional scattering centres in the added phase boundaries. Such an effect by the phase boundaries is expected to be the more pronounced the smaller the sizes of the phase grains,
Di . However, the C&J formula [27] [28] does not contain
Di .
The differences between Equation (43) and the curves “C&J” increase with increasing difference between
σA and
σB . On the other hand, for the limiting case,
σA=σB , Equation (62) and Equation (43) agree.
(2) Another striking difference between Equation (43) and Equation (62) is represented by the boundary case “
σB=0 and
σA≠0 ”, for which one obtains
Figure 2.
μ/μA and
R/RA versus
υA calculated by Equation (62) (“C&J”) and Equation (43), where
μ/μi=Rσ/Riσi . The “C&J”-curves in (a) and (c) agree with [27] (Figure 1(b), Figure 1(c) therein) and [28] (Figures.13,14 therein), where the same examples are chosen.
1R(43)=1RA(3υA−1)2, (63)
and
1RC&J=1RA(3υA+1)4, (64)
respectively, and for
σ , Equation (19) gives
σ=σA(3υA−1)2. (65)
Starting at
υA=1 , with decreasing
υA both
σ and
1R(43) decrease continuously until they vanish at
υA=1/3 . This result corresponds to the fact that for
υA<1/3 there is no longer a connected metal cluster through the composite (in correspondence with the assumption made earlier that the phase grains are spherical without preferred orientations and arranged in a symmetrical fashion). This result is, however, not reflected by Equation (64)
which gives
1RC&J>0 even for
υA<1/3 , where all the metallic granules are
separated by adjacent insulating phase regions, that is, electron transport through the sample does not happen, if additional tunneling is excluded.
These two differences, (1) and (2), suggest the fact that Equation (63) represents the physical situation better than Equation (64).
2.6. Composites with Non-Spherical Phase Grains
An essential assumption for the derivation of the EMT formulas, Equations (19), (22), (31), (32), (43) and (44) was the fact that the phase grains are spherical. In real composites this assumption is often not fulfilled, especially not for man-made composites with large phase grains. A typical feature for these composites is the fact that the transition from an infinite phase i cluster through the composite does not occur at
υi=1/3 but at larger values called the percolation edge. For such cases McLachlan and co-workers [30] [31] have modifyed the EMT formulas by introduction of percolation elements in the formulas for two-phase composites leading to the generalized form (GEMT)
υAσ1/tA−σ1/tσ1/tA+A⋅σ1/t+υBσ1/tB−σ1/tσ1/tB+A⋅σ1/t=0 (66)
and
υAκ1/tA−κ1/tκ1/tA+A⋅κ1/t+υBκ1/tB−κ1/tκ1/tB+A⋅κ1/t=0 (67)
where A is given by
A=(1−υc)/υc .
υc is the volume fraction of the phase A, where the actual percolation threshold is assumed to occur. t represents the asymmetry of the microstructure. Vaney et al. [29] have applied this idea to the EMT formula for the thermopower. With this modification Equation (32) reads now
υA(σA/SA)1/t−(σ/S)1/t(σA/SA)1/t+A⋅(σ/S)1/t+υB(σB/SB)1/t−(σ/S)1/t(σB/SB)1/t+A⋅(σ/S)1/t≈0, (68)
Vaney et al. [29] have studied real composites composed of a crystalline phase “Bi0.4Sb1.6Te3” (≡phase A) and a glassy phase “Si10As15Te75” (≡phase B) produced by spark plasma sintering (SPS) of powder mixtures with different volume fractions
υA=0.1,⋯,0.5 .
Inserting the
σi and
Si data (given by Vaney et al. [29] ) in Equation (32), it follows the concentration dependence of S on
υA as drawn in Figure 3(a). For
υA>1/3 the calculated S on
υA curve agree with the experimental data. For
υA<1/3 the EMT formulas are no longer a good description, because the
Figure 3. S versus
υA for
[Bi0.4Sb1.6Te3]υA[Si10As15Te75]1−υA composites produced by spark plasma sintering (SPS) [29] and comparison with the experimental data (full points,
T=300 K ). (a), calculated by Equation (32), (b), calculated by the GEMT equation, Equation (68), for
υc=0.16 and
υc=0.10 ;
t=2 for both curves. For the calculations the same values have been applied as given by Vaney et al. [29] :
σ1=0.040 Ω−1⋅m−1 and
SA=165 μV/K for the crystalline phase (Bi0.4Sb1.6Te3),
σ2=2.95−9 Ω−1⋅m−1 and
SB=1400 μV/K for the glassy phase (Si10As15Te75).
phase grains in the man-made composites as studied in [29] are arranged randomly (accidentally), that is, the transition from a continuous A phase cluster through the composite to separated A phase grains does not occur at
υA=1/3 , but it is shifted to smaller values of
υA . These percolation effects are considered by the generalized form (GEMT), Equation (68).
In Figure 3(b), S versus
υA is drawn calculated by Equation (68) for
t=2 and
υc=0.16 (as applied in Figure 7(a) of [29] ) as well as for
t=2 and
υc=0.10 . As can be seen in Figure 3(b), for the parameters
t=2 and
υc=0.10 there is an excellent agreement with the experimental data.
Both in [29] and in Figure 3, the values for
Si and
σi were assumed to be independent of
υi . This is, however, only an approximation, because the electrochemical potential,
˜μ , in the two phases is generally different, as long as the phases are separated from each other. In the composite, however,
˜μ must be uniform. This is ensured by an electron transfer
δn between the phases. The temperature dependent part of
δn leads to a change of
d˜μ/dT . This
provides an additional contribution to the thermopower,
ΔS=1|e|d˜μdT as will be
discussed in Section 3.1. For a calculation of
d˜μ/dT , knowledge of the band structure data of the phases is necessary. For the composite considered by Vaney et al. [29] , the contribution of
d˜μ/dT to S is surely small, because
SA and
SB have the same sign. However, for composites with different signs of
SA and
SB , the effect of
d˜μ/dT can be essentially larger. That is, a really sensitive hardness test of Equation (32), respectively Equation (68), would be a composite with
SA>0 and
SB<0 or vice versa.
3. Formulas for the Transport Coefficients of the Phases
3.1. Composites with Large Phase Grains
“Large phase grains” means that the classical transport theory can be applied separately to each of the phase grains. This fact is immediarely plausible for crystalline phases, but also for amorphous phases as argued in [11] [32] [33] . First, we focus our attention on
σ ,
κe and S. In the framework of the kinetic transport equations, for a zero-magnetic field situation the transport coefficients of the phase i can be written as [9] [10] [14] [15] [34] [35] [36] [37]
σi=e2iK11,i, (69)
κe,i=K31,i−K221,i/K11,iT, (70)
S0i=K21,i/K11,i−˜μ0ieiT. (71)
where
˜μ0i is the chemical potential for the phase i.
S0i is the Seebeck coefficient of the phase i for the (hypothetic) case that electron transfer does not happen between the different phases and that the band edge does not depend on T. Therefore an additional term,
ΔS , is to be introduced considering these effects realizing a common electrochemical potential in the composite,
˜μ=˜μi. (72)
In the composite the Seebeck coefficient of the phase i is given by [12]
Si=S0i+ΔS=S0i+1|e|d˜μdT, (73)
where
S0i is a scattering term.
1|e|d˜μdT is an additional term taking into
account the change of
˜μ with temperature. One consequence of Equation (73) is the fact that the thermopowers of the phases,
Si , dependent on
υi . For a calculation of
d˜μ/dT , knowledge of the band structure data of the phases is necessary. Assuming a two-phase composite consisting of the phase A with electron conductivity and the phase B with hole conductivity,
d˜μ/dT can be calculated by ( [38] , Section 3.3. therein)
d˜μdT=∂EC,A∂T+∂˜μ0A∂T−∂˜μ0A∂T+∂˜μ0B∂T+∂EC,A∂T−∂EV,B∂T1+υA(∂˜μ0B∂p−|e|∂φB∂nB)υB(∂˜μ0A∂n−|e|∂φA∂n), (74)
where
∂EC,A/∂T and
∂EV,B/∂T are the band edge shifts with temperature. n (
≡nA ) is the electron density in the phase A. p (
≡nB ) is the hole density in the phase B.
EC,A and
EV,B are the energies of the band edges of the conduction band (CB) and valence band (VB) in the phase A and phase B, respectively. In Equation (74),
∂φi/∂T=0 ,
∂EC,A/∂n=0 , and
∂EV,B/∂p=0 is assumed. The first assumption corresponds to the fact that the electrostatic potential does not depend on T, the second and third ones are equivalent to the assumption that
EC,A and
EV,B do not depend on occupation of the CB and VB. From Equation (73) and Equation (74) it follows immediately the fact that
Si depends on
υi .
Now,
d˜μ/dT can be calculated by Equation (74) if
∂˜μ0i∂T=−π2k2BT6EF,i, (75)
∂˜μ0i∂ni=2EF,i3ni, (76)
are taken into account and if the band edge shifts for the CB and VB,
∂EC,A/∂T and
∂EV,B/∂T , are known. Equations (75) and (76) follow from the Fermi-Dirac-statistics, where
˜μ0i is given by
˜μ0i=EF,i−π2k2BT26ddE[lnNi(E)]E=EF,i=EF,i−π2k2BT212EF,i (77)
(lowest order in the powers of
kBT/EF,i ) with
EF,i=h28mi(3π)2/3n2/3i, (78)
where
EF,i is the Fermi energy and
mi the effective mass in the phase i, respectively. h is Planck’s constant. The second equation in Equation (77) corresponds to the NFE-approximation.
For metallic phases, Equations (69)-(71) provides (NFE-approximation, [9] [10] [15] )
σi=2(π3)1/3e2hLin2/3i, (79)
κe,i=16π39miLiEF,ih3k2BT, (80)
S0i=π2k2BT(1+ri)3eiEF,i, (81)
ri characterizes the scattering mechanism and represents the energy dependence of the mean free path
Li in the phase i, according to
Li∝Eri .
For composites with semiconducting phases, the
σi ,
κe,i and
S0i are to be calculated according to the rules for semiconducting solids, in correspondence with the two-band (or multiband) model (see, e.g., Harman and Honig [9] , pp. 37 and 129), where additionally the condition Equation (72) is to be taken into account.
For metallic phases, the non-electronic contribution to
κi ,
κne,i , can often be neclected compared to
κe,i , if the carrier densities are not too small. If not,
κ is to be calculated by Equation (23) under consideration of Equation (24), where the
κne,i are to be determined separately; that becomes especially important for semiconducting phases and if the phase i does not form a macro-scopic cluster.
For the Hall coefficient of the phase i in a nonmagnetic composite,
Ri in NFE-approximation is given by
Ri=−C|e|⋅ni=μiσi, (82)
where
μi is the Hall mobility of the phase i and C is a parameter of the order of one depending slightly on the magnetic field. [14] [15] .
The volume fractions of the phases,
υi , can be calculated from the atomic concentrations of the composite and the phases, x and
xi , respectively. For a two-phase composite, the
υi can be determined by
υA=1−υB=[1+NA(x−xA)NB(xB−x)]−1, (83)
where
NA and
NB are the atomic densities in the phases A and B, respectively.
3.2. Nanocomposites
3.2.1. Electron Density in Nanocomposites
Let us consider the class of amorphous transition-metal-metalloid alloys in more detail: it was one of the most important results that amorphous transition-metal-metalloid alloys do generally not consist of randomly mixed atoms, but they form composites consisting of amourphous phases which differ regarding their short range order (SRO). Interpreting a series of transport data of amorphous transition-metal-metalloid alloys, Sonntag [32] has drawn the following conclusions:
For large ranges of concentration there is
(i) amorphous phase separation between two different amorphous phases called phase A and phase B, where each phase has its “own” short-range order (SRO),
(ii) the amorphous phase separation leads to band separation in the conduction band (CB) and valence band (VB) connected with the phases A and B, respectively, and the electrons are freely propagating and the corresponding wave functions are extended with respect to connected phase ranges.
(iii) Between the two coexisting phases there is electron redistribution (electron transfer) which can be described by
n(ζ)=nAexp(−βζ), (84)
where
ζ is the quotient of the volume or atomic fractions of the two coexisting phases.
n(ζ) is the electron density in the phase A with
nA=n(0) .
β is a constant for a given alloy, which is determined by the average potential difference between the two phases.
For a series of amorphous transition-metal-metalloid alloys, conclusion (i) and conclusion (ii) are now confirmed experimentally or supported by independent authors [2] [3] [4] [5] [6] [43] [44] [45] [46] [47] . For details see [33] (Section 1 therein). Although the conclusion (iii) is not yet confirmed or supported by independent authors, there are a series of experimental findings supporting this conclusion (iii), as demonstrated in [32] , where for some amorphous transition-metal-metalloid alloys it is shown that
ρ versus
x/(1−x) follows an exponential concentration dependence (see Figures 1-6 in [32] ). For instance, a-Cr1-xSix thin films measured by Helms et al. [39] and shown in Figure 4 follow an exponential dependence of
ρ versus
x/(1−x) . As argued in [32] , one can assume that for
x>0.25 a-Cr1-xSix thin films consist of the two amorphous phases a-Cr3Si (=phase A) and a-Si (=phase B), i.e., the silicon concentrations in the two phases are characterized by
xA=0.25 and
xB=1 , respectively, whereas for
0≤x<0.25 , the films are one-phase. In the two-phase range,
ζ (
=XB/XA ) can be calculated by
ζ=XBXA=x−xAxB−x (85)
As can be seen in Figure 4(b) (lower grafic),
ρ follows an exponential dependence on
ζ (
=XB/XA ). This result is in correspondence with Equation (84) if the effect of
σB (about Equation (19)) as well as the effect of a variation of the electronic mean free path L can be neglected. Because the scattering in a-Cr1-xSix is very strong, L is very small. However, L cannot be smaller than the average atomic distance, d; therefore one can assume that it is nearly independent of x, that is,
d≈L and
σA (
=1/ρA ) is expected to be directly proportional to the electron density, n. The effect of the second phase (a-Si = phase B) on
ρ (
=1/σ ) (about Equation (19)) can be neglected, because in a-Cr1−xSix the concentration dependence of
σB is similar to
σA . This conclusion follows from calculations done in [11] (see Figure 4(c) therein).
Support for Equation (84) comes also from Hall coefficient measurements of metal-insulator composites as shown in Figure 5. The sizes of the phase grains are of the order of some atomic distances as well. In Figure 5(a) and Figure 5(b) the experimental original data of Cu1-y(SiO2)y and Ni1-y(SiO2)y taken from [40]
Figure 4. Experimental data
ρ versus
x/(1−x) (upper diagram) and
ρ versus
XB/XA (lower diagram) for a-Cr1-xSix thin films at
T=300 K taken from Helms et al. [39] .
XB/XA is calculated by
XB/XA=(x−xA)/(xB−x) with
xA=0.25 and
xB=1.00 .
XB and
XA are the atomic fractions of the amorphous phases a-Cr3Si (=phase A) and a-Si (=phase B).
and [41] are drawn versus
y/(1−y) , and in Figure 5(c) and Figure 5(d) the corresponding
RA drawings calculated by Equation (63) (
σB=0 , because SiO2 is an insulator). The exponential concentration dependence of
RA in Figure 5(c) and Figure 5(d) reflects immediately Equation (84) because of Equation (82).
R and
RA versus
y/(1−y) follow an exponential dependence also for Ni1-y(SiO2)y, although we had assumed non-magnetic materials for the derivation of Equation (43), wherease Ni is a magnetic material. This is however not a
Figure 5. Experimental Hall coefficient data at 5 K versus
y/(1−y) for Cu1-y(SiO2)y, (a), (c), and Ni1-y (SiO2)y, (b), (d), taken from [40] (circles), [41] (triangles) and [42] (diamonds). (c), (d):
RA calculated by
RA=R(3υA−1)2 according to Equation (43), where
υA=1−y is set.
discrepance. For magnetic metal-insulator composites Equation (82) holds approximately if “=” is replaced by “
∝ ” considering the effect of the additional internal magnetic field due to the magnetization: An electron sees the effective magnet field
Hw=H+Hi , where
Hi≫H . H is the external field applied to the specimen and
Hi is the internal field produced by the quantum mechanical exchange forces ( [48] , p. 341). An electron does not distinguish between H and
Hi . It moves according to the Lorentz force determined by
Hw and the electric field E. One can assume that
Hw is nearly proportional to H as long as
Hi is nearly proportional to the magnetization produced by H.
This assumption is supported by the experimental finding by Xiong et al. [49] that (for not too small fields H), in the granular Co-Ag system, the Hall resistivity
ρxy is linearly proportional to H. If so, the measured R values differ from the calculated R values, Equation (82), only by a factor which is nearly constant. Therefore, we assume that the EMT formula for R, Equation (43), can be applied to magnetic composites as well.
If the metallic phase of a M-I composite is a noble metal, the NFE-approximation is a good one for the metallic phase, above all as the Fermi surface moves away from the Brillouin zone boundary as n decreases. For the metallic phase in Ni-SiO2 the NFE-approximation is surely also a good one, because Ni has only 0.55 4s valence electrons per Ni atom ( [48] , p. 271).
Additional support for Equation (84) comes from experimental resistivity data of amorphous (Cr1-xSix)1-zNz and (Cr1-xSix)1-yOy thin films (in preparation).
1) Amorphous Metals
The electron transfer
δn between the phases described by Equation (84) leads to a lowering of the total energy of the composite compared with a situation, where the phases exist alone. This is the reason for the fact that metallic composites with an amorphous structure can exist at all. Almost invariably, amorphous metals contain a metalloid as one of the constituents. e.g., Au-Si, Pd-Si, Fe-P-C, … [50] . That is, one-phase amorphous metals practically do not exist. If yet (for instance after an evaporation process at extremely low temperatures), then there is the strong suspicion that during the preparation process impurities as oxygen and nitrogen are incorporated in the metal film.
Of course, for the crystalline state the energy gain is surely larger compared with that of the amorphous state. However, the transition from the amorphous state to the crystalline one realized by atomic diffusions processes requires additional energy to overcome energy barriers.
2) The Giant Hall Effect (GHE)
As reasons for the very large values of the Hall coefficient in metal-insulator composites as shown in Figure 5 for the examples Cu1-y(SiO2)y and Ni1-y(SiO2)y (known as the “Giant Hall effect”), quantum size effects and quantum interference effects on the mesoscopic scale have been discussed. [40] [51] [52] [53] [54] [55] However, to our knowledge, until now, there is no explanatory model which can interpret the phenomenon of GHE. Now, we believe that Equation (84) gives a simple and plausible explanation for the GHE.
3.2.2. Electronic Transport in Nanocomposites
For many alloys with phase separation the phase grains are very small so that an application of the classical transport equations to the phase grains (Section 3.1) does not seem to be appropriate. In spite of this reservation, both the BTE and the approximation of free electrons (NFE-approximation) may be good descriptions for a phase i of a composite, as long as it forms an infinite cluster in the composite and the scattering processes are elastic. For amorphous transition-metal-metalloid alloys this condition is fulfilled for
υi>1/3 . This point of view is justified in [33] (Section 4.2 therein) and [12] (Section 5.2. therein). That is, formulas for the transport coefficients of the phases, Equations (79)-(82), may be a good description if the phase i forms a macroscopic cluster through the composite corresponding to
υi>1/3 . For this case the phase i can form a quasi-continuous energy band and the scattering is elastic. With decreasing
υi , the macroscopic i phase cluster decomposes into separate phase grains; for amorphous composites, this transition occurs very precisely at
υi=1/3 . The reasons are the following [12] [33] :
(1) As the grain diameters are very small (of the order of ~1 - 2 nm [5] [6] , the number of grains (of the same sort i) must be very large to form a cluster for which a quasi-continuous energy band is realized. Accidentially formed agglomerates of some (for instance 10 or 100 or 1000) neighbouring grains (with ~1 - 2 nm diameter) are not sufficient for forming a quasi-continuous energy band.
(2) During the film deposition of a composite, the atoms of the different atom sorts arrive at the substrate equally distributed; therefore the different phase grains (A and B) can also be assumed to be locally equally distributed in the amorphous composite, because the diffusion paths during solidification are very short, which is a prerequisite for forming amorphous composites.
A phase is an “electronic phase” determined by a solution of the Schrödinger equation; after hitting at the substrate, the atoms move locally only so long until they can form a phase which corresponds to a solution of the Schrödinger equation. That is why, the phase grains of the same sort i are also locally equally distributed, as the compositions of the different phases are very different, i.e., the local distribution of the i phase grains is not completely random, as, e.g., assumed within the framework of the classical percolation theory. For such a locally equally distributed arrangement of the i phase grains in an amorphous matrix (formed by the rest of the composite), it follows that this merging to a macroscopic cluster through the sample occurs very precisely at a specific concentration
υi ; and this specific concentration is
υi=1/3 , as follows, for instance, from Equation (19) setting
σj=0 , but
σi≠0 (
j≠i ). (For a comparison with the classical percolation theory concerning this critical value of
υi , see also [33] (Section 4.1 therein).
(3) A macroscopic i phase cluster is only realized, if all the atoms in this i phase cluster are directly connected to atoms belonging to the same phase sort i. When, for instance, two grains of the same phase sort i are separated by a monoatomic layer of a different(!) phase (j), these two i phase grains cannot be considered as (nearly) one i phase cluster, because the overlap of the electron wave functions is interrupted by this monoatomic layer. Within an i phase grain or i phase cluster the wave functions,
ψi(r) , overlap, but they do not overlap between two i phase grains or i phase clusters separated by a monoatomic layer of a different(!) phase (j).
If a sufficiently large number of i phase grains form a macroscopic cluster, the overlapping wave functions
ψi(r) form a quasi-continuous energy band, while the wave functions fall off exponentially in a very short distance
|δr| outside this macroscopic cluster. This “falling off” is comparable with the decrease of the molecular orbitals of (large) molecules at their molecule boundaries (see, e.g., [56] , pp. 409 and pp. 435, or [57] ). Because of this “non-overlapping” of the
ψ(r) belonging to different phases, they do not “penetrate” each other, i.e., there are sharp boundaries between the different (electronic) phases.
Considering the fact that the phase i does not form a quasi-continuous energy spectrum for
υi<1/3 , but there is a discrete energy spectrum typical for separate grains, then there are no electronic states immediately above and below
˜μ (within the energy range
kBT ), i.e., the electrons cannot be activated to higher energies (at the hot end of the sample) and cannot deliver energy (at the colder end), if
υi<1/3 . Under this condition it follows that (at the temperature T)
κe,i=0 for
υi<1/3 . For
Si the situation is analogous: for
υi<1/3 ,
S0i=0 , i.e., at the transition from
υi>1/3 to
υi<1/3 , both
κe,i and
Si change discontinuously. Such discontinuities are especially to be expected in composites with metallic phases.
For amorphous transition-metal-metalloid alloys, “
κe,i=0 ” and “
S0i=0 ” for
υi<1/3 is only an approximation, because, at the boundary faces between the different phases, there are p-d bonds, i.e., d orbitals of the Cr atoms (of the A phase grains) overlap with p orbitals of the boundary faces atoms on the B phase grains resulting to a p-d band, which is incompletely occupied. (For a detailed discussion see Section 2.1 of [33] ). Because of this p-d overlapping, this (incompletely occupied) p-d band exists also if
υA<1/3 : although the energy levels of the s states in the separate A phase grains are discret, the electrons can be activated to higher energies (at the hot end of the sample) and deliver energy (at the colder end) by a transition to the p-d band. That is why,
κe,A≠0 and
S0A≠0 , also for
υA<1/3 . However, the real values of
κe,A and
S0A (for
υA<1/3 ) are expected to be essentially smaller than calculated by Equations (80) and (81), because the d states are localized.
4. Discontinuity in the Concentration Dependence of the Thermopower
A specific feature of composites with
SA>0 and
SB<0 or vice versa is the fact that a discontinuity (step) in the calculated S vs.
υi can occur, which is an additional possibility to check experimentally the thermopower formula. One example for a composite with different signs of
SA and
SB is a-Cr1-xSix consisting of the two amophous phases [33]
a-Cr1−xASixA and
a-Cr1−xBSixB with
xA≈0.25 and
xB≈1 . (Amorphous phase separation is typical for many amorphous transition-metal-metalloid alloys [2] [3] [4] [5] [6] ). In [12] S versus x had been calculated for a-Cr1-xSix applying Equation (31) for
xA=0.25 and
xB=0.90 (bold line in Figure 1 of [12] ).
Now, the calculations of [12] we have done once more, but for
xA=0.25 and
xB=1.00 shown in Figure 6 commonly with the result of the calculations for
xA=0.25 and
xB=0.90 . Both the discontinuity (step) at
x=0.49 and the general trend of the calculated curves agree with the experimental data.
The discontinuity of S versus
υi for composites with
SA>0 and
SB<0 or vice versa has its origin in the mathematic structure of the formula
S(±)=4⋅κeκe,ASA(3υA−1)+κe,BSB(3υB−1)±√(κe,ASA(3υA−1)+κe,BSB(3υB−1))2+8κe,A⋅κe,BSA⋅SB
(86)
following from Equation (31) for two-phase composites, where
υB=1−υA . Equation (86) has two solutions,
S(−) and
S(+) , which both show a discontinuity (step) at the same concentration, when
S(−) and
S(+) passes the value “0” coming from negative values crossing to positive values or vice versa. The physics follows only one of them,
S(−) , as suggested by the results of [11] (compare Figure 4(g) and Figure 4(h) with Figure 5 therein).
As this discontinuity occurs at
S=0 , this phenomenon opens the possibility
Figure 6. Thermopower versus x for a-Cr1-xSix at
T=300 K calculated by Equation (32) with
xA=0.25 and
xB=0.90 (bold line) and
xA=0.25 and
xB=1.00 (dotted line) and comparison with experimental data. Details of the calculations and experimental data are described in [11] (section IVB therein).
to produce reference standards for absolute thermopower
S=0 even for
T>Tc (
Tc stands for the transition temperature of any superconductor).
5. Extension of the Classical Thermopower Formula
In the limit
υB=0 , the composite degenerates to a homogeneous alloy consisting exclusively of the phase A. On the opposite side, for
υB=1 we get a homogeneous alloy consisting exclusively of the phase B. For these two limiting cases the formulae must hold as well. Setting
υB=0 , it follows from Equation (74) that
d˜μdT=∂EC,A∂T+∂˜μ0A∂T, (87)
and taking into account Equations (73), (81), (68) and Equation (22),
S(+)=SA=π2k2BT(1+rA)3eAEF,A+1|e|(∂EC,A∂T+∂˜μ0A∂T). (88)
Analogously it follows for
υB=1 that
d˜μdT=∂EV,B∂T−∂˜μ0B∂T, (89)
S(+)=SB=π2k2BT(1+rB)3eBEF,B+1|e|(∂EV,B∂T−∂˜μ0B∂T). (90)
Inserting Equation (75) in Equations (88), (90) we get
S=−π2k2BT(z+r)3|e|EF+1|e|dECdT, (91)
S=π2k2BT(z+r)3|e|EF+1|e|dEVdT, (92)
with
z=3/2 , where the index i is omitted and
∂EC,A/∂T and
∂EV,B/∂T are replaced by
dEC/dT and
dEV/dT , respectively, because a homogeneous metal consists only of one phase. Each of the two limiting cases, Equation (91) and Equation (92), represents a homogeneous metal with electron conductivity and hole conductivity, respectively, and constant carrier density, i.e., interband transfer of electrons (in dependence on temperature) are not considered.
The term “
∂˜μ0i/∂T ” reflects the lowering of the chemical potential with T described by the Fermi-Dirac-statistics. In Equations (87)-(92) it is assumed that Equation (81) does exclusively represent the scattering part of S. It is, however, not completely clear, whether Equation (81) does indirectly contain this term “
∂˜μ0i/∂T ”, already. If so, then in Equations (88), (90)-(92) the term “
∂˜μ0i/∂T ”
(respectively “
−π2k2BT6EF,i ”) is to be deleted, and the value for z is to be replaced by
z=1 .
Equation (81) agrees with the first term in Equation (91) if
z=1 is set. It represents the contribution of the scatterring on the thermopower, whereas the second term in Equation (91) represents the effect of the temperature dependance of the band edge on the thermopower. Therefore, we can say that for normal metals positive sign of thermopower will be measured if
dEC/dT>0 and if this effect overcompensates for the influence of the first term in Equation (91). This conclusion holds exactly if
z=1 . If
z=3/2 , this fact is to be considered as a tendency.
As mentioned earlier, it is not yet completely clear whether
z=3/2 or
z=1 . We believe that for metals
z=3/2 is correct. However, the final answer depends on the question whether or not Equation (81) contains exclusively the scattering contribution [corresponding to
S0 , Equation (81)]. This question is a matter of future studies.
Simple Metals with Positive Thermopower
It is a widely forgotten mystery, why the thermopower of a metal as simple as lithium is positive [10] [58] [59] [60] [61] . Considering only the scatterring term of Equation (91) (as it was classical standard theory according to Equation (81)), positive sign was incomprehensible. However, considering the second term of Equation (91), this puzzle is now solved. Positive sign of thermopower can also be expected for crystalline metals: As a trend, positive sign of thermopower will be measured, if the conduction band edge shifts to higher values with increasing T and if this effect overcompensates for the influence the scatterring term. Thus, for the crystalline Cu, Ag, Au, and Li the experimental thermopowers also have positive sign; at 0˚C they are +1.7, +1.4, +1.1, and +11.5 μV/K, respectively ( [48] , p. 316). For more details see [38] .
6. Constraints of the Combinatorial Development of Materials
The combinatorial development of materials is a very effective method to get experimental data about a certain material system, because a large range of different compositions can be realized concurrently on one large substrat, for instance by deposition of thin films by simultaneous co-sputtering from two or three targets on a large substrat. On the other hand, this method is to be used with great caution because the results can be completely different from a situation where the different compositions are produced separated in single manufacturing processes. The reasons are the following ones.
In a composite a common electrochemical potential is realized by electron transfer between the different phases. This can lead to the fact that in one of the phases another topological structure (atomic configuration) is more favorable than for the case that this phase exists alone (as homogeneous material). And in nanocomposites the electron transfer has an essential larger effect on the carrier densities and therefore on the transport coefficients in the phases compared with composites with large phase grains.
Figure 7 shows two series of resistivity measurements of Cr1-xSix thin films
Figure 7. Experimental resistivity data of Cr1-xSix thin films cosputtered from a chromium target and a silicon target. By variation of the sputtering rates two different concentration ranges about 4 inch glass wafers are realized, 1)
0<x<0.30 , and 2)
0.22<x<0.80 . Details of film deposition, analysis and resistivity measurements are described in [62] . Additionally the resistivity data of co-sputtered Cr1-xSix thin films by Helms et al. [39] are shown (circles; already shown as a logarithmic representation in Figure 4), where the single samples were separated from each other during the deposition process.
cosputtered from a chromium target and a silicon target and deposited on 4 inch glass wafers (the deposition conditions are specified in [62] [63] ). Three parallely lying measuring lines have been measured with a distance of 4.5 mm to each other. The first series extend continuously about the concentration range
0<x<0.30 , the second one about
0.22<x<0.80 . The corresponding film thickness were 200 nm > d > 80 nm and 165 nm > d > 131 nm, respectively. Additionally resistivity data received from Helms et al. [39] are drawn. These data were produced by co-sputtering on glass as well, but with the difference that each measuring point represent a single sample; these single samples were separated from each other during the deposition process, i.e., there were no electrical connections between the single samples. The structure of all the films were detected to be X-ray amorphous, where for the chromium rich samples microcristalline precipitations has been found; number and average size of these precipitations decrease with the silicon content.
As can be seen in Figure 7, in the overlapping range between the two series “
0<x<0.30 ” and “
0.22<x<0.80 ”, there are very large
ρ differences, nearly a factor of 1.7 … 1.8. The factor between the “
0<x<0.30 ” series and the Helms samples is even 3.0 … 3.1.
Figure 8 shows experimental thermopower data of the same samples as shown in Figure 7 supplemented by a third series for the concentration range
0.49<x<0.96 with the corresponding film thickness of 70 nm < d < 260 nm [62] . The structures were detected to be X-ray amorphous. S has been measured applying the Potential-Seebeck Microprobe (PSM) from the company LOT-Oriel Group Europe [64] [65] .
The S curves for the concentration ranges
0<x<0.30 and
0.22<x<0.80 as well as for the single samples (triangles) correspond relatively well. However, between the second series and the third one there are large differences. The discontinuity at
S=0 occurring in the S curve of the single samples does not occur in the co-sputtered samples.
Figure 8. Experimental thermopower data of Cr1-xSix thin films for the same samples as shown in Figure 7 supplemented by a third series
0.49<x<0.96 . Details of film deposition, analysis and thermopower measurements are described in [62] . Additionally it is drawn the experimental S data of the Cr1-xSix thin films shown in Figure 6 already.
Figure 9 shows thermopower data of a-(Cr1-xSix)1-yOy thin films co-sputtered as described in [62] for Cr1-xSix thin films, where additionally an oxygen gradient was realized, lateral and perpendicular to the wafer. S has been measured at the complete wafer and after that the wafer was sawed producing single 8 × 2 mm2 pieces. Now S was measured once more at these pieces. The S vs.
x/(1−x) curves for the “wafer” and the “pieces” are shifted nearly parallely to each other, where the zero-crossing of S is shifted from
x=0.31 to
x=0.59 (corresponding to
x/(1−x)=0.45 to
x/(1−x)=1.45 .
The experimenal results shown in Figures 7-9 demonstrate strikingly that the combinatorial development of materials can provide results which can be different from samples produced by single procedures with a given composition.
7. Conclusions
Formulas have been presented for calculation of
σ ,
κ , S, and R in composites. We have shown that these formulas can also be applied successfully to nanocomposites if certain conditions are taken into account, especially the
Figure 9. S vs.
x/(1−x) of a a-Cr1-xSix-O thin film with an additional O gradient vertical to the wafer surface as well as in the same direction as those of silicon concentration gradient (measured for the complete wafer and the pieces after sawing to 8 × 2 mm2 pieces). The measuring positions with respect to the wafer, before and after sawing, were the same ones. This figure is taken from [62] .
phenomenon of electron transfer between the different phases in the nano-composite. It is argued that, under certain conditions, an noncritically application of the combinatorial development of materials can provide unfeasible results.
The formulas and the theory described in the present article give answers to some mysterious puzzles, for which the scientific research had no final answers:
1) Why there are simple metals with positive thermopower?
2) What is the reason for the phenomenon of the “Giant Hall effect”?
3) What is the reason for the fact that amorphous metals can exist at all?
4) Until to the end of the twentieth century amorphous metallic alloys were assumed to be a random and homogeneous distribution of the metal atoms in the amorphous matrix. On this basis, a quantitative calculation applying classic theories was not successful.
The answers given in the present article are the following:
1) The reason for positive thermopower of some metals comes from the temperature dependance of the band edge: If
dEC/dT>0 and if this effect overcompensates for the influence of the scattering term in the thermopower formula, Equation (91), then the thermopower is positive.
2) The reason for the phenomenon of the “Giant Hall effect” in metal-insulator composites is the exponential reduction of the electron density in the metallic phase due to electron transfer from the metallic phase in direction to the insulating phase, described by Equation (84). The transferred electrons are pinned at the phase boundaries between the phases.
3) The reason for the fact that amorphous metals can exist is an electron transfer between the phases described by Equation (84). This electron transfer leads to a lowering of the total energy of the composite compared with a situation, where the phases exist alone. (Amorphous metals are generally composites.)
4) Experimental and theoretical studies at amorphous transition-metal-metalloid alloys have shown that in these amorphous alloys there exists amorphous nano-scaled phase separation between two different amorphous phases, where each phase has its “own” short-range order (SRO) and each phase may be described by its “own” band structure. Within this framework the formulas for
σ ,
κ , S, and R described in this paper can be applyed.
Acknowledgements
The authors thank Mr. Alan Savan from the Institute of Materials, Faculty of Mechanical Engineering of the Ruhr-University Bochum, and Mr. Michael Kieschnik from RUBION, Ruhr-University Bochum, for co-sputtering the Cr-Si thin films and the RBS analysis, respectively, shown in Figures 7-9.