The use of non-standard CT conversion ramps for Monte Carlo verification of 6 MV prostate IMRT plans

Monte Carlo (MC) dose calculation algorithms have been widely used to verify the accuracy of intensity-modulated radiotherapy (IMRT) dose distributions computed by conventional algorithms due to the ability to precisely account for the effects of tissue inhomogeneities and multileaf collimator characteristics. Both algorithms present, however, a particular difference in terms of dose calculation and report. Whereas dose from conventional methods is traditionally computed and reported as the water-equivalent dose (D_w), MC dose algorithms calculate and report dose to medium (D_m). In order to compare consistently both methods, the conversion of MC D_m into D_w is therefore necessary.This study aims to assess the effect of applying the conversion of MC-based D_m distributions to D_w for prostate IMRT plans generated for 6 MV photon beams. MC phantoms were created from the patient CT images using three different ramps to convert CT numbers into material and mass density: a conventional four material ramp (CTCREATE) and two simplified CT conversion ramps: (1) air and water with variable densities and (2) air and water with unit density. MC simulations were performed using the BEAMnrc code for the treatment head simulation and the DOSXYZnrc code for the patient dose calculation. The conversion of D_m to D_w by scaling with the stopping power ratios of water to medium was also performed in a post-MC calculation process.The comparison of MC dose distributions calculated in conventional and simplified (water with variable densities) phantoms showed that the effect of material composition on dose-volume histograms (DVH) was less than 1% for soft tissue and about 2.5% near and inside bone structures. The effect of material density on DVH was less than 1% for all tissues through the comparison of MC distributions performed in the two simplified phantoms considering water. Additionally, MC dose distributions were compared with the predictions from an Eclipse treatment planning system (TPS), which employed a pencil beam convolution (PBC) algorithm with Modified Batho Power Law heterogeneity correction. Eclipse PBC and MC calculations (conventional and simplified phantoms) agreed well (<1%) for soft tissues. For femoral heads, differences up to 3% were observed between the DVH for Eclipse PBC and MC calculated in conventional phantoms. The use of the CT conversion ramp of water with variable densities for MC simulations showed no dose discrepancies (0.5%) with the PBC algorithm. Moreover, converting D_m to D_w using mass stopping power ratios resulted in a significant shift (up to 6%) in the DVH for the femoral heads compared to the Eclipse PBC one.Our results show that, for prostate IMRT plans delivered with 6 MV photon beams, no conversion of MC dose from medium to water using stopping power ratio is needed. In contrast, MC dose calculations using water with variable density may be a simple way to solve the problem found using the dose conversion method based on the stopping power ratio.     

Tissue composition and density impact on the clinical parameters for 125I prostate implants dosimetry

The MCNPX code was used to calculate the TG-43U1 recommended parameters in water and prostate tissue in order to quantify the dosimetric impact in 30 patients treated with ¹²⁵I prostate implants when replacing the TG-43U1 formalism parameters calculated in water by a prostate-like medium in the planning system (PS) and to evaluate the uncertainties associated with Monte Carlo (MC) calculations. The prostate density was obtained from the CT of 100 patients with prostate cancer. The deviations between our results for water and the TG-43U1 consensus dataset values were −2.6% for prostate V₁₀₀, −13.0% for V₁₅₀, and −5.8% for D₉₀; −2.0% for rectum V₁₀₀, and −5.1% for D_0.1; −5.0% for urethra D₁₀, and −5.1% for D₃₀. The same differences between our water and prostate results were all under 0.3%. Uncertainties estimations were up to 2.9% for the g_L(r) function, 13.4% for the F(r,θ) function and 7.0% for Λ, mainly due to seed geometry uncertainties. Uncertainties in extracting the TG-43U1 parameters in the MC simulations as well as in the literature comparison are of the same order of magnitude as the differences between dose distributions computed for water and prostate-like medium. The selection of the parameters for the PS should be done carefully, as it may considerably affect the dose distributions. The seeds internal geometry uncertainties are a major limiting factor in the MC parameters deduction.     

Commissioning and validation of commercial deformable image registration software for adaptive contouring

PurposeTo report the commissioning and validation of deformable image registration(DIR) software for adaptive contouring.MethodsDIR (SmartAdapt®v13.6) was validated using two methods namely contour propagation accuracy and landmark tracking, using physical phantoms and clinical images of various disease sites. Five in-house made phantoms with various known deformations and a set of 10 virtual phantoms were used. Displacement in lateral, anterio-posterior (AP) and superior-inferior (SI) direction were evaluated for various organs and compared with the ground truth. Four clinical sites namely, brain (n = 5), HN (n = 9), cervix (n = 18) and prostate (n = 23) were used. Organs were manually delineated by a radiation oncologist, compared with the deformable image registration (DIR) generated contours. 3D slicer v4.5.0.1 was used to analyze Dice Similarity Co-efficient (DSC), shift in centre of mass (COM) and Hausdorff distances Hf_95%/avg.ResultsMean (SD) DSC, Hf_95% (mm), Hf_avg (mm) and COM of all the phantoms 1–5 were 0.84 (0.2) mm, 5.1 (7.4) mm, 1.6 (2.2) mm, and 1.6 (0.2) mm respectively. Phantom-5 had the largest deformation as compared to phantoms 1–4, and hence had suboptimal indices. The virtual phantom resulted in consistent results for all the ROIs investigated. Contours propagated for brain patients were better with a high DSC score (0.91 (0.04)) as compared to other sites (HN: 0.84, prostate: 0.81 and cervix 0.77). A similar trend was seen in other indices too. The accuracy of propagated contours is limited for complex deformations that include large volume and shape change of bladder and rectum respectively. Visual validation of the propagated contours is recommended for clinical implementation.ConclusionThe DIR algorithm was commissioned and validated for adaptive contouring.     

A new cone-beam computed tomography dosimetry method providing optimal measurement positions: A Monte Carlo study

PurposeThe conventional weighted computed tomography dose index (CTDI_w) may not be suitable for cone-beam computed tomography (CBCT) dosimetry because a cross-sectional dose distribution is angularly inhomogeneous owing to partial angle irradiations. This study was conducted to develop a new dose metric (f(0)^CB_w) for CBCT dosimetry to determine a more accurate average dose in the central cross-sectional plane of a cylindrical phantom using Monte Carlo simulations.MethodsFirst, cross-sectional dose distributions of cylindrical polymethyl methacrylate phantoms over a wide range of phantom diameters (8–40 cm) were calculated for various CBCT scan protocols. Then, by obtaining linear least-squares fits of the full datasets of the cross-sectional dose distributions, the optimal radial positions, which represented measurement positions for the average phantom dose, were determined. Finally, the f(0)^CB_w method was developed by averaging point doses at the optimal radial positions of the phantoms. To demonstrate its validity, the relative differences between the average doses and each dose index value were estimated for the devised f(0)^CB_w, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ResultsThe relative differences between the average doses and each dose index value were within 4.1%, 16.7%, and 11.9% for the devised, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ConclusionsThe devised f(0)^CB_w value was calculated by averaging four “point doses” at 90° intervals and the optimal radial positions of the cylindrical phantom. The devised method can estimate the average dose more accurately than the previously developed CTDI_w methods for CBCT dosimetry.     

Monte Carlo dose in a prosthesis phantom based on exact geometry vs streak artefact contaminated CT data as benchmarked against Gafchromic film measurements

PurposeIn radiotherapy, accurate calculation of patient radiation dose is very important for good clinical outcome. In the presence of metallic implants, the dose calculation accuracy could be compromised by metal artefacts generated in computed tomography (CT) images of patients. This study investigates the influence of metal-induced CT artefacts on MC dose calculations in a pelvic prosthesis phantom.MethodsA pelvic phantom containing unilateral Ti prosthesis was CT-scanned and accurate Hounsfield unit (HU) values were assigned to known materials of the phantom as opposed to HU values produced through the artefact CT images of the phantom. Using the DOSXYZnrc MC code, dose calculations were computed in the phantom model constructed from the original CT images containing the artefacts and artefact-free images made from the exact geometry of the phantom with known materials. The dose calculations were benchmarked against Gafchromic EBT3 film measurements using 15 MeV electron and 10 MV photon beams.ResultsThe average deviations between film and MC dose data decreased from 3 ± 2% to 1 ± 1% and from about 6 ± 2% to 3 ± 1% for the artefact and artefact-free phantom models against film data for the electron and photon fields, respectively.ConclusionsFor the Ti prosthesis phantom, the presence of metal-induced CT artefacts could cause dose inaccuracies of about 3%. Construction of an artefact-free phantom model made from the exact geometry of the phantom with known materials to overcome the effect of artefacts is advantageous compared to using CT data directly of which the exact tissue composition is not well-known.     

Development of a deformable phantom for experimental verification of 4D Monte Carlo simulations in a deforming anatomy

PurposeTo verify the accuracy of 4D Monte Carlo (MC) simulations, using the 4DdefDOSXYZnrc user code, in a deforming anatomy. We developed a tissue-equivalent and reproducible deformable lung phantom and evaluated 4D simulations of delivered dose to the phantom by comparing calculations against measurements.MethodsA novel deformable phantom consisting of flexible foam, emulating lung tissue, inside a Lucite external body was constructed. A removable plug, containing an elastic tumor that can hold film and other dosimeters, was inserted in the phantom. Point dose and position measurements were performed inside and outside the tumor using RADPOS 4D dosimetry system. The phantom was irradiated on an Elekta Infinity linac in both stationary and moving states. The dose delivery was simulated using delivery log files and the phantom motion recorded with RADPOS.ResultsReproducibility of the phantom motion was determined to be within 1 mm. The phantom motion presented realistic features like hysteresis. MC calculations and measurements agreed within 2% at the center of tumor. Outside the tumor agreements were better than 5% which were within the positional/dose reading uncertainties at the measurement points. More than 94% of dose points from MC simulations agreed within 2%/2 mm compared to film measurements.ConclusionThe deformable lung phantom presented realistic and reproducible motion characteristics and its use for verification of 4D dose calculations was demonstrated. Our 4DMC method is capable of accurate calculations of the realistic dose delivered to a moving and deforming anatomy during static and dynamic beam delivery techniques.     

Advanced optimization methods for whole pelvic and local prostate external beam therapy

PurposeRadiation treatment planning inherently involves multiple conflicting planning goals, which makes it a suitable application for multicriteria optimization (MCO). This study investigates a MCO algorithm for VMAT planning (VMAT–MCO) for prostate cancer treatments including pelvic lymph nodes and uses standard inverse VMAT optimization (sVMAT) and Tomotherapy planning as benchmarks.MethodsFor each of ten prostate cancer patients, a two stage plan was generated, consisting of a stage 1 plan delivering 22 Gy to the prostate, and a stage 2 plan delivering 50.4 Gy to the lymph nodes and 56 Gy to the prostate with a simultaneous integrated boost. The single plans were generated by three planning techniques (VMAT–MCO, sVMAT, Tomotherapy) and subsequently compared with respect to plan quality and planning time efficiency.ResultsPlan quality was similar for all techniques, but sVMAT showed slightly better rectum (on average D_mean −7%) and bowel sparing (D_mean −17%) compared to VMAT–MCO in the whole pelvic treatments. Tomotherapy plans exhibited higher bladder dose (D_mean +42%) in stage 1 and lower rectum dose (D_mean −6%) in stage 2 than VMAT–MCO. Compared to manual planning, the planning time with MCO was reduced up to 12 and 38 min for stage 1 and 2 plans, respectively.ConclusionMCO can generate highly conformal prostate VMAT plans with minimal workload in the settings of prostate-only treatments and prostate plus lymph nodes irradiation. In the whole pelvic plan manual VMAT optimization led to slightly improved OAR sparing over VMAT–MCO, whereas for the primary prostate treatment plan quality was equal.     

Generating and using patient-specific whole-body models for organ dose estimates in CT with increased accuracy: Feasibility and validation

The estimation of patient dose using Monte Carlo (MC) simulations based on the available patient CT images is limited to the length of the scan. Software tools for dose estimation based on standard computational phantoms overcome this problem; however, they are limited with respect to taking individual patient anatomy into account. The purpose of this study was to generate whole-body patient models in order to take scattered radiation and over-scanning effects into account. Thorax examinations were performed on three physical anthropomorphic phantoms at tube voltages of 80 kV and 120 kV; absorbed dose was measured using thermoluminescence dosimeters (TLD). Whole-body voxel models were built as a combination of the acquired CT images appended by data taken from widely used anthropomorphic voxel phantoms. MC simulations were performed both for the CT image volumes alone and for the whole-body models. Measured and calculated dose distributions were compared for each TLD chip position; additionally, organ doses were determined.MC simulations based only on CT data underestimated dose by 8%–15% on average depending on patient size with highest underestimation values of 37% for the adult phantom at the caudal border of the image volume. The use of whole-body models substantially reduced these errors; measured and simulated results consistently agreed to better than 10%.This study demonstrates that combined whole-body models can provide three-dimensional dose distributions with improved accuracy. Using the presented concept should be of high interest for research studies which demand high accuracy, e.g. for dose optimization efforts.     

Machine log file-based dose verification using novel iterative CBCT reconstruction algorithm in commercial software during volumetric modulated arc therapy for prostate cancer patients

PurposeTo evaluate the utility of the use of iterative cone-beam computed tomography (CBCT) for machine log file-based dose verification during volumetric modulated arc therapy (VMAT) for prostate cancer patients.MethodsAll CBCT acquisition data were used to reconstruct images with the Feldkamp-Davis-Kress algorithm (FDK-CBCT) and the novel iterative algorithm (iCBCT). The Hounsfield unit (HU)-electron density curves for CBCT images were created using the Advanced Electron Density Phantom. The I’mRT and anthropomorphic phantoms were irradiated with VMAT after CBCT registration. Subsequently, fourteen prostate cancer patients received VMAT after CBCT registration. Machine log files and both CBCT images were exported to the PerFRACTION software, and a 3D patient dose was reconstructed. Mean dose for planning target volume (PTV), the bladder, and rectum and the 3D gamma analysis were evaluated.ResultsFor the phantom studies, the variation of HU values was observed at the central position surrounding the bones in FDK-CBCT. There were almost no changes in the difference of doses at the isocenter between measurement and reconstructed dose for planning CT (pCT), FDK-CBCT, and iCBCT. Mean dose differences of PTV, rectum, and bladder between iCBCT and pCT were approximately 2% lower than those between FDK-CBCT and pCT. For the clinical study, average gamma analysis for 2%/2 mm was 98.22% ± 1.07 and 98.81% ± 1.25% in FDK-CBCT and iCBCT, respectively.ConclusionsA similar machine log file-based dose verification accuracy is obtained for FDK-CBCT and iCBCT during VMAT for prostate cancer patients.     

Measuring the dose in bone for spine stereotactic body radiotherapy

PurposeCurrent quality assurance of radiotherapy involving bony regions generally utilises homogeneous phantoms and dose calculations, ignoring the challenges of heterogeneities with dosimetry problems likely occurring around bone. Anthropomorphic phantoms with synthetic bony materials enable realistic end-to-end testing in clinical scenarios. This work reports on measurements and calculated corrections required to directly report dose in bony materials in the context of comprehensive end-to-end dosimetry audit measurements (63 plans, 6 planning systems).Materials and methodsRadiochromic film and microDiamond measurements were performed in an anthropomorphic spine phantom containing bone equivalent materials. Medium dependent correction factors, k_med, were established using 6 MV and 10 MV Linear Accelerator Monte Carlo simulations to account for the detectors being calibrated in water, but measuring in regions of bony material. Both cortical and trabecular bony material were investigated for verification of dose calculations in dose-to-medium (D_m,m) and dose-to-water (D_w,w) scenarios.ResultsFor D_m,m calculations, modelled correction factors for cortical and trabecular bone in film measurements, and for trabecular bone in microDiamond measurements were 0.875(±0.1%), 0.953(±0.3%) and 0.962(±0.4%), respectively. For D_w,w calculations, the corrections were 0.920(±0.1%), 0.982(±0.3%) and 0.993(±0.4%), respectively. In the audit, application of the correction factors improves the mean agreement between treatment plans and measured microDiamond dose from −2.4%(±3.9%) to 0.4%(±3.7%).ConclusionMonte Carlo simulations provide a method for correcting the dose measured in bony materials allowing more accurate comparison with treatment planning system doses. In verification measurements, algorithm specific correction factors should be applied to account for variations in bony material for calculations based on D_m,m and D_w,w.     

Dosimetric verification of a high dose rate brachytherapy treatment planning system in homogeneous and heterogeneous media

ObjectivesTo verify the dosimetric accuracy of treatment plans in high dose rate (HDR) brachytherapy by using Gafchromic EBT2 film and to demonstrate the adequacy of dose calculations of a commercial treatment planning system (TPS) in a heterogeneous medium.MethodsAbsorbed doses at chosen points in anatomically different tissue equivalent phantoms were measured using Gafchromic EBT2 film. In one case, tandem ovoid brachytherapy was performed in a homogeneous cervix phantom, whereas in the other, organ heterogeneities were introduced in a phantom to replicate the upper thorax for esophageal brachytherapy treatment. A commercially available TPS was used to perform treatment planning in each case and the EBT2 films were irradiated with the HDR Ir-192 brachytherapy source.ResultsFilm measurements in the cervix phantom were found to agree with the TPS calculated values within 3% in the clinically relevant volume. In the thorax phantom, the presence of surrounding heterogeneities was not seen to affect the dose distribution in the volume being treated, whereas, a little dose perturbation was observed at the lung surface. Doses to the spinal cord and to the sternum bone were overestimated and underestimated by 14.6% and 16.5% respectively by the TPS relative to the film measurements. At the trachea wall facing the esophagus, a dose reduction of 10% was noticed in the measurements.ConclusionsThe dose calculation accuracy of the TPS was confirmed in homogeneous medium, whereas, it was proved inadequate to produce correct dosimetric results in conditions of tissue heterogeneity.     

Validation of a commercial TPS based on the VMC++ Monte Carlo code for electron beams: Commissioning and dosimetric comparison with EGSnrc in homogeneous and heterogeneous phantoms

The aim of the present work was the validation of the VMC⁺⁺ Monte Carlo (MC) engine implemented in the Oncentra Masterplan (OMTPS) and used to calculate the dose distribution produced by the electron beams (energy 5-12 MeV) generated by the linear accelerator (linac) Primus (Siemens), shaped by a digital variable applicator (DEVA). The BEAMnrc/DOSXYZnrc (EGSnrc package) MC model of the linac head was used as a benchmark.Commissioning results for both MC codes were evaluated by means of 1D Gamma Analysis (2%, 2 mm), calculated with a home-made Matlab (The MathWorks) program, comparing the calculations with the measured profiles. The results of the commissioning of OMTPS were good [average gamma index (γ) > 97%]; some mismatches were found with large beams (size ≥ 15 cm). The optimization of the BEAMnrc model required to increase the beam exit window to match the calculated and measured profiles (final average γ > 98%).Then OMTPS dose distribution maps were compared with DOSXYZnrc with a 2D Gamma Analysis (3%, 3 mm), in 3 virtual water phantoms: (a) with an air step, (b) with an air insert, and (c) with a bone insert.The OMTPD and EGSnrc dose distributions with the air-water step phantom were in very high agreement (γ ∼ 99%), while for heterogeneous phantoms there were differences of about 9% in the air insert and of about 10–15% in the bone region. This is due to the Masterplan implementation of VMC⁺⁺ which reports the dose as “dose to water”, instead of “dose to medium”.     

Evaluation of PENFAST – A fast Monte Carlo code for dose calculations in photon and electron radiotherapy treatment planning

The aim of the present study is to demonstrate the potential of accelerated dose calculations, using the fast Monte Carlo (MC) code referred to as PENFAST, rather than the conventional MC code PENELOPE, without losing accuracy in the computed dose. For this purpose, experimental measurements of dose distributions in homogeneous and inhomogeneous phantoms were compared with simulated results using both PENELOPE and PENFAST. The simulations and experiments were performed using a Saturne 43 linac operated at 12 MV (photons), and at 18 MeV (electrons). Pre-calculated phase space files (PSFs) were used as input data to both the PENELOPE and PENFAST dose simulations. Since depth–dose and dose profile comparisons between simulations and measurements in water were found to be in good agreement (within ±1% to 1 mm), the PSF calculation is considered to have been validated. In addition, measured dose distributions were compared to simulated results in a set of clinically relevant, inhomogeneous phantoms, consisting of lung and bone heterogeneities in a water tank. In general, the PENFAST results agree to within a 1% to 1 mm difference with those produced by PENELOPE, and to within a 2% to 2 mm difference with measured values. Our study thus provides a pre-clinical validation of the PENFAST code. It also demonstrates that PENFAST provides accurate results for both photon and electron beams, equivalent to those obtained with PENELOPE. CPU time comparisons between both MC codes show that PENFAST is generally about 9–21 times faster than PENELOPE.     

Accuracy of dose calculation on iterative CBCT for head and neck radiotherapy

PurposeTo evaluate the feasibility of the use of iterative cone-beam computed tomography (CBCT) for dose calculation in the head and neck region.MethodsThis study includes phantom and clinical studies. All acquired CBCT images were reconstructed with Feldkamp–Davis–Kress algorithm-based CBCT (FDK-CBCT) and iterative CBCT (iCBCT) algorithm. The Hounsfield unit (HU) consistency between the head and body phantoms was determined in both reconstruction techniques. Volumetric modulated arc therapy (VMAT) plans were generated for 16 head and neck patients on a planning CT scan, and the doses were recalculated on FDK-CBCT and iCBCT with Anisotropic Analytical Algorithm (AAA) and Acuros XB (AXB). As a comparison of the accuracy of dose calculations, the absolute dosimetric difference and 1%/1 mm gamma passing rate analysis were analyzed.ResultsThe difference in the mean HU values between the head and body phantoms was larger for FDK-CBCT (max value: 449.1 HU) than iCBCT (260.0 HU). The median dosimetric difference from the planning CT were <1.0% for both FDK-CBCT and iCBCT but smaller differences were found with iCBCT (planning target volume D_50%: 0.38% (0.15–0.59%) for FDK-CBCT, 0.28% (0.13–0.49%) for iCBCT, AAA; 0.14% (0.04–0.19%) for FDK-CBCT, 0.07% (0.02–0.20%) for iCBCT). The mean gamma passing rate was significantly better in iCBCT than FDK-CBCT (AAA: 98.7% for FDK-CBCT, 99.4% for iCBCT; AXB: 96.8% for FDK_CBCT, 97.5% for iCBCT).ConclusionThe iCBCT-based dose calculation in VMAT for head and neck cancer was accurate compared to FDK-CBCT.     

A Monte Carlo study on dose distribution evaluation of Flexisource 192Ir brachytherapy source

Aim

The aim of this study is to evaluate the dose distribution of the Flexisource ¹⁹²Ir source.

Background

Dosimetric evaluation of brachytherapy sources is recommended by task group number 43 (TG. 43) of American Association of Physicists in Medicine (AAPM).

Materials and methods

MCNPX code was used to simulate Flexisource ¹⁹²Ir source. Dose rate constant and radial dose function were obtained for water and soft tissue phantoms and compared with previous data on this source. Furthermore, dose rate along the transverse axis was obtained by simulation of the Flexisource and a point source and the obtained data were compared with those from Flexiplan treatment planning system (TPS).

Results

The values of dose rate constant obtained for water and soft tissue phantoms were equal to 1.108 and 1.106, respectively. The values of the radial dose function are listed in the form of tabulated data. The values of dose rate (cGy/s) obtained are shown in the form of tabulated data and figures. The maximum difference between TPS and Monte Carlo (MC) dose rate values was 11% in a water phantom at 6.0 cm from the source.

Conclusion

Based on dosimetric parameter comparisons with values previously published, the accuracy of our simulation of Flexisource ¹⁹²Ir was verified. The results of dose rate constant and radial dose function in water and soft tissue phantoms were the same for Flexisource and point sources. For Flexisource ¹⁹²Ir source, the results of TPS calculations in a water phantom were in agreement with the simulations within the calculation uncertainties. Furthermore, the results from the TPS calculation for Flexisource and MC calculation for a point source were practically equal within the calculation uncertainties.     

Phantom design and dosimetric characterization for multiple simultaneous cell irradiations with active pencil beam scanning

A new phantom was designed for in vitro studies on cell lines in horizontal particle beams. The phantom enables simultaneous irradiation at multiple positions along the beam path. The main purpose of this study was the detailed dosimetric characterization of the phantom which consists of various heterogeneous structures. The dosimetric measurements described here were performed under non-reference conditions. The experiment involved a CT scan of the phantom, dose calculations performed with the treatment planning system (TPS) RayStation employing both the Pencil Beam (PB) and Monte Carlo (MC) algorithms, and proton beam delivery. Two treatment plans reflecting the typical target location for head and neck cancer and prostate cancer treatment were created. Absorbed dose to water and dose homogeneity were experimentally assessed within the phantom along the Bragg curve with ionization chambers (ICs) and EBT3 films. LET_d distributions were obtained from the TPS. Measured depth dose distributions were in good agreement with the Monte Carlo-based TPS data. Absorbed dose calculated with the PB algorithm was 4% higher than the absorbed dose measured with ICs at the deepest measurement point along the spread-out Bragg peak. Results of experiments using melanoma (SKMel) cell line are also presented. The study suggested a pronounced correlation between the relative biological effectiveness (RBE) and LET_d, where higher LET_d leads to elevated cell death and cell inactivation. Obtained RBE values ranged from 1.4 to 1.8 at the survival level of 10% (RBE₁₀). It is concluded that dosimetric characterization of a phantom before its use for RBE experiments is essential, since a high dosimetric accuracy contributes to reliable RBE data and allows for a clearer differentiation between physical and biological uncertainties.

     

Experimental benchmarking of RayStation proton dose calculation algorithms inside and outside the target region in heterogeneous phantom geometries

PurposeThe aim of the presented study was to complement existing literature on benchmarking proton dose by comparing dose calculations with experimental measurements in heterogeneous phantom. Points of interest inside and outside the target were considered to quantify the magnitude of calculation uncertainties in current and previous proton therapy practice that might especially have an impact on the dose in organs at risk (OARs).MethodsThe RayStation treatment planning system (RaySearch Laboratories), offering two dose calculation algorithms for pencil beam scanning in proton therapy, i.e., Pencil Beam (PB) and Monte Carlo (MC), was utilized. Treatment plans for a target located behind the interface of the heterogeneous tissues were generated. Dose measurements within and behind the target were performed in a water phantom with embedded slabs of various tissue equivalent materials and 24 PinPoint ionization chambers (PTW). In total 12 test configurations encompassing two different target depths, oblique beam incidence of 30 degrees and range shifter, were considered.ResultsPB and MC calculated doses agreed equally well with the measurements for all test geometries within the target, including the range shifter (mean dose differences ± 3%). Outside the target, the maximum dose difference of 9% (19%) was observed for MC (PB) for the oblique beam incidence and inserted range shifter.ConclusionThe accuracy of MC dose algorithm was superior compared to the PB algorithm, especially outside the target volumes. MC based dose calculation should therefore be preferred in treatment scenarios with heterogeneities, especially to reduce clinically relevant uncertainties for OARs.     

Validation of irtGPUMCD,a GPU-based Monte Carlo internal dosimetry framework for radionuclide therapy

PurposeMonte Carlo (MC) simulations are highly desirable for dose treatment planning and evaluation in radiation oncology. This is true also in emerging nuclear medicine applications such as internal radiotherapy with radionuclides. The purpose of this study is the validation of irtGPUMCD, a GPU-based MC code for dose calculations in internal radiotherapy.MethodsThe female and male phantoms of the International Commission on Radiological Protection (ICRP 110) were used as benchmarking geometries for this study focused on ¹⁷⁷Lu and including ^99mTc and ¹³¹I. Dose calculations were also conducted for a real patient. For phantoms, twelve anatomical structures were considered as target/source organs. The S-values were evaluated with irtGPUMCD simulations (10⁸ photons), with gamma branching ratios of ICRP 107 publication. The ¹⁷⁷Lu electrons S-values were calculated for source organs only, based on local deposition of dose in irtGPUMCD. The S-value relative difference between irtGPUMCD and IDAC-DOSE were evaluated for all targets/sources considered. A DVHs comparison with GATE was conducted. An exponential track length estimator was introduced in irtGPUMCD to increase computational efficiency.ResultsThe relative S-value differences between irtGPUMCD and IDAC-DOSE were <5% while this comparison with GATE was <1%. The DVHs dosimetric indices comparison between GATE and irtGPUMCD for the patient led to an excellent agreement (<2%). The time required for the simulation of 10⁸ photons was 1.5 min for the female phantom, and one minute for the real patient (<1% uncertainty). These results are promising and let envision the use of irtGPUMCD for internal dosimetry in clinical applications.     

Loose versus stranded seeds in permanent prostate brachytherapy: Dosimetric comparison of intraoperative plans

PurposeTo compare target volume coverage and critical organ dosimetry of intraoperative treatment plans for loose seed (LS) and stranded seed (SS) ¹²⁵I permanent implants for low and intermediate risk prostate cancer.MethodsTwo hundred and five patients who underwent permanent seed brachytherapy were included in the study. For prostate dosimetry V90, V100, V150, V200, D90 and COIN were used. The dose to urethra and rectum was determined by the maximal dose and relative doses that cover specified volumes. Means and standard deviations were calculated and statistically compared.ResultsOn average, 54 (range, 30–78) and 48 (range, 31–67) seeds were implanted in the prostate with individual median activities of 0.62 U (range, 0.52–0.70 U) and 0.71 U (range, 0.65–0.71 U) for LS and SS technique, respectively. The target coverage was slightly better with SS (V100: 98% vs. 96%, p < 0.05; D90: 172 Gy vs. 166 Gy, p < 0.05), but more conformal dose distributions were observed with LS (COIN: 0.70 vs. 0.63, p < 0.05). The dose homogeneity did not differ significantly between the two groups. Regarding the dose to urethra and rectum all dose parameters were significantly lower with LS.ConclusionsLS resulted in less dose to the urethra and rectum compared to SS in intraoperative dosimetry. A slightly better target volume coverage with decreased conformity of dose distribution is reported with SS. More studies are necessary to determine how these results will affect postoperative dosimetry, and ultimately, clinical outcome.     

Dose calculation accuracy of different image value to density tables for cone-beam CT planning in head & neck and pelvic localizations

PurposeWe aimed to identify the most accurate combination of phantom and protocol for image value to density table (IVDT) on volume-modulated arc therapy (VMAT) dose calculation based on kV-Cone-beam CT imaging, for head and neck (H&N) and pelvic localizations.MethodsThree phantoms (Catphan^®600, CIRS^®062M (inner phantom for head and outer phantom for body), and TomoTherapy^® “Cheese” phantom) were used to create IVDT curves of CBCT systems with two different CBCT protocols (Standard-dose Head and Standard Pelvis). Hounsfield Unit (HU) time stability and repeatability for a single On-Board-Imager (OBI) and compatibility of two distinct devices were assessed with Catphan^®600. Images from the anthropomorphic phantom CIRS ATOM^® for both CT and CBCT modalities were used for VMAT dose calculation from different IVDT curves. Dosimetric indices from CT and CBCT imaging were compared.ResultsIVDT curves from CBCT images were highly different depending on phantom used (up to 1000 HU for high densities) and protocol applied (up to 200 HU for high densities). HU time stability was verified over seven weeks. A maximum difference of 3% on the dose calculation indices studied was found between CT and CBCT VMAT dose calculation across the two localizations using appropriate IVDT curves. One IVDT curve per localization can be established with a bi-monthly verification of IVDT-CBCT.ConclusionsThe IVDT-CBCT_CIRS-Head phantom with the Standard-dose Head protocol was the most accurate combination for dose calculation on H&N CBCT images. For pelvic localizations, the IVDT-CBCT_Cheese established with the Standard Pelvis protocol provided the best accuracy.