New optically stimulated luminescence dosimetry film optimized for energy dependence guided by Monte Carlo simulations.

Optically stimulated luminescence (OSL) film dosimeters, based on BaFBr:Eu2+phosphor material, have major dosimetric advantages such as dose linearity, high spatial resolution, film re-usability, and immediate film readout. However, they exhibit an energy-dependent over-response at low photon energies because they are not made of tissue-equivalent materials. In this work, the OSL energy-dependent response was optimized by lowering the phosphor grain size and seeking an optimal choice of phosphor concentration and film thickness to achieve sufficient signal sensitivity. This optimization process combines measurement-based assessments of energy response in narrow x-ray beams with various energy response calculation methods applied to different film metrics. Theoretical approaches and MC dose simulations were used for homogeneous phosphor distributions and for isolated phosphor grains of different dimensions, where the dose in the phosphor grain was calculated. In total 8 OSL films were manufactured with different BaFBr:Eu2+median particle diameters (D50): 3.2µm, 1.5µm and 230 nm and different phosphor concentrations (1.6%, 5.3% and 21.3 %) and thicknesses (from 5.2 to 49µm). Films were irradiated in narrow x-ray spectra (N60, N80, N-150 and N-300) and the signal intensity relative to the nominal dose-to-water value was normalized to Co-60. Finally, we experimentally tested the response of several films in Varian 6MV TrueBeam STx linear accelerator using the following settings: 10 × 10 cm2field, 0deggantry angle, 90 cm SSD, 10 cm depth. The x-ray irradiation experiment reported a reduced energy response for the smallest grain size with an inverse correlation between response and grain size. The N-60 irradiation showed a 43% reduction in the energy over-response when going from 3µm to 230 nm grain size for the 5% phosphor concentration. Energy response calculation using a homogeneous dispersion of the phosphor underestimated the experimental response and was not able to obtain the experimental correlation between grain size and energy response. Isolated grain size modeling combined with MC dose simulations allowed to establish a good agreement with experimental data, and enabled steering the production of optimized OSL-films. The clinical 6 MV beam test confirmed a reduction in energy dependence, which is visible in small-grain films where a decrease in out-of-field over-response was observed.

Mini-GRID radiotherapy on the CLEAR very-high-energy electron beamline: collimator optimization, film dosimetry, and Monte Carlo simulations.

Objective.Spatially-fractionated radiotherapy (SFRT) delivered with a very-high-energy electron (VHEE) beam and a mini-GRID collimator was investigated to achieve synergistic normal tissue-sparing through spatial fractionation and the FLASH effect.Approach.A tungsten mini-GRID collimator for delivering VHEE SFRT was optimized using Monte Carlo (MC) simulations. Peak-to-valley dose ratios (PVDRs), depths of convergence (DoCs, PVDR ≤ 1.1), and peak and valley doses in a water phantom from a simulated 150 MeV VHEE source were evaluated. Collimator thickness, hole width, and septal width were varied to determine an optimal value for each parameter that maximized PVDR and DoC. The optimized collimator (20 mm thick rectangular prism with a 15 mm × 15 mm face with a 7 × 7 array of 0.5 mm holes separated by 1.1 mm septa) was 3D-printed and used for VHEE irradiations with the CERN linear electron accelerator for research beam. Open beam and mini-GRID irradiations were performed at 140, 175, and 200 MeV and dose was recorded with radiochromic films in a water tank. PVDR, central-axis (CAX) and valley dose rates and DoCs were evaluated.Main results.Films demonstrated peak and valley dose rates on the order of 100 s of MGy/s, which could promote FLASH-sparing effects. Across the three energies, PVDRs of 2-4 at 13 mm depth and DoCs between 39 and 47 mm were achieved. Open beam and mini-GRID MC simulations were run to replicate the film results at 200 MeV. For the mini-GRID irradiations, the film CAX dose was on average 15% higher, the film valley dose was 28% higher, and the film PVDR was 15% lower than calculated by MC.Significance.Ultimately, the PVDRs and DoCs were determined to be too low for a significant potential for SFRT tissue-sparing effects to be present, particularly at depth. Further beam delivery optimization and investigations of new means of spatial fractionation are warranted.

Experimental validation of Monte Carlo dosimetry for therapeutic beta emitters with radiochromic film in a 3D-printed phantom.

PURPOSE: Validation of dosimetry software, such as Monte Carlo (MC) radiation transport codes used for patient-specific absorbed dose estimation, is critical prior to their use in clinical decision making. However, direct experimental validation in the clinic is generally not performed for low/medium-energy beta emitters used in radiopharmaceutical therapy (RPT) due to the challenges of measuring energy deposited by short-range particles. Our objective was to design a practical phantom geometry for radiochromic film (RF)-based absorbed dose measurements of beta-emitting radionuclides and perform experiments to directly validate our in-house developed Dose Planning Method (DPM) MC code dedicated to internal dosimetry. METHODS: The experimental setup was designed for measuring absorbed dose from beta emitters that have a range sufficiently penetrating to ∼200 µm in water as well as to capture any photon contributions to absorbed dose. Assayed 177 Lu and 90 Y liquid sources, 13-450 MBq estimated to deliver 0.5-10 Gy to the sensitive layer of the RF, were injected into the cavity of two 3D-printed half-cylinders that had been sealed with 12.7 µm or 25.4 µm thick Kapton Tape. A 3.8 × 6 cm strip of GafChromic EBT3 RF was sandwiched between the two taped half-cylinders. After 2-48 h exposures, films were retrieved and wipe tested for contamination. Absorbed dose to the RF was measured using a commercial triple-channel dosimetry optimization method and a calibration generated via 6 MV photon beam. Profiles were analyzed across the central 1 cm2 area of the RF for validation. Eleven experiments were completed with 177 Lu and nine with 90 Y both in saline and a bone equivalent solution. Depth dose curves were generated for 177 Lu and 90 Y stacking multiple RF strips between a single filled half-cylinder and an acrylic backing. All experiments were modeled in DPM to generate voxelized MC absorbed dose estimates. We extended our study to benchmark general purpose MC codes MCNP6 and EGSnrc against the experimental results as well. RESULTS: A total of 20 experiments showed that both the 3D-printed phantoms and the final absorbed dose values were reproducible. The agreement between the absorbed dose estimates from the RF measurements and DPM was on average -4.0% (range -10.9% to 3.2%) for all single film 177 Lu experiments and was on average -1.0% (range -2.7% to 0.7%) for all single film 90 Y experiments. Absorbed depth dose estimates by DPM agreed with RF on average 1.2% (range -8.0% to 15.2%) across all depths for 177 Lu and on average 4.0% (range -5.0% to 9.3%) across all depths for 90 Y. DPM absorbed dose estimates agreed with estimates from EGSnrc and MCNP across the board, within 4.7% and within 3.4% for 177 Lu and 90 Y respectively, for all geometries and across all depths. MC showed that absorbed dose to RF from betas was greater than 92% of the total (betas + other radiations) for 177 Lu, indicating measurement of dominant beta contribution with our design. CONCLUSIONS: The reproducible results with a RF insert in a simple phantom designed for liquid sources demonstrate that this is a reliable setup for experimentally validating dosimetry algorithms used in therapies with beta-emitting unsealed sources. Absorbed doses estimated with the DPM MC code showed close agreement with RF measurement and with results from two general purpose MC codes, thereby validating the use of this algorithms for clinical RPT dosimetry.

Microdosimetric modeling of the sensitometric curve of GafChromic films in the photon fields.

The sensitometric curve of EBT3 GafChromic film located at a depth 5 cm of RW3 water-equivalent phantom exposed to 6 MV X-rays is investigated. Variation of optical density for absorbed doses less than 2 Gy is determined by the experimental measurement together with the microdosimetric one-hit detector model. It is found that this model needs two fitting parameters, a maximum optical density and a saturation parameter. Both of them depend on the film structure as well as the photon spectrum. Meanwhile, the saturation parameter is a function of the microdosimetric single-event distribution of specific energy, i.e., f1(z). To calculate this distribution, irradiation of the films is simulated by the Geant4 toolkit. A sample of EBT3 film with 2 mm × 2 mm area is simulated. Active layer of the film is considered to contain 5000 cylindrical sensitive targets (SV) with 9.4 µm length and 1.62 µm diameter, located in random positions with different axial directions. The results obtained show that below an absorbed dose 1 Gy the maximum difference between the measured and the calculated optical densities is about 11%, while for the doses above 1 Gy the discrepancy is at most 3%. Eventually, it can be concluded that the sensitometric curve of EBT3 GafChromic film can satisfactorily be determined by the microdosimetric one-hit detector model, especially in the dose range above 1 Gy.

Electron Scattering in Conventional Cell Flask Experiments and Dose Distribution Dependency.

Electron beam therapy (EBT) is commonly used for treating superficial and subdermal tumors. Previous cellular radiosensitivity research using EBT may be underestimating the contribution from flask wall scattering and the corresponding dose distribution. Single cell suspensions of Chinese hamster ovary (CHO) cells were plated on flasks and irradiated with 3, 4, 7, 9, and 18 MeV energy electron beams from two different institutions, and the spatial locations of surviving colonies were recorded. Gafchromic film dosimetry and Monte Carlo simulations were carried out to determine the spatial electron scattering contribution from the flask walls. Low electron irradiation resulted in an uneven surviving colony distribution concentrated near the periphery of the flasks, while spatial colony formation was statistically uniform at energies above 7 MeV. Our data demonstrates that without proper dosimetric corrections, studies using low energy electrons can lead to misinterpretations of energy dependent cellular radiosensitivity in culture vessels, and radiotherapeutic applications.

Verification of a Monte Carlo dose calculation engine in proton minibeam radiotherapy in a passive scattering beamline for preclinical trials.

OBJECTIVES: Proton minibeam radiation therapy (pMBRT) is a novel therapeutic strategy that combines the benefits of proton therapy with the remarkable normal tissue preservation observed with the use of submillimetric spatially fractionated beams. This promising technique has been implemented at the Institut Curie-Proton therapy centre (ICPO) using a first prototype of a multislit collimator. The purpose of this work was to develop a Monte Carlo-based dose calculation engine to reliably guide preclinical studies at ICPO. METHODS: The whole "Y1"-passive beamline at the ICPO, including pMBRT implementation, was modelled using the Monte Carlo GATE v. 7.0 code. A clinically relevant proton energy (100 MeV) was used as starting point. Minibeam generation by means of the brass collimator used in the first experiments was modelled. A virtual source was modelled at the exit of the beamline nozzle and outcomes were compared with dosimetric measurements performed with EBT3 gafchromic films and a diamond detector in water. Dose distributions were recorded in a water phantom and in rat CT images (7-week-old male Fischer rats). RESULTS: The dose calculation engine was benchmarked against experimental data and was then used to assess dose distributions in CT images of a rat, resulting from different irradiation configurations used in several experiments. It reduced computational time by an order of magnitude. This allows us to speed up simulations for in vivo trials, where we obtained peak-to-valley dose ratios of 1.20 ± 0.05 and 6.1 ± 0.2 for proton minibeam irradiations targeting the tumour and crossing the rat head. Tumour eradication was observed in the 67 and 22% of the animals treated respectively. CONCLUSION: A Monte Carlo dose calculation engine for pMBRT implementation with mechanical collimation has been developed. This tool can be used to guide and interpret the results of in vivo trials. ADVANCES IN KNOWLEDGE: This is the first Monte Carlo dose engine for pMBRT that is being used to guide preclinical trials in a clinical proton therapy centre.

A new Monte Carlo model of a Cyberknife® system for the precise determination of out-of-field doses.

In a previous work, a PENELOPE Monte Carlo model of a Cyberknife system equipped with fixed collimator was developed and validated for in-field dose evaluation. The aim of this work is to extend it to evaluate peripheral doses and to determine the precision of the treatment planning system (TPS) Multiplan in evaluating the off-axis doses. The Cyberknife® head model was completed with surrounding components based on manufacturer drawings. The contribution of the different head parts on the out-of-field dose was studied. To model the attenuation and the modification of particle energy caused by components not modelled, the photon transport was modified in one of the added components. The model was iteratively adjusted to fit dose profiles measured with EBT3 films and an ionization chamber for several collimator sizes. Finally, dose profiles were calculated using the two Multiplan TPS algorithms and were compared to our simulations. The contributions to out-of-field dose were identified as scattered radiation from the phantom and head leakage and scatter originating at the secondary collimator level. Particle transport in the additional pieces was modified to model this radiation. The maximum differences between simulated and measured doses are of 20.4%. Regarding the detector responses away from axis, EBT3 films and the Farmer chamber give similar response (less than 20% difference). The TPS Monte Carlo algorithm underestimates the doses away from axis more importantly for the smaller field sizes (up to 98%). Besides, RayTracing simplifies peripheral dose to a constant value with no inclusion of particle transport. A Monte Carlo model of a Cyberknife system for the determination of out-of-field doses up to 14 cm off-axis was successfully developed and validated for different depths and field sizes in comparison with measurements. This study also confirms that TPS algorithms do not model peripheral dose properly.

Extended field radiotherapy measurements in a single shot using a BaFBr-based OSL-film.

This work evaluated the use of a class solution specific calibration for an extra-large BaFBr-based optically stimulated luminescence film (OSL; 43 × 35 cm2; Z eff = 4.55). The clinical need for such large dosimeters follows from the increased use of extended-field radiation therapy (EFRT). E.g. for prostate cancer EFRT is currently used in the first prospective trial investigating the benefit of adding elective irradiation of the para-aortic lymph nodes in pN1 prostate cancer. The full extent of these EFRT dose distributions is not covered by the well-established standard sized radiochromic film or 2D detector arrays. Here we investigate an OSL calibration methodology, that tackles BaFBr-based OSL's inherent energy dependence by a class solution specific calibration. 10 EFRT treatment plans used in the PART trial were investigated. One plan was used to build a class solution specific bilinear calibration model, that distinguishes between in-field and penumbra dose contributions. The effect of this calibration was evaluated with respect to a standard linear calibration, using standard IMRT patterns, the nine remaining patient plans, and to smaller prostate treatment plans. A single OSL-dosimeter could be reused for all measurements. The dosimeter captured the full extent of the dose distributions (maximum EFRT field size = 33.5 cm). The bilinear correction reduced the residual dose differences from above 10% to an average of 0.7% (max 3.6%) in comparison with a Monte Carlo simulation. Consequently global gamma agreement scores (3%-3 mm) of 95.5% ± 2.7% were reached. A more strict local evaluation resulted in an average gamma-agreement score of 93.3% ± 3.2%. The BaFBr-based OSL film, with reduced Z eff requires a class-solution specific correction. The current work shows that such a correction can be as simple as a bilinear residual dose correction driven by the measured signal. As far as we know this is the first 2D dosimeter combining reusability, a sub-mm resolution, and a size covering the typical EFRT treatment plans.

Radiological tissue equivalence of deformable silicone-based chemical radiation dosimeters (FlexyDos3D).

FlexyDos3D, a silicone-based chemical radiation dosimeter, has great potential to serve as a three-dimensional (3D) deformable dosimetric tool to verify complex dose distributions delivered by modern radiotherapy techniques. To facilitate its clinical application, its radiological tissue needs to be clarified. In this study we investigated its tissue-equivalence in comparison with water and Solid Water (RMI457). We found that its effective and mean atomic numbers were 40% and 20% higher and the total interaction probabilities for kV x-ray photons were larger than those of water respectively. To assess the influence of its over-response to kV photons, its HU value was measured by kV computed tomography (CT) and was found higher than all the soft-tissue substitutes. When applied for dose calculation without correction, this effect led to an 8% overestimation in electron density via HU-value mapping and 0.65% underestimation in target dose. Furthermore, depth dose curves (PDDs) and off-axis ratios (profiles) at various beam conditions as well as the dose distribution of a full-arc VMAT plan in FlexyDos3D and reference materials were simulated by Monte Carlo, where the results showed great agreement. As indicated, FlexyDos3D exhibits excellent radiological water-equivalence for clinical MV x-ray dosimetry, while its nonwater-equivalent effect for low energy x-ray dosimetry requires necessary correction. The key findings of this study provide pertinent reference for further FlexyDos3D characterization research.

Investigation of TLD and EBT3 performance under the presence of 1.5T, 0.35T, and 0T magnetic field strengths in MR/CT visible materials.

PURPOSE: The aim of this study was to investigate thermoluminescent dosimeters (TLD) and radiochromic EBT3 film inside MR/CT visible geometric head and thorax phantoms in the presence of: 0, 0.35, and 1.5 T magnetic fields. METHODS: Thermoluminescent Dosimeters reproducibility studies were examined by irradiating IROC-Houston's TLD acrylic block five times under 0 and 1.5 T configurations of Elekta's Unity system and three times under 0 and 0.35 T configurations of ViewRay's MRIdian Cobalt-60 (60 Co) system. Both systems were irradiated with an equivalent 10 × 10 cm2 field size, and a prescribed dose of 3 Gy to the maximum depth deposition (dmax). EBT3 film and TLDs were investigated using two geometrical Magnetic Resonance (MR)-guided Radiation Therapy (MRgRT) head and thorax phantoms. Each geometrical phantom had eight quadrants that combined to create a centrally located rectangular tumor (3 × 3 × 5 cm3 ) surrounded by tissue to form a 15 × 15 × 15 cm3 cubic phantom. Liquid polyvinyl chloride plastic and Superflab were used to simulate the tumor and surrounding tissue in the head phantom, respectively. Synthetic ballistic gel and a heterogeneous in-house mixture were used to construct the tumor and surrounding tissue in the thorax phantom, respectively. EBT3 and double-loaded TLDs were used in the phantoms to compare beam profiles and point dose measurements with and without magnetic fields. GEANT4 Monte Carlo simulations were performed to validate the detectors for both Unity 0 T/1.5 T and MRIdian 0 T/0.35 T configurations. RESULTS: Average TLD block measurements which, compared the magnetic field effects (magnetic field vs 0 T) on the Unity and MRIdian systems, were 0.5% and 0.6%, respectively. The average ratios between magnetic field effects for the geometric thorax and head phantoms under the Unity system were -0.2% and 1.6% and for the MRIdian system were 0.2% and -0.3%, respectively. Beam profiles generated with both systems agreed with Monte Carlo measurements and previous literature findings. CONCLUSIONS: TLDs and EBT3 film dosimeters could potentially be used in MR/CT visible tissue equivalent phantoms that will experience a magnetic field environment.

A linear relationship for the LET-dependence of Gafchromic EBT3 film in spot-scanning proton therapy.

Radiochromic film (RCF) is a valuable dosimetric tool, primarily due to its sub-millimeter spatial resolution. For accurate proton dosimetry, the dependence of film response on linear energy transfer (LET) must be characterized and calibrated. In this work, we characterized film under-response, or 'quenching', as a function of dose-weighted linear energy transfer (LETd) in several proton fields and established a simple, linear relationship with LETd. We performed measurements as a function of depth in a PMMA phantom irradiated by a spot-scanning proton beam. The fields had energies of 71.3 MeV, 71.3 MeV with filter, and 159.9 MeV. At each depth (measurements taken in depth step sizes of 0.5-1 mm in the Bragg peak), we measured dose with a PTW Markus ion chamber and EBT3 RCF. EBT3 under-response was characterized by the ratio of dose measured with film to that with ion chamber. LETd values for our experimental setup were calculated using in-house clinical Monte Carlo code. Measured film under-response increased with LETd, from approximately 10% under-response for LETd = 5 keV µm-1 to approximately 20% for LETd = 8 keV µm-1. The under-response for all measurements was plotted versus LETd. A linear fit to the data was performed, yielding a function for under-response, [Formula: see text], with respect to LETd: [Formula: see text]. Finally, the linear under-response relationship was applied to a film measurement within a spread-out Bragg peak (SOBP). Without correcting for LETd-dependence in the SOBP measurement, the discrepancy between film and Monte Carlo profiles was greater than 15% at the distal edge. With correction, the corrected film profile was within 2% and 1 mm of the Monte Carlo profile. RCF response depends on LETd, potentially under-responding by >15% in clinically-relevant scenarios. Therefore, it is insufficient to perform only a dose calibration; LET calibration is also necessary for accurate proton film dosimetry. The LETd-dependence of EBT3 can be described by a single, linear function over a range of clinically-relevant proton therapy LETd values.

Technical Note: Manufacturing of a realistic mouse phantom for dosimetry of radiobiology experiments.

PURPOSE: The goal of this work was to design a realistic mouse phantom as a useful tool for accurate dosimetry in radiobiology experiments. METHODS: A subcutaneous tumor-bearing mouse was scanned in a microCT scanner, its organs manually segmented and contoured. The resulting geometries were converted into a stereolithographic file format (STL) and sent to a multimaterial 3D printer. The phantom was split into two parts to allow for lung excavation and 3D-printed with an acrylic-like material and consisted of the main body (mass density ρ=1.18 g/cm3 ) and bone (ρ=1.20 g/cm3 ). The excavated lungs were filled with polystyrene (ρ=0.32 g/cm3 ). Three cavities were excavated to allow the placement of a 1-mm diameter plastic scintillator dosimeter (PSD) in the brain, the center of the body and a subcutaneous tumor. Additionally, a laser-cut Gafchromic film can be placed in between the two phantom parts for 2D dosimetric evaluation. The expected differences in dose deposition between mouse tissues and the mouse phantom for a 220-kVp beam delivered by the small animal radiation research platform (SARRP) were calculated by Monte Carlo (MC). RESULTS: MicroCT scans of the phantom showed excellent material uniformity and confirmed the material densities given by the manufacturer. MC dose calculations revealed that the dose measured by tissue-equivalent dosimeters inserted into the phantom in the brain, abdomen, and subcutaneous tumor would be underestimated by 3-5%, which is deemed to be an acceptable error assuming the proposed 5% accuracy of radiobiological experiments. CONCLUSIONS: The low-cost mouse phantom can be easily manufactured and, after a careful dosimetric characterization, may serve as a useful tool for dose verification in a range of radiobiology experiments.

In-phantom neutron dose measurement using Gafchromic film dosimeter for QA of BNCT beams.

In order to improve the quality assurance (QA) procedure for beams of boron neutron capture therapy (BNCT), this study introduced using the Gafchromic film dosimeter for neutron dose measurement of BNCT beams. The crucial part of this study was investigating an approach to employ the Gafchromic film dosimeter placed inside a PMMA phantom irradiated by a BNCT beam. The spatial distribution of neutron dose of the film was determined using measurements and Monte Carlo calculations. By employing the present approach, the two-dimensional distributions of the neutron dose component of the film at specific depths in the phantom were successfully obtained. The determined neutron dose profiles were in good agreement with the calculated ones. This study also confirmed the finding that the film dosimeter is sensitive to thermal neutrons by comparing the depth-capture-reaction-rate and depth-dose distributions. Results of this work not only proved the feasibility of using the proposed method for the QA measurement of beam delivery but also revealed the advantages of easy-handling and remarkable spatial resolution of the film dosimeter when applied to BNCT fields. The present work can help to verify the dose uniformity and output stability of BNCT beams prior to clinical irradiation.

Technical Note: Reconstruction of physical and biological dose distributions of carbon-ion beam through deconvolution of longitudinal dosimeter responses.

PURPOSE: This is a theoretical simulation study for proof of concept of radiochromic film dosimetry to measure physical and biological doses without plan-based quenching correction for patient-specific quality assurance of carbon-ion radiotherapy. METHODS: We took a layer-stacking carbon-ion beam comprised of range-shifted beamlets. The dosimeter response was simulated according to an experimental quenching model. The beam model followed a treatment planning system. The beam was decomposed into finely arranged beamlets with weights estimated by deconvolution of longitudinal dosimeter responses. The distributions of physical and biological doses were reconstructed from the estimated weights and were compared with the plan. We also evaluated the sensitivity to measurement errors and to erratic delivery with an undelivered beamlet. RESULTS: The reconstructed physical and biological doses accurately reproduced the simulated delivery with errors approximately corresponding to the measurement errors. The erratic beam delivery was easily detectable by comparison of biological dose distribution to the plan. CONCLUSIONS: We have developed a method to measure physical and biological doses by longitudinal dosimetry of quenched response without using plan data. The method only involves a general optimization algorithm, a radiobiology model, and experimental beamlet data, and requires no extra corrections. Theoretically, this approach is applicable to various dosimeters and to proton and ion beams of any delivery method, regardless of quenching or biological effectiveness.

A novel rectal applicator for contact radiotherapy with HDR 192Ir sources.

PURPOSE: Dose escalation to rectal tumors leads to higher complete response rates and may thereby enable omission of surgery. Important advantages of endoluminal boosting techniques include the possibility to apply a more selective/localized boost than using external beam radiotherapy. A novel brachytherapy (BT) rectal applicator with lateral shielding was designed to be used with a rectoscope for eye-guided positioning to deliver a dose distribution similar to the one of contact x-ray radiotherapy devices, using commonly available high-dose-rate 192Ir BT sources. METHODS AND MATERIALS: A cylindrical multichannel BT applicator with lateral shielding was designed by Monte Carlo modeling, validated experimentally with film dosimetry and compared with results found in the literature for the Papillon 50 (P50) contact x-ray radiotherapy device regarding rectoscope dimensions, radiation beam shape, dose fall-off, and treatment time. RESULTS: The multichannel applicator designed is able to deliver 30 Gy under 13 min with a 20350 U (5 Ci) source. The use of multiple channels and lateral shielding provide a uniform circular treatment surface with 22 mm in diameter. The resulting dose fall-off is slightly steeper (maximum difference of 5%) than the one generated by the P50 device with the 22 mm applicator. CONCLUSIONS: A novel multichannel rectal applicator for contact radiotherapy with high-dose-rate 192Ir sources that can be integrated with commercially available treatment planning systems was designed to produce a dose distribution similar to the one obtained by the P50 device.

Monte Carlo optimization of a microbeam collimator design for use on the small animal radiation research platform (SARRP).

Microbeam radiation therapy (MRT) is a pre-clinical, spatially-fractionated treatment modality noted for its ability to achieve a large differential response between normal and tumoral tissues. In the present study, TOPAS Monte Carlo (MC) simulations were used to optimize the design of a compact, affordable multi-slit collimator (MSC) suitable for use with the small animal radiation research platform (SARRP). MRT dose distributions in a (1 × 1 × 3)cm3 water phantom were simulated for a tungsten MSC using different focal spot sizes (0.4, 3 mm), beam energies (40, 80, 220 kVp), slit widths (100, 125, 150, 175, 200 µm), collimator thicknesses (1.5, 2.5, 3 cm) and collimator-to-surface distances (CSD of 1 and 3 cm). Key MRT figures of merit, namely the peak-to-valley dose ratio (PVDR), full-width at half-maximum and peak dose rate were determined. Use of the small focal spot maximized the PVDR (~40 at surface) and reduced the system's sensitivity to changes in CSD, but decreased the collimated beam output to 55.2 cGy min-1. The large focal spot was ill-suited for large CSD irradiations, but increased the beam output by a factor of 2.8, to 153.0 cGy min-1, and decreased the sensitivity to changes in slit width. A modular MSC, using divergent plastic spacer materials in place of excavated slits, was also investigated. Polypropylene and polyethylene terephthalate material spacers were considered and while neither reduced the PVDR compared to air slits, the dose rate was reduced by 37% and 47%, respectively. Lastly, a steel parallel-slit MSC was used in a preliminary test of MRT delivery using the SARRP. Discrepancies between the results of film dosimetry and the corresponding MC simulations highlight the need to fabricate a more well-defined collimator for use in future validation and radiobiological work. The simulated results of this study are being used to inform the design of such a collimator, which will additionally boast a high degree of modularity at reasonable cost.

Dosimetric verification and commissioning for a small animal image-guided irradiator.

In recent years, small animal image-guided irradiators have been widely utilized in preclinical studies involving rodent models. However, the dosimetry commissioning of such equipment involving kilovoltage small-field dosimetry has not received as much interest as the megavoltage small-field dosimetry used clinically. To date, a paucity of measured kilovoltage beam data, especially for field sizes less than 3 mm, can be found in the literature. For improvement of rodent treatments in the future, this work aims to provide comprehensive and accurate beam data for the small animal radiation research platform (SARRP, Xstrahl) using EBT3 Gafchromic films and Monte Carlo calculation, with submillimeter resolution and accuracy. This work includes three primary tasks: (1) establish an optimized film measurement protocol for small field dosimetry of kilovoltage photon beam. (2) Acquire dosimetric data including (a) depth dose curves from the surface to 6 cm depth (b) beam profiles, (c) penumbra, (d) cone factors and (e) 2D dose distribution. These tasks were undertaken for a 220 kVp photon beam with five different small field widths and 33 cm source to surface distance (0.5 mm and 1 mm circular fields, 3 × 3 mm2, 5 × 5 mm2, 10 × 10 mm2 square fields). Beam data was measured with EBT3 films. (3) Provide comparative dosimetry for film measurements, Monte Carlo calculations, and the dose calculations performed with the SARRP treatment planning system, Muriplan. For the majority of parameters, film measurement agreed with Monte Carlo simulation within 1%. There were, however, discrepancies between measured beam data and Muriplan treatment planning data. Specifically, for PDD, Muriplan underestimates the dose for field sizes of 0.5 mm and 1 mm. For beam profiles comparisons, the calculation from Muriplan predicts a smaller lateral distance between the 50% isodose lines compared to film measurement. There is a difference of 0.18, 0.72, 0.6 mm between Muriplan and film for field sizes of 3, 5, 10 mm, respectively. This work demonstrates that accurate and precise kilovoltage small-field dosimetry can be conducted using EBT3 Gafchromic film with an optimized protocol. In addition, discrepancies between measured beam data and Muriplan were identified.

Monte Carlo uncertainty analysis of dose estimates in radiochromic film dosimetry with single-channel and multichannel algorithms.

PURPOSE: To provide a multi-stage model to calculate uncertainty in radiochromic film dosimetry with Monte-Carlo techniques. This new approach is applied to single-channel and multichannel algorithms. MATERIAL AND METHODS: Two lots of Gafchromic EBT3 are exposed in two different Varian linacs. They are read with an EPSON V800 flatbed scanner. The Monte-Carlo techniques in uncertainty analysis provide a numerical representation of the probability density functions of the output magnitudes. From this numerical representation, traditional parameters of uncertainty analysis as the standard deviations and bias are calculated. Moreover, these numerical representations are used to investigate the shape of the probability density functions of the output magnitudes. Also, another calibration film is read in four EPSON scanners (two V800 and two 10000XL) and the uncertainty analysis is carried out with the four images. RESULTS: The dose estimates of single-channel and multichannel algorithms show a Gaussian behavior and low bias. The multichannel algorithms lead to less uncertainty in the final dose estimates when the EPSON V800 is employed as reading device. In the case of the EPSON 10000XL, the single-channel algorithms provide less uncertainty in the dose estimates for doses higher than four Gy. CONCLUSION: A multi-stage model has been presented. With the aid of this model and the use of the Monte-Carlo techniques, the uncertainty of dose estimates for single-channel and multichannel algorithms are estimated. The application of the model together with Monte-Carlo techniques leads to a complete characterization of the uncertainties in radiochromic film dosimetry.

Validation of GPU-accelerated superposition-convolution dose computations for the Small Animal Radiation Research Platform.

PURPOSE: The Small Animal Radiation Research Platform (SARRP) has been developed for conformal microirradiation with on-board cone beam CT (CBCT) guidance. The graphics processing unit (GPU)-accelerated Superposition-Convolution (SC) method for dose computation has been integrated into the treatment planning system (TPS) for SARRP. This paper describes the validation of the SC method for the kilovoltage energy by comparing with EBT2 film measurements and Monte Carlo (MC) simulations. METHODS: MC data were simulated by EGSnrc code with 3 × 108 -1.5 × 109 histories, while 21 photon energy bins were used to model the 220 kVp x-rays in the SC method. Various types of phantoms including plastic water, cork, graphite, and aluminum were used to encompass the range of densities of mouse organs. For the comparison, percentage depth dose (PDD) of SC, MC, and film measurements were analyzed. Cross beam (x,y) dosimetric profiles of SC and film measurements are also presented. Correction factors (CFz) to convert SC to MC dose-to-medium are derived from the SC and MC simulations in homogeneous phantoms of aluminum and graphite to improve the estimation. RESULTS: The SC method produces dose values that are within 5% of film measurements and MC simulations in the flat regions of the profile. The dose is less accurate at the edges, due to factors such as geometric uncertainties of film placement and difference in dose calculation grids. CONCLUSION: The GPU-accelerated Superposition-Convolution dose computation method was successfully validated with EBT2 film measurements and MC calculations. The SC method offers much faster computation speed than MC and provides calculations of both dose-to-water in medium and dose-to-medium in medium.

Magnetic fields are causing small, but significant changes of the radiochromic EBT3 film response to 6 MV photons.

The optical density (OD) of EBT3 radiochromic films (Ashland Specialty Ingredients, Bridgewater, NJ, USA) exposed to absorbed doses to water up to D = 20 Gy in magnetic fields of B = 0.35 and 1.42 T was measured in the three colour channels of an Epson Expression 10000XL flatbed scanner. A 7 cm wide water phantom with fixed film holder was placed between the pole shoes of a constant-current electromagnet with variable field strength and was irradiated by a 6 MV photon beam whose axis was directed at right angles with the field lines. The doses at the film position at water depth 5 cm were measured with a calibrated ionization chamber when the magnet was switched off and were converted to the doses in presence of the magnetic field via the monitor units and by a Monte Carlo-calculated correction accounting for the slight change of the depth dose curves in magnetic fields. In the presence of the 0.35 and 1.42 T fields small negative changes of the OD values at given absorbed doses to water occurred and just significantly exceeded the uncertainty margin given by the stochastic and the uncorrected systematic deviations. This change can be described by a +2.1% change of the dose values needed to produce a given optical density in the presence of a 1.42 T field. The thereby modified OD versus D function remained unchanged irrespective of whether the original short film side-the preference direction of the monomer crystals of the film-was directed parallel or orthogonal to the magnetic field. The 'orientation effect', the difference between the optical densities measured in the 'portrait' or 'landscape' film positions on the scanner bed caused by the reflection of polarised light in the scanner's mirror system, remained unaltered after EBT3 film exposure in magnetic fields. An independent optical bench investigation of EBT3 films exposed to doses of 10 and 20 Gy at 0.35 and 1.42 T showed that the direction of the electric vector of polarised light experiencing the largest transmission through EBT3 films remained unaltered after film exposure in the magnetic fields. The observed small modification of the OD versus D curve of the radiochromic film EBT3 in the range up to 20 Gy and 1.42 T, hardly exceeding the experimental uncertainty margin, numerically confirms other recent studies on EBT3 film. A stronger magnetic field effect had been observed with the previous product EBT2 exposed to 60Co gamma radiation at 0.35 T.