The charge of Task Group 186 (TG-186) is to provide guidance for early adopters of model-based dose calculation algorithms (MBDCAs) for brachytherapy (BT) dose calculations to ensure practice ...uniformity. Contrary to external beam radiotherapy, heterogeneity correction algorithms have only recently been made available to the BT community. Yet, BT dose calculation accuracy is highly dependent on scatter conditions and photoelectric effect cross-sections relative to water. In specific situations, differences between the current water-based BT dose calculation formalism (TG-43) and MBDCAs can lead to differences in calculated doses exceeding a factor of 10. MBDCAs raise three major issues that are not addressed by current guidance documents: (1) MBDCA calculated doses are sensitive to the dose specification medium, resulting in energy-dependent differences between dose calculated to water in a homogeneous water geometry (TG-43), dose calculated to the local medium in the heterogeneous medium, and the intermediate scenario of dose calculated to a small volume of water in the heterogeneous medium. (2) MBDCA doses are sensitive to voxel-by-voxel interaction cross sections. Neither conventional single-energy CT nor ICRU/ICRP tissue composition compilations provide useful guidance for the task of assigning interaction cross sections to each voxel. (3) Since each patient-source-applicator combination is unique, having reference data for each possible combination to benchmark MBDCAs is an impractical strategy. Hence, a new commissioning process is required. TG-186 addresses in detail the above issues through the literature review and provides explicit recommendations based on the current state of knowledge. TG-43-based dose prescription and dose calculation remain in effect, with MBDCA dose reporting performed in parallel when available. In using MBDCAs, it is recommended that the radiation transport should be performed in the heterogeneous medium and, at minimum, the dose to the local medium be reported along with the TG-43 calculated doses. Assignments of voxel-by-voxel cross sections represent a particular challenge. Electron density information is readily extracted from CT imaging, but cannot be used to distinguish between different materials having the same density. Therefore, a recommendation is made to use a number of standardized materials to maintain uniformity across institutions. Sensitivity analysis shows that this recommendation offers increased accuracy over TG-43. MBDCA commissioning will share commonalities with current TG-43-based systems, but in addition there will be algorithm-specific tasks. Two levels of commissioning are recommended: reproducing TG-43 dose parameters and testing the advanced capabilities of MBDCAs. For validation of heterogeneity and scatter conditions, MBDCAs should mimic the 3D dose distributions from reference virtual geometries. Potential changes in BT dose prescriptions and MBDCA limitations are discussed. When data required for full MBDCA implementation are insufficient, interim recommendations are made and potential areas of research are identified. Application of TG-186 guidance should retain practice uniformity in transitioning from the TG-43 to the MBDCA approach.
To quantify differences that exist between dosimetry models used for
Y selective internal radiation therapy (SIRT).
Retrospectively, 37 tumors were delineated on 19 post-therapy quantitative
Y single ...photon emission computed tomography/computed tomography scans. Using matched volumes of interest (VOIs), absorbed doses were reported using 3 dosimetry models: glass microsphere package insert standard model (SM), partition model (PM), and Monte Carlo (MC). Univariate linear regressions were performed to predict mean MC from SM and PM. Analysis was performed for 2 subsets: cases with a single tumor delineated (best case for PM), and cases with multiple tumors delineated (typical clinical scenario). Variability in PM from the ad hoc placement of a single spherical VOI to estimate the entire normal liver activity concentration for tumor (T) to nontumoral liver (NL) ratios (TNR) was investigated. We interpreted the slope of the resulting regression as bias and the 95% prediction interval (95%PI) as uncertainty. MC
represents MC absorbed doses to the NL for the single tumor patient subset; other combinations of calculations follow a similar naming convention.
SM was unable to predict MC
or MC
(p>.12, 95%PI >±177 Gy). However, SM
was able to predict (p<.012) MC
, albeit with large uncertainties; SM
and SM
yielded biases of 0.62 and 0.71, and 95%PI of ±40 and ± 32 Gy, respectively. PM
and PM
predicted (p<2E-6) MC
and MC
with biases of 0.52 and 0.54, and 95%PI of ±38 and ± 111 Gy, respectively. The TNR variability in PM
increased the 95%PI for predicting MC
(bias = 0.46 and 95%PI = ±103 Gy). The TNR variability in PM
modified the bias when predicting MC
(bias = 0.32 and 95%PI = ±110 Gy).
The SM is unable to predict mean MC tumor absorbed dose. The PM is statistically correlated with mean MC, but the resulting uncertainties in predicted MC are large. Large differences observed between dosimetry models for
Y SIRT warrant caution when interpreting published SIRT absorbed doses. To reduce uncertainty, we suggest the entire NL VOI be used for TNR estimates when using PM.
Purpose:
The novel deterministic radiation transport algorithm, Acuros XB (AXB), has shown great potential for accurate heterogeneous dose calculation. However, the clinical impact between AXB and ...other currently used algorithms still needs to be elucidated for translation between these algorithms. The purpose of this study was to investigate the impact of AXB for heterogeneous dose calculation in lung cancer for intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT).
Methods:
The thorax phantom from the Radiological Physics Center (RPC) was used for this study. IMRT and VMAT plans were created for the phantom in the Eclipse 11.0 treatment planning system. Each plan was delivered to the phantom three times using a Varian Clinac iX linear accelerator to ensure reproducibility. Thermoluminescent dosimeters (TLDs) and Gafchromic EBT2 film were placed inside the phantom to measure delivered doses. The measurements were compared with dose calculations from AXB 11.0.21 and the anisotropic analytical algorithm (AAA) 11.0.21. Two dose reporting modes of AXB, dose-to-medium in medium (D
m,m) and dose-to-water in medium (D
w,m), were studied. Point doses, dose profiles, and gamma analysis were used to quantify the agreement between measurements and calculations from both AXB and AAA. The computation times for AAA and AXB were also evaluated.
Results:
For the RPC lung phantom, AAA and AXB dose predictions were found in good agreement to TLD and film measurements for both IMRT and VMAT plans. TLD dose predictions were within 0.4%–4.4% to AXB doses (bothD
m,m and D
w,m); and within 2.5%–6.4% to AAA doses, respectively. For the film comparisons, the gamma indexes (±3%/3 mm criteria) were 94%, 97%, and 98% for AAA, AXB_D
m,m, and AXB_D
w,m, respectively. The differences between AXB and AAA in dose–volume histogram mean doses were within 2% in the planning target volume, lung, heart, and within 5% in the spinal cord. However, differences up to 8% between AXB and AAA were found at lung/soft tissue interface regions for individual IMRT fields. AAA was found to be 5–6 times faster than AXB for IMRT, while AXB was 4–5 times faster than AAA for VMAT plan.
Conclusions:
AXB is satisfactorily accurate for the dose calculation in lung cancer for both IMRT and VMAT plans. The differences between AXB and AAA are generally small except in heterogeneous interface regions. AXBD
w,m and D
m,m calculations are similar inside the soft tissue and lung regions. AXB can benefit lung VMAT plans by both improving accuracy and reducing computation time.
To investigate the dosimetric impact of the heterogeneity dose calculation Acuros (Transpire Inc., Gig Harbor, WA), a grid-based Boltzmann equation solver (GBBS), for brachytherapy in a cohort of ...cervical cancer patients.
The impact of heterogeneities was retrospectively assessed in treatment plans for 26 patients who had previously received (192)Ir intracavitary brachytherapy for cervical cancer with computed tomography (CT)/magnetic resonance-compatible tandems and unshielded colpostats. The GBBS models sources, patient boundaries, applicators, and tissue heterogeneities. Multiple GBBS calculations were performed with and without solid model applicator, with and without overriding the patient contour to 1 g/cm(3) muscle, and with and without overriding contrast materials to muscle or 2.25 g/cm(3) bone. Impact of source and boundary modeling, applicator, tissue heterogeneities, and sensitivity of CT-to-material mapping of contrast were derived from the multiple calculations. American Association of Physicists in Medicine Task Group 43 (TG-43) guidelines and the GBBS were compared for the following clinical dosimetric parameters: Manchester points A and B, International Commission on Radiation Units and Measurements (ICRU) report 38 rectal and bladder points, three and nine o'clock, and (D2cm3) to the bladder, rectum, and sigmoid.
Points A and B, D(2) cm(3) bladder, ICRU bladder, and three and nine o'clock were within 5% of TG-43 for all GBBS calculations. The source and boundary and applicator account for most of the differences between the GBBS and TG-43 guidelines. The D(2cm3) rectum (n = 3), D(2cm3) sigmoid (n = 1), and ICRU rectum (n = 6) had differences of >5% from TG-43 for the worst case incorrect mapping of contrast to bone. Clinical dosimetric parameters were within 5% of TG-43 when rectal and balloon contrast were mapped to bone and radiopaque packing was not overridden.
The GBBS has minimal impact on clinical parameters for this cohort of patients with unshielded applicators. The incorrect mapping of rectal and balloon contrast does not have a significant impact on clinical parameters. Rectal parameters may be sensitive to the mapping of radiopaque packing.
Purpose:
To evaluate the dose distributions of an
source (model VS2000) in homogeneous water geometry calculated using a deterministic grid‐based Boltzmann transport equation solver (GBBS) in the ...commercial treatment planning system (TPS) (
BRACHYVISION‐ACUROS
v8.8).
Methods:
Using percent dose differences
, the GBBS (BV‐
ACUROS
) was compared to the (1) published TG‐43 data, (2)
MCNPX
Monte Carlo (MC) simulations of the
source centered in a 15 cm radius water sphere, and (3) TG‐43 output from the TPS using vendor supplied (BV‐TG43‐vendor) and user extended (BV‐TG43‐extended) 2D anisotropy functions
. BV‐
ACUROS
assumes 1 mm of NiTi cable, while the TPS TG‐43 algorithm uses data based on a 15 cm cable. MC models of various cable lengths were simulated.
Results:
The MC simulations resulted in
dose deviations along the cable for 1, 2, and 3 mm cable lengths relative to 15 cm. BV‐
ACUROS
comparisons with BV‐TG43‐vendor and BV‐TG43‐extended yielded magnitude of differences, consistent with those seen in MC simulations. However, differences
extended further
when using the vendor supplied anisotropy function
. These differences were also seen in comparisons of
derived from the TPS output.
Conclusions:
The results suggest that
near the cable region is larger than previously estimated. The spatial distribution of the dose deviation is highly dependent on the reference TG‐43 data used to compare to GBBS. The differences observed, while important to realize, should not have an impact on clinical dosimetry in homogeneous water.
To investigate the potential of a novel deterministic solver, Attila, for external photon beam radiotherapy dose calculations.
Two hypothetical cases for prostate and head-and-neck cancer photon beam ...treatment plans were calculated using Attila and EGSnrc Monte Carlo simulations. Open beams were modeled as isotropic photon point sources collimated to specified field sizes. The sources had a realistic energy spectrum calculated by Monte Carlo for a Varian Clinac 2100 operated in a 6-MV photon mode. The Attila computational grids consisted of 106,000 elements, or 424,000 spatial degrees of freedom, for the prostate case, and 123,000 tetrahedral elements, or 492,000 spatial degrees of freedom, for the head-and-neck cases.
For both cases, results demonstrate excellent agreement between Attila and EGSnrc in all areas, including the build-up regions, near heterogeneities, and at the beam penumbra. Dose agreement for 99% of the voxels was within the 3% (relative point-wise difference) or 3-mm distance-to-agreement criterion. Localized differences between the Attila and EGSnrc results were observed at bone and soft-tissue interfaces and are attributable to the effect of voxel material homogenization in calculating dose-to-medium in EGSnrc. For both cases, Attila calculation times were <20 central processing unit minutes on a single 2.2-GHz AMD Opteron processor.
The methods in Attila have the potential to be the basis for an efficient dose engine for patient-specific treatment planning, providing accuracy similar to that obtained by Monte Carlo.
Abstract Purpose For the custom‐built construction of eye plaques, the iodine (I‐125) seeds of different source strengths are recycled in our eye plaque program. To return I‐125 seeds to the correct ...lot, we developed a novel 3D‐printed conical plaque QA holder for relative assay for eye plaques. Materials and methods A universal 3D‐printed conical plaque holder was designed to accommodate six plaque sizes and fit reproducibly in a well‐type dose calibrator. A reproducibility test was used to compare the plaque placement consistency in the holder versus without the holder. Plaque assays were performed for assembled plaques both before implant and after explant. The explant reading was compared with the implant reading adjusted for decay, and the relative error was calculated. The plaque response fraction (PRF) is defined as the fraction of well chamber implant reading over the total seed strength for a plaque. The PRF was aggregated for each individual plaque to confirm the seed lot before implant. Results The reproducibility test showed the chamber reading's relative standard deviation of 0.40% with the QA holder compared to 0.68% without it. The batch relative assay was performed for 251 plaques. The absolute value of measurement deviation between explant and decay‐corrected implant readings is 0.89% ± 0.86% (mean ± standard deviation). The PRFs for individual plaques range from 36.49% to 49.87%, with a maximum standard deviation of 2%. Conclusions This novel 3D‐printed QA holder provides reproducible setup for assaying assembled eye plaques in a well chamber. Batch relative assay can validate the seed batch used and plaque integrity during the implant without assaying individual seeds, saving valuable physicist time and radiation exposure from seed handling.
With the recent introduction of heterogeneity correction algorithms for brachytherapy, the AAPM community is still unclear on how to commission and implement these into clinical practice. The ...recently‐published AAPM TG‐186 report discusses important issues for clinical implementation of these algorithms. A charge of the AAPM‐ESTRO‐ABG Working Group on MBDCA in Brachytherapy (WGMBDCA) is the development of a set of well‐defined test case plans, available as references in the software commissioning process to be performed by clinical end‐users. In this practical medical physics course, specific examples on how to perform the commissioning process are presented, as well as descriptions of the clinical impact from recent literature reporting comparisons of TG‐43 and heterogeneity‐based dosimetry.
Learning Objectives:
1.Identify key clinical applications needing advanced dose calculation in brachytherapy.
2.Review TG‐186 and WGMBDCA guidelines, commission process, and dosimetry benchmarks.
3.Evaluate clinical cases using commercially available systems and compare to TG‐43 dosimetry
Purpose:
The deterministic Acuros XB (AXB) algorithm was recently implemented in the Eclipse treatment planning system. The goal of this study was to compare AXB performance to Monte Carlo (MC) and ...two standard clinical convolution methods: the anisotropic analytical algorithm (AAA) and the collapsed-cone convolution (CCC) method.
Methods:
Homogeneous water and multilayer slab virtual phantoms were used for this study. The multilayer slab phantom had three different materials, representing soft tissue, bone, and lung. Depth dose and lateral dose profiles from AXB v10 in Eclipse were compared to AAA v10 in Eclipse, CCC in Pinnacle3, and EGSnrc MC simulations for 6 and 18 MV photon beams with open fields for both phantoms. In order to further reveal the dosimetric differences between AXB and AAA or CCC, three-dimensional (3D) gamma index analyses were conducted in slab regions and subregions defined by AAPM Task Group 53.
Results:
The AXB calculations were found to be closer to MC than both AAA and CCC for all the investigated plans, especially in bone and lung regions. The average differences of depth dose profiles between MC and AXB, AAA, or CCC was within 1.1, 4.4, and 2.2%, respectively, for all fields and energies. More specifically, those differences in bone region were up to 1.1, 6.4, and 1.6%; in lung region were up to 0.9, 11.6, and 4.5% for AXB, AAA, and CCC, respectively. AXB was also found to have better dose predictions than AAA and CCC at the tissue interfaces where backscatter occurs. 3D gamma index analyses (percent of dose voxels passing a 2%/2 mm criterion) showed that the dose differences between AAA and AXB are significant (under 60% passed) in the bone region for all field sizes of 6 MV and in the lung region for most of field sizes of both energies. The difference between AXB and CCC was generally small (over 90% passed) except in the lung region for 18 MV 10 × 10 cm2 fields (over 26% passed) and in the bone region for 5 × 5 and 10 × 10 cm2 fields (over 64% passed). With the criterion relaxed to 5%/2 mm, the pass rates were over 90% for both AAA and CCC relative to AXB for all energies and fields, with the exception of AAA 18 MV 2.5 × 2.5 cm2 field, which still did not pass.
Conclusions:
In heterogeneous media, AXB dose prediction ability appears to be comparable to MC and superior to current clinical convolution methods. The dose differences between AXB and AAA or CCC are mainly in the bone, lung, and interface regions. The spatial distributions of these differences depend on the field sizes and energies.
Purpose:
To evaluate the dose distributions of an
I
192
r
source (model VS2000) in homogeneous water geometry calculated using a deterministic grid-based Boltzmann transport equation solver (GBBS) in ...the commercial treatment planning system (TPS) (BRACHYVISION-ACUROS v8.8).
Methods:
Using percent dose differences
(
%
Δ
D
)
, the GBBS (BV-ACUROS) was compared to the (1) published TG-43 data, (2) MCNPX Monte Carlo (MC) simulations of the
I
192
r
source centered in a 15 cm radius water sphere, and (3) TG-43 output from the TPS using vendor supplied (BV-TG43-vendor) and user extended (BV-TG43-extended) 2D anisotropy functions
F
(
r
,
θ
)
. BV-ACUROS assumes 1 mm of NiTi cable, while the TPS TG-43 algorithm uses data based on a 15 cm cable. MC models of various cable lengths were simulated.
Results:
The MC simulations resulted in
>
20
%
dose deviations along the cable for 1, 2, and 3 mm cable lengths relative to 15 cm. BV-ACUROS comparisons with BV-TG43-vendor and BV-TG43-extended yielded magnitude of differences, consistent with those seen in MC simulations. However, differences
>
20
%
extended further
(
θ
≤
10
°
)
when using the vendor supplied anisotropy function
F
ven
(
r
,
θ
)
. These differences were also seen in comparisons of
F
(
r
,
θ
)
derived from the TPS output.
Conclusions:
The results suggest that
%
Δ
D
near the cable region is larger than previously estimated. The spatial distribution of the dose deviation is highly dependent on the reference TG-43 data used to compare to GBBS. The differences observed, while important to realize, should not have an impact on clinical dosimetry in homogeneous water.