forked from maherharb/Diffraction
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Main_Form_Powder_Diffraction.m
47 lines (34 loc) · 1.33 KB
/
Main_Form_Powder_Diffraction.m
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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
% EXAMPLE SCRIPT TO GENERATE DIFFRACTION PATTERN FROM NONCOPLANAR DIFFRACTION GEOMETRY
% Last Updated 5/26/17 - Jacob Bolduc
% CLEAN WORKSPACE
home
clear
close all
% PATH TO DIFFWIZ LIBRARY
addpath(genpath('C:\Users\maher\Google drive\DIFFWIZ\'))
% DEFINE LATTICE OR LOAD LATTICE FROM EXISTING STRUCTURE LIBRARY
load graphite
% When making the material file, keep the space group in mind. For
% instance, you might have 5 atoms in your stoichiometry with a 1:1:3
% ration, but because of the space group, the number of atoms and positions
% in the unit cell is dictated by an algorithm specific to the group. That
% is why, if the material requires 4 positions for the same atom in the
% unit cell due to its symmetry type, then from 5 atoms you will end with
% 20 atoms and you will need 20 inputs in Lattice.AtomicPositions and
% Lattice.AtomicTypes.
% DEFINE DIRECTION OF CRYSTAL NORMAL
Lattice.Normal = [0 0 1];
% DEFINE X-RAYS
Probe.Type = 'xrays';
Probe.Energy = 8048.3; % [eV]
% DEFINE DETECTOR
% Detector is not used for Powder Diffraction.
% MILLER INDICES TO LOOP OVER
hkl = 6;
%Ignores all peaks which have intensities under 1.
Threshold = 1;
%Resolution between peaks.
%Defines width of the peak. Decrease to improve resolution.
Separation = 0.1;
% MAIN FUNCTION
I = GenerateIntensity2theta(Lattice, Probe,1,hkl,Threshold, Separation)