summaryrefslogtreecommitdiff
path: root/635/CH8/EX8.12/Ch08Ex12.sci
blob: 783ae3008b715eb935f03c4d4d3261ed91e22325 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
// Scilab code Ex8.12: Reciprocal lattice parameters from 2-D direct lattice parameters Page 277 (2010)
a = 3e-010;    // First lattice parameter of direct lattice
b = 5e-010;    // Second lattice parameter of direct lattice
theta = 60;    // Angle between two lattice vectors of the direct lattice
// if a_prime and b_prime are the lattice vectors for the reciprocal lattice, then
// a_prime*a = 2*%pi and a_prime*b = 0 
// Similarly, b_prime*b = 2*%pi and b_prime*a = 0
// Solving for a_prime and b_prime, we have
a_prime = 2*%pi/(a*cosd(90-theta)); // Lattice vector for reciprocal lattice, per metre
b_prime = 2*%pi/(b*cosd(90-theta)); // Lattice vector for reciprocal lattice, per metre
printf("\nThe reciprocal lattice vectors are:\n a_prime = %5.2f per angstrom and b_prime = %5.2f per angstrom", a_prime*1e-010, b_prime*1e-010);
// Result
// The reciprocal lattice vectors are:
// a_prime =  2.42 per angstrom and b_prime =  1.45 per angstrom