Toolkit implementation to exchange phase-space files between IAEA and MCNP6 monte Carlo code format

PURPOSE

Some Monte Carlo simulation codes can read and write phase space files in IAEA format, which are used to characterize accelerators, brachytherapy seeds and other radiation sources. Moreover, as the format has been standardized, these files can be used with different simulation codes. However, MCNP6 has not still implemented this capability, which complicate the studies involving this kind of sources and the reproducibility of results among independent researchers. Therefore, the purpose of this work is to develop a tool to perform conversions between IAEA and MCNP6 phase space files formats, to be used for Monte Carlo simulations.

MATERIALS AND METHODS

This paper presents a toolkit written in C language that uses the IAEA libraries to convert phase space files between IAEA and MCNP6 format and vice versa. To test the functionality of the provided tool, a set of verification tests has been carried out. In addition, a linear accelerator treatment has been simulated with the PENELOPE library using the PenEasy framework, which is already capable to read and write IAEA phase space files, and MCNP6 using the developed tools.

RESULTS

Both codes show compatible depth dose curves and profiles in a water tank, demonstrating that the conversion tools work properly. Moreover, the phase space file formats have been converted from IAEA to MCNP6 format and back again to IAEA format, reproducing the very same results.

CONCLUSION

The toolkit developed in this work offers MCNP6 scientific community an external and validated program able to convert phase space files in IAEA format to MCNP6 internal format and use them for Monte Carlo applications. Furthermore, the developed tools provide also the reverse conversion, which allow sharing MCNP6 results with users of other Monte Carlo codes. This capability in the MCNP6 ecosystem provides to the scientific community the ability not only to share radiation sources, but also to facilitate the reproducibility among different groups using different codes via the standard format specified by the IAEA.

Study of the photoneutron generation caused by a LinAc Beryllium window with a 6 MeV treatment beam

In most conventional radiation therapy treatments, special attention is payed for neutron contamination when working with energy beams above 8 MeV and generally it is only considered for shielding requirements, not for dose study in patients or employees. The present work is focused on studying the unwanted generated photoneutrons in a Medical Linear Accelerator (LinAc) Varian TrueBeam using a 6 MeV radiation treatment beam. To that, Monte Carlo (MC) simulation code MCNP6.1.1 was used. This version of the code allows the use of unstructured mesh geometries as a novelty, offering more reliable results and higher speed computation. The particularity of the studied LinAc is the presence of a beryllium filter at the treatment head. Since Beryllium causes photonuclear reactions (γ, n) at energies much lower than other LinAc composing materials, this work aims to analyze if this type of units, when using low energy treatment beams (6 MeV), produce neutron pollution and to ensure that this unwanted radiation can be considered negligible.