sc_paracrystal¶
Simple cubic lattice with paracrystalline distortion
Parameter |
Description |
Units |
Default value |
---|---|---|---|
scale |
Scale factor or Volume fraction |
None |
1 |
background |
Source background |
cm-1 |
0.001 |
dnn |
Nearest neighbor distance |
Å |
220 |
d_factor |
Paracrystal distortion factor |
None |
0.06 |
radius |
Radius of sphere |
Å |
40 |
sld |
Sphere scattering length density |
10-6Å-2 |
3 |
sld_solvent |
Solvent scattering length density |
10-6Å-2 |
6.3 |
theta |
c axis to beam angle |
degree |
0 |
phi |
rotation about beam |
degree |
0 |
psi |
rotation about c axis |
degree |
0 |
The returned value is scaled to units of cm-1 sr-1, absolute scale.
Definition
Calculates the scattering from a simple cubic lattice with paracrystalline distortion. Thermal vibrations are considered to be negligible, and the size of the paracrystal is infinitely large. Paracrystalline distortion is assumed to be isotropic and characterized by a Gaussian distribution.
The scattering intensity \(I(q)\) is calculated as
where scale is the volume fraction of crystal in the sample volume, \(V_\text{lattice}\) is the volume fraction of spheres in the crystal, \(V_p\) is the volume of the primary particle, \(P(q)\) is the form factor of the sphere (normalized), and \(Z(q)\) is the paracrystalline structure factor for a simple cubic structure.
Note
At this point the GUI does not return \(V_\text{lattice}\) separately so that the user will need to calculate it from the equation given and the appropriate returned parameters.
Warning
As per the equations below, this model will return I(q)=0 for all q if the distortion factor is equal to 0. The model is not meant to support perfect crystals.
Following the derivation from reference [1], as corrected in reference [2], and based on the above figure, the primitive unit cell vectors are \(\vec{a_1},\vec{a_2}\), and \(\vec{a_3}\) which in this case are the same as the conventional unit cell vectors (\(\vec{b_1}\), \(\vec{b_2}\), and \(\vec{b_3}\)) so that
where \(a\) is the lattice parameter which is also in this case the nearest neighbor distance \(D\).
The volume fraction (\(V_{lattice}\)) of spherical particles with radius \(R\) sitting on the sc lattice is then given by:
Now, continuing to follow [1], the structure (lattice) factor \(Z(\vec{q})\) for a 3D paracrystal can be written as:
with
and where \(F_k(\vec{q})\) is the structure factor of the primitive unit cell defined as:
Here, \(\vec{a_k}\) are the primitive unit cell vectors \(\vec{a_1}\), \(\vec{a_2}\), and \(\vec{a_3}\). Furthermore, \(\Delta a_k\) is the isotropic distortion of the lattice point from its ideal position and can be defined by a constant factor \(g=\Delta a / |\vec{a_1}| = \Delta a / |\vec{a_2}| = \Delta a / |\vec{a_3}|=\Delta a/D\).
Finally, assuming the definitions presented in this document, the authors of reference [1] have derived the lattice factors which, substituting \(D\) for the lattice parameter \(a\), are given by:
Finally, the position of the Bragg peaks for the sc lattice are indexed by (reduced q-values):
where for a simple cubic lattice there are no selection rules for h, k, l so that all permutations of them give constructive interference. Thus the peak positions correspond to (just the first 5)
Note
The calculation of \(Z(q)\) is a double numerical integral that must be carried out with a high density of points to properly capture the sharp peaks of the paracrystalline scattering. So be warned that the calculation is slow. Fitting of any experimental data must be resolution smeared for any meaningful fit. This makes a triple integral which may be very slow. If a double-precision GPU with OpenCL support is available this may improve the speed of the calculation.
The 2D (Anisotropic model) is based on the reference below where I(q) is approximated for 1d scattering. Thus the scattering pattern for 2D may not be accurate particularly at low \(q\). For general details of the calculation and angular dispersions for oriented particles see Oriented Particles. Note that we are not responsible for any incorrectness of the 2D model computation.
Source
sc_paracrystal.py
\(\ \star\ \) sc_paracrystal.c
\(\ \star\ \) gauss150.c
\(\ \star\ \) sphere_form.c
\(\ \star\ \) sas_3j1x_x.c
Reference
Authorship and Verification
Author: NIST IGOR/DANSE Date: pre 2010
Last Modified by: Paul Butler Date: Oct 27, 2022
Last Reviewed by: Jonathan Gaudet Date: Nov 2, 2022