Showing papers on "Dosimetry published in 2013"

In vivo dosimetry in external beam radiotherapy

Abstract: In vivo dosimetry (IVD) is in use in external beam radiotherapy (EBRT) to detect major errors, to assess clinically relevant differences between planned and delivered dose, to record dose received by individual patients, and to fulfill legal requirements. After discussing briefly the main characteristics of the most commonly applied IVD systems, the clinical experience of IVD during EBRT will be summarized. Advancement of the traditional aspects of in vivo dosimetry as well as the development of currently available and newly emerging noninterventional technologies are required for large-scale implementation of IVD in EBRT. These new technologies include the development of electronic portal imaging devices for 2D and 3D patient dosimetry during advanced treatment techniques, such as IMRT and VMAT, and the use of IVD in proton and ion radiotherapy by measuring the decay of radiation-induced radionuclides. In the final analysis, we will show in this Vision 20/20 paper that in addition to regulatory compliance and reimbursement issues, the rationale for in vivo measurements is to provide an accurate and independent verification of the overall treatment procedure. It will enable the identification of potential errors in dose calculation, data transfer, dose delivery, patient setup, and changes in patient anatomy. It is the authors’ opinion that all treatments with curative intent should be verified through in vivo dose measurements in combination with pretreatment checks.

Impact of plan parameters on the dosimetric accuracy of volumetric modulated arc therapy.

Abstract: Purpose: To evaluate the effect of plan parameters on volumetric modulated arc therapy (VMAT) dosimetric accuracy, together with the possibility of scoring plan complexity. Methods: 142 clinical VMAT plans initially optimized using a 4° control point (CP) separation were evaluated. All plans were delivered by a 6 MV Linac to a biplanar diode array for patient-specific quality assurance (QA). Local Γ index analysis (3%, 3 mm and 2%, 2 mm) enabled the comparison between delivered and calculated dose. The following parameters were considered for each plan: average leaf travel (LT), modulation complexity score applied to VMAT (MCSv), MU value, and a multiplicative combination of LT and MCSv (LTMCS). Pearson's correlation analysis was performed between Γ passing rates and each parameter. The effects of CP angular separation on VMAT dosimetric accuracy were also analyzed by focusing on plans with high LT values. Forty out of 142 plans with LT above 350 mm were further optimized using a finer angle spacing (3° or 2°) and Γ analysis was performed. The average Γ passing rates obtained at 4° and at 3°/2° sampling were compared. A further correlation analysis between all parameters and the Γ pass-rates was performed on 142 plans, but including the newly optimized 40 plans (CP every 3° or 2°) in place of the old ones (CP every 4°). Results: A moderate significant (p < 0.05) correlation between each examined parameter and Γ passing rates was observed for the original 142 plans at 4° CP discretization. A negative correlation was found for LT with Pearson's r absolute values above 0.6, suggesting that a lower dosimetric accuracy may be expected for higher LT values when a 4° CP sampling is used. A positive correlation was observed for MCSv and LTMCS with r values above 0.5. In order to score plan complexity, threshold values of LTMCS were defined. The average Γ passing rates were significantly higher for the plans created using the finer CP spacing (3°/2°) compared to the plans optimized using the standard 4° spacing (Student t-test p < 0.05). The correlation between LT and passing rates was strongly diminished when plans with finer angular separations were considered, yielding Pearson's r absolute values below 0.45. Conclusions: At 4° CP sampling, LT, MCSv, and LTMCS were found to be significantly correlated with VMAT dosimetric accuracy, expressed as Γ pass-rates. These parameters were found to be possible candidates for scoring plan complexity using threshold values. A finer CP separation (3°/2°) led to a significant increase in dosimetric accuracy for plans with high leaf travel values, and to a decrease in correlation between LT and Γ passing rates. These results indicated that the influence of LT on VMAT dosimetric accuracy can be controlled by reducing CP separation. CP spacing for all plans requiring large leaf motion should not exceed 3°. The reported data were integrated to optimize our clinical workflow for plan creation, optimization, selection among rival plans, and patient-specific QA of VMAT treatments.

A dual model HU conversion from MRI intensity values within and outside of bone segment for MRI-based radiotherapy treatment planning of prostate cancer.

Abstract: Purpose: The lack of electron density information in magnetic resonance images (MRI) poses a major challenge for MRI-based radiotherapy treatment planning (RTP). In this study the authors convert MRI intensity values into Hounsfield units (HUs) in the male pelvis and thus enable accurate MRI-based RTP for prostate cancer patients with varying tissue anatomy and body fat contents. Methods: T{sub 1}/T{sub 2}*-weighted MRI intensity values and standard computed tomography (CT) image HUs in the male pelvis were analyzed using image data of 10 prostate cancer patients. The collected data were utilized to generate a dual model HU conversion technique from MRI intensity values of the single image set separately within and outside of contoured pelvic bones. Within the bone segment local MRI intensity values were converted to HUs by applying a second-order polynomial model. This model was tuned for each patient by two patient-specific adjustments: MR signal normalization to correct shifts in absolute intensity level and application of a cutoff value to accurately represent low density bony tissue HUs. For soft tissues, such as fat and muscle, located outside of the bone contours, a threshold-based segmentation method without requirements for any patient-specific adjustments was introduced to convert MRI intensity values intomore » HUs. The dual model HU conversion technique was implemented by constructing pseudo-CT images for 10 other prostate cancer patients. The feasibility of these images for RTP was evaluated by comparing HUs in the generated pseudo-CT images with those in standard CT images, and by determining deviations in MRI-based dose distributions compared to those in CT images with 7-field intensity modulated radiation therapy (IMRT) with the anisotropic analytical algorithm and 360° volumetric-modulated arc therapy (VMAT) with the Voxel Monte Carlo algorithm. Results: The average HU differences between the constructed pseudo-CT images and standard CT images of each test patient ranged from −2 to 5 HUs and from 22 to 78 HUs in soft and bony tissues, respectively. The average local absolute value differences were 11 HUs in soft tissues and 99 HUs in bones. The planning target volume doses (volumes 95%, 50%, 5%) in the pseudo-CT images were within 0.8% compared to those in CT images in all of the 20 treatment plans. The average deviation was 0.3%. With all the test patients over 94% (IMRT) and 92% (VMAT) of dose points within body (lower than 10% of maximum dose suppressed) passed the 1 mm and 1% 2D gamma index criterion. The statistical tests (t- and F-tests) showed significantly improved (p ≤ 0.05) HU and dose calculation accuracies with the soft tissue conversion method instead of homogeneous representation of these tissues in MRI-based RTP images. Conclusions: This study indicates that it is possible to construct high quality pseudo-CT images by converting the intensity values of a single MRI series into HUs in the male pelvis, and to use these images for accurate MRI-based prostate RTP dose calculations.« less

Effectiveness of robust optimization in intensity-modulated proton therapy planning for head and neck cancers

Abstract: Purpose: Intensity-modulated proton therapy (IMPT) is highly sensitive to uncertainties in beam range and patient setup. Conventionally, these uncertainties are dealt using geometrically expanded planning target volume (PTV). In this paper, the authors evaluated a robust optimization method that deals with the uncertainties directly during the spot weight optimization to ensure clinical target volume (CTV) coverage without using PTV. The authors compared the two methods for a population of head and neck (H&N) cancer patients. Methods: Two sets of IMPT plans were generated for 14 H&N cases, one being PTV-based conventionally optimized and the other CTV-based robustly optimized. For the PTV-based conventionally optimized plans, the uncertainties are accounted for by expanding CTV to PTV via margins and delivering the prescribed dose to PTV. For the CTV-based robustly optimized plans, spot weight optimization was guided to reduce the discrepancy in doses under extreme setup and range uncertainties directly, while delivering the prescribed dose to CTV rather than PTV. For each of these plans, the authors calculated dose distributions under various uncertainty settings. The root-mean-square dose (RMSD) for each voxel was computed and the area under the RMSD-volume histogram curves (AUC) was used to relatively compare plan robustness. Data derived from the dose volume histogram in the worst-case and nominal doses were used to evaluate the plan optimality. Then the plan evaluation metrics were averaged over the 14 cases and were compared with two-sided paired t tests. Results: CTV-based robust optimization led to more robust (i.e., smaller AUCs) plans for both targets and organs. Under the worst-case scenario and the nominal scenario, CTV-based robustly optimized plans showed better target coverage (i.e., greater D95%), improved dose homogeneity (i.e., smaller D5% − D95%), and lower or equivalent dose to organs at risk. Conclusions: CTV-based robust optimization provided significantly more robust dose distributions to targets and organs than PTV-based conventional optimization in H&N using IMPT. Eliminating the use of PTV and planning directly based on CTV provided better or equivalent normal tissue sparing.

Post-radioembolization yttrium-90 PET/CT - part 2: dose-response and tumor predictive dosimetry for resin microspheres.

Abstract: Background Coincidence imaging of low-abundance yttrium-90 (90Y) internal pair production by positron emission tomography with integrated computed tomography (PET/CT) achieves high-resolution imaging of post-radioembolization microsphere biodistribution. Part 2 analyzes tumor and non-target tissue dose-response by 90Y PET quantification and evaluates the accuracy of tumor 99mTc macroaggregated albumin (MAA) single-photon emission computed tomography with integrated CT (SPECT/CT) predictive dosimetry.

A comparative dosimetric study for treating left-sided breast cancer for small breast size using five different radiotherapy techniques: conventional tangential field, filed-in-filed, Tangential-IMRT, Multi-beam IMRT and VMAT

Abstract: To compare the dosimetry for the left-sided breast cancer treatment using five different radiotherapy techniques. Twenty patients with left sided breast cancer were treated with conservative surgery followed by radiotherapy. They were planned using five different radiotherapy techniques, including: 1) conventional tangential wedge-based fields (TW); 2) field-in-field (FIF) technique; 3) tangential inverse planning intensity-modulated radiation therapy (T-IMRT); 4) multi-field IMRT (M-IMRT); and 5) volumetric modulated arc therapy (VMAT). The CTV, PTV and OARs including the heart, the regions of coronary artery (CA), the contralateral breast, the left and right lung were delineated. The PTV dose was prescribed 50Gy and V47.5≥95%. Same dose constraint was used for all five plans. The planned volumetric dose of PTV and PRV-OARs were compared and analyzed. Except VMAT (Average V47.5 was 94.72%±1.2%), all the other four plans were able to meet the V95% (V47.5) requirement. T-IMRT plan improved the PTV dose homogeneity index (HI) by 0.02 and 0.03 when compared to TW plan and VMAT plan, and decreased the V5, V10 and V20 of all PRV-OARs. However, the high dose volume (≥ 30Gy) of the PRV-OARs in T-IMRT plan had no statistically significant difference compared with the other two inverse plans. In all five plans, the dose volume of coronary artery area showed a strong correlation to the dose volume of the heart (the correlation coefficients were 0.993, 0.996, 1.000, 0.995 and 0.986 respectively). Compared to other techniques, the T-IMRT technology reduced radiation dose exposure to normal tissues and maintained reasonable target homogeneity, VMAT is not recommended for left-sided breast cancer treatment. In five techniques, the dose-volume histogram (DVH) of the heart can be used to predict the dose-volume histogram (DVH) of the coronary artery.

Mapping motion from 4D-MRI to 3D-CT for use in 4D dose calculations: a technical feasibility study

Abstract: Purpose: Target sites affected by organ motion require a time resolved (4D) dose calculation. Typical 4D dose calculations use 4D-CT as a basis. Unfortunately, 4D-CT images have the disadvantage of being a 'snap-shot' of the motion during acquisition and of assuming regularity of breathing. In addition, 4D-CT acquisitions involve a substantial additional dose burden to the patient making many, repeated 4D-CT acquisitions undesirable. Here the authors test the feasibility of an alternative approach to generate patient specific 4D-CT data sets. Methods: In this approach motion information is extracted from 4D-MRI. Simulated 4D-CT data sets [which the authors call 4D-CT(MRI)] are created by warping extracted deformation fields to a static 3D-CT data set. The employment of 4D-MRI sequences for this has the advantage that no assumptions on breathing regularity are made, irregularities in breathing can be studied and, if necessary, many repeat imaging studies (and consequently simulated 4D-CT data sets) can be performed on patients and/or volunteers. The accuracy of 4D-CT(MRI)s has been validated by 4D proton dose calculations. Our 4D dose algorithm takes into account displacements as well as deformations on the originating 4D-CT/4D-CT(MRI) by calculating the dose of each pencil beam based on an individual time stamp of whenmore » that pencil beam is applied. According to corresponding displacement and density-variation-maps the position and the water equivalent range of the dose grid points is adjusted at each time instance. Results: 4D dose distributions, using 4D-CT(MRI) data sets as input were compared to results based on a reference conventional 4D-CT data set capturing similar motion characteristics. Almost identical 4D dose distributions could be achieved, even though scanned proton beams are very sensitive to small differences in the patient geometry. In addition, 4D dose calculations have been performed on the same patient, but using 4D-CT(MRI) data sets based on variable breathing patterns to show the effect of possible irregular breathing on active scanned proton therapy. Using a 4D-CT(MRI), including motion irregularities, resulted in significantly different proton dose distributions. Conclusions: The authors have demonstrated that motion information from 4D-MRI can be used to generate realistic 4D-CT data sets on the basis of a single static 3D-CT data set. 4D-CT(MRI) presents a novel approach to test the robustness of treatment plans in the circumstance of patient motion.« less

EANM Dosimetry Committee Series on Standard Operational Procedures for Pre-Therapeutic Dosimetry II. Dosimetry prior to radioiodine therapy of benign thyroid diseases

Abstract: The EANM Dosimetry Committee Series “Standard Operational Procedures for Pre-Therapeutic Dosimetry” (SOP) provides advice to scientists and clinicians on how to perform patient-specific absorbed dose assessments. This particular SOP describes how to tailor the therapeutic activity to be administered for radioiodine therapy of benign thyroid diseases such as Graves’ disease or hyperthyroidism. Pretherapeutic dosimetry is based on the assessment of the individual 131I kinetics in the target tissue after the administration of a tracer activity. The present SOP makes proposals on the equipment to be used and guides the user through the measurements. Time schedules for the measurement of the fractional 131I uptake in the diseased tissue are recommended and it is shown how to calculate from these datasets the therapeutic activity necessary to administer a predefined target dose in the subsequent therapy. Potential sources of error are pointed out and the inherent uncertainties of the procedures depending on the number of measurements are discussed. The theoretical background and the derivation of the listed equations from compartment models of the iodine kinetics are explained in a supplementary file published online only.

Comparative Biodistribution and Radiation Dosimetry of 68Ga-DOTATOC and 68Ga-DOTATATE in Patients with Neuroendocrine Tumors

Abstract: 68Ga-DOTATOC and 68Ga-DOTATATE are 2 radiolabeled somatostatin analogs for in vivo diagnosis of neuroendocrine tumors with PET. The aim of the present work was to measure their comparative biodistribution and radiation dosimetry. Methods: Ten patients diagnosed with neuroendocrine tumors were included. Each patient underwent a 45-min dynamic and 3 whole-body PET/CT scans at 1, 2, and 3 h after injection of each tracer on consecutive days. Absorbed doses were calculated using OLINDA/EXM 1.1. Results: Data from 9 patients could be included in the analysis. Of the major organs, the highest uptake at 1, 2, and 3 h after injection was observed in the spleen, followed by kidneys and liver. For both tracers, the highest absorbed organ doses were seen in the spleen and urinary bladder wall, followed by kidney, adrenals, and liver. The absorbed doses to the liver and gallbladder wall were slightly but significantly higher for 68Ga-DOTATATE. The total effective dose was 0.021 ± 0.003 mSv/MBq for both tracers. Conclusion: The effective dose for a typical 100-MBq administration of 68Ga-DOTATATE and 68Ga-DOTATOC is 2.1 mSv for both tracers. Therefore, from a radiation dosimetry point of view, there is no preference for either tracer for PET/CT evaluation of somatostatin receptor–expressing tumors.

Experimental validation of the TOPAS Monte Carlo system for passive scattering proton therapy.

Abstract: Purpose: TOPAS (TOol for PArticle Simulation) is a particle simulation code recently developed with the specific aim of making Monte Carlo simulations user-friendly for research and clinical physicists in the particle therapy community. The authors present a thorough and extensive experimental validation of Monte Carlo simulations performed with TOPAS in a variety of setups relevant for proton therapy applications. The set of validation measurements performed in this work represents an overall end-to-end testing strategy recommended for all clinical centers planning to rely on TOPAS for quality assurance or patient dose calculation and, more generally, for all the institutions using passive-scattering proton therapy systems. Methods: The authors systematically compared TOPAS simulations with measurements that are performed routinely within the quality assurance (QA) program in our institution as well as experiments specifically designed for this validation study. First, the authors compared TOPAS simulations with measurements of depth-dose curves for spread-out Bragg peak (SOBP) fields. Second, absolute dosimetry simulations were benchmarked against measured machine output factors (OFs). Third, the authors simulated and measured 2D dose profiles and analyzed the differences in terms of field flatness and symmetry and usable field size. Fourth, the authors designed a simple experiment using a half-beam shifter to assess the effects of multiple Coulomb scattering, beam divergence, and inverse square attenuation on lateral and longitudinal dose profiles measured and simulated in a water phantom. Fifth, TOPAS’ capabilities to simulate time dependent beam delivery was benchmarked against dose rate functions (i.e., dose per unit time vs time) measured at different depths inside an SOBP field. Sixth, simulations of the charge deposited by protons fully stopping in two different types of multilayer Faraday cups (MLFCs) were compared with measurements to benchmark the nuclear interaction models used in the simulations. Results: SOBPs’ range and modulation width were reproduced, on average, with an accuracy of +1, −2 and ±3 mm, respectively. OF simulations reproduced measured data within ±3%. Simulated 2D dose-profiles show field flatness and average field radius within ±3% of measured profiles. The field symmetry resulted, on average in ±3% agreement with commissioned profiles. TOPAS accuracy in reproducing measured dose profiles downstream the half beam shifter is better than 2%. Dose rate function simulation reproduced the measurements within ∼2% showing that the four-dimensional modeling of the passively modulation system was implement correctly and millimeter accuracy can be achieved in reproducing measured data. For MLFCs simulations, 2% agreement was found between TOPAS and both sets of experimental measurements. The overall results show that TOPAS simulations are within the clinical accepted tolerances for all QA measurements performed at our institution. Conclusions: Our Monte Carlo simulations reproduced accurately the experimental data acquired through all the measurements performed in this study. Thus, TOPAS can reliably be applied to quality assurance for proton therapy and also as an input for commissioning of commercial treatment planning systems. This work also provides the basis for routine clinical dose calculations in patients for all passive scattering proton therapy centers using TOPAS.

Dosimetry of 223Ra-chloride: dose to normal organs and tissues

Abstract: 223Ra-Chloride (also called Alpharadin®) targets bone metastases with short range alpha particles. In recent years several clinical trials have been carried out showing, in particular, the safety and efficacy of palliation of painful bone metastases in patients with castration-resistant prostate cancer using 223Ra-chloride. The purpose of this work was to provide a comprehensive dosimetric calculation of organ doses after intravenous administration of 223Ra-chloride according to the present International Commission on Radiological Protection (ICRP) model for radium. Absorbed doses were calculated for 25 organs or tissues. Bone endosteum and red bone marrow show the highest dose coefficients followed by liver, colon and intestines. After a treatment schedule of six intravenous injections with 0.05 MBq/kg of 223Ra-chloride each, corresponding to 21 MBq for a 70 kg patient, the absorbed alpha dose to the bone endosteal cells is about 16 Gy and the corresponding absorbed dose to the red bone marrow is approximately 1.5 Gy. The comprehensive list of dose coefficients presented in this work will assist in comparing and evaluating organ doses from various therapy modalities used in nuclear medicine and will provide a base for further development of patient-specific dosimetry.

A comparative study of small field total scatter factors and dose profiles using plastic scintillation detectors and other stereotactic dosimeters: the case of the CyberKnife.

Abstract: Purpose: Small-field dosimetry is challenging, and the main limitations of most dosimeters are insufficient spatial resolution, water nonequivalence, and energy dependence. The purpose of this study was to compare plastic scintillation detectors (PSDs) to several commercial stereotactic dosimeters by measuring total scatter factors and dose profiles on a CyberKnife system. Methods: Two PSDs were developed, having sensitive volumes of 0.196 and 0.785 mm{sup 3}, and compared with other detectors. The spectral discrimination method was applied to subtract Cerenkov light from the signal. Both PSDs were compared to four commercial stereotactic dosimeters by measuring total scatter factors, namely, an IBA dosimetry stereotactic field diode (SFD), a PTW 60008 silicon diode, a PTW 60012 silicon diode, and a microLion. The measured total scatter factors were further compared with those of two independent Monte Carlo studies. For the dose profiles, two commercial detectors were used for the comparison, i.e., a PTW 60012 silicon diode and Gafchromics EBT2. Total scatter factors for a CyberKnife system were measured in circular fields with diameters from 5 to 60 mm. Dose profiles were measured for the 5- and 60-mm cones. The measurements were performed in a water tank at a 1.5-cm depth and an 80-cm source-axis distance.more » Results: The total scatter factors measured using all the detectors agreed within 1% with the Monte Carlo values for cones of 20 mm or greater in diameter. For cones of 10-20 mm in diameter, the PTW 60008 silicon diode was the only dosimeter whose measurements did not agree within 1% with the Monte Carlo values. For smaller fields (<10 mm), each dosimeter type showed different behaviors. The silicon diodes over-responded because of their water nonequivalence; the microLion and 1.0-mm PSD under-responded because of a volume-averaging effect; and the 0.5-mm PSD was the only detector within the uncertainties of the Monte Carlo simulations for all the cones. The PSDs, the PTW 60012 silicon diode, and the Gafchromics EBT2 agreed within 2% and 0.2 mm (gamma evaluation) for the measured dose profiles except in the tail of the 60-mm cone. Conclusions: Silicon diodes can be used to accurately measure small-field dose profiles but not to measure total scatter factors, whereas PSDs can be used to accurately measure both. The authors' measurements show that the use of a 1.0-mm PSD resulted in a negligible volume-averaging effect (under-response of Almost-Equal-To 1%) down to a field size of 5 mm. Therefore, PSDs are strong candidates to become reference radiosurgery detectors for beam characterization and quality assurance measurements.« less

Monte Carlo investigation of the increased radiation deposition due to gold nanoparticles using kilovoltage and megavoltage photons in a 3D randomized cell model

Abstract: Purpose: Investigation of increased radiation dose deposition due to gold nanoparticles (GNPs) using a 3D computational cell model during x-ray radiotherapy. Methods: Two GNP simulation scenarios were set up in Geant4; a single 400 nm diameter gold cluster randomly positioned in the cytoplasm and a 300 nm gold layer around the nucleus of the cell. Using an 80 kVp photon beam, the effect of GNP on the dose deposition in five modeled regions of the cell including cytoplasm, membrane, and nucleus was simulated. Two Geant4 physics lists were tested: the default Livermore and custom built Livermore/DNA hybrid physics list. 106 particles were simulated at 840 cells in the simulation. Each cell was randomly placed with random orientation and a diameter varying between 9 and 13 μm. A mathematical algorithm was used to ensure that none of the 840 cells overlapped. The energy dependence of the GNP physical dose enhancement effect was calculated by simulating the dose deposition in the cells with two energy spectra of 80 kVp and 6 MV. The contribution from Auger electrons was investigated by comparing the two GNP simulation scenarios while activating and deactivating atomic de-excitation processes in Geant4. Results: The physical dose enhancement ratio (DER) of GNP was calculated using the Monte Carlo model. The model has demonstrated that the DER depends on the amount of gold and the position of the gold cluster within the cell. Individual cell regions experienced statistically significant (p < 0.05) change in absorbed dose (DER between 1 and 10) depending on the type of gold geometry used. The DER resulting from gold clusters attached to the cell nucleus had the more significant effect of the two cases (DER ∼ 55). The DER value calculated at 6 MV was shown to be at least an order of magnitude smaller than the DER values calculated for the 80 kVp spectrum. Based on simulations, when 80 kVp photons are used, Auger electrons have a statistically insignificant (p < 0.05) effect on the overall dose increase in the cell. The low energy of the Auger electrons produced prevents them from propagating more than 250–500 nm from the gold cluster and, therefore, has a negligible effect on the overall dose increase due to GNP. Conclusions: The results presented in the current work show that the primary dose enhancement is due to the production of additional photoelectrons.

Commissioning dose computation models for spot scanning proton beams in water for a commercially available treatment planning system.

Abstract: dated by measurements. MC-generated IDDs were converted to units of Gymm 2 /MU using the measured IDDs at a depth of 2 cm employing the largest commercially available parallel-plate ionization chamber. The sensitive area of the chamber was insufficient to fully encompass the entire lateral dose deposited at depth by a pencil beam (spot). To correct for the detector size, correction factors as a function of proton energy were defined and determined using MC. The fluence of individual spots was initially modeled as a single Gaussian (SG) function and later as a double Gaussian (DG) function. The DG fluence model was introduced to account for the spot fluence due to contributions of large angle scattering from the devices within the scanning nozzle, especially from the spot profile monitor. To validate the DG fluence model, we compared calculations and measurements, including doses at the center of spread out Bragg peaks (SOBPs) as a function of nominal field size, range, and SOBP width, lateral dose profiles, and depth doses for different widths of SOBP. Dose models were validated extensively with patient treatment field-specific measurements. Results: We demonstrated that the DG fluence model is necessary for predicting the field size dependence of dose distributions. With this model, the calculated doses at the center of SOBPs as a function of nominal field size, range, and SOBP width, lateral dose profiles and depth doses for rectangular target volumes agreed well with respective measured values. With the DG fluence model for our scanning proton beam line, we successfully treated more than 500 patients from March 2010 through June 2012 with acceptable agreement between TPS calculated and measured dose distributions. However, the current dose model still has limitations in predicting field size dependence of doses at some intermediate depths of proton beams with high energies. Conclusions: We have commissioned a DG fluence model for clinical use. It is demonstrated that the DG fluence model is significantly more accurate than the SG fluence model. However, some deficiencies in modeling the low-dose envelope in the current dose algorithm still exist. Further improvements to the current dose algorithm are needed. The method presented here should be useful for commissioning pencil beam dose algorithms in new versions of TPS in the future. © 2013 American Association of Physicists in Medicine .[ http://dx.doi.org/10.1118/1.4798229]

Kidney Dosimetry in 177Lu and 90Y Peptide Receptor Radionuclide Therapy: Influence of Image Timing, Time-Activity Integration Method, and Risk Factors

Abstract: Kidney dosimetry in 177Lu and 90Y PRRT requires 3 to 6 whole-body/SPECT scans to extrapolate the peptide kinetics, and it is considered time and resource consuming. We investigated the most adequate timing for imaging and time-activity interpolating curve, as well as the performance of a simplified dosimetry, by means of just 1-2 scans. Finally the influence of risk factors and of the peptide (DOTATOC versus DOTATATE) is considered. 28 patients treated at first cycle with 177Lu DOTATATE and 30 with 177Lu DOTATOC underwent SPECT scans at 2 and 6 hours, 1, 2, and 3 days after the radiopharmaceutical injection. Dose was calculated with our simplified method, as well as the ones most used in the clinic, that is, trapezoids, monoexponential, and biexponential functions. The same was done skipping the 6 h and the 3 d points. We found that data should be collected until 100 h for 177Lu therapy and 70 h for 90Y therapy, otherwise the dose calculation is strongly influenced by the curve interpolating the data and should be carefully chosen. Risk factors (hypertension, diabetes) cause a rather statistically significant 20% increase in dose (-test, ), with DOTATATE affecting an increase of 25% compared to DOTATOC (-test, ).

Quantitative Monte Carlo–Based 90Y SPECT Reconstruction

Abstract: The evaluation of radiation absorbed doses in tumorous and healthy tissues is of increasing interest for 90Y microsphere radioembolization of liver malignancies. The objectives of this work were to introduce and validate a new reconstruction method for quantitative 90Y bremsstrahlung SPECT to improve posttreatment dosimetry. Methods: A fast Monte Carlo simulator was adapted for 90Y and incorporated into a statistical reconstruction algorithm (SPECTMC). Photon scatter and attenuation for all photons sampled from the full 90Y energy spectrum were modeled during reconstruction by Monte Carlo simulations. The energy- and distance-dependent collimator–detector response was modeled with precalculated convolution kernels. The National Electrical Manufacturers Association 2007/International Electrotechnical Commission 2008 image quality phantom was used to quantitatively evaluate the performance of SPECT-MC in comparison with those of state-of-the-art clinical SPECT reconstruction and PET. The liver radiation absorbed doses estimated by SPECT, PET, and SPECT-MC were evaluated in 5 patients consecutively treated with radioembolization. Results: In comparison with state-of-the-art clinical 90 Y SPECT reconstruction, SPECT-MC substantially improved image contrast (e.g., from 25% to 88% for the 37-mm sphere) and decreased the mean residual count error in the lung insert (from 73% to 15%) at the cost of higher image noise. Image noise and the mean count error were lower for SPECT-MC than for PET. Image contrast was higher in the larger spheres (diameter of $28 mm) but lower in the smaller spheres (#22 mm) for SPECT-MC than for PET. In the clinical study, mean absorbed dose estimates in liver regions with high absorbed doses were consistently higher for SPECT-MC than for SPECT (P 5 0.0625) and consistently higher for SPECT-MC than for PET (P 5 0.0625). Conclusion: The quantitative accuracy of 90Y bremsstrahlung SPECT is substantially improved by Monte Carlo–based modeling of the image-degrading factors. Consequently, 90Y bremsstrahlung SPECT may be used as an alternative to 90Y PET.

Experimentally studied dynamic dose interplay does not meaningfully affect target dose in VMAT SBRT lung treatments.

Abstract: Purpose: The effects of respiratory motion on the tumor dose can be divided into the gradient and interplay effects. While the interplay effect is likely to average out over a large number of fractions, it may play a role in hypofractionated [stereotactic body radiation therapy (SBRT)] treatments. This subject has been extensively studied for intensity modulated radiation therapy but less so for volumetric modulated arc therapy (VMAT), particularly in application to hypofractionated regimens. Also, no experimental study has provided full four-dimensional (4D) dose reconstruction in this scenario. The authors demonstrate how a recently described motion perturbation method, with full 4D dose reconstruction, is applied to describe the gradient and interplay effects during VMAT lung SBRT treatments. Methods: VMAT dose delivered to a moving target in a patient can be reconstructed by applying perturbations to the treatment planning system-calculated static 3D dose. Ten SBRT patients treated with 6 MV VMAT beams in five fractions were selected. The target motion (motion kernel) was approximated by 3D rigid body translation, with the tumor centroids defined on the ten phases of the 4DCT. The motion was assumed to be periodic, with the period T being an average from the empirical 4DCT respiratory trace. The real observed tumor motion (total displacement ≤8 mm) was evaluated first. Then, the motion range was artificially increased to 2 or 3 cm. Finally, T was increased to 60 s. While not realistic, making T comparable to the delivery time elucidates if the interplay effect can be observed. For a single fraction, the authors quantified the interplay effect as the maximum difference in the target dosimetric indices, most importantly the near-minimum dose (D99%), between all possible starting phases. For the three- and five-fractions, statistical simulations were performed when substantial interplay was found. Results: For the motion amplitudes and periods obtained from the 4DCT, the interplay effect is negligible (<0.2%). It is also small (0.9% average, 2.2% maximum) when the target excursion increased to 2–3 cm. Only with large motion and increased period (60 s) was a significant interplay effect observed, with D99% ranging from 16% low to 17% high. The interplay effect was statistically significantly lower for the three- and five-fraction statistical simulations. Overall, the gradient effect dominates the clinical situation. Conclusions: A novel method was used to reconstruct the volumetric dose to a moving tumor during lung SBRT VMAT deliveries. With the studied planning and treatment technique for realistic motion periods, regardless of the amplitude, the interplay has nearly no impact on the near-minimum dose. The interplay effect was observed, for study purposes only, with the period comparable to the VMAT delivery time.

Superficial dosimetry imaging based on Čerenkov emission for external beam radiotherapy with megavoltage x-ray beam.

Abstract: Purpose: Cerenkov radiation emission occurs in all tissue, when charged particles (either primary or secondary) travel at velocity above the threshold for the Cerenkov effect (about 220 KeV in tissue for electrons). This study presents the first examination of optical Cerenkov emission as a surrogate for the absorbed superficial dose for MV x-ray beams. Methods: In this study, Monte Carlo simulations of flat and curved surfaces were studied to analyze the energy spectra of charged particles produced in different regions near the surfaces when irradiated by MV x-ray beams. Cerenkov emission intensity and radiation dose were directly simulated in voxelized flat and cylindrical phantoms. The sampling region of superficial dosimetry based on Cerenkov radiation was simulated in layered skin models. Angular distributions of optical emission from the surfaces were investigated. Tissue mimicking phantoms with flat and curved surfaces were imaged with a time domain gating system. The beam field sizes (50 × 50–200 × 200 mm2), incident angles (0°–70°) and imaging regions were all varied. Results: The entrance or exit region of the tissue has nearly homogeneous energy spectra across the beam, such that their Cerenkov emission is proportional to dose. Directly simulated local intensity of Cerenkov and radiation dose in voxelized flat and cylindrical phantoms further validate that this signal is proportional to radiation dose with absolute average discrepancy within 2%, and the largest within 5% typically at the beam edges. The effective sampling depth could be tuned from near 0 up to 6 mm by spectral filtering. The angular profiles near the theoretical Lambertian emission distribution for a perfect diffusive medium, suggesting that angular correction of Cerenkov images may not be required even for curved surface. The acquisition speed and signal to noise ratio of the time domain gating system were investigated for different acquisition procedures, and the results show there is good potential for real-time superficial dose monitoring. Dose imaging under normal ambient room lighting was validated, using gated detection and a breast phantom. Conclusions: This study indicates that Cerenkov emission imaging might provide a valuable way to superficial dosimetry imaging in real time for external beam radiotherapy with megavoltage x-ray beams.

Dosimetric impact of Acuros XB deterministic radiation transport algorithm for heterogeneous dose calculation in lung cancer.

Abstract: 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 (Dm,m) and dose-to-water in medium (Dw,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 (both Dm,m and Dw,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_Dm,m, and AXB_Dw,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. AXB Dw,m and Dm,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.

On the validity of 3D polymer gel dosimetry: III. MRI-related error sources.

Abstract: In MRI (PAGAT) polymer gel dosimetry, there exists some controversy on the validity of 3D dose verifications of clinical treatments The relative contribution of important sources of uncertainty in MR scanning to the overall accuracy and precision of 3D MRI polymer gel dosimetry is quantified in this study The performance in terms of signal-to-noise and imaging artefacts was evaluated on three different MR scanners (two 15 T and a 3 T scanner) These include: (1) B0-field inhomogeneity, (2) B1-field inhomogeneity, (3) dielectric effects (losses and standing waves) and (4) temperature inhomogeneity during scanning B0-field inhomogeneities that amount to maximum 5 ppm result in dose deviations of up to 43% and deformations of up to 5 pixels Compensation methods are proposed B1-field inhomogeneities were found to induce R2 variations in large anthropomorphic phantoms both at 15 and 3 T At 15 T these effects are mainly caused by the coil geometry resulting in dose deviations of up to 25% After the correction of the R2 maps using a heuristic flip angle–R2 relation, these dose deviations are reduced to 24% At 3 T, the dielectric properties of the gel phantoms are shown to strongly influence B1-field homogeneity, hence R2 homogeneity, especially of large anthropomorphic phantoms The low electrical conductivity of polymer gel dosimeters induces standing wave patterns resulting in dose deviations up to 50% Increasing the conductivity of the gel by adding NaCl reduces the dose deviation to 25% after which the post-processing is successful in reducing the remaining inhomogeneities caused by the coil geometry to within 24% The measurements are supported by computational modelling of the B1-field Finally, temperature fluctuations of 1 °C frequently encountered in clinical MRI scanners result in dose deviations up to 15% It is illustrated that with adequate temperature stabilization, the dose uncertainty is reduced to within 258%

Superficial dosimetry imaging of Čerenkov emission in electron beam radiotherapy of phantoms

Abstract: Cerenkov emission is generated from ionizing radiation in tissue above 264 keV energy. This study presents the first examination of this optical emission as a surrogate for the absorbed superficial dose. Cerenkov emission was imaged from the surface of flat tissue phantoms irradiated with electrons, using a range of field sizes from 6 cm × 6 cm to 20 cm × 20 cm, incident angles from 0° to 50°, and energies from 6 to 18 MeV. The Cerenkov images were compared with the estimated superficial dose in phantoms from direct diode measurements, as well as calculations by Monte Carlo and the treatment planning system. Intensity images showed outstanding linear agreement (R(2) = 0.97) with reference data of the known dose for energies from 6 to 18 MeV. When orthogonal delivery was carried out, the in-plane and cross-plane dose distribution comparisons indicated very little difference (± 2-4% differences) between the different methods of estimation as compared to Cerenkov light imaging. For an incident angle 50°, the Cerenkov images and Monte Carlo simulation show excellent agreement with the diode data, but the treatment planning system had a larger error (OPT = ± 1~2%, diode = ± 2~3%, TPS = ± 6-8% differences) as would be expected. The sampling depth of superficial dosimetry based on Cerenkov radiation has been simulated in a layered skin model, showing the potential of sampling depth tuning by spectral filtering. Taken together, these measurements and simulations indicate that Cerenkov emission imaging might provide a valuable method of superficial dosimetry imaging from incident radiotherapy beams of electrons.

NAIRAS aircraft radiation model development, dose climatology, and initial validation

Abstract: [1] The Nowcast of Atmospheric Ionizing Radiation for Aviation Safety (NAIRAS) is a real-time, global, physics-based model used to assess radiation exposure to commercial aircrews and passengers. The model is a free-running physics-based model in the sense that there are no adjustment factors applied to nudge the model into agreement with measurements. The model predicts dosimetric quantities in the atmosphere from both galactic cosmic rays (GCR) and solar energetic particles, including the response of the geomagnetic field to interplanetary dynamical processes and its subsequent influence on atmospheric dose. The focus of this paper is on atmospheric GCR exposure during geomagnetically quiet conditions, with three main objectives. First, provide detailed descriptions of the NAIRAS GCR transport and dosimetry methodologies. Second, present a climatology of effective dose and ambient dose equivalent rates at typical commercial airline altitudes representative of solar cycle maximum and solar cycle minimum conditions and spanning the full range of geomagnetic cutoff rigidities. Third, conduct an initial validation of the NAIRAS model by comparing predictions of ambient dose equivalent rates with tabulated reference measurement data and recent aircraft radiation measurements taken in 2008 during the minimum between solar cycle 23 and solar cycle 24. By applying the criterion of the International Commission on Radiation Units and Measurements (ICRU) on acceptable levels of aircraft radiation dose uncertainty for ambient dose equivalent greater than or equal to an annual dose of 1 mSv, the NAIRAS model is within 25% of the measured data, which fall within the ICRU acceptable uncertainty limit of 30%. The NAIRAS model predictions of ambient dose equivalent rate are generally within 50% of the measured data for any single-point comparison. The largest differences occur at low latitudes and high cutoffs, where the radiation dose level is low. Nevertheless, analysis suggests that these single-point differences will be within 30% when a new deterministic pion-initiated electromagnetic cascade code is integrated into NAIRAS, an effort which is currently underway.

Performance parameters of a liquid filled ionization chamber array

Abstract: Purpose: In this work, the properties of the two-dimensional liquid filled ionization chamber array Octavius 1000SRS (PTW-Freiburg, Germany) for use in clinical photon-beam dosimetry are investigated. Methods: Measurements were carried out at an Elekta Synergy and Siemens Primus accelerator. For measurements of stability, linearity, and saturation effects of the 1000SRS array a Semiflex 31013 ionization chamber (PTW-Freiburg, Germany) was used as a reference. The effective point of measurement was determined by TPR measurements of the array in comparison with a Roos chamber (type 31004, PTW-Freiburg, Germany). The response of the array with varying field size and depth of measurement was evaluated using a Semiflex 31010 ionization chamber as a reference. Output factor measurements were carried out with a Semiflex 31010 ionization chamber, a diode (type 60012, PTW-Freiburg, Germany), and the detector array under investigation. The dose response function for a single detector of the array was determined by measuring 1 cm wide slit-beam dose profiles and comparing them against diode-measured profiles. Theoretical aspects of the low pass properties and of the sampling frequency of the detector array were evaluated. Dose profiles measured with the array and the diode detector were compared, and an intensity modulated radiation therapy (IMRT) field was verified using the Gamma-Index method and the visualization of line dose profiles. Results: The array showed a short and long term stability better than 0.1% and 0.2%, respectively. Fluctuations in linearity were found to be within ±0.2% for the vendor specified dose range. Saturation effects were found to be similar to those reported in other studies for liquid-filled ionization chambers. The detector's relative response varied with field size and depth of measurement, showing a small energy dependence accounting for maximum signal deviations of ±2.6% from the reference condition for the setup used. The σ-values of the Gaussian dose response function for a single detector of the array were found to be (0.72 ± 0.25) mm at 6 MV and (0.74 ± 0.25) mm at 15 MV and the corresponding low pass cutoff frequencies are 0.22 and 0.21 mm−1, respectively. For the inner 5 × 5 cm2 region and the outer 11 × 11 cm2 region of the array the Nyquist theorem is fulfilled for maximum sampling frequencies of 0.2 and 0.1 mm−1, respectively. An IMRT field verification with a Gamma-Index analysis yielded a passing rate of 95.2% for a 3 mm/3% criterion with a TPS calculation as reference. Conclusions: This study shows the applicability of the Octavius 1000SRS in modern dosimetry. Output factor and dose profile measurements illustrated the applicability of the array in small field and stereotactic dosimetry. The high spatial resolution ensures adequate measurements of dose profiles in regular and intensity modulated photon-beam fields.