THERMAL: A routine designed to calculate neutron thermal scattering. Revision 1

THERMAL is designed to calculate neutron thermal scattering that is elastic and isotropic in the center of mass system. At low energy thermal motion will be included. At high energies the target nuclei are assumed to be stationary. The point of transition between low and high energies has been defined to insure a smooth transition. It is assumed that at low energy the elastic cross section is constant in the relative system. At high energy the cross section can be of any form. You can use this routine for all energies where the elastic scattering is isotropic in the center of mass system. In most materials this will be a fairly high energy, e.g., the keV energy range. The THERMAL method is simple, clean, easy to understand, and most important very efficient; on a SUN SPARC-10 workstation, at low energies with thermal scattering it can do almost 6 million scatters a minute and at high energy over 13 million. Warning: This version of THERMAL completely supersedes the original version described in the same report number, dated February 24, 1995. The method used in the original code is incorrect, as explained in this report