Single-crystal dispersion: S(Q,E) = S(Q) delta(E - E(Q)) with analytic E(Qx,Qy,Qz).
Author: Fahima Islam (McStas port of the MCViNE E_vQ_Kernel; MCViNE by J. Y. Y. Lin et al.)
Origin: MCViNE, https://github.com/mcvine/mcvine (mccomponents/lib/kernels/sample/E_vQ_Kernel.icc)
Date: 2026-09-29
Excitation with dispersion E(Qx,Qy,Qz) [meV] and intensity S(Qx,Qy,Qz), given as expressions of the Q vector in the sample (component) frame [Å^-1]. For a random final direction, all solutions kf of Ei - Ef(kf) = E(ki - kf) with Ef in [Ei-Emax, Ei] are found (nsteps sub-intervals + Ridders’ method) and one is picked; the weight contains the Jacobian of the delta function.
Functions are given as strings evaluated at run time (MCViNE uses fparser): + - * / ^ (or **), % , comparisons, && ||, if(c,a,b), sin cos tan asin acos atan sinh cosh tanh exp log log10 log2 sqrt abs floor ceil int sign cbrt, pow atan2 min max hypot fmod, constants pi and e.
Transport (same as MCViNE HomogeneousNeutronScatterer): on the first passage the neutron is forced to scatter at a uniformly chosen depth, weighted by path length, attenuation exp(-(mu+sigma)x) and the scattering coefficient; on exit it is attenuated by exp(-(mu+sigma)L_out). With order>1 further scatterings are sampled analogically (truncated exponential). With p_transmit>0 that fraction of events is kept as the attenuated direct beam. mu(v) = absorption_coefficient*2200/v. Events for which the kernel cannot scatter (kinematically forbidden) are absorbed.
Geometry: box (xwidth,yheight,zdepth), cylinder (radius,yheight), hollow cylinder (+thickness) or sphere (radius only). All vectors (Q, targets, reciprocal vectors, atom positions) are in the component’s local frame.
Uses the shared runtime share/mcvine-lib.h/.c.
Example: MCViNE_E_vQ(E_Q="20*(sin(Qx*1.57)^2+sin(Qy*1.57)^2+sin(Qz*1.57)^2)", S_Q="1", Emax=60, scattering_coefficient=10, xwidth=0.02, yheight=0.02, zdepth=0.02)
Parameters in boldface are required; the others are optional.
|
Name |
Unit |
Description |
Default |
|
|
|||
|
E_Q |
str |
E(Qx,Qy,Qz) expression [meV] |
"10" |
|
S_Q |
str |
S(Qx,Qy,Qz) expression |
"1" |
|
Emax |
meV |
Maximum energy transfer considered |
10 |
|
nsteps |
1 |
Number of sub-intervals for the kf root search (MCViNE: 1000) |
1000 |
|
xacc |
Å\(^{-1}\) |
Root accuracy in kf (MCViNE: 1e-3) |
1e-3 |
|
absorption_coefficient |
m\(^{-1}\) |
Absorption coefficient at 2200 m/s (MCViNE convention; scales as 1/v) |
0 |
|
scattering_coefficient |
m\(^{-1}\) |
Scattering coefficient |
0 |
|
sigma_abs |
barn |
Alternative: absorption cross section per unit cell at 2200 m/s (used when Vc>0) |
0 |
|
sigma_scat |
barn |
Alternative: scattering cross section per unit cell (used when Vc>0) |
0 |
|
Vc |
Å\(^{3}\) |
Unit cell volume; when >0, coefficients are computed from sigma_abs, sigma_scat |
0 |
|
radius |
m |
Outer radius of a cylinder (with yheight) or of a sphere (yheight=0) |
0 |
|
xwidth |
m |
Width of a box sample |
0 |
|
yheight |
m |
Height of a box or cylinder sample |
0 |
|
zdepth |
m |
Depth of a box sample |
0 |
|
thickness |
m |
Wall thickness of a hollow cylinder (0: filled) |
0 |
|
pack |
1 |
Packing factor (scales absorption and scattering coefficients) |
1 |
|
p_transmit |
1 |
Monte Carlo fraction of events kept as unscattered transmitted beam (0: always scatter) |
0 |
|
order |
1 |
Maximum number of scattering events per neutron (1: single scattering) |
1 |
|
|
|||
|
|
|||
|
|
|||
Component source code found in file MCViNE_E_vQ.comp.
MCViNE documentation: https://mcvine.github.io