octave_packages/secs2d-0.0.8/Utilities/constants.m

   1 % This file is part of
   2 %
   3 %            SECS2D - A 2-D Drift--Diffusion Semiconductor Device Simulator
   4 %         -------------------------------------------------------------------
   5 %            Copyright (C) 2004-2006  Carlo de Falco
   6 %
   7 %
   8 %
   9 %  SECS2D is free software; you can redistribute it and/or modify
  10 %  it under the terms of the GNU General Public License as published by
  11 %  the Free Software Foundation; either version 2 of the License, or
  12 %  (at your option) any later version.
  13 %
  14 %  SECS2D is distributed in the hope that it will be useful,
  15 %  but WITHOUT ANY WARRANTY; without even the implied warranty of
  16 %  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  17 %  GNU General Public License for more details.
  18 %
  19 %  You should have received a copy of the GNU General Public License
  20 %  along with SECS2D; If not, see <http://www.gnu.org/licenses/>.
  21
  22 Kb           = 1.3806503e-23;
  23 q            = 1.602176462e-19;
  24 e0           = 8.854187817e-12;
  25 esir         = 11.7;
  26 esio2r       = 3.9;
  27 esi          = e0 * esir;
  28 esio2        = e0 * esio2r;
  29 hplanck      = 6.626e-34;
  30 hbar         = ( hplanck/ (2*pi));
  31 mn0          = 9.11e-31;
  32 mn           = 0.26*mn0;
  33 mh           = 0.18*mn0;
  34
  35
  36 qsue         = q / esi;
  37 T0           = 300 ;
  38 Vth          = Kb * T0 / q;
  39 un           = 1417e-4;
  40 up           = 480e-4;
  41
  42 vsatn0        = 1.07e5;
  43 vsatp0        = 8.37e4;
  44 vsatnexp        = 0.87;
  45 vsatpexp        = 0.52;
  46 vsatn           = vsatn0*(300/T0).^vsatnexp;
  47 vsatp           = vsatp0*(300/T0).^vsatpexp;
  48
  49 mubn0         = 1.109;
  50 mubp0         = 1.213;
  51 mubnexp       = 0.66;
  52 mubpexp       = 0.17;
  53 mubn          = mubn0*(T0/300)^mubnexp;
  54 mubp          = mubp0*(T0/300)^mubpexp;
  55 mudopparn     = [   52.2e-4 %mumin1
  56                     52.2e-4 %mumin2
  57                     43.4e-4 %mu1
  58                     0e+6    %Pc
  59                     9.68e22 %Cr
  60                     3.43e26 %Cs
  61                     0.680   %alpha
  62                     2.0     %beta
  63                 ];
  64 mudopparp     = [   44.9e-4 %mumin1
  65                     0.00e-4 %mumin2
  66                     29.0e-4 %mu1
  67                     9.23e22 %Pc
  68                     2.23e23 %Cr
  69                     6.10e26 %Cs
  70                     0.719   %alpha
  71                     2.0     %beta
  72                 ];
  73
  74
  75 tp           = 1e-7;
  76 tn           = 1e-7;
  77 twp          = .1e-12;
  78 twn          = .1e-12;
  79
  80
  81 mnl          = 0.98*mn0;
  82 mnt          = 0.19*mn0;
  83 mndos        = (mnl*mnt*mnt)^(1/3);
  84
  85 mhh             = 0.49*mn0;
  86 mlh             = 0.16*mn0;
  87 mhdos           = (mhh^(3/2)+mlh^(3/2))^(2/3);
  88
  89 rn              = 3;
  90 aleph           = hbar^2/(4*rn*q*mn);
  91 alephn          = aleph;
  92 rp              = .1;
  93 alephp          = hbar^2/(4*rp*q*mh);
  94
  95 Nc              = (6/4)*(2*mndos*Kb*T0/(hbar^2*pi))^(3/2);
  96 Nv              = (1/4)*(2*mhdos*Kb*T0/(hbar^2*pi))^(3/2);
  97 Eg0             = 1.16964*q;
  98 alfaEg          = 4.73e-4*q;
  99 betaEg          = 6.36e2;
 100 Egap            = Eg0-alfaEg*((T0^2)/(T0+betaEg));
 101
 102 ni              = sqrt(Nc*Nv)*exp(-Egap/(2*(Kb * T0)));
 103 Phims           = -.79;%- Egap /(2*q);
 104
 105 kappaSi         = 154.86;
 106
 107 % Last Revision:
 108 % $Author: adb014 $
 109 % $Date: 2008-02-04 16:26:27 +0100 (man, 04 feb 2008) $
 110
 111