Empirical description and Monte Carlo simulation of fast neutron pencil beams as basis of a treatment planning system.

Abstract

The fast neutron beam, used for fast neutron therapy in Essen, is produced by the nuclear reaction of a 14 MeV cyclotron-based deuteron beam on a thick beryllium target. The resulting neutron beam has a continuous energy spectrum with a mean and a maximum energy equal to 5.5 and 18 MeV, respectively. The dose delivered to the patient is computed by a treatment planning system (TPS) based on an empirical model, in which the dose components (neutron and photon) are described by analytical functions. In order to improve the dose calculation, and thus to use the fast neutron beam for other applications (e.g., Boron Neutron Capture Enhancement of Fast Neutron Therapy), in this work we aim to develop a new TPS. For this purpose, a model based on pencil beams of mono-energetic neutrons has been created. The neutron energy ranged from 0.25 MeV up to 17.25 MeV by steps of 0.5 MeV in order to cover the energy range of the Essen facility. The Monte Carlo method was then used to simulate the transport of neutrons within such pencil beams in a homogeneous water phantom. By using Monte Carlo techniques, it is possible to distinguish the energy deposition due to a primary collision in water to that due to scattered neutrons. The energy deposition due to pencil beams of 2.224 MeV photons, coming from hydrogen neutron capture reaction in the phantom or in the collimator, was also determined. In order to complete this work, air filled cylinders have been introduced in the water phantom. It is shown that the resulting depth dose curves for primary neutrons can be easily derived using the homogeneous phantom, and that the description of the effect on scattered neutron dose distribution is more complex. In this work we demonstrate the relevance of Monte Carlo simulations of mono-energetic neutron pencil beams for purposes of neutron treatment planning. Some additional work is still required to describe a clinical situation (continuous energy neutron spectrum) as well as to experimentally validate the method described here.