Scattering kernels for fast neutron therapy treatment planning

The University of Washington (UW) Clinical Neutron Therapy System (CNTS) has been used to treat over 3300 patients. Treatment planning for these patients is currently performed using an MV x-ray model in Pinnacle® adapted to fit measurements of fast neutron output factors, wedge factors, depth-dose and lateral profiles. While this model has provided an adequate representation of the CNTS for 3D conformal treatment planning, later versions of Pinnacle did not allow for isocentric gantry rotation machines with a source-to-axis distance of 150 cm. This restriction limited the neutron model to version 9.0 while the photon and electron treatment planning at the UW had moved on to newer versions. Also, intensity modulated neutron therapy (IMNT) is underdevelopment at the UW, and the Pinnacle neutron model developed cannot be used for inverse treatment planning. We have used the MCNP6 Monte Carlo code system to develop Collapsed Cone (CC) and Singular Value Decomposition (SVD) neutron scattering kernels suitable for integration into the RayStation treatment planning system. Kernels were generated for monoenergetic neutrons with energies ranging from 1 keV to 51 MeV, i.e. the energy range most relevant to the CNTS. Percent depth dose (PDD) profiles computed in RayStation for the CC and SVD kernels are in excellent agreement with each other for depths beyond the beam's dose build-up region (depths greater than about 1.7 cm) for open 2.8 × 2.8 cm², 10.3 × 10.3 cm², and 28.8 × 32.8 cm² fields. Lateral profiles at several depths, as well as the PDD, calculated using the CC kernels in RayStation for the full CNTS energy spectrum pass a 3%/3 mm gamma test for field sizes of 2.8 × 2.8 cm², 10.0 × 10.3 cm², and 28.8 × 32.8 cm².