Rejection-based sampling of inelastic neutron scattering

Distributions of inelastically scattered neutrons can be quantum dynamically described by a scattering kernel. We present an accurate and computationally effcient rejection method for sampling a given scattering kernel of any isotropic material. The proposed method produces continuous neutron energy and angular distributions, typically using just a single interpolation per sampling.

We benchmark the results of this method against those from accurate analytical models and one of the major neutron transport codes. We also show the results of applying this method to the conventional discrete double differential cross sections.