blob: 15df6a9d918c1acb76ae4dad4493f6867ee07880 (
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
26
|
// Scilab Code Ex4.9a: Page-141 (2006)
clc; clear;
E_F = 1; // For simplicity assume Fermi energy to be unity, eV
k = 1.38e-023; // Boltzmann constant, J/mol/K
e = 1.6e-019; // Energy equivalent of 1 eV, J/eV
dE = 0.2; // Exces energy above Fermi level, eV
T = 0+273; // Room temperature, K
E = E_F + dE; // Energy of the level above Fermi level, eV
f_E = 1/(exp((E-E_F)*e/(k*T))+1); // Occupation probability of the electron at 0.1 eV above E_F
printf("\nAt 273 K:");
printf("\n=========");
printf("\nThe occupation probability of electron at %3.1f eV above Fermi energy = %4.2e", dE, f_E);
T = 100+273; // Given temperature of 100 degree celsius, K
f_E = 1/(exp((E-E_F)*e/(k*T))+1); // Occupation probability of the electron at 0.1 eV below E_F
printf("\n\nAt 373 K:");
printf("\n=========");
printf("\nThe occupation probability of electron at %3.1f eV above Fermi energy = %4.2e", dE, f_E);
// Result
// At 273 K:
// =========
// The occupation probability of electron at 0.2 eV above Fermi energy = 2.05e-004
// At 373 K:
// =========
// The occupation probability of electron at 0.2 eV above Fermi energy = 1.99e-003
|
