Instituto de Seguridad Industrial, Radiofísica y Medioambiental (ISIRYM), Universitat Politècnica de València, València, Spain.
Int J Radiat Biol. 2023;99(3):373-383. doi: 10.1080/09553002.2022.2110296. Epub 2023 Jan 17.
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.
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.
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.
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.
一些蒙特卡罗模拟代码可以读取和写入国际原子能机构(IAEA)格式的相空间文件,这些文件用于描述加速器、近距离治疗放射源和其他辐射源。此外,由于该格式已经标准化,这些文件可以与不同的模拟代码一起使用。然而,MCNP6 尚未实现此功能,这使得涉及此类源的研究以及独立研究人员之间结果的可重复性变得复杂。因此,这项工作的目的是开发一种工具,用于在 IAEA 和 MCNP6 相空间文件格式之间进行转换,以便用于蒙特卡罗模拟。
本文介绍了一个用 C 语言编写的工具包,该工具包使用 IAEA 库在 IAEA 和 MCNP6 格式之间转换相空间文件,反之亦然。为了测试提供的工具的功能,已经进行了一组验证测试。此外,使用 PenEasy 框架,已使用 PENELOPE 库对线性加速器治疗进行了模拟,PenEasy 框架已经能够读取和写入 IAEA 相空间文件,并且可以使用开发的工具读取和写入 MCNP6 相空间文件。
两种代码在水箱中都显示出兼容的深度剂量曲线和剖面图,表明转换工具运行正常。此外,已将相空间文件从 IAEA 格式转换为 MCNP6 格式,然后再转换回 IAEA 格式,重现了完全相同的结果。
这项工作中开发的工具包为 MCNP6 科学界提供了一个外部且经过验证的程序,能够将 IAEA 格式的相空间文件转换为 MCNP6 内部格式,并将其用于蒙特卡罗应用。此外,开发的工具还提供了反向转换,允许与其他蒙特卡罗代码的用户共享 MCNP6 结果。MCNP6 生态系统中的此功能不仅为科学界提供了共享辐射源的能力,还通过 IAEA 规定的标准格式促进了不同群体之间使用不同代码的可重复性。