Demonstration of three-dimensional deterministic radiation transport theory dose distribution analysis for boron neutron capture therapy

The Monte Carlo stochastic simulation technique has traditionally been the only well-recognized method for computing three-dimensional radiation dose distributions in connection with boron neutron capture therapy (BNCT) research. A deterministic approach to this problem would offer some advantages over the Monte Carlo method. This paper describes an application of a deterministic method to analytically simulate BNCT treatment of a canine head phantom using the epithermal neutron beam at the Brookhaven medical research reactor (BMRR). Calculations were performed with the TORT code from Oak Ridge National Laboratory (ORNL), an implementation of the discrete ordinates, or Sn method. Calculations were from first principles and used no empirical correction factors. The phantom surface was modeled by flat facets of approximately 1 cm2. The phantom interior was homogeneous. Energy-dependent neutron and photon scalar fluxes were calculated on a 32 x 16 x 22 mesh structure with 96 discrete directions in angular phase space. The calculation took 670 min on an Apollo DN10000 workstation. The results were subsequently integrated over energy to obtain full three-dimensional dose distributions. Isodose contours and depth-dose curves were plotted for several separate dose components of interest. Phantom measurements were made by measuring neutron activation (and therefore neutron flux) as a function of depth in copper-gold alloy wires that were inserted through catheters placed in holes drilled in the phantom. Measurements agreed with calculations to within about 15%. The calculations took about an order of magnitude longer than comparable Monte Carlo calculations but provided various conveniences, as well as a useful check.