On the perturbation effect and LET dependence of beam quality correction factors in carbon ion beams

BACKGROUND

In a recent study, we reported beam quality correction factors, f_Q , in carbon ion beams using Monte Carlo (MC) methods for a cylindrical and a parallel-plate ionization chamber (IC). A non-negligible perturbation effect was observed; however, the magnitude of the perturbation correction due to the specific IC subcomponents was not included. Furthermore, the stopping power data presented in the International Commission on Radiation Units and Measurements (ICRU) report 73 were used, whereas the latest stopping power data have been reported in the ICRU report 90.

PURPOSE

The aim of this study was to extend our previous work by computing f_Q correction factors using the ICRU 90 stopping power data and by reporting IC-specific perturbation correction factors. Possible energy or linear energy transfer (LET) dependence of the f_Q correction factor was investigated by simulating both pristine beams and spread-out Bragg peaks (SOBPs).

METHODS

The TOol for PArticle Simulation (TOPAS)/GEANT4 MC code was used in this study. A 30 × 30 × 50 cm³ water phantom was simulated with a uniform 10 × 10 cm² parallel beam incident on the surface. A Farmer-type cylindrical IC (Exradin A12) and two parallel-plate ICs (Exradin P11 and A11) were simulated in TOPAS using the manufacturer-provided geometrical drawings. The f_Q correction factor was calculated in pristine carbon ion beams in the 150-450 MeV/u energy range at 2 cm depth and in the middle of the flat region of four SOBPs. The k_Q correction factor was calculated by simulating the f_Qo correction factor in a ⁶⁰ Co beam at 5 cm depth. The perturbation correction factors due to the presence of the individual IC subcomponents, such as the displacement effect in the air cavity, collecting electrode, chamber wall, and chamber stem, were calculated at 2 cm depth for monoenergetic beams only. Additionally, the mean dose-averaged and track-averaged LET was calculated at the depths at which the f_Q was calculated.

RESULTS

The ICRU 90 f_Q correction factors were reported. The p_dis correction factor was found to be significant for the cylindrical IC with magnitudes up to 1.70%. The individual perturbation corrections for the parallel-plate ICs were <1.0% except for the A11 p_cel correction at the lowest energy. The f_Q correction for the P11 IC exhibited an energy dependence of >1.00% and displayed differences up to 0.87% between pristine beams and SOBPs. Conversely, the f_Q for A11 and A12 displayed a minimal energy dependence of <0.50%. The energy dependence was found to manifest in the LET dependence for the P11 IC. A statistically significant LET dependence was found only for the P11 IC in pristine beams only with a magnitude of <1.10%.

CONCLUSIONS

The perturbation and k_Q correction factor should be calculated for the specific IC to be used in carbon ion beam reference dosimetry as a function of beam quality.

A Monte Carlo simulation for Moving Light Source in Intracavity PDT

We developed a simulation method for modeling the light fluence delivery in intracavity Photodynamic Therapy (icav-PDT) for pleural lung cancer using a moving light source. Due to the large surface area of the pleural lung cavity, the light source needs to be moved to deliver a uniform dose around the entire cavity. While multiple fixed detectors are used for dosimetry at a few locations, an accurate simulation of light fluence and fluence rate is still needed for the rest of the cavity. We extended an existing Monte Carlo (MC) based light propagation solver to support moving light sources by densely sampling the continuous light source trajectory and assigning the proper number of photon packages launched along the way. The performance of Simphotek GPU CUDA-based implementation of the method - PEDSy-MC - has been demonstrated on a life-size lung-shaped phantom, custom printed for testing icav-PDT navigation system at the Perlman School of Medicine (PSM) - calculations completed under a minute (for some cases) and within minutes have been achieved. We demonstrate results within a 5% error of the analytic solution for multiple detectors in the phantom. PEDSy-MC is accompanied by a dose-cavity visualization tool that allows real-time inspection of dose values of the treated cavity in 2D and 3D, which will be expanded to ongoing clinical trials at PSM. PSM has developed a technology to measure 8-detectors in a pleural cavity phantom using Photofrin-mediated PDT that has been used during validation.

Assessment of EMF Human Exposure Levels Due to Wearable Antennas at 5G Frequency Band

(1) Background: This work aims to assess human exposure to EMF due to two different wearable antennas tuned to two 5G bands. (2) Methods: The first one was centered in the lower 5G band, around f = 3.5 GHz, whereas the second one was tuned to the upper 5G band, at 26.5 GHz. Both antennas were positioned on the trunk of four simulated human models. The exposure assessment was performed by electromagnetic numerical simulations. Exposure levels were assessed by quantifying the specific absorption rate averaged on 10 g of tissue (SAR_10g) and the absorbed power density (S_ab), depending on the frequency of the wearable antenna. (3) Results: the higher exposure values that resulted were always mainly concentrated in a superficial area just below the antenna itself. In addition, these resulting distributions were narrowed around their peak values and tended to flatten toward lower values in farther anatomical body regions. All the exposure levels complied with ICNIRP guidelines when considering realistic input power. (4) Conclusions: This work highlights the importance of performing an exposure assessment when the antenna is placed on the human wearer, considering the growth of wearable technology and its wide variety of application, particularly regarding future 5G networks.

Exploring the Physical and Biological Aspects of BNCT with a Carboranylmethylbenzo[b]acridone Compound in U87 Glioblastoma Cells

Boron neutron capture therapy (BNCT) is a re-emerging technique for selectively killing tumor cells. Briefly, the mechanism can be described as follows: after the uptake of boron into cells, the thermal neutrons trigger the fission of the boron atoms, releasing the α-particles and recoiling lithium particles and high-energy photons that damage the cells. We performed a detailed study of the reactor dosimetry, cellular dose assessment, and radiobiological effects induced by BNCT in glioblastoma (GBM) cells. At maximum reactor power, neutron fluence rates were ϕ0 = 6.6 × 107 cm−2 s−1 (thermal) and θ = 2.4 × 104 cm−2 s−1 with a photon dose rate of 150 mGy·h−1. These values agreed with simulations to within 85% (thermal neutrons), 78% (epithermal neutrons), and 95% (photons), thereby validating the MCNPX model. The GEANT4 simulations, based on a realistic cell model and measured boron concentrations, showed that >95% of the dose in cells was due to the BNC reaction. Carboranylmethylbenzo[b]acridone (CMBA) is among the different proposed boron delivery agents that has shown promising properties due to its lower toxicity and important cellular uptake in U87 glioblastoma cells. In particular, the results obtained for CBMA reinforce radiobiological effects demonstrating that damage is mostly induced by the incorporated boron with negligible contribution from the culture medium and adjacent cells, evidencing extranuclear cell radiosensitivity.

Addressing transcranial electrical stimulation variability through prospective individualized dosing of electric field strength in 300 participants across two samples: the 2-SPED approach

Objective. Transcranial electrical stimulation (tES) is a promising method for modulating brain activity and excitability with variable results to date. To minimize electric (E-)field strength variability, we introduce the 2-sample prospective E-field dosing (2-SPED) approach, which uses E-field strengths induced by tES in a first population to individualize stimulation intensity in a second population.Approach. We performed E-field modeling of three common tES montages in 300 healthy younger adults. First, permutation analyses identified the sample size required to obtain a stable group average E-field in the primary motor cortex (M1), with stability being defined as the number of participants where all group-average E-field strengths ± standard deviation did not leave the population's 5-95 percentile range. Second, this stable group average was used to individualize tES intensity in a second independent population (n = 100). The impact of individualized versus fixed intensity tES on E-field strength variability was analyzed.Main results. In the first population, stable group average E-field strengths (V/m) in M1 were achieved at 74-85 participants, depending on the tES montage. Individualizing the stimulation intensity (mA) in the second population resulted in uniform M1 E-field strength (all p < 0.001) and significantly diminished peak cortical E-field strength variability (all p < 0.01), across all montages.Significance. 2-SPED is a feasible way to prospectively induce more uniform E-field strengths in a region of interest. Future studies might apply 2-SPED to investigate whether decreased E-field strength variability also results in decreased physiological and behavioral variability in response to tES.

Comparison of organ and effective dose estimations from different Monte Carlo simulation-based software methods in infant CT and comparison with direct phantom measurements

Purpose

Computational dosimetry software is routinely used to evaluate the organ and effective doses from computed tomography (CT) examinations. Studies have shown a significant variation in dose estimates between software in adult cohorts, and few studies have evaluated software for pediatric dose estimates. This study aims to compare the primary organ and effective doses estimated by four commercially available CT dosimetry software to thermoluminescent dosimeter (TLD) measurements in a 1‐year‐old phantom.

Methods

One hundred fifteen calibrated LiF (Mg, Cu, P)‐TLD 100‐H chips were embedded within an anthropomorphic phantom representing a 1‐year‐old child at positions that matched the approximate location of organs within an infant. The phantom was scanned under three protocols, each with whole‐body coverage. The mean absorbed doses from 25 radiosensitive organs and skeletal tissues were determined from the TLD readings. Effective doses for each of the protocols were subsequently calculated using ICRP 103 formalism. Dose estimates by the four Monte Carlo–based dose calculation systems were determined and compared to the directly measured doses.

Results

Most organ doses determined by computation dosimetry software aligned to phantom measurements within 20%. Additionally, comparisons between effective doses are calculated using computational and direct measurement methods aligned within 20% across the three protocols. Significant variances were found in bone surface dose estimations among dosimetry methods, likely caused by differences in bone tissue modeling.

Conclusion

All four‐dosimetry software evaluated in this study provide adequate primary organ and effective dose estimations. Users should be aware, however, of the possible estimated uncertainty associated with each of the programs.

Virtual estimation of effective dose in neutron fields

The PODIUM project aims to provide real-time assessments of occupationally exposed workers by tracking their motion and combining this with a simulation of the radiation field. The present work describes the approach that would be taken in mixed neutron-gamma fields, and details the methods for generating and applying an effective dose rate map; the required fluence to effective dose conversion coefficients at intercardinal angles are also presented. A proof-of-concept of the approach is demonstrated using a simple simulated workplace field within a calibration laboratory, with corroborative comparisons made against survey instrument measurements generally confirming good agreement. Simulated tracking of an individual within the facility was performed, recording a 1.25μSv total effective dose and accounting for dose rates as low as 0.5 nSv h^-1, which is much lower than anything that could be accurately measured by physical neutron dosemeters in such a field.

The impact of CT image parameters and skull heterogeneity modeling on the accuracy of transcranial focused ultrasound simulations

Objective. Low-intensity transcranial ultrasound stimulation (TUS) is a promising non-invasive brain stimulation (NIBS) technique. TUS can reach deeper areas and target smaller regions in the brain than other NIBS techniques, but its application in humans is hampered by the lack of a straightforward and reliable procedure to predict the induced ultrasound exposure. Here, we examined how skull modeling affects computer simulations of TUS.Approach. We characterized the ultrasonic beam after transmission through a sheep skull with a hydrophone and performed computed tomography (CT) image-based simulations of the experimental setup. To study the skull model's impact, we varied: CT acquisition parameters (tube voltage, dose, filter sharpness), image interpolation, segmentation parameters, acoustic property maps (speed-of-sound, density, attenuation), and transducer-position mismatches. We compared the impact of modeling parameter changes on model predictions and on measurement agreement. Spatial-peak intensity and location, total power, and the Gamma metric (a measure for distribution differences) were used as quantitative criteria. Modeling-based sensitivity analysis was also performed for two human head models.Main results. Sheep skull attenuation assignment and transducer positioning had the most important impact on spatial peak intensity (overestimation up to 300%, respectively 30%), followed by filter sharpness and tube voltage (up to 20%), requiring calibration of the mapping functions. Positioning and skull-heterogeneity-structure strongly affected the intensity distribution (gamma tolerances exceeded in>80%, respectively>150%, of the focus-volume in water), necessitating image-based personalized modeling. Simulation results in human models consistently demonstrate a high sensitivity to the skull-heterogeneity model, attenuation tuning, and transducer shifts, the magnitude of which depends on the underlying skull structure complexity.Significance. Our study reveals the importance of properly modeling the skull-heterogeneity and its structure and of accurately reproducing the transducer position. The results raise red flags when translating modeling approaches among clinical sites without proper standardization and/or recalibration of the imaging and modeling parameters.

USE OF ANTHROPOMORPHIC HETEROGENEOUS PHYSICAL PHANTOMS FOR VALIDATION OF COMPUTATIONAL DOSIMETRY OF MEDICAL PERSONNEL AND PATIENTS

In the dosimetry of ionizing radiation, the phantoms of the human body, which are used as a replacement for thehuman body in physical measurements and calculations, play an important, but sometimes underestimated, role.There are physical phantoms used directly for measurements, and mathematical phantoms for computationaldosimetry. Their complexity varies from simple geometry applied for calibration purposes up to very complex, whichsimulates in detail the shapes of organs and tissues of the human body. The use of physical anthropomorphic phantoms makes it possible to effectively optimize radiation doses by adjusting the parameters of CT-scanning (computed tomography) in accordance with the characteristics of the patient without compromising image quality. The useof phantoms is an indispensable approach to estimate the actual doses to the organs or to determine the effectivedose of workers - values that are regulated, but cannot be directly measured.The article contains an overview of types, designs and the fields of application of anthropomorphic heterogeneousphysical phantoms of a human with special emphasis on their use for validation of models and methods of computational dosimetry.