Existence of Weak Solutions to a Class of Degenerate Semiconductor Equations Modeling Avalanche Generation
Bin Wu1,2
Academic Editor: Shijun Liao
Received22 Feb 2009
Accepted03 Jun 2009
Published30 Jul 2009
Abstract
We consider the drift-diffusion model with avalanche generation for evolution
in time of electron and hole densities , coupled with the electrostatic potential in a semiconductor device. We also assume that the diffusion term is degenerate. The existence of local weak solutions to this Dirichlet-Neumann mixed boundary value problem is obtained.
1. Introduction
In this paper, we consider the following degenerate semiconductor equations modeling avalanche generation:
with initial and boundary conditions
Here the unknowns and denote the electrostatic potential, the electron density, and the hole density, respectively. The boundary consists of two disjoint subsets and . The carrier densities and the potential are fixed at (Ohmic contacts), whereas models the union of insulating boundary segments. represents the electron current, and is the analogously defined physical quantity of the positively charged holes. Function denotes the doping profile (fixed charged background ions) characterizing the semiconductor under consideration, while the term models the effect of impact ionization (avalanche generation of charged particles) (cf. [1, 2] for details). is the net recombination-generation rate, where characterizes the mechanism of particle transition. The constant is the adiabatic or isothermal (if ) exponent. The regime describes a fast diffusion process in the electron (hole) density, whereas is related to slow diffusion.
The standard drift-diffusion model corresponding to has been mathematically and numerically investigated in many papers (see [3β6]). Existence and uniqueness of weak solutions have been shown. The standard model can be derived from Boltzmann's equation once assumed that the semiconductor device is in the low injection regime, that is, for small absolute values of the applied voltage. In [7] JΓΌngel showed that in the high-injection regime diffusion terms are no longer linear. A useful choice for is . In this case, the parabolic equations (1.2) and (1.3) become of degenerate type, and existence of solutions does not follow from standard theory. Recently, many authors [8β10] have studied the existence and uniqueness of weak solutions of this type of degenerate semiconductor equations without avalanche generation term. In [9], the degenerate semiconductor equations based on Fermi-Dirac statistics were introduced by JΓΌngel for the first time. The existence and uniqueness results are shown under the assumption that the solution of Poisson equation with Dirichlet-Neumann mixed boundary conditions had the regularity , this amounts to a geometric condition on , for example and . ([11, Theorem 3.29]). Then Guan and Wu [8] obtain similar results without the assumption above.
There are some papers concerning the semiconductor equations modeling avalanche generation. For instance, the existence of weak solutions of nondegenerate stationary problem has been investigated in [12, 13]. When , that is, the diffusion term is not degenerate, the authors [14] obtained the existence of local weak solutions of problem (1.1)β(1.6).
Our main goal in this paper is to study the existence of weak solutions of problem (1.1)β(1.6). In contrast to the above works, degeneration of diffusion term we are going to study introduces significant new technical difficulties to estimate the avalanche term.
We make the following assumptions:
(H1) is bounded and , whose outward normal vector is and , ;(H2);(H3) is a locally Lipschitz continuous function defined for and ;(H4), =ββconst., ;(H5), , and , in ;(H6) and a.e. in .
Let
Definition 1.1. is called the weak solution to the problem (1.1)β(1.6) if , , , , and there hold
Our main result in this paper is as follows.
Theorem 1.2. Under hypotheses (H1)β(H6), there exists at least one local weak solution to the problem (1.1)β(1.6).
2. Approximate Problem
For simplicity, we assume that . As [14], we first construct the following bound approximate sequence of avalanche generation term :
here . Obviously, .
Now we introduce the following approximate problem with the initial and boundary conditions (1.4)β(1.6):
This section is devoted to the proof of global existence of weak solutions to the above approximate problem (2.2)β(2.4), (1.4)β(1.6). We will prove the following existence theorem.
Theorem 2.1. Under hypotheses (H1)β(H6), there exists at least one global weak solution to the problem (2.2)β(2.4), (1.4)β(1.6).
The proof is based on Schauder's fixed pointed theorem. The main difficulty in the proof is that problem (2.2)β(2.4), (1.4)β(1.6) is degenerated at points where . This difficulty leads us to consider the following auxiliary regularized problem with the initial and boundary conditions (1.4)β(1.6):
where and
for any and .
Let
and . It is obvious that is a closed convex set and weakly compact in .
The theory of linear elliptic boundary value problems [15] gives a unique such that
where is dependent on and the norms for and , but not on .
Next, for the unique weak solution to problem (2.10)β(2.11), we consider the following problem:
Lemma 2.2. Under hypotheses (H1)β(H6), there exists one global and unique weak solution to the problem (2.13)β(2.16). Further there are bounds on and which depend on and the known data, but not on . Similar estimates also hold for .
Proof. We begin by choosing a constant such that
and observe that satisfies (2.13)β(2.16) if and only if satisfies
or if and only if satisfies
where
Clearly, . First, we prove the uniqueness of weak solution to the problem (2.13)β(2.16) which is equivalent to (2.18)β(2.21). Let be two weak solutions to the problem (2.18)β(2.21), then satisfies
here
such that
Take as test functions in (2.27), (2.28), respectively. By (2.32) and HΓΆlder inequality,
Thus the uniqueness is established by Gronwall's inequality. We are now in a position to prove the existence result. Define , . Clearly, is a Hilbert space with respect to the scalar product
Set, for ,
The operator is well defined and bounded (because ). To prove the existence result by using [16, Theorem 30.A], it suffices to verify that the operator is hemicontinuous, monotone, and coercive. Note that
for any The hemicontinuity of is easily obtained by the standard method. For the monotone, we first notice that
Hence
By the choice of , we can easily obtain the monotone of . Moreover, from (2.38) we also know that the operator is coercive. Therefore, there exists a unique with such that in and
Especially,
That is, is a weak solution to the problem (2.22)β(2.25). Finally, noting that we can easily establish the bounds on and by the standard energy estimate.
Lemma 2.3. Under hypotheses (H1)β(H6), there exists at least one global weak solution to the problem (2.5)β(2.7), (1.4)β(1.6).
Proof. We define the mapping as
with solution of (2.13)β(2.16). From Lemma 2.2, we know that is well defined and compact. Indeed, lies in a bounded sunsets of and lies in a bounded subset of . Since the injection is compact, we conclude from Aubin's lemma that is relatively compact in . And for given and , holds if we choose large enough. To apply Schauder's fixed point theorem, we still need to prove that the mapping is continuous. Consider any sequence strongly in and let . Since is relatively compact in and bounded in , we can extract subsequences such that
We only have to show . To do this, we only need to prove
where is a test function, and are solutions of (2.10)β(2.11) corresponding to , respectively. The reminder of convergence proof is standard (details see [8] or [9]). Use as test function in a modification of (2.13) in which the functions have been replaced by respectively. Then we have
A similar estimate holds for . Then adding the two inequalities and using strongly in (details see [9]) and (2.42), we conclude that This implies that a.e. in . Then we can easily prove (2.43) by using Vitali's theorem. Now existence of a fixed point of follows which is a solution of (2.5)β(2.7), (1.4)β(1.6).
To obtain the existence result of problem (2.2)β(2.4), (1.4)β(1.6), the following estimates on uniformly in are necessary.
Lemma 2.4. The solutions of problem (2.5)β(2.7), (1.4)β(1.6) satisfy the estimates
where is dependent on and the known data, but not on .
Proof. By taking as test function in (2.6), we have
By taking into account in and the nonnegativity of and we obtain
and thus a.e. in . Similarly, we have a.e. in . To obtain the upper bound set
and use as test function in (2.6), then
Adding the equality (2.49) for and a similar inequality for , we get
Noticing that
and applying HΓΆlder inequality, we further have
where is independent of and . Gronwall's inequality then implies that
for all and . Since the right-hand side of this equality does not depend on , we can let and then to obtain the desired upper bound.
Thus, taking large enough, we see that solves
subject to the initial and boundary conditions (1.4)β(1.6).
Proof of Theorem 2.1. Noticing that the function is bound, we can obtain the following convergence properties by using the same method as the proof in [8, Theorem 1.1] and Lemma 3.2:
In addition, a standard elliptic estimate gives
from which we conclude
and furthermore
Next, using as test function in (2.55), we get
as , where we have used
and (2.57)β(2.60). The same argument shows that as . Thus, there exists a subsequence (not relabeled) such that almost everywhere in as . Then it follows from Vitali's theorem that
Now we can conclude that is the solution of the problem (2.2)β(2.4), (1.4)β(1.6) from the above convergence by standard method and then complete the proof of Theorem 2.1.
3. Proof of the Main Result
In the last section, we prove that there is at least one global weak solution to the problem (2.2)β(2.4), (1.4)β(1.6) for every given . In the following what we need to do is to prove that the limit of is a solution of (1.1)β(1.6). To this end, we first give some uniform estimates for the problem (2.2)β(2.4), (1.4)β(1.6). For simplicity, we drop the subscript of and set
Lemma 3.1. , there holds
where , satisfy
This is the well-known Gagliardo-Nirenberg Inequality [15].
Lemma 3.2. If is the solution of the problem (2.2)β(2.4), (1.4)β(1.6), the following estimate holds:
where is a sufficiently small constant.
Lemma 3.3. If is the solution of the problem (2.2)β(2.4), (1.4)β(1.6), there holds that
for some sufficiently small ; here positive constant is independent of .
Proof. Without loss of generality, we assume and the or case is easier.Case 1. . In this case, from HΓΆlder inequality and (3.3), we have
Noting that , we can choose a constant such that . Then by taking and as test function in (2.3) and (2.4), respectively, and adding them together, we get
We estimate the right-hand side term by term. Due to the equation of and (3.6), we obtain
where we use that and when . and can be bounded as
where we have used the nonnegativity of and . Inserting (3.13)β(3.15) into (3.12) we conclude that
Similarly, taking and as test function in (2.3) and (2.4), we conclude that
Choosing sufficiently small, and then summing (3.16) and (3.17), we have
Now we estimate the term (or ). Due to Lemma 2.2 we get
where satisfy
Consequently we obtain, taking into account (3.11) and the choice of ,
for sufficiently small and such that , where only depends on . This proves that
for some sufficiently small by Lemma 2.2 of the appendix in [14] , and thus
Using the above estimates and (3.4)β(3.6), we can obtain (3.10).Case 2. . In this case, the estimates on are the same as Case 1, and
Then we have
Next, by a small change in (3.17), we obtain
Choose such that
and reset
By a discuss similar to (3.19), we obtain
Finally, noting , from (3.25)β(3.29) and estimate (3.3) we deduce that (3.21) also holds for . Therefore, we can also obtain estimates (3.8)β(3.10) when .
The following lemma is indispensable to prove the estimates of and uniformly in .
Lemma 3.4. If is the solution of the problem (2.2)β(2.4), (1.4)β(1.6), there holds that
Proof. Denote by a Banach space in which function satisfies
The following proof is base on Moser's iteration technique [15]. We insert into (2.3), into (2.4), integrate over the interval and add the equations to obtain
Since is negative, and choosing , we have
for all . Applying HΓΆlder inequality, we obtain
where denotes measure, and . Choose sufficiently small, together with (3.33)β(3.35) and , then
On the other hand, we have
Let . From (3.36) and (3.37), we obtain
which proves (3.30) by [15, Lemma 5.2 in Chapter 2].
Lemma 3.5. If is the solution of the problem (2.2)β(2.4), (1.4)β(1.6), there holds that
Proof. We only prove (3.39). The proof of (3.40) is completely the same as that for (3.39). Take test function for (2.3). Then by Lemma 3.2β3.4 we conclude from
that . The proof is complete.
Proof of Theorem 1.2. By passing to a subsequence if necessary, from Lemma 3.2β3.5, together with compact lemma of Lions [18], we infer that
Then by the same argument as the proof of Theorem 2.1, we can conclude that is the solution of the problem (1.1)β(1.6) from the convergence of (3.42) and then complete the proof of Theorem 1.2.
Acknowledgment
This work was supported by the Science Research Foundation of NUIST (20080295).
References
S. Seldman, Analysis and Simulation of Semiconductor Devices, Springer, New York, NY, USA, 1984.
P. A. Markowich, The Stationary Semiconductor Device Equations, Computational Microelectronics, Springer, Vienna, Austria, 1986.
W. Fang and K. Ito, βOn the time-dependent drift-diffusion model for semiconductors,β Journal of Differential Equations, vol. 117, no. 2, pp. 245β280, 1995.
W. Fang and K. Ito, βGlobal solutions of the time-dependent drift-diffusion semiconductor equations,β Journal of Differential Equations, vol. 123, no. 2, pp. 523β566, 1995.
H. Gajewski, βOn existence, uniqueness and asymptotic behavior of solutions of the basic equations for carrier transport in semiconductors,β Zeitschrift fΓΌr Angewandte Mathematik und Mechanik, vol. 65, no. 2, pp. 101β108, 1985.
A. JΓΌngel, βAsymptotic analysis of a semiconductor model based on Fermi-Dirac statistics,β Mathematical Methods in the Applied Sciences, vol. 19, no. 5, pp. 401β424, 1996.
P. Guan and B. Wu, βExistence of weak solutions to a degenerate time-dependent semiconductor equations with temperature effect,β Journal of Mathematical Analysis and Applications, vol. 332, no. 1, pp. 367β380, 2007.
A. JΓΌngel, βOn the existence and uniqueness of transient solutions of a degenerate nonlinear drift-diffusion model for semiconductors,β Mathematical Models & Methods in Applied Sciences, vol. 4, no. 5, pp. 677β703, 1994.
J. I. DΓaz, G. Galiano, and A. JΓΌngel, βOn a quasilinear degenerate system arising in semiconductors theory. I. Existence and uniqueness of solutions,β Nonlinear Analysis: Real World Applications, vol. 2, no. 3, pp. 305β336, 2001.
J. Frehse and J. Naumann, βStationary semiconductor equations modeling avalanche generation,β Journal of Mathematical Analysis and Applications, vol. 198, no. 3, pp. 685β702, 1996.
J. Frehse and J. Naumann, βOn the existence of weak solutions to a system of stationary semiconductor equations with avalanche generation,β Mathematical Models & Methods in Applied Sciences, vol. 4, no. 2, pp. 273β289, 1994.
J. Naumann, βSemiconductor equations modelling nonstationary avalanche generation in 3D,β Differential and Integral Equations, vol. 7, no. 5-6, pp. 1511β1526, 1994.