Add a useful package (from Source forge) for octave

[CreaPhase.git] / octave_packages / secs2d-0.0.8 / Utilities / constants.m

diff --git a/octave_packages/secs2d-0.0.8/Utilities/constants.m b/octave_packages/secs2d-0.0.8/Utilities/constants.m
new file mode 100644 (file)
index 0000000..72be361
--- /dev/null
+++ b/octave_packages/secs2d-0.0.8/Utilities/constants.m
@@ -0,0 +1,111 @@
+% This file is part of 
+%
+%            SECS2D - A 2-D Drift--Diffusion Semiconductor Device Simulator
+%         -------------------------------------------------------------------
+%            Copyright (C) 2004-2006  Carlo de Falco
+%
+%
+%
+%  SECS2D is free software; you can redistribute it and/or modify
+%  it under the terms of the GNU General Public License as published by
+%  the Free Software Foundation; either version 2 of the License, or
+%  (at your option) any later version.
+%
+%  SECS2D is distributed in the hope that it will be useful,
+%  but WITHOUT ANY WARRANTY; without even the implied warranty of
+%  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+%  GNU General Public License for more details.
+%
+%  You should have received a copy of the GNU General Public License
+%  along with SECS2D; If not, see <http://www.gnu.org/licenses/>.
+
+Kb           = 1.3806503e-23;
+q            = 1.602176462e-19;
+e0           = 8.854187817e-12;
+esir        = 11.7;
+esio2r              = 3.9;
+esi         = e0 * esir;
+esio2       = e0 * esio2r;
+hplanck             = 6.626e-34;
+hbar         = ( hplanck/ (2*pi));
+mn0          = 9.11e-31;
+mn           = 0.26*mn0;
+mh           = 0.18*mn0;
+
+
+qsue         = q / esi;
+T0           = 300 ;
+Vth         = Kb * T0 / q;
+un           = 1417e-4;
+up           = 480e-4;
+
+vsatn0        = 1.07e5;
+vsatp0        = 8.37e4;
+vsatnexp        = 0.87;
+vsatpexp        = 0.52;
+vsatn           = vsatn0*(300/T0).^vsatnexp;
+vsatp           = vsatp0*(300/T0).^vsatpexp;
+
+mubn0         = 1.109;
+mubp0         = 1.213;
+mubnexp       = 0.66;
+mubpexp       = 0.17;
+mubn          = mubn0*(T0/300)^mubnexp;
+mubp          = mubp0*(T0/300)^mubpexp;
+mudopparn     = [   52.2e-4 %mumin1
+                    52.2e-4 %mumin2
+                    43.4e-4 %mu1
+                    0e+6    %Pc
+                    9.68e22 %Cr
+                    3.43e26 %Cs
+                    0.680   %alpha
+                    2.0     %beta
+                ];
+mudopparp     = [   44.9e-4 %mumin1
+                    0.00e-4 %mumin2
+                    29.0e-4 %mu1
+                    9.23e22 %Pc
+                    2.23e23 %Cr
+                    6.10e26 %Cs
+                    0.719   %alpha
+                    2.0     %beta
+                ];
+                    
+
+tp           = 1e-7;
+tn           = 1e-7;
+twp          = .1e-12;
+twn          = .1e-12;
+
+
+mnl          = 0.98*mn0;
+mnt          = 0.19*mn0;
+mndos        = (mnl*mnt*mnt)^(1/3); 
+
+mhh             = 0.49*mn0;
+mlh             = 0.16*mn0;
+mhdos           = (mhh^(3/2)+mlh^(3/2))^(2/3);
+
+rn              = 3;
+aleph           = hbar^2/(4*rn*q*mn);
+alephn          = aleph;
+rp              = .1;
+alephp          = hbar^2/(4*rp*q*mh);
+
+Nc              = (6/4)*(2*mndos*Kb*T0/(hbar^2*pi))^(3/2);   
+Nv              = (1/4)*(2*mhdos*Kb*T0/(hbar^2*pi))^(3/2);
+Eg0             = 1.16964*q;
+alfaEg          = 4.73e-4*q;
+betaEg          = 6.36e2;
+Egap            = Eg0-alfaEg*((T0^2)/(T0+betaEg));
+
+ni              = sqrt(Nc*Nv)*exp(-Egap/(2*(Kb * T0)));
+Phims           = -.79;%- Egap /(2*q);
+
+kappaSi         = 154.86;
+
+% Last Revision:
+% $Author: adb014 $
+% $Date: 2008-02-04 16:26:27 +0100 (man, 04 feb 2008) $
+
+