Arti run this code and you will see the graph:
Nde = 2.7*10^19;
Nsdp = 1*10^20;
lateralstraggle = 2*10^-9;
seff = log(Nde/Nsdp)*(-2*(lateralstraggle)^2);
Lg=70*10^-9;
Leff=Lg-(2*seff);
Deff=seff+Leff;
k=1;
T=20;
q=1.6*10^-19;
Vt=(k*T)/q;
ni=1;
Eg=1.17*1.6*10^-19;
nieff=(ni^2*exp(Eg/k*T))^1/2;
Nde=2.7*10^19;
Na=56;
Vbi=Vt*log(Nde*Na/(nieff)^2);
epsilon=32;
epsilonox=4;
tsi=20*10^-9;
tox=1.2*10^-9;
lambda=(epsilon*tsi*tox)/(2*epsilonox);
vgs=2;
vfb=1;
q=1.6*10^-19;
Na=56;
Nsdp= 1*10^20;
x=4.6*(1.6*10^-19);
lateralstraggle=7:70; % Here many values for lateralstraggle to see multiple values on graph
Nsdx=Nsdp*exp((-x^2)./(2*lateralstraggle.^2));
Lg=4;
Sd=4;
Nde=2.7*10^19;
Ei0=0.054*(1.6*10^-19);
Ei=Ei0*(1-(Nsdx/Nde).^1/3);
Eg=34;
eg=24;
Egeff=Eg-eg;
Ed= Egeff-Ei;
Ef=8;
Eg=3;
ni=8;
nieff=(ni^2*exp(Eg/k*T))^1/2;
Nsdx=Nsdp*exp((-x^2)./(2*lateralstraggle.^2))+Nsdp*exp((-(Lg-x)/2*lateralstraggle.^2))/(1+Sd*exp((Ef-Ed)/k*T));
PI=(vgs-vfb)-lambda^2*q*(Na-Nsdx)/epsilon;
Vds=1.3607e+003;
c2=((Vbi-PI)*(1-(exp(-Deff/lambda))/(exp(-seff/lambda))+Vds))/(exp(Deff/lambda))-(exp(seff/lambda))/(exp(-seff/lambda))*(exp(-Deff/lambda));
c1=(Vbi-c2*exp(seff/lambda)/exp(-seff/lambda))*exp(-Deff/lambda);
si=c1*exp(-x/lambda)+c2*exp(x/lambda)+PI;
Ni=2;
fi=Vbi*log(Na/Ni);
Na=3;
Vth=vfb+(2*fi)-(c1*exp(-x/lambda))-(c2*exp(x/lambda))+(lambda^2)*q*(Na-Nsdx)/epsilon;
plot(lateralstraggle,Vth)
ylabel('Vth, V')
xlabel('lateralstraggle,nm')
Is it what you are looking for?
