1## Copyright (C) 2004-2008  Carlo de Falco
2##
3## SECS1D - A 1-D Drift--Diffusion Semiconductor Device Simulator
4##
5## SECS1D is free software; you can redistribute it and/or modify
6## it under the terms of the GNU General Public License as published by
7## the Free Software Foundation; either version 2 of the License, or
8## (at your option) any later version.
9##
10## SECS1D is distributed in the hope that it will be useful,
11## but WITHOUT ANY WARRANTY; without even the implied warranty of
12## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
13## GNU General Public License for more details.
14##
15## You should have received a copy of the GNU General Public License
16## along with SECS1D; If not, see <http://www.gnu.org/licenses/>.
17##
18## author: Carlo de Falco <cdf _AT_ users.sourceforge.net>
19
20## -*- texinfo -*-
21##
22## @deftypefn {Script File} constants
23##
24## Compute global constants needed for Drift-Diffusion simulation
25##
26## @end deftypefn
27
28
29
30Kb           = 1.3806503e-23;
31q            = 1.602176462e-19;
32e0           = 8.854187817e-12;
33esir 	     = 11.7;
34esio2r 	     = 3.9;
35esi 	     = e0 * esir;
36esio2 	     = e0 * esio2r;
37hplanck	     = 6.626e-34;
38hbar         = ( hplanck/ (2*pi));
39mn0          = 9.11e-31;
40mn           = 0.26*mn0;
41mh           = 0.18*mn0;
42
43
44qsue         = q / esi;
45T0           = 300 ;
46Vth 	     = Kb * T0 / q;
47un           = 1417e-4;
48up           = 480e-4;
49tp           = 1e-7;
50tn           = 1e-7;
51
52mnl          = 0.98*mn0;
53mnt          = 0.19*mn0;
54mndos        = (mnl*mnt*mnt)^(1/3);
55
56mhh             = 0.49*mn0;
57mlh             = 0.16*mn0;
58mhdos           = (mhh^(3/2)+mlh^(3/2))^(2/3);
59
60rn              = .1;
61aleph           = hbar^2/(4*rn*q*mn);
62alephn          = aleph;
63rp              = .1;
64alephp          = hbar^2/(4*rp*q*mh);
65
66Nc              = (6/4)*(2*mndos*Kb*T0/(hbar^2*pi))^(3/2);
67Nv              = (1/4)*(2*mhdos*Kb*T0/(hbar^2*pi))^(3/2);
68Eg0             = 1.16964*q;
69alfaEg          = 4.73e-4*q;
70betaEg          = 6.36e2;
71Egap            = Eg0-alfaEg*((T0^2)/(T0+betaEg));
72
73ni              = sqrt(Nc*Nv)*exp(-Egap/(2*(Kb * T0)));
74Phims           = - Egap /(2*q);
75