Search

Abstract 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 16O, 12C and 14N 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.

The interaction of ionizing radiation with matter is a stochastic process and statistical analysis of such a process would be a crucial step in understanding radioactivity. Geiger–Müller (GM) counter is a widely used radiation detector used in nuclear radiation surveying, which produces counts upon exposure to a radioactive source. There are a variety of multi-purpose software that can be used to perform statistical analysis of measured counts from a GM counter. However, statistical analysis is a lengthy, error prone and time-consuming process, which gets more tedious when the number of measurements increases. In the present work, we have developed an open-source and easy-to-use graphical user interface (GUI) computer program named RadStat for statistical analysis of counts measured by a GM counter. RadStat has its own scripting syntaxes and bundled with gnuplot for quick visualization of output results. We believe the present open-source GUI program would be a useful tool for research and teaching of nuclear radiation physics.

The Monte Carlo method is employed in this study to simulate the proton irradiation of a water-gel phantom. Positron-emitting radionuclides such as 11C, 15O, and 13N are scored using the Particle and Heavy Ion Transport Code System Monte Carlo code package. Previously, it was reported that as a result of 16O(p,2p2n)13N nuclear reaction, whose threshold energy is relatively low (5.660 MeV), a 13N peak is formed near the actual Bragg peak. Considering the generated 13N peak, we obtain offset distance values between the 13N 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 13N 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 16O(p,2p2n)13N 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 13N 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 13N radionuclide distribution in the modeled phantom. The obtained results and the developed methodologies will facilitate future investigations into proton range monitoring for therapeutic applications.

Abstract A calibration method was proposed in the present work to determine the medium-thickness-dependent proton doses absorbed in cellular components (i.e., cellular cytoplasm and nucleus) in radiobiological experiments. Consideration of the dependency on medium thickness was crucial as the linear energy transfer (LET) of protons could rise to a sharp peak (known as the Bragg peak) towards the end of their ranges. Relationships between the calibration coefficient R vs medium-layer thickness were obtained for incident proton energies of 10, 15, 20, 25, 30 and 35 MeV, and for various medium thicknesses up to 5000 μm, where R was defined as the ratio D A /D E , D A was the absorbed proton dose in cellular components, and D E was the absorbed proton dose in a separate radiation detector. In the present work, D A and D E were determined using the MCNPX (Monte Carlo N-Particle eXtended) code version 2.4.0. For lower incident proton energies (i.e., 10, 15 and 20 MeV), formation of Bragg-peak-like features were noticed in their R-vs-medium-layer-thickness relationships, and large R values of >7 and >6 were obtained for cytoplasm and nucleus of cells, respectively, which highlighted the importance of careful consideration of the medium thickness in radiobiological experiments.

The Monte Carlo method was employed to simulate realistic treatment situations for photon and proton radiation therapy for a set of Oak Ridge National Laboratory (ORNL) pediatric phantoms for 15, 10, 5 and 1-year olds as well as newborns. Complete radiotherapy situations were simulated using the previously developed NRUrad input code for Monte Carlo N-Particle (MCNP) code package. Each pediatric phantom was irradiated at five different positions, namely, the testes, colon, liver, left lung and brain, and the doses in targeted organs (Dt) were determined using the track length estimate of energy. The dispersed photon and proton doses in non-targeted organs (Dd), namely, the skeleton, skin, brain, spine, left and right lungs were computed. The conversion coefficients (F = Dd/Dt) of the dispersed doses were used to study the dose dispersion in different non-targeted organs for phantoms for 15, 10, 5 and 1-year olds as well as newborns. In general, the F values were larger for younger patients. The F values for non-targeted organs for phantoms for 1-year olds and newborns were significantly larger compared to those for other phantoms. The dispersed doses from proton radiation therapy were also found to be significantly lower than those from conventional photon radiation therapy. For example, the largest F values for the brain were 65.6% and 0.206% of the dose delivered to the left lung (P4) for newborns during photon and proton radiation therapy, respectively. The present results demonstrated that dispersion of photons and generated electrons significantly affected the absorbed doses in non-targeted organs during pediatric photon therapy, and illustrated that proton therapy could in general bring benefits for treatment of pediatric cancer patients.

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 137Cs 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.

X-ray and γ-ray photons have been widely used for studying radiobiological effects of ionizing radiations. Photons are indirectly ionizing radiations so they need to set in motion electrons (which are a directly ionizing radiation) to perform the ionizations. When the photon dose decreases to below a certain limit, the number of electrons set in motion will become so small that not all cells in an "exposed" cell population can get at least one electron hit. When some cells in a cell population are not hit by a directly ionizing radiation (in other words not irradiated), there will be rescue effect between the irradiated cells and non-irradiated cells, and the resultant radiobiological effect observed for the "exposed" cell population will be different. In the present paper, the mechanisms underlying photon interactions in radiobiological experiments were studied using our developed NRUphoton computer code, which was benchmarked against the MCNP5 code by comparing the photon dose delivered to the cell layer underneath the water medium. The following conclusions were reached: (1) The interaction fractions decreased in the following order: 16O > 12C > 14N > 1H. Bulges in the interaction fractions (versus water medium thickness) were observed, which reflected changes in the energies of the propagating photons due to traversals of different amount of water medium as well as changes in the energy-dependent photon interaction cross-sections. (2) Photoelectric interaction and incoherent scattering dominated for lower-energy (10 keV) and high-energy (100 keV and 1 MeV) incident photons. (3) The fractions of electron ejection from different nuclei were mainly governed by the photoelectric effect cross-sections, and the fractions from the 1s subshell were the largest. (4) The penetration fractions in general decreased with increasing medium thickness, and increased with increasing incident photon energy, the latter being explained by the corresponding reduction in interaction cross-sections. (5) The areas under the angular distribution curves of photons exiting the medium layer and subsequently undergoing interactions within the cell layer became smaller for larger incident photon energies. (6) The number of cells suffering at least one electron hit increased with the administered dose. For larger incident photon energies, the numbers of cells suffering at least one electron hit became smaller, which was attributed to the reduction in the photon interaction cross-section. These results highlighted the importance of the administered dose in radiobiological experiments. In particular, the threshold administered doses at which all cells in the exposed cell array suffered at least one electron hit might provide hints on explaining the intriguing observation that radiation-induced cancers can be statistically detected only above the threshold value of ~100 mSv, and thus on reconciling controversies over the linear no-threshold model.

Monte Carlo method was used to study the characteristics of neutron interactions with cells underneath a water medium layer with varying thickness. The following results were obtained. (1) The fractions of neutron interaction with 1H, 12C, 14N and 16O nuclei in the cell layer were studied. The fraction with 1H increased with increasing medium thickness, while decreased for 12C, 14N and 16O nuclei. The bulges in the interaction fractions with 12C, 14N and 16O nuclei were explained by the resonance spikes in the interaction cross-section data. The interaction fraction decreased in the order: 1H > 16O > 12C > 14N. (2) In general, as the medium thickness increased, the number of "interacting neutrons" which exited the medium and then further interacted with the cell layer increased. (3) The area under the angular distributions for "interacting neutrons" decreased with increasing incident neutron energy. Such results would be useful for deciphering the reasons behind discrepancies among existing results in the literature.

Radiotherapy is a common cancer treatment module, where a certain amount of dose will be delivered to the targeted organ. This is achieved usually by photons generated by linear accelerator units. However, radiation scattering within the patient's body and the surrounding environment will lead to dose dispersion to healthy tissues which are not targets of the primary radiation. Determination of the dispersed dose would be important for assessing the risk and biological consequences in different organs or tissues. In the present work, the concept of conversion coefficient (F) of the dispersed dose was developed, in which F = (Dd/Dt), where Dd was the dispersed dose in a non-targeted tissue and Dt is the absorbed dose in the targeted tissue. To quantify Dd and Dt, a comprehensive model was developed using the Monte Carlo N-Particle (MCNP) package to simulate the linear accelerator head, the human phantom, the treatment couch and the radiotherapy treatment room. The present work also demonstrated the feasibility and power of parallel computing through the use of the Message Passing Interface (MPI) version of MCNP5.

Monte Carlo method has been used to determine the efficiency for proton production and to study the energy and angular distributions of the generated protons. The ENDF library of cross sections is used to simulate the interactions between the neutrons and the atoms in a polyethylene (PE) layer, while the ranges of protons with different energies in PE are determined using the Stopping and Range of Ions in Matter (SRIM) computer code. The efficiency of proton production increases with the PE layer thickness. However the proton escaping from a certain polyethylene volume is highly dependent on the neutron energy and target thickness, except for a very thin PE layer. The energy and angular distributions of protons are also estimated in the present paper, showing that, for the range of energy and thickness considered, the proton flux escaping is dependent on the PE layer thickness, with the presence of an optimal thickness for a fixed primary neutron energy.