Validation of PRIMO Monte Carlo Model of Clinac®iX 6MV Photon Beam

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.

Dosimetric accuracy of Acuros® XB and AAA algorithms for stereotactic body radiotherapy (SBRT) lung treatments: evaluation with PRIMO Monte Carlo code

Purpose:

The study aimed to compare the dosimetric performance of Acuros® XB (AXB) and anisotropic analytical algorithm (AAA) for lung SBRT plans using Monte Carlo (MC) simulations.

Methods:

We compared the dose calculation algorithms AAA and either of the dose reporting modes of AXB (dose to medium (AXB-Dm) or dose to water (AXB-Dw)) algorithms implemented in Eclipse® (Varian Medical Systems, Palo Alto, CA) Treatment planning system (TPS) with MC. PRIMO code was used for the MC simulations. The TPS-calculated dose profiles obtained with a multi-slab heterogeneity phantom were compared to MC. A lung phantom with a tumour was used to validate TPS algorithms using different beam delivery techniques. 2D gamma values obtained from Gafchromic film measurements in the tumour isocentre plane were compared with TPS algorithms and MC. Ten VMAT SBRT plans generated in TPS with each algorithm were recalculated with a PRIMO MC system for identical beam parameters for the clinical plan validation. A dose–volume histogram (DVH) based plan comparison and a 3D global gamma analysis were performed.

Results:

AXB demonstrated better agreement with MC and film measurements in the lung phantom validation, with good agreement in PDD, profiles and gamma analysis. AAA showed an overestimated PDD, a significant difference in dose profiles and a lower gamma pass rate near the field borders. With AAA, there was a dose overestimation at the periphery of the tumour. For clinical plan validation, AXB demonstrated higher agreement with MC than AAA.

Conclusions:

AXB provided better agreement with MC than AAA in the phantom and clinical plan evaluations.

Hybrid Cerenkov-scintillation detector validation using Monte Carlo simulations

Objective.This study aimed at investigating through Monte Carlo simulations the limitations of a novel hybrid Cerenkov-scintillation detector and the associated method for irradiation angle measurements.Approach.Using Monte Carlo simulations, previous experimental irradiations of the hybrid detector with a linear accelerator were replicated to evaluate its general performances and limitations. Cerenkov angular calibration curves and irradiation angle measurements were then compared. Furthermore, the impact of the Cerenkov light energy dependency on the detector accuracy was investigated using the energy spectra of electrons travelling through the detector.Main results.Monte Carlo simulations were found to be in good agreement with experimental values. The irradiation angle absolute mean error was found to be less than what was obtained experimentally, with a maximum value of 1.12° for the 9 MeV beam. A 0.4% increase of the ratio of electrons having an energy below 1 MeV to the total electrons was found to impact the Cerenkov light intensity collected as a function of the incident angle. The effect of the Cerenkov intensity variation on the measured angle was determined to vary according to the slope of the angular calibration curve. While the contribution of scattered electrons with a lower energy affects the detector accuracy, the greatest discrepancies result from the limitations of the calculation method and the calibration curve itself.Significance.A precise knowledge of the limitations of the hybrid detector and the irradiation angle calculation method is crucial for a clinical implementation. Moreover, the simulations performed in this study also corroborate hypotheses made regarding the relations between multiple Cerenkov dependencies and observations from the experimental measurements.

Monte Carlo simulation of linac using PRIMO

Background

Monte Carlo simulation is considered as the most accurate method for dose calculation in radiotherapy. PRIMO is a Monte-Carlo program with a user-friendly graphical interface.

Material and method

A VitalBeam with 6MV and 6MV flattening filter free (FFF), equipped with the 120 Millennium multileaf collimator was simulated by PRIMO. We adjusted initial energy, energy full width at half maximum (FWHM), focal spot FWHM, and beam divergence to match the measurements. The water tank and ion-chamber were used in the measurement. Percentage depth dose (PDD) and off axis ratio (OAR) were evaluated with gamma passing rates (GPRs) implemented in PRIMO. PDDs were matched at different widths of standard square fields. OARs were matched at five depths. Transmission factor and dose leaf gap (DLG) were simulated. DLG was measured by electronic portal imaging device using a sweeping gap method.

Result

For the criterion of 2%/2 mm, 1%/2 mm and 1%/1 mm, the GPRs of 6MV PDD were 99.33–100%, 99–100%, and 99–100%, respectively; the GPRs of 6MV FFF PDD were 99.33–100%, 98.99–99.66%, and 97.64–98.99%, respectively; the GPRs of 6MV OAR were 96.4–100%, 90.99–100%, and 85.12–98.62%, respectively; the GPRs of 6MV FFF OAR were 95.15–100%, 89.32–100%, and 87.02–99.74%, respectively. The calculated DLG matched well with the measurement (6MV: 1.36 mm vs. 1.41 mm; 6MV FFF: 1.07 mm vs. 1.03 mm, simulation vs measurement). The transmission factors were similar (6MV: 1.25% vs. 1.32%; 6MV FFF: 0.8% vs. 1.12%, simulation vs measurement).

Conclusion

The calculated PDD, OAR, DLG and transmission factor were all in good agreement with measurements. PRIMO is an independent (with respect to analytical dose calculation algorithm) and accurate Monte Carlo tool.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02149-5.

Comparison of PRIMO Monte Carlo code and Eclipse treatment planning system in calculation of dosimetric parameters in brain cancer radiotherapy

Background

It is important to evaluate the dose calculated by treatment planning systems (TPSs) and dose distribution in tumor and organs at risk (OARs). The aim of this study is to compare dose calculated by the PRIMO Monte Carlo code and Eclipse TPS in radiotherapy of brain cancer patients.

Materials and methods

PRIMO simulation code was used to simulate a Varian Clinac 600C linac. The simulations were validated for the linac by comparison of the simulation and measured results. In the case of brain cancer patients, the dosimetric parameters obtained by the PRIMO code were compared with those calculated by Eclipse TPS. Gamma function analysis with 3%, 3 mm criteria was utilized to compare the dose distributions. The evaluations were based on the dosimetric parameters for the planning target volume (PTV) and OAR including D min, D mean, and D max, homogeneity index (HI), and conformity index (CI).

Results

The gamma function analysis showed a 98% agreement between the results obtained by the PRIMO code and measurement for the percent depth dose (PDD) and dose profiles. The corresponding value in comparing the dosimetric parameters from PRIMO code and Eclipse TPS for the brain patients was 94%, on average. The results of the PRIMO simulation were in good agreement with the measured data and Eclipse TPS calculations.

Conclusions

Based on the results of this study, the PRIMO code can be utilized to simulate a medical linac with good accuracy and to evaluate the accuracy of treatment plans for patients with brain cancer.

Determination and validation of the initial beam parameters of Elekta Agility collimator head by Monte Carlo simulations

The availability of geometrical, physical, and initial beam parameters for Monte Carlo (MC) simulations of the Elekta Agility collimator head has become very difficult due to the proprietary nature of this data. This study presents strategies to independently determine the geometrical and physical properties of the components and initial beam parameters of the Agility collimator head for full beam simulations and postulates a benchmarking process using the EGSnrc MC toolkit. Target material of W (90%) and Re (10%) of 0.09 cm thickness, flattening filter of 1.77 cm thick stainless steel placed on 0.5 cm Al disc, and primary and secondary collimators of Tungsten alloy have been found to best fit the Agility head. The initial beam energy of 6.0 MeV with a radial distribution given as full-width half maxima (FWHM) of 0.301 cm (crossline) × 0.201 cm (inline) for 6 MV beam with a mean angular spread of 1.34° has been found best fitting the model. Variations of 0.29% and 0.59% have been noted in the measured and calculated values of TPR_20,10 and D₁₀ respectively. More than 90% dose points for all simulations passed the 2D gamma criteria of 3% DD, 3 mm DTA. This MC model of the Agility head can be used for dose calculation and validation of advanced treatment techniques.

Determination of dose distributions by monte-carlo simulation of 6 MV photon beam of varian vitalbeam accelerator using geant4 multithreaded code

Background:

Accuracy of dose delivery in radiation therapy is a primary requirement for effective cancer treatment. In practice, dose delivery accuracy of ±5% is desired. To achieve this accuracy level, an accurate method for calculating the dose distributions in the tumor volume is required. Monte-Carlo method is one of the methods considered to be the most accurate for calculating dose distributions.

Materials and Methods:

G4 linac-MT code was used to simulate a 6 MV photon beam. The initial electron beam parameters were tuned to validate the beam modeling from depth doses and beam profile. The dose distributions measured in water phantom were compared to the calculated dose distributions based on gamma index criterion.

Results:

The beam tuning showed the initial electron energy, sigma and full width at half maximum of 6.2 MeV, 0.8 MeV, and 1.18 mm, respectively, best match the measured dose distributions. The gamma index tests showed the calculated depth doses and beam profile were generally comparable with measurements, passing the standard acceptance criterion of 2%/2 mm. The simulated photon beam was justified by the index of beam quality, which showed excellent agreement with measured doses with a discrepancy of 0.1%.

Conclusion:

The observed agreement confirm the accuracy of the simulated 6 MV photon beam. It can therefore be used as radiation source for calculating dose distributions and further investigations aimed at improving dose delivery and planning in cancer patients.

Validation of Monte carlo Geant4 multithreading code for a 6 MV photon beam of varian linac on the grid computing

Aim

To evaluate the computation time efficiency of the multithreaded code (G4Linac-MT) in the dosimetry application, using the high performance of the HPC-Marwan grid to determine with high accuracy the initial parameters of the 6 MV photon beam of Varian CLINAC 2100C.

Background

The difficulty of Monte Carlo methods is the long computation time, this is one of the disadvantages of the Monte Carlo methods.

Materials and methods

Calculations are performed by the multithreaded code G4Linac-MT and Geant4.10.04.p02 using the HPC-Marwan computing grid to evaluate the computing speed for each code. The multithreaded version is tested in several CPUs to evaluate the computing speed according to the number of CPUs used. The results were compared to the measurements using different types of comparisons, TPR_20.10, penumbra, mean dose error and gamma index.

Results

The results obtained for this work indicate a much higher computing time saving for the G4Linac-MT version compared to the Geant4.10.04 version, the computing time decreases with the number of CPUs used, can reach about 12 times if 64CPUs are used. After optimization of the initial electron beam parameters, the results of the dose simulations obtained for this work are in very good agreement with the experimental measurements with a mean dose error of up to 0.41% on the PDDs and 1.79% on the lateral dose.

Conclusions

The gain in computation time leads us to perform Monte Carlo simulations with a large number of events which gives a high accuracy of the dosimetry results obtained in this work.