Dosimetric performance of the Elekta Unity MR-linac system: 2D and 3D dosimetry in anthropomorphic inhomogeneous geometry.

Following the clinical introduction of the Elekta Unity MR-linac, there is an urgent need for development of dosimetry protocols and tools, not affected by the presence of a magnetic field. This work presents a benchmarking methodology comprising 2D/3D passive dosimetry and involving on-couch adaptive treatment planning, a unique step in MR-linac workflows. Two identical commercially available 3D-printed head phantoms (featuring realistic bone anatomy and MR/CT contrast) were employed. One phantom incorporated a film dosimetry insert, while the second was filled with polymer gel. Gel dose-response characteristics were evaluated under the Unity irradiation and read-out conditions, using vials and a cubic container filled with gel from the same batch. Treatment plan for the head phantoms involved a hypothetical large C-shape brain lesion, partly surrounding the brainstem. An IMRT step-and-shoot 7-beam plan was employed. Pre-treatment on-couch MR-images were acquired in order for the treatment planning system to calculate the virtual couch shifts and perform adaptive planning. Absolute 2D and relative 3D measurements were compared against calculations related to both adapted and original plans. Real-time dose accumulation monitoring in the gel-filled phantom was also performed. Results from the vials and cubic container suggest that gel dose-response is linear in the dose range investigated and signal integrity is mature at the read-out timings considered. Head phantom 2D and 3D measurements agreed well with calculations with 3D gamma index passing rates above 90% in all cases, even with the most stringent criteria used (2 mm/2%). By exploiting the 3D information provided by the gel, comparison also involved DVHs, dose-volume and plan quality metrics, which also reflected the agreement between adapted and delivered plans within ±4%. No considerable discrepancies were detected between adapted and original plans. A novel methodology was developed and implemented, suitable for QA procedures in Unity. TPS calculations were validated within the experimental uncertainties involved.

Characterization of a novel 3D printed patient specific phantom for quality assurance in cranial stereotactic radiosurgery applications.

In single-isocenter stereotactic radiosurgery/radiotherapy (SRS/SRT) intracranial applications, multiple targets are being treated concurrently, often involving non-coplanar arcs, small photon beams and steep dose gradients. In search for more rigorous quality assurance protocols, this work presents and evaluates a novel methodology for patient-specific pre-treatment plan verification, utilizing 3D printing technology. In a patient's planning CT scan, the external contour and bone structures were segmented and 3D-printed using high-density bone-mimicking material. The resulting head phantom was filled with water while a film dosimetry insert was incorporated. Patient and phantom CT image series were fused and inspected for anatomical coherence. HUs and corresponding densities were compared in several anatomical regions within the head. Furthermore, the level of patient-to-phantom dosimetric equivalence was evaluated both computationally and experimentally. A single-isocenter multi-focal SRS treatment plan was prepared, while dose distributions were calculated on both CT image series, using identical calculation parameters. Phantom- and patient-derived dose distributions were compared in terms of isolines, DVHs, dose-volume metrics and 3D gamma index (GI) analysis. The phantom was treated as if the real patient and film measurements were compared against the patient-derived calculated dose distribution. Visual inspection of the fused CT images suggests excellent geometric similarity between phantom and patient, also confirmed using similarity indices. HUs and densities agreed within one standard deviation except for the skin (modeled as 'bone') and sinuses (water-filled). GI comparison between the calculated distributions resulted in passing rates better than 97% (1%/1 mm). DVHs and dose-volume metrics were also in satisfying agreement. In addition to serving as a feasibility proof-of-concept, experimental absolute film dosimetry verified the computational study results. GI passing rates were above 90%. Results of this work suggest that employing the presented methodology, patient-equivalent phantoms (except for the skin and sinuses areas) can be produced, enabling literally patient-specific pre-treatment plan verification in intracranial applications.

A Holistic Approach to Assessment of Population Exposure to Radiation: Challenges and Initiatives of a Regulatory Authority.

A regulatory authority for radiation safety should continuously evaluate and improve the national safety framework, in line with current requirements and standards. In this context, the Greek Atomic Energy Commission initiated a series of concerted actions. The radiation dose to the population due to public and medical exposures was assessed. The assessment of dose due to public exposure was based on measurements of radon concentrations in dwellings, radionuclide concentrations in environmental samples, and air dose rates; the assessment of dose due to medical exposure was based on dose measurements for typical examinations or procedures and data on their frequency. The mean effective dose to a member of the population was found to be 4.5 mSv (1.8 mSv and 2.7 mSv from medical and public exposures, respectively). Regarding occupational exposure, aircrew dose assessment, eye lens monitoring, and the national dose registry were significantly improved. With respect to artificial tanning (sun beds), the ultraviolet radiation produced was assessed and the practices followed were observed. Results demonstrated exceedance of the 0.3 W m erythema effective irradiance limit set in European Union standards by 63.5% of the sun beds measured, along with general noncompliance with standards. An overarching activity was the upgrade of the Greek Atomic Energy Commission information system in order to collect and disseminate radiation data electronically, launch a networking strategy for interaction with stakeholders, and facilitate the process of regulatory control. In response to the above findings, regulatory actions have been initiated.

On the total system error of a robotic radiosurgery system: phantom measurements, clinical evaluation and long-term analysis.

The total system error (TSE) of a CyberKnife® system was measured using two phantom-based methods and one patient-based method. The standard radiochromic film (RCF) end-to-end (E2E) test using an anthropomorphic head and neck phantom and isocentric treatment delivery was used with the 6Dskull, Fiducial and Xsight® spine (XST) tracking methods. More than 200 RCF-based E2E results covering the period from installation in 2006 until 2017 were analyzed with respect to tracking method, system hardware and software versions, secondary collimation system, and years since installation. An independent polymer gel E2E method was also applied, involving a 3D printed head phantom and multiple spherical target volumes widely distributed within the brain. Finally, the TSE was assessed by comparing the delineated target in the planning computed tomography images of a patient treated for a thalamic functional target with the radiation-induced lesion defined on the six-month follow-up magnetic resonance (MR) images. Statistical analysis of the RCF-based TSE results showed mean ± standard deviation values of 0.40 ± 0.18 mm, 0.40 ± 0.19 mm, and 0.55 ± 0.20 mm for the 6Dskull, Fiducial, and XST tracking methods, respectively. Polymer gel TSE values smaller than 0.66 mm were found for seven targets distributed within the brain, showing that the targeting accuracy of the system is sustained even for targets situated up to 80 mm away from the center of the skull. An average clinical TSE value of 0.87 ± 0.25 mm was also measured using the FSE T2 and FLAIR post-treatment MR image data. Analysis of the long-term RCF-based E2E tests showed no changes of TSE over time. This study is the first to report long-term (>10 years) analysis of TSE, TSE measurement for targets positioned at large distances from the virtual machine isocenter, or a clinical assessment of TSE for the CyberKnife system. All of these measurements demonstrate TSE consistently < 1 mm.

A retrospective dosimetric comparison of TG43 and a commercially available MBDCA for an APBI brachytherapy patient cohort.

PURPOSE: To compare dosimetry using a contemporary model based dose calculation algorithm (MBDCA) following TG186 recommendations, and the conventional TG43 method in an (192)Ir high dose rate (HDR) accelerated partial breast irradiation (APBI) patient cohort. METHODS: Data of 38 APBI patients were studied. Dosimetry for the treatment plans was performed using both the TG43 and TG186 dose calculation methods of the Oncentra Brachy v4.4 treatment planning system (TPS). Analysis included indices of clinical interest for the planning target volume (PTV coverage, dose homogeneity, conformity) as well as dose volume histograms (DVH) for the breast, lung, heart, rib and skin. Significance testing of observed differences between TG43 and TG186 results was carried out and the effect of target location to these differences was studied. RESULTS: Statistically significant differences were observed in the values of clinically relevant DVH parameters for the PTV and the organs at risk (OAR), except for the heart. Differences for the PTV are relatively small (<1% for coverage, on the order of 2% for homogeneity and conformity) with a slight TG43 overestimation except for the dose homogeneity. Percentage differences are larger for the rib and lung (on the order of 4% for Dmax and 5% for V10Gy, respectively) and maximum for the skin (on the order of 6% for D10cc), with a correlation of the observed differences with target location. CONCLUSION: While the MBDCA option of the TPS appears to improve dosimetric accuracy, differences from TG43 do not appear to warrant dose prescription changes or treatment protocol amendment..

Brachytherapy structural shielding calculations using Monte Carlo generated, monoenergetic data.

PURPOSE: To provide a method for calculating the transmission of any broad photon beam with a known energy spectrum in the range of 20-1090 keV, through concrete and lead, based on the superposition of corresponding monoenergetic data obtained from Monte Carlo simulation. METHODS: MCNP5 was used to calculate broad photon beam transmission data through varying thickness of lead and concrete, for monoenergetic point sources of energy in the range pertinent to brachytherapy (20-1090 keV, in 10 keV intervals). The three parameter empirical model introduced by Archer et al. ["Diagnostic x-ray shielding design based on an empirical model of photon attenuation," Health Phys. 44, 507-517 (1983)] was used to describe the transmission curve for each of the 216 energy-material combinations. These three parameters, and hence the transmission curve, for any polyenergetic spectrum can then be obtained by superposition along the lines of Kharrati et al. ["Monte Carlo simulation of x-ray buildup factors of lead and its applications in shielding of diagnostic x-ray facilities," Med. Phys. 34, 1398-1404 (2007)]. A simple program, incorporating a graphical user interface, was developed to facilitate the superposition of monoenergetic data, the graphical and tabular display of broad photon beam transmission curves, and the calculation of material thickness required for a given transmission from these curves. RESULTS: Polyenergetic broad photon beam transmission curves of this work, calculated from the superposition of monoenergetic data, are compared to corresponding results in the literature. A good agreement is observed with results in the literature obtained from Monte Carlo simulations for the photon spectra emitted from bare point sources of various radionuclides. Differences are observed with corresponding results in the literature for x-ray spectra at various tube potentials, mainly due to the different broad beam conditions or x-ray spectra assumed. CONCLUSIONS: The data of this work allow for the accurate calculation of structural shielding thickness, taking into account the spectral variation with shield thickness, and broad beam conditions, in a realistic geometry. The simplicity of calculations also obviates the need for the use of crude transmission data estimates such as the half and tenth value layer indices. Although this study was primarily designed for brachytherapy, results might also be useful for radiology and nuclear medicine facility design, provided broad beam conditions apply.

Dosimetric accuracy of a deterministic radiation transport based (192)Ir brachytherapy treatment planning system. Part III. Comparison to Monte Carlo simulation in voxelized anatomical computational models.

PURPOSE: To compare TG43-based and Acuros deterministic radiation transport-based calculations of the BrachyVision treatment planning system (TPS) with corresponding Monte Carlo (MC) simulation results in heterogeneous patient geometries, in order to validate Acuros and quantify the accuracy improvement it marks relative to TG43. METHODS: Dosimetric comparisons in the form of isodose lines, percentage dose difference maps, and dose volume histogram results were performed for two voxelized mathematical models resembling an esophageal and a breast brachytherapy patient, as well as an actual breast brachytherapy patient model. The mathematical models were converted to digital imaging and communications in medicine (DICOM) image series for input to the TPS. The MCNP5 v.1.40 general-purpose simulation code input files for each model were prepared using information derived from the corresponding DICOM RT exports from the TPS. RESULTS: Comparisons of MC and TG43 results in all models showed significant differences, as reported previously in the literature and expected from the inability of the TG43 based algorithm to account for heterogeneities and model specific scatter conditions. A close agreement was observed between MC and Acuros results in all models except for a limited number of points that lay in the penumbra of perfectly shaped structures in the esophageal model, or at distances very close to the catheters in all models. CONCLUSIONS: Acuros marks a significant dosimetry improvement relative to TG43. The assessment of the clinical significance of this accuracy improvement requires further work. Mathematical patient equivalent models and models prepared from actual patient CT series are useful complementary tools in the methodology outlined in this series of works for the benchmarking of any advanced dose calculation algorithm beyond TG43.

On the output factor measurements of the CyberKnife iris collimator small fields: Experimental determination of the k(Q(clin),Q(msr) ) (f(clin),f(msr) ) correction factors for microchamber and diode detectors.

PURPOSE: To measure the output factors (OFs) of the small fields formed by the variable aperture collimator system (iris) of a CyberKnife (CK) robotic radiosurgery system, and determine the k(Q(clin),Q(msr) ) (f(clin),f(msr) ) correction factors for a microchamber and four diode detectors. METHODS: OF measurements were performed using a PTW PinPoint 31014 microchamber, four diode detectors (PTW-60017, -60012, -60008, and the SunNuclear EDGE detector), TLD-100 microcubes, alanine dosimeters, EBT films, and polymer gels for the 5 mm, 7.5 mm, 10 mm, 12.5 mm, and 15 mm iris collimators at 650 mm, 800 mm, and 1000 mm source to detector distance (SDD). The alanine OF measurements were corrected for volume averaging effects using the 3D dose distributions registered in polymer gel dosimeters. k(Q(clin),Q(msr) ) (f(clin),f(msr) ) correction factors for the PinPoint microchamber and the diode dosimeters were calculated through comparison against corresponding polymer gel, EBT, alanine, and TLD results. RESULTS: Experimental OF results are presented for the array of dosimetric systems used. The PinPoint microchamber was found to underestimate small field OFs, and a k(Q(clin),Q(msr) ) (f(clin),f(msr) ) correction factor ranging from 1.127 ± 0.022 (for the 5 mm iris collimator) to 1.004 ± 0.010 (for the 15 mm iris collimator) was determined at the reference SDD of 800 mm. The PinPoint k(Q(clin),Q(msr) ) (f(clin),f(msr) ) correction factor was also found to increase with decreasing SDD; k(Q(clin),Q(msr) ) (f(clin),f(msr) ) values equal to 1.220 ± 0.028 and 1.077 ± 0.016 were obtained for the 5 mm iris collimator at 650 mm and 1000 mm SDD, respectively. On the contrary, diode detectors were found to overestimate small field OFs and a correction factor equal to 0.973 ± 0.006, 0.954 ± 0.006, 0.937 ± 0.007, and 0.964 ± 0.006 was measured for the PTW-60017, -60012, -60008 and the EDGE diode detectors, respectively, for the 5 mm iris collimator at 800 mm SDD. The corresponding correction factors for the 15 mm iris collimator were found equal to 0.997 ± 0.010, 0.994 ± 0.009, 0.988 ± 0.010, and 0.986 ± 0.010, respectively. No correlation of the diode k(Q(clin),Q(msr) ) (f(clin),f(msr) ) correction factors with SDD was observed. CONCLUSIONS: This work demonstrates an experimental procedure for the determination of the k(Q(clin),Q(msr) ) (f(clin),f(msr) ) correction factors required to obtain small field OF results of increased accuracy.

Dosimetric accuracy of a deterministic radiation transport based 192Ir brachytherapy treatment planning system. Part II: Monte Carlo and experimental verification of a multiple source dwell position plan employing a shielded applicator.

PURPOSE: The aim of this work is the dosimetric validation of a deterministic radiation transport based treatment planning system (BRACHYVISION v. 8.8, referred to as TPS in the following) for multiple 192Ir source dwell position brachytherapy applications employing a shielded applicator in homogeneous water geometries. METHODS: TPS calculations for an irradiation plan employing seven VS2000 192Ir high dose rate (HDR) source dwell positions and a partially shielded applicator (GM11004380) were compared to corresponding Monte Carlo (MC) simulation results, as well as experimental results obtained using the VIP polymer gel-magnetic resonance imaging three-dimensional dosimetry method with a custom made phantom. RESULTS: TPS and MC dose distributions were found in agreement which is mainly within +/- 2%. Considerable differences between TPS and MC results (greater than 2%) were observed at points in the penumbra of the shields (i.e., close to the edges of the "shielded" segment of the geometries). These differences were experimentally verified and therefore attributed to the TPS. Apart from these regions, experimental and TPS dose distributions were found in agreement within 2 mm distance to agreement and 5% dose difference criteria. As shown in this work, these results mark a significant improvement relative to dosimetry algorithms that disregard the presence of the shielded applicator since the use of the latter leads to dosimetry errors on the order of 20%-30% at the edge of the "unshielded" segment of the geometry and even 2%-6% at points corresponding to the potential location of the target volume in clinical applications using the applicator (points in the unshielded segment at short distances from the applicator). CONCLUSIONS: Results of this work attest the capability of the TPS to accurately account for the scatter conditions and the increased attenuation involved in HDR brachytherapy applications employing multiple source dwell positions and partially shielded applicators.

Accuracy of multislice CT angiography for the assessment of in-stent restenoses in the iliac arteries at reduced dose: a phantom study.

OBJECTIVE: We investigated the potential of low-dose CT angiography for accurate assessment of in-stent restenoses (ISRs) of the iliac artery. METHOD: A Rando anthropomorphic phantom (Alderson Research Labs, Stanford, CA), custom-made wax simulating hyperplastic tissue and a nitinol stent were used to simulate a patient with clinically relevant iliac artery ISRs. The cylindrical lumen was filled with a solution of iodine contrast medium diluted in saline, representing a patient's blood during CT angiography. The phantom was subjected to standard- and low-dose angiographic exposures using a modern multidetector (MD) CT scanner. The percentage of ISR was determined using the profile along a line normal to the lumen axis on reconstructed images of 2 and 5 mm slice thickness. Percentage ISRs derived using the standard- and low-dose protocols were compared. In a preliminary study, seven patients with stents were subjected to standard- and low-dose MDCT angiography during follow-up. The resulting images were assessed and compared by two experienced radiologists. RESULTS: The accuracy in measuring the percentage ISR was found to be better than 12% for all simulated stenoses. The differences between percentage ISRs measured on images obtained at 120 kVp/160 mAs and 80 kVp/80 mAs were below 6%. Patient image sets acquired using low-exposure factors were judged to be of satisfactory diagnostic quality. The assessment of ISR did not differ significantly between image sets acquired using the standard factors and those acquired using the low-exposure factors, although the mean reduction in patient effective dose was 48%. CONCLUSION: A reduction in exposure factors during MDCT angiography of the iliac artery is possible without affecting the accuracy in the determination of ISRs.

On the implementation of a recently proposed dosimetric formalism to a robotic radiosurgery system.

PURPOSE: The aim of this work is to implement a recently proposed dosimetric formalism for nonstandard fields to the calibration and small field output factor measurement of a robotic stereotactic radiosurgery system. METHODS: Reference dosimetry measurements were performed in the nonstandard, 60 mm diameter machine specific reference (msr) field using a Farmer ion chamber, five other cylindrical chambers with cavity lengths ranging from 16.25 down to 2.7 mm, and alanine dosimeters. Output factor measurements were performed for the 5, 7.5, 10, and 15 mm field sizes using microchambers, diode detectors, alanine dosimeters, TLD microcubes, and EBT Gafchromic films. Measurement correction factors as described in the proposed formalism were calculated for the ion chamber and diode detector output factor measurements based on published Monte Carlo data. Corresponding volume averaging correction factors were calculated for the alanine output factor measurements using 3D dose distributions, measured with polymer gel dosimeters. RESULTS: Farmer chamber and alanine reference dosimetry results were found in close agreement, yielding a correction factor of k(Q(msr),Q)(f(msr),f(ref)) = 0.999 +/- 0.016 for the chamber readings. These results were also found to be in agreement within experimental uncertainties with corresponding results obtained using the shorter cavity length ionization chambers. The mean measured dose values of the latter, however, were found to be consistently greater than that of the Farmer chamber. This finding, combined with an observed inverse relationship between the mean measured dose and chamber cavity length that follows the trend predicted by theoretical volume averaging calculations in the msr field, implies that the Farmer k(Q(msr),Q)(f(msr),f(ref)) correction is greater than unity. Regarding the output factor results, deviations as large as 33% were observed between the different dosimeters used. These deviations were substantially decreased when appropriate correction factors were applied to the measured microchamber, diode, and alanine values. After correction, all diode and microchamber measured output factors agreed within 1.6% with the corresponding alanine measurements, and within 3.1% with the TLD measurements. The weighted mean output factors were 0.681 +/- 0.001, 0.824 +/- 0.001, 0.875 +/- 0.001, and 0.954 +/- 0.001 for the 5, 7.5, 10, and 15 mm beams, respectively. CONCLUSIONS: The comparison of Farmer chamber measurements versus alanine reference dosimetry validates the use of the former for dosimetry in the msr field of this treatment delivery system. The corresponding results of this work obtained using chambers with different cavity lengths, combined with previous literature findings, suggest that a k(Q(msr),Q)(f(msr),f(ref)) Farmer chamber dose response correction factor of 1.01 may improve calibration measurement accuracy when using the proposed dosimetric formalism. The k(Q(msr),Q)(f(msr),f(ref)) correction factor is within 0.5% from unity for ion chambers with cavity lengths less than 10 mm. Substantial improvements in small field output factor measurement accuracy can be obtained when using microchambers and diodes by applying appropriately calculated correction factors to the detector measurements according to the proposed dosimetric formalism, and their routine use is therefore recommended.

Dosimetric accuracy of a deterministic radiation transport based 192Ir brachytherapy treatment planning system. Part I: single sources and bounded homogeneous geometries.

PURPOSE: The aim of this work is to validate a deterministic radiation transport based treatment planning system (TPS) for single 192Ir brachytherapy source dosimetry in homogeneous water geometries. METHODS: TPS results were obtained using the deterministic radiation transport option of a BRACHYVISION v. 8.8 system for three characteristic source designs (VS2000, GMPlus HDR, and GMPlus PDR) with each source either centered in a 15 cm radius spherical water phantom, or positioned at varying distance away from the phantom center. Corresponding MC simulations were performed using the MCNPX code v.2.5.0 and source geometry models prepared using information provided by the manufacturers. RESULTS: Comparison in terms of the AAPM TG-43 dosimetric formalism quantities, as well as dose rate distributions per unit air kerma strength with a spatial resolution of 0.1 cm, yielded close agreement between TPS and MC results for the sources centered in the phantom. Besides some regions close to the source longitudinal axes where discrepancies could be characterized as systematic, overall agreement for all three sources studied is comparable to the statistical (type A) uncertainty of MC simulations (1% at the majority of points in the geometry increasing to 2%-3% at points lying both away from the source center and close to the source longitudinal axis). A corresponding good agreement was also found between TPS and MC results for the sources positioned away from the phantom center. CONCLUSIONS: Results of this work attest the capability of the TPS to accurately account for the scatter conditions regardless of the size or shape of a given geometry of dosimetric interest, and the position of a source within it. This is important since, as shown in the literature and summarized also in this work, these factors could introduce a significant dosimetric effect that is currently ignored in clinical treatment planning. It is concluded that the implementation of the deterministic radiation transport option of the BRACHYVISION v. 8.8 system for 192Ir brachytherapy dosimetry in homogeneous water geometries yields results of comparable accuracy to the golden standard of Monte Carlo simulation, in clinically viable calculation times.