Initial Guess

Procedures

Initial guess assumes charge neutrality

For Scharfetter-Gummel (Elec Hole DevPsi) solution

proc InitialGuess {Doping} { 
	sel z= {(Doping>0.0) 
		?  ( 0.025*log( (Doping+1.0e10) / 1.0e10))
		:  (-0.025*log(-(Doping+1.0e10) / 1.0e10))} name = DevPsi
	sel z=1.0e10*exp(DevPsi/0.025) name=Elec
	sel z=1.0e10*exp(-DevPsi/0.025) name=Hole
	sel z=0.0 name=Qfn
	sel z=0.0 name=Qfp
} 

For Quasi-Fermi Solution (Qfn Qfp DevPsi) solution

proc InitialGuessQF {Doping} {
	sel z= {(Doping>0.0) 
		?  ( 0.025*log( (Doping+1.0e10) / 1.0e10))
		:  (-0.025*log(-(Doping+1.0e10) / 1.0e10))} name = DevPsi store
	sel z=1.0e10*exp((DevPsi)/0.025) name=Elec store
	sel z=1.0e10*exp((-DevPsi)/0.025) name=Hole store
	sel z=0.0 name=Qfn store
	sel z=0.0 name=Qfp store
}

For Slotboom (SBn SBp DevPsi) solution

proc InitialGuessSB {Doping} {
	sel z= {(Doping>0.0) 
		?  ( 0.025*log( (Doping+1.0e10) / 1.0e10))
		:  (-0.025*log(-(Doping+1.0e10) / 1.0e10))} name = DevPsi store
	sel z=1.0e10*exp((DevPsi)/0.025) name=Elec 
	sel z=1.0e10*exp((-DevPsi)/0.025) name=Hole 
	sel z=1.0 name=SBn store
	sel z=1.0 name=SBp store
}