[ Identification | Description | Input parameters | Output parameters | Links ]
Isotropic_Sqw
ComponentAn isotropic sample handling multiple scattering and including as input the dynamic structure factor of the chosen sample (e.g. from Molecular Dynamics). Handles elastic/inelastic, coherent and incoherent scattering - depending on the input S(Q,w) - with multiple scattering and absorption. Only the norm of Q is handled (not the vector), and thus suitable for liquids, gazes, amorphous and powder samples. If incoherent S(Q,w) file is specified as empty (0 or "") then the scattering is constant isotropic (Vanadium like). In case you only have one S(q,w) data containing both coherent and incoherent contributions you should e.g. use 'file_coh' and set 'sigma_coh' to the total scattering cross section. The S(q,w) are used as probability tables, and the intensity is weighted to the S(q,w) integral value. This integral value may be set from: norm_sqw=-1 compute |S| from the density (refer to the Component Manual) norm_sqw>0 use |S| = norm_sqw norm_sqw=0 compute |S| = \int S(q,w) dq dw from raw data files The default and recommanded normalization is 'norm_sqw=-1'. Focusing on the relevant [q,w] data range corresponding to the instrument setting can speed-up substantially the computation and prevent some neutrons to be removed when energy transfert is higher than actual neutron energy. An automatic such process may be activated when 'auto_qw=1' (recommanded). Additionally, for single order scattering (order=1), you may restrict the vertical spreading of the scattering area using d_phi parameter. If you use this component and produce valuable scientific results, please cite authors with references bellow (in Links). Sample shape: Sample shape may be a cylinder, a sphere, a box or a hollow cylinder/sphere. box/plate: xwidth x yheight x zthick (thickness=0) hollow box/plate:xwidth x yheight x zthick and thickness>0 cylinder: radius_o x yheight (radius_i=0) hollow cylinder: radius_o x yheight and radius_i>0 or thickness>0 sphere: radius_o (yheight=0 radius_i=0) hollow sphere: radius_o and radius_i>0 or thickness>0 (yheight=0) Concentric components: This component has the ability to contain other components when used in hollow cylinder geometry (namely sample environment, e.g. cryostat and furnace structure). Such component 'shells' should be split into input and output side surrounding the 'inside' components. First part must then use 'concentric=1' flag to enter the inside part. The component itself must be repeated to mark the end of the concentric zone. The number of concentric shells and number of components inside is not limited. COMPONENT S_in = Isotropic_Sqw(Sqw_coh="Al.laz", concentric=1, ...) AT (0,0,0) RELATIVE sample_position COMPONENT something_inside ... // e.g. the sample itself or other walls ... COMPONENT S_out = Isotropic_Sqw(Sqw_coh="Al.laz", ...) AT (0,0,0) RELATIVE sample_position Sqw file format: File format for S(Q,w) (coherent and incoherent) should contain 3 numerical blocks, defining q axis values (vector), then energy axis values (vector), then a matrix with one line per q axis value, containing Sqw values for each energy axis value. Comments (starting with '#') and non numerical lines are ignored and used to separate blocks. Sampling must be regular. Example: # q axis values # vector of m values in Angstroem-1 0.001000 .... 3.591000 # w axis values # vector of n values in meV 0.001391 ... 1.681391 # sqw values (one line per q axis value) # matrix of S(q,w) values (m rows x n values), one line per q value, 9.721422 10.599145 ... 0.000000 10.054191 11.025244 ... 0.000000 ... 0.000000 ... 3.860253 See for instance file He4.sqw. Such files may be obtained from e.g. INX, Nathan, Lamp and IDA softwares, as well as Molecular Dynamics. Powder file format: Files for coherent elastic powder scattering may also be used. Format specification follows the same principle as in the PowderN component, with parameters: powder_format=Crystallographica or powder_format=Fullprof or powder_format=Lazy or powder_format={j,d,F2,DW,Delta_d/d,1/2d,q,F} (column indexes 1:n) or column indexes (starting from 1) given as comments in the file header (e.g. '#column_j 4'). Refer to the PowderN component for more details. Delta_d/d and Debye-Waller factor may be specified for all lines with the 'powder_Dd' and 'powder_DW' parameters. Additionally a special [q,Sq] format is also defined with: powder_format=qSq for which column 1 is 'q' and column 2 is 'S(q)'. Examples: 1- Vanadium-like incoherent elastic scattering Isotropic_Sqw(V_rho=1/13.827, sigma_abs=5.08, sigma_inc=4.935, sigma_coh=0) 2- liq-4He parameters Isotropic_Sqw(..., Sqw_coh="He4.sqw", T=10) 3- powder sample Isotropic_Sqw(..., Sqw_coh="Al.laz", save_sqw=1) %BUGS: When used in concentric mode, multiple bouncing scattering (traversing the hollow part) is not taken into account. %VALIDATION For Vanadium incoherent scattering mode, V_sample, Single_crystal and Isotropic_Sqw produce equivalent results, eventhough the two later are more accurate (geometry, multiple scattering).
Name | Unit | Description | Default |
---|---|---|---|
Sqw_coh | str | Name of the file containing the values of Q, w and S(Q,w) Coherent part; Q in Angs-1, E in meV, S(q,w) in meV-1 Use 0 or "" to disable. | 0 |
Sqw_inc | str | Name of the file containing the values of Q, w and S(Q,w) Incoherent (self) part. Use 0 or "" to scatter isotropically (V-like). | 0 |
powder_format | no quotes | name or definition of column indexes in file | Undefined |
radius_i | m | Inner radius of sample in (x,z) plane. cylinder/sphere. | 0 |
radius_o | m | Outer radius of sample in (x,z) plane. cylinder/sphere. | 0 |
thickness | m | Thickness of hollow sample (overrides radius_i). | 0 |
xwidth | m | width for a box sample shape | 0 |
yheight | m | Height of sample in vertical direction for box/cylinder shapes | 0 |
zthick | m | thickness for a bulk box sample shape | 0 |
qmin | Angs-1 | Minimum Q value to use in S(q,w). | 0 |
qmax | Angs-1 | Maximum Q value to use in S(q,w). | 0 |
wmin | meV | Minimum Energy value to use in S(q,w). | 0 |
wmax | meV | Maximum Energy value to use in S(q,w). | 0 |
auto_qw | 1 | When set to 1, the [q,w] range will automatically be tuned to optimal setting whenever required. | 0 |
threshold | 1 | Value under which S(Q,w) is not accounted for. to set according to the S(Q,w) values ,i.e. not too low. | 1e-10 |
order | 1 | Limit multiple scattering up to given order | 0 |
T | K | Temperature of sample, detailed balance | 0 |
verbose | 1 | Verbosity level (0:silent, 1:normal, 2:verbose, 3:debug). | 1 |
d_phi | deg | scattering vertical angular spreading (usually the heigh of the next component/detector). Use 0 for full space. This is only relevant for single scattering (order=1). | 0 |
concentric | 1 | indicates that this component in hollow cylinder geometry may contain other components. It should be then duplicated after the inside part. See description for an example. | 0 |
V_rho | AA-3 | Density of atoms (nb atoms/unit cell V_0). | 0 |
sigma_abs | barns | Absorption cross-section at 2200 m/s. 0 to disables | 0 |
sigma_coh | barns | Coherent Scattering cross-section. 0 to disable | 0 |
sigma_inc | barns | Incoherent Scattering cross-section. 0 to disable | 0 |
save_sqw | 1 | When set to 1, saves S(q), DOS and S(q,w) as monitors. | 0 |
norm_sqw | 1 | When positive, the value is used as the integral of S(q,w) Null value uses computed integral. Any negative value will normalize on the density. | -1 |
powder_Dd | 1 | global Delta_d/d spreading, or 0 if ideal. | 0 |
powder_DW | 1 | global Debey-Waller factor, if not in |F2| or 1. | 0 |
powder_Vc | AA^3 | volume of the unit cell | 0 |
density | g/cm^3 | density of material. V_rho=density/weight/1e24*N_A | 0 |
weight | g/mol | atomic/molecular weight of material | 0 |
Name | Unit | Description | Default |
---|---|---|---|
VarSqw | - | .type interaction type of event 'c' (coherent), 't' (transmitted) 'i' (incoherent) 'v' (isotropic incoherent, Vanadium-like). | |
columns |
Isotropic_Sqw.comp
.
[ Identification | Description | Input parameters | Output parameters | Links ]
Generated automatically by McDoc, Peter Willendrup <peter.willendrup@risoe.dk> / Sun Jul 11 17:36:12 2010