Microdosimetric Analysis for Boron Neutron Capture Therapy via Monte Carlo Track Structure Simulation with Modified Lithium Cross-sections.

Boron neutron capture therapy (BNCT) is a cellular-level hadron therapy achieving therapeutic effects via the synergistic action of multiple particles, including Lithium, alpha, proton, and photon. However, evaluating the relative biological effectiveness (RBE) in BNCT remains challenging. In this research, we performed a microdosimetric calculation for BNCT using the Monte Carlo track structure (MCTS) simulation toolkit, TOPAS-nBio. This paper reports the first attempt to derive the ionization cross-sections of low-energy (>0.025 MeV/u) Lithium for MCTS simulation based on the effective charge cross-section scalation method and phenomenological double-parameter modification. The fitting parameters λ1=1.101,λ2=3.486 were determined to reproduce the range and stopping power data from the ICRU report 73. Besides, the lineal energy spectra of charged particles in BNCT were calculated, and the influence of sensitive volume (SV) size was discussed. Condensed history simulation obtained similar results with MCTS when using Micron-SV while overestimating the lineal energy when using Nano-SV. Furthermore, we found that the microscopic boron distribution can significantly affect the lineal energy for Lithium, while the effect for alpha is minimal. Similar results to the published data by PHITS simulation were observed for the compound particles and monoenergetic protons when using micron-SV. Spectra with nano-SV reflected that the different track densities and absorbed doses in the nucleus together result in the dramatic difference in the macroscopic biological response of BPA and BSH. This work and the developed methodology could impact the research fields in BNCT where understanding radiation effects is crucial, such as the treatment planning system, source evaluation, and new boron drug development.

An integrated Monte Carlo track-structure simulation framework for modeling inter and intra-track effects on homogenous chemistry.

Objective. The TOPAS-nBio Monte Carlo track structure simulation code, a wrapper of Geant4-DNA, was extended for its use in pulsed and longtime homogeneous chemistry simulations using the Gillespie algorithm.Approach. Three different tests were used to assess the reliability of the implementation and its ability to accurately reproduce published experimental results: (1) a simple model with a known analytical solution, (2) the temporal evolution of chemical yields during the homogeneous chemistry stage, and (3) radiolysis simulations conducted in pure water with dissolved oxygen at concentrations ranging from 10µM to 1 mM with [H2O2] yields calculated for 100 MeV protons at conventional and FLASH dose rates of 0.286 Gy s-1and 500 Gy s-1, respectively. Simulated chemical yield results were compared closely with data calculated using the Kinetiscope software which also employs the Gillespie algorithm.Main results. Validation results in the third test agreed with experimental data of similar dose rates and oxygen concentrations within one standard deviation, with a maximum of 1% difference for both conventional and FLASH dose rates. In conclusion, the new implementation of TOPAS-nBio for the homogeneous long time chemistry simulation was capable of recreating the chemical evolution of the reactive intermediates that follow water radiolysis.Significance. Thus, TOPAS-nBio provides a reliable all-in-one chemistry simulation of the physical, physico-chemical, non-homogeneous, and homogeneous chemistry and could be of use for the study of FLASH dose rate effects on radiation chemistry.