General powder/polycrystalline sample with neutronmatter interaction based
on neutron cross section calculations of a unit cell
Identification
Author: Mirko Boin
Origin: HelmholtzZentrum Berlin fuer Materialien und Energie (Germany)
Date: Jan 2019
Version: 1.3
Description
Features:
 coherent and incoherent scattering, absorption and transmission
 multiple scattering
 wavelengthdependent neutron cross section calculation
 unit cell definition (via input file)
 focussing option for scattered neutrons (via d_phi)
 strain gradient simulation for axial symmetric samples (...to be extended)
Geometry is a powder filled cylinder or a box defined by radius or xwidth,
yheight, zthick respectively. The component handles coherent and incoherent
scattering (also multiple scattering), absorption and transmission. Hence, it
is suitable for diffraction (scattering) and imaging (transmission) instruments
at the same time. In order to improve the neutron statistics at the detector
these individual features can be enabled/disabled using TransOnly, IncohScat,
MultiScat.
For example, if scattering shall not be monitored (e.g. imaging), then TransOnly
can be set to 1 to let the component calculate the neutron transmission through
the sample only. Likewise, incoherent scattering can be switched on and off.
The decision of whether a neuron shall be scattered, absorbed or transmitted is
based on the wavelengthdependent neutron cross sections (nxs), i.e. the neutron
velocity is used here. These cross sections represent the individual scattering
and absorption probabilities for each individual neutron. However, the calculation
of neutron cross sections depends on the material defined for this sample component.
Therefore, several input parameters such as the crystal structure and the atoms
involved (isotope mass, coherent scattering length, ...). This information is
provided by an input file with the following format:
Al.nxs
# define the unit cell parameters:
# space_group  the space group number or Hermann or Hall symbol [string]
# lattice_a, ...b, ...c  the lattice spacings a,b,c [angstrom]
# lattice_alpha, ...beta, ...gamma  the lattice angles alpha,beta,gamma [degree]
# debye_temp  the Debye temperature [K]
space_group = F 4 2 3 # space group number is also allowed (= 225)
lattice_a = 4.049
lattice_b = 4.049
lattice_c = 4.049
lattice_alpha = 90
lattice_beta = 90
lattice_gamma = 90
debye_temp = 429.0
# add atoms to the unit cell:
# notation is "atom_number = name b_coh sigma_inc sigma_abs_2200 molar_mass x y z"
# name  labels the current atom/isotope [string]
# b_coh  the coherent scattering length [fm]
# sigma_inc  the incoherent scattering cross section [barns]
# sigma_abs_2200  the absorption cross sect. at 2200 m/s [barns]
# molar_mass  the Molar mass [g/mol]
# debye_temp  the Debye temperature [K]
# x y z  the Wyckoff postion of the atom inside the unit cell
The above example defines a pure (fcc) aluminium. Other example files are provided.
Alternatively, a fallback mechanism exists to simulate an alphaFe (bcc) sample
if an input file is missing. But have a look for warning messages during simulation to
make sure you are not running the fallback alphaFe simulation accidently!
CHANGELOG:
v1.0  first official release with McStas 2.0
v1.1  enabled strain distribution
v1.2  use updated libnxs (global Debye temp > parameter files changed)
 corrected sigma incoherent
 updated towards McStas 2.1 conventions
v1.3  bugfixes
For details, please contact the author and visit the nxs website and/or the repository.
NXS Website: http://www.helmholtzberlin.de/people/boin/nxs_en.html
NXS Repository: http://bitbucket.org/mirkoboin/nxs
WARNING: This is a contributed Component.
Input parameters
Parameters in boldface are required;
the others are optional.
Name
Unit
Description
Default
TransOnly

enable/disable (=1/=0) the transmission only option
0
IncohScat

enable/disable (=1/=0) incoherent scattering
1
MultiScat

enable/disable (=1/=0) multiple scattering
1
xwidth
m
horizontal dimension of sample, as a width
0
yheight
m
vertical dimension of sample, as a height
0
zthick
m
thickness of sample
0
radius
m
radius of sample in (x,z) plane
0
thickness
0
geometry
0
nxsFileName

filename of the unit cell definition
""
max_hkl

maximum (hkl) indices to consider for the nxs calculation, e.g. max_hkl=3 means (hkl)=333
8
d_phi
deg,0180
=0: Angle corresponding to the vertical angular range
to focus to, e.g. detector height. 0 for no focusing