Closure relations problem of hydrodynamical models in semiconductors is considered by expressing third- and fourth-order closure relations for the moments of the distribution function
in terms of second-order Lagrange multipliers using a generalized Maxwell-Boltzmann distribution
function within information theory. Calculation results are commented and compared with
others to justify the accuracy of the approach developed in this paper. The comparison involves,
in the first part with good agreements, the closure relations results obtained within extended
thermodynamics which were checked by means of Monte Carlo simulations, in the second part,
the results obtained by Grad's method which expands the distribution function up to fourth-order
in Hermite polynomials. It is seen that the latter method cannot give any restriction on closure
relations for higher-order moments, within the same conditions proposed in our approach. The
important role of Lagrange multipliers for the determination of all closure relations is asserted.
1. Introduction
The analysis of transport in small semiconductor devices is essential for the optimization of their functioning. Such transport could in principle be described by means of Boltzmann transport equation (BTE) for charge carriers. However, in small devices the electric fields are extremely large, and therefore nonlinear effects are unavoidable [1, 2] which leads to insurmountable difficulties to obtain solutions. Notwithstanding this, BTE contains more information than needed in practical applications. It is common in practice to consider only the lowest-order moments of the distribution function, which are directly related to density, charge flux, kinetic energy, heat flux, and so on. These variables are measured and controlled. This kind of approach is called a hydrodynamical approach [3–5].
The basic model, in which the various steps and approximations are derived and discussed in detail, is due to Blotekjaer [4]. So as to close the set of balance equations considered by Blotekjaer, one assumes that higher-order moments have the value appropriate for a displaced Maxwellian. A slightly different model has been suggested by Hänsch and Miura-Mattausch [6]. In their model, the distribution function is expanded in Legendre polynomials and only the first two terms in the expansion are retained. Only the five balance equations for particle number, momentum, and energy are considered; then the closure is accomplished by means of the Wiedemann-Franz law for heat flux. Both models [4, 6] are then further simplified in order to provide a manageable set of equations appropriate for the device simulation. However, for more accurate results, the full models must be registered. In approaches similar to those cited in references [4, 6], Woolard et al. [7] and Thoma et al. [8] have proposed models taking into account the nonparabolicity of the band structure of the crystal. All these approaches have in common the assumption that some higher-order moments can be calculated by means of a displaced Maxwellian, such an approximation is rather rough and imprecise and its range of validity needs to be assessed [9].
Another method suitable for deriving hydrodynamical-like equations is Grad’s method of moments [10]. This method yields, with an appropriate truncation, a set of evolution equations for the thirteen fields comprising, beside the five balance laws corresponding to particle number, momentum, and energy, rate-type equations for heat flux and anisotropic stress. These equations are known to describe dilute gases only near thermal equilibrium and fail in nonequilibrium situations. Note that Banach and Larecki [11] proposed a generalized Grad’s method, which begins by expanding the phase density about an anisotropic Planck function, in order to close systems of moment equations in nonequilibrium situations.
Extended irreversible thermodynamics [1, 2] and extended thermodynamics [12, 13] are relatively recent approaches to nonequilibrium thermodynamics phenomena [14, 15], which at variance with classical irreversible thermodynamics [16] incorporate higher-order moments (to be interpreted as flux of fluxes) in the thermodynamic state description of system. Though both approaches coincide at least up to second-order in the fluxes, we adopt in this paper extended thermodynamics in order to derive hydrodynamical equations for carrier transport in semiconductors and try to close the system by means of a generalized distribution function within information theory (IT) [17]. Furthermore, our approach deals with the higher-order moments based models which still require a lot of fine tuning and a detailed understanding of the underlying physical phenomena as mentioned by Grasser et al. [18]. Furthermore, two recent contributions based on the maximum entropy principle (MEP) and validated through numerical results are to underline at this stage; the first one concerns an exact closure obtained of the 8-moment model for silicon semiconductors [19], and, the second one is a technique for the construction of realizable 5-moment closures whose extension to a fully three-dimensional gas was not necessarily simple. In fact, the cost of the accurate numerical integration of multidimensional distribution functions required for the resynchronization step was expected to be overwhelming [20].
Some rudimentary results of our approach were published previously [21] which contain some difficulties in comparing with other works based on Monte Carlo simulations, specially [22]. Here, we revisit it completely to overcome difficulties and to provide accurate results. It will be done by assuming nonvanishing higher-order Lagrange multipliers (L-Multipliers). Such hypothesis is not the first one to this field but there is some relevant contributions in quantum regime due to Trovato and Reggiani [23–26], in which they used quantum maximum entropy principle (QMEP) and dealt with closure problems by assuming that L-Multipliers could be expanded in powers of ℏ2; ℏ being the reduced Planck constant, in order to determine the reduced Wigner function [27] for equilibrium and nonequilibrium conditions. Thus, we compare and comment briefly, the aims and the main results of the recent unified physical approach [25] with those of the present paper. Finally, an attempt to solve the closure problem by means of Grad’s method will show within the same conditions proposed to deal with closure relations (CR) in our approach.
The plan of this paper is as follows: in Section 2, we recall the basic formalism of the hydrodynamic description in order to write the CR, that is, third- and fourth-order moment tensors in terms of lower-order quantities. In Section 3, we summarize some foregoing works dealing with the CR problem. Section 4 presents the treatment of the closure problem, in the case of an ideal gas under heat flux or subject to both heat flux and viscous pressure by expanding a generalized distribution function up to second-order. Furthermore, we compare our results with other works. In Section 5, we summarize for comparison some results of the unified physical approach underlined above, particularly those concerning nonlinear CR. Section 6 is devoted to concluding remarks. The detail of some expressions used in this paper and some useful integrals are included as Appendices A, B, and C.
2. Hydrodynamic Description
The evolution equations for the moments of the distribution function are directly obtained from the semiclassical BTE for charge carriers [3]. In the case of electrons in the conduction band of a semiconductor it is written as
(1)∂f∂t+v(k)·∇f-eE·∇kf=P,
with f(x,k,t) the distribution function, k the electron momentum, v(k) the electron group velocity given byv(k)=∇kε, where ε is the energy, e the absolute value of electron charge, E the electric field, and P the collision term.
In the effective mass approximation [28], the energy is given by
(2)ε(k)=k22m*withv(k)=km*,
where m* is the effective electron mass which in silicon is m*=0.26me with me the electron mass and the reduced Planck’s constant is taken as ℏ=1 for convenience.
By multiplying (1) by several products of the components of the momentum vector k and integrating, we obtain a hierarchy of equations for the different moments of the distribution function. Thus, we obtain different hydrodynamical models which differ in the choice of selected variables and on the level at which the truncation is made in the hierarchy [1, 2]. We will summarize these models in Section 3. For instance, let us write this hierarchy of equations [9, 14, 15, 29].
For consistency, the boundary of the first Brillouin zone is moved toward infinity, the particle density n(x,t) and the mean velocity u(x,t) are defined, respectively, as n(x,t)=∫dkf(x,t,k) and u(x,t)=∫dkv(k)f(x,t,k), where J=nu defines the particle flux. Assuming as usual that f(x,t,k) vanishes sufficiently fast at infinity and integrating (1) in k-space, we obtain the particle continuity equation
(3)∂n∂t+∇·(nu)=0.
Multiplying (1) by ki and integrating yield the momentum balance equation
(4)∂(nui)∂t+∂θij∂xj+neEim*=Pi.
Multiplying (1) by kikj and integrating, we find the expression of stress equation for the traceless part of θ^〈ij〉(5)∂θ^〈ij〉∂t+∂θ^〈ij〉r∂xr+2neE〈iuj〉m*=P〈ij〉.
We denote the completely symmetric and traceless part of a tensor Aij as
(6)A〈ij〉=12[Aij+Aji-23δijAkk].
From the trace of θ^ij, we obtain the balance equation for the energy W, with PW the energy production,
(7)∂W∂t+∇·S+neE·u=PW.
The energy flux equation for Si is governed by the following equation, with Pi′ as the energy-flux production:
(8)∂Si∂t+∂Sij∂xj+e[Eiθ^ij+1m*WEi]=Pi′.
Here Sij denotes the flux of energy flux which is written as
(9)Sij=Wuiuj+2nkBTuiuj+12u2nkBTδij+m*2(2urujθ^〈ir〉+2uiurθ^〈jr〉+u2θ^〈ij〉)+ujQi+uiQj+urθ^ijr+θ^ijrr,
with Qi the heat flux. θ^ijr and θ^ijrr are, respectively, the random parts of the third- and fourth-order moment tensors, (10).
3. Closure Problem and Extended Thermodynamics
As we have mentioned above, different hydrodynamical models have tried to close the hierarchy of equations for the different moments of the distribution function. The variables appearing in the thirteen-moment equations (3)–(8) are the particle number density n, the momentum nui, the temperature T, the stress tensor θ^〈ij〉, the heat flux Qi, and plus the random parts of the third- and fourth-order moment tensors θ^ijr and θ^ijrr given by
(10)θ^ijr=m*2∫fcicjcrdk,θ^ijrs=m*2∫fcicjcrcsdk,
where the random component c of the moment vector k is used for convenience and which is defined by k=m*(u+c) with u the mean velocity.
So as to close the system at this level of approximation, we need explicit expressions for θ^ijr and θ^ijrr. One way is based on the assumption that the distribution function could be approximated by a drifted-Maxwellian [3, 4, 8, 30]. This approach implies the following closure relations for the random parts of the third-order and fourth-order moment tensors, (10),
(11)θ^ijr=15(Qiδjr+Qjδir+Qkδij),(12)θ^ijrr=52n(kBT)2m*δij.
However, this model presents some limitations which makes a sounder analysis of the possible truncation of equations to optimize the description of the system desirable [1, 2]. Therefore, CR related to the hierarchy of moment equations of the BTE take an essential role in building hydrodynamical models of carrier transport in semiconductors, and it deserves a thorough investigation starting from first principles [12, 13]. This is done in framework of extended thermodynamics [1, 9, 12–15, 29]. The critical assumption is that θ^ijr and θ^ijrr can be considered as functions of the lower-order moments of particle density, momentum, energy, components of the pressure tensor, and heat flux. Within the definition and the justification of the partial thermal equilibrium state, Anile and Muscato [29] derived the results (11) and (12) as θ^〈ijr〉=0 and θ^ijrr=(5/2)((n(kBT)2)/m*)δij which were obtained with the drift-Maxwellian assumption. The last closure led to a flux-limited expression for the heat flow by using the Schwartz’s inequality as
(13)|Q|≤Qmax=3(52)1/2nm*(kBT)3/2.
Let us mention that other attempts to investigate the flux-limited phenomena are well known, for instance, in radiations hydrodynamics [9, 31, 32], in plasma physics [33, 34] and in electronics [1, 2, 29, 35–38]. Furthermore, Anile et al. have applied [9, 14, 15, 29, 39, 40] a method based on extended thermodynamics [1, 2, 13] in order to establish expressions corresponding to third- and fourth-order moment tensors, that is, θ^ijr and θ^ijrr, in the second-order of Qi and θ^ij. They obtained the following constitutive equations up to second-order about partial thermal equilibrium:
(14)θ^〈ijk〉=25m*nkBT(Qiθ^jk+Qjθ^ki+Qkθ^ij215iiiiiiiiiiiiiiiiiii-215Ql(θ^〈li〉δjk+θ^〈lj〉δki+θ^〈lk〉δij)),(15)θ^illj=(5n(kBT)22m*+(σMkBTL-25nkBT)Q2)δij+2σnkBTLQiQj+7kBT2θ^〈ij〉+m*nθ^〈il〉θ^〈jl〉,
where TL is the lattice temperature.
In the constitutive equation (15), only the free parameter σ can be determined by comparison with Monte Carlo results [22]. In order to check CR, (14) and (15), the authors evaluate the quantities n, ui, T, Qi, and θ^〈ij〉 by Monte Carlo simulations and by using a suitable fitting. Accordingly, they determine the free parameter σ appearing in (15). In fact, this free parameter σ is in the range −0.15–+0.15 and it gives a maximal error less than 7% with respect to the Monte Carlo calculations as the authors described. Furthermore, one finds that in (14), the error increases with electrical field E. Indeed, for E≃104 V/cm the error is 0.01% and for E≈105 V/cm, the error is 1%. In (15), the error is 7% for E≈102 V/cm and reduces to 2% for E between 103 V/cm and 104 V/cm.
Note that Romano [41] has dealt with another method based on a finite difference scheme of Scharfetter-Gummel type in order to simulate a consistent energy-transport model for electron transport in semiconductors devices. This last method is free of any fitting parameters. Otherwise, it is known that the maximum entropy moment systems of the gas-dynamical Boltzmann equation suffer from severe disadvantages which are related to the nonsolvability of an underlying maximum entropy moment problem unless restrictions on the choice of the macroscopic variables are made. Thus, Junk and Romano [42] showed that no such difficulties appeared in the semiconductor case if Kane’s dispersion relation is used for the energy band of electrons. In addition, a hydrodynamic subband model for semiconductors had been formulated by closing the moment system derived from the Schrödinger-Poisson-Boltzmann equations on the basis of the MEP where explicit closure relations for the fluxes and the production terms were obtained taking into account scattering of electrons with acoustic and nonpolar optical phonons, as well as surface scattering. For this model, a suitable numerical scheme was presented together with simulations of a nanoscale silicon diode [43]. Evenly, Camiola et al. [44] proposed an expression of the entropy combines quantum effects and semiclassical transport by weighting the contribution of each subband with the square modulus of the envelope functions arising from the Schrödinger-Poisson subsystem. The simulations shown that their model was able to capture the relevant confining and transport features, and, assessed the robustness of the numerical scheme. By taking into account nonparabolic energy bands of Kane’s type, Mascali and Romano [45] obtained explicit closure relations for fluxes and production terms, including scattering of electrons with acoustic and nonpolar optical phonons and surface scattering. Numerical simulations of a quantum diode showed the feasibility of their model and the importance of the nonparabolicity was assessed. Likewise, a recent unified physical approach developed by Trovato and Reggiani [23–26] deal with the same CR problem in quantum scheme, we will report in Section 5 some results of this latter for a comparison asserting that the knowledge of the L-Multipliers implies the determination of all CR for the system.
4. Information Theory and Closure Relations
In the microscopic description, several methods (kinetic theory and information theory) coincide in the first-order in the fluxes. Here, we want to check disagreements in higher-orders (second-order for instance). We think that production terms and closure relations are good fields to do that. In this work, we consider only the latter ones. We apply IT [17] to deal with closure relations of a gas in nonequilibrium steady state subject to heat flux and viscous pressure. The treatment proposed here differs from that we have presented in the precedent section. Indeed, there are two modifications: (i) from a microscopic point of view, it is based on a second-order expansion of a generalized distribution function obtained from IT and (ii) from a macroscopic point of view, it is based on extended thermodynamics but with a slightly different interpretation of the inequalities arising from the second law of thermodynamics.
Thus, we consider a gas in a nonequilibrium steady state characterized by numerical density n, internal energy, or local equilibrium temperature T and by a heat flux Q and viscous pressure θ^ij. The conditions to be imposed on the nonequilibrium distribution function fneq are
(16)∫fneqdk=∫feqdk=n,(17)∫fneqcdk=∫feqcdk=0,(18)∫fneq(12m*c2)dk=∫feq(12m*c2)dk=32nkBT,(19)∫fneq(12m*c2)cidk=Qi,(20)∫fneqc〈icj〉dk=θ^〈ij〉,
where feq is the local equilibrium distribution function corresponding to the number density n and the internal energy (3/2)nkBT. Since fneq differs from feq, the higher-order moments of fneq, as for instance corresponding to nonconserved quantities, will in general differ from the corresponding moments (19) and (20) of feq [46].
4.1. System under Heat Flux and Flux Limiters
When we take into account the heat flux as the only dissipative source, the maximum entropy solution for the distribution function has the form
(21)fneq=Aexp(-β12m*c2-(12m*c2-52β)γ·c),
where β and γ are the L-Multipliers to be determined from the constraints expressed by (18) and (19). A is a normalization constant to be obtained from (16). The Lagrange multiplier corresponding to the velocity has been written as -(5/2β)γ, a relation which, as it is known, satisfies the restriction (17). Therefore, the term (5/2β)γ·c follows from the requirement that the system is at rest; that is, 〈c〉=0 [47]; it behaves near equilibrium as (5/2β)c≈((5kBT)/2)c.
By using the shorthand notation α=(1/2)m*β and B(c)=((1/2)m*c2-(5/2β)), we expand (21) up to the second-order in γ as
(22)fneq=Aexp(-αc2)[1-B(c)γ·c+12B2(c)(γ·c)2].
Thus, if one restricts the expansion to the first-order in γ, one would obtain the Grad distribution function [10] corresponding to a steady heat flux without viscous pressure. Also this second-order expansion in fneq (which corresponds to a fourth-order expansion in c) is not exactly equal to the fourth-order expansion in the Grad’s scheme as it will be commented below. Note finally that according to the usual definition of the entropy in terms of the velocity distribution function f(23)ρs=-kB∫fln(f)dc.
We obtain, according to the restrictions imposed on f(24)ds=kBβdu+kBγ·dQ.
Therefore, the entropy is no longer the local equilibrium entropy [1, 2], but it depends also on the heat flux Q (through the Lagrange multiplier γ). Later, we will include also the viscous pressure θ^ij as a further independent quantity.
Introduction of (22) into (10) yields up to second-order in L-Multipliers corresponding to heat flux
(25)θ^ijk=15(Qiδjk+Qjδik+Qkδij),(26)θ^ijrr=5n2m*β2δij+C2ij(β,n)γiγj,
where C2ij(β,n)=A∫exp(-β((m*c2)/2))(B2(c)/2)c2ci2cj2dc.
We conclude that relation (25) coincides with (11) and by eliminating the second-order correction term in the right-hand side of (26) we find the same result as (12). Indeed, we can justify the second-order correction appearing in (26) by using a generalized distribution function.
Introduction of (22) into (16)–(19) yields the following simplified expressions:
(27)A=nm*(απ)3/2[1+58γ2αβ2]-1,(28)32nkBT=32nβ[1+158γ2αβ2][1+58γ2αβ2]-1,(29)Qi=-52nβ3m*γi[1+58γ2αβ2]-1.
Then the relation between Qi and γi can be written as
(30)Qi=-A′γi1+B′γi2,
where A′ and B′ are functions of m* and β and n are obtained directly from (28) and are written as A′=(5/2)(n/β3m*) and B′=(5/8)(1/αβ2).
Thus, the expression (30) yields a flux limiter for the heat flux because of the relation
(31)Q≤A′2B′=(54)1/2nm*(β)-(3/2)≈(54)1/2nm*(kBT)3/2.
This saturation behavior of heat flux is attained when γ=1/B′. So as to have a more accurate expression of heat flux limiter, we should take higher-order terms in (22). Then, we conclude that within this approach, the heat flow is limited (31) which could described a flux-limited phenomena as mentioned below (13).
4.2. System Subject to Reduced Heat Flux and Viscous Pressure
The generalized distribution function obtained from IT which describes a priori the nonequilibrium steady states of a system subjected to heat flux Q and viscous pressure θ^ij takes the following form:
(32)fneq=Aexp(-β12m*c2-(12m*c2-52β)ciγiiiiiiiiiiiiiiiii12-Γ〈jk〉:m*c〈jck〉).
The expansion of this latter up to the second-order in the L-Multipliers, β, γi, and Γ〈jk〉 conjugated, respectively, to energy, heat flux, and viscous pressure, is written with a trace-less symmetric tensor Λ as
(33)fneq=Aexp(-αc2)[1-B(c)(γici)-m*(Γ〈jk〉:Λ)cjckm*22iiiiiiiiiiiiiiiiiiiiiiiiii+γi22B2(c)ci2+m*22(Γ〈jk〉:Λ)2cj2ck2iiiiiiiiiiiiiiiiiiiiiiiiiim*22+m*γi(Γ〈jk〉:Λ)B(c)cicjck],
where A is the normalization constant.
Denoting that the condition imposed on the distribution function to give the reduced heat flux is written as
(34)∫fneqB(c)cidk=Qi.
Hence, in order to incorporate nonlinear effects in the CR we propose up to the second-order in reduced heat flux Qi and viscous pressure tensor θ^〈ij〉, the following expressions for the third- and fourth-order moments:
(35)θ^ijk=15(Qiδjk+Qjδik+Qkδij)+C^1(Qiθ^jk+Qjθ^ki+Qkθ^ij)+C^2Ql(θ^〈li〉δjk+θ^〈lj〉δki+θ^〈lk〉δij),(36)θ^ijrr=(5n2m*β2+C^3Q2)δij+C^4QiQj+C^5θ^〈ij〉+C^6θ^〈il〉θ^〈jl〉,
where C^i(i=1,…,6) are coefficients to be determined.
Note that all these coefficients are equal to zero in the simplified hypotheses which recover (11) and (12). In the following subsection, we will derive the coefficients C^i(i=1,…,6) within IT [17] by imposing some restrictions on Qi, θ^〈ij〉, θ^ijk, and θ^ijrr.
4.2.1. Calculations of C^i
The first two hypotheses imposed over Qi and θ^jk to evaluate C^1 (resp. C^2) are Q1≠0 and θ^23≠0 (resp. Q1≠0 and θ^11≠0), and, the other components of Qi and θ^ij are equal to zero. After carrying out the calculations, (35) yields the expressions of C^1 and C^2 up to the second-order in L-Multipliers
(37)C^1=θ^123Q1θ^23=25m*βn[1+58γ12αβ2+12(Γ〈23〉:Λ)2β2],(38)3C^1+2C^2=5θ^111-3Q15Q1θ^11=-6925m*βn[(Γ〈11〉:Λ)β]Φ,
with
(39)Φ=(1-(Γ〈11〉:Λ)β+58γ12α2β2+32(Γ〈11〉:Λ)2β2)×((1-275(Γ〈11〉:Λ)β)152(Γ〈11〉:Λ)2β2iiiiiiii×(1-3(Γ〈11〉:Λ)β+278γ12α2β2152(Γ〈11〉:Λ)2β2iiiiiiiiiiiiiii+152(Γ〈11〉:Λ)2β2))-1.
Note that we let (Γ〈11〉:Λ)/β term without simplification by β because it will play an important role as we will see later in this paper.
The second two hypotheses imposed over Qi and θ^ij to evaluate C^3 (resp. C^4) are Q1≠0 (resp. Q1≠0 and Q2≠0), and, the other components of Qi and θ^ij are equal to zero. After integration, (36) gives the following expressions for C^3and C^4 up to second-order in the L-Multipliers
(40)C^3+C^4=θ^11rr-5n/2m*β2Q12=16625βn[1+58γ12αβ2],(41)C^4=θ^12rrQ1Q2=11925βn[1+58γ12+γ22αβ2].
The third two hypotheses imposed over Qi and θ^ij to evaluate C^5 (resp. C^6) are θ^12≠0 (resp. θ^11≠0), and, the other components of Qi and θ^ij are equal to zero. Thus (36) yields the following expressions for C^5and C^6 up to second-order in the L-Multipliers:
(42)C^5=θ^12rrθ^〈12〉=721β,(43)C^6=θ^11rr-5n/2m*β2-C^5θ^〈11〉θ^〈11〉2=-214m*n([1+87(Γ〈11〉:Λ)β-757(Γ〈11〉:Λ)2β2]iiiiiiiiiiiiiiiii×[1-(Γ〈11〉:Λ)β+32(Γ〈11〉:Λ)2β2])×([1-3(Γ〈11〉:Λ)β+152(Γ〈11〉:Λ)2β2]2)-1.
Here, we should underline that the coefficients C^i(i=1,…,6) depend, in our case, not only on m*, n, and T as in [22], but on the L-Multipliers β, γi, and Γ〈jk〉 by taking into account only the second-order in the expansion of the distribution function. Note that explicit expressions of the generalized distribution function corresponding to each hypothesis are given in Appendix A.
4.2.2. Comparison with Anile and Muscato [22]
In zero-order in L-Multipliers [21] or in first-order in the L-Multipliers with (Γ〈11〉:Λ)/β is equal to zero, the coefficients C^i(i=1,…,6), that is, (37)–(43), become C^1≈(2/5)(m*β/n), C^2≈-(1/5)(m*β/n), C^3≈(47/25)(β/n), C^4≈(119/25)(β/n), C^5≈(7/2)(1/β), and C^6≈-(21/4)(m*/n). At this stage, the expressions for CR near equilibrium state take the forms
(44)θ^〈ijk〉≈25m*βn(Qiθ^jk+Qjθ^ki+Qkθ^ij12iiiiiiiiiiiiiiiii-12Ql(θ^〈li〉δjk+θ^〈lj〉δki+θ^〈lk〉δij)),(45)θ^ijrr≈(5n2m*β2+4725βnQ2)δij+11925βnQiQj+721βθ^〈ij〉-214m*nθ^〈il〉θ^〈jl〉.
When we neglect the contribution of γi and Γ〈jk〉, we can recover easily (14) and (15) established by Anile and Muscato [22], with a slight suitable difference concerning the definitions of Q, θ^ij, θ^ijk, and θ^ijrr. In our opinion, the main difference is related to the use of different microscopic definitions of the heat flow and a vanishing expression of (Γ〈11〉:Λ)/β. We will study later in this paper a nonvanishing expression of (Γ〈11〉:Λ)/β which will give us a solution to overcome the slight differences among our expressions and those established by Anile and Muscato [22].
Otherwise to evaluate CR, Anile and Muscato [22] used for the heat-flux the expression Q=(1/2)m*c2c but in our case we have used Q=((1/2)m*c2-(5/2β))c. Thus, the two definitions gave the same average value for Q, as we have described above, but they yield different values for higher-order moments. The adopted definition in our calculation corresponds to what one calls the reduced heat flow, which it is used in fluctuation-dissipation theory where moments of Q play an important role.
Furthermore, if we simply compare (14) and (15) with those established in our approach (44) and (45), we can estimate the unknown free parameter σ. In fact, a definite value of this latter is provided by the maximum entropy approach or alternatively it may be determined by Monte Carlo simulations [1]. This is what we deal with in our approach based on IT. It is clear from (15) and (45) that the unknown parameter σ is connected to C^3 (resp. C^4). It is a nondimensional free parameter and a function of the lattice temperature TL (300 K) and the electron temperature T as σ≈(57/25)(MTL/nT) (resp. σ≈(119/50)(TL/T)) with β≈(1/kBT) in the case of partial thermal equilibrium.
So as to compare our results with those obtained in Monte Carlo simulations, we plot in Figure 1 for the same values of TL and T versus field E. The data of these latters (see Table 1) are extracted from [22]. Thus, we choose the second expression of σ related to C^4 because the authors were silent about the expression of the parameter M used in expression (15)
(46)σ≈11950TLT.
Estimates of the free parameter σ, (46), by using data extracted from [22].
Electron temperature: T (K)
275.27
281.77
348.27
449.46
1694.33
3281.89
5123.61
Electric field: E (V/cm)
1109
1449
5518
10768
41000
61230
80000
Free parameter: σ
2.594
2.534
2.050
1.589
0.421
0.218
0.139
The free parameter σ from (46) versus electric field E.
As shown in Figure 1, we find that free parameter σ which is obtained by using a fitting of Monte Carlo simulation data [22] becomes approximately constant at high field more than 60000 V/cm and σ≤0.139.
In Figure 2, we plot the impact of electron temperature T on the nondimensional free parameter σ. It is clear from Figure 2 that at high electron temperature T is more than 3500 K and σ≤0.13. This last σ-value is obtained without using any fitting. Hence, we can use the expression (46) to evaluate the value of the electric field which could be applied in future experiences dealing with electron transport in semiconductors and to compare with Monte Carlo simulation results, vice versa to have electron temperature T if we know electric field E.
Determination of the free parameter σ from (46) by varying electron temperature T, without using a fit nor data of Monte Carlo simulation.
4.2.3. Overcoming Slight Differences in CR with Anile and Muscato [22]
In order to overcome the slight difference between our results and those performed by Anile and Muscato [22], we have to take into account, in the first-order in L-Multipliers, a nonvanishing expression of (Γ〈11〉:Λ)/β. Denoting firstly by C^iA(i=1,…,6), the coefficients related to (14) and (15) which are written as C^1A≈(2/15)(m*β/n), C^2A≈-(4/75)(m*β/n), C^3A≈((σ/MkBTL)-(2/5)(β/n)), C^4A=(2σ/nkBTL), C^5A≈(7/2)(1/β), and C^6A=(m*/n) with (1/kBT)≈β. Thus, by subtracting from (38) three times (37) and adding (38) to (37), that is, (38) −3 × (37) and (38) + (37), we obtain the coefficients C^2 and C^1 in the first-order in L-Multipliers with a nonvanishing (Γ〈11〉:Λ)/β(47)C^2=C^2A[2078(Γ〈11〉:Λ)β+454],(48)C^1=-15C^1A[2078(Γ〈11〉:Λ)β+454]-m*βn[69100(Γ〈11〉:Λ)β+220].
We find from (47) that for (Γ〈11〉:Λ)/β=-(82/207) we have C^2=C^2A. This last equality into (48) yields C^1=C^1A. Thus, introducing the latter nonvanishing expression of (Γ〈11〉:Λ)/β into (43) without any trancation in L-Multipliers, it gives C^6=0.857C6A. Then, the closure relations (35) and (36) become
(49)θ^〈ijk〉≈25m*βn(Qiθ^jk+Qjθ^ki+Qkθ^ij215iiiiiiiiiiiiiiiii-215Ql(θ^〈li〉δjk+θ^〈lj〉δki+θ^〈lk〉δij)),θ^ijrr≈(5n2m*β2+4725βnQ2)δij+11925βnQiQj+721βθ^〈ij〉+0.857m*nθ^〈il〉θ^〈jl〉.
By comparing (49) with (14) and (15), we conclude that we have overcome the silght difference arising at zero-order and at first-order in L-Multipliers where we have used a vanishing expression of (Γ〈11〉:Λ)/β. Thus, the accuracy of our approach checked against the results of Anile and Muscato [22], involving the percent error between C^i and C^iA is asserted such as C^1/C^1A=100%, C^2/C^2A=100%, C^5/C^5A=100%, and C^6/C^6A=85.7%. We recall, as we have processed above that C^4 and C^4A (resp. C^3 and C^3A) allow us to find the nondimensional free parameter σ. It is clear now that taking into account higher-order in L-Multipliers will give better accurate results as it is seen from the calculation of coefficient C^6 determined without any truncation in L-Multipliers.
4.2.4. Closure Relations and Grad’s Method [10]
The most well-known assumed form for the distribution function is the Grad type which is based on an expansion of the distribution function in Hermite polynomials [10]. This decouples the dependence of the closure coefficients and greatly simplifies the derivation. In fact, Grad considered both 13- and 20-moment closures, extensions to many moments have been considered by others [13, 48]. However, it is an unfortunate fact that members of the Grad hierarchy suffer from several problems: (i) the distribution function is not always positive; it is therefore not a properly defined probability density function and (ii) for modest departures from local equilibrium, it is possible for the resulting moment equations to become nonhyperbolic [49–51]. Here, we assert that Grad’s method fall in nonequilibrium situation. Indeed, we expand the distribution function fneq up to the fourth-order in the power of the velocity (this is the same order used in the precedent section). In the notation of Hermite polynomials where Hi(n) is a tensor of order n, Grad’s distribution function writes as
(50)fG=feq∑n=041n!ai(n)Hi(n),
where the coefficient ai(n) is also a tensor of order n.
The first four closure relations are written as
(51)θ^ijk=∑n=14θ^ijkn,θ^ijrr=∑n=14θ^ijrrn.
In order to derive the coefficients C^i appearing in (35) and (36), we should use the same conditions as in Section 4.2.1; that is, Q1≠0 and θ^23≠0, and, the other components are equal to zero. Thus, up to the third-order Grad’s distribution function takes the form
(52)fG=feq[1+12p23pc2c3+16Q1c1pRT(c2-5)].
According to the last equation we observe that θ^123=0, which means that C^1=0. Grad’s distribution function can be generalized by including the full third-order contributions
(53)fG=feq[1+θ23pc2c3+Q1c1pRT(c2-5)iiiiiiiiiii+θ123pRTc1c2c3+⋯],
where θ^123 is arbitrary.
Hence, if one imposes only conditions on p23 and on Q1 but not on the additional third-order terms, we observe that C^1 is arbitrary too, since C^1 is given by C^1=θ^123/(Q1θ^23). We checked also the same method for the other coefficients C^2,C^3,…,C^6 in (35) and (36). We concluded that at higher-order, if conditions are imposed only on heat flux and viscous pressure, Grad’s method cannot give any restrictions in the closure relations (35) and (36). We can understand the difference between this model and our approach based on IT. Consequently, if we compare fneq with fG used to provid θ^123, we will see clearly that in our approach θ^123 is connected in a unique manner to γ1 and Γ〈23〉 but in Grad’s method it is an independent coefficient. In fact, the distribution function of Grad assumes to θ^123=0 in the case, where only Q1≠0 and θ^23≠0. Thus, this freedom of coefficients in Grad’s method is present also for the others coefficients C^2,C^3,…,C^6.
5. Nonlinear Closure Relations in Quantum Scheme
The present section is devoted to compare the aims and the results developed in Section 4 with those of the work of Trovato and Reggiani [23–26]. The authors proposed MEP in its local semiclassical and nonlocal quantum formulation as the basis of an unified physical approach. In fact, the local semiclassical case investigated MEP both in the usual kinetic theory of degenerate gases and in transport phenomena for hot carriers in nondegenerate semiconductors and mostly overviewed results obtained in the last decade. The nonlocal quantum case mostly summarized recent results and presented a series of original findings. Otherwise, Romano [52] proposed quantum corrections to the hydrodynamical model of semiconductors based on MEP were obtained at ℏ2 order with a Chapman-Enskog expansion in the high field approximation, modeling the ℏ2 part of the collision term in a relaxation form. The author deduced limiting energy-transport and drift-diffusion models.
As constraints in QMEP approach, Trovato and Reggiani [24] have determined a closed quantum hydrodynamic (QHD) system for macroscopic variables. By showing that only a higher-order expansion of the distribution function can be fruitfully applied to describe transport phenomena under conditions extremely far from thermodynamic equilibrium in the presence of very high electric fields and associated gradients such as those occurring in submicron devices. Therefore, to describe accurately the transport phenomena in submicron semiconductor devices, they have considered a nonlinear expansion of the distribution function in terms of the first thirteen relevant macrovariables involved physical interpretation. For validation purposes, the approach was applied to many one-dimensional n+nn+ submicron Si structures by using different band structure models, different doping profiles, different lattice temperatures, and different applied biases. From the numerical results, it is shown that only a strongly nonlinear expansion provided a more accurate description of all moments [25]. We report below an abstract of those works which concern only nonlinear extended thermodynamics with thirteen moments and we comment briefly the main characteristics of this recent unified physical approach.
5.1. Unified Physical Approach due to Trovato and Reggiani [25]
In the local semiclassical case, Trovato and Reggiani [25] introduced BTE for a fermion or bose gas by developing a local theory with a temperature scheme to describe the state of a gas in which the moments of the distribution function f(r,u,t) are taken as basic fields. Through a formal expansion around the local equilibrium configuration, they explained the extended thermodynamic theories of N moments and degree α, by determining a set of closed hydrodynamic systems for bosons and/or fermions, in the frame work of the usual gas dynamics. By introducing the kinetic field
(54)Ψ{n,mui,mu2,mu〈iuj〉,mu2ui},
where u denotes velocity vector; they obtained the corresponding first thirteen macroscopic variables, that is, the masse density n, the moment density mui, the energy density mu2, the traceless momentum flux density mu〈iuj〉, and the energy flux density mu2ui that admit a direct physical interpretation. In the same way, by defining the higher-order moments, the external field productions and collisional productions, they wrote the balance equations for usual first thirteen moments.
Thus, by considering the decomposition of a moment in its convective and central parts, they formally obtained the new set of variables {n,P,M〈ij〉,Qi}, n being the numerical density, P the pressure, M〈ij〉 the stress deviator, and Qi the heat flux. Analogously, they obtained the following constitutive equations for the third- and fourth-order moments tensors θ^ijr and θ^illk as
(55)θ^〈ijk〉=Mijr+3M(〈ij〉vk)-25vk(M〈ri〉δjk+M〈rj〉δki+M〈rk〉δij)+mnv(ivjvk),θ^illk=Mijr+2vrM〈rik〉+7Pvivk+285Q(ivk)+45(Qrvr)δik+v2M〈ik〉+Pv2δik+2vr(M〈ri〉vk+M〈rk〉vi)+mnv2vivk,
where vi is the mean velocity.
The central moments Mi1⋯ink are expressed as
(56)Mi1⋯isk=∫mu~i⋯u~isu~kf(r,u,t)du,
where u~=ui-vi the peculiar velocity and s=0,1,…,M by construction with arbitrary values for the integer M.
The higher-order moments θ^〈ijk〉 and θ^illk (55) are determined by decomposing the distribution function and L-Multipliers into equilibrium and nonequilibrium parts where L-Multipliers are formulated versus bose and fermi integrals. In fact, by determining an analytic expression for the L-Multipliers versus n, P, Qi, M〈ij〉 and the fermi and bose integral functions, Trovato and Reggiani [25] obtained an explicit distribution function up to the second-order of the local equilibrium. Thus, all the results obtained for the constitutive relations are determined in terms of the fermi and/or bose functions. Therefore, only by means of an efficient evaluation of these integrals, it is possible to analyse the gas behavior for various degeneracy levels.
5.2. Comparison with Trovato and Reggiani [25]
Within the QMEP Trovato and Reggiani [25] introduced a quantum entropy functional of the reduced density matrix; the principle of quantum maximum entropy is then asserted as fundamental principle of quantum statistical mechanics. Accordingly, they developed a comprehensive theoretical formalism to construct rigorously a closed QHD transport within a Wigner function approach [25]. Their approach is formulated both in thermodynamic equilibrium and nonequilibrium conditions, and the quantum contributions are obtained by only assuming that the L-Multipliers can be expanded in powers of ℏ2, while the classical results are recovered in the limit ℏ→0. The authors asserted that the knowledge of the L-Multipliers involves the determination of all CR for the system, up to the first quantum correction [24]. Thus, to prove the validity of the QMEP results in the case of nonlinear extended thermodynamics with thirteen moments, Trovato and Reggiani compared the analytic results existing in the literature with some recent numerical methods for the evaluation of bose and fermi integrals.
Hence it is clear that the L-Multipliers play an important role in the conception dealing with CR in hydrodynamical models both in the relevant contributions of Trovato and Reggiani and in the present paper. In our case, we have assumed some restrictions on L-Multipliers where third- and fourth-order CR for the moments of the distribution function are expressed merely in terms of second-order L-Multipliers: (i) to recover the formulation of flux-limiters (30) and (ii) to overcome difficulties in comparing our results with those of Anile and Muscato [22]. In fact, we have used on the one hand a nonvanishing expression of (Γ〈11〉:Λ)/β at the first-order in L-Multipliers to find the accurate coefficients C^1 and C^2, and on the other hand the fall expression of C^6 without any truncation in L-Multipliers (49). Also, the fact that third- and fourth-order moments are related to the numerical values of bose and fermi integrals, in the unified physical approach due to Trovato and Reggiani, does not allow us to reach one of the main objectives of the present paper which is to find an analytical expression of the nondimensional free parameter σ (46).
6. Concluding Remarks
In this paper, we have studied within extended thermodynamics implications of higher-order moments, in particular second-order moments, so as to solve the closure problem presented in hydrodynamical models. This is done by means of a generalized distribution function which is obtained from maximum-entropy argument. Thus, we have compared the results of our approach with those obtained by both models based on the assumption of a drift-Maxwellian distribution function, and, on extended thermodynamics and entropy principle; then we have found a good agreement. Further, we have applied Grad’s method to solve the same nonlinear CR. In the first part, we have restricted ourself using only heat flux as dissipative source. The CR found are the same as those obtained by a drift-Maxwellian assumption of the distribution function. We have been led to an expression for the heat flux describing the flux limited phenomena occurring in the submicron-electron devices. In the second part, we have considered a generalized distribution function for a nonequilibrium system under heat flux and viscous pressure in order to obtain constitutive equations of the third-order and the fourth-order moments. Indeed, we have found with a slight difference similar constitutive equations as those established by Anile and Muscato [22]. This slight difference which was published previously arised from the use of a vanishing first-order in L-Multipliers. Here, it was overcome by using a nonvanishing one.
Furthermore, we conclude that considerations proposed in this paper, which concern the definition of the reduced heat flux, allow us to express the nondimensional free parameter σ appearing in (15) as function of the lattice temperature TL and the electron temperature T. In fact, the σ-parameter appearing in fourth-order moment tensor, for partial thermal equilibrium, is related both to the Monte Carlo simulation and to a fitting as described by Anile and coworkers. Finally, we have used Grad’s method to calculate the CR for higher-order moments in a simple case. We have found that this method could not give any restrictions about higher-order moments. Also, we assert as Trovato and Reggiani [23–26] that the knowledge of the L-Multipliers involves the determination of all CR for the system.
In brief, we have constructed a generalized distribution function in the frame work of IT. This latter allowed us to reproduce the constitutive equations of the third-order and the fourth-order moments in terms of second-order L-Multipliers and to take into account the nonlinear effects in transport-hydrodynamical models of submicron-electron devices. We should not forgot that the generalized distribution function proposed in this paper is one of several contributions dealing with the extension of the canonical distribution function by adding to the equilibrium Hamiltonian an effective Hamiltonian due to the flow [1, 2, 53]. Thus, we have considered in (21) heat flux as a supplementary constraint, and in (32) both heat flux and viscous pressure as supplementary constraints in order to conceive the generalized distribution function used in our approach. At the end, for a full exploitation of the results mentioned in the present paper involving the nonlinear CR and the maximum-entropy should be tested against the Monte Carlo simulation incorporating second-order L-Multipliers.
AppendicesA. Explicit Expressions of fneq Used in Section 4.1
The present appendix contains the explicit expressions of the generalized distribution function corresponding to each hypothesis proposed in our paper in order to derive the coefficients Ci, that is, (37), (38), (40), (41), (42), and (43).
For the first tow hypotheses, the distribution function is written, respectively,
(A.1)fneq=Aexp(-αc2)[1-B(c)γ1c1-m*(Γ〈23〉:Λ)c2c3γ122iiiiiiiiiiiiiiiiiiiiiiiiiii+γ122B2(c)c12+m*22(Γ〈23〉:Λ)2c22c32iiiiiiiiiiiiiiiiiiiiiiiiiiiγ122+m*γ1(Γ〈23〉:Λ)B(c)c1c2c3],fneq=Aexp(-αc2)[1-B(c)γ1c1-m*(Γ〈11〉:Λ)c12γ122iiiiiiiiiiiiiiiiiiiiiiiiiii+γ122B2(c)c12+m*22(Γ〈11〉:Λ)2c14iiiiiiiiiiiiiiiiiiiiiiiiiiiγ122+m*γ1(Γ〈11〉:Λ)B(c)c13].
In the conditions of the second two hypotheses, the distribution function is written, respectively,
(A.2)fneq=Aexp(-αc2)[1-B(c)γ1c1-m*(Γ〈23〉:Λ)c2c3γ122iiiiiiiiiiiiiiiiiiiiiiiiiii+γ122B2(c)c12+m*22(Γ〈23〉:Λ)2c22c32iiiiiiiiiiiiiiiiiiiiiiiiiiiγ122+m*γ1(Γ〈23〉:Λ)B(c)c1c2c3],fneq=Aexp(-αc2)[1-B(c)γ1c1-m*(Γ〈11〉:Λ)c12γ122iiiiiiiiiiiiiiiiiiiiiiiiiii+γ122B2(c)c12+m*22(Γ〈11〉:Λ)2c14iiiiiiiiiiiiiiiiiiiiiiiiiiiγ122+m*γ1(Γ〈11〉:Λ)B(c)c13].
The distribution functions related to the third two hypotheses are written, respectively,
(A.3)fneq=Aexp(-αc2)[1-B(c)γ1c1+γ122B2(c)c12],fneq=Aexp(-αc2)[1-B(c)γ1c1-B(c)γ2c2γ122iiiiiiiiiiiiiiiiiiiiiiiiiii+γ122B2(c)c12+γ222B2(c)c22].
B. Some Equations Used in Grad’s Method
We underline here some equations used in Grad’s method to deal with CR. Indeed, the first four Hermite polynomials Hi(n) appeared in the distribution function of Grad (50) [10] may be written as
(B.1)H(0)=1,Hi(1)=ci,Hij(2)=cicj-δij,Hijk(3)=cicjck-(ciδjk+cjδik+ckδij),Hijkl(4)=cicjckcl-(cicjδkl+cickδjl+ciclδjk+cjckδil+ciclδikiiiiii+ckclδij)+(δijδkl+δikδjl+δilδjk),
and the coefficient tensors ai(n) are expressed in terms of the elementary moments as
(B.2)a(0)=1,ai(1)=0,aij(2)=pijp,aijk(3)=sijkpRT,aijkl(4)=MijklpRT-1p(pijδkl+pikδjl+pilδjk+pjkδiliiiiiiiiiiiiiiiiii+pjlδik+pklδij)-(δijδkl+δikδjl+δilδjk),
where pmn, smnl, and Mijkl take the form
(B.3)pij=∫cicjfdξ,sijk=∫cicjckfdξ,Mijkl=∫cicjckclfdξ.
Then, the simplified expressions of the closure relations (51) are given by
(B.4)θ^ijk(0)=Am*22∫exp(-αc2)cicjckdc=0,θ^ijk(1)=Am*22∫cicjckexp(-αc2)(0cm)dc=0,θ^ijk(2)=Am*2pmn4p:∫exp(-αc2)(cmcn-δmn)cicjckckdc=0,θ^ijk(3)=Am*2smnl12pRT:∫exp(-αc2)×(cmcncl-(cmδnl+cnδml+clδmn))cicjckdc,θ^ijk(4)=0,θ^ijrr(0)=m*22∫feqcicjc2dc,θ^ijrr(1)=0,θ^ijrr(2)=m*2pmn4p:∫feq(cmcn-δmn)cicjc2dc,θ^ijrr(3)=0,θ^ijrr(4)=m*2214!amnlr(4):∫feqHmnlr(4)cicjc2dc.
C. Useful Integrals
This appendix contains some useful integrals that are needed to calculate equations presented in this paper. Let F(c) be any scalar function of c such that the integrals appearing below converge, and let cx and cy be the components of the relative velocity c [1]. Then
(C.1)∫F(c)cx2dc=13∫F(c)c2dc,∫F(c)cx4dc=15∫F(c)c4dc,∫F(c)cx2cy2dc=115∫F(c)c4dc.
The following definite integrals are also useful:
(C.2)∫0+∞exp(-αc2)crdc=π2123252⋯r-12α-(r+1)/2,iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiireven,∫0+∞exp(-αc2)crdc=12(r-12)!α-(r+1)/2,rodd.
Conflict of Interests
The author declares that there is no conflict of interests in this paper.
Acknowledgments
The author thanks Professor Jou, D. from Grup de Fsica Estatistica Departament de Fsica Universitat Autonóma de Barcelona (Spain), for his precious remarks, careful reading of the paper, and helpful discussions. The author thanks also Professor Trovato, M. from Dipartimento di Matematica Università di Catania (Italy) for giving his relevant contributions [23–26] and to Professor Bennai, M. for hosting this work at LPMC. A. Salhoumi’s current address is Faculté des Sciences Juridiques Economiques et Sociales Ain Sebaa, Universitxé Hassan II-Mohammedia, Ain sebaa Morocco while his hosting address is Laboratoire de Physique de la Matière Condensée (LPMC), Faculté des Sciences Ben M’sik, Université Hassan II-Mohammedia, Morocco.
JouD.Casas-VazquezJ.LebonG.2010Berlin, GermanySpringerLebonG.JouD.Casas-VazquezJ.2008Berlin, GermanySpringerHänschW.1991Berlin, GermanySpringerBlotekjaerB. K.Transport equations for electrons in two-valley semiconductors197017138472-s2.0-0014705867BaccaraniG.WordemanM. R.An investigation on steady-state velocity overshoot in silicon198229970977HänschW.Miura-MattauschM.The hot-electron problem in small semiconductor devices19866026506562-s2.0-000026362210.1063/1.337408WoolardD. L.TianH.TrewR. J.LittlejohnM. A.KimK. W.Hydrodynamic electron-transport model: nonparabolic corrections to the streaming terms1991442011119111322-s2.0-000099731810.1103/PhysRevB.44.11119ThomaR.EmundsA.MeinerzhagenB.PeiferH.-J.EnglW. L.Hydrodynamic equations for semiconductors with nonparabolic band structure1991386134313532-s2.0-002617080210.1109/16.81625AnileA. M.PennisiS.Thermodynamic derivation of the hydrodynamical model for charge transport in semiconductors1992462013186131932-s2.0-000129710110.1103/PhysRevB.46.13186GradH.FlüggeS.Principles of the kinetic theory of gases in thermodynamics of gases195812Berlin, GermanySpringerBanachZ.LareckiW.Nine-moment phonon hydrodynamics based on the modified Grad-type approach: formulation20043741980598292-s2.0-754423417110.1088/0305-4470/37/41/013MullerI.RuggeriT.1993New York, NY, USASpringerMullerI.RuggeriT.1998New York, NY, USASpringerAnileA. M.PennisiS.Extended thermodynamics of the Blotekjaer hydrodynamical model for semiconductors1992431871972-s2.0-1544437418010.1007/BF01130290AnileA. M.PennisiS.MuscatoM.DonatoA.OliveriF.Extended thermodynamics of charge carrier transport in semiconductors199343Braunschweig, GermanyViewegNotes on Numerical Fluid Mechanicsde GrootS. R.MazurP.1985New York, NY, USADover PublicationGrayR. M.2010Berlin, GermanySpringerGrasserT.TangT.KosinaH.SelberherrS.A review of hydrodynamic and energy-transport models for semiconductor device simulation20039122512742-s2.0-3364689348510.1109/JPROC.2002.808148la RosaS.MascaliG.RomanoV.Exact maximum entropy closure of the hydrodynamical model for sisemiconductors: the 8-moment case20097037107342-s2.0-6834912496110.1137/080714282McDonaldJ. G.GrothC. P. T.Towards realizable hyperbolic moment closures for viscous heat-conducting gas ows based on a maximum-entropy distribution201325557360310.1007/s00161-012-0252-ySalhoumiA.ZakariM.BoughalebY.Nonlinear efiects in 3rd and 4th order moments of transporthydrodynamical models in semiconductors within extended thermodynamics199927175AnileA. M.MuscatoO.Extended thermodynamics tested beyond the linear regime: the case of electron transport in silicon semiconductors1996831311422-s2.0-0030548243TrovatoM.ReggianiL.Quantum hydrodynamic models from a maximum entropy principle201043102-s2.0-7714913272210.1088/1751-8113/43/10/102001102001TrovatoM.ReggianiL.Quantum maximum entropy principle for a system of identical particles20108122-s2.0-7724917569910.1103/PhysRevE.81.021119021119TrovatoM.ReggianiL.Maximum entropy principle and hydrodynamic models in statistical mechanics20123599266TrovatoM.ReggianiL.Quantum maximum entropy principle for fractional exclusion statistics20131102020404WignerE.On the quantum correction for thermodynamic equilibrium19324057497592-s2.0-3374501474210.1103/PhysRev.40.749KittelC.1963New York, NY, USAJohn Wiley & SonsAnileA. M.MuscatoO.Improved hydrodynamical model for carrier transport in semiconductors1995512316728167402-s2.0-000050577910.1103/PhysRevB.51.16728StettlerM. A.AlamM. A.LundstromM. S.A critical examination of the assumptions underlying macroscopic transport equations for silicon devices19934047337402-s2.0-002757457710.1109/16.202785AnileA. M.RomanoV.Covariant flux-limited diffusion theories199238613253292-s2.0-11944254746LevermoreC. D.Relating Eddington factors to flux limiters19843121491602-s2.0-48749142061BellA. R.EvansR. G.NicholasD. J.Elecron energy transport in steep temperature gradients in laser-produced plasmas19814642432462-s2.0-000110484110.1103/PhysRevLett.46.243EpperleinE. M.RickardG. J.BellA. R.Two-dimensional nonlocal electron transport in laser-produced plasmas19886121245324562-s2.0-394311272210.1103/PhysRevLett.61.2453ZakariM.A continued-fraction expansion for flux limiters19972403-46766842-s2.0-0031169721SalhoumiA.ZakariM.JouD.Casa-VazquezJ.BoughalebY.A brief glimpse of extended irreversible thermodynamics theory and some of its applications199816270SalhoumiA.ZakariM.BoughalebY.Another contribution to calculate ballistic velocity by means of continued fraction expansion within extended irreversible thermodynamics20002755-64864922-s2.0-003470648410.1016/S0375-9601(00)00626-5VazquezF.Size efiects on heat transport in small systems: dynamical phase transition from difiusive to ballistic regime20091056491564919AnileA. M.AllegrettoW.RinghoferC.20031823Berlin, GermanySpringerLecture Notes in MathematicsAnileA. M.MarroccoA.RomanoV.SellierJ. M.2D numerical simulation of the MEP energy-transport model with a mixed finite elements scheme200543-42312592-s2.0-3144444253710.1007/s10825-005-5039-yRomanoV.2D numerical simulation of the MEP energy-transport model with a finite difference scheme200722124394682-s2.0-3384639076410.1016/j.jcp.2006.06.028JunkM.RomanoV.Maximum entropy moment system of the semiconductor Boltzmann equation using Kane's dispersion relation20051732472672-s2.0-2324443745910.1007/s00161-004-0201-5MascaliG.RomanoV.Numerical simulation of a hydrodynamic subband model for semiconductors based on the maximum entropy principle201033155163CamiolaV. D.MascaliG.RomanoV.Numerical simulation of a double-gate MOSFET with a subband model for semiconductors based on the maximum entropy principle2012244174362-s2.0-8205517603910.1007/s00161-011-0217-6MascaliG.RomanoV.A non parabolic hydrodynamical subband model for semiconductors based on the maximum entropy principle2012553-4100310202-s2.0-8485519159310.1016/j.mcm.2011.09.026Casas-VázquezJ.JouD.Nonequilibrium temperature versus local-equilibrium temperature1994492104010482-s2.0-3375041060310.1103/PhysRevE.49.1040DomínguezR.JouD.Thermodynamic pressure in nonequilibrium gases1995511158163StruchtrupH.2005Berlin, GermanySpringerBriniF.Hyperbolicity region in extended thermodynamics with 14 moments2001131182-s2.0-003553041010.1007/s001610100036TorrilhonM.Characteristic waves and dissipation in the 13-moment-case20001252893012-s2.0-0034344280TorrilhonM.Hyperbolic moment equations in kinetic gas theory based on multi-variate Pearson-IV-distributions2010746396732-s2.0-7794975822410.4208/cicp.2009.09.049RomanoV.Quantum corrections to the semiclassical hydrodynamical model of semiconductors based on the maximum entropy principle200748122-s2.0-3764901133410.1063/1.2819600123504JouD.Casas-VázquezJ.Criado-SanchoM.Information theory and thermodynamics of polymer solutions under flow19992621-269752-s2.0-0032761524