Can anyone please help me in correcting the code for ns vs thickness AlN

I want to plot sheet charge density 'ns' vs AlN spacer thickness 'tAlN ' using eq. 9 but I am not getting correct results as I should get the plot which I have attached . Please check my code and modify it so that I can understand my mistakes. As in this eqn 'tAlN ' will vary but here used tAlN for 'detaEceff' and ns both.

clear all;
q=1.6e-19;
E0=8.85e-12;
Ealgan=10.31;
Ealn=10.78;
fib=(5.1*1.6e-19);
sigmaaln=3.38e17;
sigmaalgan=1.3e17;
detaec=(0.862*1.6e-19);
talgan=23e-9;
d0=26e-9;
p=(q^2/Ealn*E0);
taln=0:0.1:2.5;
m=length(taln);
for i=1:m
    detaec2(i)=[detaec+(p*sigmaaln*taln(i))];
end
ns(i)=(sigmaaln*taln-(Ealgan*E0*fib/q)+(Ealn*E0*detaec2(i)/q^2))/(talgan+taln+d0);
plot(taln, ns)

Answer 1

[Abderrahim. B]

Hi!

Please learn how to use for loops ...

doc for

The code below plots a graph similar to the one you attached, bmaybe still need some to be reviewed.

clear;

q=1.6e-19;
E0=8.85e-12;
Ealgan=10.31;
Ealn=10.78;
fib=(5.1*1.6e-19);
sigmaaln=3.38e17;
sigmaalgan=1.3e17;
detaec=(0.862*1.6e-19);
talgan=23e-9;
d0=26e-9;
p=((q^2)/Ealn)*E0;
taln=10:10:100;
m=length(taln);


for i=1:m
    detaec2 = detaec + (p*sigmaaln*taln(i));
    ns(i) = ((sigmaaln*taln(i) - (Ealgan*E0)*(fib/q) + ((Ealn*E0*detaec2)/(q^2)) ))/(talgan +taln(i) +d0);
end

plot(taln, ns, 'k-o')

[Nudrat Sufiyan]

I am very much thankful to you for your response as here you have taken tAlN range from 10to 100 but I need tAlN from 0.8 to 10 nm. How I'll improve this

[Abderrahim. B]

clear;
format long
q=1.6e-19;
E0=8.85e-12;
Ealgan=10.31;
Ealn=10.78;
fib=(5.1*1.6e-19);
sigmaaln=3.38e17;
sigmaalgan=1.3e17;
detaec=(0.862*1.6e-19);
talgan=23e-9;
d0=26e-9;
p=((q^2)/Ealn)*E0;
taln=0.8:0.1:10;
m=length(taln);
for i=1:m
    detaec2 = detaec + (p*sigmaaln*taln(i));
    ns(i) = ((sigmaaln*taln(i) - (Ealgan*E0)*(fib/q) + ((Ealn*E0*detaec2)/(q^2)) ))/(talgan +taln(i) +d0);
end
plot(taln, ns, 'k-o')
xlabel('taln (nm)')

[Nudrat Sufiyan]

Thankyou sir. This is really helpful for me.

[Abderrahim. B]

My pleasure!