diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /2309/CH5 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '2309/CH5')
32 files changed, 674 insertions, 0 deletions
diff --git a/2309/CH5/EX5.1/Ex5_1.sce b/2309/CH5/EX5.1/Ex5_1.sce new file mode 100755 index 000000000..bf59b4231 --- /dev/null +++ b/2309/CH5/EX5.1/Ex5_1.sce @@ -0,0 +1,21 @@ +// Chapter 5 Example 1
+//==============================================================================
+clc;
+clear;
+
+//input data
+//Copper has FCC structure
+
+r = 1.273; // Atomic radius in angstrom
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 63.5; // Atomic weight of copper in grams
+n = 4; // No. of atoms per unit cell for FCC
+
+//Calculations
+r1 = r*10^-10; // Radius conversion from angstrom to m
+a = (4*r1)/sqrt(2); // lattice parameter for FCC
+p = (n*A)/(N*a^3); // Density of copper
+
+//Output
+
+mprintf('Lattice Constant a = %3.1e m\n Density of copper = %3.1f kg/m^3',a,p);
diff --git a/2309/CH5/EX5.10/Ex5_10.sce b/2309/CH5/EX5.10/Ex5_10.sce new file mode 100755 index 000000000..9846717dd --- /dev/null +++ b/2309/CH5/EX5.10/Ex5_10.sce @@ -0,0 +1,21 @@ +// Chapter 5 Example 10
+//==============================================================================
+clc;
+clear;
+
+// input data
+// FCC structured crystal
+
+p = 6250; // Density of crystal in kg/m^3
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 60.2; // molecular weight
+n = 4; // No. of atoms per unit cell for FCC
+
+//Calculations
+
+a = ((n*A)/(N*p))^(1/3);
+
+//Output
+
+mprintf('Lattice Constant a = %3.1e m ',a);
+//==============================================================================
diff --git a/2309/CH5/EX5.11/Ex5_11.sce b/2309/CH5/EX5.11/Ex5_11.sce new file mode 100755 index 000000000..71b35a3bd --- /dev/null +++ b/2309/CH5/EX5.11/Ex5_11.sce @@ -0,0 +1,19 @@ +// Chapter 5 Example 11
+//==============================================================================
+clc;
+clear;
+
+//input data
+// (321) plane in simple cubic lattice
+h = 3; // miller indice
+k = 2; // miller indice
+l = 1; // miller indice
+a = 4.12 // inter atomic space Å
+
+// Calculations
+dhkl = a/sqrt((h^2)+(k^2)+(l^2)); // interplanar distance
+
+// Output
+mprintf('d = %3.2f Å',dhkl);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.12/Ex5_12.sce b/2309/CH5/EX5.12/Ex5_12.sce new file mode 100755 index 000000000..b10562ac6 --- /dev/null +++ b/2309/CH5/EX5.12/Ex5_12.sce @@ -0,0 +1,22 @@ +// Chapter 5 Example 12
+//==============================================================================
+clc;
+clear;
+
+// input data
+// BCC structured crystal
+
+p = 7860; // Density of iron in kg/m^3
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 55.85; // Atomic weight
+n = 2; // No. of atoms per unit cell for BCC
+
+//Calculations
+
+a = ((n*A)/(N*p))^(1/3); //lattice constant
+
+//Output
+
+mprintf('Lattice Constant of Fe = %3.3f Å \n',a*10^10);
+mprintf(' Note: density of iron is taken as 7.86 instead of 7860 in calculation')
+//==============================================================================
diff --git a/2309/CH5/EX5.14/Ex5_14.sce b/2309/CH5/EX5.14/Ex5_14.sce new file mode 100755 index 000000000..4395408db --- /dev/null +++ b/2309/CH5/EX5.14/Ex5_14.sce @@ -0,0 +1,15 @@ +// Chapter 5 Example 14
+//==============================================================================
+clc;
+clear;
+
+// input data
+r = 0.123*10^-10; // Radius of the atom
+
+// Calculations
+a = (4*r)/sqrt(3); // Lattice constant in m For a BCC structure
+V = a*a*a; // Volume of BCC
+
+// Output
+mprintf('Volume of the unit cell = %3.4e m^3',V);
+//==============================================================================
diff --git a/2309/CH5/EX5.15/Ex5_15.sce b/2309/CH5/EX5.15/Ex5_15.sce new file mode 100755 index 000000000..70ae2432f --- /dev/null +++ b/2309/CH5/EX5.15/Ex5_15.sce @@ -0,0 +1,24 @@ +// Chapter 5 Example 15
+//==============================================================================
+clc;
+clear;
+
+// input data
+a = 0.05; // unit cell edge of an orthorhombic crystal in nm
+b = 0.05; // unit cell edge of an orthorhombic crystal in nm
+c = 0.03; // unit cell edge of an orthorhombic crystal in nm
+Ia = 0.025 // intercept on 'a' in nm
+Ib = 0.02 // intercept on 'b' in nm
+Ic = 0.01 // intercept on 'c' in nm
+
+//Calculations
+
+h = a/Ia; // miller indice h
+k = b/Ib; // miller indice k
+l = c/Ic // miller indice l
+
+// Output
+mprintf('Miller indices (h k l) = (%d %d %d)',h,k,l);
+//==============================================================================
+
+
diff --git a/2309/CH5/EX5.16/Ex5_16.sce b/2309/CH5/EX5.16/Ex5_16.sce new file mode 100755 index 000000000..3645662ba --- /dev/null +++ b/2309/CH5/EX5.16/Ex5_16.sce @@ -0,0 +1,20 @@ +// Chapter 5 Example 16
+//==============================================================================
+clc;
+clear;
+// Magnesium has HCP structure
+// for HCF(Hexagonal closed packed structure) consider the relation between 'c' and 'a';
+// c/a = sqrt(8/3) = 1.6329
+//input data
+r = 0.1605*10^-9; // radius of magnesium atom in m
+
+// Calculations
+
+a = 2*r // lattice constant of HCP
+c = a*sqrt(8/3); // relation b/w c and a in HCP
+V = (3*3^0.5)*(a*a*c)/2; //Volume of unit cell in m^3
+
+// Output
+mprintf('Volume of the unit cell of magnesium = %3.3e m^3',V);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.17/Ex5_17.sce b/2309/CH5/EX5.17/Ex5_17.sce new file mode 100755 index 000000000..ff8cdc1d2 --- /dev/null +++ b/2309/CH5/EX5.17/Ex5_17.sce @@ -0,0 +1,24 @@ +// Chapter 5 Example 17
+//==============================================================================
+clc;
+clear;
+
+//input data
+// (101),(221) planes in simple cubic lattice
+h1 = 1; // miller indice
+k0 = 0; // miller indice
+l1 = 1; // miller indice
+h2 = 2; // miller indice
+k2 = 2; // miller indice
+l1 = 1; // miller indice
+a = 4.2 // inter atomic space Å
+
+// Calculations
+d101 = a/sqrt((h1^2)+(k0^2)+(l1^2)); // interplanar distance
+d221 = a/sqrt((h2^2)+(k2^2)+(l1^2)); // interplanar distance
+
+
+// Output
+mprintf('d(101) = %3.4f Å\n d(221) = %3.1f Å ',d101,d221);
+//=============================================================================
+
diff --git a/2309/CH5/EX5.2/Ex5_2.sce b/2309/CH5/EX5.2/Ex5_2.sce new file mode 100755 index 000000000..6c19c2acd --- /dev/null +++ b/2309/CH5/EX5.2/Ex5_2.sce @@ -0,0 +1,26 @@ +// Chapter 5 Example 1
+//==============================================================================
+clc;
+clear;
+
+//input data
+//given intercepts 3,4 and ∞, the recipocals of intercepts is
+// (1/3):(1/4):(1/∞)
+// LCM = 12
+// multiplying by LCM we get miller indices
+// miller indices of a plane are the smallest integers of the reciprocals of its intercerpts
+// therefore miller indices(h k l) is (4 3 0);
+
+h = 4; // miller indice
+k = 3; // miller indice
+l = 0; // miller indice
+a = 2; // primitive vector of lattice in angstrom
+
+//Calculations
+
+dhkl = a/sqrt((h^2)+(k^2)+(l^2)); // interplanar distance
+
+//Output
+mprintf('Miller indices = (4 3 0)\n');
+mprintf(' The interplanar distance d = %3.1f Å',dhkl);
+//==============================================================================
diff --git a/2309/CH5/EX5.3/Ex5_3.sce b/2309/CH5/EX5.3/Ex5_3.sce new file mode 100755 index 000000000..881fa4018 --- /dev/null +++ b/2309/CH5/EX5.3/Ex5_3.sce @@ -0,0 +1,26 @@ +// Chapter 5 Example 3
+//==============================================================================
+clc;
+clear;
+
+//input data
+//α-Iron solidifies to BCC structure
+
+r = 1.273; // Atomic radius in angstrom
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 55.85; // Atomic weight of α-Iron in kilograms
+n = 2; // No. of atoms per unit cell for BCC
+p = 7860; // density in kg/m^-3
+
+//Calculations
+
+// p = (n*A)/(N*a^3); density
+
+a = ((n*A)/(N*p))^(1/3); // lattice constant
+a1 = a*10^10; // m to angstrom conversion
+r = (a1*sqrt(3))/4 // atomic radius for BCC
+
+//Output
+mprintf('The Radius of the atom = %3.5f Å\n',r);
+mprintf(' Note : atomic wt taken as 55.58*10^-3 instead of 55.85 in calculation')
+//==============================================================================
diff --git a/2309/CH5/EX5.4/Ex5_4.sce b/2309/CH5/EX5.4/Ex5_4.sce new file mode 100755 index 000000000..c5e8434f5 --- /dev/null +++ b/2309/CH5/EX5.4/Ex5_4.sce @@ -0,0 +1,24 @@ +// Chapter 5 Example 4
+//==============================================================================
+clc;
+clear;
+
+//input data
+lamda = 1.5418; // wavelength in Å
+h = 1; // miller indice
+k = 1; // miller indice
+l = 1; // miller indice
+n = 1; // given first order
+theta = 30; // diffraction angle in degrees
+
+// Calculations
+theta1 = theta*%pi/180; // degree to radian conversion
+// d = (n*lamda)/(2*sinθ); by Braggs law ------------- 1
+// d = a/sqrt((h^2)+(k^2)+(l^2)); interplanar distance ------------ 2
+// equating 1 and 2
+
+a = (n*lamda*sqrt((h^2)+(k^2)+(l^2))/(2*sin(theta1)))
+
+// Output
+mprintf('Interatomic spacing a = %f Å',a);
+//==============================================================================
diff --git a/2309/CH5/EX5.5/Ex5_5.sce b/2309/CH5/EX5.5/Ex5_5.sce new file mode 100755 index 000000000..5ed750009 --- /dev/null +++ b/2309/CH5/EX5.5/Ex5_5.sce @@ -0,0 +1,28 @@ +// Chapter 5 Example 5
+//==============================================================================
+clc;
+clear;
+
+//input data
+h1 = 1; // miller indice
+k1 = 1; // miller indice
+l1 = 1; // miller indice
+h0 = 0; // miller indice
+k0 = 0; // miller indice
+l0 = 0; // miller indice
+
+// calculations
+// dhkl = a/sqrt((h^2)+(k^2)+(l^2)); // interplanar distance
+// assume a = 1(constant) for easier calculation in scilab
+
+a = 1;
+d100 = a/sqrt((h1^2)+(k0^2)+(l0^2)); // interplanar distance
+d110 = a/sqrt((h1^2)+(k1^2)+(l0^2)); // interplanar distance
+d111 = a/sqrt((h1^2)+(k1^2)+(l1^2)); // interplanar distance
+
+// Output
+mprintf('d100 : d110 : d111 = %d : %3.2f : %3.2f',d100,d110,d111);
+
+//==============================================================================
+
+
diff --git a/2309/CH5/EX5.6/Ex5_6.sce b/2309/CH5/EX5.6/Ex5_6.sce new file mode 100755 index 000000000..26c15f557 --- /dev/null +++ b/2309/CH5/EX5.6/Ex5_6.sce @@ -0,0 +1,19 @@ +// Chapter 5 Example 6
+//==============================================================================
+clc;
+clear;
+
+// input data
+// Aluminium is FCC
+a = 0.405*10^-9; // lattice constant of aluminium
+t = 0.005*10^-2; // thickness of aluminium foil in m
+s = 25*10^-2; // side of square in m
+
+//Calculations
+VUC = a^3; // volume of unit cell
+Val = (s^2)*t // volume of aluminium foil (area*thickness)
+N = Val/VUC // Number if unit cells
+
+//Output
+mprintf('Number of unit cells = %3.3e',N);
+//==============================================================================
diff --git a/2309/CH5/EX5.7/Ex5_7.sce b/2309/CH5/EX5.7/Ex5_7.sce new file mode 100755 index 000000000..7abef37bc --- /dev/null +++ b/2309/CH5/EX5.7/Ex5_7.sce @@ -0,0 +1,21 @@ +// Chapter 5 Example 7
+//==============================================================================
+clc;
+clear;
+
+// input data
+// metallic iron changes from BCC to FCC form at 910 degress
+rb = 0.1258*10^-9; // atomic radius of BCC iron atom
+rf = 0.1292*10^-9; // atomic radius of FCC iron atom
+
+// Calculations
+
+ab = (4*rb)/(sqrt(3)); // lattice constant for BCC
+Vbcc = (ab^3)/2; // volume occupied by one BCC atom
+af = (4*rf)/(sqrt(2)) // lattice constant for FCC
+Vfcc = (af^3)/4; // volume occupied by one FCC atom
+dv = ((Vbcc-Vfcc)/Vbcc)*100 // percentage change in volume
+
+// output
+mprintf('During the structural change the percentage change in volume = %3.4f',dv);
+//==============================================================================
diff --git a/2309/CH5/EX5.8/Ex5_8.sce b/2309/CH5/EX5.8/Ex5_8.sce new file mode 100755 index 000000000..98a450dca --- /dev/null +++ b/2309/CH5/EX5.8/Ex5_8.sce @@ -0,0 +1,22 @@ +// Chapter 5 Example 8
+//==============================================================================
+clc;
+clear;
+
+//input data
+//Copper Crystallines in FCC structure
+
+p = 8960; // Density of copper in kg/m^3
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 63.5; // Atomic weight of copper in kg/mol
+n = 4; // No. of atoms per unit cell for FCC
+
+//Calculations
+
+a = ((n*A)/(N*p))^(1/3);
+
+//Output
+
+mprintf('Lattice Constant a = %3.4f Å\n',a*10^10);
+mprintf(' atomic wt of copper is taken as 63.5*10^-3 instead of 63.5 in textbook')
+//==============================================================================
diff --git a/2309/CH5/EX5.9/Ex5_9.sce b/2309/CH5/EX5.9/Ex5_9.sce new file mode 100755 index 000000000..c9da54ece --- /dev/null +++ b/2309/CH5/EX5.9/Ex5_9.sce @@ -0,0 +1,19 @@ +// Chapter 5 Example 9
+//==============================================================================
+clc;
+clear;
+
+//input data
+// (100) planes in rock crystal
+h = 1; // miller indice
+k = 0; // miller indice
+l = 0; // miller indice
+a = 2.814 // lattice constant in Å
+
+// Calculations
+dhkl = a/sqrt((h^2)+(k^2)+(l^2)); // interplanar distance
+
+// Output
+mprintf('d-spacing for (100) plane in rock salt = %3.3f Å',dhkl);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.a.1/A_Ex5_1.sce b/2309/CH5/EX5.a.1/A_Ex5_1.sce new file mode 100755 index 000000000..dafec357c --- /dev/null +++ b/2309/CH5/EX5.a.1/A_Ex5_1.sce @@ -0,0 +1,17 @@ +// Chapter 5 additional Example 1
+//==============================================================================
+clc;
+clear;
+
+// input data
+// Copper has FCC structure
+a = 3.6; // lattice parameter of copper in Å
+
+// Calculations
+
+r = a*sqrt(2)/4; // atomic radius of copper
+
+// Output
+mprintf('Atomic Radius of copper = %3.3f Å',r);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.a.11/A_Ex5_11.sce b/2309/CH5/EX5.a.11/A_Ex5_11.sce new file mode 100755 index 000000000..5356c7731 --- /dev/null +++ b/2309/CH5/EX5.a.11/A_Ex5_11.sce @@ -0,0 +1,19 @@ +// Chapter 5 additional Example 11
+//==============================================================================
+clc;
+clear;
+
+//input data
+// (311) plane in simple cubic lattice
+h = 3; // miller indice
+k = 1; // miller indice
+l = 1; // miller indice
+a = 2.109*10^-10 // lattice constant in m
+
+// Calculations
+dhkl = a/sqrt((h^2)+(k^2)+(l^2)); // interplanar distance
+
+// Output
+mprintf('d = %3.3e m',dhkl);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.a.12/A_Ex5_12.sce b/2309/CH5/EX5.a.12/A_Ex5_12.sce new file mode 100755 index 000000000..882807a08 --- /dev/null +++ b/2309/CH5/EX5.a.12/A_Ex5_12.sce @@ -0,0 +1,19 @@ +// Chapter 5 additional Example 12
+//==============================================================================
+clc;
+clear;
+
+//input data
+
+h = 1; // miller indice
+k = 1; // miller indice
+l = 0; // miller indice
+d = 2.86*10^-10 // interplanar distance in m
+
+// Calculations
+a = d*sqrt((h^2)+(k^2)+(l^2)); // interplanar distance
+
+// Output
+mprintf('Lattice constant a = %3.3e m',a);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.a.13/A_Ex5_13.sce b/2309/CH5/EX5.a.13/A_Ex5_13.sce new file mode 100755 index 000000000..04fe09817 --- /dev/null +++ b/2309/CH5/EX5.a.13/A_Ex5_13.sce @@ -0,0 +1,21 @@ +// Chapter 5 Additional Example 13
+//==============================================================================
+clc;
+clear;
+
+h1 = 1;
+h0 = 0;
+k0 = 0;
+l0 = 0;
+l1 = 1;
+// calculations
+
+// we know that dhkl = a/sqrt( h^2 + k^2 + l^2)
+// let sqrt( h^2 + k^2 + l^2) = p
+p101 = sqrt( h1^2 + k0^2 + l1^2);
+p100 = sqrt( h1^2 + k0^2 + l0^2);
+p001 = sqrt( h0^2 + k0^2 + l1^2);
+
+// output
+mprintf('d101 : d100 : d001 :: a/%3.4f : a/%d : a/%d ',p101,p100,p001);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.14/A_Ex5_14.sce b/2309/CH5/EX5.a.14/A_Ex5_14.sce new file mode 100755 index 000000000..0917dceef --- /dev/null +++ b/2309/CH5/EX5.a.14/A_Ex5_14.sce @@ -0,0 +1,14 @@ +// Chapter 5 additional Example 14
+//==============================================================================
+clc;
+clear;
+
+// if a plane cut intercepts of lengths l1,l2,l3 the on three crystal axes ,then
+// l1 : l2 : l3 = pa : pq :rc
+// where a,b and c are primitive vectors of the unit cell and p,q and r are numbers related to miller indices (hkl) of plane by relation
+// 1/p : 1/q : 1/r = h : k : l
+//since, the crystal is simple cubic a = b = c and given that h = 1, k = 1 and l = 1
+// p : q : r = 1/h : 1/k : 1/l = 1/1 : 1/1 : 1/1
+// p : q : r = 1 : 1 : 1
+//similarly l1 : l2 : l3 = 1a : 1a : 1a
+mprintf('ratio of intercepts on the three axes by (111) plane is l1 : l2 : l3 = 1 : 1 : 1');
diff --git a/2309/CH5/EX5.a.15/A_Ex5_15.sce b/2309/CH5/EX5.a.15/A_Ex5_15.sce new file mode 100755 index 000000000..369fd9d14 --- /dev/null +++ b/2309/CH5/EX5.a.15/A_Ex5_15.sce @@ -0,0 +1,24 @@ +// Chapter 5 additional Example 15
+//==============================================================================
+clc;
+clear;
+
+//input data
+r = 1.246*10^-10; // atomic radius in m
+h1 = 1 // miller indice
+h2 = 2 // miller indice
+k0 = 0 // miller indice
+k1 = 1 // miller indice
+k2 = 2 // miller indice
+l0 = 0 // miller indice
+l1 = 1 // miller indice
+
+// Calculations
+a = (4*r)/sqrt(2); // lattice constant
+d111 = a/sqrt((h1^2)+(k1^2)+(l1^2)); // interplanar distance
+d200 = a/sqrt((h2^2)+(k0^2)+(l0^2)); // interplanar distance
+d220 = a/sqrt((h2^2)+(k2^2)+(l0^2)); // interplanar distance
+
+// Output
+mprintf('d111 = %3.3e m\n d200 = %3.4e m\n d220 = %3.3e m\n',d111,d200,d220');
+//==============================================================================
diff --git a/2309/CH5/EX5.a.16/A_Ex5_16.sce b/2309/CH5/EX5.a.16/A_Ex5_16.sce new file mode 100755 index 000000000..3e594a689 --- /dev/null +++ b/2309/CH5/EX5.a.16/A_Ex5_16.sce @@ -0,0 +1,27 @@ +// Chapter 5 additional Example 16
+//==============================================================================
+clc;
+clear;
+
+//input data
+// the intercept along X-axis be c1 = a
+// the intercept along Y-axis be c2 = b/2 and
+// the intercept along Z-axis be c3 = 3c
+// Therefore, p = c1/a = a/a = 1
+// q = c2/b = (b/2)/b = 1/2
+// r = c3/c = (3c)/c = 3
+// therefore h = 1/p = 1
+// k = 1/q = 2
+// l = 1/r = 1/3
+// lcm of 1 1 and 3 = 3
+h = 1
+k = 2
+l = 1/3
+p = [1 1 3]
+s = lcm(p);
+h1= s*h
+k1= s*k
+l1= s*l;
+// Output
+mprintf('(h k l) = (%d %d %d)',h1,k1,l1);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.17/A_Ex5_17.sce b/2309/CH5/EX5.a.17/A_Ex5_17.sce new file mode 100755 index 000000000..d9c917736 --- /dev/null +++ b/2309/CH5/EX5.a.17/A_Ex5_17.sce @@ -0,0 +1,19 @@ +// Chapter 5 Additional Example 17
+//==============================================================================
+clc;
+clear;
+
+//input data
+
+d = 1.3*10^-10 // interplanar distance
+n = 1; // given first order
+theta = 23; // Bragg reflection angle in degrees
+
+// Calculations
+theta1 = theta*%pi/180; // degree to radian conversion
+// d = (n*lamda)/(2*sinθ); by Braggs law ------------- 1
+lamda = (2*d*sin(theta1)/n)
+
+// Output
+mprintf('Wavelength of X-ray = %3.4f Å',lamda*10^10);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.2/A_Ex5_2.sce b/2309/CH5/EX5.a.2/A_Ex5_2.sce new file mode 100755 index 000000000..8842b1458 --- /dev/null +++ b/2309/CH5/EX5.a.2/A_Ex5_2.sce @@ -0,0 +1,22 @@ +// Chapter 5 additional Example 2
+//==============================================================================
+clc;
+clear;
+
+// input data
+// Copper has FCC structure
+
+r = 1.278; // Atomic radius in angstrom
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 63.54; // Atomic weight of copper
+n = 4; // No. of atoms per unit cell for FCC
+
+//Calculations
+r1 = r*10^-10; // Radius conversion from angstrom to m
+a = (4*r1)/sqrt(2); // lattice parameter for FCC
+p = (n*A)/(N*a^3); // Density of copper
+
+//Output
+
+mprintf(' Density of copper = %3.2f kg/m^3',p);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.3/A_Ex5_3.sce b/2309/CH5/EX5.a.3/A_Ex5_3.sce new file mode 100755 index 000000000..a3aebbff0 --- /dev/null +++ b/2309/CH5/EX5.a.3/A_Ex5_3.sce @@ -0,0 +1,23 @@ +// Chapter 5 additional Example 3
+//==============================================================================
+clc;
+clear;
+
+// input data
+// NaCl has FCC structure
+
+ANa = 23; // atomic wt of sodiim
+ACl = 35.45 // atomic wt of chlorine
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+n = 4 // No. of atoms per unit cell for FCC
+p = 2180; // density in kg/m^-3
+
+// Calculations
+
+// p = (n*A)/(N*a^3); density
+A = ANa+ACl; // atomic wt of NaCl
+a = ((n*A)/(N*p))^(1/3); // lattice constant
+r = a/2 // Distance b/w two adjacent atoms
+//Output
+mprintf('Distance between two adjacent atoms is r = %3.2e m',r);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.4/A_Ex5_4.sce b/2309/CH5/EX5.a.4/A_Ex5_4.sce new file mode 100755 index 000000000..c457edfcb --- /dev/null +++ b/2309/CH5/EX5.a.4/A_Ex5_4.sce @@ -0,0 +1,25 @@ +// Chapter 5 additional Example 4
+//==============================================================================
+clc;
+clear;
+
+// input data
+// iron has BCC structure
+
+r = 1.273; // Atomic radius in angstrom
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 55.85 ; // Atomic weight of Fe
+n = 2; // No. of atoms per unit cell for BCC
+p = 7860; // density in kg/m^-3
+
+//Calculations
+
+// p = (n*A)/(N*a^3); density
+
+a = ((n*A)/(N*p))^(1/3); // lattice constant
+a1 = a*10^10; // m to angstrom conversion
+r = (a1*sqrt(3))/4 // atomic radius for BCC
+
+//Output
+mprintf('The Radius of the Fe = %3.3f Å',r);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.5/A_Ex5_5.sce b/2309/CH5/EX5.a.5/A_Ex5_5.sce new file mode 100755 index 000000000..53e645382 --- /dev/null +++ b/2309/CH5/EX5.a.5/A_Ex5_5.sce @@ -0,0 +1,23 @@ +// Chapter 5 additional Example 5
+//==============================================================================
+clc;
+clear;
+
+// input data
+// KBr has FCC structure
+
+N = 6.023*10^26; // Avagadros number in atoms/kilomole
+A = 119; // Atomic weight of pottasium bromide
+n = 4; // No. of atoms per unit cell for FCC
+p = 2700; // density in kg/m^-3
+
+//Calculations
+
+// p = (n*A)/(N*a^3); density
+
+a = ((n*A)/(N*p))^(1/3); // lattice constant
+a1 = a*10^10; // m to angstrom conversion
+
+//Output
+mprintf('Lattice constant = %3.1f Å',a1);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.6/A_Ex5_6.sce b/2309/CH5/EX5.a.6/A_Ex5_6.sce new file mode 100755 index 000000000..1092c4743 --- /dev/null +++ b/2309/CH5/EX5.a.6/A_Ex5_6.sce @@ -0,0 +1,18 @@ +// Chapter 5 additional Example 6
+//==============================================================================
+clc;
+clear;
+// input data
+a = 4.3*10^-10; // Lattice constant in Å
+p = 960; // Density of crystal in kg/m^3
+A = 23; // Atomic wt
+N = 6.023*10^26; // avogadros no in atoms/kilomole
+
+//Calculations
+
+n = (p*N*(a^3))/A; // No. of atoms per unit cell
+
+// Output
+mprintf('No. of atoms per unit cell = %3.0f (BCC)',n);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.a.7/A_Ex5_7.sce b/2309/CH5/EX5.a.7/A_Ex5_7.sce new file mode 100755 index 000000000..7492185b5 --- /dev/null +++ b/2309/CH5/EX5.a.7/A_Ex5_7.sce @@ -0,0 +1,16 @@ +// Chapter 5 additional Example 7
+//==============================================================================
+clc;
+clear;
+// input data
+// given crystal has BCC structure
+r = 1.2*10^-10; // atomic radius in m
+
+// Calculations
+
+a = (4*r)/sqrt(3); // lattice constant
+V = a^3; // volume of cell
+
+//Output
+mprintf('Volume of the cell = %3.3e m^3',V);
+//==============================================================================
diff --git a/2309/CH5/EX5.a.8/A_Ex5_8.sce b/2309/CH5/EX5.a.8/A_Ex5_8.sce new file mode 100755 index 000000000..c7444c683 --- /dev/null +++ b/2309/CH5/EX5.a.8/A_Ex5_8.sce @@ -0,0 +1,21 @@ +// Chapter 5 additional Example 8
+//==============================================================================
+clc;
+clear;
+// input data
+a = 4*10^-10; // lattice constant of the crystal
+h = 1 // miller indice
+k = 0 // miller indice
+l = 0 // miller indice
+
+//Calculations
+
+// in fig consider (100) plane. the no of atoms in plane ABCD
+N = 4*(1/4); // Number of atoms
+p = N/(a*a); // planar atomic density in atoms/m^2
+p1 = p*10^-6 // planar atomic density in atoms/mm^2
+
+//Output
+mprintf('planar atomic density = %3.2e atoms/mm^2',p1);
+//==============================================================================
+
diff --git a/2309/CH5/EX5.a.9/A_Ex5_9.sce b/2309/CH5/EX5.a.9/A_Ex5_9.sce new file mode 100755 index 000000000..b49c9923b --- /dev/null +++ b/2309/CH5/EX5.a.9/A_Ex5_9.sce @@ -0,0 +1,15 @@ +// Chapter 5 additional Example 9
+//==============================================================================
+clc;
+clear;
+// input data
+// in fig 5(b) the given plane is parallel to X and Z axes.Thus,its numerical intercepts on these axes is infinity
+//The numerical intercept on y axis is 1/2. Thus the numerical intercepts of plane is (∞ 1/2 ∞)
+mprintf('Miller indices of plane shown in fig 5.(b) = (0 2 0)\n');
+// in fig 5(c) the given plane is parallel to Z axis.Thus its numerical intercept on z axis is infinity
+// The numerical intercept on x axis is 1 and y axis is 1/2. this numerical intercepts on plane is (1 1/2 ∞ )
+mprintf(' Miller indices of plane shown in fig 5.(c) = (1 2 0)\n')
+// in fig 5(d) the given plane is parallel to Z axis.Thus its numerical intercept on z axis is infinity
+// The numerical intercept on x axis is 1/2 and y axis is 1/2. this numerical intercepts on plane is (1/2 1/2 ∞ )
+mprintf(' Miller indices of plane shown in fig 5.(d) = (2 2 0)\n')
+//==============================================================================
|