Monte Carlo Simulation of Fast Neutron Treatments for Head and Neck Cancers
G Moffitt1*, R Stewart2, G Sandison3, R Emery4, T Jevremovic5, (1) University of Utah, Salt Lake City, Utah (2) University of Washington, Seattle, WA, (3) University of Washington, Seattle, WA, (4) University of Washington, Seattle, Washington, (5) University of Utah, Salt Lake City, UtahSU-E-T-470 Sunday 3:00PM - 6:00PM Room: Exhibit Hall
To quantify the physical characteristics of a therapeutic fast neutron beam at the University of Washington (UW) using a Monte Carlo model validated against measured percentage depth-dose (PDD) and lateral profiles.
MCNPX modeling and simulation of the UW neutron therapy system's treatment head is divided into three stages. In Stage I, 50.5 MeV protons are directed at a Be target. Neutrons produced in the target pass through the primary collimator and flattening filter; photons and neutrons produced in Stage I are saved to a 100x100 cm² phase space file located just below the flattening filter. In Stage II, neutrons and photons are transported through the MLC and saved to a second phase space file. Automatic adjustment of the 40 individual MLC to create regular or irregular field shapes with or without a wedge is controlled by an in-house developed python script. In Stage III, neutral particles are transported from the phase space file at the bottom of the MLC into a 40x40x40 cm³ water phantom (150 cm source to surface distance). PDD curves and lateral dose profiles are computed from dose tallies in rectangular (0.1 cm³) volumes.
PDD simulation for an open 10.4x10.3 cm² field agrees with measurement to within 4% up to 17 cm and within 8% up to 30 cm. A lateral dose profile at a depth of 1.7 cm matches well the overall shape of the profile. A representative simulation (10¹¹ source particles) requires about 150 hours of CPU time on a dual, 2.4 GHz 8-core Intel Xeon server.
The developed MCNPX model is capable of reproducing measured data with sufficient accuracy for clinical applications. Efforts to improve agreement in the penumbra of lateral profiles are underway as well as tests of the model for various sizes of open and wedged fields.