diff options
Diffstat (limited to '1580/CH3')
-rwxr-xr-x | 1580/CH3/EX3.1/Ch03Ex1.sce | 18 | ||||
-rwxr-xr-x | 1580/CH3/EX3.10/Ch03Ex10.sce | 11 | ||||
-rwxr-xr-x | 1580/CH3/EX3.11/Ch03Ex11.sce | 18 | ||||
-rwxr-xr-x | 1580/CH3/EX3.2/Ch03Ex2.sce | 9 | ||||
-rwxr-xr-x | 1580/CH3/EX3.3/Ch03Ex3.sce | 15 | ||||
-rwxr-xr-x | 1580/CH3/EX3.4/Ch03Ex4.sce | 16 | ||||
-rwxr-xr-x | 1580/CH3/EX3.5/Ch03Ex5.sce | 14 | ||||
-rwxr-xr-x | 1580/CH3/EX3.6/Ch03Ex6.sce | 13 | ||||
-rwxr-xr-x | 1580/CH3/EX3.7/Ch03Ex7.sce | 14 | ||||
-rwxr-xr-x | 1580/CH3/EX3.8/Ch03Ex8.sce | 15 | ||||
-rwxr-xr-x | 1580/CH3/EX3.9/Ch03Ex9.sce | 14 |
11 files changed, 157 insertions, 0 deletions
diff --git a/1580/CH3/EX3.1/Ch03Ex1.sce b/1580/CH3/EX3.1/Ch03Ex1.sce new file mode 100755 index 000000000..5744f67b1 --- /dev/null +++ b/1580/CH3/EX3.1/Ch03Ex1.sce @@ -0,0 +1,18 @@ +// Scilab Code Ex3.1 : Page-3.4 (2004)
+ // In a SC structure number of planes are having three arrangement (100),(110) and (111)
+ clc;clear;
+a = 1; // For simplicity lattice constant is taken to be unity
+A_100 = a^2; // Area of the plane (100), mm^2
+N_100 = 1/A_100; // Number of atoms along (100) plane, atoms per square mm
+A_110 = sqrt(2)*a^2; // Area of the plane (110), mm^2
+N_110 = 1/A_110; // // Number of atoms along (110) plane, atoms per square mm
+A_111 = 1/2*a*sqrt(2)*sqrt(2)*a^2*cosd(30); // Area of the plane (110), mm^2
+A_111t = 0.5; // Total no of atoms in (111) plane
+N_111 = A_111t/A_111; // // Number of atoms along (110) plane, atoms per square mm
+printf("\nNumber of atoms along (100) plane= %d /a^2 atoms per square mm", N_100);
+printf("\nNumber of atoms along (110) plane= %f atoms per square mm", N_110);
+printf("\nNumber of atoms along (111) plane= %5.3f /a^2 atoms per square mm", N_111);
+// Result
+// Number of atoms along (100) plane= 1 /a^2 atoms per square mm
+// Number of atoms along (110) plane= 0.707107 atoms per square mm
+// Number of atoms along (111) plane= 0.577 /a^2 atoms per square mm
diff --git a/1580/CH3/EX3.10/Ch03Ex10.sce b/1580/CH3/EX3.10/Ch03Ex10.sce new file mode 100755 index 000000000..8ce6694d0 --- /dev/null +++ b/1580/CH3/EX3.10/Ch03Ex10.sce @@ -0,0 +1,11 @@ +// Scilab Code Ex3.10 :Page-3.16 (2004)
+T = 273+25; // Temperature , K
+r = 2.82e-10; // Interionic distance, m
+N = 4/((2*r)^3); // Density of ion pairs, ion pairs
+k = 8.625e-5; // Boltzmann constant, J/K
+n = 5e+11; // Density od Schottky effects, per unit volume
+E_s = 2*k*T*2.303*log10(N/n); // Average energy required to creat Schottky defect
+printf("\nAverage energy required to create Schottky defect = %5.3f eV", E_s);
+
+// Result
+// Average energy required to create Schottky defect = 1.971 eV
diff --git a/1580/CH3/EX3.11/Ch03Ex11.sce b/1580/CH3/EX3.11/Ch03Ex11.sce new file mode 100755 index 000000000..7e73fe5f0 --- /dev/null +++ b/1580/CH3/EX3.11/Ch03Ex11.sce @@ -0,0 +1,18 @@ +// Scilab Code Ex3.11 : Ratio of vacancies in metal to create Frenkel defect:Page-3.18 (2004)
+N = 1; // For simplicity assume total number of metal ions to be unity
+Ni = 1; // For simplicity assume total number of metal ions to be unity
+k = 8.625e-5; // Boltzmann constant, J/K
+T1 = 273+20; // First temperature for metal, K
+T2 = 300+273; // Second temperature for metal, K
+E_v = 1.4; // Average energy required to create a vacancy in metal, eV
+n_293 = N*exp(-E_v/(2*k*T1)); // Number of vacancies at 500 K
+n_573 = N*exp(-E_v/(2*k*T2)); // Number of vacancies at 500 K
+n_ratio1 = n_573/n_293; // Ratio of vacancies in metal
+n_ratio2 = n_293/n_573; // Ratio of vacancies in metal
+
+printf("\nThe ratio 1 of vacancies in metal to create Frenkel defect = %5.3e", n_ratio1);
+printf("\nThe ratio 2 of vacancies in metal to create Frenkel defect = %5.3e", n_ratio2);
+
+// Result
+// The ratio 1 of vacancies in metal to create Frenkel defect = 7.558e+05
+// The ratio 2 of vacancies in metal to create Frenkel defect = 1.323e-06
diff --git a/1580/CH3/EX3.2/Ch03Ex2.sce b/1580/CH3/EX3.2/Ch03Ex2.sce new file mode 100755 index 000000000..d237f0ef2 --- /dev/null +++ b/1580/CH3/EX3.2/Ch03Ex2.sce @@ -0,0 +1,9 @@ +// Scilab Code Ex3.2 : Page-3.5(2004)
+clc;clear;
+r = 1; // For simplicity assume radius of atom to be unity, unit
+a = 4*r/sqrt(3); // Lattice constant, unit
+R = (a/2)-r; // R be the radius of interstitial sphere that can fit into void, unit
+printf ("\nMaximum Radius of sphere that can fit into BCC = %5.3fr", R);
+
+// Result
+// Maximum Radius of sphere that can fit into BCC = 0.155r
diff --git a/1580/CH3/EX3.3/Ch03Ex3.sce b/1580/CH3/EX3.3/Ch03Ex3.sce new file mode 100755 index 000000000..2ec85cbfb --- /dev/null +++ b/1580/CH3/EX3.3/Ch03Ex3.sce @@ -0,0 +1,15 @@ +// Scilab Code Ex3.3 : Page-3.6 (2010)
+clc;clear;
+r1 = 1.258e-10; // Atomic radius in BCC, metre
+a1 = 4*r1/sqrt(3); // Lattice constant for BCC, metre
+V1 = a1^3; // Volume of unit cell in BCC, metre cube
+Vpa = V1/2; // Volume occupied by one atom in BCC, metre cube
+r2 = 1.292e-10; // Atomic radius in FCC, metre
+a2 = 2*r2*sqrt(2); // Lattice constant for F CC, cube
+V2 = a2^3; // Volume of unit cell in FCC, meter cube
+Vpa1 = V2/4; // Volume occupied by one atom in FCC, metre cube
+dV = (Vpa-Vpa1)/Vpa*100; // Change in volume, percentage
+printf("\nChange in volume in percentage = %4.3f percentage", dV);
+
+// Result
+// Change in volume in percentage = 0.493 percentage
diff --git a/1580/CH3/EX3.4/Ch03Ex4.sce b/1580/CH3/EX3.4/Ch03Ex4.sce new file mode 100755 index 000000000..846e4555c --- /dev/null +++ b/1580/CH3/EX3.4/Ch03Ex4.sce @@ -0,0 +1,16 @@ +// Scilab Code Ex3.4 : Page-3.7 (2010)
+clc;clear;
+a = 0.27e-9; // Lattice constant for BCC, metre
+c = 0.494e-9; // Height of the unit cell, metre
+M = 65.37; // Atomic weight of zn, kg
+N = 6.02e+26; // Avogadro number per k mol
+m = 6*M/N; // Mass per unit cell in HCP structure, kg
+V = 3*sqrt(3)*a^2*c/2; // Volume of unit cell in HCP, metre cube
+rho = m/V; // Density of HCP Zn structure, kg per metrecube
+
+printf("\nVolume of HCP Zn structure = %4.3e metrecube", V);
+printf("\nDensity of HCP Zn structure = %4.0f kg per metrecube", rho);
+
+// Result
+// Volume of HCP Zn structure = 9.356e-29 metrecube
+// Density of HCP Zn structure = 6963 kg per metrecube
diff --git a/1580/CH3/EX3.5/Ch03Ex5.sce b/1580/CH3/EX3.5/Ch03Ex5.sce new file mode 100755 index 000000000..08973effd --- /dev/null +++ b/1580/CH3/EX3.5/Ch03Ex5.sce @@ -0,0 +1,14 @@ +/// Scilab Code Ex3.5 : Page-3.5 (2004)
+clc;clear;
+r = 0.1278; // Atomic radius, nm
+a = 4*r/sqrt(2); // Lattice constant, nm
+h1 = 1, k1 = 1, l1 = 0; // Miller Indices of (110) planes
+d_110 = a/sqrt(h1^2 + k1^2 + l1^2); // Interplanar spacing for (110) planes, nm
+h2 = 2, k2 = 1, l2 = 2; // Indices of third set of parallel planes
+d_212 = a/sqrt(h2^2 + k2^2 + l2^2); // Interplanar spacing for (111) planes, nm
+printf("\nInterplanar spacing for (110) planes = %6.4f nm", d_110);
+printf("\nInterplanar spacing for (212) planes = %6.4f nm", d_212);
+
+// Result
+// Interplanar spacing for (110) planes = 0.2556 nm
+// Interplanar spacing for (212) planes = 0.1205 nm
diff --git a/1580/CH3/EX3.6/Ch03Ex6.sce b/1580/CH3/EX3.6/Ch03Ex6.sce new file mode 100755 index 000000000..123fe17c1 --- /dev/null +++ b/1580/CH3/EX3.6/Ch03Ex6.sce @@ -0,0 +1,13 @@ +// Scilab Code Ex3.6 : Page-3.8 (2004)
+clc;clear;
+a = 1; // For simplicity we assume a to be unity, unit
+h1 = 1, k1 = 0, l1 = 0; // Indices of first set of parallel planes
+d_100 = a/sqrt(h1^2 + k1^2 + l1^2); // Interplanar spacing for (100) planes, unit
+h2 = 1, k2 = 1, l2 = 0; // Indices of second set of parallel planes
+d_110 = a/sqrt(h2^2 + k2^2 + l2^2); // Interplanar spacing for (110) planes, unit
+h3 = 1, k3 = 1, l3 = 1; // Indices of third set of parallel planes
+d_111 = a/sqrt(h3^2 + k3^2 + l3^2); // Interplanar spacing for (111) planes, unit
+printf("\nd_100 : d_110 : d_111 = %1d : %4.2f : %4.2f", d_100, d_110, d_111);
+
+// Result
+// d_100 : d_110 : d_111 = 1 : 0.71 : 0.58
diff --git a/1580/CH3/EX3.7/Ch03Ex7.sce b/1580/CH3/EX3.7/Ch03Ex7.sce new file mode 100755 index 000000000..e375812e8 --- /dev/null +++ b/1580/CH3/EX3.7/Ch03Ex7.sce @@ -0,0 +1,14 @@ +// Scilab Code Ex3.7 : Page-3.8 (2004)
+clc;clear;
+m = 1; n = 1/2; p = 3; // Coefficients of intercepts along three axes
+m_inv = 1/m; // Reciprocate the first coefficient
+n_inv = 1/n; // Reciprocate the second coefficient
+p_inv = 1/p; // Reciprocate the third coefficient
+mul_fact = double(lcm(int32([1, 1, 3]))); // Find l.c.m. of 1, 1 and 3
+m1 = m_inv*mul_fact; // Clear the first fraction
+m2 = n_inv*mul_fact; // Clear the second fraction
+m3 = p_inv*mul_fact; // Clear the third fraction
+printf("\nThe required miller indices are : (%d %d %d) ", m1,m2,m3);
+
+// Result
+// The required miller indices are : (3 6 1)
diff --git a/1580/CH3/EX3.8/Ch03Ex8.sce b/1580/CH3/EX3.8/Ch03Ex8.sce new file mode 100755 index 000000000..ef639e9b4 --- /dev/null +++ b/1580/CH3/EX3.8/Ch03Ex8.sce @@ -0,0 +1,15 @@ +// Scilab Code Ex3.8 : Page-3.13 (2004)
+clc;clear;
+N = 1; // For simplicity assume total number of metal ions to be unity
+e = 1.6e-019; // Electronic charge, C
+k = 1.38e-023/e; // Boltzmann constant, eV/K
+T1 = 500; // First temperature for metal, K
+T2 = 1000; // Second temperature for metal, K
+E_v = 1; // Average energy required to create a vacancy in metal, eV
+n_500 = N*exp(-E_v/(k*T1)); // Number of vacancies at 500 K
+n_1000 = N*exp(-E_v/(k*T2)); // Number of vacancies at 500 K
+n_ratio = n_1000/n_500; // Ratio of vacancies in metal
+printf("\nThe ratio of vacancies in metal = %5.3e", n_ratio);
+
+// Result
+// The ratio of vacancies in metal = 1.085e+05
diff --git a/1580/CH3/EX3.9/Ch03Ex9.sce b/1580/CH3/EX3.9/Ch03Ex9.sce new file mode 100755 index 000000000..f8eb30ab7 --- /dev/null +++ b/1580/CH3/EX3.9/Ch03Ex9.sce @@ -0,0 +1,14 @@ +// Scilab Code Ex3.9 : Page-3.14 (2004)
+clc;clear;
+T1 = 500+273; // First temperature for metal, K
+T2 = 1000+273; // Second temperature for metal, K
+frac_vac = 1e-010; // n1/N, the fraction of vacancy sites at 500 degree celsius
+e = 1.6e-019; // Electronic charge, C
+k = 1.38e-023/e; // Boltzmann constant, eV/K
+// n1 = N*exp(-E_f/(k*T1)); // Number of vacancies at 500 K
+// n2 = N*exp(-E_f/(k*T2)); // Number of vacancies at 500 K, solving for n2/N = x
+x = exp((T1/T2)*log(frac_vac));
+printf("\nThe fraction of vacancy sites in metal = %6.4e", x);
+
+// Result
+// The fraction of vacancy sites in metal = 8.4670e-07
|