Effect of bolus materials on dose deposition in deep tissues during electron beam radiotherapy

Several materials are utilized in the production of bolus, which is essential for superficial tumor radiotherapy. This research aimed to compare the variations in dose deposition in deep tissues during electron beam radiotherapy when employing different bolus materials. Specifically, the study developed general superficial tumor models (S-T models) and postoperative breast cancer models (P-B models). Each model comprised a bolus made of water, polylactic acid (PLA), polystyrene, silica-gel or glycerol. Geant4 was employed to simulate the transportation of electron beams within the studied models, enabling the acquisition of dose distributions along the central axis of the field. A comparison was conducted to assess the dose distributions in deep tissues. In regions where the percentage depth dose (PDD) decreases rapidly, the relative doses (RDs) in the S-T models with silica-gel bolus exhibited the highest values. Subsequently, RDs for PLA, glycerol and polystyrene boluses followed in descending order. Notably, the RDs for glycerol and polystyrene boluses were consistently below 1. Within the P-B models, RDs for all four bolus materials are consistently below 1. Among them, the smallest RDs are observed with the glycerol bolus, followed by silica-gel, PLA and polystyrene bolus in ascending order. As PDDs are ~1-3% or smaller, the differences in RDs diminish rapidly until are only around 10%. For the S-T and P-B models, polystyrene and glycerol are the most suitable bolus materials, respectively. The choice of appropriate bolus materials, tailored to the specific treatment scenario, holds significant importance in safeguarding deep tissues during radiotherapy.

Dosimetry of Large Field Valencia applicators for Cobalt-60-based brachytherapy

BACKGROUND

Non-melanoma skin cancer is one of the most common types of cancer and one of the main approaches is brachytherapy. For small lesions, the treatment of this cancer with brachytherapy can be done with two commercial applicators, one of these is the Large Field Valencia Applicators (LFVA).

PURPOSE

The aim of this study is to test the capabilities of the LFVA to use clinically 60Co sources instead of the 192Ir ones. This study was designed for the same dwell positions and weights for both sources.

METHODS

The Penelope Monte Carlo code was used to evaluate dose distribution in a water phantom when a 60Co source is considered. The LFVA design and the optimized dwell weights reported for the case of 192Ir are maintained with the only exception of the dwell weight of the central position, that was increased. 2D dose distributions, field flatness, symmetry and the leakage dose distribution around the applicator were calculated.

RESULTS

When comparing the dose distributions of both sources, field flatness and symmetry remain unchanged. The only evident difference is an increase of the penumbra regions for all depths when using the 60Co source. Regarding leakage, the maximum dose within the air volume surrounding the applicator is in the order of 20% of the prescription dose for the 60Co source, but it decreases to less than 5% at about 1 cm distance.

CONCLUSIONS

Flatness and symmetry remains unaltered as compared with 192Ir sources, while an increase in leakage has been observed. This proves the feasibility of using the LFVA in a larger range of clinical applications.

Quantifying accuracy on distance estimation tasks: A Monte Carlo study

An important aspect of perceptual learning involves understanding how well individuals can perceive distances, sizes, and time-to-contact. Oftentimes, the primary goal in these experiments is to assess participants' errors (i.e., how accurately participants perform these tasks). However, the manner in which researchers have quantified error, or task accuracy, has varied. The use of different measures of task accuracy, to include error scores, ratios, and raw estimates, indicates that the interpretation of findings depends on the measure of task accuracy utilized. In an effort to better understand this issue, we used a Monte Carlo simulation to evaluate five dependent measures of accuracy: raw distance judgments, a ratio of true to estimated distance judgments, relative error, signed error, and absolute error. We simulated data consistent with prior findings in the distance perception literature and evaluated how findings and interpretations vary as a function of the measure of accuracy used. We found there to be differences in both statistical findings (e.g., overall model fit, mean square error, Type I error rate) and the interpretations of those findings. The costs and benefits of utilizing each accuracy measure for quantifying accuracy in distance estimation studies are discussed.

First Monte Carlo beam model for ultra-high dose rate radiotherapy with a compact electron LINAC

BACKGROUND

FLASH radiotherapy based on ultra-high dose rate (UHDR) is actively being studied by the radiotherapy community. Dedicated UHDR electron devices are currently a mainstay for FLASH studies.

PURPOSE

To present the first Monte Carlo (MC) electron beam model for the UHDR capable Mobetron (FLASH-IQ) as a dose calculation and treatment planning platform for preclinical research and FLASH-radiotherapy (RT) clinical trials.

METHODS

The initial beamline geometry of the Mobetron was provided by the manufacturer, with the first-principal implementation realized in the Geant4-based GAMOS MC toolkit. The geometry and electron source characteristics, such as energy spectrum and beamline parameters, were tuned to match the central-axis percentage depth dose (PDD) and lateral profiles for the pristine beam measured during machine commissioning. The thickness of the small foil in secondary scatter affected the beam model dominantly and was fine tuned to achieve the best agreement with commissioning data. Validation of the MC beam modeling was performed by comparing the calculated PDDs and profiles with EBT-XD radiochromic film measurements for various combinations of applicators and inserts.

RESULTS

The nominal 9 MeV electron FLASH beams were best represented by a Gaussian energy spectrum with mean energy of 9.9 MeV and variance (σ) of 0.2 MeV. Good agreement between the MC beam model and commissioning data were demonstrated with maximal discrepancy < 3% for PDDs and profiles. Hundred percent gamma pass rate was achieved for all PDDs and profiles with the criteria of 2 mm/3%. With the criteria of 2 mm/2%, maximum, minimum and mean gamma pass rates were (100.0%, 93.8%, 98.7%) for PDDs and (100.0%, 96.7%, 99.4%) for profiles, respectively.

CONCLUSIONS

A validated MC beam model for the UHDR capable Mobetron is presented for the first time. The MC model can be utilized for direct dose calculation or to generate beam modeling input required for treatment planning systems for FLASH-RT planning. The beam model presented in this work should facilitate translational and clinical FLASH-RT for trials conducted on the Mobetron FLASH-IQ platform.

Health risk assessment of heavy metals based on source analysis and Monte Carlo in the downstream basin of the Zishui

In this study, stone coal mines in the lower reaches of the Zijiang River were adopted as the research object. To analyze the spatial distribution, sources, and health risks of heavy metals in the surrounding soil of stone coal mines, 82 topsoil samples were collected in the study area, and the contents of 8 heavy metals including Cd, Hg, As, Cr, Pb, Cu, Ni, and Zn were determined. The spatial distribution of heavy metals was analyzed using ArcGIS, and the pollution sources of heavy metals were identified using Positive matrix factorization (PMF). Then, Monte Carlo and health risk assessment models were used to evaluate the health risks of different populations. The results showed that the average content of heavy metals followed the order of Zn > Cr > Pb > Cu > Ni > As > Cd > Hg, and the contents of all heavy metals were higher than the soil background values of Hunan Province. The high-value areas of heavy metals content were mostly concentrated in the central region close to areas with a notable concentration of stone coal mines. PMF identified four pollution sources, namely, mining activities (26.9%), atmospheric deposition (18.8%), natural sources (32.8%) and agricultural sources (21.5%). The carcinogenic and non-carcinogenic risks for children were higher than those for adults, with As and Cd posing higher carcinogenic risks to children. Based on the source of health risks, it was determined that the health risks could be primarily attributed to agricultural sources, and As was the main heavy metal causing health risks. This study provides theoretical support for treating heavy metal pollution in mining basins.

Shape-Interaction Dualism: Unraveling Complex Phase Behavior in Triangular Particle Monolayers

This paper examines the effect of finite attractive and repulsive interactions on the self-assembly of triangular-shaped particles on a triangular lattice. The ground state analysis of the lattice model has revealed an infinite sequence of ordered structures, a phenomenon referred to as the "devil's staircase" of phase transitions. The model has been studied at finite temperatures using both the transfer-matrix and tensor renormalization group methods. The concurrent use of these two methods lends credibility to the obtained results. It has been demonstrated that the initial ordered structures of the "devil's staircase" persist at non-zero temperatures. Further increase of the attraction between particles or a decrease of the temperature induces the appearance of subsequent ordered structures of the "devil's staircase". The corresponding phase diagram of the model has been calculated. The phase behavior of our model agrees qualitatively with the phase behavior of trimesic acid adsorption layer on single crystal surfaces. .

Angle-dependent light scattering in tissue phantoms for the case of thin bone layers with predominant forward scattering

The cochlea forms a key element of the human auditory system in the temporal bone. Damage to the cochlea continues to produce significant impairment for sensory reception of environmental stimuli. To improve this impairment, the optical cochlear implant forms a new research approach. A prerequisite for this method is to understand how light propagation, as well as scattering, reflection, and absorption, takes place within the cochlea. We offer a method to study the light distribution in the human cochlea through phantom materials which have the objective to mimic the optical behavior of bone and Monte-Carlo simulations. The calculation of an angular distribution after scattering requires a phase function. Often approximate functions like Henyey-Greenstein, two-term Henyey-Greenstein or Legendre polynomial decompositions are used as phase function. An alternative is to exactly calculate a Mie distribution for each scattering event. This method provides a better fit to the data measured in this work.

Beam quality correction factors for ionization chambers in a 0.35 T magnetic resonance (MR)-linac - A Monte Carlo study

PURPOSE

The purpose of this study was to directly calculate [Formula: see text] correction factors for four cylindrical ICs for a 0.35 T MR-linac using the Monte Carlo (MC) method.

METHODS

A previously-validated TOPAS/GEANT4 MC head model of the 0.35 T MR-linac was employed. The MR-compatible Exradin A12, A1SL, A26, and A28 cylindrical ICs were modeled considering the dead volume in the air cavity. The [Formula: see text] correction factor was determined for initial electron energies of 5-7 MeV. The correction factor was calculated for all four angular orientations in the lateral plane. The impact of the 0.35 T magnetic field on the IC response was also investigated.

RESULTS

The maximum beam quality dependence in the [Formula: see text] exhibited by the A12, A1SL, A26, and A28 ICs was 1.10 %, 2.17 %, 0.81 %, and 1.75 %, respectively, considering all angular orientations. The magnetic field dependence was < 1 % and the maximum [Formula: see text] correction was < 2 % when the detector was aligned along the direction of the magnetic field at 0° and 180° angles. The A12 IC over-responded up to 5.40 % for the orthogonal orientation. An asymmetry in the response of up to 8.30 % was noted for the A28 IC aligned at 90° and 270° angles.

CONCLUSIONS

A parallel orientation for the IC, with respect to the magnetic field, is recommended for reference dosimetry in MRgRT. Both over and under-response in the IC signal was noted for the orthogonal orientations, which is highly dependent on the cavity diameter, cavity length, and the dead volume.

MCGPU-PET: An Open-Source Real-Time Monte Carlo PET Simulator

Monte Carlo (MC) simulations are commonly used to model the emission, transmission, and/or detection of radiation in Positron Emission Tomography (PET). In this work, we introduce a new open-source MC software for PET simulation, MCGPU-PET, which has been designed to fully exploit the computing capabilities of modern GPUs to simulate the acquisition of more than 100 million coincidences per second from voxelized sources and material distributions. The new simulator is an extension of the PENELOPE-based MCGPU code previously used in cone-beam CT and mammography applications. We validated the accuracy of the accelerated code by comparing it to GATE and PeneloPET simulations achieving an agreement within 10 percent approximately. As an example application of the code for fast estimation of PET coincidences, a scan of the NEMA IQ phantom was simulated. A fully 3D sinogram with 6382 million true coincidences and 731 million scatter coincidences was generated in 54 seconds in one GPU. MCGPU-PET provides an estimation of true and scatter coincidences and spurious background (for positron-gamma emitters such as 124I) at a rate 3 orders of magnitude faster than CPU-based MC simulators. This significant speed-up enables the use of the code for accurate scatter and prompt-gamma background estimations within an iterative image reconstruction process.

Possible improvements in effective fill factor using X-ray fluorescent interpixel reflectors

BACKGROUND

The spatial resolution of energy-integrating diagnostic CT scanners is limited by interpixel reflectors on the detector, which optically isolate pixels but create dead space. Because the width of the reflector cannot easily be decreased, fill factor diminishes as resolution increases.

PURPOSE

We propose loading (or mixing) a high-Z element into the reflectors, causing the reflectors to be X-ray fluorescent. Re-emitted characteristic X-rays could be detected in adjacent pixels, increasing the effective fill factor and compensating for fill factor loss with higher-resolution detectors. The purpose of this work is to understand the physical principles of this approach and to analyze its effectiveness using Monte Carlo simulations.

METHODS

Detector pixels were modeled using the GEANT4 Monte Carlo package. The width of the reflector was kept constant at 0.1 mm throughout, and we considered pixel pitches between 0.5 and 1 mm. The pixelated scintillator material was gadolinium oxysulfide, 3 mm thick. The baseline reflector material was chosen to be acrylic, and varying concentrations of a high-Z element were loaded into the material. We assumed that the optical characteristics of pixels were ideal (no absorption within pixels, perfect reflection at boundaries). The detector was irradiated uniformly with 10,000 X-ray photons to estimate its spectral response. The figure of merit was the variance of the detector signal at zero frequency normalized to that of an ideal single-bin photon-counting detector with 100% fill factor. Sensitivity analyses were conducted to understand the effect of varying the high-Z element concentration and the spectrum.

RESULTS

Initial simulations suggested that a k-edge near 50 keV would be ideal. Gd was therefore selected as the high-Z material. The relative variances for a conventional energy integrating detector without Gd at 1 mm pixel pitch (81% fill factor) and 0.5 mm pixel pitch (64% fill factor) were 1.38 and 1.74, compared to 1.00 for an ideal photon counting detector, implying a 26% variance penalty for 0.5 mm pitch. When 1 g/cm3 Gd was loaded into the interpixel reflector, the relative variance improved to 1.27 and 1.43, respectively, implying that the variance penalty for including Gd together with 0.5 mm pitch is only 4%. Performance was nearly maximized at 1.0 g/cm3 of Gd, but a concentration of 0.5 g/cm3 of Gd showed most of the benefit. Improvements depend weakly on kV, with lower kV associated with higher improvements. An external anti-scatter grid was not modeled in our simulations and would reduce the expected benefit, depending greatly on the pitch and dimensionality of the anti-scatter grid.

CONCLUSIONS

The losses in fill factor associated with smaller pixel pitch can be reduced if Gd or a similar element could be loaded into the interpixel reflector. These improvements in noise efficiency are yet to be verified experimentally.

Comparison of CBCT radiation doses with conventional radiographs in TMJ imaging using Monte Carlo simulations

Three-dimensional imaging methods are widely used for evaluation of bony changes of temporomandibular joint (TMJ). Besides, lateral and posterio-anterior TMJ projections in both closed- and open-jaw positions for each temporomandibular joint are used as two-dimensional diagnostic tools. The purpose of the present study was to compare effective and mean organ absorbed doses of plain radiography techniques with those of different modalities of cone beam computed tomography (CBCT) scanning of an adult's temporomandibular joint. PCXMC 2.0 software was used to calculate mean organ and effective doses. A NewTom CBCT device (Newtom 5G XL; QR systems; Verona, Italy) was simulated at 360° rotation using a 6 × 6 cm2 FOV in standard and high-resolution modes. Lateral and posterio-anterior TMJ plain projections were simulated according to recommendations of the manufacturer of the Planmeca ProMax® 2D S3 device. Doses for both projections were simulated with Monte Carlo methods and the International Commission on radiological protection adult reference computational phantoms. The highest mean organ absorbed doses occurred in bone surfaces, salivary glands, and skull for posterio-anterior TMJ and lateral TMJ, and for CBCT scanning in all examinations. The effective doses of posterio-anterior and lateral TMJ plain radiographs were found to be higher than those of the Standard Mode-Eco Scan CBCT. Therefore, the lowest effective dose was calculated in Standard Mode-Eco Scan CBCT. It is concluded that NewTom 5G XL Standard Mode-Eco Scan CBCT can be used instead of plain radiographs (lateral and posterio-anterior TMJ) in temporomandibular joint imaging as it allows visualizing the three-dimensional structure of the temporomandibular joint as an advantage.

Beam quality and the mystery behind the lower percentage depth dose in grid radiation therapy

Grid therapy recently has been picking momentum due to favorable outcomes in bulky tumors. This is being termed as Spatially Fractionated Radiation Therapy (SFRT) and lattice therapy. SFRT can be performed with specially designed blocks made with brass or cerrobend with repeated holes or using multi-leaf collimators where dosimetry is uncertain. The dosimetric challenge in grid therapy is the mystery behind the lower percentage depth dose (PDD) in grid fields. The knowledge about the beam quality, indexed by TPR20/10 (Tissue Phantom Ratio), is also necessary for absolute dosimetry of grid fields. Since the grid may change the quality of the primary photons, a new [Formula: see text] should be evaluated for absolute dosimetry of grid fields. A Monte Carlo (MC) approach is provided to resolving the dosimetric issues. Using 6 MV beam from a linear accelerator, MC simulation was performed using MCNPX code. Additionally, a commercial grid therapy device was used to simulate the grid fields. Beam parameters were validated with MC model for output factor, depth of maximum dose, PDDs, dose profiles, and TPR20/10. The electron and photon spectra were also compared between open and grid fields. The dmax is the same for open and grid fields. The PDD with grid is lower (~ 10%) than the open field. The difference in TPR20/10 of open and grid fields is observable (~ 5%). Accordingly, TPR20/10 is still a good index for the beam quality in grid fields and consequently choose the correct [Formula: see text] in measurements. The output factors for grid fields are 0.2 lower compared to open fields. The lower depth dose with grid therapy is due to lower depth fluence with scatter radiation but it does not impact the dosimetry as the calibration parameters are insensitive to the effective beam energies. Thus, standard dosimetry in open beam based on international protocol could be used.

The Influence of B4C Film Density on Damage Threshold Based on Monte Carlo Method for X-ray Mirror

The uniformity and consistency of X-ray mirror film materials prepared by experimental methods are difficult to guarantee completely. These factors directly affect the service life of free electron laser devices in addition to its own optical properties. Therefore, the quality of the film material, especially the density, has a critical effect on its application. Boron carbide film and monocrystalline silicon substrate were suitable examples to explore their influence of density on the damage threshold based on Monte Carlo and heat-conduction methods. Through simulation results, it was found that the change in film density could affect the energy deposition depth and damage threshold. When the film density was 2.48 g/cm3, it had relatively high damage threshold in all energy ranges. And then the specific incident parameter for practical application was investigated. It was found that the damage mechanism of the B4C/Si was the melting of the interface. And the damage threshold was also higher with the film density of 2.48 g/cm3. Therefore, it was recommended to maintain the density at this value as far as possible when preparing the film, and to ensure the uniformity and consistency of the film material.

The Effect of Polymer-Solvent Interaction on the Swelling of Polymer Matrix Tablets: A Magnetic Resonance Microscopy Study Complemented by Bond Fluctuation Model Simulations

Polymer matrix tablets are an important drug-delivery system widely used for oral drug administration. Understanding the tablet hydration process, both experimentally and theoretically, is, thus, very important for the development of drug delivery systems that exhibit high drug loading capacity and controlled release potential. In this study, we used magnetic resonance microscopy (MRM) to nondestructively and dynamically analyze the water hydration process of xanthan-based tablets. The swelling process was characterized by well-resolved fronts of erosion, swelling, and penetration. The experimental results were complemented by numerical simulations of the polymer matrix hydration process. In the simulations, the polymer tablet matrix was modeled as an assembly of interacting chains with embedded drug particles, while its hydration process was mediated by interaction with solvent particles. The swelling dynamics were modeled within a Monte Carlo-based bond fluctuation model (BFM) that elegantly accounted for steric and nearest-neighbor interactions. This study provides an efficient experimental-theoretical approach for the study of polymer matrix swelling processes.

Assessing the impact of PV panel climate-based degradation rates on inverter reliability in grid-connected solar energy systems

This paper provides an evaluation of a 4-kW grid-connected full-bridge PV inverter under three different scenarios to assess its reliability with a fixed PV degradation rate, with a climate-based degradation rate, and without considering PV degradation. The climate-based degradation rates are estimated using a physics-based model that considers the different parameters influencing the PV reliability. Three different locations representing three different climate zones (hot and dry, hot and humid, and moderate climates) have been chosen in this study. The estimated lifetime of the IGBT, the switching device in the PV inverter, varies depending on the location, with the inclusion of fixed and climate-based degradation rates extending the lifespan of the PV inverter in the examined locations. The results demonstrate the significant impact of PV climate-based degradation rates on power electronics' reliability assessment and the importance of considering various factors in predicting device failures. To ensure the PV inverter's lifespan over the desired period in areas with high solar irradiation rates and extremely hot climates, the design parameters should be slightly elevated above the standard value.

Voxel-level dosimetry for combined Iodine-131 radiopharmaceutical therapy and external beam radiotherapy treatment paradigms of head and neck cancer

PURPOSE

Targeted radiopharmaceutical therapy (RPT) in combination with external beam radiotherapy (EBRT) shows promise as a method to increase tumor control and mitigate potential high-grade toxicities associated with re-treatment for patients with recurrent head and neck cancer. This work establishes a patient-specific dosimetry framework that combines Monte Carlo based dosimetry from the two radiation modalities at the voxel level using deformable image registration (DIR) and radiobiological constructs for patients enrolled in a phase I clinical trial combining EBRT and RPT.

METHODS

Serial SPECT/CT patient scans performed at approximately 24, 48, 72, and 168 hours post-injection of 577.2 MBq/m2 (15.6 mCi/m2) iodine-131 containing RPT agent called XXX 131. Clinical EBRT treatment plans were created on a treatment planning CT (TPCT) using RayStation; SPECT/CT images were deformably registered to the TPCT using the Elastix DIR module in 3D Slicer and assessed by measuring mean activity concentrations and absorbed doses. Monte Carlo EBRT dosimetry was computed using EGSnrc. RPT dosimetry was conducted using a GEANT4 based RPT dosimetry platform named XXX. Radiobiological metrics (BED, EQD2) were utilized to combine the two radiation modalities.

RESULTS

The DIR method provided good agreement for the activity concentrations and calculated absorbed dose in the tumor volumes for the SPECT/CT and TPCT; the maximum mean absorbed dose difference was -11.2%. Based on the RPT absorbed dose calculations, two to four EBRT fractions were removed from patients' EBRT treatments. From the combined treatment, the absorbed dose to target volumes ranged from 57.14 - 75.02 Gy. When including partial volume corrections (PVC), the mean EQD2 to the PTV from EBRT+RPT was -3.11% to 1.40% different compared to EBRT alone.

CONCLUSION

This work demonstrated the clinical feasibility of performing combined EBRT+RPT dosimetry on TPCTs. Dosimetry guides treatment decisions for EBRT and this work provides a bridge for the same paradigm to be implemented within the rapidly emerging clinical RPT space.

Radiation shielding assessment for interventional radiology personnel: Geant4 dosimetry of lead-free compositions

The inherent biological hazards associated with ionizing radiation necessitate the implementation of effective shielding measures, particularly in medical applications. Interventional radiology, in particular, poses a unique challenge as it often exposes medical personnel to prolonged periods of high x-ray doses. Historically, lead and lead-based compounds have been the primary materials employed for shielding against photons. However, the drawbacks of lead, including its substantial weight causing personnel's inflexibility and its toxicity, have raised concerns regarding its long-term impact on both human health and the environment. Barium tantalate has emerged as a promising alternative, due to its unique attenuation properties against low-energy x-rays, specifically targeting the weak absorption area of lead. In the present study, we employ the Geant4 Monte Carlo simulation tool to investigate various formulations of barium tantalate doped with rare earth elements. The aim is to identify the optimal composition for shielding x-rays in the context of interventional radiology. To achieve this, we employ a reference x-ray spectrum typical of interventional radiology procedures, with energies extending up to 90 keV, within a carefully designed simulation setup. Our primary performance indicator is the reduction in air kerma transmission. Furthermore, we assess the absorbed doses to critical organs at risk within a standard human body phantom protected by the shield. Our results demonstrate that specific concentrations of the examined rare earth impurities can enhance the shielding performance of barium tantalate. To mitigate x-ray exposure in interventional radiology, our analysis reveals that the most effective shielding performance is achieved when using barium tantalate compositions containing 15% Erbium or 10% Samarium by weight. These findings suggest the possibility of developing lead-free shielding solutions or apron for interventional radiology personnel, offering a remarkable reduction in weight (exceeding 30%) while maintaining shielding performance at levels comparable to traditional lead-based materials.

Application of NEMA protocols to verify GATE models based on the Digital Biograph Vision and the Biograph Vision Quadra scanners

PURPOSE

The Monte Carlo method is an effective tool to simulate and verify PET systems. Furthermore, it can help in the design and optimization of new medical imaging devices and algorithms. In this framework, the goal of this work is to verify the GATE toolkit performance when applied to simulate two Siemens Healthineers PET scanners: a standard axial field-of-view Biograph Vision scanner and the new long axial field-of-view Biograph Vision Quadra scanner.

METHODS

The simulation toolkit GATE is based on GEANT4, comprising its main functionalities and a set of domain-specific features in the field of medical physics. To accomplish our purpose, the guidelines described in the NEMA NU 2-2018 protocol are reproduced. Then the simulated results are compared to experimental data available in the literature for both PET scanners. The assessment of the models includes different studies of sensitivity, count rate performances, spatial resolution and image quality. These tests are intended to evaluate the image quality of PET devices.

RESULTS

In the spatial resolution test, relative errors lower than 8% are obtained between the experiments and GATE models. The sensitivity is 17.2 cps/kBq (Vision) and 175.9 cps/kBq (Quadra), representing relative differences with the experiment of 6% and 0.3%, respectively. Deviations from peak NECR are less than 9%. In the Image Quality test, the contrast recovery coefficient for hot spheres, with 8 iterations and 5 subsets, ranges between 57-83% for Vision and 54-86% for Quadra. These values represent a maximum deviation between the simulations and the experiments of 10% for the Quadra scanner. In the case of the Vision scanner, the highest difference is observed for the 10 mm sphere (∼38%) due to the higher contrast recovery of the experiment caused by the Gibbs artifact from the use of PSF reconstruction.

CONCLUSIONS

The results of the simulations have provided evidence of a good agreement between the experimental data and the results obtained with GATE. Thus, this work supports the capability of this MC toolkit to accurately simulate the models of the Vision and Quadra scanners. This study has laid the basis for further research in this field and has identified several areas that could be explored.

Integrating microdosimetricin vitroRBE models for particle therapy into TOPAS MC using the MicrOdosimetry-based modeliNg for RBE ASsessment (MONAS) tool

Objective.In this paper, we present MONAS (MicrOdosimetry-based modelliNg for relative biological effectiveness (RBE) ASsessment) toolkit. MONAS is a TOPAS Monte Carlo extension, that combines simulations of microdosimetric distributions with radiobiological microdosimetry-based models for predicting cell survival curves and dose-dependent RBE.Approach.MONAS expands TOPAS microdosimetric extension, by including novel specific energy scorers to calculate the single- and multi-event specific energy microdosimetric distributions at different micrometer scales. These spectra are used as physical input to three different formulations of themicrodosimetric kinetic model, and to thegeneralized stochastic microdosimetric model(GSM2), to predict dose-dependent cell survival fraction and RBE. MONAS predictions are then validated against experimental microdosimetric spectra andin vitrosurvival fraction data. To show the MONAS features, we present two different applications of the code: (i) the depth-RBE curve calculation from a passively scattered proton SOBP and monoenergetic12C-ion beam by using experimentally validated spectra as physical input, and (ii) the calculation of the 3D RBE distribution on a real head and neck patient geometry treated with protons.Main results.MONAS can estimate dose-dependent RBE and cell survival curves from experimentally validated microdosimetric spectra with four clinically relevant radiobiological models. From the radiobiological characterization of a proton SOBP and12C fields, we observe the well-known trend of increasing RBE values at the distal edge of the radiation field. The 3D RBE map calculated confirmed the trend observed in the analysis of the SOBP, with the highest RBE values found in the distal edge of the target.Significance.MONAS extension offers a comprehensive microdosimetry-based framework for assessing the biological effects of particle radiation in both research and clinical environments, pushing closer the experimental physics-based description to the biological damage assessment, contributing to bridging the gap between a microdosimetric description of the radiation field and its application in proton therapy treatment with variable RBE.

Multiscale Monte Carlo simulations for dosimetry in x-ray breast imaging: Part I - Macroscopic scales

BACKGROUND

X-ray breast imaging modalities are commonly employed for breast cancer detection, from screening programs to diagnosis. Thus, dosimetry studies are important for quality control and risk estimation since ionizing radiation is used.

PURPOSE

To perform multiscale dosimetry assessments for different breast imaging modalities and for a variety of breast sizes and compositions. The first part of our study is focused on macroscopic scales (down to millimeters).

METHODS

Nine anthropomorphic breast phantoms with a voxel resolution of 0.5 mm were computationally generated using the BreastPhantom software, representing three breast sizes with three distinct values of volume glandular fraction (VGF) for each size. Four breast imaging modalities were studied: digital mammography (DM), contrast-enhanced digital mammography (CEDM), digital breast tomosynthesis (DBT) and dedicated breast computed tomography (BCT). Additionally, the impact of tissue elemental compositions from two databases were compared. Monte Carlo (MC) simulations were performed with the MC-GPU code to obtain the 3D glandular dose distribution (GDD) for each case considered with the mean glandular dose (MGD) fixed at 4 mGy (to facilitate comparisons).

RESULTS

The GDD within the breast is more uniform for CEDM and BCT compared to DM and DBT. For large breasts and high VGF, the ratio between the minimum/maximum glandular dose to MGD is 0.12/4.02 for DM and 0.46/1.77 for BCT; the corresponding results for a small breast and low VGF are 0.35/1.98 (DM) and 0.63/1.42 (BCT). The elemental compositions of skin, adipose and glandular tissue have a considerable impact on the MGD, with variations up to 30% compared to the baseline. The inclusion of tissues other than glandular and adipose within the breast has a minor impact on MGD, with differences below 2%. Variations in the final compressed breast thickness alter the shape of the GDD, with a higher compression resulting in a more uniform GDD.

CONCLUSIONS

For a constant MGD, the GDD varies with imaging modality and breast compression. Elemental tissue compositions are an important factor for obtaining MGD values, being a source of systematic uncertainties in MC simulations and, consequently, in breast dosimetry.

How to check a simulation study

Simulation studies are powerful tools in epidemiology and biostatistics, but they can be hard to conduct successfully. Sometimes unexpected results are obtained. We offer advice on how to check a simulation study when this occurs, and how to design and conduct the study to give results that are easier to check. Simulation studies should be designed to include some settings in which answers are already known. They should be coded in stages, with data-generating mechanisms checked before simulated data are analysed. Results should be explored carefully, with scatterplots of standard error estimates against point estimates surprisingly powerful tools. Failed estimation and outlying estimates should be identified and dealt with by changing data-generating mechanisms or coding realistic hybrid analysis procedures. Finally, we give a series of ideas that have been useful to us in the past for checking unexpected results. Following our advice may help to prevent errors and to improve the quality of published simulation studies.

Radiation Shielding efficiency of lead-tungsten-boron glasses with Sb, Al, and Bi against gamma, neutron and charge particles

We report on newly developed nuclear shielding glass system based on lead-tungsten-boron (PWB) for radiation applications against photon, neutron and charge particles. This newly developed system contains also different additions, in low concentrations, such as Sb, Al and Bi. The gamma/photon shielding performance was tested by using FLUKA Monte Carlo. Moreover, the shielding efficiency of the present system is examined against charged particles (light and heavy ones) and neutrons. The highest gamma/photons attenuation is observed in the lowest incident energy and this is at the region of the photoelectric absorption. We also observe that the values of effective atomic number (Zeff) show a peak at 100 keV incident energy. The reduction of these values is higher for photon energy range 0.1-1 MeV than below 80 keV energies. The lowest half value layer (d1/2), reflecting the best shielding efficiency, is recorded for the PWB-Bi system. The PWB-Bi system demonstrates promising performance better than many of commercial and standard systems and heavy concretes.

Study of the radionuclide deposition in the radioactive ion line of the Selective Production of Exotic Species (SPES) facility

SPES (Selective Production of Exotic Species) is a second generation facility for the production of radioactive ion beams that is going to be commissioned at the Laboratori Nazionali di Legnaro of INFN at Legnaro, Padua, Italy. Radioactive neutron-rich isotopes are expected to be produced by nuclear fission induced by a 40 MeV, 200 μA primary proton beam impinging on a 238UCx target. The expected reaction rate is about 1013 fission/s. Radioactive ion beams are produced using the isotope separation on-line technique. The production of such an amount of radioactive species raises radiological issues throughout the life cycle of the facility. A study of the radioactive contamination of the components of the radioactive ion beam line is performed with the FLUKA Monte Carlo simulation code, under realistic hypotheses for the produced isobaric beams. The present results complete previous studies focused on the radiological impact of the production target irradiation, the residual activation of the primary proton beam line and the radioactive contamination of the ion source complex. The overall ambient dose equivalent rate due to the different radiation sources is calculated at several positions inside the production bunker and at different times after a typical one-year operating period of the facility with the 238UCx target at full power. The obtained results and the developed methodology provide the guidelines and the needed tools to plan ordinary and extraordinary interventions as well as final decommissioning of the SPES facility.

Retrieval of chromophore concentration changes in a digital human head model using analytical mean partial pathlengths of photons

Significance

Continuous-wave functional near-infrared spectroscopy has proved to be a valuable tool for assessing hemodynamic activity in the human brain in a non-invasively and inexpensive way. However, most of the current processing/analysis methods assume the head is a homogeneous medium, and hence do not appropriately correct for the signal coming from the scalp. This effect can be reduced by considering light propagation in a layered model of the human head, being the Monte Carlo (MC) simulations the gold standard to this end. However, this implies large computation times and demanding hardware capabilities.

Aim

In this work, we study the feasibility of replacing the homogeneous model and the MC simulations by means of analytical multilayered models, combining in this way, the speed and simplicity of implementation of the former with the robustness and accuracy of the latter.

Approach

Oxy- and deoxyhemoglobin (HbO and HbR, respectively) concentration changes were proposed in two different layers of a magnetic resonance imaging (MRI)-based meshed model of the human head, and then these changes were retrieved by means of (i) a typical homogeneous reconstruction and (ii) a theoretical layered reconstruction.

Results

Results suggest that the use of analytical models of light propagation in layered models outperforms the results obtained using traditional homogeneous reconstruction algorithms, providing much more accurate results for both, the extra- and the cerebral tissues. We also compare the analytical layered reconstruction with MC-based reconstructions, achieving similar degrees of accuracy, especially in the gray matter layer, but much faster (between 4 and 5 orders of magnitude).

Conclusions

We have successfully developed, implemented, and validated a method for retrieving chromophore concentration changes in the human brain, combining the simplicity and speed of the traditional homogeneous reconstruction algorithms with robustness and accuracy much more similar to those provided by MC simulations.

Validation of an in vivo transit dosimetry algorithm using Monte Carlo simulations and ionization chamber measurements

PURPOSE

Transit dosimetry is a safety tool based on the transit images acquired during treatment. Forward-projection transit dosimetry software, as PerFRACTION, compares the transit images acquired with an expected image calculated from the DICOM plan, the CT, and the structure set. This work aims to validate PerFRACTION expected transit dose using PRIMO Monte Carlo simulations and ionization chamber measurements, and propose a methodology based on MPPG5a report.

METHODS

The validation process was divided into three groups of tests according to MPPG5a: basic dose validation, IMRT dose validation, and heterogeneity correction validation. For the basic dose validation, the fields used were the nine fields needed to calibrate PerFRACTION and three jaws-defined. For the IMRT dose validation, seven sweeping gaps fields, the MLC transmission and 29 IMRT fields from 10 breast treatment plans were measured. For the heterogeneity validation, the transit dose of these fields was studied using three phantoms: 10 , 30 , and a 3 cm cork slab placed between 10 cm of solid water. The PerFRACTION expected doses were compared with PRIMO Monte Carlo simulation results and ionization chamber measurements.

RESULTS

Using the 10 cm solid water phantom, for the basic validation fields, the root mean square (RMS) of the difference between PerFRACTION and PRIMO simulations was 0.6%. In the IMRT fields, the RMS of the difference was 1.2%. When comparing respect ionization chamber measurements, the RMS of the difference was 1.0% both for the basic and the IMRT validation. The average passing rate with a γ(2%/2 mm, TH = 20%) criterion between PRIMO dose distribution and PerFRACTION expected dose was 96.0% ± 5.8%.

CONCLUSION

We validated PerFRACTION calculated transit dose with PRIMO Monte Carlo and ionization chamber measurements adapting the methodology of the MMPG5a report. The methodology presented can be applied to validate other forward-projection transit dosimetry software.