blob: bf0d8d95e8e5c53f631ca5ff7a861c1d31cfee4f (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
//Example 2_8_u2
clc();
clear;
//TO Find the energy at which probability of occupancy is 0.9 and density of states and the population density
fe=0.9
k=(1/fe)-1
logk=log(k)
kb=8.62*10^-5
t=1000 //units in K
E=logk*kb*t //units in eV
ef=7.06 //units in eV
energy=E+ef //units in eV
printf("The energy of this state is E=%.2f eV",energy)
n=9.1*10^-31
EE=energy //units in eV
h=4.14*10^-15
ZE=(8*sqrt(2)*%pi*n^1.5*sqrt(EE))/h^3 //units in met^-3 (eV)^-1
ZE=ZE*1.56*10^28 //units met^-3 (eV)^-1
printf("\nThe density of the states for this energy is Z(E)=")
disp(ZE)
printf("met^-3 (eV)^-1")
ne=ZE*fe //units in met^-3 (eV)^-1
printf("\nThe population density for this energy is N(E)=")
disp(ne)
printf("met^-3 (eV)^-1")
|