Monte Carlo model of the Studsvik BNCT clinical beam: description and validation

The neutron beam at the Studsvik facility for boron neutron capture therapy (BNCT) and the validation of the related computational model developed for the MCNP-4B Monte Carlo code are presented. Several measurements performed at the epithermal neutron port used for clinical trials have been made in order to validate the Monte Carlo computational model. The good general agreement between the MCNP calculations and the experimental results has provided an adequate check of the calculation procedure. In particular, at the nominal reactor power of 1 MW, the calculated in-air epithermal neutron flux in the energy interval between 0.4 eV-10 keV is 3.24 x 10(9) n cm(-2) s(-1) (+/- 1.2% 1 std. dev.) while the measured value is 3.30 x 10(9) n cm(-20 s(-1) (+/- 5.0% 1 std. dev.). Furthermore, the calculated in-phantom thermal neutron flux, equal to 6.43 x 10(9) n cm(-2) s(-1) (+/- 1.0% 1 std. dev.), and the corresponding measured value of 6.33 X 10(9) n cm(-2) s(-1) (+/- 5.3% 1 std. dev.) agree within their respective uncertainties. The only statistically significant disagreement is a discrepancy of 39% between the MCNP calculations of the in-air photon kerma and the corresponding experimental value. Despite this, a quite acceptable overall in-phantom beam performance was obtained, with a maximum value of the therapeutic ratio (the ratio between the local tumor dose and the maximum healthy tissue dose) equal to 6.7. The described MCNP model of the Studsvik facility has been deemed adequate to evaluate further improvements in the beam design as well as to plan experimental work.