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.
- Thermal neutron scattering
- Monte Carlo simulation
- Rejection method
- Isotropic material
Cai, X-X., Kittelmann, T., Klinkby, E., & Damián, J. I. M. (2019). Rejection-based sampling of inelastic neutron scattering. Journal of Computational Physics, 380, 400-407. https://doi.org/10.1016/j.jcp.2018.11.043