1 # Created by Octave 3.6.1, Thu Mar 22 16:18:37 2012 UTC <root@t61>
10 DDGelectron_driftdiffusion
13 # name: <cell-element>
17 n=DDGelectron_driftdiffusion(psi,x,ng,p)
18 Solves the continuity equation for electrons
19 input: psi electric potential
21 ng initial guess and BCs for electron density
22 p hole density (for SRH recombination)
23 output: n updated electron density
27 # name: <cell-element>
31 n=DDGelectron_driftdiffusion(psi,x,ng,p)
32 Solves the continuity equation fo
36 # name: <cell-element>
43 # name: <cell-element>
49 DDGgummelmap (x,idata,toll,maxit,ptoll,pmaxit,verbose)
51 Solves the scaled stationary bipolar DD equation system
52 using Gummel algorithm
55 idata.D doping profile
56 idata.p initial guess for hole concentration
57 idata.n initial guess for electron concentration
58 idata.V initial guess for electrostatic potential
59 idata.Fn initial guess for electron Fermi potential
60 idata.Fp initial guess for hole Fermi potential
61 idata.l2 scaled electric permittivity (diffusion coefficient in Poisson equation)
62 idata.un scaled electron mobility
63 idata.up scaled electron mobility
64 idata.nis scaled intrinsic carrier density
65 idata.tn scaled electron lifetime
66 idata.tp scaled hole lifetime
67 toll tolerance for Gummel iterarion convergence test
68 maxit maximum number of Gummel iterarions
69 ptoll tolerance for Newton iterarion convergence test for non linear Poisson
70 pmaxit maximum number of Newton iterarions
71 verbose verbosity level: 0,1,2
73 output: odata.n electron concentration
74 odata.p hole concentration
75 odata.V electrostatic potential
76 odata.Fn electron Fermi potential
77 odata.Fp hole Fermi potential
78 it number of Gummel iterations performed
79 res total potential increment at each step
83 # name: <cell-element>
92 # name: <cell-element>
96 DDGhole_driftdiffusion
99 # name: <cell-element>
103 p=DDGhole_driftdiffusion(psi,x,pg,n)
104 Solves the continuity equation for holes
105 input: psi electric potential
107 pg initial guess and BCs for hole density
108 n electron density (to compute SRH recombination)
109 output: p updated hole density
113 # name: <cell-element>
117 p=DDGhole_driftdiffusion(psi,x,pg,n)
118 Solves the continuity equation for ho
122 # name: <cell-element>
129 # name: <cell-element>
133 phin = DDGn2phin (V,n);
134 computes the qfl for electrons using Maxwell-Boltzmann
139 # name: <cell-element>
143 phin = DDGn2phin (V,n);
144 computes the qfl for electrons using Maxwell-B
148 # name: <cell-element>
155 # name: <cell-element>
159 [V,n,p,res,niter] = DDGnlpoisson (x,sinodes,Vin,nin,...
160 pin,Fnin,Fpin,D,l2,toll,maxit,verbose)
161 Solves the non linear Poisson equation
162 $$ - lamda^2 *V'' + (n(V,Fn) - p(V,Fp) -D)=0 $$
163 input: x spatial grid
164 sinodes index of the nodes of the grid which are in the
165 semiconductor subdomain
166 (remaining nodes are assumed to be in the oxide subdomain)
167 Vin initial guess for the electrostatic potential
168 nin initial guess for electron concentration
169 pin initial guess for hole concentration
170 Fnin initial guess for electron Fermi potential
171 Fpin initial guess for hole Fermi potential
173 l2 scaled electric permittivity (diffusion coefficient)
174 toll tolerance for convergence test
175 maxit maximum number of Newton iterations
176 verbose verbosity level: 0,1,2
177 output: V electrostatic potential
178 n electron concentration
180 res residual norm at each step
181 niter number of Newton iterations
185 # name: <cell-element>
189 [V,n,p,res,niter] = DDGnlpoisson (x,sinodes,Vin,nin,.
193 # name: <cell-element>
200 # name: <cell-element>
204 phip = DDGp2phip (V,p);
205 computes the qfl for holes using Maxwell-Boltzmann statistics
209 # name: <cell-element>
213 phip = DDGp2phip (V,p);
214 computes the qfl for holes using Maxwell-Boltzmann sta
218 # name: <cell-element>
225 # name: <cell-element>
229 n = DDGphin2n (V,phin);
230 computes the electron density using Maxwell-Boltzmann
235 # name: <cell-element>
239 n = DDGphin2n (V,phin);
240 computes the electron density using Maxwell-Bo
244 # name: <cell-element>
251 # name: <cell-element>
255 p = DDGphip2p (V,phip);
256 computes the hole density using Maxwell-Boltzmann
261 # name: <cell-element>
265 p = DDGphip2p (V,phip);
266 computes the hole density using Maxwell-Boltzm
270 # name: <cell-element>
277 # name: <cell-element>
281 DDGplotresults(x,n,p,V,Fn,Fp);
285 # name: <cell-element>
289 DDGplotresults(x,n,p,V,Fn,Fp);