A clinical study of lung cancer dose calculation accuracy with Monte Carlo simulation

Correction method for ionization chamber dosimetry in flattening filter free radiotherapy based on Monte Carlo simulation

BACKGROUND

The clinical use of flattening filter free (FFF) radiotherapy has significantly increased in recent years due to its effective enhancement of dose rates and reduction of scatter dose. A proposal has been made to adjust the incident electron angle of the accelerator to expand the application of FFF beams in areas such as large planning target volumes (PTVs). However, the inherent softening characteristics and non-uniformity of lateral dose distribution in FFF beams inevitably lead to increased dosimetry errors, especially for ionization chambers widely used in clinical practice, which may result in serious accidents during FFF radiotherapy.

PURPOSE

This study constructs a comprehensive Monte Carlo model that encompasses not only conventional FFF beams but also incorporates FFF beams with varying incident electron angles, to investigate dosimetry errors and correction methods in FFF radiotherapy.

METHODS

We have innovatively introduced a FFF output correction factor (k Q F F F , Q W F F ${k}_{{Q}_{FFF},{Q}_{WFF}}$ ) to address dosimetry errors in various ionization chambers under different incident electron angle conditions in FFF beams. The primary variations ink Q F F F , Q W F F ${k}_{{Q}_{FFF},{Q}_{WFF}}$ were analytically determined to result from changes ins w , a i r ${s}_{w,air}$ and the perturbation correction terms of the ionization chamber.

RESULTS

Ionization chambers with smaller sensitive volumes typically exhibit reduced dosimetry errors. Our findings indicate that for ionization chambers with sensitive volumes ranging from 0.016  to 0.125 cm3, the dosimetry error under various FFF beam conditions consistently remains below 1.15%. This study provides crucial guidance for selecting appropriate ionization chambers in FFF radiotherapy.

CONCLUSION

A correlation was established between the absorbed dose to water in beams with a flattening filter (WFF) and those without (FFF), defined by the FFF output factor (O F Q F F F , Q W F F $O{F}_{{Q}_{FFF},{Q}_{WFF}}$ ). Using the proposed Monte Carlo model, theO F Q F F F , Q W F F $O{F}_{{Q}_{FFF},{Q}_{WFF}}$ can be derived and applied to theoretically calculate the absorbed dose to water in FFF beams at varying incident electron angles, with a relative standard uncertainty of 0.2. This study provides a valuable reference for clinical dose measurements and crucial support for establishing dose calibration standards in FFF radiotherapy.

Dosimetric Impact of Prescription Point Placement in Heterogeneous Medium for Conformal Radiotherapy Dose Calculation with Various Algorithms

Objective

The aim of the study is to compare the accuracy of dose calculation for different dose calculation algorithms with different prescription points (air, tissue, air-tissue interface in carcinoma lung patients and bone, tissue, and bone-tissue interface in carcinoma buccal Mucosa tumors).

Materials and Methods

Forty-one patients with carcinoma lung and buccal mucosa were retrospectively selected for this study. A three-dimensional conformal radiotherapy reference plan was created using the prescription point in the tissue with Monte Carlo (MC) algorithms for both the groups of patients. The reference plan was modified by changing the prescription point and algorithms in the tissue, air, air-tissue interface for lung patients and tissue, bone, and bone-tissue interface for buccal mucosa patients. The dose received by the target volume and other organs at risk (OAR) structures was compared. To find out the statistical difference between different prescription points and algorithms, the statistical tests were performed with repeated measures ANOVA.

Results

The target volume receiving 95% dose coverage in lung patients decreased to -3.08%, -5.75%, and -1.87% in the dose prescription point at the air-tissue interface with the dose calculation algorithms like MC, collapsed cone (CC), and pencil beam (PB), respectively, compared to that of the MC tissue. Spinal cord dose was increased in the CC and PB algorithms in all prescription points in patients with lung and buccal mucosa. OAR dose calculated by PB in all prescription points showed a significant deviation compared to MC tissue prescription point.

Conclusion

This study will help demonstrate the accuracy of dose calculation for the different dose prescription points with the different treatment algorithms in radiotherapy treatment planning.

Application of a linear interpolation algorithm in radiation therapy dosimetry for 3D dose point acquisition

Air-vented ion chambers are generally used in radiation therapy dosimetry to determine the absorbed radiation dose with superior precision. However, in ion chamber detector arrays, the number of array elements and their spacing do not provide sufficient spatial sampling, which can be overcome by interpolating measured data. Herein, we investigated the potential principle of the linear interpolation algorithm in volumetric dose reconstruction based on computed tomography images in the volumetric modulated arc therapy (VMAT) technique and evaluated how the ion chamber spacing and anatomical mass density affect the accuracy of interpolating new data points. Plane measurement doses on 83 VMAT treatment plans at different anatomical sites were acquired using Octavius 729, Octavius1500, and MatriXX ion chamber detector arrays, followed by the linear interpolation to reconstruct volumetric doses. Dosimetric differences in planning target volumes (PTVs) and organs at risk (OARs) between treatment planning system and reconstruction were evaluated by dose volume histogram metrics. The average percentage dose deviations in the mean dose (D_mean) of PTVs reconstructed by 729 and 1500 arrays ranged from 4.7 to 7.3% and from 1.5 to 2.3%, while the maximum dose (D_max) counterparts ranged from 2.3 to 5.5% and from 1.6 to 7.6%, respectively. The average percentage dose/volume deviations of mixed PTVs and OARs in the abdomen/gastric and pelvic sites were 7.6%, 3.5%, and 7.2%, while mediastinum and lung plans showed slightly larger values of 8.7%, 5.1%, and 8.9% for 729, 1500, and MatriXX detector arrays, respectively. Our findings indicated that the smaller the spacing between neighbouring detectors and the more ion chambers present, the smaller the error in interpolating new data points. Anatomical regions with small local mass density inhomogeneity were associated with superior dose reconstruction. Given a large mass density difference in the various human anatomical structures and the characteristics of the linear interpolation algorithm, we suggest that an alternative data interpolation method should be used in radiotherapy dosimetry.

A Monte Carlo study on dose distribution of an orthovoltage radiation therapy system

It is important to plan radiotherapy treatment and establish optimal dose distribution to reduce the chances of side effects and injury. Because there are no commercially available tools for calculating dose distribution in orthovoltage radiotherapy in companion animals, we developed an algorithm to accomplish this and verified its characteristics using tumor disease cases. First, we used the Monte Carlo method to develop an algorithm to calculate the dose distribution of orthovoltage radiotherapy (280 kVp; MBR-320, Hitachi Medical Corporation, Tokyo, Japan) using BEAMnrc at our clinic. Using development of Monte Carlo method, dose distribution for tumor and normal organs were evaluated in brain tumors, squamous cell carcinomas of the head, and feline nasal lymphomas. In all cases of brain tumors, the mean dose delivered to the GTV ranged from 36.2 to 76.1% of the prescribed dose due to the decrease through the skull. In the nasal lymphoma in cats, the eyes with covered a 2 mm-thick lead plate, the respective average dose to the eyes was 71.8% and 89.9% less than that to the uncovered eyes. The findings may be useful for informed decision making in orthovoltage radiotherapy with more effective and targeted irradiation and data collection allowing detailed informed consent.

Monte Carlo evaluation of target dose coverage in lung stereotactic body radiation therapy with flattening filter-free beams

Aim

Previous studies showed that replacing conventional flattened beams (FF) with flattening filter-free (FFF) beams improves the therapeutic ratio in lung stereotactic body radiation therapy (SBRT), but these findings could have been impacted by dose calculation uncertainties caused by the heterogeneity of the thoracic anatomy and by respiratory motion, which were particularly high for target coverage. In this study, we minimized such uncertainties by calculating doses using high-spatial-resolution Monte Carlo and four-dimensional computed tomography (4DCT) images. We aimed to evaluate more reliably the benefits of using FFF beams for lung SBRT.

Materials and methods

For a cohort of 15 patients with early stage lung cancer that we investigated in a previous treatment planning study, we recalculated dose distributions with Monte Carlo using 4DCT images. This included fifteen FF and fifteen FFF treatment plans.

Results

Compared to Monte Carlo, the treatment planning system (TPS) over-predicted doses in low-dose regions of the planning target volume. For most patients, replacing FF beams with FFF beams improved target coverage, tumor control, and uncomplicated tumor control probabilities.

Conclusions

Monte Carlo tends to reveal deficiencies in target coverage compared to coverage predicted by the TPS. Our data support previously reported benefits of using FFF beams for lung SBRT.

Experimental validation of a linac head Geant4 model under a grid computing environment

Background and purpose: This work aims to present the strategy to simulate a clinical linear accelerator based on the geometry provided by the manufacturer and summarize the corresponding experimental validation. Simulations were performed with the Geant4 Monte Carlo code under a grid computing environment. The objective of this contribution is reproducing therapeutic dose distributions in a water phantom with an accuracy less than 2%. Materials and methods: A Geant4 Monte Carlo model of an Elekta Synergy linear accelerator has been established, the simulations were launched in a large grid computing platform. Dose distributions are calculated for a 6 MV photon beam with treatment fields ranging from 5 × 5 cm2 to 20 × 20 cm 2 at a source - surface distance of 100 cm. Results: A high degree of agreement is achieved between the simulation results and the measured data, with dose differences of about 1.03% and 1.96% for the percentage depth dose curves and lateral dose profiles, respectively. This agreement is evaluated by the gamma index comparisons. Over 98% of the points for all simulations meet the restrictive acceptability criteria of 2%/2 mm. Conclusion: We have demonstrated the possibility to establish an accurate linac head Monte Carlo model for dose distribution simulations and quality assurance. Percentage depth dose curves and beam quality indices are in perfect agreement with the measured data with an accuracy of better than 2%.

Positron emission tomography guided dose painting by numbers of lung cancer: Alanine dosimetry in an anthropomorphic phantom

•

DPBN was delivered to a phantom based on the anatomy of a lung cancer patient examined by FDG PET/CT prior to radiotherapy.

•

Alanine dosimetry showed that DPBN can be delivered with high accuracy to the tumour in the anthropomorphic phantom.

•

For regions outside the tumour, high correspondence between planned and delivered doses were also found.

•

Positioning errors can lead to large deviations and potentially sub-optimal tumor doses.

Background and purpose

Dose painting by numbers (DPBN) require a high degree of dose modulation to fulfill the image-based voxel wise dose prescription. The aim of this study was to assess the dosimetric accuracy of ¹⁸F-fluoro-2-deoxy-glucose positron emission tomography(¹⁸F-FDG-PET)-based DPBN in an anthropomorphic lung phantom using alanine dosimetry.

Materials and methods

A linear dose prescription based on ¹⁸F-FDG-PET image intensities within the gross tumor volume (GTV) of a lung cancer patient was employed. One DPBN scheme with low dose modulation (Scheme A; minimum/maximum fraction dose to the GTV 2.92/4.26 Gy) and one with a high modulation (Scheme B; 2.81/4.52 Gy) were generated. The plans were transferred to a computed tomograpy (CT) scan of a thorax phantom based on CT images of the patient. Using volumetric modulated arc therapy (VMAT), DPBN was delivered to the phantom with embedded alanine dosimeters. A plan was also delivered to an intentionally misaligned phantom. Absorbed doses at various points in the phantom were measured by alanine dosimetry.

Results

A pointwise comparison between GTV doses from prescription, treatment plan calculation and VMAT delivery showed high correspondence, with a mean and maximum dose difference of <0.1 Gy and 0.3 Gy, respectively. No difference was found in dosimetric accuracy between scheme A and B. The misalignment caused deviations up to 1 Gy between prescription and delivery.

Conclusion

DPBN can be delivered with high accuracy, showing that the treatment may be applied correctly from a dosimetric perspective. Still, misalignment may cause considerable dosimetric erros, indicating the need for patient immobilization and monitoring.

Clinical verification of treatment planning dose calculation in lung SBRT with GATE Monte Carlo simulation code

PURPOSE

This study aims to use GATE/Geant4 simulation code to evaluate the performance of dose calculations with Anisotropic Analytical Algorithm (AAA) in the context of lung SBRT for complex treatments considering images of patients.

METHODS

Four cases of non-small cell lung cancer treated with SBRT were selected for this study. Irradiation plans were created with AAA and recalculated end to end using Monte Carlo (MC) method maintaining field configurations identical to the original plans. Each treatment plan was evaluated in terms of PTV and organs at risk (OARs) using dose-volume histograms (DVH). Dosimetric parameters obtained from DVHs were used to compare AAA and MC.

RESULTS

The comparison between the AAA and MC DVH using gamma analysis with the passing criteria of 3%/3% showed an average passing rate of more than 90% for the PTV structure and 97% for the OARs. Tightening the criteria to 2%/2% showed a reduction in the average passing rate of the PTV to 86%. The agreement between the AAA and MC dose calculations for PTV dosimetric parameters (V₁₀₀; V₉₀; Homogeneity index; maximum, minimum and mean dose; CI_Paddick and D_2cm) was within 18.4%. For OARs, the biggest differences were observed in the spinal cord and the great vessels.

CONCLUSIONS

In general, we did not find significant differences between AAA and MC. The results indicate that AAA could be used in complex SBRT cases that involve a larger number of small treatment fields in the presence of tissue heterogeneities.

The status of medical physics in radiotherapy in China

PURPOSE

To present an overview of the status of medical physics in radiotherapy in China, including facilities and devices, occupation, education, research, etc. MATERIALS AND METHODS: The information about medical physics in clinics was obtained from the 9-th nationwide survey conducted by the China Society for Radiation Oncology in 2019. The data of medical physics in education and research was collected from the publications of the official and professional organizations.

RESULTS

By 2019, there were 1463 hospitals or institutes registered to practice radiotherapy and the number of accelerators per million population was 1.5. There were 4172 medical physicists working in clinics of radiation oncology. The ratio between the numbers of radiation oncologists and medical physicists is 3.51. Approximately, 95% of medical physicists have an undergraduate or graduate degrees in nuclear physics and biomedical engineering. 86% of medical physicists have certificates issued by the Chinese Society of Medical Physics. There has been a fast growth of publications by authors from mainland of China in the top international medical physics and radiotherapy journals since 2018.

CONCLUSIONS

Demand for medical physicists in radiotherapy increased quickly in the past decade. The distribution of radiotherapy facilities in China became more balanced. High quality continuing education and training programs for medical physicists are deficient in most areas. The role of medical physicists in the clinic has not been clearly defined and their contributions have not been fully recognized by the community.

Monte Carlo as quality control tool of stereotactic body radiation therapy treatment plans

PURPOSE/OBJECTIVE

The objective of this study was to verify the accuracy of treatment plans of stereotactic body radiation therapy (SBRT) and to verify the feasibility of the use of Monte Carlo (MC) as quality control (QC) on a daily basis.

MATERIAL/METHODS

Using EGSnrc, a MC model of Agility™ linear accelerator was created. Various measurements (Percentage depth dose (PDD), Profiles and Output factors) were done for different fields sizes from 1x1 up to 40x40 (cm²). An iterative model optimization was performed to achieve adequate parameters of MC simulation. 40 SBRT patient's dosimetry plans were calculated by Monaco™ 3.1.1. CT images, RT-STRUCT and RT-PLAN files from Monaco™ being used as input for Moderato MC code. Finally, dose volume histogram (DVH) and paired t-tests for each contour were used for dosimetry comparison of the Monaco™ and MC.

RESULTS

Validation of MC model was successful, as <2% difference comparing to measurements for all field's sizes. The main energy of electron source incident on the target was 5.8 MeV, and the full width at half maximum (FWHM) of Gaussian electron source were 0.09 and 0.2 (cm) in X and Y directions, respectively. For 40 treatment plan comparisons, the minimum absolute difference of mean dose of planning treatment planning (PTV) was 0.1% while the maximum was 6.3%. The minimum absolute difference of Max dose of PTV was 0.2% while the maximum was 8.1%.

CONCLUSION

SBRT treatment plans of Monaco agreed with MC results. It possible to use MC for treatment plans verifications as independent QC tool.

Development of a Tongue Immobilization Device Using a 3D Printer for Intensity Modulated Radiation Therapy of Nasopharyngeal Cancer Patients

Purpose

This study aimed to reduce radiation doses to the tongue, a patient-specific semi-customized tongue immobilization device (SCTID) was developed using a 3D printer for helical tomotherapy (HT) of nasopharyngeal cancer (NPCa). Dosimetric characteristics and setup stability of the SCTID were compared with those of a standard mouthpiece (SMP).

Materials and Methods

For displacement and robust immobilization of the tongue, the SCTID consists of four parts: upper and lower tooth stoppers, tongue guider, tongue-tip position guide bar, and connectors. With the SCTID and SMP, two sets of planning computed tomography and HT plans were obtained for 10 NPCa patients. Dosimetric and geometric characteristics were compared. Position reproducibility of the tongue with SCTID was evaluated by comparing with planned dose and adaptive accumulated dose of the tongue and base of the tongue based on daily setup mega-voltage computed tomography.

Results

Using the SCTID, the tongue was effectively displaced from the planning target volume compared to the SMP. The median mucosa of the tongue (M-tongue) dose was significantly reduced (20.7 Gy vs. 27.8 Gy). The volumes of the M-tongue receiving a dose of 15 Gy, 30 Gy, and 45 Gy and the volumes of the mucosa of oral cavity and oropharynx (M-OC/OP) receiving a dose of 45 Gy and 60 Gy were significantly lower than using the SMP. No significant differences was observed between the planned dose and the accumulated adaptive dose in any dosimetric characteristics of the tongue and base of tongue.

Conclusion

SCTID can not only reduce the dose to the M-tongue and M-OC/OP dramatically, when compared to SMP, but also provide excellent reproducibility and easy visual verification.

Evaluation of the Lung Dose in Three-dimensional Conformal Radiation Therapy of Left-Sided Breast Cancer: A Phantom Study

Background:

Three-dimensional 3D-CRT: conformal radiation therapy is a selective modality in many radiotherapy centers for the treatment of breast cancer. One of the most common side effects of this method is radiation lung injury. Considering such an injury, lung dose deserves to be studied in depth.

Methods:

Computed tomography scan of a node-positive left-sided breast cancer woman was used for generating a thorax phantom. Ten thermoluminescent dosimeters (TLDs) were distributed evenly in the left lung of the phantom, and the phantom was scanned. The optimal plan, including supraclavicular and tangential fields, was created by the treatment planning system (TPS). The results of TLD dose measurements at the selected points in the phantom were compared to TPS dose calculations.

Results:

Lung doses calculated by TPS are significantly different from those measured by the TLDs (P = 0.007). The minimum and maximum differences were −0.91% and 4.46%, respectively. TLDs that were on the inner margin of the lung and breast tissue showed higher dose differences than the TLDs in the lung.

Conclusion:

The results of this study showed that TPS generally overestimated doses compared to TLD measurements due to incorrect beam modeling caused by contaminated electrons in the lung.

Pre-treatment 2D and 3D dosimetric verification of volumetric arc therapy. A correlation study between gamma index passing rate and clinical dose volume histogram

OBJECTIVES

To evaluate methods for the pre-treatment verification of volumetric modulated arc therapy (VMAT) based on the percentage gamma passing rate (%GP) and its correlation and sensitivity with percentage dosimetric errors (%DE).

METHODS

A total of 25 patients with prostate cancer and 15 with endometrial cancer were analysed. The %GP values of 2D and 3D verifications with different acceptance criteria (1%/1 mm, 2%/2 mm, and 3%/3 mm) were obtained using OmniPro and Compass. The %DE was calculated using a planned dose volume histogram (DVH) created in Monaco's treatment planning system (TPS), which relates radiation dose to tissue and the patient's predicted dose volume histogram in Compass. Statistical correlation between %GP and %DE was verified using Pearson's correlation coefficient. Sensitivity was calculated based on the receiver operating characteristics (ROC) curve. Plans were calculated using Collapsed Cone Convolution and the Monte Carlo algorithm.

RESULTS

The t-test results of the planned and estimated DVH showed that the mean values were comparable (P > 0.05). For the 3%/3 mm criterion, the average %GP was acceptable for the prostate and endometrial cancer groups, with average rates of 99.68 ± 0.49% and 99.03 ± 0.59% for 2D and 99.86 ± 0.39% and 99.53 ± 0.44% for 3D, respectively. The number of correlations was poor for all analysed data. The mean Pearson's R-values for prostate and endometrial cancer were < 0.45 and < 0.43, respectively. The area under the ROC curve for the prostate and endometrial cancer groups, was lower than 0.667.

CONCLUSIONS

Analysis of the %GP versus %DE values revealed only weak correlations between 2D and 3D verifications. DVH results obtained using the Compass system will be helpful in confirming that the analysed plans respect dosimetric constraints.

Development of a semi-customized tongue displacement device using a 3D printer for head and neck IMRT

Purpose

To reduce radiation doses to the tongue, a patient-specific semi-customized tongue displacement device (SCTDD) was developed using a 3D printer (3DP) for head and neck (H&N) radiation therapy (RT). Dosimetric characteristics of the SCTDD were compared with those of a standard mouthpiece (SMP).

Materials and methods

The SCTDD consists of three parts: a mouthpiece, connector with an immobilization mask, and tongue displacer, which can displace the tongue to the contralateral side of the planning target volume. Semi-customization was enabled by changing the thickness and length of the SCTDD. The instrument was printed using a 3DP with a biocompatible material. With the SCTDD and SMP, two sets of planning computed tomography (CT) and tomotherapy plans were obtained for seven H&N cancer patients. Dosimetric and geometric characteristics were compared.

Results

Using the SCTDD, the tongue was effectively displaced from the planning target volume without significant tongue volume change compared to the SMP. The median tongue dose was significantly reduced (29.6 Gy vs. 34.3 Gy). The volumes of the tongue receiving a dose of 15 Gy, 30 Gy, 35 Gy, 45 Gy, and 60 Gy were significantly lower than using the SMP.

Conclusion

The SCTDD significantly decreased the radiation dose to the tongue compared to the SMP, which may potentially reduce RT-related tongue toxicity.

Dosimetric impact of statistical uncertainty on Monte Carlo dose calculation algorithm in volumetric modulated arc therapy using Monaco TPS for three different clinical cases

Aim

To study the dosimetric impact of statistical uncertainty (SU) per plan on Monte Carlo (MC) calculation in Monaco™ treatment planning system (TPS) during volumetric modulated arc therapy (VMAT) for three different clinical cases.

Background

During MC calculation SU is an important factor to decide dose calculation accuracy and calculation time. It is necessary to evaluate optimal acceptance of SU for quality plan with reduced calculation time.

Materials and methods

Three different clinical cases as the lung, larynx, and prostate treated using VMAT technique were chosen. Plans were generated with Monaco™ V5.11 TPS with 2% statistical uncertainty. By keeping all other parameters constant, plans were recalculated by varying SU, 0.5%, 1%, 2%, 3%, 4%, and 5%. For plan evaluation, conformity index (CI), homogeneity index (HI), dose coverage to PTV, organ at risk (OAR) dose, normal tissue receiving dose ≥5 Gy and ≥10 Gy, integral dose (NTID), calculation time, gamma pass rate, calculation reproducibility and energy dependency were analyzed.

Results

CI and HI improve as SU increases from 0.5% to 5%. No significant dose difference was observed in dose coverage to PTV, OAR doses, normal tissue receiving dose ≥5 Gy and ≥10 Gy and NTID. Increase of SU showed decrease in calculation time, gamma pass rate and increase in PTV max dose. No dose difference was seen in calculation reproducibility and dependent on energy.

Conclusion

For VMAT plans, SU can be accepted from 1% to 3% per plan with reduced calculation time without compromising plan quality and deliverability by accepting variations in point dose within the target.

The dosimetric effect of electron density overrides in 3DCRT Lung SBRT for a range of lung tumor dimensions

The combined effects of lung tumor motion and limitations of treatment planning system dose calculations in lung regions increases uncertainty in dose delivered to the tumor and surrounding normal tissues in lung stereotactic body radiotherapy (SBRT). This study investigated the effect on plan quality and accuracy when overriding treatment volume electron density values. The QUASAR phantom with modified cork cylindrical inserts, each containing a simulated spherical tumor of 15‐mm, 22‐mm, or 30‐mm diameter, was used to simulate lung tumor motion. Using Monaco 5.1 treatment planning software, two standard plans (50% central phase (50%) and average intensity projection (AIP)) were compared to eight electron density overridden plans that focused on different target volumes (internal target volume (ITV), planning target volume (PTV), and a hybrid plan (HPTV)). The target volumes were set to a variety of electron densities between lung and water equivalence. Minimal differences were seen in the 30‐mm tumor in terms of target coverage, plan conformity, and improved dosimetric accuracy. For the smaller tumors, a PTV override showed improved target coverage as well as better plan conformity compared to the baseline plans. The ITV plans showed the highest gamma pass rate agreement between treatment planning system (TPS) and measured dose (P < 0.040). However, the low electron density PTV and HPTV plans also showed improved gamma pass rates (P < 0.035, P < 0.011). Low‐density PTV overrides improved the plan quality and accuracy for tumor diameters less than 22 mm only. Although an ITV override generated the most significant increase in accuracy, the low‐density PTV plans had the additional benefit of plan quality improvement. Although this study and others agreed that density overrides improve the treatment of SBRT, the optimal density override and the conditions under which it should be applied were found to be department specific, due to variations in commissioning and calculation methods.

Dosimetric comparison of pencil beam and Monte Carlo algorithms in conformal lung radiotherapy

PURPOSE

In this study, lung radiotherapy target volumes as well as critical organs such as the lungs, spinal cord, esophagus, and heart doses calculated using pencil beam (PB) and Monte Carlo (MC) algorithm-based treatment planning systems (TPSs) were compared. The main aim was the evaluation of calculated dose differences between the PB and MC algorithms in a highly heterogeneous medium.

METHODS

A total of 6 MV photon energy conformal treatment plans were created for a RANDO lung phantom using one PB algorithm-based Precise Plan Release 2.16 TPS and one MC algorithm-based Monaco TPS. Thermoluminescence dosimeters (TLDs) were placed into appropriate slices within the RANDO phantom and then irradiated with an Elekta-Synergy® Linear Accelerator for dose verification. Doses were calculated for the V5, V10, V20, and mean lung doses (MLDs) in bilateral lungs and D50, D98, D2, and mean doses in the target volume (planning target volume, PTV).

RESULTS

The minimum, maximum, and mean doses of the target volumes and critical organs in two treatment plans were compared using dose volume histograms (DVHs). The mean dose difference between the PB and MC algorithms for the PTV was 0.3%, whereas the differences in V5, V10, V20, and MLD were 12.5%, 15.8%, 14.4%, and 9.1%, respectively. The differences in PTV coverage between the two algorithms were 0.9%, 2.7% and 0.7% for D50, D98 and D2, respectively.

CONCLUSIONS

A comparison of the dose data acquired in this study reveals that the MC algorithm calculations are closer to the 60 Gy prescribed dose for PTV, while the difference between the PB and MC algorithms was found to be non-significant. Because of the major difference arising from the dose calculation techniques by TPS that was observed in the MLD with significant medium heterogeneity, we recommend the use of the MC algorithm in such heterogeneous sites.

CPU time optimization and precise adjustment of the Geant4 physics parameters for a VARIAN 2100 C/D gamma radiotherapy linear accelerator simulation using GAMOS

We have verified the GAMOS/Geant4 simulation model of a 6 MV VARIAN Clinac 2100 C/D linear accelerator by the procedure of adjusting the initial beam parameters to fit the percentage depth dose and cross-profile dose experimental data at different depths in a water phantom. Thanks to the use of a wide range of field sizes, from 2  ×  2 cm² to 40  ×  40 cm², a small phantom voxel size and high statistics, fine precision in the determination of the beam parameters has been achieved. This precision has allowed us to make a thorough study of the different physics models and parameters that Geant4 offers. The three Geant4 electromagnetic physics sets of models, i.e. Standard, Livermore and Penelope, have been compared to the experiment, testing the four different models of angular bremsstrahlung distributions as well as the three available multiple-scattering models, and optimizing the most relevant Geant4 electromagnetic physics parameters. Before the fitting, a comprehensive CPU time optimization has been done, using several of the Geant4 efficiency improvement techniques plus a few more developed in GAMOS.

Normal lung sparing Tomotherapy technique in stage III lung cancer

PURPOSE

Radiation pneumonitis (RP) has been a challenging obstacle in treating stage III lung cancer patients. Beam angle optimization (BAO) technique for Tomotherapy was developed to reduce the normal lung dose for stage III non-small cell lung cancer (NSCLC). Comparative analyses on plan quality by 3 different Intensity-modulated radiation therapy (IMRT) methods with BAO were done.

MATERIALS AND METHODS

Ten consecutive stage IIIB NSCLC patients receiving linac-based static IMRT (L-IMRT) with total 66 Gy in 33 fractions to the PTV were selected. Two additional Tomotherapy-based IMRT plans (helical beam (TH-IMRT) and static beam (TD-IMRT)) were generated on each patient. To reduce the normal lung dose, Beam angles were optimized by using complete and directional block functions in Tomotherapy based on knowledge based statistical analysis. Plan quality was compared with target coverage, normal organ sparing capability, and normal tissue complication probability (NTCP). Actual beam delivery times and risk of RP related with planning target volume (PTV) were also evaluated.

RESULTS

The best PTV coverage measured by conformity index and homogeneity index was achievable by TH-IMRT (0.82 and 1.06), followed by TD-IMRT (0.81 and 1.07) and L-IMRT (0.75 and 1.08). Mean lung dose was the lowest in TH-IMRT plan followed by TD-IMRT and L-IMRT, all of which were ≤20 Gy. TH-IMRT plan could significantly lower the lung volumes receiving low to medium dose levels: V5~30 when compared to L-IMRT plan; and V5~20 when compared to TD-IMRT plan, respectively. TD-IMRT plan was significantly better than L-IMRT with respects to V20 and V30 and there was no significant difference with respect to V40 among three plans. The NTCP of the lung was the lowest in TH-IMRT plan, followed by TD-IMRT and L-IMRT (6.42% vs. 6.53% vs. 8.11%). Beam delivery time was the shortest in TD-IMRT plan followed by L-IMRT. As PTV length increased, NTCP and Mean lung dose proportionally increased significantly in all three plans.

CONCLUSION

Advantageous profiles by TH-IMRT could be achieved by BAO by complete and directional block functions. Current observation could help radiation oncologists to make wise selection of IMRT method for stage IIIB NSCLC.

A clinical database to assess action levels and tolerances for the ongoing use of Mobius3D

In radiation therapy, calculation of dose within the patient contains inherent uncertainties, inaccuracies, limitations, and the potential for random error. Thus, point dose‐independent verification of such calculations is a well‐established process, with published data to support the setting of both action levels and tolerances. Mobius3D takes this process one step further with a full independent calculation of patient dose and comparisons of clinical parameters such as mean target dose and voxel‐by‐voxel gamma analysis. There is currently no published data to directly inform tolerance levels for such parameters, and therefore this work presents a database of 1000 Mobius3D results to fill this gap. The data are tested for normality using a normal probability plot and found to fit this distribution for three sub groups of data; Eclipse,iPlan and the treatment site Lung. The mean (μ) and standard deviation (σ) of these sub groups is used to set action levels and tolerances at μ ± 2σ and μ ± 3σ, respectively. A global (3%, 3 mm) gamma tolerance is set at 88.5%. The mean target dose tolerance for Eclipse data is the narrowest at ± 3%, whilst iPlan and Lung have a range of −5.0 to 2.2% and −1.8 to 5.0%, respectively. With these limits in place, future results failing the action level or tolerance will fall within the worst 5% and 1% of historical results and an informed decision can be made regarding remedial action prior to treatment.

Target dose conversion modeling from pencil beam (PB) to Monte Carlo (MC) for lung SBRT

BACKGROUND

A challenge preventing routine clinical implementation of Monte Carlo (MC)-based lung SBRT is the difficulty of reinterpreting historical outcome data calculated with inaccurate dose algorithms, because the target dose was found to decrease to varying degrees when recalculated with MC. The large variability was previously found to be affected by factors such as tumour size, location, and lung density, usually through sub-group comparisons. We hereby conducted a pilot study to systematically and quantitatively analyze these patient factors and explore accurate target dose conversion models, so that large-scale historical outcome data can be correlated with more accurate MC dose without recalculation.

METHODS

Twenty-one patients that underwent SBRT for early-stage lung cancer were replanned with 6MV 360° dynamic conformal arcs using pencil-beam (PB) and recalculated with MC. The percent D95 difference (PB-MC) was calculated for the PTV and GTV. Using single linear regression, this difference was correlated with the following quantitative patient indices: maximum tumour diameter (MaxD); PTV and GTV volumes; minimum distance from tumour to soft tissue (dmin); and mean density and standard deviation of the PTV, GTV, PTV margin, lung, and 2 mm, 15 mm, 50 mm shells outside the PTV. Multiple linear regression and artificial neural network (ANN) were employed to model multiple factors and improve dose conversion accuracy.

RESULTS

Single linear regression with PTV D95 deficiency identified the strongest correlation on mean-density (location) indices, weaker on lung density, and the weakest on size indices, with the following R(2) values in decreasing orders: shell2mm (0.71), PTV (0.68), PTV margin (0.65), shell15mm (0.62), shell50mm (0.49), lung (0.40), dmin (0.22), GTV (0.19), MaxD (0.17), PTV volume (0.15), and GTV volume (0.08). A multiple linear regression model yielded the significance factor of 3.0E-7 using two independent features: mean density of shell2mm (P = 1.6E-7) and PTV volume (P = 0.006). A 4-feature ANN model slightly improved the modeling accuracy.

CONCLUSION

Quantifiable density features were proposed, replacing simple central/peripheral location designation, which showed strong correlations with PB-to-MC target dose conversion magnitude, followed by lung density and target size. Density in the immediate outer and inner areas of the PTV showed the strongest correlations. A multiple linear regression model with one such feature and PTV volume established a high significance factor, improving dose conversion accuracy.

Dose differences in intensity-modulated radiotherapy plans calculated with pencil beam and Monte Carlo for lung SBRT

For patients with medically inoperable early‐stage non‐small cell lung cancer (NSCLC) treated with stereotactic body radiation therapy, early treatment plans were based on a simpler dose calculation algorithm, the pencil beam (PB) calculation. Because these patients had the longest treatment follow‐up, identifying dose differences between the PB calculated dose and Monte Carlo calculated dose is clinically important for understanding of treatment outcomes. Previous studies found significant dose differences between the PB dose calculation and more accurate dose calculation algorithms, such as convolution‐based or Monte Carlo (MC), mostly for three‐dimensional conformal radiotherapy (3D CRT) plans. The aim of this study is to investigate whether these observed dose differences also exist for intensity‐modulated radiotherapy (IMRT) plans for both centrally and peripherally located tumors. Seventy patients (35 central and 35 peripheral) were retrospectively selected for this study. The clinical IMRT plans that were initially calculated with the PB algorithm were recalculated with the MC algorithm. Among these paired plans, dosimetric parameters were compared for the targets and critical organs. When compared to MC calculation, PB calculation overestimated doses to the planning target volumes (PTVs) of central and peripheral tumors with different magnitudes. The doses to 95% of the central and peripheral PTVs were overestimated by 9.7%±5.6% and 12.0%±7.3%, respectively. This dose overestimation did not affect doses to the critical organs, such as the spinal cord and lung. In conclusion, for NSCLC treated with IMRT, dose differences between the PB and MC calculations were different from that of 3D CRT. No significant dose differences in critical organs were observed between the two calculations.

PACS number: 87.53.Ly