Hydrodynamical modeling of charge transport in semiconductors

Enhanced functional integration in modern electron devices requires an accurate modeling of energy transport in semiconductors in order to describe high-field phenomena such as hot electron propagation, impact ionization and heat generation in the bulk material. The standard drift-diffusion models cannot cope with high-field phenomena because they do not comprise energy as a dynamical variable. Furthermore for many applications in optoelectronics one needs to describe the transient interaction of electromagnetic radiation with carriers in complex semiconductor materials and since the characteristic times are of order of the electron momentum or energy flux relaxation times, some higher moments of the distribution function must be necessarily involved. Therefore these phenomena cannot be described within the framework of the drift-diffusion equations (which are valid only in the quasi-stationary limit). Therefore generalizations of the drift-diffusion equations have been sought which would incorporate energy as a dynamical variable and also would not be restricted to quasi-stationary situations. These models are loosely speaking called hydrodynamical models. One of the earliest hydrodynamical models currently used in applications was originally put forward by Blotekjaer [1] and subsequently investigated by Baccarani and Wordeman [2] and by other authors [3]. Eventually other models have also been investigated, some including also non-parabolic effects [4-6, 8-20]. Most of the implemented hydrodynamical models suffer from serious theoretical drawbacks due to the ad hoc treatment of the closure problem (lacking a physically convincing motivation) and the modeling of the production terms (usually assumed to be of the relaxation type and this, as we shall see, leads to serious inconsistencies with the Onsager reciprocity relations). In these lectures we present a general overview of the theory underlying hydrodynamical models. In particular we investigate in depth both the closure problem and the modeling of the production terms and present a recently introduced approach based on the maximum entropy principle (physically set in the framework of extended thermodynamics [21, 22]). The considerations and the results reported in the paper are exclusively concerned with silicon.