Water Phantom Research Articles

A quantification of the electron return effect using Monte Carlo simulations and experimental measurements for the MRI-linac.

The integration of magnetic resonance (MR) imaging and linear accelerators into hybrid treatment systems has made MR-guided radiation therapy a clinical reality. This work aims to evaluate the influence of the Electron Return Effect (ERE) on the dose distributions.
This study was conducted using MRIdian (ViewRay, Cleveland, Ohio) system. Monte-Carlo simulations (MCs) and experimental measurements with EBT3 Gafchromic films were performed to investigate the dose distribution in a slab water phantom with and without a 2-cm air gap. Plus, MCs took into account different field sizes and a lung gap. A gamma analysis compared calculated versus measured dose distributions.
The MCs have shown an increase of the ERE with the radiation field size both in Percent Depth Dose (PDD) and crossline direction. As concerns to the PDD direction, the smallest field for which there was a significant dose accumulation was 4.15x4.15 cm2 both for air-gap (13.5%) and lung-gap (3.3%). The largest field for which there was a significant dose accumulation was 24.07x24.07 cm2 both for air-gap (39.7%) and lung-gap (4.9%). Instead for the crossline direction, the smallest field for which there was a significant dose accumulation was 2.49x2.49 cm2 both for air-gap (8.6% ) and lung-gap (0.5%). The largest field for which there was a significant dose accumulation was 24.07x24.07 cm2 both for air-gap (46.2%) and lung-gap (4.2%).
PDD and crossline profiles showed good agreement with a gamma-passing rate higher than 91.15% for 2%/2mm. The ERE can be adequately calculated by MC dose calculation platform available in the MRIdian Treatment Planning System. The MCs show an increase of the ERE directly proportional with the radiation field size. Good agreement was observed between the experimental measurements and calculated dose distributions.

Improved heterogeneity handling in the collapsed cone dose engine for brachytherapy.

Model-based dose calculation algorithms (MBDCA), such as the Advanced Collapsed cone Engine (ACE) in Oncentra Brachy® can be used to overcome the limitations of the TG-43 formalism. ACE is a point kernel superposition algorithm that calculates the total dose separated into primary, first-scatter, and multiple-scatter dose. Albeit ACE yields accurate results under most circumstances, several studies have reported underestimations of the dose to cortical bone. These underestimations are likely caused by approximations in the handling of multiple-scatter dose for non-water media. Such would result in noticeable deviations where the multiple-scatter is a considerable fraction of the total dose, that is, at greater distances from the source. To improve and test the accuracy of the multiple-scatter dose component in the ACE algorithm to remedy its inaccuracy for non-water geometries. A careful analysis of the transport and absorption of the multiple-scatter energy fluence revealed an inconsistency in the scaling of energy absorption ratios for non-water media of the multiple-scatter kernel. We implemented an updated algorithm version, ACEcorr, and tested it for three different geometries. All had a single 192Ir-source at the center of a cubic water phantom with a box-shaped heterogeneity of either cortical bone or air, positioned at different distances from the source. Dose distributions for the three cases were calculated with ACE and ACEcorr and compared to Monte Carlo simulations, using the percentage dose difference ratio as figure-of-merit. All dose calculation methods scored separately the dose deposited by primary, first-scattered, and multiple-scattered photons. The accuracy of the updated algorithm ACEcorr was superior to ACE. In the cortical bone heterogeneity, the mean percentage dose difference ratio for the total dose improved from to (in the worst case) by our update. Less impact was seen in the air heterogeneity, where both ACE and ACEcorr deviated less than 2% from the Monte Carlo results. The algorithm update mainly concerns the multiple-scattered dose component, but an accompanying data processing update also had a small effect ( 0.5% difference) on the primary and first-scattered dose. The calculation times were not affected. The updates to ACE improved the accuracy of multiple-scatter dose calculation for non-water media, without increasing calculation times.

Study on the impact of bowtie filter on photon-counting CT imaging

Objective.The aim of this study was to investigate the impact of the bowtie filter on the image quality of the photon-counting detector (PCD) based CT imaging.Approach.Numerical simulations were conducted to investigate the impact of bowtie filters on image uniformity using two water phantoms, with tube potentials ranging from 60 to 140 kVp with a step of 5 kVp. Subsequently, benchtop PCD-CT imaging experiments were performed to verify the observations from the numerical simulations. Additionally, various correction methods were validated through these experiments.Main results.It was found that the use of a bowtie filter significantly alters the uniformity of PCD-CT images, depending on the size of the object and the x-ray spectrum. Two notable effects were observed: the capping effect and the flattening effect. Furthermore, it was demonstrated that the conventional beam hardening correction method could effectively mitigate such non-uniformity in PCD-CT images, provided that dedicated calibration parameters were used.Significance.It was demonstrated that the incorporation of a bowtie filter results in varied image artifacts in PCD-CT imaging under different conditions. Certain image correction methods can effectively mitigate and reduce these artifacts, thereby enhancing the overall quality of PCD-CT images.

Geant4 Monte Carlo simulation of human exposure to indoor 222Rn from building materials.

The present study aimed to develop a Monte Carlo model to estimate the annual effective dose due to radon exposure sourced by radon gas in the walls and floor of a standard model room. With the purpose of developing a tool for radon level assessment in dwellings and workplaces, Geant4 toolkit was used to simulate the energy deposited by gamma rays emitted by radioactive radon progeny in a water phantom positioned at three different locations within the model room. The energy deposition was then used to estimate the annual effective dose through a deterministic approach. The simulation outcomes showed good agreement with experimental data, with the ratio between the simulated and the experimental data displaying the overestimation by a factor of approximately 1.09. Both simulation and experimental data fell within the same range, with a relative deviation of 7.7%. Additionally, the influence of various parameters, such as receptor position in the room, wall, and floor thicknesses, wall cover, and building material bulk density, on the annual effective dose due to radon inhalation in the room was evaluated. Geant4 Monte Carlo toolkit proved to be a reliable tool for radon modeling in real exposure situations.

The dosimetric accuracy of a commercial model-based dose calculation algorithm in modeling a six-groove direction modulated brachytherapy tandem applicator

Objective.With advancements in high-dose rate brachytherapy, the clinical translation of intensity modulated brachytherapy (IMBT) innovations necessitates utilization of model-based dose calculation algorithms (MBDCA) for accurate and rapid dose calculations. This study uniquely benchmarks a commercial MBDCA, BrachyVision ACUROSTM(BVA), against Monte Carlo (MC) simulations, evaluating dose distributions for a novel IMBT applicator, termed as thesix-grooveDirection Modulated Brachytherapy (DMBT) tandem, expanding beyond previous focus on partially shielded vaginal cylinder applicators, through a novel methodology.Approach.The DMBT tandem applicator, made of a tungsten alloy with six evenly spaced grooves, was simulated using the GEANT4 MC code. Subsequently, two main scenarios were created using the BVA and reproduced by the MC simulations: 'Source at the Center of the Water Phantom (SACWP)' and 'Source at the Middle of the Applicator (SAMA)' for three cubical virtual water phantoms (20 cm)3, (30 cm)3, and (40 cm)3. A track length estimator was utilized for dose calculation and 2D/3D scoring were performed. The difference in isodose surfaces/lines (i.e. coverage) at each voxel,ΔDIsodose Levels/Lines, was thus calculated for relevant normalization points (rref).Results.The coverage was comparable, based on 2D scoring, between the BVA and MC isodose surfaces/lines for the region of clinical relevance, (i.e. within 5 cm radius from the source) withΔDIsodose Lines(rref: 1 cm from the source) falling within 2% for the two scenarios for all phantom sizes. For the phantom (20 cm)3,ΔDIsodose Levels(3D scoring) recorded the range [-3.0% +6.5%] ([-7.4% +7.3%]) for 95% of the voxels of the same scoring volume for the SACWP (SAMA) scenario.Significance.The results indicated that the BVA could render comparable coverage as compared to the MC simulations in the region of clinical relevance for different phantom sizes.ΔDIsodose Linesmay offer an advantageous metric for evaluation of MBDCAs in clinical setting.

Computed tomography employing sensing of high energy particle current

Objective. We demonstrate High Energy Current Computed Tomography (HEC-CT) employing megavoltage linac x-rays.Approach. Using deterministic radiation transport we simulate two-parameter HEC-CT projections and using inverse Fourier transform we reconstruct two distinct material parameters for water phantoms with ICRP tissue inserts and materials of different atomic number Z. The HEC-CT projections are obtained by beam scanning and rotating the object.Main Results. The first HEC-CT material parameter is alike the standard attenuation coefficient with dependence on atomic number and material density similar to the Hounsfield Units. The second material parameter has opposite trends and does not find any analogy in the standard CT framework.Significance. New CT method has been invented for medical imaging or non-destructive testing. The key feature of the technique is a two-value CT reconstruction based on particle current instead of transmitted dose.

A Fast 3D Range-Modulator Delivery Approach: Validation of the FLUKA Model on a Varian ProBeam System Including a Robustness Analysis

A 3D range-modulator (RM), optimized for a single energy and a specific target shape, is a promising and viable solution for the ultra-fast dose delivery in particle therapy. The aim of this work was to investigate the impact of potential beam and modulator misalignments on the dose distribution. Moreover, the FLUKA Monte Carlo model, capable of simulating 3D RMs, was adjusted and validated for the 250 MeV single-energy proton irradiation from a Varian ProBeam system. A 3D RM was designed for a cube target shape rotated 45° around two axes using a Varian-internal research version of the Eclipse treatment planning software, and the resulting dose distribution was simulated in a water phantom. Deviations from the ideal alignment were introduced, and the dose distributions from the modified simulations were compared to the original unmodified one. Finally, the FLUKA model and the workflow were validated with base-line data measurements and dose measurements of the manufactured modulator prototype at the HollandPTC facility in Delft. The adjusted FLUKA model, optimized particularly in the scope of a single-energy FLASH irradiation with a PMMA pre-absorber, demonstrated very good agreement with the measured dose distribution resulting from the 3D RM. Dose deviations resulting from modulator-beam axis misalignments depend on the specific 3D RM and its shape, pin aspect ratio, rotation angle, rotation point, etc. A minor modulator shift was found to be more relevant for the distal dose distribution than for the spread-out Bragg Peak (SOBP) homogeneity. On the other hand, a modulator tilt (rotation away from the beam axis) substantially affected not only the depth dose profile, transforming a flat SOBP into a broad, Gaussian-like distribution with increasing rotation angle, but also shifted the lateral dose distribution considerably. This work strives to increase awareness and highlight potential pitfalls as the 3D RM method progresses from a purely research concept to pre-clinical studies and human trials. Ensuring that gantry rotation and the combined weight of RM, PMMA, and aperture do not introduce alignment issues is critical. Given all the other range and positioning uncertainties, etc., not related to the modulator, the RM must be aligned with an accuracy below 1° in order to preserve a clinically acceptable total uncertainty budget. Careful consideration of critical parameters like the pin aspect ratio and possibly a novel robust modulator geometry optimization are potential additional strategies to mitigate the impact of positioning on the resulting dose. Finally, even the rotated cube 3D modulator with high aspect ratio pin structures (~80 mm height to 3 mm pin base width) was found to be relatively robust against a slight misalignment of 0.5° rotation or a 1.5 mm shift in one dimension perpendicular to the beam axis. Given a reliable positioning and QA concept, the additional uncertainties introduced by the 3D RM can be successfully managed adopting the concept into the clinical routine.

Deep learning architecture for scatter estimation in cone-beam computed tomography head imaging with varying field-of-measurement settings.

X-ray scatter causes considerable degradation in the cone-beam computed tomography (CBCT) image quality. To estimate the scatter, deep learning-based methods have been demonstrated to be effective. Modern CBCT systems can scan a wide range of field-of-measurement (FOM) sizes. Variations in the size of FOM can cause a major shift in the scatter-to-primary ratio in CBCT. However, the scatter estimation performance of deep learning networks has not been extensively evaluated under varying FOMs. Therefore, we train the state-of-the-art scatter estimation neural networks for varying FOMs and develop a method to utilize FOM size information to improve performance. We used FOM size information as additional features by converting it into two channels and then concatenating it to the encoder of the networks. We compared our approach for a U-Net, Spline-Net, and DSE-Net, by training them with and without the FOM information. We utilized a Monte Carlo-simulated dataset to train the networks on 18 FOM sizes and test on 30 unseen FOM sizes. In addition, we evaluated the models on the water phantoms and real clinical CBCT scans. The simulation study demonstrates that our method reduced average mean-absolute-percentage-error for U-Net by 38%, Spline-Net by 40%, and DSE-net by 33% for the scatter estimation in the 2D projection domain. Furthermore, the root-mean-square error on the 3D reconstructed volumes was improved for U-Net by 43%, Spline-Net by 30%, and DSE-Net by 23%. Furthermore, our method improved contrast and image quality on real datasets such as water phantom and clinical data. Providing additional information about FOM size improves the robustness of the neural networks for scatter estimation. Our approach is not limited to utilizing only FOM size information; more variables such as tube voltage, scanning geometry, and patient size can be added to improve the robustness of a single network.

Characterization of novel 3D-printed metal shielding for brachytherapy applicators.

To characterize 3D-printed stainless steel metal samples in the presence of an Iridium-192 source for organ-at-risk sparing in gynecologic brachytherapy. Samples of 3D-printed stainless steel (5.5×5.5 cm2, thickness range 1-5mm) were embedded in a solid water phantom at varying distances from source catheters. An Ir-192 brachytherapy source was passed through the phantom and the dose was measured using EBT3 Gafchromic film. The film was initially positioned in the sagittal plane 2cm away from the catheters, with the metal directly below and then 1cm from the film. A uniform dose was delivered at the film plane. A second setup measured a depth dose curve in solid water with film in the transverse plane directly above the metal samples. This setup was recreated using Monte Carlo simulations (EGSnrc egs_brachy). Validation between methods was performed with unshielded (solid water only) measurements. The planar dose passing through the metal samples (thickness 1-5mm) at the midpoint between the film and catheters, decreased compared to solid water by 7.4%±6.9% to 26.5%±5.5%. Dose enhancement on the order of 5% was noted when metal was directly adjacent to the film. The average decrease in depth dose from a single dwell position ranged from 10.0%±5.9% (1mm) to 21.1%±5.3% (5mm) as measured with film, and from 3.8%±0.9% (1mm) to 16.3%±0.9% (5mm) using MC simulation. The average depth dose values were measured using a line width of 2.5mm for film, and 3mm for MC simulation, and the measurements generally agree within standard error. The 3D-printed metal samples show potential for personalized applicators. Maximum dose reduction of 26.5%±5.5% compared to solid water was measured at 2cm from the source using the 5mm sample. An outer layer of solid water could potentially be used to reduce dose enhancement due to increased scatter near the metal.

Effect of neutron beam properties on dose distributions in a water phantom for boron neutron capture therapy.

From the viewpoints of the advantage depths (ADs), peak tumor dose and skin dose, we evaluated the effect on the dose distribution of neutron beam properties, namely the ratio between thermal and epithermal neutron fluxes (thermal/epithermal ratio), fast neutron component and γ-ray component. Several parameter surveys were conducted with respect to the beam properties of neutron sources for boron neutron capture therapy assuming boronophenylalanine as the boron agent using our dose calculation tool, called SiDE. The ADs decreased by 3% at a thermal/epithermal ratio of 20-30% compared with the current recommendation of 5%. The skin dose increased with the increasing thermal/epithermal ratio, reaching a restricted value of 14 Gyeq at a thermal/epithermal ratio of 48%. The fast neutron component was modified using two different models, namely the 'linear model', in which the fast neutron intensity decreases log-linearly with the increasing neutron energy, and the 'moderator thickness (MT) model', in which the fast neutron component is varied by adjusting the MT in a virtual beam shaping assembly. Although a higher fast neutron component indicated a higher skin dose, the increment was <10% at a fast neutron component of <1 × 10-12Gycm2 for both models. Furthermore, in the MT model, the epithermal neutron intensity at a fast neutron component of 6.8 × 10-13Gycm2 was 41% higher compared with that of 2 × 10-13Gycm2. The γ-ray component also caused no significant disadvantages up to several times larger compared with the current recommendation.

Dosimetric evaluation of Rhizophora spp. particleboard bonded with animal protein adhesive at X-ray energies below 100 kVp

This research developed a bio-based adhesive (AP) derived from industrial slaughterhouse waste, comprising over 85% protein. The adhesive was characterized by a melting point of 193.14 °C, a neutral pH of 7, and a viscosity comparable to common wood adhesives such as urea-formaldehyde and phenol-formaldehyde. Utilizing this adhesive, a Rhizophora spp. particleboard phantom was produced, featuring wood particles of ≤149 μm, an adhesive concentration of 12%, and a target density of 1 g/cm³, adhering to a standard phantom dimension of 30 cm × 30 cm × 30 cm. The dosimetric properties of this particleboard phantom were subsequently compared with those of water and Perspex phantoms within an X-ray energy range of 60–100 kVp, employing high-sensitivity thermo-luminescent dosimeters (TL–100H). The findings indicated that the percentage depth dose (PDD) values of the AP-Rhizophora spp. particleboard were closely aligned with those of the Perspex and water phantoms, with the greatest discrepancy observed at 60 kVp. Additionally, the half-value layer (HVL) of the particleboard was similar to those of Perspex and water, particularly at diagnostic X-ray energies. These results demonstrate that the AP adhesive is effective for creating Rhizophora spp. particleboard phantoms, exhibiting dosimetric properties comparable to tissue-equivalent materials.

Proof-of-concept of real-time electromagnetic guidance for gynecologic interstitial catheters in high dose rate brachytherapy

Background and PurposeThe addition of interstitial needles to intracavitary gynecologic (GYN) high dose rate (HDR) brachytherapy has been shown to improve target coverage and organs-at-risk (OAR) sparing. However, no commercial solution allows real-time guidance of interstitial catheter placement. This phantom study aimed to evaluate the feasibility of an electromagnetic (EM) tracking system guidance workflow for GYN HDR brachytherapy treatment in a magnetic resonance imaging (MRI) and real-time transrectal ultrasound (TRUS) fusion scenario. Materials and MethodsA clinical investigational system combining a treatment planning system and the EM tracking technology was used. The 3D T2 weighted magnetic resonance (MR) image set of a patient treated with intracavitary and interstitial HDR brachytherapy was retrospectively chosen. The MR image set was used to delineate the target and the OARs. A preplan was generated to determine needles positions in advance. The implant was reproduced in a water phantom. A 3D TRUS scan was acquired, and a rigid registration between the MR and the TRUS images was performed. ResultsThe accuracy of the EM tracking system was < 1 mm for both the sagittal and the transverse modes of the TRUS probe. Contours that were delineated on the MRI were propagated on the TRUS images after the rigid registration. Needle insertion was successfully guided in real time with the EM tracking system on the TRUS live image using the MRI contours for guidance. ConclusionBased on this proof-of-concept, real-time EM-guidance of interstitial needle for GYN HDR brachytherapy appears to be feasible.

Comprehensive analysis of peripheral dose in electron beam therapy with a Varian TrueBeam Linac

ObjectiveTo investigate the characteristics of peripheral doses outside electron-beam applicators in Varian TrueBeam linacs. MethodPeripheral doses outside the electron applicator were measured for 6-, 9- and 12-MeV beams at the maximum dose depth (Dmax) for each energy source and at a source-to-surface distance of 100 cm. Measurements were performed using EBT3 films in solid water phantoms. The impact of field size on the penumbra width and peripheral doses was studied using various cutouts, including 3 cm × 3 cm, 6 cm × 6 cm, and 10 cm × 10 cm in a 10 cm × 10 cm applicator with the gantry and collimator at 0°. The influence of the applicator size was investigated using a circular cutout of 5 cm in diameter for various applicator sizes, including 6 cm × 6 cm, 10 cm × 10 cm, 15 cm × 15 cm, 20 × 20 cm, and 25 cm × 25 cm, at Dmax for each energy, while keeping the gantry and collimator angle at 0°. The measured dose profiles were compared with the Eclipse treatment planning system (TPS) predicted dose profiles. The effect of varying gantry angles (0°, 90°, and 270°) for a 3 cm × 3 cm cutout in a 10 cm × 10 cm applicator for each energy source and varying collimator angles (0°, 90°, and 270°) for a 10 cm × 10 cm field were investigated to determine their effects on the penumbra widths and peripheral doses. ResultsBoth the penumbra width and peripheral dose values increased with energy across different field sizes, gantry angles, collimator angles, and applicator sizes. Root Mean Square Deviation (RMSD) analysis indicated minimal differences between the measured profiles and TPS data. Peripheral doses remained below 5% of the maximum dose approximately 10–15 mm away from the field edges, suggesting the potential for implementing additional shielding where required. ConclusionsThis study highlights the importance of considering peripheral doses in electron radiotherapy. It is important to note the impact on healthy tissues beyond the treatment area to ensure patient safety and prevent the long-term side effects of treatment. These findings emphasize the necessity of implementing appropriate measures to minimize peripheral doses.

Monte Carlo-based dosimetry of proposed bi-radionuclide (125I and 106Ru/106Rh) eye plaque: A feasibility study.

Combining the sharp dose fall off feature of beta-emitting 106Ru/106Rh radionuclide with larger penetration depth feature of photon-emitting125I radionuclide in a bi-radionuclide plaque, prescribed dose to the tumor apex can be delivered while maintaining the tumor dose uniformity and sparing the organs at risk. The potential advantages of bi-radionuclide plaque could be of interest in context of ocular brachytherapy. The aim of the study is to evaluate the dosimetric advantages of a proposed bi-radionuclide plaque for two different designs, consisting of indigenous 125I seeds and 106Ru/106Rh plaque, using Monte Carlo technique. The study also explores the influence of other commercial 125I seed models and presence or absence of silastic/acrylic seed carrier on the calculated dose distributions. The study further included the calculation of depth dose distributions for the bi-radionuclide eye plaque for which experimental data are available. The proposed bi-radionuclide plaque consists of a 1.2-mm-thick silver (Ag) spherical shell with radius of curvature of 12.5mm, 20µm-thick-106Ru/106Rh encapsulated between 0.2mm Ag disk, and a 0.1-mm-thick Ag window, and water-equivalent gel containing 12 symmetrically arranged 125I seeds. Two bi-radionuclide plaque models investigated in the present study are designated as Design I and Design II. In Design I, 125I seeds are placed on the top of the plaque, while in Design II 106Ru/106Rh source is positioned on the top of the plaque. In Monte Carlo calculations, the plaque is positioned in a spherical water phantom of 30cm diameter. The proposed bi-radionuclide eye plaque demonstrated superior dose distributions as compared to 125I or 106Ru plaque for tumor thicknesses ranges from 5 to 10mm. Amongst the designs, dose at a given voxel for Design I is higher as compared to the corresponding voxel dose for Design II. This difference is attributed to the higher degree of attenuation of 125I photons in Ag as compared to beta particles. Influence of different 125I seed models on the normalized lateral dose profiles of Design I (in the absence of carrier) is negligible and within 5% on the central axis depth dose distribution as compared to the corresponding values of the plaque that has indigenous 125I seeds. In the presence of a silastic/acrylic seed carrier, the normalized central axis dose distributions of Design I are smaller by 3%-12% as compared to the corresponding values in the absence of a seed carrier. For the published bi-radionuclide plaque model, good agreement is observed between the Monte Carlo-calculated and published measured depth dose distributions for clinically relevant depths. Regardless of the type of 125I seed model utilized and whether silastic/acrylic seed carrier is present or not, Design I bi-radionuclide plaque offers superior dose distributions in terms of tumor dose uniformity, rapid dose fall off and lesser dose to nearby critical organs at risk over the Design II plaque. This shows that Design I bi-radionuclide plaque could be a promising alternative to 125I plaque for treatment of tumor sizes in the range 5 to 10mm.

Verification of mean dose in a tumor using 3D-printed tumor-model scintillators in anthropomorphic phantoms and a spherical phantom

To enhance treatment outcomes of radiotherapy, accurate verification of dosimetric parameters such as the mean, minimum, and maximum doses is essential. This study aimed to measure the mean absorbed dose in a tumor-shaped target and compare the results with the calculated value of the treatment planning system. Tumor model scintillators (TMSs) mimicking brain tumors were fabricated with a commercial 3D printer based on stereolithography files of tumors. The TMSs were irradiated following treatment plans made with the Leksell Gamma Plan® (Elekta AB, Stockholm, Sweden), and the mean dose in each TMS was measured three times and compared with the calculated value. Eleven TMSs were set at the center of a spherical solid water phantom. Two additional TMSs were fabricated according to images of vestibular schwannomas and were set at the tumor location in head model phantoms following each patient's head shape. The mean relative differences between the measured and calculated mean doses were less than one percent. This method provides an independent parameter for evaluating the accuracy of the treatment planning system.

Dosimetry evaluation of photon beam profile characteristics for different treatment parameters of quality assurance.

To compare beam profiles of MatriXX scanning system and water phantom for different treatment parameters. The cross-sectional study was conducted at Al-Amal National Hospital for Cancer Treatment, Baghdad, Iraq, from November 2020 to March 2021. Beam data for 6MV and 10MV photon beams generated from the linear accelerator was utilised at field sizes 20×20cm2, 15×15 cm2, 10×10cm2 and 5×5cm2 at depth 10 and source-to-skin distance 100cm. Data was obtained for both water phantom and MatriXX system. The dose distribution for the two systems were compared. Data was analysed using SPSS 24. The 32 measures taken were all related to symmetry and flatness. Flatness data indicated that all measurements were within tolerance except for cross line plane variations in 10x10cm2 field size with 6MV energy (-3.81%) and 5x5cm2 field size with 10MV energy (-3.01).Symmetry data revealed all measurement differences were within tolerance. MatriXX system could also be used for routine photon profile measurements as a substitute for water phantom.

Validation of GAMOS Monte Carlo simulation for 131Cs brachytherapy source in water and dosimetric characterization of solid water, bone, and brain tissues

PurposeThe treatment of brain tumours using permanently implanted 131Cs has become more frequent since GammaTile brachytherapy platform was cleared by the U.S. Food and Drug Administration (FDA). Current treatment planning systems (TPS) (BrachyVision v11.0.47, Varian Medical Systems, Palo Alto, CA, USA) used clinically do not account for anatomical heterogeneities. GAMOS Monte Carlo simulation was validated for the dose deposited by 131Cs in water and used to accurately determine the dose deposited in simulated human tissues, including bone and brain tissue, as well as in Solid Water (SW) phantom material. MethodGAMOS was validated and benchmarked for dose deposited by 131Cs (CS-1 Rev2) in water. The evaluation performed included air-kerma strength, radial dose function, anisotropy function, and dose rate constant based on the recommendation of TG-43U1S2 protocol. Dosimetric simulation of the 131Cs brachytherapy source was also performed in SW, bone, and brain tissue using GAMOS. Conversion factors were calculated by the ratio of the dose rate constant from water to SW, bone, and brain tissue, respectively. ResultResults for dose in water were benchmarked with MCNP simulation and experiential data published in the TG-43U1S2. The dose rate constant in water was 1.039 ± 0.017 cGy h−1 U−1 which is consistent with the results used to determine the TG-43U1S2 consensus data. The conversion factor of water to SW, water to the bone, and water to brain tissue were 0.936, 1.15, and 0.964, respectively. ConclusionOther Monte Carlo (MC) simulation codes, such as MCNP, have been previously used in the literature to determine the dose deposited by 131Cs. In this work, doses deposited by 131Cs radioactive sources were successfully simulated using GAMOS MC simulation. To our knowledge, this is the first time that GAMOS was used to simulate the dose deposited by 131Cs seed. MC simulation of dose for brachytherapy sources is especially important to accurately estimate the dose deposited in high Z materials such as skull bone, which may be widely different than the dose calculated by the TG-43U1S2 formalism.

Estimation of Dose Perturbation Due to the Presence of Metal Hip Prostheses in Radiotherapy with Electron and Photon Beams: A Monte Carlo Study

Purpose: The impact of various hip prosthesis materials on the amount of dose perturbation generated by 9 MeV electron and 18 MV photon beams during pelvic radiation therapy is analyzed in this research. Materials and Methods: The Varian 2100 C/D LINAC head for the electron (9 MeV) and photon (18 MV) modes and a water phantom with a realistic hip prosthesis were modeled using the Monte Carlo (MC) code; MCNPX (Ver. 2.6.0) and were benchmarked for measurement. Four different materials, including Cobalt–Chromium–Molybdenum-alloy (CCM), Stainless Steel (SS), Titanium, and Titanium-alloy (Ti-alloy) were evaluated. The changes in electron and photon fluences and dose perturbations due to the presence of the hip prostheses were investigated. Results: An increased dose of 13.29%, 13.77%, 6.16%, and 5.93% for 9 MeV and 30.43%, 33.05%, 10.89%, and 11.27% for 18 MV was calculated for Ti-alloy, Ti, CCM, and SS prosthesis, respectively. At 0.5 mm distance from the prosthesis, the Electron Backscatter Factor (EBF) of 1.31, 1.33, 1.15, and 1.14 for 9 MeV and 1.32, 1.34, 1.11, and 1.12 for 18 MV was calculated for Ti-alloy, Ti, CCM, and SS prosthesis, respectively. Dose perturbation is higher at a near distance from the prosthesis; by reducing the distance from the Ti-alloy prosthesis (1.5 to 0.5 mm), an increased dose of 7.70% and 5.54% resulted in 18 MV and 9 MeV, respectively. The dose decreases of up to 21% behind the hip prosthesis were calculated for 18 MV. Conclusion: It is essential to consider the dose perturbation due to the presence of a hip prosthesis to achieve the optimal treatment planning in radiotherapy.

Long-term beam output stability of an accelerator-based boron neutron capture therapy system.

Accelerator-based boron neutron capture therapy (AB-BNCT) systems are becoming commercially available and are expected to be widely used in hospitals. To ensure the safety of BNCT, establishing a quality assurance (QA) program and properly managing the stability of the system are necessary. In particular, a high level of beam output stability is required to avoid accidents because beam output is a major factor in patient dose. However, no studies have analyzed the long-term beam output stability of AB-BNCT systems. This study aimed to retrospectively analyze the long-term stability of the beam output by statistical process control (SPC) based on the QA results over 3 years. The data analyzed are the results of daily QA (DQA) and weekly QA (WQA) in an AB-BNCT system and were taken between June 2020 and September 2023. The evaluation of the stability of the beam output was based on the reaction rate between gold and neutrons calculated using the activation foil method using a gold foil. In DQA, which can be performed in a short time, the gold foil was applied directly to the beam irradiation aperture in air. In WQA, measurements were performed at the phantom surface, 2-cm depth, and 6-cm depth using a dedicated water phantom. The acquired data were retrospectively analyzed by individuals and a moving range chart (I-MR chart), exponentially weighted moving average control chart (EWMA chart), and several process capability indexes (PCIs). Over 99% of the DQA I-MR chart results were within control limits, whereas the WQA I-MR chart results showed that 1.8%, 4.1%, and 2.0% of the measurements exceeded the control limits at the surface, 2-cm depth, and 6-cm depth, respectively. The variation in the reaction rate of the gold foil before and after the replacement of the target was<0.5%. The EWMA chart results revealed no significant beam output drift for either DQA or WQA. Most measured data were normal based on the results of the Anderson-Darling test and met the requirements for PCI evaluation; most PCI values were>1.0; however, the Cpmk of DQA and the 2- and 6-cm depth WQAs between August 2021 and November 2022 in treatment course 2 were 0.83, 0.77, and 0.87, respectively, which were<1.0. The long-term stability of beam output was confirmed using SPC in an AB-BNCT system. The results of the control chart revealed no significant variation or drift in the beam output, and the quantitative evaluation using PCI revealed high stability. A routine QA program will enable us to provide safe BNCT.

Monte Carlo calculated absorbed-dose energy dependence of EBT3 and EBT4 films for 5-200 MeV electrons and 100 keV-15 MeV photons.

To use Monte Carlo simulations to study the absorbed-dose energy dependence of GAFChromic EBT3 and EBT4 films for 5-200MeV electron beams and 100keV-15MeV photon beams considering two film compositions: a previous EBT3 composition (Bekerat etal.) and the final composition of EBT3/current composition of EBT4 (Palmer etal.). A water phantom was simulated with films at 5-50mm depth in 5mm intervals. The water phantom was irradiated with flat, monoenergetic 5-200MeV electron beams and 100 and 150keV kilovoltage and 1-15MeV megavoltage photon beams and the dose to the active layer of the films was scored. Simulations were rerun with the films defined as water to compare the absorbed-dose response of film to water, . For electrons, the Bekerat etal. composition had variations in of up to from 5to200MeV. Similarly, the Palmer etal. composition had differences in up to from 5to200MeV. For photons, varied up to and from 100keV to 15MeV for the Bekerat etal. and Palmer etal. compositions, respectively. The depth of films did not appear to significantly affect for photons at any energy and for electrons at energies 50MeV. However, for 5 and 10MeV electrons, decreases of up to in were seen due to stacked films and increased beam attenuation in films compared to water. The up to and variations in for electrons and photons, respectively, across the energies considered in this study indicate the importance of calibrating films with the energy intended for measurement. Additionally, this work emphasizes potential issues with stacking films to measure depth dose curves, particularly for electron beams with energies 10MeV.