Determination of intrinsic energy dependence of point-like inorganic scintillation detector in brachytherapy.

BACKGROUND: Inorganic scintillation detectors (ISDs) are promising for in vivo dosimetry in brachytherapy (BT). ISDs have fast response, providing time resolved dose rate information, and high sensitivity, attributed to high atomic numbers. However, the conversion of the detector signal to absorbed dose-to-water is highly dependent on the energy spectrum of the incident radiation. This dependence is comprised of absorbed dose energy dependence, obtainable with Monte Carlo (MC) simulation, and the absorbed dose-to-signal conversion efficiency or intrinsic energy dependence requiring measurements. Studies have indicated negligible intrinsic energy dependence of ZnSe:O-based ISDs in Ir-192 BT. A full characterization has not been performed earlier. PURPOSE: This study characterizes the intrinsic energy dependence of ZnSe:O-based ISDs for kV X-ray radiation qualities, with energies relevant for BT. METHODS: Three point-like ISDs made from fiber-coupled cuboid ZnSe:O-based scintillators were calibrated at the Swedish National Metrology Laboratory for ionizing radiation. The calibration was done in terms of air kerma free-in-air, Kair , in 13 X-ray radiation qualities, Q, from 25 to 300 kVp (CCRI 25-250 kV and ISO 4037 N-series), and in terms of absorbed dose to water, Dw , in a Co-60 beam, Q0 . The mean absorbed dose to the ISDs, relative to Kair and Dw , were obtained with the MC code TOPAS (Geant4) using X-ray spectra obtained with SpekPy software and laboratory filtration data and a generic Co-60 source. The intrinsic energy dependence was determined as a function of effective photon energy, E e f f ${E}_{eff}$ , (relative to Co-60). The angular dependence of the ISD signal was measured in a 25 kVp (0.20 mm Al HVL) and 135 kVp beam (0.48 mm Cu HVL), by rotating the ISDs 180° around the fiber's longitudinal axis (perpendicular to the beam). A full 360° was not performed due to setup limitations. The impact of detector design was quantified with MC simulation. RESULTS: Above 30 keV E e f f ${E}_{eff}$ the intrinsic energy dependence varied with less than 5 ± 4% from unity for all detectors (with the uncertainty expressed as the mean of all expanded measurement uncertainties for individual E e f f ${E}_{eff}$ above 30 keV, k = 2). Below 30 keV, it decreased with up to 17% and inter-detector variations of 13% were observed, likely due to differences in detector geometry not captured by the simulations using nominal geometry. In the 25 kVp radiation quality, the ISD signal varied with 24% over a ∼45° rotation. For 135 kVp, the corresponding variation was below 3%. Assuming a 0.05 mm thicker layer of reflective paint around the sensitive volume changed the absorbed dose with 6.3% at the lowest E e f f ${E}_{eff}$ , and with less than 2% at higher energies. CONCLUSION: The study suggests that the ISDs have an intrinsic energy dependence relative to Co-60 lower than 5 ± 4% in radiation qualities with E e f f ${E}_{eff}\ $ > 30 keV. Therefore, they could in principle be calibrated in a Co-60 beam quality and transferred to such radiation qualities with correction factors determined only by the absorbed dose energy dependence obtained from MC simulations. This encourages exploration of the ISDs' applications in intensity modulated BT with Yb-169 or other novel intermediate energy isotopes.

AAPM WGDCAB Report 372: A joint AAPM, ESTRO, ABG, and ABS report on commissioning of model-based dose calculation algorithms in brachytherapy.

The introduction of model-based dose calculation algorithms (MBDCAs) in brachytherapy provides an opportunity for a more accurate dose calculation and opens the possibility for novel, innovative treatment modalities. The joint AAPM, ESTRO, and ABG Task Group 186 (TG-186) report provided guidance to early adopters. However, the commissioning aspect of these algorithms was described only in general terms with no quantitative goals. This report, from the Working Group on Model-Based Dose Calculation Algorithms in Brachytherapy, introduced a field-tested approach to MBDCA commissioning. It is based on a set of well-characterized test cases for which reference Monte Carlo (MC) and vendor-specific MBDCA dose distributions are available in a Digital Imaging and Communications in Medicine-Radiotherapy (DICOM-RT) format to the clinical users. The key elements of the TG-186 commissioning workflow are now described in detail, and quantitative goals are provided. This approach leverages the well-known Brachytherapy Source Registry jointly managed by the AAPM and the Imaging and Radiation Oncology Core (IROC) Houston Quality Assurance Center (with associated links at ESTRO) to provide open access to test cases as well as step-by-step user guides. While the current report is limited to the two most widely commercially available MBDCAs and only for 192 Ir-based afterloading brachytherapy at this time, this report establishes a general framework that can easily be extended to other brachytherapy MBDCAs and brachytherapy sources. The AAPM, ESTRO, ABG, and ABS recommend that clinical medical physicists implement the workflow presented in this report to validate both the basic and the advanced dose calculation features of their commercial MBDCAs. Recommendations are also given to vendors to integrate advanced analysis tools into their brachytherapy treatment planning system to facilitate extensive dose comparisons. The use of the test cases for research and educational purposes is further encouraged.

Dosimetric impact of a robust optimization approach to mitigate effects from rotational uncertainty in prostate intensity-modulated brachytherapy.

BACKGROUND: Intensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT). PURPOSE: We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria. METHODS: A previously suggested prototype for platinum-shielded prostate 169 Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional 192 Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes. RESULTS: We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra ( D 0.1 c c $D_{0.1cc}$ ) and rectum ( D 0.1 c c $D_{0.1cc}$ and D 1 c c $D_{1cc}$ ) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors. CONCLUSIONS: We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.

Technical note: Evaluation of a spatial optimization model for prostate high dose-rate brachytherapy in a clinical treatment planning system.

BACKGROUND: Spatial properties of a dose distribution, such as volumes of contiguous hot spots, are of clinical importance in treatment planning for high dose-rate brachytherapy (HDR BT). We have in an earlier study developed an optimization model that reduces the prevalence of contiguous hot spots by modifying a tentative treatment plan. PURPOSE: The aim of this study is to incorporate the correction of hot spots in a standard inverse planning workflow and to validate the integrated model in a clinical treatment planning system. The spatial function is included in the objective function for the inverse planning, as opposed to in the previous study where it was applied as a separate post-processing step. Our aim is to demonstrate that fine-adjustments of dose distributions, which are often performed manually in today's clinical practice, can be automated. METHODS: A spatial optimization function was introduced in the treatment planning system RayStation (RaySearch Laboratories AB, Stockholm, Sweden) via a research interface. A series of 10 consecutive prostate patients treated with HDR BT was retrospectively replanned with and without the spatial function. RESULTS: Optimization with the spatial function decreased the volume of the largest contiguous hot spot by on average 31%, compared to if the function was not included. The volume receiving at least 200% of the prescription dose decreased by on average 11%. Target coverage, measured as the fractions of the clinical target volume (CTV) and the planning target volume (PTV) receiving at least the prescription dose, was virtually unchanged (less than a percent change for both metrics). Organs-at-risk received comparable or slightly decreased doses if the spatial function was included in the optimization model. CONCLUSIONS: Optimization of spatial properties such as the volume of contiguous hot spots can be integrated in a standard inverse planning workflow for brachytherapy, and need not be conducted as a separate post-processing step.

An audit of high dose-rate prostate brachytherapy treatment planning at six Swedish clinics.

PURPOSE: High dose-rate prostate brachytherapy has been implemented in Sweden in the late 1980s and early 1990s in six clinics using the same schedule: 20 Gy in two fractions combined with 50 Gy in 25 fractions with external beam radiation therapy. Thirty years have passed and during these years, various aspects of the treatment process have developed, such as ultrasound-guided imaging and treatment planning system. An audit was conducted, including a questionnaire and treatment planning, which aimed to gather knowledge about treatment planning methods in Swedish clinics. MATERIAL AND METHODS: A questionnaire and a treatment planning case (non-anatomical images) were sent to six Swedish clinics, in which high-dose-rate prostate brachytherapy is performed. Treatment plans were compared using dosimetric indices and equivalent 2 Gy doses (EQD2). Treatment planning system report was used to compare dwell positions and dwell times. RESULTS: For all the clinics, the planning aim for the target was 10.0 Gy, but the volume to receive the dose differed from 95% to 100%. Dose constraints for organs at risk varied with up to 2 Gy. The dose to 90% of target volume ranged from 10.0 Gy to 11.1 Gy, equivalent to 26.0 Gy EQD2 and 31.3 Gy EQD2, respectively. Dose non-homogeneity ratio differed from 0.18 to 0.32 for clinical target volume (CTV) in treatment plans and conformity index ranged from 0.52 to 0.59 for CTV. CONCLUSIONS: Dose constraints for the organs at risk are showing a larger variation than that reflected in compared treatments plans. In all treatment plans in our audit, at least 10 Gy was administered giving a total treatment of 102 Gy EQD2, which is in the upper part of the prescription doses published in the GEC/ESTRO recommendations.

SEMI-AUTOMATED 3D SEGMENTATION OF PELVIC REGION BONES IN CT VOLUMES FOR THE ANNOTATION OF MACHINE LEARNING DATASETS.

Automatic segmentation of bones in computed tomography (CT) images is used for instance in beam hardening correction algorithms where it improves the accuracy of resulting CT numbers. Of special interest are pelvic bones, which-because of their strong attenuation-affect the accuracy of brachytherapy in this region. This work evaluated the performance of the JJ2016 algorithm with the performance of MK2014v2 and JS2018 algorithms; all these algorithms were developed by authors. Visual comparison, and, in the latter case, also Dice similarity coefficients derived from the ground truth were used. It was found that the 3D-based JJ2016 performed better than the 2D-based MK2014v2, mainly because of the more accurate hole filling that benefitted from information in adjacent slices. The neural network-based JS2018 outperformed both traditional algorithms. It was, however, limited to the resolution of 1283 owing to the limited amount of memory in the graphical processing unit (GPU).

Optimization in treatment planning of high dose-rate brachytherapy - Review and analysis of mathematical models.

Treatment planning in high dose-rate brachytherapy has traditionally been conducted with manual forward planning, but inverse planning is today increasingly used in clinical practice. There is a large variety of proposed optimization models and algorithms to model and solve the treatment planning problem. Two major parts of inverse treatment planning for which mathematical optimization can be used are the decisions about catheter placement and dwell time distributions. Both these problems as well as integrated approaches are included in this review. The proposed models include linear penalty models, dose-volume models, mean-tail dose models, quadratic penalty models, radiobiological models, and multiobjective models. The aim of this survey is twofold: (i) to give a broad overview over mathematical optimization models used for treatment planning of brachytherapy and (ii) to provide mathematical analyses and comparisons between models. New technologies for brachytherapy treatments and methods for treatment planning are also discussed. Of particular interest for future research is a thorough comparison between optimization models and algorithms on the same dataset, and clinical validation of proposed optimization approaches with respect to patient outcome.

Development and assessment of a quality assurance device for radiation field-light field congruence testing in diagnostic radiology.

Purpose: Existing methods for checking the light field-radiation field congruence on x-ray equipment either do not fully meet the conditions of various quality control standards regarding inherent uncertainty requirements or contain subjective steps, further increasing the uncertainty of the end result. The aim of this work was to develop a method to check the light field-radiation field congruence on all x-ray equipment. The result should have a low uncertainty which is accomplished by eliminating most subjective user steps in the method. A secondary aim was to maintain the same level of usability as of comparable methods but still able to store the result. Approach: A new device has been developed where the light field and corresponding radiation field are monitored through measurements of the field edge locations (in total: 2 × 4 edges ). The maximum field size location deviation between light field and radiation field in the new method is constrained by the physical limitations of the sensors used in various versions of the prototype: linear image sensors (LISs) of 25 to 29 mm active sensor length. The LISs were sensitized to x-rays by applying a phosphor strip of Gd 2 O 2 S : Tb covering the light sensor input area. Later prototypes of the completed LIS device also have the option of a Bluetooth (100-m range standard) connection, thus increasing the mobility. Results: The developed device has a special feature of localization a field edge without any prior, subjective, alignment procedure of the user, i.e., the signals produced were processed by software storing the associated field edge profiles, localizing the edges in them, and finally displaying the calculated deviation. The uncertainty in field edge location difference was estimated to be < 0.1 mm ( k = 2 ). The calculated uncertainty is lower than for other, commercially available, methods for light field-radiation field congruence also presented in this work. Conclusions: A method to check the light field-radiation field congruence of x-ray systems was developed to improve the limitations found in existing methods, such as device detector resolution, subjective operator steps, or the lack of storing results for later analysis. The development work overcame several challenges including mathematically describing real-life edges of light and radiation fields, noise reduction of radiation edges, and mapping/quantification of the rarely observed phenomenon of focal spot wandering. The assessment of the method showed that the listed limitations were overcome, and the aims were accomplished. It is therefore believed that the device can improve the work in quality controls of x-ray systems.

A mathematical optimization model for spatial adjustments of dose distributions in high dose-rate brachytherapy.

High dose-rate brachytherapy is a modality of radiation therapy used for cancer treatment, in which the radiation source is placed within the body. The treatment goal is to give a high enough dose to the tumour while sparing nearby healthy tissue and organs (organs-at-risk). The most common criteria for evaluating dose distributions are dosimetric indices. For the tumour, such an index is the portion of the volume that receives at least a specified dose level (e.g. the prescription dose), while for organs-at-risk it is instead the portion of the volume that receives at most a specified dose level. Dosimetric indices are aggregate criteria and do not consider spatial properties of the dose distribution. Further, there are neither any established evaluation criteria for characterizing spatial properties, nor have such properties been studied in the context of mathematical optimization of brachytherapy. Spatial properties are however of clinical relevance and therefore dose plans are sometimes adjusted manually to improve them. We propose an optimization model for reducing the prevalence of contiguous volumes with a too high dose (hot spots) or a too low dose (cold spots) in a tentative dose plan. This model is independent of the process of constructing the tentative plan. We conduct computational experiments with tentative plans obtained both from optimization models and from clinical practice. The objective function considers pairs of dose points and each pair is given a distance-based penalty if the dose is either too high or too low at both dose points. Constraints are included to retain dosimetric indices at acceptable levels. Our model is designed to automate the manual adjustment step in the planning process. In the automatic adjustment step large-scale optimization models are solved. We show reductions of the volumes of the largest hot and cold spots, and the computing times are feasible in clinical practice.

Testing the methodology for a dosimetric end-to-end audit of IMRT/VMAT: results of IAEA multicentre and national studies.

Introduction: Within an International Atomic Energy Agency (IAEA) co-ordinated research project (CRP), a remote end-to-end dosimetric quality audit for intensity modulated radiation therapy (IMRT)/ volumetric arc therapy (VMAT) was developed to verify the radiotherapy chain including imaging, treatment planning and dose delivery. The methodology as well as the results obtained in a multicentre pilot study and national trial runs conducted in close cooperation with dosimetry audit networks (DANs) of IAEA Member States are presented.Material and methods: A solid polystyrene phantom containing a dosimetry insert with an irregular solid water planning target volume (PTV) and organ at risk (OAR) was designed for this audit. The insert can be preloaded with radiochromic film and four thermoluminescent dosimeters (TLDs). For the audit, radiotherapy centres were asked to scan the phantom, contour the structures, create an IMRT/VMAT treatment plan and irradiate the phantom. The dose prescription was to deliver 4 Gy to the PTV in two fractions and to limit the OAR dose to a maximum of 2.8 Gy. The TLD measured doses and film measured dose distributions were compared with the TPS calculations.Results: Sixteen hospitals from 13 countries and 64 hospitals from 6 countries participated in the multicenter pilot study and in the national runs, respectively. The TLD results for the PTV were all within ±5% acceptance limit for the multicentre pilot study, whereas for national runs, 17 participants failed to meet this criterion. All measured doses in the OAR were below the treatment planning constraint. The film analysis identified seven plans in national runs below the 90% passing rate gamma criteria.Conclusion: The results proved that the methodology of the IMRT/VMAT dosimetric end-to-end audit was feasible for its intended purpose, i.e., the phantom design and materials were suitable; the phantom was easy to use and it was robust enough for shipment. Most importantly the audit methodology was capable of identifying suboptimal IMRT/VMAT delivery.

END-TO-END AUDIT: COMPARISON OF TLD AND LITHIUM FORMATE EPR DOSIMETRY.

The aim of this study was to test two different solid state dosimetry systems for the purpose of end-to-end audits of radiotherapy volumetric modulated arc therapy (VMAT) technique; a lithium formate electron paramagnetic resonance system and a lithium fluoride thermoluminescent dosimetry system. As a complement to the solid state systems, ion chamber measurements were performed. A polystyrene phantom with a planning target volume (PTV) and an organ at risk (OAR) structure was scanned using CT. A VMAT dose plan was optimized to deliver 2 Gy to the target volume and to minimize the dose to the OAR. The different detectors were inserted into the phantom and the planned dose distribution was delivered. The measured doses were compared to the treatment planning system (TPS) calculated doses. Good agreement was found between the TPS calculated and the measured doses, well accepted for the dose determinations in remote dosimetry audits of VMAT treatment technique.

A generic TG-186 shielded applicator for commissioning model-based dose calculation algorithms for high-dose-rate 192 Ir brachytherapy.

PURPOSE: A joint working group was created by the American Association of Physicists in Medicine (AAPM), the European Society for Radiotherapy and Oncology (ESTRO), and the Australasian Brachytherapy Group (ABG) with the charge, among others, to develop a set of well-defined test case plans and perform calculations and comparisons with model-based dose calculation algorithms (MBDCAs). Its main goal is to facilitate a smooth transition from the AAPM Task Group No. 43 (TG-43) dose calculation formalism, widely being used in clinical practice for brachytherapy, to the one proposed by Task Group No. 186 (TG-186) for MBDCAs. To do so, in this work a hypothetical, generic high-dose rate (HDR) 192 Ir shielded applicator has been designed and benchmarked. METHODS: A generic HDR 192 Ir shielded applicator was designed based on three commercially available gynecological applicators as well as a virtual cubic water phantom that can be imported into any DICOM-RT compatible treatment planning system (TPS). The absorbed dose distribution around the applicator with the TG-186 192 Ir source located at one dwell position at its center was computed using two commercial TPSs incorporating MBDCAs (Oncentra® Brachy with Advanced Collapsed-cone Engine, ACE™, and BrachyVision ACUROS™) and state-of-the-art Monte Carlo (MC) codes, including ALGEBRA, BrachyDose, egs_brachy, Geant4, MCNP6, and Penelope2008. TPS-based volumetric dose distributions for the previously reported "source centered in water" and "source displaced" test cases, and the new "source centered in applicator" test case, were analyzed here using the MCNP6 dose distribution as a reference. Volumetric dose comparisons of TPS results against results for the other MC codes were also performed. Distributions of local and global dose difference ratios are reported. RESULTS: The local dose differences among MC codes are comparable to the statistical uncertainties of the reference datasets for the "source centered in water" and "source displaced" test cases and for the clinically relevant part of the unshielded volume in the "source centered in applicator" case. Larger local differences appear in the shielded volume or at large distances. Considering clinically relevant regions, global dose differences are smaller than the local ones. The most disadvantageous case for the MBDCAs is the one including the shielded applicator. In this case, ACUROS agrees with MC within [-4.2%, +4.2%] for the majority of voxels (95%) while presenting dose differences within [-0.12%, +0.12%] of the dose at a clinically relevant reference point. For ACE, 95% of the total volume presents differences with respect to MC in the range [-1.7%, +0.4%] of the dose at the reference point. CONCLUSIONS: The combination of the generic source and generic shielded applicator, together with the previously developed test cases and reference datasets (available in the Brachytherapy Source Registry), lay a solid foundation in supporting uniform commissioning procedures and direct comparisons among treatment planning systems for HDR 192 Ir brachytherapy.

Collapsed cone dose calculations for heterogeneous tissues in brachytherapy using primary and scatter separation source data.

BACKGROUND AND OBJECTIVE: Brachytherapy is a form of radiation therapy using sealed radiation sources inserted within or in the vicinity of the tumor of, e.g., gynecological, prostate or head- and neck cancers. Accurate dose calculation is a crucial part of the treatment planning. Several reviews have called for clinical software with model-based algorithms that better take into account the effects of patient individual distribution of tissues, source-channel and shielding attenuation than the commonly employed TG-43 formalism which simply map homogeneous water dose distributions onto the patient. In this paper we give a comprehensive and thorough derivation of such an algorithm based on collapsed cone point-kernel superposition, and describe details of its implementation into a commercial treatment planning system for clinical use. METHODS: A brachytherapy version of the collapsed-cone algorithm using analytical raytraces of the primary photon radiation followed by successive scattering dose calculation for once- and multiply scattered photons is described in detail, including derivation of the corresponding set of recursive equations for energy transport along cone axes/transport lines and the coupling to clinical source modeling. Specific implementation issues for setting up of the calculation grid, handling of intravoxel gradients and voxels partly containing non-patient applicator material are given. RESULTS: Sample runs for two clinical cases are shown, one being a gynecological application with a tungsten-shielded applicator and one a breast implant. These two cases demonstrate the impact of improved dose calculation versus TG-43 formalism. CONCLUSIONS: Use of model-based dose calculation algorithms for brachytherapy taking the three-dimensional treatment geometry into account increases the dosimetric accuracy in planning and follow up of treatments. The comprehensive description and derivations provided gives a rigid background for further clinical, educational and research applications.

The collapsed cone algorithm for (192)Ir dosimetry using phantom-size adaptive multiple-scatter point kernels.

The aim of this work was to investigate how dose distributions calculated with the collapsed cone (CC) algorithm depend on the size of the water phantom used in deriving the point kernel for multiple scatter. A research version of the CC algorithm equipped with a set of selectable point kernels for multiple-scatter dose that had initially been derived in water phantoms of various dimensions was used. The new point kernels were generated using EGSnrc in spherical water phantoms of radii 5 cm, 7.5 cm, 10 cm, 15 cm, 20 cm, 30 cm and 50 cm. Dose distributions derived with CC in water phantoms of different dimensions and in a CT-based clinical breast geometry were compared to Monte Carlo (MC) simulations using the Geant4-based brachytherapy specific MC code Algebra. Agreement with MC within 1% was obtained when the dimensions of the phantom used to derive the multiple-scatter kernel were similar to those of the calculation phantom. Doses are overestimated at phantom edges when kernels are derived in larger phantoms and underestimated when derived in smaller phantoms (by around 2% to 7% depending on distance from source and phantom dimensions). CC agrees well with MC in the high dose region of a breast implant and is superior to TG43 in determining skin doses for all multiple-scatter point kernel sizes. Increased agreement between CC and MC is achieved when the point kernel is comparable to breast dimensions. The investigated approximation in multiple scatter dose depends on the choice of point kernel in relation to phantom size and yields a significant fraction of the total dose only at distances of several centimeters from a source/implant which correspond to volumes of low doses. The current implementation of the CC algorithm utilizes a point kernel derived in a comparatively large (radius 20 cm) water phantom. A fixed point kernel leads to predictable behaviour of the algorithm with the worst case being a source/implant located well within a patient/phantom for which low doses at phantom edges can be overestimated by 2-5 %. It would be possible to improve the situation by using a point kernel for multiple-scatter dose adapted to the patient/phantom dimensions at hand.

Dose specification for ¹9²Ir high dose rate brachytherapy in terms of dose-to-water-in-medium and dose-to-medium-in-medium.

Dose calculation in high dose rate brachytherapy with (192)Ir is usually based on the TG-43U1 protocol where all media are considered to be water. Several dose calculation algorithms have been developed that are capable of handling heterogeneities with two possibilities to report dose: dose-to-medium-in-medium (Dm,m) and dose-to-water-in-medium (Dw,m). The relation between Dm,m and Dw,m for (192)Ir is the main goal of this study, in particular the dependence of Dw,m on the dose calculation approach using either large cavity theory (LCT) or small cavity theory (SCT). A head and neck case was selected due to the presence of media with a large range of atomic numbers relevant to tissues and mass densities such as air, soft tissues and bone interfaces. This case was simulated using a Monte Carlo (MC) code to score: Dm,m, Dw,m (LCT), mean photon energy and photon fluence. Dw,m (SCT) was derived from MC simulations using the ratio between the unrestricted collisional stopping power of the actual medium and water. Differences between Dm,m and Dw,m (SCT or LCT) can be negligible (<1%) for some tissues e.g. muscle and significant for other tissues with differences of up to 14% for bone. Using SCT or LCT approaches leads to differences between Dw,m (SCT) and Dw,m (LCT) up to 29% for bone and 36% for teeth. The mean photon energy distribution ranges from 222 keV up to 356 keV. However, results obtained using mean photon energies are not equivalent to the ones obtained using the full, local photon spectrum. This work concludes that it is essential that brachytherapy studies clearly report the dose quantity. It further shows that while differences between Dm,m and Dw,m (SCT) mainly depend on tissue type, differences between Dm,m and Dw,m (LCT) are, in addition, significantly dependent on the local photon energy fluence spectrum which varies with distance to implanted sources.

A linear programming model for optimizing HDR brachytherapy dose distributions with respect to mean dose in the DVH-tail.

PURPOSE: Recent research has shown that the optimization model hitherto used in high-dose-rate (HDR) brachytherapy corresponds weakly to the dosimetric indices used to evaluate the quality of a dose distribution. Although alternative models that explicitly include such dosimetric indices have been presented, the inclusion of the dosimetric indices explicitly yields intractable models. The purpose of this paper is to develop a model for optimizing dosimetric indices that is easier to solve than those proposed earlier. METHODS: In this paper, the authors present an alternative approach for optimizing dose distributions for HDR brachytherapy where dosimetric indices are taken into account through surrogates based on the conditional value-at-risk concept. This yields a linear optimization model that is easy to solve, and has the advantage that the constraints are easy to interpret and modify to obtain satisfactory dose distributions. RESULTS: The authors show by experimental comparisons, carried out retrospectively for a set of prostate cancer patients, that their proposed model corresponds well with constraining dosimetric indices. All modifications of the parameters in the authors' model yield the expected result. The dose distributions generated are also comparable to those generated by the standard model with respect to the dosimetric indices that are used for evaluating quality. CONCLUSIONS: The authors' new model is a viable surrogate to optimizing dosimetric indices and quickly and easily yields high quality dose distributions.

Specification of absorbed dose to water using model-based dose calculation algorithms for treatment planning in brachytherapy.

Model-based dose calculation algorithms (MBDCAs), recently introduced in treatment planning systems (TPS) for brachytherapy, calculate tissue absorbed doses. In the TPS framework, doses have hereto been reported as dose to water and water may still be preferred as a dose specification medium. Dose to tissue medium Dmed then needs to be converted into dose to water in tissue Dw,med. Methods to calculate absorbed dose to differently sized water compartments/cavities inside tissue, infinitesimal (used for definition of absorbed dose), small, large or intermediate, are reviewed. Burlin theory is applied to estimate photon energies at which cavity sizes in the range 1 nm-10 mm can be considered small or large. Photon and electron energy spectra are calculated at 1 cm distance from the central axis in cylindrical phantoms of bone, muscle and adipose tissue for 20, 50, 300 keV photons and photons from (125)I, (169)Yb and (192)Ir sources; ratios of mass-collision-stopping powers and mass energy absorption coefficients are calculated as applicable to convert Dmed into Dw,med for small and large cavities. Results show that 1-10 nm sized cavities are small at all investigated photon energies; 100 µm cavities are large only at photon energies <20 keV. A choice of an appropriate conversion coefficient Dw, med/Dmed is discussed in terms of the cavity size in relation to the size of important cellular targets. Free radicals from DNA bound water of nanometre dimensions contribute to DNA damage and cell killing and may be the most important water compartment in cells implying use of ratios of mass-collision-stopping powers for converting Dmed into Dw,med.

Response of LiF:Mg,Ti thermoluminescent dosimeters at photon energies relevant to the dosimetry of brachytherapy (<1 MeV).

PURPOSE: High energy photon beams are used in calibrating dosimeters for use in brachytherapy since absorbed dose to water can be determined accurately and with traceability to primary standards in such beams, using calibrated ion chambers and standard dosimetry protocols. For use in brachytherapy, beam quality correction factors are needed, which include corrections for differences in mass energy absorption properties between water and detector as well as variations in detector response (intrinsic efficiency) with radiation quality, caused by variations in the density of ionization (linear energy transfer (LET) -distributions) along the secondary electron tracks. The aim of this work was to investigate experimentally the detector response of LiF:Mg,Ti thermoluminescent dosimeters (TLD) for photon energies below 1 MeV relative to (60)Co and to address discrepancies between the results found in recent publications of detector response. METHODS: LiF:Mg,Ti dosimeters of formulation MTS-N Poland were irradiated to known values of air kerma free-in-air in x-ray beams at tube voltages 25-250 kV, in (137)Cs- and (60)Co-beams at the Swedish Secondary Standards Dosimetry Laboratory. Conversions from air kerma free-in-air into values of mean absorbed dose in the dosimeters in the actual irradiation geometries were made using EGSnrc Monte Carlo simulations. X-ray energy spectra were measured or calculated for the actual beams. Detector response relative to that for (60)Co was determined at each beam quality. RESULTS: An increase in relative response was seen for all beam qualities ranging from 8% at tube voltage 25 kV (effective energy 13 keV) to 3%-4% at 250 kV (122 keV effective energy) and (137)Cs with a minimum at 80 keV effective energy (tube voltage 180 kV). The variation with effective energy was similar to that reported by Davis et al. [Radiat. Prot. Dosim. 106, 33-43 (2003)] with our values being systematically lower by 2%-4%. Compared to the results by Nunn et al. [Med. Phys. 35, 1861-1869 (2008)], the relative detector response as a function of effective energy differed in both shape and magnitude. This could be explained by the higher maximum read-out temperature (350 °C) used by Nunn et al. [Med. Phys. 35, 1861-1869 (2008)], allowing light emitted from high-temperature peaks with a strong LET dependence to be registered. Use of TLD-100 by Davis et al. [Radiat. Prot. Dosim. 106, 33-43 (2003)] with a stronger super-linear dose response compared to MTS-N was identified as causing the lower relative detector response in this work. CONCLUSIONS: Both careful dosimetry and strict protocols for handling the TLDs are required to reach solid experimental data on relative detector response. This work confirms older findings that an over-response relative to (60)Co exists for photon energies below 200-300 keV. Comparison with the results from the literature indicates that using similar protocols for annealing and read-out, dosimeters of different makes (TLD-100, MTS-N) differ in relative detector response. Though universality of the results has not been proven and further investigation is needed, it is anticipated that with the use of strict protocols for annealing and read-out, it will be possible to determine correction factors that can be used to reduce uncertainties in dose measurements around brachytherapy sources at photon energies where primary standards for absorbed dose to water are not available.

Characterization of a 60Co unit at a secondary standard dosimetry laboratory: Monte Carlo simulations compared to measurements and results from the literature.

PURPOSE: To compare a Monte Carlo (MC) characterization of a 60Co unit at the Swedish Secondary Standard Dosimetry Laboratory (SSDL) with the results of both measurements and literature with the aims of (1) resolving a change in the ratio of air-kerma free in air Kair and absorbed dose to water Dw in a water phantom noted experimentally after a source exchange in the laboratory and (2) reviewing results from the literature on similar MC simulations. Although their use in radiotherapy is decreasing, the characteristics of 60Co beams are of interest since 60Co beams are utilized in calibrating ionization chambers for the absolute dosimetry of radiotherapy beams and as reference radiation quality in evaluating the energy dependence of radiation detectors and in studies on radiobiological effectiveness. METHODS: The BEAMnrc MC code was used with a detailed geometrical model of the treatment head and two models of the 60Co source representing the sources used before and after source exchange, respectively. The active diameters of the 60Co sources were 1.5 cm in pellet form and 2.0 cm in sintered form. Measurements were performed on the actual unit at the Swedish SSDL. RESULTS: Agreement was obtained between the MC and the measured results within the estimated uncertainties for beam profiles, water depth-dose curve, relative air-kerma output factors, and for the ratios of Kair/Dw before and after source exchange. The on-axis energy distribution of the photon fluence free in air for the unit loaded with its present (1.5 cm in diameter) source agreed closely with the results from the literature in which a source of the same make and active diameter, inside a different treatment head, was simulated. The spectrum for the larger (2.0 cm in diameter) source was in close agreement with another published spectrum, also modeling a 60Co source with an active diameter of 2.0 cm inside a different treatment head. CONCLUSIONS: The reduction in the value of Kair/Dw following source exchange was explained by the spectral differences between the two sources that were larger in the free in-air geometry used for Kair calibrations than at 5 g/cm2 depth in the water phantom used for Dw calibrations. Literature review revealed differences between published in-air 60Co spectra derived for sources of different active diameters, and investigators in need of an accurately determined 60Co in-air spectrum should be aware of differences due to source active diameter.