Monte Carlo studies on photon interactions in radiobiological experiments

Development of DynamicMC for PHITS Monte Carlo package

Previously, we have developed DynamicMC for modeling relative movement of Oak Ridge National Laboratory phantom in a radiation field for the Monte Carlo N-Particle package (Health Physics. 2023,124(4):301-309). Using this software, three-dimensional dose distributions in a phantom irradiated by a certain mono-energetic (Mono E) source can be deduced through its graphical user interface. In this study, we extended DynamicMC to be used in combination with the Particle and Heavy Ion Transport code System (PHITS) by providing it with a higher flexibility for dynamic movement for an anthropomorphic phantom. For this purpose, we implemented four new functions into the software, which are (1) to generate not only Mono E sources but also those having an energy spectrum of an arbitrary radioisotope (2) to calculate the absorbed doses for several radiologically important organs (3) to automatically average the calculated absorbed doses along the path of the phantom and (4) to generate user-defined slab shielding materials. The first and third items utilize the PHITS-specific modalities named radioisotope-source and sumtally functions, respectively. The computational cost and complexity can be dramatically reduced with these features. We anticipate that the present work and the developed open-source tools will be in the interest of nuclear radiation physics community for research and teaching purposes.

On the effectiveness of proton boron fusion therapy (PBFT) at cellular level

The present work introduced a framework to investigate the effectiveness of proton boron fusion therapy (PBFT) at the cellular level. The framework consisted of a cell array generator program coupled with PHITS Monte Carlo package with a dedicated terminal-based code editor that was developed in this work. The framework enabled users to model large cell arrays with normal, all boron, and random boron filled cytoplasm, to investigate the underlying mechanism of PBFT. It was found that alpha particles and neutrons could be produced in absence of boron mainly because of nuclear reaction induced by proton interaction with ¹⁶O, ¹²C and ¹⁴N nuclei. The effectiveness of PBFT is highly dependent on the incident proton energy, source size, cell array size, buffer medium thickness layer, concentration and distribution of boron in the cell array. To quantitatively assess the effectiveness of PBFT, of the total energy deposition by alpha particle for different cases were determined. The number of alpha particle hits in cell cytoplasm and nucleus for normal and 100 ppm boron were determined. The obtained results and the developed tools would be useful for future development of PBFT to objectively determine the effectiveness of this treatment modality.

Chemical Overview of Gel Dosimetry Systems: A Comprehensive Review

Advances in radiotherapy technology during the last 25 years have significantly improved both dose conformation to tumors and the preservation of healthy tissues, achieving almost real-time feedback by means of high-precision treatments and theranostics. Owing to this, developing high-performance systems capable of coping with the challenging requirements of modern ionizing radiation is a key issue to overcome the limitations of traditional dosimeters. In this regard, a deep understanding of the physicochemical basis of gel dosimetry, as one of the most promising tools for the evaluation of 3D high-spatial-resolution dose distributions, represents the starting point for developing new and innovative systems. This review aims to contribute thorough descriptions of the chemical processes and interactions that condition gel dosimetry outputs, often phenomenologically addressed, and particularly formulations reported since 2017.

Development of PHITS graphical user interface for simulation of positron emitting radioisotopes production in common biological materials during proton therapy

The Monte Carlo (MC) method is a powerful tool for modeling nuclear radiation interaction with matter. A variety of MC software packages has been developed, especially for applications in radiation therapy. Most widely used MC packages require users to write their own input scripts for their systems, which can be a time consuming and error prone process and requires extensive user experience. In the present work, we have developed a graphical user interface (GUI) bundled with a custom-made 3D OpenGL visualizer for PHITS MC package. The current version focuses on modeling proton induced positron emitting radioisotopes, which in turn can be used for verification of proton ranges in proton therapy. The developed GUI program does not require extensive user experience. The present open-source program is distributed under GPLv3 license that allows users to freely download, modify, recompile and redistribute the program.

Proton range monitoring using 13N peak for proton therapy applications

The Monte Carlo method is employed in this study to simulate the proton irradiation of a water-gel phantom. Positron-emitting radionuclides such as ¹¹C, ¹⁵O, and ¹³N are scored using the Particle and Heavy Ion Transport Code System Monte Carlo code package. Previously, it was reported that as a result of ¹⁶O(p,2p2n)¹³N nuclear reaction, whose threshold energy is relatively low (5.660 MeV), a ¹³N peak is formed near the actual Bragg peak. Considering the generated ¹³N peak, we obtain offset distance values between the ¹³N peak and the actual Bragg peak for various incident proton energies ranging from 45 to 250 MeV, with an energy interval of 5 MeV. The offset distances fluctuate between 1.0 and 2.0 mm. For example, the offset distances between the ¹³N peak and the Bragg peak are 2.0, 2.0, and 1.0 mm for incident proton energies of 80, 160, and 240 MeV, respectively. These slight fluctuations for different incident proton energies are due to the relatively stable energy-dependent cross-section data for the ¹⁶O(p,2p2n)¹³N nuclear reaction. Hence, we develop an open-source computer program that performs linear and non-linear interpolations of offset distance data against the incident proton energy, which further reduces the energy interval from 5 to 0.1 MeV. In addition, we perform spectral analysis to reconstruct the ¹³N Bragg peak, and the results are consistent with those predicted from Monte Carlo computations. Hence, the results are used to generate three-dimensional scatter plots of the ¹³N radionuclide distribution in the modeled phantom. The obtained results and the developed methodologies will facilitate future investigations into proton range monitoring for therapeutic applications.

MCHP (Monte Carlo + Human Phantom): Platform to facilitate teaching nuclear radiation physics

Some concepts in nuclear radiation physics are abstract and intellectually demanding. In the present paper, an “MCHP platform” (MCHP was an acronym for Monte Carlo simulations + Human Phantoms) was proposed to provide assistance to the students through visualization. The platform involved Monte Carlo simulations of interactions between ionizing radiations and the Oak Ridge National Laboratory (ORNL) adult male human phantom. As an example to demonstrate the benefits of the proposed MCHP platform, the present paper investigated the variation of the absorbed photon dose per photon from a ¹³⁷Cs source in three selected organs, namely, brain, spine and thyroid of an adult male for concrete and lead shields with varying thicknesses. The results were interesting but not readily comprehensible without direct visualization. Graphical visualization snapshots as well as video clips of real time interactions between the photons and the human phantom were presented for the involved cases, and the results were explained with the help of such snapshots and video clips. It is envisaged that, if the platform is found useful and effective by the readers, the readers can also propose examples to be gradually added onto this platform in future, with the ultimate goal of enhancing students’ understanding and learning the concepts in an undergraduate nuclear radiation physics course or a related course.

Advances on inorganic scintillator-based optic fiber dosimeters

This article presents a new perspective on the development of inorganic scintillator-based fiber dosimeters (IOSFDs) for medical radiotherapy dosimetry (RTD) focusing on real-time in vivo dosimetry. The scintillator-based optical fiber dosimeters (SFD) are compact, free of electromagnetic interference, radiation-resistant, and robust. They have shown great potential for real-time in vivo RTD. Compared with organic scintillators (OSs), inorganic scintillators (IOSs) have larger X-ray absorption and higher light output. Variable IOSs with maximum emission peaks in the red part of the spectrum offer convenient stem effect removal. This article outlines the main advantages and disadvantages of utilizing IOSs for SFD fabrication. IOSFDs with different configurations are presented, and their use for dosimetry in X-ray RT, brachytherapy (BT), proton therapy (PT), and boron neutron capture therapy (BNCT) is reviewed. Challenges including the percentage depth dose (PDD) deviation from the standard ion chamber (IC) measurement, the angular dependence, and the Cherenkov effect are discussed in detail; methods to overcome these problems are also presented.

Radiation-induced rescue effect

Radiation-induced rescue effect (RIRE) refers to the phenomenon in which detrimental effects in targeted irradiated cells are reduced upon receiving feedback signals from partnered non-irradiated bystander cells, or from the medium previously conditioning these partnered non-irradiated bystander cells. For convenience, in the current review we define two types of RIRE: (i) Type 1 RIRE (reduced detrimental effects in targeted cells upon receiving feedback signals from bystander cells) and (ii) Type 2 RIRE (exacerbated detrimental effects in targeted cells upon receiving feedback signals from bystander cells). The two types of RIRE, as well as the associated mechanisms and chemical messengers, have been separately reviewed. The recent report on the potential effects of RIRE on the traditional colony-formation assays has also been reviewed. Finally, future priorities and directions for research into RIRE are discussed.