Verification of the McStas code using two double axis neutron diffractometers

Abstract

The Monte Carlo ray-tracing simulation package McStas is verified by modeling two double axis neutron diffractometers, KARL and KANDI-II, located at the Israel Research reactors IRR-1 and IRR-2, respectively. The advantage of simulating these instruments is their simple design and small number of components. We focus on the simulation of a Soller slit collimator as well as Pyrolytic Graphite and Copper single crystal monochromators. The simulation is compared to neutron flux measurements, performed via activation analysis along the neutron beam as well as observed diffraction patterns from different powder samples. Calculated diffraction patterns are found to agree with observations to within 20%–50% in reflection peak width, while calculated attenuation of the neutron beam agrees to within 20% when the entire instrument is considered, and can show a discrepancy of up to 60% for individual components. Finally, a possible optimization of the neutron flux in the KARL diffractometer is presented.