- Open Access
Clinical evaluation for the difference of absorbed doses calculated to medium and calculated to water by Monte Carlo method
Radiation Oncology volume 13, Article number: 137 (2018)
To evaluate the difference of absorbed doses calculated to medium and to water by a Monte Carlo (MC) algorithm based treatment planning system (TPS), and to assess the potential clinical impact to dose prescription.
Thirty patients, 10 nasopharyngeal cancer (NPC), 10 lung cancer and 10 bone metastases cases, were selected for this study. For each case, the treatment plan was generated using a commercial MC based TPS and dose was calculated to medium (Dm). The plan was recalculated for dose to water (Dw) using the same Monitor Units (MU) and control points. The differences between Dm and Dw were qualitatively evaluated by dose-volume parameters and by the plan subtraction method. All plans were measured using the MapCheck2, and gamma passing rates were calculated.
For NPC and Lung cases, the mean differences between Dw and Dm for the targets were less than 2% and the maximum difference was 3.9%. The maximum difference of D2% for the organs at risk (OARs) was 6.7%. The maximum differences between Dw and Dm were as high as 10% in certain high density regions. For bone metastases cases, the mean differences between Dw and Dm for the targets were more than 2.2% and the maximum difference was 7.1%. The differences between Dw and Dm for the OARs were basically negligible. At 3%&3 mm criterion, the gamma passing rate of Dw plan and Dm plan were close (> 94%).
The differences between Dw and Dm has little clinical impact for most clinical cases. In bony structures the differences may become clinically significant if the target/OAR is receiving doses close to its tolerance limit which can potentially influence the selection or rejection of a particular plan.
Absorbed dose is an important parameter in characterizing the effect of radiation therapy for the efficacy of tumor eradication and protection from unacceptable damage to normal organs . For historical reasons, in terms of dose, Dw has been assumed for reporting the dose to various media. However, human body is not only composed of water. Many tissues in the body have different densities than water, especially the bones and lung. For radiation therapy the dose absorbed to water cannot accurately represent the actual dose absorbed in different tissues. In practice, traditional treatment planning system (TPS) typically takes the effect of different tissue densities with attenuation and scatter into considerations but reports the dose at each location as the dose to water. Monte Carlo (MC) algorithm is the most accurate algorithm for dose calculation in that it simulates the transport properties of various particles in various media in the region of interest and scores the dose contribution locally to the medium with its assigned chemical composition as well as density. The resulting dose distributions may be different from those calculated by traditional dose calculation algorithms, especially for tissues of heterogeneity [2,3,4]. In recent years, MC has been increasingly adopted in clinical application [5,6,7]. There are a number of reasons for using Dw for reporting of MC calculated doses. Two major ones are that it has been used in decades of clinical studies for outcome correlation with the dose, and that the calibration protocols are all referenced to water. A technical issue related to dose calibration is that an MC based TPS could model the chemical composition of various biological tissues by approximation as a function of Computed Tomography (CT) numbers based on data of the human body (reference International Commission on Radiation Units & Measurements reports 44 and 46). Such an approximation may not perform well for non-biological materials like in a quality assurance (QA) phantom. MC based dose calculations typically report absorbed dose to media (Dm). Therefore there is a need to convert between Dm and Dw, and, as Siebers JV et al.  argued, MC is capable of doing the conversion. Siebers et al. presented a method to calculate the difference between Dm and Dw by applying the Bragg-Gray cavity theory, and their results showed a difference exceeding 10% in cortical bones.
Currently there is no consensus regarding whether Dm or Dw should be used for an MC based TPS [9, 10]. When it comes to clinical application, the difference between Dw and Dm will affect interpretation of dose distribution and perhaps the value of prescription dose, leading to differences in plan evaluation, dose reporting, and dose verification. In this work, Dm and Dw were both calculated using Monaco TPS for 10 nasopharyngeal cancer (NPC) cases, 10 lung cancer cases and 10 bone target cases, in order to investigate the issue in two common clinical sites in which differences of dose distributions may be highlighted. Dose Volume Histogram (DVH) was used to analyze dose parameters in the target and organ at risk (OAR), and three dimensional dose difference distributions between Dm and Dw were calculated. Gamma passing rates (measurement results vs Dm/Dw plans) were calculated at different QA criteria to evaluate the dose accuracy.
Dm plan originally created for treatment
Ten NPC cases in stage T3 or T4, 10 lung cancer cases and 10 bone target cases (7 cases of lumbar vertebra metastasis, 3 cases of thoracic vertebra metastasis) treated at Sun Yat-sen University Cancer Center were retrospectively chosen in this study. The gross tumor volumes (GTVs) and clinical tumor volume (CTV) were contoured by experienced radiation oncologists according to definitions in the ICRU 50 and ICRU 62 reports [11, 12], and the planning target volume (PTV) were generated following a set of physician prescribed margins that were consistent with departmental protocols specific to the disease sites. Monaco TPS (Version 5.0, Elekta) was used to create the treatment plans for step-and-shoot IMRT with an Elekta Synergy linac, and MC calculated Dm was chosen for dose reporting. Nine equally spaced fields were used for NPC cases. The prescription of NPC cases and Lung cancer cases were 70 Gy (32 or 33 fractions, 5 days/week) and 65 Gy (26 fractions, 5 days/week) respectively. The main planning objectives for NPC are PTV V100% > 98% and PTV V110% < 10% (Vx%, is the percentage volume of reign of interest (ROI) that receives at least x% prescription dose), spinal cord D2% < 45Gy, brain stem D2% < 54Gy, parotid gland D50% < 30Gy, optical nerve D2% < 54Gy, and the dose to lens as low as possible. For lung IMRT cases 5–7 fields were used. The planning objectives are PTV V100% > 95% and PTV V110% < 2%, spinal cord D2% < 45Gy, normal lung V20 Gy < 35% (VD Gy, is the percentage volume of ROI that receives at least absorbed dose D) and normal lung mean dose <19Gy, heart V30 Gy < 40%, and the maximum esophagus dose <65Gy. For bone target cases, 5–7 fields were used. The prescription of bone target cases was 25 Gy (5Gy/fractions, 5 days/week). The main planning objectives are for PTV, V100% > 95% and V110% < 10%, for spinal cord Dmax < 26 Gy, for lung V10Gy < 15%, and the maximum esophagus dose < 26 Gy.
The MC algorithm in the Monaco TPS used for this study, called XVMC, calculates dose based upon mass density. A technical issue of dose calculation with MC in treatment planning is how to obtain the density and chemical composition data for the patient model from the CT. An approximation is made by assigning a voxel to certain type of tissue in the human body based on its Hounsfield unit (HU) in a certain range, and the mass density and composition data can be looked up in the International Commission on Radiation Units & Measurements Reports No. 46 . XVMC algorithm converts CT numbers to ED numbers using the user-defined CT-to-ED calibration table and takes with a fit function that maps continuously the electron density to mass density for matching a tissue with approximating cross section and attenuation coefficient data .
The conversion to Dw can be calculated based on the distribution of Dm plan according to the Bragg-Gray cavity theory:
where sw,med is the mean unconstrained mass stop power ratio of water to media of primary electron spectrum, and Dw is understood as the dose to the voxel replacement of water embedded to the actual media. Theoretically mass stop power ratio can be calculated by the following formula :
where (S/ρ)w and (S/ρ)med are the unconstrained mass stop power of water and media, respectively. (ΦE)m is the primary electron fluence in the medium and Emax is the maximum energy in the (ΦE)m distribution. The stopping power ratio in Moncao was pre-calculated by approximation for tissue-like media.
The conversion from Dm to Dw in Monaco with a clinically accepted plan involved a simple recalculation with exactly the same set of plan parameters (all the geometric parameters and monitor units (MU)) retained. The stopping power ratios dependent of mass density were applied voxel by voxel. The matrix of dose calculation grid was 0.3 cm × 0.3 cm × 0.3 cm, and the Monte Carlo statistical uncertainty was set at 3% per control point.
Dm and Dw dose verification
All the plans were measured with MapCHECK2 (Sun Nuclear, Florida, USA) to verify the dose distribution. MacpCHECK2 was mounted in a water-equivalent phantom (MapPHAN) with a 5 cm equivalent depth from the surface to the detectors. The TPS planed dose was calculated on the real phantom CT images without overriding the density. The measured dose distributions of composite fields were compared with the corresponding planned dose distributions (Dm or Dw), and the local dose normalization gamma (γ) passing rates were calculated at the setting dose difference (DD) and distance to agreement (DTA). In order to eliminate dose in the out-of-field region where a large relative dose difference can be calculated and hence skew theγ result, a lower dose threshold (10%) was set and below the threshold theγ result was ignored. Using 3%&3 mm, 2%&2 mm and 1%&1 mm tolerances, the gamma passing rates were calculated to find how the pass rates change with reduction of dose difference and DTA limits.
According to the ICRU 83 report, the volume-dose is recommended to describe the dose information in the ROIs, as Dx% to note the dose that X% of volume of ROI receives . For example, D98% means 98% of volume received the dose at specified value such as 65Gy. These DVH parameters were used for statistical analysis of Dw and Dm dose distributions. The bin width of the DVHs was 1 cGy, and the resolution for DVH sampling was 0.1 cm. The difference between the Dw and Dm was calculated by:
The plan subtraction method was used to evaluate the spatial dose difference distribution of Dw and Dm.
Paired t-tests were performed using the SPSS software (Version 19, SPSS, Inc., USA) to determine the statistical significance of the difference between Dw and Dm, with a p-value < 0.05 as the threshold for consideration as statistically significant.
Dw and Dm for NPC cases
Figure 1 shows the comparison of the DVH results with Dw and Dm for a typical NPC treatment plan. There were small but systematic deviations from Dm to Dw in the planning target volumes (PTVs). Table 1 shows the mean and difference in dose-volume indices calculated with MC, evaluated for 10 NPC cases. Except for the D50% and D2% of PTV66, and D98% of PTV54, all DVH indices for all PTVs were different with statistical significance (p < 0.05), including D98%, D50%, and D2% (Dx%, the minimum dose that x% of the volume of the organ receives from the cumulative DVH). The possible reason for PTV66 behaved differently from the others may be that PTV66 is the lymph gland target, small in size and relatively variable in location among different patients. For the D2% of PTV70, PTV66, PTV60 and PTV54, the values of the Dm plan are less than that of Dw, and the mean deviation was 1.9 ± 1.1%, 0.4 ± 1.0%, 1.7 ± 1.0% and 1.3 ± 0.7%, respectively. The difference between Dw and Dm in the mean dose of PTVs were within 1%.
As for the OARs, the D50% increased when Dm was converted to Dw, and this was a statistically significant result except for the optic nerve and parotid gland. The median dose of T-M joints and mandibular in the Dm plans were at least 5% less than that in the Dw plans. The D2% of spinal cord, brain stem, parotid gland, lens, optic nerves, temporal lobe, and tongue increased by less than 1% from Dm to Dw. However, the D2% of T-M joints and mandibular suffered about 5% change from Dm to Dw.
Dw and Dm for lung cancer cases
Figure 2 shows that, for lung cancer cases, the difference between Dw and Dm is less obvious than in the NPC cases. Table 2 shows that the D2% of PTV65 and the D98% of PTV50 were statistically significant (p < 0.05), and the mean deviation were 0.3 ± 0.4% and − 0.3 ± 0.3%, respectively. There were no other statistically significant differences for other DVH indices evaluated for PTVs. All deviations were with 1%. For the OARs, the median dose D50% of spinal cord and heart were slightly increased from Dm to Dw with the mean deviation at 0.3 ± 0.3% and 1.1 ± 0.5%, respectively, and this was statistically significant. There were no statistically significant differences between Dw and Dm in lung and esophagus. For the D2% of spinal cord, lung, esophagus and heart, there were statistically significant differences between Dw and Dm, and the mean deviation were 0.3 ± 0.4%, − 0.6 ± 0.5%, − 0.7 ± 0.5%, and 0.6 ± 0.6%, respectively. All the differences in the DVH indices evaluated were within 2%.
Dw and Dm for bone target cases
Figure 3 shows that, for bone metastases cases, the differences between Dw and Dm for PTV targets are more obvious than those in the NPC cases and lung cases. From Table 3, all DVH indices for the PTVs were different with statistical significance (p < 0.01). The D98%, D50%, and D2% deviation of PTV25 were 3.0 ± 1.2%, 3.5 ± 1.4% and 4.4 ± 1.9%, respectively. For the PTV20, D98%, D50%, and D2% deviations were 2.2 ± 0.7%, 2.8 ± 0.7% and 3.8 ± 1.7%, respectively. There were basically negligible differences between Dw and Dm in spinal, lung and esophagus. All the differences in the DVH indices evaluated for OARs were within 0.6%.
Dose difference distribution maps
By subtracting the re-calculated Dw plan and original Dm plans, the dose difference of three-dimensional distribution can be obtained. The dose difference (diff) is defined by diff (%) = (Dw - Dm)/ Dp × 100, where Dp is the prescription dose. Figure 4 shows the difference distribution in three-dimensions of a typical NPC case between Dw and Dm. A typical case of lung cancer is shown in Fig. 5 and a case of bone metastasis is shown in Fig. 6. The blue to purple gradient legend represented the dose difference values ranging from 0 to 10%. It can be seen from Fig. 4 and Fig. 5 that the difference between Dw and Dm could be higher than 5% in bone, while the differences between Dw and Dm in soft tissues were less obvious (usually smaller than 3%). From Fig. 6 the differences between Dw and Dm in thoracic vertebra bone were about 3–8%, a little lower than the result in head bone in Fig. 3. It’s probably because the bone density of the thoracic vertebra is different from that of the head bone.
At normal QA criterion, 3% dose difference and 3 mm distance to agreement, the gamma pass rates of Dw plan and Dm plan are all above 94% and very close. But when the tolerances become stricter, the gamma passing rates decreases dramatically, and Dw plans gamma pass rates become better than the Dm plans (Table 4).
With the application of MC algorithm for dose calculation in radiation therapy, whether the dose should be calculated to medium or to water has been an unsettled debate [9, 10, 16]. The arguments that support Dw include that beam data was measured in water, that the beam output was calibrated in water, and that most clinical experience were based on dose to water, etc. However, the compelling argument to support the use of Dm is that it represents the true dose at each location of specific medium. It is the unique advantage of Monte Carlo in that Dm can be calculated directly, but Dm to Dw using stopping power ratios may involve an uncertainty . In reality, different TPS use different dose calculation algorithms to produce Dw, from direct calculation to applying conversion factors. According to the AAPM TG 105 report , when the element components are considered in dose calculation, both Dm and Dw should be available for evaluation. When comes to a specific clinical situation, the difference between Dm and Dw should be known. N Dogan et al.  showed that converting Dm to Dw in EGS4 MC-calculated IMRT treatment plans introduces a systematic error in target and critical structure DVHs, and this systematic error may reach up to 5.8% for H&N and 8.0% for prostate cases when the hard-bone-containing structures such as femoral heads are present.
From our work using Monaco for NPC and lung cancer, Dm was less than Dw. The mean deviation for soft tissues was within 2%. For T-M joints and mandibular, the mean deviation was greater than 5%, and in regions of unspecified normal bone the difference could reach 10%. Our results agreed nicely with the work by Siebers et al. . It is interesting to find, based on our study, that there was hardly any difference between Dw and Dm in low density regions. Although the stopping power ratio for both cortical bone and air can be above 1.10, the stopping power ratio is close to 1 for low density tissues like lung. For this reason, the issue with using Dw or Dm may have a minimal effect for majority of clinical situations.
The dose difference between Dw and Dm in bony structures may become clinically significant if the OAR is receiving doses close to its tolerance dose limit which can influence selection or rejection of a particular plan. The dose calculated by MC may need to be carefully evaluated in certain situations, e.g. bone metastasis, bone tumor, or constraining a hot spot in bone that becomes a limiting factor in plan optimization. From the Fig. 3, for PTV of the bone target cases, though the target dose coverages (the target volume (%) received the prescription dose) of Dm and converted Dw plan were similar, the mean median dose of Dw plan increased by 3.5% comparing with that of Dm plan (Table 3). That means the dose prescription for bone target could be about 3.5% higher than that of using Dw dose, and their treatment response and outcome may need further study in the future.
Previous studies [16, 20] using EGS4/MCSIM Monte Carlo and AXB dose calculations proved that conventional model based algorithms predicted dose distributions in bone that were closer to Dm distributions than to Dw distributions. It is therefore better to use Dm for consistency with previous radiation therapy experience. Our measurements showed that at widely used reference standard, 3% dose difference and 3 mm DTA, the Dm and Dw plan gamma passing rates were very close, but when the gamma calculation standard became stricter, the Dw was closer to the result of measurement than the Dm. That’s because the MapCheck2 CT images without forcing density were used to calculate the planned dose distribution, where the MapCheck2 detectors are made of high density metallic elements and the detectors are always calibrated by Dw. The CT scanner used for acquisition of patient simulation images has the limitation of scanning high density material such as the diode and the TPS also has limitation while accepting CT images with high density material. In our practice, Dm is used for treatment planning, and physicians and physicists will be consulted in case conversion to Dw in bone may affect the decisions to choose the appropriate dose distribution for treatment.
Conversion to Dw may be necessary for dose verification in the quality assurance phantom. If a water phantom is used, the difference between Dm and Dw can be ignored. Kan MW et al.  showed that for a heterogeneous phantom with high density materials contained the difference between Dm and Dw has an effect on the passing rate of QA measurement. Our results (Table 4) showed there were obvious differences between the Dm and Dw plan gamma passing rates when the QA criteria became strict. A simple method to bypass the problem is to assign a uniform density to the phantom and calculate to either Dm or Dw in a consistent manner. The choice of an appropriate density needs to be validated by an independent method such as point dose measurement.
Overall, the dose differences between Dm and Dw calculated by MC algorithm in Monaco are small in regions that have densities close or low to water. Our results show that dose calculated to medium by Monaco can be used clinically. In high density regions like cortical bone, the difference was 5 to 10%, and this may have a clinical consequence and needs to be carefully considered in certain clinical situations.
Clinical target volume
- Dm :
Dose to media
Distance to agreement
Dose volume histogram
- Dw :
Dose to water
Gross tumor volume
Intensity modulated radiation therapy
Organ at risk
Planning target volume
Region of interest
Treatment planning system
International Commission on Radiation Units and Measurements (ICRU). ICRU Report No.33: Radiation quantities and units. Washington: ICRU; 1980.
Papanikolaou N, Battista J, Boyer A, Kappas C, Kein E, Mackie TR, Sharpe M, Dyk JV. AAPM Report No.85: Tissue inhomogeneity corrections for megavoltage photon beams. AAPM Report No. 85. Madison: Medical Physics Publishing; 2004. p. 1–135.
Fogliata A, Nicolini G, Vanetti E, Clivio A, Winkler P, Cozzi L. The impact of photon dose calculation algorithms on expected dose distributions in lungs under different respiratory phases. Phys Med Biol. 2008;53:2375–90.
Kathirvel M, Subramanian S, Clivio A, Arun G, Fogliata A, Nicolini G, Subramani V, Thirumalai Swamy S, Vanetti E, Cozzi L. Critical appraisal of the accuracy of Acuros-XB and anisotropic analytical algorithm compared to measurement and calculations with the compass system in the delivery of RapidArc clinical plans. Radiat Oncol. 2013;8(1):140.
Hartmann Siantar CL, Walling RS, Daly TP, Faddegon B, Albright N, Bergstrom P, Beilajew AF, Chuang C, Garrett D, RK H, Knapp D, Wieczorek DJ, Verhey LJ. Description and dosimetric verification of the PEREGRINE Monte Carlo dose calculation system for photon beams incident on a water phantom. Med Phys. 2001;28:1322–37.
Mika S, Christ G. Experimental validation of a Monte Carlo-based treatment-planning system for Electron beams. Strahlenther Onkol. 2007;3:150–6.
Li J, Doemer A, Harrison AS, Galvin JM, Xiao Y. Dose calculation using Monte Carlo algorithm of Monaco treatment planning system. Int J of Radiat Oncol Bio Phys. 2008;72:S674.
Siebers JV, Keall PJ, Nahum AE, Mohan M. Converting absorbed dose to medium to absorbed dose to water for Monte Carlo based photon beam dose calculations. Phys Med Biol. 2000;45:983–95.
Liu H, Keall P. Dm rather than Dw should be used in Monte Carlo treatment planning. Med Phys. 2002;29:922–4.
Andreo P. Dose to ‘water-like’ media or dose to tissue in MV photons radiotherapy treatment planning: still a matter of debate. Phys Med Biol. 2015;60:309–37.
International Commission on Radiation Units and Measurements (ICRU). Report 50: Prescribing, recording, and reporting photon beam therapy. Bethesda: ICRU; 1993.
International Commission on Radiation Units and Measurements (ICRU). Report 62: Prescribing, recording and reporting photon beam therapy (supplement to ICRU report 50). Bethesda: ICRU; 1999.
International Commission on Radiation Units and Measurements (ICRU). Report 46: Photon, electron, proton and neutron interaction data for body issues. Bethesda: ICRU; 1992.
Fippel M. Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm. Med Phys. 1999;26:1466-75.
International Commission on Radiation Units and Measurements (ICRU). Report 83: Prescribing, recording and reporting photon-beam Intensity-Modulated Radiation Therapy (IMRT). Journal of ICRU, 2010; https://doi.org/10.1093/jicru/ndq002.
Ma CM, Li J. Dose specification for radiation therapy: dose to water or dose to medium? Phys Med Biol. 2011;56:3073–89.
Fippel M, Nusslin F. Comments on ‘ converting absorbed dose to medium to absorbed dose to water for Monte Carlo based photon beam dose calculations’. Phys Med Biol. 2000;45:L17.
Chetty IJ, Curran B, Cygler JE, DeMarco JJ, Ezzel G, Faddegon BA, Kawrakow I, Keall PJ, Liu H, Ma CM, Rogers DW, Seuntjens J, Bagheri DS, Siebers JV. Report of the AAPM Task Group No.105: Issues associated with clinical implementation of Monte Carl-based photon and electron external beam treatment planning. Med Phys. 2007;34:4818–53.
Dogam N, Siebers JV, Keall PJ. Clinical comparison of head and neck and prostate IMRT plans using absorbed dose to medium and absorbed dose to water. Phys Med Biol. 2006;51:4967–80.
Kan MW, Leung LH, So RW, Yu PK. Experimental verification of the Acuros XB and AAA dose calculation adjacent to heterogeneous media for IMRT and RapidArc of nasopharygeal carcinoma. Med Phys. 2013;40:031714. 1–19
The authors would like to thank Dr. Yan Chen for many very helpful discussions.
Availability of data and materials
The datasets are backed up on the Research Data Deposit (RDD, http://www.researchdata.org.cn, approval number: RDDA2017000317) and are available on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, L., Huang, B., Huang, X. et al. Clinical evaluation for the difference of absorbed doses calculated to medium and calculated to water by Monte Carlo method. Radiat Oncol 13, 137 (2018) doi:10.1186/s13014-018-1081-3
- Monte Carlo dose calculation
- Absorbed dose to medium
- Absorbed dose to water
- Nasopharyngeal cancer
- Gamma analysis