Technical note: Temporal and thermal stability of optical response for silicone-based 3D radiochromic dosimeters.

BACKGROUND: Radiochromic silicone-based dosimeters are flexible 3D dosimeters, which at appropriate concentration of leucomalachite green (LMG) and curing agent are dose-rate independent for clinical photon beams. However, their dose response is based on chemical processes that can be influenced by temporal and thermal conditions, impacting measurement stability. PURPOSE: The aim of this study was to investigate the temporal stability of the dose response of radiochromic dosimeters for different curing times and post-irradiation storage temperatures. METHODS: Six cylindrical dosimeters (5 cm diameter, 5 cm length) were produced in a single batch and separated into two groups that were irradiated 72 and 118 h after production. The same photon plan, consisting of two 10 × 1.6 cm2 opposing fields, was delivered to all dosimeters. After irradiation, the dosimeters were separated into three groups, stored at 5°C, 15°C, and 20°C, and read out for five consecutive days. RESULTS: Storage temperature influenced the measurement stability, and changes in the optical response with time differed between irradiated and non-irradiated parts of the dosimeters. The relative change between signal and background was greater than 10% for all measurements performed 24 h or more after irradiation, except for dosimeters stored at 5°C, which changed by 2%-5% after 24 h. The dosimeter temporal stability was not influenced by curing time. CONCLUSIONS: For room temperature storage (15°C and 20°C), readout should take place as soon as possible after irradiation since the background color increased rapidly for both curing times (72 and 118 h), whereas the dosimeters are stored at 5°C, readout can be performed up to 24 h after.

Empirical quenching correction in radiochromic silicone-based three-dimensional dosimetry of spot-scanning proton therapy.

BACKGROUND AND PURPOSE: Three-dimensional dosimetry of proton therapy (PT) with chemical dosimeters is challenged by signal quenching, which is a lower dose-response in regions with high ionization density due to high linear-energy-transfer (LET) and dose-rate. This study aimed to assess the viability of an empirical correction model for 3D radiochromic silicone-based dosimeters irradiated with spot-scanning PT, by parametrizing its LET and dose-rate dependency. MATERIALS AND METHODS: Ten cylindrical radiochromic dosimeters (Ø50 and Ø75 mm) were produced in-house, and irradiated with different spot-scanning proton beam configurations and machine-set dose rates ranging from 56 to 145 Gy/min. Beams with incident energies of 75, 95 and 120 MeV, a spread-out Bragg peak and a plan optimized to an irregular target volume were included. Five of the dosimeters, irradiated with 120 MeV beams, were used to estimate the quenching correction factors. Monte Carlo simulations were used to obtain dose and dose-averaged-LET (LETd) maps. Additionally, a local dose-rate map was estimated, using the simulated dose maps and the machine-set dose-rate information retrieved from the irradiation log-files. Finally, the correction factor was estimated as a function of LETd and local dose-rate and tested on the different fields. RESULTS: Gamma-pass-rates of the corrected measurements were >94% using a 3%-3 mm gamma analysis and >88% using 2%-2 mm, with a dose deviation of <5.6 ± 1.8%. Larger dosimeters showed a 20% systematic increase in dose-response, but the same quenching in signal when compared to the smaller dosimeters. CONCLUSION: The quenching correction model was valid for different dosimeter sizes to obtain relative dosimetric maps of complex dose distributions in PT.

Impact of interfractional target motion in locally advanced cervical cancer patients treated with spot scanning proton therapy using an internal target volume strategy.

BACKGROUND AND PURPOSE: The more localized dose deposition of proton therapy (PT) compared to photon therapy might allow a reduction in treatment-related side effects but induces additional challenges to address. The aim of this study was to evaluate the impact of interfractional motion on the target and organs at risk (OARs) in cervical cancer patients treated with spot scanning PT using an internal target volume (ITV) strategy. METHODS AND MATERIALS: For ten locally advanced cervical cancer patients, empty and full bladder planning computed tomography (pCT) as well as 25 daily cone beam CTs (CBCTs) were available. The Clinical Target Volume (CTV), the High Risk CTV (CTVHR) (gross tumor volume and whole cervix), the non-involved uterus as well as the OARs (bowel, bladder and rectum) were contoured on the daily CBCTs and transferred to the pCT through rigid bony match. Using synthetic CTs derived from pCTs, four-beam spot scanning PT plans were generated to target the patient-specific ITV with 45 Gy(RBE) in 25 fractions. This structure was defined based on pre-treatment MRI and CT to anticipate potential target motion throughout the treatment. D98% of the targets and V40Gy(RBE) of the OARs were extracted from the daily anatomies, accumulated and analyzed. In addition, the impact of bladder volume deviations from planning values on target and bowel dose was investigated. RESULTS: The ITV strategy ensured a total accumulated dose >42.75 Gy(RBE) to the CTVHR for all ten patients. Two patients with large bladder-related uterus motion had accumulated dose to the non-involved uterus of 35.7 Gy(RBE) and 41.1 Gy(RBE). Variations in bowel V40Gy(RBE) were found to be correlated (Pearson r = -0.55; p-value <0.0001) with changes in bladder volume during treatment. CONCLUSION: The ITV concept ensured adequate dose to the CTVHR, but was insufficient for the non-involved uterus of patients subject to large target interfractional motion. CBCT monitoring and occasional replanning is recommended along the same lines as with photon radiotherapy in cervical cancer.

Hepatic regeneration following radiation-induced liver injury is associated with increased hepatobiliary secretion measured by PET in Göttingen minipigs.

Normal liver tissue is highly vulnerable towards irradiation, which remains a challenge in radiotherapy of hepatic tumours. Here, we examined the effects of radiation-induced liver injury on two specific liver functions and hepatocellular regeneration in a minipig model. Five Göttingen minipigs were exposed to whole-liver stereotactic body radiation therapy (SBRT) in one fraction (14 Gy) and examined 4-5 weeks after; five pigs were used as controls. All pigs underwent in vivo positron emission tomography (PET) studies of the liver using the conjugated bile acid tracer [N-methyl-11C]cholylsarcosine ([11C]CSar) and the galactose-analogue tracer [18F]fluoro-2-deoxy-D-galactose ([18F]FDGal). Liver tissue samples were evaluated histopathologically and by immunohistochemical assessment of hepatocellular mitosis, proliferation and apoptosis. Compared with controls, both the rate constant for secretion of [11C]CSar from hepatocytes into intrahepatic bile ducts as well as back into blood were doubled in irradiated pigs, which resulted in reduced residence time of [11C]CSar inside the hepatocytes. Also, the hepatic systemic clearance of [18F]FDGal in irradiated pigs was slightly increased, and hepatocellular regeneration was increased by a threefold. In conclusion, parenchymal injury and increased regeneration after whole-liver irradiation was associated with enhanced hepatobiliary secretion of bile acids. Whole-liver SBRT in minipigs ultimately represents a potential large animal model of radiation-induced liver injury and for testing of normal tissue protection methods.

Evaluation of an a priori scatter correction algorithm for cone-beam computed tomography based range and dose calculations in proton therapy.

BACKGROUND AND PURPOSE: Scatter correction of cone-beam computed tomography (CBCT) projections may enable accurate online dose-delivery estimations in photon and proton-based radiotherapy. This study aimed to evaluate the impact of scatter correction in CBCT-based proton range/dose calculations, in scans acquired in both proton and photon gantries. MATERIAL AND METHODS: CBCT projections of a Catphan and an Alderson phantom were acquired on both a proton and a photon gantry. The scatter corrected CBCTs (corrCBCTs) and the clinical reconstructions (stdCBCTs) were compared against CTs rigidly registered to the CBCTs (rigidCTs). The CBCTs of the Catphan phantom were segmented by materials for CT number analysis. Water equivalent path length (WEPL) maps were calculated through the Alderson phantom while proton plans optimized on the rigidCT and recalculated on all CBCTs were compared in a gamma analysis. RESULTS: In medium and high-density materials, the corrCBCT CT numbers were much closer to those of the rigidCT than the stdCBCTs. E.g. in the 50% bone segmentations the differences were reduced from above 300 HU (with stdCBCT) to around 60-70 HU (with corrCBCT). Differences in WEPL from the rigidCT were typically well below 5 mm for the corrCBCTs, compared to well above 10 mm for the stdCBCTs with the largest deviations in the head and thorax regions. Gamma pass rates (2%/2mm) when comparing CBCT-based dose re-calculations to rigidCT calculations were improved from around 80% (with stdCBCT) to mostly above 90% (with corrCBCT). CONCLUSION: Scatter correction leads to substantial artefact reductions, improving accuracy of CBCT-based proton range/dose calculations.

Comparison of single and dual energy CT for stopping power determination in proton therapy of head and neck cancer.

BACKGROUND AND PURPOSE: Patients with head and neck (HN) cancer may benefit from proton therapy due to the potential for sparing of normal tissue. For planning of proton therapy, dual-energy CT (DECT) has been shown to provide superior stopping power ratio (SPR) determination in phantom materials and organic tissue samples, compared to single-energy CT (SECT). However, the benefit of DECT in HN cancer patients has not yet been investigated. This study therefore compared DECT- and SECT-based SPR estimation for HN cancer patients. MATERIALS AND METHODS: Fourteen HN cancer patients were DECT scanned. Eight patients were scanned using a dual source DECT scanner and six were scanned with a conventional SECT scanner by acquiring two consecutive scans. SECT image sets were computed as a weighted summation of the low and high energy DECT image sets. DECT- and SECT-based SPR maps were derived. Water-equivalent path lengths (WEPLs) through the SPR maps were compared in the eight cases with dual source DECT scans. Mean SPR estimates over region-of-interests (ROIs) in the cranium, brain and eyes were analyzed for all patients. RESULTS: A median WEPL difference of 1.9 mm (1.5%) was found across the eight patients. Statistically significant SPR differences were seen for the ROIs in the brain and eyes, with the SPR estimates based on DECT overall lower than for SECT. CONCLUSIONS: Clinically relevant WEPL and SPR differences were found between DECT and SECT, which could imply that the accuracy of treatment planning for proton therapy would benefit from DECT-based SPR estimation.

Impact of bowel gas and body outline variations on total accumulated dose with intensity-modulated proton therapy in locally advanced cervical cancer patients.

INTRODUCTION: Density changes occurring during fractionated radiotherapy in the pelvic region may degrade proton dose distributions. The aim of the study was to quantify the dosimetric impact of gas cavities and body outline variations. MATERIAL AND METHODS: Seven patients with locally advanced cervical cancer (LACC) were analyzed through a total of 175 daily cone beam computed tomography (CBCT) scans. Four-beams intensity-modulated proton therapy (IMPT) dose plans were generated targeting the internal target volume (ITV) composed of: primary tumor, elective and pathological nodes. The planned dose was 45 Gy [Relative-Biological-Effectiveness-weighted (RBE)] in 25 fractions and simultaneously integrated boosts of pathologic lymph nodes were 55-57.5 Gy (RBE). In total, 475 modified CTs were generated to evaluate the effect of: 1/gas cavities, 2/outline variations and 3/the two combined. The anatomy of each fraction was simulated by propagating gas cavities contours and body outlines from each daily CBCT to the pCT. Hounsfield units corresponding to gas and fat were assigned to the propagated contours. D98 (least dose received by the hottest 98% of the volume) and D99.9 for targets and V43Gy(RBE) (volume receiving ≥43 Gy(RBE)) for organs at risk (OARs) were recalculated on each modified CT, and total dose was evaluated through dose volume histogram (DVH) addition across all fractions. RESULTS: Weight changes during radiotherapy were between -3.1% and 1.2%. Gas cavities and outline variations induced a median [range] dose degradation for ITV45 of 1.0% [0.5-3.5%] for D98 and 2.1% [0.8-6.4%] for D99.9. Outline variations had larger dosimetric impact than gas cavities. Worst nodal dose degradation was 2.0% for D98 and 2.3% for D99.9. The impact on bladder, bowel and rectum was limited with V43Gy(RBE) variations ≤3.5 cm3. CONCLUSION: Bowel gas cavities and outline variations had minor impact on accumulated dose in targets and OAR of four-field IMPT in a LACC population of moderate weight changes.

Functional image-guided dose escalation in gliomas using of state-of-the-art photon vs. proton therapy.

BACKGROUND: Recurrences of glioma are usually local, suggesting the need for higher tumor dose. We investigated the boundaries for dose escalation of an 18F-fluoro-ethyl-tyrosine positron emission tomography defined target by intensity-modulated photon therapy (IMRT), volumetric modulated arc therapy (VMAT) and intensity-modulated proton therapy (IMPT). MATERIALS AND METHODS: Standard dose (60 Gy) and dose-escalated plans were calculated for seven patients using IMRT, VMAT and IMPT. The achieved boost dose, the dose to the organs at risk (OAR), the dose homogeneity (defined as overdose volume, ODV) and the ratio of the 30 Gy isodose curve and the boost volume (R30) were compared. The risk of radionecrosis was estimated using the ratio of the dose volume histograms of the brain (range 30-60 Gy). RESULTS: The mean boost dose was 77.1 Gy for IMRT, 79.2 Gy for VMAT and 85.1 GyE for IMPT. Compared with the standard plan, the ODV was unchanged and the R30 increased (17%) for IMRT. For VMAT, the ODV decreased (7%) and the R30 was unchanged whereas IMPT substantially decreased ODV (61%), R30 (22%), OAR doses as well as the risk of radionecrosis. CONCLUSIONS: Dose escalation can be achieved with IMRT, VMAT and IMPT while respecting normal tissue constraints, yet with IMPT being most favorable.

Technical Note: Improving proton stopping power ratio determination for a deformable silicone-based 3D dosimeter using dual energy CT.

PURPOSE: The aim of this study was to investigate whether the stopping power ratio (SPR) of a deformable, silicone-based 3D dosimeter could be determined more accurately using dual energy (DE) CT compared to using conventional methods based on single energy (SE) CT. The use of SECT combined with the stoichiometric calibration method was therefore compared to DECT-based determination. METHODS: The SPR of the dosimeter was estimated based on its Hounsfield units (HUs) in both a SECT image and a DECT image set. The stoichiometric calibration method was used for converting the HU in the SECT image to a SPR value for the dosimeter while two published SPR calibration methods for dual energy were applied on the DECT images. Finally, the SPR of the dosimeter was measured in a 60 MeV proton by quantifying the range difference with and without the dosimeter in the beam path. RESULTS: The SPR determined from SECT and the stoichiometric method was 1.10, compared to 1.01 with both DECT calibration methods. The measured SPR for the dosimeter material was 0.97. CONCLUSIONS: The SPR of the dosimeter was overestimated by 13% using the stoichiometric method and by 3% when using DECT. If the stoichiometric method should be applied for the dosimeter, the HU of the dosimeter must be manually changed in the treatment planning system in order to give a correct SPR estimate. Using a wrong SPR value will cause differences between the calculated and the delivered treatment plans.

A simulation study on proton computed tomography (CT) stopping power accuracy using dual energy CT scans as benchmark.

BACKGROUND: Accurate stopping power estimation is crucial for treatment planning in proton therapy, and the uncertainties in stopping power are currently the largest contributor to the employed dose margins. Dual energy x-ray computed tomography (CT) (clinically available) and proton CT (in development) have both been proposed as methods for obtaining patient stopping power maps. The purpose of this work was to assess the accuracy of proton CT using dual energy CT scans of phantoms to establish reference accuracy levels. MATERIAL AND METHODS: A CT calibration phantom and an abdomen cross section phantom containing inserts were scanned with dual energy and single energy CT with a state-of-the-art dual energy CT scanner. Proton CT scans were simulated using Monte Carlo methods. The simulations followed the setup used in current prototype proton CT scanners and included realistic modeling of detectors and the corresponding noise characteristics. Stopping power maps were calculated for all three scans, and compared with the ground truth stopping power from the phantoms. RESULTS: Proton CT gave slightly better stopping power estimates than the dual energy CT method, with root mean square errors of 0.2% and 0.5% (for each phantom) compared to 0.5% and 0.9%. Single energy CT root mean square errors were 2.7% and 1.6%. Maximal errors for proton, dual energy and single energy CT were 0.51%, 1.7% and 7.4%, respectively. CONCLUSION: Better stopping power estimates could significantly reduce the range errors in proton therapy, but requires a large improvement in current methods which may be achievable with proton CT.

An adaptive radiotherapy planning strategy for bladder cancer using deformation vector fields.

PURPOSE: Adaptive radiotherapy (ART) has considerable potential in treatment of bladder cancer due to large inter-fractional changes in shape and size of the target. The aim of this study was to compare our clinically applied method for plan library creation that involves manual bladder delineations (Clin-ART) with a method using the deformation vector fields (DVFs) resulting from intensity-based deformable image registrations (DVF-based ART). MATERIALS AND METHODS: The study included thirteen patients with urinary bladder cancer who had daily cone beam CTs (CBCTs) acquired for set-up. In both ART strategies investigated, three plan selection volumes were generated using the CBCTs from the first four fractions; in Clin-ART boolean combinations of delineated bladders were used, while the DVF-based strategy applied combinations of the mean and standard deviation of patient-specific DVFs. The volume ratios (VRs) of the course-averaged PTV for the two ART strategies relative the non-adaptive PTV were calculated. RESULTS: Both Clin-ART and DVF-based ART considerably reduced the course-averaged PTV, compared to non-adaptive RT. The VR for DVF-based ART was lower than for Clin-ART (0.65 vs. 0.73; p<0.01). CONCLUSIONS: DVF-based ART for bladder irradiation has a considerable normal tissue sparing potential surpassing our already highly conformal clinically applied ART strategy.

Improved proton computed tomography by dual modality image reconstruction.

PURPOSE: Proton computed tomography (CT) is a promising image modality for improving the stopping power estimates and dose calculations for particle therapy. However, the finite range of about 33 cm of water of most commercial proton therapy systems limits the sites that can be scanned from a full 360° rotation. In this paper the authors propose a method to overcome the problem using a dual modality reconstruction (DMR) combining the proton data with a cone-beam x-ray prior. METHODS: A Catphan 600 phantom was scanned using a cone beam x-ray CT scanner. A digital replica of the phantom was created in the Monte Carlo code Geant4 and a 360° proton CT scan was simulated, storing the entrance and exit position and momentum vector of every proton. Proton CT images were reconstructed using a varying number of angles from the scan. The proton CT images were reconstructed using a constrained nonlinear conjugate gradient algorithm, minimizing total variation and the x-ray CT prior while remaining consistent with the proton projection data. The proton histories were reconstructed along curved cubic-spline paths. RESULTS: The spatial resolution of the cone beam CT prior was retained for the fully sampled case and the 90° interval case, with the MTF = 0.5 (modulation transfer function) ranging from 5.22 to 5.65 linepairs/cm. In the 45° interval case, the MTF = 0.5 dropped to 3.91 linepairs/cm For the fully sampled DMR, the maximal root mean square (RMS) error was 0.006 in units of relative stopping power. For the limited angle cases the maximal RMS error was 0.18, an almost five-fold improvement over the cone beam CT estimate. CONCLUSIONS: Dual modality reconstruction yields the high spatial resolution of cone beam x-ray CT while maintaining the improved stopping power estimation of proton CT. In the case of limited angles, the use of prior image proton CT greatly improves the resolution and stopping power estimate, but does not fully achieve the quality of a 360° proton CT scan.