Perfect, laue crystal with common cubic structures (diamond, fcc, or bcc, and others if symmetry form factor multipliers provided explicitly)
Author: Marcus H Mendenhall, NIST
Origin: Marcus H. Mendenhall, NIST, Gaithersburg, MD, USA
Date: June, 2017
NIST_Laue_crystal_BC.comp is written by Marcus H. Mendenhall, NIST, Gaithersburg, MD, USA
It is based on the full vector math and exact solution of the dispersion relation in
Batterman and Cole, Reviews of Modern Physics 36 number 3, page 681, July 1964
Contributed by the National Institute of Standards and Technology; not subject to copyright in the United States.
This is not an official contribution, in that the results are in no way certified by NIST.
Perfect crystal with common cubic structures (diamond, fcc, or bcc, and others if symmetry form factor multipliers provided explicitly)
Reads atomic form factors from a data input file.
The Laue_Crystal code reflects rays in an ideal geometry, does not include surface imperfections or mosaicity.
The crystal is positioned such that the long axis of the crystal surface coincides with
z-axis and the outer normal to the crystal surface is along +y.
The ratio of the transmitted beam and forward-diffracted Borrman-effect beam is a hack. The sum of the
two is exactly right, but the actual ratio depends critically on geometry, and I just put in a wild estimate
to allow one to demonstrate what the Borrmann effect looks like. If this is turned on, the displacement of the
transmitted beam and forward diffracted beam at the back side of the crystal will be correctly computed.
This displacement is only exact for symmetrical Laue; asymmetrical computation requires more effort,
and is probably not worth it.
The sampling of these processes are controlled by the 3 variables transmission_sampling, forward_diffraction_sampling,
and laue_sampling. Since 99% of uses of this will have the transmitted beam turned off, and use laue diffraction mode,
the values should be just 0,0,1. If the general behavior of the transmitted beams
is interesting, use 1,1,1 which samples all beams equally. Results weights are adjusted for this,
so computed intensities won't be affected.
WARNING: This is a contributed Component.
Parameters in boldface are required;
the others are optional.
zdepth (length) of the crystal.
width of the crystal.
along y-axis, the surface normal
[m] thickness of crystal
unit cell volume
"FormFactors.txt" from McXtrace install, usually
[ ] Si ("Si.txt"), Ge
x component of normal (unit vector) to crystal planes. Vector is usually [0,0,1] for symmetric Laue. Crystal surface itself has normal [0,1,0].
y component of normal (unit vector) to crystal planes.
z component of normal (unit vector) to crystal planes.
Debye-Waller temperature factor, M=B*(sin(theta)/lambda)^2*(2/3), default=silicon at room temp, 0.4632
1 => Mx_crystal_explicit: provide explicit real and imaginary form factor multipliers structure_factor_scale_r, structure_factor_scale_i,
2 => Mx_crystal_diamond: diamond
3 => Mx_crystal_fcc: fcc
4 => Mx_crystal_fcc: bcc
1st Miller index of reflection
2nd Miller index of reflection
3rd Miller index of reflection
real part of complex explicit override of structure factor multiplier for crystal structure if Bragg_crystal_explicit
imaginary part of complex explicit override of structure factor multiplier for crystal structure if Bragg_crystal_explicit