18 files changed, 311 insertions, 0 deletions
diff --git a/635/CH1/EX1.1/Ch01Ex1.sci b/635/CH1/EX1.1/Ch01Ex1.sci
new file mode 100755
index 000000000..42263bbe7
--- /dev/null
+++ b/635/CH1/EX1.1/Ch01Ex1.sci
@@ -0,0 +1,13 @@
+// Scilab Code Ex1.1 Relationship among cyrstal elements: Page-2 (2010)
+f = 18;    // Number of faces of the quartz crystal
+c = 14;    // Number of angles in the quratz crystal
+// The relationship amongst the crystal elements can be 
+// expressed by the following formula:
+//    f + c = e + 2;
+// Solving for e 
+e = f + c - 2;
+disp (e, "The number of edges of the quartz crystal is : ");
+
+// Result 
+// The number of edges of the quartz crystal is : 
+//    30
+\ No newline at end of file
diff --git a/635/CH1/EX1.17/Ch01Ex17.sci b/635/CH1/EX1.17/Ch01Ex17.sci
new file mode 100755
index 000000000..b8924e70c
--- /dev/null
+++ b/635/CH1/EX1.17/Ch01Ex17.sci
@@ -0,0 +1,13 @@
+// Scilab Code Ex 1.17 Angle between two crystal directions: Page-23 (2010)
+h1 = 1;k1 = 1;l1 = 1; // Miller indices of first set of planes
+h2 = 0;k2 = 0;l2 = 1; // Miller indices of second set of planes
+// We know that
+// cos(theta) = (h1*h2+k1*k2+l1*l2)/(sqrt(h1^2+k1^2+l1^2)*sqrt(h1^2+k1^2+l1^2))
+// Solving for theta
+theta = acos((h1*h2+k1*k2+l1*l2)/(sqrt(h1^2+k1^2+l1^2)*sqrt(h2^2+k2^2+l2^2)));
+printf("\nThe angle between [%d%d%d] and [%d%d%d] directions in the cubic crystal, in degrees, is : %4.2f", h1,k1,l1,h2,k2,l2, theta*180/%pi);
+
+// Result 
+// The angle between [111] and [001] directions in the cubic crystal, in degrees, is :
+//     54.74
+
diff --git a/635/CH1/EX1.18/Ch01Ex18.sci b/635/CH1/EX1.18/Ch01Ex18.sci
new file mode 100755
index 000000000..191fb004e
--- /dev/null
+++ b/635/CH1/EX1.18/Ch01Ex18.sci
@@ -0,0 +1,13 @@
+// Scilab Code Ex 1.18 Angle between two directions of cubic crystal: Page-23(2010)
+h1 = 1; k1 = 1; l1 = 1 // Miller indices for first set of planes
+h2 = -1; k2 = -1; l2 = 1;    // Miller indices for second set of planes
+// We know that
+// cos(theta) = (h1*h2+k1*k2+l1*l2)/(sqrt(h1^2+k1^2+l1^2)*sqrt(h2^2+k2^2+l2^2))
+// Solving for theta
+theta = acos((h1*h2+k1*k2+l1*l2)/(sqrt(h1^2+k1^2+l1^2)*sqrt(h2^2+k2^2+l2^2)));
+printf("\nThe angle between [%d%d%d] and [%d %d %d] directions in the cubic crystal, in degrees, is : %4.1f", h1,k1,l1,h2,k2,l2, theta*180/%pi);
+
+// Result 
+// The angle between [111] and [-1-1 1] directions in the cubic crystal, in degrees, is :
+//     109.5
+
diff --git a/635/CH1/EX1.19/Ch01Ex19.sci b/635/CH1/EX1.19/Ch01Ex19.sci
new file mode 100755
index 000000000..c578f60f8
--- /dev/null
+++ b/635/CH1/EX1.19/Ch01Ex19.sci
@@ -0,0 +1,13 @@
+// Scilab Code Ex 1.19 Miller indices of the crystal plane: Page-25 (2010)
+m = 2; n = 3; p = 6; // 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([m,n,p]))); // Find l.c.m. of m,n and p
+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 2 1)
+\ No newline at end of file
diff --git a/635/CH1/EX1.2/Ch01Ex2.sci b/635/CH1/EX1.2/Ch01Ex2.sci
new file mode 100755
index 000000000..50be14759
--- /dev/null
+++ b/635/CH1/EX1.2/Ch01Ex2.sci
@@ -0,0 +1,14 @@
+// Scilab Code Ex1.2 Primitive unit cell: Page-4 (2010)
+a = 3, b = 3;     // Lattice translation vectors along X and Y direction, angstrom
+c_bar = 3;    // Assumed translation vector along Z direction, angstrom
+c = 1.5*(a+b+c_bar);    // Real translation vector along Z direction, angstrom
+printf("\n%3.1f is the body centered position of a cubic unit cell defined by the primitive translation vectors a, b and c_bar.", c);
+V_con = a^3;    // Volume of conventional unit cell, metre cube
+V_primitive = 1/2*V_con;    // Volume of primitive unit cell, metre cube
+printf("\nThe volume of conventional unit cell: %2d angstrom cube", V_con);
+printf("\nThe volume of primitive unit cell: %4.1f angstrom cube", V_primitive);
+
+// Result 
+// 13.5 is the body centered position of a cubic unit cell defined by the primitive translation vectors a, b and c_bar.
+// The volume of conventional unit cell: 27 angstrom cube
+// The volume of primitive unit cell: 13.5 angstrom cube 
+\ No newline at end of file
diff --git a/635/CH1/EX1.20/Ch01Ex20.sci b/635/CH1/EX1.20/Ch01Ex20.sci
new file mode 100755
index 000000000..33e9f1f9b
--- /dev/null
+++ b/635/CH1/EX1.20/Ch01Ex20.sci
@@ -0,0 +1,16 @@
+// Scilab Code Ex 1.20 Indices of lattice plane: Page-25 (2010)
+m = 10000; // Coefficient of intercept along x-axis, can be taken as some large value
+n = 2;    // Coefficient of intercept along y-axis
+p = 1/2;    // Coefficient of intercept along z-axis
+m_inv = 1/m;    // Reciprocate m
+n_inv = 1/n;    // Reciprocate n
+p_inv = 1/p;    // Reciprocate p
+mul_fact = n;    // multiplicative factor
+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 :
+//     0, 1, 4
diff --git a/635/CH1/EX1.21/Ch01Ex21.sci b/635/CH1/EX1.21/Ch01Ex21.sci
new file mode 100755
index 000000000..4dbe299a1
--- /dev/null
+++ b/635/CH1/EX1.21/Ch01Ex21.sci
@@ -0,0 +1,19 @@
+// Scilab Code Ex 1.21 Length of the intercepts: Page-26 (2010)
+a = 1.21D-10;    // Lattice parameter of the unit cell, m    
+b = 1.84D-10;    // Lattice parameter of the unit cell, m
+c = 1.97D-10;    // Lattice parameter of the unit cell, m
+p = 1/2;    // Reciprocal of miller index on x-axis
+q = 1/3;    // Reciprocal of miller index on y-axis
+r = 1/(-1);    // Reciprocal of miller index on z-axis
+l1 = 1.21D-10;    // Actual length of the intercept along x-axis, m
+mul_fact = l1/(p*a); // Calculate multiplication factor
+l2 = mul_fact*q*b;    // Actual length of the interceptalong y-axis, m
+l3 = mul_fact*r*c;    // Actual length of the intercept along z-axis, m
+disp(l2/1D-10, "The length of the intercept along y-axis, in angstrom, is : ");
+disp(l3/1D-10, "The length of the intercept along z-axis, in angstrom, is : ");
+
+// Result 
+// The length of the intercept along y-axis, in angstrom, is :
+//     1.2266667
+// The length of the intercept along z-axis, in angstrom, is :
+//   - 3.94
diff --git a/635/CH1/EX1.22/Ch01Ex22.sci b/635/CH1/EX1.22/Ch01Ex22.sci
new file mode 100755
index 000000000..31b8ccdae
--- /dev/null
+++ b/635/CH1/EX1.22/Ch01Ex22.sci
@@ -0,0 +1,22 @@
+// Scilab Code Ex 1.22 Miller indices of lattice plane: Page-26 (2010)
+a = 4;    // Lattice parameter of the unit cell    
+b = 3;    // Lattice parameter of the unit cell
+c = 2;    // Lattice parameter of the unit cell
+l1 = 2;    // Length of the intercept along x-axis, m
+l2 = 3;    // Length of the intercept along y-axis, m
+l3 = 4;    // Length of the intercept along z-axis, m
+l = l1/a;    // Intercept per unit translation along x-axis
+m = l2/b;    // Intercept per unit translation along y-axis
+n = l3/c;    // Intercept per unit translation along z-axis
+r1 = 1/l;    // Reciprocal of l
+r2 = 1/m;    // Reciprocal of m
+r3 = 1/n;    // Reciprocal of n
+m1 = 2*r1;    // miller index along x-axis
+m2 = 2*r2;    // miller index along y-axis
+m3 = 2*r3;    // miller index along z-axis
+printf("The required miller indices of the plane are : %d %d %d", m1, m2, m3);
+
+// Result 
+// The required miller indices of the plane are :
+//     4, 2, 1
+
diff --git a/635/CH1/EX1.23/Ch01Ex23.sci b/635/CH1/EX1.23/Ch01Ex23.sci
new file mode 100755
index 000000000..a3a69a500
--- /dev/null
+++ b/635/CH1/EX1.23/Ch01Ex23.sci
@@ -0,0 +1,23 @@
+// Scilab Code Ex 1.23 Indices of tetragonal lattice: Page-26 (2010)
+// For a tetragonal system we have a = b
+a = 1;    // Lattice parameter of the unit cell along x-axis    
+b = 1;    // Lattice parameter of the unit cell along y-axis
+c = 1.5;    // Lattice parameter of the unit cell along z-axis
+l1 = 3;    // Length of the intercept along x-axis, angstrom
+l2 = 4;    // Length of the intercept along y-axis, angstrom
+l3 = 3;    // Length of the intercept along z-axis, angstrom
+l = l1/a;    // Intercept per unit translation along x-axis
+m = l2/b;    // Intercept per unit translation along y-axis
+n = l3/c;    // Intercept per unit translation along z-axis
+r1 = 1/l;    // Reciprocal of l
+r2 = 1/m;    // Reciprocal of m
+r3 = 1/n;    // Reciprocal of n
+mul_fact = double(lcm(int32([l,m,n])));
+m1 = mul_fact*r1;    // miller index along x-axis
+m2 = mul_fact*r2;    // miller index along y-axis
+m3 = mul_fact*r3;    // miller index along z-axis
+printf("The required miller indices of the plane are : %d %d %d", m1, m2, m3);
+
+// Result 
+// The required miller indices of the plane are : 4 3 6
+
diff --git a/635/CH1/EX1.24/Ch01Ex24.sci b/635/CH1/EX1.24/Ch01Ex24.sci
new file mode 100755
index 000000000..cb9f21a91
--- /dev/null
+++ b/635/CH1/EX1.24/Ch01Ex24.sci
@@ -0,0 +1,18 @@
+// Scilab Code Ex 1.24 Miller-Bravias indices for Miller indices: Page-29 (2010)
+function [i] = f(h,k)
+    i = -(h + k);
+endfunction     
+h1 = 1; k1 = 1; l1 = 0 ;  // First set of Miller indices
+h2 = 1; k2 = -1; l2 = 0;    // Second set of miller indices
+h3 = 3; k3 = 4; l3 = 5;    // Third set of miller indices
+h4 = 3; k4 = -4; l4 = 5;    // Fourth set of miller indices
+printf("\nThe Miller-Bravias indices corresponding to the miller indices (%d %d %d), = (%d %d %d %d)", h1, k1, l1, h1, k1, f(h1,k1), l1);
+printf("\nThe Miller-Bravias indices corresponding to the miller indices (%d %d %d), = (%d %d %d %d)", h2, k2, l2, h2, k2, f(h2,k2), l2);
+printf("\nThe Miller-Bravias indices corresponding to the miller indices (%d %d %d), = (%d %d %d %d)", h3, k3, l3, h3, k3, f(h3,k3), l3);
+printf("\nThe Miller-Bravias indices corresponding to the miller indices (%d %d %d), = (%d %d %d %d)", h4, k4, l4, h4, k4, f(h4,k4), l4);
+
+// Result 
+// The Miller-Bravias indices corresponding to the miller indices (1 1 0), = (1 1 -2 0)
+// The Miller-Bravias indices corresponding to the miller indices (1 -1 0), = (1 -1 0 0)
+// The Miller-Bravias indices corresponding to the miller indices (3 4 5), = (3 4 -7 5)
+// The Miller-Bravias indices corresponding to the miller indices (3 -4 5), = (3 -4 1 5)
+\ No newline at end of file
diff --git a/635/CH1/EX1.25/Ch01Ex25.sci b/635/CH1/EX1.25/Ch01Ex25.sci
new file mode 100755
index 000000000..cf7c74355
--- /dev/null
+++ b/635/CH1/EX1.25/Ch01Ex25.sci
@@ -0,0 +1,51 @@
+// Scilab Code Ex 1.25 Miller Bravias indices of lattice planes: Page-30 (2010)
+function [h] = fh(H,K)    // Function for calculating (2H-K)/3
+    h = (2*H - K)/3;    
+endfunction     
+
+function [k] = fk(H,K)    // Function for calculating (2K-H)/3
+    k = (2*K - H)/3;
+endfunction
+
+function [i] = f(h,k)    // Function for calculating i
+    i = -(h + k);
+endfunction
+
+function [l] = fl(L)    // Function for calculating l
+    l = L;
+endfunction
+
+H1 = 1; K1 = 0; L1 = 0 ;  // First set of Miller indices
+H2 = 0; K2 = 1; L2 = 0;    // Second set of miller indices
+H3 = 1; K3 = 1; L3 = 0;    // Third set of miller indices
+
+h1 = fh(H1,K1)*3;    // Call function fh
+k1 = fk(H1,K1)*3;    // Call function fk
+l1 = fl(L1)*3;       // Call function fl
+i1 = f(h1,k1);       // Call function 
+
+h2 = fh(H2,K2)*3;    // Call function fh
+k2 = fk(H2,K2)*3;    // Call function fk
+l2 = fl(L2)*3;        // Call function l2
+i2 = f(h2,k2);        // Call function f
+
+h3 = fh(H3,K3)*3;    // Call function fh
+k3 = fk(H3,K3)*3;    // Call function fk
+l3 = fl(L3)*3;       // Call function l3
+i3 = f(h3,k3);        // Call function f
+
+printf("\n The Miller Bravias indices of [%d%d%d] are [%d %d %d %d]", H1, K1, L1, h1, k1,i1,l1);
+printf("\n The Miller Bravias indices of [%d%d%d] are [%d %d %d %d]", H2, K2, L2, h2, k2,i2,l2);
+printf("\n The Miller Bravias indices of [%d%d%d] are [%d %d %d %d]", H3, K3, L3, h3, k3,i3,l3);
+
+// Result
+// The Miller Bravias indices of [100] are [2 -1 -1 0]
+// The Miller Bravias indices of [010] are [-1 2 -1 0]
+// The Miller Bravias indices of [110] are [1 1 -2 0] 
+
+
+
+
+
+
+
diff --git a/635/CH1/EX1.26/Ch01Ex26.sci b/635/CH1/EX1.26/Ch01Ex26.sci
new file mode 100755
index 000000000..e22c6795b
--- /dev/null
+++ b/635/CH1/EX1.26/Ch01Ex26.sci
@@ -0,0 +1,12 @@
+// Scilab Code Ex 1.26 Lattice parameter of a cubic crystal: Page-33 (2010)
+h = 1; k = 1; l = 1; // Miller Indices for planes in a cubic crystal
+d = 2D-10;    // Interplanar spacing, m
+// For cubic crystals, the interplanar spacing is given by
+// d = a/(h^2+k^2+l^2)^1/2;
+// Solving for a
+a = (h^2+k^2+l^2)^(1/2)*d;   // lattice parameter of cubic crystal, m 
+disp(a/1D-10, "The lattice parameter of the cubic crystal, in angstrom, is :");
+
+// Result
+// The lattice parameter of the cubic crystal, in angstrom, is :
+// 3.4641016
diff --git a/635/CH1/EX1.27/Ch01Ex27.sci b/635/CH1/EX1.27/Ch01Ex27.sci
new file mode 100755
index 000000000..91f6b1047
--- /dev/null
+++ b/635/CH1/EX1.27/Ch01Ex27.sci
@@ -0,0 +1,9 @@
+// Scilab Code Ex 1.27 Interplanar spacing in tetragonal crystal: Page-33 (2010)
+h = 1; k = 0; l = 1; // Miller Indices for planes in a cubic crystal
+a = 2.42D-10;   b = 2.42D-10;   c = 1.74D-10; // Lattice parameters of a tetragonal crystal, each in m
+d = [(h^2+k^2)/a^2 + l^2/c^2]^(-1/2); // The interplanar spacing for cubic crystal, m  
+disp(d/1D-10, "The interplanar spacing between consecutive (101) planes : in angstrom, is :"); 
+
+// Result
+// The interplanar spacing between consecutive (101) planes : in angstrom, is :
+// 1.4127338
+\ No newline at end of file
diff --git a/635/CH1/EX1.28/Ch01Ex28.sci b/635/CH1/EX1.28/Ch01Ex28.sci
new file mode 100755
index 000000000..34fad3969
--- /dev/null
+++ b/635/CH1/EX1.28/Ch01Ex28.sci
@@ -0,0 +1,9 @@
+// Scilab Code Ex 1.28 Interplanar spacing in cubic crystal: Page-36 (2010)
+h = 3; k = 2; l = 1; // Miller Indices for planes in a cubic crystal
+a = 4.21D-10;    // Interatomic spacing, m
+d = a/(h^2+k^2+l^2)^(1/2);  // The interplanar spacing for cubic crystals, m
+disp(d/1D-10, "The interplanar spacing between consecutive (321) planes : in angstrom, is :");
+
+// Result
+// The interplanar spacing between consecutive (321) planes : in angstrom, is :
+// 1.1251698
diff --git a/635/CH1/EX1.3/Ch01Ex3.sci b/635/CH1/EX1.3/Ch01Ex3.sci
new file mode 100755
index 000000000..c6394ab66
--- /dev/null
+++ b/635/CH1/EX1.3/Ch01Ex3.sci
@@ -0,0 +1,14 @@
+// Scilab Code Ex1.3 Number of Lattice points per unit cell Page-9 (2010)
+a = 3.60D-10;    // Lattice parameter, m:
+M = 63.6;        // Atomic weight, gram per mole
+d = 8960D+03;    // Density of copper, g per metre cube
+N = 6.023D+23;    // Avogadro's No.
+// Volume of the unit cell is given by
+// a^3 = M*n/(N*d)
+// Solving for n
+n = a^3*d*N/M; // Number of lattice points per unit cell
+disp (n, "The number of atoms per unit cell for an fcc lattice of copper crystal is :");
+
+// Result 
+// The number of atoms per unit cell for an fcc lattice of copper crystal
+// 3.9588702
+\ No newline at end of file
diff --git a/635/CH1/EX1.4/Ch01Ex4.sci b/635/CH1/EX1.4/Ch01Ex4.sci
new file mode 100755
index 000000000..0fe424fc6
--- /dev/null
+++ b/635/CH1/EX1.4/Ch01Ex4.sci
@@ -0,0 +1,14 @@
+// Scilab Code Ex 1.4 Lattice constant of a unit cell: Page-9 (2010)
+M = 58.5;        // Atomic weight of NaCl, gram per mole
+d = 2180D+03;    // Density of rock salt, per metre cube
+n = 4;    // No. of atoms per unit cell for an fcc lattice of NaCl crystal
+N = 6.023D+23;    // Avogadro's No.
+// Volume of the unit cell is given by
+// a^3 = M*n/(N*d)
+// Solving for a
+a = (n*M/(d*N))^(1/3);    // Lattice constant of unit cell of NaCl
+disp (a/1D-10, "Lattice constant for the rock salt (NaCl) crystal, in angstrom, is : ");
+
+// Result 
+// Lattice constant for the rock salt (NaCl) crystal, in angstrom, is :
+// 5.6275
+\ No newline at end of file
diff --git a/635/CH1/EX1.5/Ch01Ex5.sci b/635/CH1/EX1.5/Ch01Ex5.sci
new file mode 100755
index 000000000..d58bd7e43
--- /dev/null
+++ b/635/CH1/EX1.5/Ch01Ex5.sci
@@ -0,0 +1,17 @@
+// Scilab Code Ex 1.5 Density of diamond: Page-9 (2010)
+a = 3.57D-10; // Lattice parameter of a diamond crystal
+M = 12D-03;     // Atomic weight of diamond, kg per mole
+n1 = 8;    // No. of corner atoms in the diamond cubic unit cell
+n2 = 6;    // No. of face centered atoms in the diamond cubic unit cell
+n3 = 4;    // No. of atoms completely within the unit cell
+n = 1/8*n1+1/2*n2+1*n3;    // No. of atoms per unit cell for an fcc lattice of NaCl crystal
+N = 6.023D+23;    // Avogadro's No.
+// Volume of the unit cell is given by
+// a^3 = M*n/(N*d)
+// Solving for d
+d = M*n/(N*a^3);    // Density of diamond cubic unit cell
+disp (round(d), "Density of diamond cubic unit cell, in kg per metre cube, is : ");
+
+// Result 
+// Density of diamond cubic unit cell, in kg per metre cube, is :
+// 3503
+\ No newline at end of file
diff --git a/635/CH1/EX1.6/Ch01Ex6.sci b/635/CH1/EX1.6/Ch01Ex6.sci
new file mode 100755
index 000000000..de265e9cd
--- /dev/null
+++ b/635/CH1/EX1.6/Ch01Ex6.sci
@@ -0,0 +1,21 @@
+// Scilab Code Ex 1.6 Calculating Unit cell dimensions: Page-9 (2010)
+d = 2.7D+03;    // Density of fcc structure of aluminium, kg per metre cube
+M = 26.98D-03;     // Atomic weight of aluminium, kg per mole
+n = 4;       // No. of atoms per unit cell of fcc lattice structure of  aluminium
+N = 6.023D+23;    // Avogadro's No.
+// Volume of the unit cell is given by
+// a^3 = M*n/(N*d)
+// Solving for a
+a = ((M*n)/(N*d))^(1/3);    // Lattice parameter of alumitnium unit cell
+// For an fcc cryatal lattice, 
+// 2^(1/2) = 4R = 2D
+// Solving for D
+D = (a/2^(1/2)); // Diameter of aluminium atom
+disp (a/1D-10, "The Lattice parameter of aluminium, in angstrom, is : ");
+disp(D/1D-10, "The diameter of aluminium atom, in angstrom, is : ");
+
+// Result 
+// The Lattice parameter of aluminium, in angstrom, is :
+//     4.0486332
+// The diameter of aluminium atom, in angstrom, is :
+//     2.862816  
+\ No newline at end of file