Phase Space Files Research Articles

Electron and Photon Energy Spectra Outside of 6 MV X-ray Small Radiotherapy Field Edges Produced by a Varian iX Linac

Due to the increase in the survival probability for patients treated with modern radiotherapy techniques to live enough for experimenting the late radiation effect, low dose outside the treatment volume becomes a concern. However, besides the absorbed dose, the beam quality outside the field edge should be taken into account. This work aimed at investigating the photon and electron fluence spectra outside the field edges for several small radiotherapy fields for determining the quality of the beams in order to better evaluate the secondary effect after modern radiotherapy treatments. Phase-space files of a 6 MV X-ray beam produced by a Varian iX linac for eight small fields of 0.7 × 0.7 cm2, 0.9 × 0.9 cm2, 1.8 × 1.8 cm2, 2.2 × 2.2 cm2, 2.7 × 2.7 cm2, 3.1 × 3.1 cm2, 3.6 × 3.6 cm2, and 4.5 × 4.5 cm2 and for the reference 10 × 10 cm2 field at SSD = 100 cm were generated using the BEAMnrc code. The photon and electron fluences in each field were calculated at 0.15, 1.35, and 9.85 cm water depth and several off-axis distances using FLURZnrc. The number of low-energy electrons between 1 and 10 keV at 2 cm outside the field edge increases by 60% compared to the central axis. Due to the relatively high linear energy transfer (LET) of these electrons, the results of this work should help to better evaluate the possible late effect of secondary radiation on healthy organs close to the tumor volume after radiotherapy treatment. We also observed high-energy electrons outside the field edge that are attributed to the leakage of the primary electron beam from the head of the linac. From a standpoint of radiological protection, these electrons should be taken into account when evaluating the dose delivered to the patient’s skin.

PENELOPE simulations and experiment for 6 MV clinac iX accelerator for standard and small static fields

The goal of this work was to produce accurate data for use as a ‘gold standard’ and a valid tool for measurements in reference dosimetry for standard/small static field sizes from 0.5 × 0.5 to 10 × 10 cm2. It is based on the accuracy of the phase space files (PSFs) as a key quantity. Because the IAEA general public database provides few PSFs for the Varian iX, we simulated the head through Monte Carlo (MC) simulations and calculated validated PSFs for 12 square field sizes including seven for small static fields. The resulting dosimetric calculations allowed us to reach a good level of agreement in comparison to our relative and absolute dose measurements performed on a Varian iX in water phantom. Measured and MC calculated output factors were investigated for different detectors. Based on the TRS 483 formalism and MC (PENELOPE/penEasy), we calculated output correction factors for the unshielded Diode-E (T60017) and the PinPoint-3D (T31016) micro-chamber according to manufacturers' blueprints. Our MC results were in agreement with the recommended data; they compete with recent measurements and MC simulations and in particular the TRS 483 MC data obtained from similar simulations. Moreover, our MC results provide supplemental data in comparison to TRS 483 data in particular for the PinPoint-3D (T31016). We suggest our MC output correction factors as new datasets for future TRS compilations. The work was substantial, used different robust MC strategies depending on the scoring regions, and led in most cases to uncertainties of less than 1%.

The relation between XR-QA2 and RT-QA2 GafchromicTM film optical density and absorbed dose in water produced by radionuclides

Purpose. In this study, Monte Carlo (MC) simulations were done to relate the dose-response of the film to that in water. The effect of backscattering materials (PMMA, lead, polystyrene, and air) was investigated on its influence on film density for radionuclides including Am-241, Tc-99m, I-131, Cs-137. Methods. A BEAMnrc MC simulation was designed to score a phase-space file (PSF) below the container of the radionuclide under consideration to use as an input file for the subsequent DOSXYZnrc MC simulation. The geometry of the container holding the radionuclide was built using the component modules available in BEAMnrc. BEAMDP was used to investigate the container effect on the radionuclide spectrum as well as the fluence. The DOSXYZnrc simulation produced the absorbed dose in XR-QA2 and RT-QA2 GafchromicTM films. The DOSXYZnrc simulations were repeated for the GafchromicTM film now replaced with water to get the absorbed dose in water. From these results, conversion factors for the dose in water to the film dose for the different radionuclides, Am-241, Tc-99m, I-131, and Cs-137 were obtained. The actual film dose was calculated using the specific gamma exposure constant (Γ) at a distance of 50 cm for a point source approximation. From the BEAMnrc simulations, the particle fluence was extracted from PSFs to correct for the fluence at 0.1 cm below the sources from the fluence 50 cm away since the inverse square law will not apply to finite-size sources. The absorbed dose profiles in the film were compared to the absorbed dose profiles from the MC simulations. Results. A fitting function based on the neutron depletion model fits the optical density versus absorbed film dose data well and can be used as a calibration tool to obtain the film dose from its optical density. Lead as a backscatter material results in a higher optical density change but a lower absorbed dose. The XR-QA2 GafchromicTM film is more sensitive than the RT-QA2 GafchromicTM film, showing a more responsive optical density (OD) change in the energy range of radionuclides used in this study. Conversion factors were determined to convert the dose in water to the dose in GafchromicTM film. The Am-241 and I-131 simulated absorbed dose in the film to dose in water does not fluctuate as much as the simulated absorbed dose in film and water when using Tc-99m and Cs-137. Validation was shown for the comparison of the film and MC simulation absorbed dose profiles. Conclusions. MC BEAMnrc simulations are useful to simulate radionuclides and their containers. BEAMDP extracted energy spectra showed that the radionuclide containers produced a Compton effect on the energy spectra and added filtration on the lower spectral photon components. Extracted fluence ratios from PSFs were used to calculate the absorbed dose value at 0.1 cm distance from the source. By using the fit function, the dose in the film can be determined for known optical density values. The effect of the backscatter materials showed that the XR-QA2 GafchromicTM film results in higher optical density values than the RT-QA2 GafchromicTM film. The absorbed dose in both the films is comparable but not for a radionuclide such as Am-241 with an activity of 74MBq. The lead backscatter material showed to be the most prominent in optical density enhancement, and the air equivalent material was the least prominent. The XR-QA2 GafchromicTM film is the most sensitive and will be the best option if working with low energies. The absorbed dose in the XR-QA2 GafchromicTM film also showed a good comparison to the absorbed dose in water for the Am-241 radionuclide with an activity of 74MBq. The absorbed dose in the films compares well to the MC simulated doses.

Single pencil beam benchmark of a module for Monte Carlo simulation of proton transport in the PENELOPE code.

PENH is a recently coded module for simulation of proton transport in conjunction with the Monte Carlo code PENELOPE. PENELOPE applies class II simulation to all type of interactions, in particular, to elastic collisions. PENH uses calculated differential cross sections for proton elastic collisions that include electron screening effects as well as nuclear structure effects. Proton-induced nuclear reactions are simulated from information in the ENDF-6 database or from alternative nuclear databases in ENDF format. The purpose of this work is to benchmark this module by simulating absorbed dose distributions from a single finite spot size proton pencil beam in water. Monte Carlo simulations with PENH are compared with simulation results from TOPAS Monte Carlo (v3.1p2) and RayStation Monte Carlo (v6). Different beam models are examined in terms of mean energy and energy spread to match the measured profiles. The phase-space file is derived from experimental measurements. Simulated absorbed dose distributions are compared to experimental data obtained with the ionization chamber array MatriXX 2D detector (IBA Dosimetry) in a water tank. The experiments were conducted with a clinical IBA pencil beam scanning dedicated nozzle. In all simulations a Fermi-Eyges phase-space representation of a single finite spot size proton pencil beam is used. In general, there is a good agreement between simulated results and experimental data up to a distance of 3cm from the central axis. In the core region (region where the dose is more than 10% of the maximum dose) PENH shows, overall, the smallest deviations from experimental data, with the largest radial rms (root mean square) smaller than 0.2. The results achieved by TOPAS and RayStation in that region are very close to those of PENH. For the halo region, that is the area of the dose distribution outside the core region reaching down to 0.01% of the maximum intensity, the largest rms achieved by TOPAS is always smaller than 0.5, yielding better results than the rest of the codes. The physics modeling of the PENELOPE/PENH code yields results consistent with measurements in the dose range relevant for proton therapy. The discrepancies between PENH appearing at distances larger than 3cm from the central-beam axis are accountable to the lack of neutron simulation in this code. In contradistinction, TOPAS has a better agreement with experimental data at large distances from the central-beam axis because of the simulation of neutrons.

GEANT4 simulation study of over-response phenomenon of fiber x-ray sensor* *Project supported by the Natural Science Foundation of Heilongjiang Province, China (Grant No. ZD2019H003), the Joint Research Fund in Astronomy under Cooperative Agreement Between the National Natural Science Foundation of China and Chinese Academy of Sciences (Grant Nos. U1631239 and U1931206), the 111 Project, China (Grant No. B13015), and the

The purpose of this article is to explore the cause of the over-response phenomenon of fiber x-ray sensor. The sensor is based on a length of PMMA fiber, whose end is filled with the scintillation material Gd2O2S:Tb. The Monte Carlo simulation software GEANT4 uses the phase space file provided by the International Atomic Energy Agency (IAEA), by irradiating the fiber x-ray sensor in the water phantom, counting the fluorescence signal of the optical fiber x-ray sensor after propagation through the fiber. In addition, the number of Cerenkov photons propagating through the fiber is also counted. Comparing this article with previous research, we believe that one of the reasons for the over-response of the fiber x-ray sensor is the non-linear response of the deposition energy of the scintillator to the fluorescence. By establishing a region of interest and counting the x-rays in this region, the simulation results show that the counted number of x-rays that may affect the fiber x-ray sensor is the biggest in the area of interest at a water depth of 5 cm. This result is close to the maximum dose point of the experimental and simulated percentage depth dose (PDD) curve of fiber x-ray sensor. Therefore, the second reason of the over-response phenomenon is believed to be fact that the inorganic materials such as Gd2O2S:Tb have larger effective atomic numbers, so the fiber x-ray sensors will cause more collisions with x-ray in a low energy region of 0.1 MeV–1.5 MeV.

Initial Steps Towards a Clinical FLASH Radiotherapy System: Pediatric Whole Brain Irradiation with 40 MeV Electrons at FLASH Dose Rates.

In this work, we investigated the delivery of a clinically acceptable pediatric whole brain radiotherapy plan at FLASH dose rates using two lateral opposing 40-MeV electron beams produced by a practically realizable linear accelerator system. The EGSnrc Monte Carlo software modules, BEAMnrc and DOSXYZnrc, were used to generate whole brain radiotherapy plans for a pediatric patient using two lateral opposing 40-MeV electron beams. Electron beam phase space files were simulated using a model of a diverging beam with a diameter of 10 cm at 50 cm SAD (defined at brain midline). The electron beams were collimated using a 10-cm-thick block composed of 5 cm of aluminum oxide and 5 cm of tungsten. For comparison, a 6-MV photon plan was calculated with the Varian AAA algorithm. Electron beam parameters were based on a novel linear accelerator designed for the PHASER system and powered by a commercial 6-MW klystron. Calculations of the linear accelerator's performance indicated an average beam current of at least 6.25 µA, providing a dose rate of 115 Gy/s at isocenter, high enough for cognition-sparing FLASH effects. The electron plan was less homogenous with a homogeneity index of 0.133 compared to the photon plan's index of 0.087. Overall, the dosimetric characteristics of the 40-MeV electron plan were suitable for treatment. In conclusion, Monte Carlo simulations performed in this work indicate that two lateral opposing 40-MeV electron beams can be used for pediatric whole brain irradiation at FLASH dose rates of >115 Gy/s and serve as motivation for a practical clinical FLASH radiotherapy system, which can be implemented in the near future.

Monitor unit calculation in electron therapy using Monte Carlo Simulation: a GUI for the phase-space field trimming

Purpose: We developed a graphical user interface (GUI) for electron phase-space field trimming using Monte Carlo simulations. This GUI can be used for monitor unit (MU) calculation in electron therapy. Methods: The GUI and electron field trimming algorithm were developed using MATLAB and C code. Phase-space files for the electron fields were generated using the EGSnrc code based on a Varian 21EX Linac with variables of applicator size, field size and energy. Verification of the alogrithm was carried out by comparing the relative output factor, which was used for MU calculation, predicted by Monte Carlo simulations and from actual measurements. Results: Our electron field trimming algorithm was found to be about five times faster than the original Monte Carlo simulation. Clinically, the GUI performed best when using voxel size ≥ 0.3 × 0.3 × 0.3 cm3, and field size larger than 2 cm in radius based on an acceptable deviation of 2%. Conclusion: A GUI for generating irregular field for MU calculation using Monte Carlo simulations was created as a user-friendly tool in electron therapy.

Dosimetry modeling of focused kV x‐ray radiotherapy for wet age‐related macular degeneration

Wet (neovascular) age-related macular degeneration (AMD) is the leading cause of blindness in the United States. The mainstay treatment requires monthly intravitreal injection of anti-vascular endothelial growth factor (anti-VEGF) drugs, associated with multiple visits, high cost, and the risk of procedural injury and infection. Anti-VEGF drugs inhibit the formation of neovasculature but do not directly attack it. Radiotherapy can destroy neovasculature and potentially also inhibit wet-AMD associated inflammation and fibrosis not addressed by VEGF inhibitors. However, the current collimation-based radiotherapy device uses fixed 4mm beams, which are prone to overtreat or undertreat the choroidal neovascularization (CNV) lesions because of their various sizes and shapes. This simulation study evaluates personalized conformal treatment with focused kV radiation using cutting-edge polycapillary x-ray optics. Simulation of the polycapillary optics was achieved via Monte Carlo (MC)-based three-dimensional (3D) geometric ray tracing. Phase-space files modeling the focused photons were generated. The method was previously verified by phantom measurements. The ultrasmall ~0.2mm beam focal spot perpendicular to the beam direction enables spatially fractionated grid therapy, which has been shown to preferentially damage abnormal neovascular blood vessels vs normal ones. Geant4-based MC simulations of scanning while rotating beam delivery were performed to conformally treat three clinical cases of large, medium, and small CNV lesions with regular and grid deliveries. Dose delivery uncertainties due to positioning errors were analyzed, including ±0.75mm displacement in the three orthogonal directions and ±5° vertical/horizontal rotation of the eyeball. The simulated CNV treatments by 60-kVp focused x-ray beams show highly conformal delivery of dose to the lesion plus margin (0.75mm) with sharp dose fall-offs and controllable spatial modulation patterns. The 90%-10% isodose penumbra is <0.5mm. With a prescription dose of 16Gy to the lesions, the critical structure doses are well below the tolerance. The average CNV dose varies within 10% (mostly within 4%) due to 0.75-mm linear displacements and 5-degree gaze angle rotation of the eyeball. Focused kV technique allows personalized treatment of CNV lesions and reduces unwanted radiation to adjacent healthy tissue. The simulated dose distribution is superior to currently available techniques.

Validation of histogram‐based virtual source models for different IGRT kV‐imaging systems

Image-guided radiotherapy (IGRT) improves tumor control but its intensive use may entrain late side effects caused by the additional imaging doses. There is a need to better quantify the additional imaging doses, so they can be integrated in the therapeutic workflow. Currently, no dedicated software enables to compute patient-specific imaging doses on a wide range of systems and protocols. As a first step toward this objective, we propose a common methodology to model four different kV-imaging systems used in radiotherapy (Varian's OBI, Elekta's XVI, Brainlab's ExacTrac, and Accuray's Cyberknife) using a new type of virtual source model based on Monte Carlo calculations. We first describe our method to build a simplified description of the photon output, or virtual source models (VSMs), of each imaging system. Instead of being constructed using measurement data, as it is most commonly the case, our VSM is used as the summary of the phase-space files (PSFs) resulting from a first Monte Carlo simulation of the considered x-ray tube. Second, the VSM is used as a photon generator for a second MC simulation in which we compute the dose. Then, the proposed VSM is thoroughly validated against standard MC simulation using PSFs on the XVI system. Last, each modeled system is compared to profiles and depth-dose-curve measurements performed in homogeneous phantom. Comparisons between PSF-based and VSM-based calculations highlight that VSMs could provide equivalent dose results (within 1% of difference) than PSFs inside the imaging field-of-view (FOV). In contrast, VSMs tend to underestimate (for up to 20%) calculated doses outside of the imaging FOV due to the assumptions underlying the VSM construction. In addition, we showed that the use of VSMs allows reducing calculation time by at least a factor of 2.8. Indeed, for identical simulation times, statistical uncertainties on dose distributions computed using VSMs were much lower than those obtained from PSF-based calculations. For each of the four imaging systems, VSMs were successfully validated against measurements in homogeneous phantoms, and are therefore ready to be used for future preclinical studies in heterogeneous or anthropomorphic phantoms. The cross system modeling methodology developed here should enable, later on, to estimate precisely and accurately patient-specific 3D dose maps delivered during a large range of kV-imaging procedures.

A method for generating intensity-modulated radiation therapy fields for small animal irradiators utilizing 3D-printed compensator molds.

The purpose of this study was to investigate the feasibility of using fused deposition modeling (FDM) three-dimensional (3D) printer to generate radiation compensators for high-resolution (~1mm) intensity-modulated radiation therapy (IMRT) for small animal radiation treatment. We propose a novel method incorporating 3D-printed compensator molds filled with NaI powder. The inverse planning module of the computational environment for radiotherapy research (CERR) software was adapted to simulate the XRAD-225Cx irradiator, both geometry and kV beam quality (the latter using a phase space file provided for XRAD-225Cx). A nine-field IMRT treatment was created for a scaled-down version of the imaging and radiation oncology core (IROC) Head and Neck IMRT credentialing test, recreated on a 2.2-cm-diameter cylindrical phantom. Optimized fluence maps comprising nine fields and a total of 2564 beamlets were calculated at resolution of 1.25×1.25mm2 . A hollow compensator mold was created (using in-house software and algorithm) for each field using 3D printing with polylactic acid (PLA) filaments. The molds were then packed with sodium iodide powder (NaI, measured density ρNaI =2.062g/cm3 ). The mounted compensator mold thickness was limited to 13.8mm due to clearance issues with couch collision. At treatment delivery, each compensator was manually mounted to a customized block tray attached to the reference 40×40mm2 collimator. Compensator reproducibility among three repeated 3D-printed molds was measured with Radiochromic EBT2 film. The two-dimensional (2D) dose distributions of the nine fields were compared to calculated 2D doses from CERR using gamma comparisons with distance-to-agreement criteria of 0.5-0.25mm and dose difference criteria of 3-5%. Good reproducibility of 3D-printed compensator manufacture was observed with mean error of ±0.024Gy and relative dose error of ±4.2% within the modulated part of the beam. Within the limit of 13.8mm compensator height, a maximum radiation blocking efficiency of 91.5% was achieved. Per field, about 45.5g of NaI powder was used. Gamma analysis on each of the nine delivered IMRT fields using radiochromic films resulted in eight of nine treatment fields with >90% pass rate with 5%/0.5mm tolerances. However, low gamma passing rate of 49-66% (3%/0.25mm to 5%/0.5mm) was noted in one field, attributed to fabrication errors resulting in over-filling the mold. The nine-field treatment plan was delivered in under 30min with no mechanical or collisional issues. We show the feasibility of high spatial resolution IMRT treatment on a small animal irradiator utilizing 3D-printed compensator shells packed with NaI powder. Using the PLA mold with NaI powder was attractive due to the ease of 3D printing a PLA mold at high geometric resolution and the well-balanced attenuation properties of NaI powders that prevented the mold from becoming too bulky. IMRT fields with 1.25-mm resolution are capable with significant fluence modulation with relative dose accuracy of ±4.2%.

Validation of a Lévy electron energy straggling model for an Elekta Synergy® linear accelerator

AimIn this study, an EGSnrc based Monte Carlo electron model was validated for an Elekta Synergy® 160-leaf Agility™ linear accelerator. A previously reported electron energy straggling model based on a Lévy distribution was tested against water tank measurements and a specially designed heterogeneous multi-layered phantom. This included PDD, beam profile, and relative output factor (ROF) comparison. All data passed a 2%/2 mm gamma criterion with the exception of some ROF data, which showed discrepancies of up to 2.7%. MethodsBEAMnrc was used to accurately model the linac that included the improved exit electron energy spectrum based on a Lévy distribution. The resulting BEAMnrc phase space files were used as sources in DOSXYZnrc for water tank dose distribution simulations consisting of 6 electron beam energies, 11 field sizes, and source-to-surface distances (SSDs) of 95 and 100 cm. Evaluation parameters included PDD, dose profiles, and relative output factors, as well as phantom PDD and dose profile measurements with EBT3 gafchromic film. ResultsThe improved exit electron beam energy spectrum caused simulated data to comply with measured data (PPD's and dose profiles) with a 100% pass rate using a 2%/2 mm criterion except for some relative output factors that deviated by 2.7% from measured ones in water. This was observed for both 95 and 100 cm SSD data. Good agreement was obtained between film and simulation data within 2% in more than 90% of PDD and profile measurements. ConclusionsThe Lévy based energy straggling model for electron beams allowed for accurate electron beam characterization in water tank and phantom measurements.

Technical note: development and validation of a Monte Carlo tool for analysis of patient-generated photon scatter

Scattered radiation unavoidably generated in the patient will negatively impact both kilovoltage (KV) and megavoltage (MV) imaging applications. Recently, ‘hybrid’ methods (i.e. combining analytical and Monte Carlo (MC) techniques) are being investigated as a solution to accurately yet quickly calculate the scattered contribution for both KV and MV images. We have developed a customized MC simulation user code for investigating the individual components of patient-scattered photon fluence, which serves as a valuable tool in this area of research. The MC tool is based on the EGSnrc/DOSXYZnrc user code. The IAUSFL flag options associated with subroutine AUSGAB, combined with LATCH tracking, are used to classify the various interactions of particles with the media. Photons are grouped into six different categories: primary, 1st Compton scatter, 1st Rayleigh scatter, multiple scatter, bremsstrahlung, and positron annihilation. We take advantage of the geometric boundary check in DOSXYZnrc, to write exiting photon particle information to a phase-space file. The tool is validated using homogeneous and heterogeneous phantom configurations with monoenergetic and polyenergetic beams under parallel and divergent beam geometry, comparing MC-simulated exit primary fluence and singly-scattered fluence to corresponding analytical calculations. This MC tool has been validated to separately score the primary and scatter fluence components of the KV and MV imaging applications in the field of radiation therapy. The results are acceptable for the various configurations and beam energies tested here. Overall, the mean percentage differences are less than 0.2% and standard deviations less than 1.6%. This will be a critical test instrument for research in photon scatter applications and particularly for the development of hybrid methods, and is freely available from the authors for research purposes.

Small field output correction factors of the microSilicon detector and a deeper understanding of their origin by quantifying perturbation factors.

PurposeThe aim of this study is the experimental and Monte Carlo‐based determination of small field correction factors for the unshielded silicon detector microSilicon for a standard linear accelerator as well as the Cyberknife System. In addition, a detailed Monte Carlo analysis has been performed by modifying the detector models stepwise to study the influences of the detector’s components.MethodsSmall field output correction factors have been determined for the new unshielded silicon diode detector, microSilicon (type 60023, PTW Freiburg, Germany) as well as for the predecessors Diode E (type 60017, PTW Freiburg, Germany) and Diode SRS (type 60018, PTW Freiburg, Germany) for a Varian TrueBeam linear accelerator at 6 MV and a Cyberknife system. For the experimental determination, an Exradin W1 scintillation detector (Standard Imaging, Middleton, USA) has been used as reference. The Monte Carlo simulations have been performed with EGSnrc and phase space files from IAEA as well as detector models according to manufacturer blueprints. To investigate the influence of the detector’s components, the detector models have been modified stepwise.ResultsThe correction factors for the smallest field size investigated at the TrueBeam linear accelerator (equivalent dosimetric square field side length Sclin = 6.3 mm) are 0.983 and 0.939 for the microSilicon and Diode E, respectively. At the Cyberknife system, the correction factors of the microSilicon are 0.967 at the smallest 5‐mm collimator compared to 0.928 for the Diode SRS. Monte Carlo simulations show comparable results from the measurements and literature.ConclusionThe microSilicon (type 60023) detector requires less correction than its predecessors, Diode E (type 60017) and Diode SRS (type 60018). The detector housing has been demonstrated to cause the largest perturbation, mainly due to the enhanced density of the epoxy encapsulation surrounding the silicon chip. This density has been rendered more water equivalent in case of the microSilicon detector to minimize the associated perturbation. The sensitive volume itself has been shown not to cause observable field size‐dependent perturbation except for the volume‐averaging effect, where the slightly larger diameter of the sensitive volume of the microSilicon (1.5 mm) is still small at the smallest field size investigated with corrections <2%. The new microSilicon fulfils the 5% correction limit recommended by the TRS 483 for output factor measurements at all conditions investigated in this work.

Evaluation of dosimetric parameters of small fields of 6 MV flattening filter free photon beam measured using various detectors against Monte Carlo simulation

AbstractPurpose:This study aims to evaluate dosimetric parameters like percentage depth dose, dosimetric field size, depth of maximum dose surface dose, penumbra and output factors measured using IBA CC01 pinpoint chamber, IBA stereotactic field diode (SFD), PTW microDiamond against Monte Carlo (MC) simulation for 6 MV flattening filter-free small fields.Materials and Methods:The linear accelerator used in the study was a Varian TrueBeam® STx. All field sizes were defined by jaws. The required shift to effective point of measurement was given for CC01, SFD and microdiamond for depth dose measurements. The output factor of a given field size was taken as the ratio of meter readings normalised to 10 × 10 cm2 reference field size without applying any correction to account for changes in detector response. MC simulation was performed using PRIMO (PENELOPE-based program). The phase space files for MC simulation were adopted from the MyVarian Website.Results and Discussion:Variations were seen between the detectors and MC, especially for fields smaller than 2 × 2 cm2 where the lateral charge particle equilibrium was not satisfied. Diamond detector was seen as most suitable for all measurements above 1 × 1 cm2. SFD was seen very close to MC results except for under-response in output factor measurements. CC01 was observed to be suitable for field sizes above 2 × 2 cm2. Volume averaging effect for penumbra measurements in CC01 was observed. No detector was found suitable for surface dose measurement as surface ionisation was different from surface dose due to the effect of perturbation of fluence. Some discrepancies in measurements and MC values were observed which may suggest effects of source occlusion, shift in focal point or mismatch between real accelerator geometry and simulation geometry.Conclusion:For output factor measurement, TRS483 suggested correction factor needs to be applied to account for the difference in detector response. CC01 can be used for field sizes above 2 × 2 cm2 and microdiamond detector is suitable for above 1 × 1 cm2. Below these field sizes, perturbation corrections and volume averaging corrections need to be applied.

Investigating the impact of collimator size variation on the single beam radiation of Gamma Knife PerfexionTM based on Monte Carlo simulation

The Gamma Knife PerfexionTM (Elekta Oncology Systems) is a specialized stereotactic radiotherapy device that uses 192 Cobalt-60 radioactive sources. A dynamic collimator system enables three circular field sizes of 4 mm, 8 mm, and 16 mm diameter to be delivered. This study investigates the dosimetric characteristics of the three different sized collimated beams using Monte Carlo (MC) simulations. The MC simulation made use of previously calculated phase space from a Monte-Carlo model of the encapsulated Cobalt-60 sources. The collimators were modelled using the FLATFILT component module of the BEAMnrc user code. The output of the single beam radiation was stored in a circular scoring plane with a radius of 1.4 cm at source surface distance (SSD) 17.20 cm in the phase space file. The simulation parameters were used 2.1 x 109 particles, electron cut off energy (ECUT) of 0.7 MeV, and photon cut off energy (PCUT) of 0.01 MeV. To investigate the physical full width at half maximum (FWHM) and penumbra width (80-20%), the single beam profiles and mean energy were obtained using BEAMnrc/BEAMDP code. Based on the results that there was an increase of the full width at half maximum (FWHM) with the size of the collimator on the beam profile and the mean energy. The FHWM of the energy fluence profile for the collimator 4, 8, and 16 mm are 2.22, 3.90, and 7.80 mm, respectively. The FHWM of the fluence profile for the collimator sizes 4, 8, and 16 mm are 1.87, 3.87, and 7.80 mm. The FHWM of mean energy as a function of position for the collimator sizes 4, 8, and 16 mm are 2.56, 4.20, and 7.94 mm, respectively. The maximum mean energy for the collimators decreases as much as 0.0185 MeV.

Dose-volume Histograms Comparison From Two Different Beam Configurations In The Monte Carlo Simulation of 6MV Clinac2100 Using PRIMO Program

The Monte Carlo (MC) method is widely considered to provide the most accurate dose distribution in external beam radiotherapy. A self-contained, full MC linac simulator and dose calculator, PRIMO is a windows-based, freely-distributed software which allows the importation of external compliant phase-space files among its many capabilities. Using PRIMO, a Varian Clinac2100 is simulated at 6 MV nominal energy for 108 number of histories. The dose calculations are compared using two different initial electron beam configurations. The tuned-beam profile has initial electron beam energy of 6.26 MeV, 0.150 MeV full-width-half-maximum (FWHM), 0.150 cm focal spot FWHM and 2 degrees beam divergence. The default configuration given in PRIMO has 5.40 MeV as initial electron beam energy and zero values for all the other beam parameters. A brain computerized tomography (CT) volume is imported into PRIMO and using its contouring tools, these structures are delineated: a hypothetical grosstumor volume (GTV), the brain, the left and right eye lens. The CT volume is irradiated with two parallel-opposed fields of size 10 cm by 10 cm with the brain conformed using the 52-leaf MLC. The dose-volume histogram (DVH) of the GTV and the brain for the tuned-beam profile has a larger extent into the high dose region compared to the DVHs of the default beam configuration. The dose given to 95% of the volume (D95) and the percentage of the volume receiving a dose equal to 95% (V95) are also compared. The GTV for the default beam profile has D95 = 86.31% with V95 = 0.44%, while for the tuned-beam profile, D95 = 88.30% with V95 = 0.95%. For the brain, D95 = 2.55% with V95 = 1.13% for the default profile while the tuned beam profile gives D95 = 2.43% with V95 = 2.33%. For the left [right] lens, D95 = 0.15% [D95 = 0.11%] for the default profile and D95 = 0.16% [D95 = 0.081%] for the tuned-beam profile. Both beam profiles give V95 = 0% for the left- and right- eye lens.

Automatic phase space generation for Monte Carlo calculations of intensity modulated particle therapy

Monte Carlo (MC) is generally considered as the most accurate dose calculation tool for particle therapy. However, a proper description of the beam particle kinematics is a necessary input for a realistic simulation. Such a description can be stored in phase space (PS) files for different beam energies. A PS file contains kinetic information such as energies, positions and travelling directions for particles traversing a plane perpendicular to the beam direction. The accuracy of PS files plays a critical role in the performance of the MC method for dose calculations. A PS file can be generated with a set of parameters describing analytically the beam kinematics. However, determining such parameters can be tedious and time consuming. Thus, we have developed an algorithm to obtain those parameters automatically and efficiently. In this paper, we presented such an algorithm and compared dose calculations using PS automatically generated for the Shanghai Proton and Heavy Ion Center (SPHIC) with measurements. The gamma-index for comparing calculated depth dose distributions (DDD) with measurements are above 96.0% with criterion 0.6%/0.6 mm. For each single energy, the mean difference percentage between calculated lateral spot sizes at 5 different locations along beam direction and measurements are below 3.5%.

Эффективный способ моделирования системы Leksell Gamma Knife Perfexion методом поворота в файле фазового пространства

Purpose: To develop an effective method of Monte Carlo simulation of the GammaKnife Perfexion system by rotating particles in the phase space file (PSF). This method does not require simulating of all 192 sources that are distributed in the conical form of the Perfexion collimator. The simulation was performed only for 5 out of 192 sources for each collimator size.&#x0D; Material and methods: Monte Carlo simulation of dose distribution for previous models of GammaKnife system requires phase space file for only one source, since this phase space is identical for all the 201 sources. The Perfexion model is more complex due to the non-coaxial positions of the sources and the complexity of the collimator system itself.&#x0D; In this work, we present an effective method to simulate the Perfexion model using a phase space file. Penelope Monte Carlo code was used to perform this simulation. In this method, the PSF was obtained for one source in each ring, resulting in five files for each collimator size. PSF for other sources were created by azimuthal redistribution of particles, in the obtained PSF, by rotation around the Z-axis. The phase space files of the same ring were then stored together in a single file.&#x0D; Results: The paper presented MC simulation using the azimuthal redistribution of particles in the phase space file by rotation around the Z-axis. The simulation has been validated comparing the dose profiles and output factors with the data of the algorithm TMR10 planning system Leksell Gamma Plan (LGP) in a homogeneous environment. The acceptance criterion between TMR10 and Monte Carlo calculations for the profiles was based on the gamma index (GI). Index values more than one were not detected in all cases, which indicates a good agreement of results. The differences between the output factors obtained in this work and the TMR10 data for collimators 8 mm and 4 mm are 0.74 and 0.73 %, respectively.&#x0D; Conclusion: In this work successfully implemented an effective method of simulating the Leksell Gamma knife Perfexion system. The presented method does not require modeling for all 192 sources distributed in the conical form of the Perfexion collimator. The simulation was performed for only five sources for each collimator and their files PSF were obtained. These files were used to create the PSF files for other sources by azimuthal redistribution of particles, in these files, by rotation around the Z-axis providing correct calculations of dose distributions in a homogeneous medium for 16, 8 and 4 mm collimators.

Monte Carlo studies on dose conversion factors from graphite to water for high energy X-ray beams

An evaluation of the dose conversion factor from graphite to water is essential to measure the absorbed dose to water using a graphite calorimeter. However, it is difficult to obtain the conversion factors from measurements. In the present study, Monte Carlo (MC) studies were performed to evaluate the absorbed dose conversion factor from graphite to water. Based on the Electron Gamma Shower MC codes from National Research Council Canada, the modeling of the medical linear accelerator gantry at the Korea Research Institute of Standards and Science (KRISS) was optimized to obtain phase space files for evaluation of the absorbed dose conversion factor from graphite to water. The evaluated dose conversion factors for the KRISS graphite calorimeter were in the range of 1.040–1.085 for beams of 6–25 MV. A linear correlation between the dose conversion factors and tissue phantom ratio in water at depths of 20 cm and 10 cm (TPR20,10) of each X-ray beam was observed. To understand sensitivity of the obtained dose conversion factors, case studies were performed under different conditions in the simulation, such as different electron stopping power, assumed graphite bulk density, and different MC code. The uncertainties of the dose conversion factors will be finalized considering the results of sensitivity studies in the present study and additional studies in the future.

Dosimetric evaluation of the Leksell GammaPlan™ Convolution dose calculation algorithm

The dosimetric accuracy of the Leksell GammaPlan Convolution calculation algorithm was evaluated through comparison with corresponding Monte Carlo (MC) dosimetric results. MC simulations were based on generated sector phase space files for the 4 mm, 8 mm and 16 mm collimator sizes, using a previous comprehensive Gamma Knife Perfexion™ source model and validated using film dosimetry. Test cases were designed for the evaluation of the Convolution algorithm involving irradiation of homogeneous and inhomogeneous phantom geometries mimicking clinical cases, with radiation fields created using one sector (single sector), all sectors with the same (single shot) or different (composite shot) collimator sizes. Dose calculations using the Convolution algorithm were found to be in excellent agreement (gamma pass rate greater than 98%, applying 1%/1 mm local dose difference and distance agreement criteria), with corresponding MC calculations, indicating the accuracy of the Convolution algorithm in homogeneous and heterogeneous model geometries. While of minor clinical importance, large deviations were observed for the voxels laying inside air media. The calculated beam on times using the Convolution algorithm were found to increase (up to 7%) relative to the TMR 10 algorithm currently used in clinical practice, especially in a test case mimicking a brain metastasis close to the skull, in excellent agreement with corresponding MC calculations.