X-ray-Sensitive Doped CaF2-Based MRI Contrast Agents for Local Radiation Dose Measurement

Ionizing radiation has become widely used in medicine, with application in diagnostic techniques, such as computed tomography (CT) and radiation therapy (RT), where X-rays are used to diagnose and treat tumors. The X-rays used in CT and, in particular, in RT can have harmful side effects; hence, an accurate determination of the delivered radiation dose is of utmost importance to minimize any damage to healthy tissues. For this, medical specialists mostly rely on theoretical predictions of the delivered dose or external measurements of the dose. To extend the practical use of ionizing radiation-based medical techniques, such as magnetic resonance imaging (MRI)-guided RT, a more precise measurement of the internal radiation dose internally is required. In this work, a novel approach is presented to measure dose in liquids for potential future in vivo applications. The strategy relies on MRI contrast agents (CAs) that provide a dose-sensitive signal. The demonstrated materials are (citrate-capped) CaF2 nanoparticles (NPs) doped with Eu3+ or Fe2+/Fe3+ ions. Free electrons generated by ionizing radiation allow the reduction of Eu3+, which produces a very small contrast in MRI, to Eu2+, which induces a strong contrast. Oxidative species generated by high-energy X-rays can be measured indirectly using Fe2+ because it oxidizes to Fe3+, increasing the contrast in MRI. Notably, in the results, a strong increase in the proton relaxation rates is observed for the Eu3+-doped NPs at 40 kV. At 6 MV, a significant increase in proton relaxation rates is observed using CaF2 NPs doped with Fe2+/Fe3+ after irradiation. The presented concept shows great promise for use in the clinic to measure in vivo local ionizing radiation dose, as these CAs can be intravenously injected in a saline solution.

Material decomposition with a prototype photon-counting detector CT system: expanding a stoichiometric dual-energy CT method via energy bin optimization and K-edge imaging

Objective.Computed tomography (CT) has advanced since its inception, with breakthroughs such as dual-energy CT (DECT), which extracts additional information by acquiring two sets of data at different energies. As high-flux photon-counting detectors (PCDs) become available, PCD-CT is also becoming a reality. PCD-CT can acquire multi-energy data sets in a single scan by spectrally binning the incident x-ray beam. With this, K-edge imaging becomes possible, allowing high atomic number (high-Z) contrast materials to be distinguished and quantified. In this study, we demonstrated that DECT methods can be converted to PCD-CT systems by extending the method of Bourqueet al(2014). We optimized the energy bins of the PCD for this purpose and expanded the capabilities by employing K-edge subtraction imaging to separate a high-atomic number contrast material.Approach.The method decomposes materials into their effective atomic number (Zeff) and electron density relative to water (ρe). The model was calibrated and evaluated using tissue-equivalent materials from the RMI Gammex electron density phantom with knownρevalues and elemental compositions. TheoreticalZeffvalues were found for the appropriate energy ranges using the elemental composition of the materials.Zeffvaried slightly with energy but was considered a systematic error. Anex vivobovine tissue sample was decomposed to evaluate the model further and was injected with gold chloride to demonstrate the separation of a K-edge contrast agent.Main results.The mean root mean squared percent errors on the extractedZeffandρefor PCD-CT were 0.76% and 0.72%, respectively and 1.77% and 1.98% for DECT. The tissue types in theex vivobovine tissue sample were also correctly identified after decomposition. Additionally, gold chloride was separated from theex vivotissue sample with K-edge imaging.Significance.PCD-CT offers the ability to employ DECT material decomposition methods, along with providing additional capabilities such as K-edge imaging.

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.

Dosimetric characterization of a novel UHDR megavoltage X-ray source for FLASH radiobiological experiments

A first irradiation platform capable of delivering 10 MV X-ray beams at ultra-high dose rates (UHDR) has been developed and characterized for FLASH radiobiological research at TRIUMF. Delivery of both UHDR (FLASH mode) and low dose-rate conventional (CONV mode) irradiations was demonstrated using a common source and experimental setup. Dose rates were calculated using film dosimetry and a non-intercepting beam monitoring device; mean values for a 100 μA pulse (peak) current were nominally 82.6 and 4.40 × 10-2 Gy/s for UHDR and CONV modes, respectively. The field size for which > 40 Gy/s could be achieved exceeded 1 cm down to a depth of 4.1 cm, suitable for total lung irradiations in mouse models. The calculated delivery metrics were used to inform subsequent pre-clinical treatments. Four groups of 6 healthy male C57Bl/6J mice were treated using thoracic irradiations to target doses of either 15 or 30 Gy using both FLASH and CONV modes. Administration of UHDR X-ray irradiation to healthy mouse models was demonstrated for the first time at the clinically-relevant beam energy of 10 MV.

Non-coplanar lung SABR treatments delivered with a gantry-mounted x-ray tube

Objective.To create two non-coplanar, stereotactic ablative radiotherapy (SABR) lung patient treatment plans compliant with the radiation therapy oncology group (RTOG) 0813 dosimetric criteria using a simple, isocentric, therapy with kilovoltage arcs (SITKA) system designed to provide low cost external radiotherapy treatments for low- and middle-income countries (LMICs).Approach.A treatment machine design has been proposed featuring a 320 kVp x-ray tube mounted on a gantry. A deep learning cone-beam CT (CBCT) to synthetic CT (sCT) method was employed to remove the additional cost of planning CTs. A novel inverse treatment planning approach using GPU backprojection was used to create a highly non-coplanar treatment plan with circular beam shapes generated by an iris collimator. Treatments were planned and simulated using the TOPAS Monte Carlo (MC) code for two lung patients. Dose distributions were compared to 6 MV volumetric modulated arc therapy (VMAT) planned in Eclipse on the same cases for a Truebeam linac as well as obeying the RTOG 0813 protocols for lung SABR treatments with a prescribed dose of 50 Gy.Main results.The low-cost SITKA treatments were compliant with all RTOG 0813 dosimetric criteria. SITKA treatments showed, on average, a 6.7 and 4.9 Gy reduction of the maximum dose in soft tissue organs at risk (OARs) as compared to VMAT, for the two patients respectively. This was accompanied by a small increase in the mean dose of 0.17 and 0.30 Gy in soft tissue OARs.Significance.The proposed SITKA system offers a maximally low-cost, effective alternative to conventional radiotherapy systems for lung cancer patients, particularly in low-income countries. The system's non-coplanar, isocentric approach, coupled with the deep learning CBCT to sCT and GPU backprojection-based inverse treatment planning, offers lower maximum doses in OARs and comparable conformity to VMAT plans at a fraction of the cost of conventional radiotherapy.

Characterization of a 0.8 mm3Medscint plastic scintillator detector system for small field dosimetry

Objective. Plastic scintillator detectors (PSDs) have demonstrated ability to meet requirements of small field dosimetry. Medscint developed a 1 mm long, 1 mm diameter cylindrical PSD with effective volume of 0.8 mm3. Clinically relevant, small field dosimetric properties of this detector, combined with a novel scintillation dosimetry system-HYPERSCINT RP-200, and HYPERDOSE analysis software were evaluated in this study.Approach. This novel scintillator-based dosimetry system was characterized with 6 MV-WFF and 10 MV-FFF x-ray beams delivered by Varian TrueBeamTMlinear accelerator. The detector was characterized for leakage, short-term repeatability, dose response linearity, angular response, dose rate response, and field size dependence for radiation field sizes of 0.25 × 0.25 to 10 × 10 cm2. Measured detector specific output ratios were compared with microDiamond output factors to determine small field output correction factors,kQclin,Qmsrfclin,fmsr.Main results. The dosimetry system showed excellent short-term repeatability with standard deviation of only 0.04 ± 0.01%. It demonstrated good dose linearity with variations less than 1.0% for 14.4 cGy and above. The dosimetry system was found to be independent of dose rate and angle of irradiation, with deviations for both below 0.5%. Leakage was found to be comparable to background readings. For 6 MV-WFF energy beams, detector specific output ratios for field sizes down to 1 × 1 cm2agreed with output factors measured with PTW TN60019 microDiamond, thus,kQclin,Qmsrfclin,fmsrequates to unity for these field sizes. For 10 MV-FFF energy beams, detector specific output ratios for field sizes down to 2 × 2 cm2agreed with PTW TN60019 microDiamond output factors, thus,kQclin,Qmsrfclin,fmsrequates to unity for these field sizes.kQclin,Qmsrfclin,fmsrfor field sizes down to 0.5 × 0.5 cm2were determined to be within 6% of unity for both 6 MV-WFF and 10 MV-FFF energy beams.Significance. The HYPERSCINT RP-200 dosimetry system coupled with a 0.8 mm3PSD showed excellent dosimetric properties and was found to be clinically relevant for relative dosimetry down to field sizes of 0.5 × 0.5 cm2and potentially smaller.

A feasibility study of ultra-high dose rate mini-GRID therapy using very-high-energy electron beams for a simulated pediatric brain case

Ultra-high dose rate (UHDR, >40 Gy/s), spatially-fractionated minibeam GRID (mini-GRID) therapy using very-high-energy electrons (VHEE) was investigated using Monte Carlo simulations. Multi-directional VHEE treatments with and without mini-GRID-fractionation were compared to a clinical 6 MV volumetric modulated arc therapy (VMAT) plan for a pediatric glioblastoma patient using dose-volume histograms, volume-averaged dose rates in critical patient structures, and planning target volume D98s. Peak-to-valley dose ratios (PVDRs) and dose rates in organs at risk (OARs) were evaluated due to their relevance for normal-tissue sparing in FLASH and spatially-fractionated techniques. Depths of convergence, defined where the PVDR is first ≤1.1, and depths at which dose rates fall below the UHDR threshold were also evaluated. In a water phantom, the VHEE mini-GRID treatments presented a surface (5 mm depth) PVDR of (51±2) and a depth of convergence of 42 mm at 150 MeV and a surface PVDR of (33±1) with a depth of convergence of 57 mm at 250 MeV. For a pediatric GBM case, VHEE treatments without mini-GRID-fractionation produced 25% and 22% lower volume-averaged doses to OARs compared to the 6 MV VMAT plan and 8/9 and 9/9 of the patient structures were exposed to volume-averaged dose rates >40 Gy/s for the 150 MeV and 250 MeV plans, respectively. The 150 MeV and 250 MeV mini-GRID treatments produced 17% and 38% higher volume-averaged doses to OARs and 3/9 patient structures had volume-averaged dose rates above 40 Gy/s. VHEE mini-GRID plans produced many comparable dose metrics to the clinical VMAT plan, encouraging further optimization.

Framework for Quality Assurance of Ultrahigh Dose Rate Clinical Trials Investigating FLASH Effects and Current Technology Gaps

FLASH radiation therapy (FLASH-RT), delivered with ultrahigh dose rate (UHDR), may allow patients to be treated with less normal tissue toxicity for a given tumor dose compared with currently used conventional dose rate. Clinical trials are being carried out and are needed to test whether this improved therapeutic ratio can be achieved clinically. During the clinical trials, quality assurance and credentialing of equipment and participating sites, particularly pertaining to UHDR-specific aspects, will be crucial for the validity of the outcomes of such trials. This report represents an initial framework proposed by the NRG Oncology Center for Innovation in Radiation Oncology FLASH working group on quality assurance of potential UHDR clinical trials and reviews current technology gaps to overcome. An important but separate consideration is the appropriate design of trials to most effectively answer clinical and scientific questions about FLASH. This paper begins with an overview of UHDR RT delivery methods. UHDR beam delivery parameters are then covered, with a focus on electron and proton modalities. The definition and control of safe UHDR beam delivery and current and needed dosimetry technologies are reviewed and discussed. System and site credentialing for large, multi-institution trials are reviewed. Quality assurance is then discussed, and new requirements are presented for treatment system standard analysis, patient positioning, and treatment planning. The tables and figures in this paper are meant to serve as reference points as we move toward FLASH-RT clinical trial performance. Some major questions regarding FLASH-RT are discussed, and next steps in this field are proposed. FLASH-RT has potential but is associated with significant risks and complexities. We need to redefine optimization to focus not only on the dose but also on the dose rate in a manner that is robust and understandable and that can be prescribed, validated, and confirmed in real time. Robust patient safety systems and access to treatment data will be critical as FLASH-RT moves into the clinical trials.

Publicly available framework for simulating and experimentally validating clinical PET systems

BACKGROUND

Monte Carlo (MC) simulations are a powerful tool to model medical imaging systems. However, before simulations can be considered the ground truth, they have to be validated with experiments.

PURPOSE

To provide a pipeline that models a clinical positron emission tomography (PET)/CT system using MC simulations after extensively validating the results against experimental measurements.

METHODS

A clinical four-ring PET imaging system was modeled using Geant4 application for tomographic emission (v. 9.0). To validate the simulations, PET images were acquired of a cylindrical phantom, point source, and image quality phantom with the modeled system and the simulations of the experimental procedures. For the purpose of validating the quantification capabilities and image quality provided by the simulation pipeline, the simulations were compared against the measurements in terms of their count rates and sensitivity as well as their image uniformity, resolution, recovery coefficients (RCs), coefficients of variation, contrast, and background variability.

RESULTS

When compared to the measured data, the number of true detections in the MC simulations was within 5%. The scatter fraction was found to be 30.0% ± 2.2% and 28.8% ± 1.7% in the measured and simulated scans, respectively. Analyzing the measured and simulated sinograms, the sensitivities were found to be 8.2 and 7.8 cps/kBq, respectively. The fraction of random coincidences were 19% in the measured data and 25% in the simulation. When calculating the image uniformity within the axial slices, the measured image exhibited a uniformity of 0.015 ± 0.005, whereas the simulated image had a uniformity of 0.029 ± 0.011. In the axial direction, the uniformity was measured to be 0.024 ± 0.006 and 0.040 ± 0.015 for the measured and simulated data, respectively. Comparing the image resolution, an average percentage difference of 2.9% was found between the measurements and simulations. The RCs calculated in both the measured and simulated images were found to be within the EARL ranges, except for that of the simulation of the smallest sphere. The coefficients of variation for the measured and simulated images were found to be 12% and 13%, respectively. Lastly, the background variability was consistent between the measurements and simulations, whereas the average percentage difference in the sphere contrasts was found to be 8.8%.

CONCLUSION

The clinical PET/CT system was modeled and validated to provide a simulation pipeline for the community. The pipeline and the validation procedures have been made available (https://github.com/teaghan/PET_MonteCarlo).

Design and CT imaging of casper, an anthropomorphic breathing thorax phantom

The goal of this work was to build an anthropomorphic thorax phantom capable of breathing motion with materials mimicking human tissues in x-ray imaging applications. The thorax phantom, named Casper, was composed of resin (body), foam (lungs), glow polyactic acid (bones) and natural polyactic acid (tumours placed in the lungs). X-ray attenuation properties of all materials prior to manufacturing were evaluated by means of photon-counting computed tomography (CT) imaging on a table-top system. Breathing motion was achieved by a scotch-yoke mechanism with diaphragm motion frequencies of 10-20 rpm and displacements of 1 to 2 cm. Casper was manufactured by means of 3D printing of moulds and ribs and assembled in a complex process. The final phantom was then scanned using a clinical CT scanner to evaluate material CT numbers and the extent of tumour motion. Casper CT numbers were close to human CT numbers for soft tissue (46 HU), ribs (125 HU), lungs (-840 HU) and tumours (-45 HU). For a 2 cm diaphragm displacement the largest tumour displacement was 0.7 cm. The five tumour volumes were accurately assessed in the static CT images with a mean absolute error of 4.3%. Tumour sizes were either underestimated for smaller tumours or overestimated for larger tumours in dynamic CT images due to motion blurring with a mean absolute difference from true volumes of 10.3%. More Casper information including a motion movie and manufacturing data can be downloaded from http://web.uvic.ca/~bazalova/Casper/.

Next generation high resolution perovskite direct conversion detector: Monte Carlo design optimisation and virtual clinical trial

We implement virtual clinical integration of next-generation perovskite detectors into common x-ray imaging devices. This was achieved by performing Monte Carlo (MC) optimisation of the design and benchmarking of low cost, high spatial resolution, direct conversion perovskite crystal x-ray flat panel imagers for a next generation of breast-, MV-, and kV-cone beam CT detectors. Semiconductor methylammonium lead bromide perovskite crystals energy deposition efficiencies calculated in TOPAS were benchmarked against four common detector materials for twelve detector crystal thicknesses between 40 to 15 mm and ten beam energies ranging from 20 keV to 6 MeV. Based on these simulations, Koning's dedicated breast CT, and Varian's Truebeam kV- and MV-cone beam CT systems were designated as suitable applications for perovskite detectors. System specific Fastcat hybrid MC cone beam CT image simulation was subsequently used to optimise the perovskite detector design and conduct virtual clinical trials. Device-specific optimal perovskite crystal thicknesses were calculated to be 0.30, 0.86, and 1.99 mm for Koning breast CT and Truebeam kV- and MV-cone beam CT systems, respectively. Replacing the current detectors on these machines with low cost perovskite crystal detectors could be advantageous as it would simultaneously yield 12.1%, 9.5% and 86.1% improvements in detective quantum efficiency as well as increases in contrast to noise ratio in brain, lung, and bone tissues.

Metal artifact correction in photon-counting detector computed tomography: metal trace replacement using high-energy data

BACKGROUND

Metal artifacts have been an outstanding issue in computed tomography (CT) since its first uses in the clinic and continue to interfere. Metal artifact reduction (MAR) methods continue to be proposed and photon-counting detectors (PCDs) have recently been the subject of research toward this purpose. PCDs offer the ability to distinguish the energy of incident x-rays and sort them in a set number of energy bins. High-energy data captured using PCDs have been shown to reduce metal artifacts in reconstructions due to reduced beam hardening.

PURPOSE

High-energy reconstructions using PCD-CT have their drawbacks, such as reduced image contrast and increased noise. Here, we demonstrate a MAR algorithm, trace replacement MAR (TRMAR), in which the data corrupted by metal artifacts in full energy spectrum projections are corrected using the high-energy data captured during the same scan. The resulting reconstructions offer similar MAR to that seen in high-energy reconstructions, but with improved image quality.

METHODS

Experimental data were collected using a bench-top PCD-CT system with a cadmium zinc telluride PCD. Simulations were performed to determine the optimal high-energy threshold and to test TRMAR in simulations using the XCAT phantom and a biological sample. For experiments a 100-mm diameter cylindrical phantom containing vials of water, two screws, various densities of Ca(ClO₄ )₂ , and a spatial resolution phantom was imaged with and without the screws. The screws were segmented in the initial reconstruction and forward projected to identify them in the sinogram space in order to perform TRMAR. The resulting reconstructions were compared to the control and to reconstructions corrected using normalized metal artifact reduction (NMAR). Additionally, a beef short rib was imaged with and without metal to provide a more realistic phantom.

RESULTS

XCAT simulations showed a reduction in the streak artifact from -978 HU in uncorrected images to -10 HU with TRMAR. The magnitude of the metal artifact in uncorrected images of the 100-mm phantom was -442 HU, compared to the desired -81 HU with no metal. TRMAR reduced the magnitude of the artifact to -142 HU, with NMAR reducing the magnitude to -96 HU. Relative image noise was reduced from 176% in the high-energy image to 56% using TRMAR. Density quantification was better with NMAR, with the Ca(ClO₄ )₂ vial affected most by metal artifacts showing 0.8% error compared to 2.1% with TRMAR. Small features were preserved to a greater extent with TRMAR, with the limiting spatial frequency at 20% of the MTF fully maintained at 1.31 lp/mm, while with NMAR it was reduced to 1.22 lp/mm. Images of the beef short rib showed better delineation of the shape of the metal using TRMAR.

CONCLUSIONS

NMAR offers slightly better performance compared to TRMAR in streak reduction and image quality metrics. However, TRMAR is less susceptible to metal segmentation errors and can closely approximate the reduction in the streak metal artifact seen in NMAR at 1/3 the computation time. With the recent introduction of PCD-CT into the clinic, TRMAR offers notable potential for fast, effective MAR.

Monte Carlo optimization of a GRID collimator for preclinical megavoltage ultra-high dose rate spatially-fractionated radiation therapy

Objective. A 2-dimensional pre-clinical SFRT (GRID) collimator was designed for use on the ultra-high dose rate (UHDR) 10 MV ARIEL beamline at TRIUMF. TOPAS Monte Carlo simulations were used to determine optimal collimator geometry with respect to various dosimetric quantities. Approach. The GRID-averaged peak-to-valley dose ratio (PVDR) and mean dose rate of the peaks were investigated with the intent of maximizing both values in a given design. The effects of collimator thickness, focus position, septal width, and hole width on these metrics were found by testing a range of values for each parameter on a cylindrical GRID collimator. For each tested collimator geometry, photon beams with energies of 10, 5, and 1 MV were transported through the collimator and dose rates were calculated at various depths in a water phantom located 1.0 cm from the collimator exit. Main results. In our optimization, hole width proved to be the only collimator parameter which increased both PVDR and peak dose rates. From the optimization results, it was determined that our optimized design would be one which achieves the maximum dose rate for a PVDR ≥ 5 at 10 MV. Ultimately, this was achieved using a collimator with a thickness of 75 mm, 0.8 mm septal and hole widths, and a focus position matched to the beam divergence. This optimized collimator maintained the PVDR of 5 in the phantom between water depths of 0–10 cm at 10 MV and had a mean peak dose rate of 3.06 ± 0.02 Gy s − 1 at 0–1 cm depth. Significance. We have investigated the impact of various GRID-collimator design parameters on the dose rate and spatial fractionation of 10, 5, and 1 MV photon beams. The optimized collimator design for the 10 MV ultra-high dose rate photon beam could become a useful tool for radiobiology studies synergizing the effects of ultra-high dose rate (FLASH) delivery and spatial fractionation.

Lead-doped scintillator dosimeters for detection of ultrahigh dose-rate x-rays

Objective. Lead-doped scintillator dosimeters may be well suited for the dosimetry of FLASH-capable x-ray radiotherapy beams. Our study explores the dose rate dependence and temporal resolution of scintillators that makes them promising in the accurate detection of ultrahigh dose-rate (UHDR) x-rays. Approach. We investigated the response of scintillators with four material compositions to UHDR x-rays produced by a conventional x-ray tube. Scintillator output was measured using the HYPERSCINT-RP100 dosimetry research platform. Measurements were acquired at high frame rates (400 fps) which allowed for accurate dose measurements of sub-second radiation exposures from 1 to 100 ms. Dose-rate dependence was assessed by scaling tube current of the x-ray tube. Scintillator measurements were validated against Monte Carlo simulations of the probe geometries and UHDR x-ray system. Calibration factors converting dose-to-medium to dose-to-water were obtained from simulation data of plastic and lead-doped scintillator materials. Main Results. The results of this work suggest that lead-doped scintillators were dose-rate independent for UHDR x-rays from 1.1 to 40.1 Gy s−1 and capable of measuring conventional radiotherapy dose-rates (0.1 Gy s−1) at extended distance from the x-ray focal spot. Dose-to-water measured with a 5% lead-doped scintillator detector agreed with simulations within 0.6%. Significance. Lead-doped scintillators may be a valuable tool for the accurate real-time dosimetry of FLASH-capable UHDR x-ray beams.

Design optimization of an electron-to-photon conversion target for ultra-high dose rate x-ray (FLASH) experiments at TRIUMF

Objective. To develop a bremsstrahlung target and megavoltage (MV) x-ray irradiation platform for ultrahigh dose-rate (UHDR) irradiation of small-animals on the Advanced Rare Isotope Laboratory (ARIEL) electron linac (e-linac) at TRIUMF. Approach. An electron-to-photon converter design for UHDR radiotherapy (RT) was centered around optimization of a tantalum–aluminum (Ta–Al) explosion-bonded target. Energy deposition within a homogeneous water-phantom and the target itself were evaluated using EGSnrc and FLUKA MC codes, respectively, for various target thicknesses (0.5–1.5 mm), beam energies (E e− = 8, 10 MeV) and electron (Gaussian) beam sizes ( 2 σ = 2–10 mm). Depth dose-rates in a 3D-printed mouse phantom were also calculated to infer the compatibility of the 10 MV dose distributions for FLASH-RT in small-animal models. Coupled thermo-mechanical FEA simulations in ANSYS were subsequently used to inform the stress–strain conditions and fatigue life of the target assembly. Main results. Dose-rates of up to 128 Gy s−1 at the phantom surface, or 85 Gy s−1 at 1 cm depth, were obtained for a 1 × 1 cm2 field size, 1 mm thick Ta target and 7.5 cm source-to-surface distance using the FLASH-mode beam (E e− = 10 MeV, 2 σ = 5 mm, P = 1 kW); furthermore, removal of the collimation assembly and using a shorter (3.5 cm) SSD afforded dose-rates >600 Gy s−1, albeit at the expense of field conformality. Target temperatures were maintained below the tantalum, aluminum and cooling-water thresholds of 2000 °C, 300 °C and 100 °C, respectively, while the aluminum strain behavior remained everywhere elastic and helped ensure   the converter survives its prescribed 5 yr operational lifetime. Significance. Effective design iteration, target cooling and failure mitigation have culminated in a robust target compatible with intensive transient (FLASH) and steady-state (diagnostic) applications. The ARIEL UHDR photon source will facilitate FLASH-RT experiments concerned with sub-second, pulsed or continuous beam irradiations at dose rates in excess of 40 Gy s−1.

Dose calculations for pre-clinical radiobiology experiments conducted with single-field cabinet irradiators

PURPOSE

To provide percentage depth-dose (PDD) data along the central axis for dosimetry calculations in small-animal radiation biology experiments performed in cabinet irradiators. The PDDs are provided as a function of source-to-surface distance (SSD), field size and animal size.

METHODS

The X-ray tube designs for four biological cabinet irradiators, the RS2000, RT250, MultiRad350 and XRAD320, were simulated using the BEAMnrc Monte Carlo code to generate 160, 200, 250 and 320 kVp photon beams, respectively. The 320 kVp beam was simulated with two filtrations: a soft F1 aluminium filter and a hard F2 thoraeus filter made of aluminium, tin and copper. Beams were collimated into circular fields with diameters of 0.5 - 10 cm at SSDs of 10 - 60 cm. Monte Carlo dose calculations in 1 - 5-cm diameter homogeneous (soft tissue) small-animal phantoms as well as in heterogeneous phantoms with 3-mm diameter cylindrical lung and bone inserts (rib and cortical bone) were performed using DOSXYZnrc. The calculated depth doses in three test-cases were estimated by applying SSD, field size and animal size correction factors to a reference case (40 cm SSD, 1 cm field and 5 cm animal size) and these results were compared with the specifically simulated (i.e., expected) doses to assess the accuracy of this method. Dosimetry for two test-case scenarios of 160 and 250 kVp beams (representative of end-user beam qualities) was also performed, whereby the simulated PDDs at two different depths were compared with the results based on the interpolation from reference data.

RESULTS

The depth doses for three test-cases calculated at 200, 320 kVp F1 and 320 kVp F2, with half value layers (HVL) ranging from ∼0.6 mm to 3.6 mm Cu, agreed well with the expected doses, yielding dose differences of 1.2, 0.1 and 1.0%, respectively. The two end-user test-cases for 160 and 250 kVp beams with respective HVLs of ∼0.8 and 1.8 mm Cu yielded dose differences of 1.4 and 3.2% between the simulated and the interpolated PDDs. The dose increase at the bone-tissue proximal interface ranged from 1.2 to 2.5 times the dose in soft tissue for rib and 1.3 to 3.7 times for cortical bone. The dose drop-off at 1-cm depth beyond the bone ranged from 1.3 - 6.0% for rib and 3.2 - 11.7% for cortical bone. No drastic dose perturbations occurred in the presence of lung, with lung-tissue interface dose of >99% of soft tissue dose and <3% dose increase at 1-cm depth beyond lung.

CONCLUSIONS

The developed dose estimation method can be used to translate the measured dose at a point to dose at any depth in small-animal phantoms, making it feasible for pre-clinical calculation of dose distributions in animals irradiated with cabinet-style irradiators. The dosimetric impact of bone must be accurately quantified as dramatic dose perturbations at and beyond the bone interfaces can occur due to the relative importance of the photoelectric effect at kilovoltage energies. These results will help improve dosimetric accuracy in pre-clinical experiments. This article is protected by copyright. All rights reserved.

Single-Crystal Bismuth Thiophosphate, BiPS4, as a Nontoxic and Mechanically Robust X-ray Detector

Bismuth thiophosphate, BiPS₄, is a promising nontoxic, high-density ternary chalcogenide semiconductor. Polycrystalline BiPS₄ was synthesized from the stoichiometric melt of Bi, P, and S. Phosphorus was purified via an in-situ sublimation method. Single crystals of BiPS₄ were grown using a modified vertical Bridgman method with a thermal gradient of 18 °C/cm. The material exhibits an electrical resistivity of 2 × 10⁹ Ω·cm. The Knoop hardness of the single crystals is 128 ± 0.8 kg mm^-2. A blocking contact behavior was observed with asymmetric contacts of Ga/BiPS₄/Ag. A clear photocurrent response was observed from a BiPS₄ crystal under an electric field as low as 1.14 V mm^-1. Using a tungsten X-ray source, an X-ray sensitivity of 52 μ Gy^-1 cm^-2 was measured under an electric field of 80 V mm^-1. When a single-crystal BiPS₄ radiation detector device was used in a pulse-height radiation detection system, a clear charge collection response was observed under a ²⁴¹Am α-particle radiation source.

Experimental validation of Fastcat kV and MV cone beam CT (CBCT) simulator

PURPOSE

To experimentally validate the Fastcat cone beam computed tomography (CBCT) simulator against kV and MV CBCT images acquired with a Varian Truebeam linac.

METHODS

kV and MV CBCT images of a Catphan 504 phantom were acquired using a 100 kVp beam with the on-board imager (OBI) and a 6 MV treatment beam with the electronic portal imaging device (EPID), respectively. The kV Fastcat simulation was performed using detailed models of the x-ray source, bowtie filter, a high resolution voxelized virtual Catphan phantom, anti-scatter grid, and the CsI scintillating detector. Likewise, an MV Fastcat CBCT was simulated with detailed models for the beam energy spectrum, flattening filter, a high-resolution voxelized virtual Catphan phantom, and the gadolinium oxysulfide (GOS) scintillating detector. Experimental and simulated CBCT images of the phantom were compared with respect to HU values, contrast to noise ratio (CNR), and dose linearity. Detector modulation transfer function (MTF) for the two detectors were also experimentally validated. Fastcat's dose calculations were compared to Monte Carlo (MC) dose calculations performed with Topas.

RESULTS

For the kV and MV simulations, respectively: Contrast agreed within 14 and 9 HUs and detector MTF agreed within 4.2% and 2.5%. Likewise, CNR had a root mean squared error (RMSE) of 2.6% and 1.4%. Dose agreed within 2.4% and 1.6% of MC values. The kV and MV CBCT images took 71 and 72 s to simulate in Fastcat with 887 and 493 projections, respectively.

CONCLUSIONS

We present a multienergy experimental validation of a fast and accurate CBCT simulator against a commercial linac. The simulator is open source and all models found in this work can be downloaded from https://github.com/jerichooconnell/fastcat.git.

Characterization of an x-ray tube-based ultrahigh dose-rate system for in vitro irradiations

PURPOSE

To present an x-ray tube system capable of in vitro ultrahigh dose-rate (UHDR) irradiation of small < 0.3 mm samples and to characterize it by means of a plastic scintillation detector (PSD).

METHODS AND MATERIALS

A conventional x-ray tube was modified for the delivery of short UHDR irradiations. A beam shutter system with a sample holder was designed and installed in a close proximity of an x-ray tube window to enable <1 s irradiations at UHDR. The dosimetry was performed with a small 0.5-mm long 0.5-mm in diameter PSD irradiated with 80, 100, and 120 kVp beams and beam currents of 1-37.5 mA. The PSD signal was recorded at frame rates of 20 and 50 fps for shutter exposure between 100 and 1125 ms. Irradiation reproducibility was studied with the PSD. The x-ray tube irradiation setup was modeled with Monte Carlo (MC) and dose on a surface of a phantom was also measured with films. The effect of dose delivery uncertainty to 300-μm spheroids due to positioning and spheroid size was evaluated.

RESULTS

MC simulations showed good agreement with PSD measurements acquired at both frame rates of 20 and 50 fps in terms of beam temporal profile. PSD-measured dose exhibited excellent linearity as a function of instantaneous dose rate from 3.1 to 118.0 Gy/s as well as shutter exposure time from 100 and 1125 ms for all investigated beam energies. PSD absorbed dose for the 80, 100, and 120 kVp beams agreed with MC simulations to within 5%. The total delivered doses ranged from 0.4 Gy for a 1-mA, 80 kVp beam, and 100 ms shutter exposure to 166.9 Gy for a 37.5-mA, 80 kVp beam, and a 1125 ms exposure. PSD irradiation reproducibility was < 0.5%. Simulated and measured dose fall off agreed and it was steep along the axis of the shutter slit (1%/0.1 mm) and with depth (2%/0.1 mm at 1-mm depth). Spheroid positioning uncertainty of 300 μm resulted in dose difference of < 3% for x and y shifts but up to 7% uncertainty for a z-shift parallel to the beam axis. A 16% difference in spheroid size resulted in <5% dose difference in spheroid absorbed dose.

CONCLUSIONS

We have presented a cost-effective x-ray tube-based system with a beam shutter designed for in vitro UHDR delivery and reaching dose rates of up to 118.0 Gy/s. The described shutter system can be easily implemented at other institutions, which might enable new researchers to investigate the radiobiology of UHDR irradiations in vitro.

FLASH radiotherapy with photon beams

Ultra-high-dose rate "FLASH" radiotherapy (FLASH-RT) has been shown to drastically reduce normal tissue toxicities while being as efficacious as conventional dose rate radiotherapy to treat tumors. A large number of preclinical studies describing this so-called FLASH effect have led to the clinical translation of FLASH-RT using ultra-high-dose rate electron and proton beams. Although the vast majority of radiation therapy treatments are delivered using X-rays, few preclinical data using ultra-high-dose rate X-ray irradiation have been published. This review focuses on different methods that can be used to generate ultra-high-dose rate X-rays and their beam characteristics along with their effect on the biological tissues and the perspectives for the development of FLASH-RT with X-rays.

A detective quantum efficiency for spectroscopic X-ray imaging detectors

PURPOSE

Spectroscopic X-ray detectors (SXDs) are under development for X-ray imaging applications. Recent efforts to extend the detective quantum efficiency (DQE) to SXDs impose a barrier to experimentation and/or do not provide a task-independent measure of detector performance. The purpose of this article is to define a task-independent DQE for SXDs that can be measured using a modest extension of established DQE-metrology methods.

METHODS

We defined a task-independent spectroscopic DQE and performed a simulation study to determine the relationship between the zero-frequency DQE and the ideal-observer signal-to-noise ratio (SNR) of low-frequency soft-tissue, bone, iodine, and gadolinium signals. In our simulations, we used calibrated models of the spatioenergetic response of cadmium telluride (CdTe) and cadmium-zinc-telluride (CdZnTe) SXDs. We also measured the zero-frequency DQE of a CdTe detector with two energy bins and of a CdZnTe detector with up to six energy bins for an RQA9 spectrum and compared with model predictions.

RESULTS

The spectroscopic DQE accounts for spectral distortions, energy-bin-dependent spatial resolution, interbin spatial noise correlations, and intrabin spatial noise correlations; it is mathematically equivalent to the squared SNR per unit fluence of the generalized least-squares estimate of the height of an X-ray impulse in a uniform noisy background. The zero-frequency DQE has a strong linear relationship with the ideal-observer SNR of low-frequency soft-tissue, bone, iodine, and gadolinium signals, and can be expressed in terms of the product of the quantum efficiency and a Swank noise factor that accounts for DQE degradation due to, for example, charge sharing (CS) and electronic noise. The spectroscopic Swank noise factor of the CdTe detector was measured to be 0.81 ± 0.04 and 0.83 ± 0.04 with and without anticoincidence logic for CS suppression, respectively. The spectroscopic Swank noise factor of the CdZnTe detector operated with four energy bins was measured to be 0.82 ± 0.02 which is within 5% of the theoretical value.

CONCLUSIONS

The spectroscopic DQE defined here is (1) task-independent, (2) can be measured using a modest extension of existing DQE-metrology methods, and (3) is predictive of the ideal-observer SNR of soft-tissue, bone, iodine, and gadolinium signals. For CT applications, the combination of CS and electronic noise in CdZnTe spectroscopic detectors will degrade the zero-frequency DQE by 10 %-20 % depending on the electronic noise level and pixel size.

Dosimetry of a novel converging X-ray source for kilovoltage radiotherapy

PURPOSE

The objective of this work was to evaluate phantom dosimetry of a novel kilovoltage (kV) X-ray source, which employs a stationary tungsten anode and a linearly swept scanning electron beam. The source utilizes converging X-ray collimation along with orthogonal mechanical rotation to distribute surface flux over large area. In this study, this was investigated as a potential solution to fast-falloff limitations expected with kV radiotherapy. This was done with the aim of future clinical development of a lower cost radiotherapy alternative to megavoltage (MV) linac systems.

METHODS

Radiochromic film was employed for dosimetry on the kV X-ray source of the linear-converging radiotherapy system (LCRS). The source utilizes charge particle optics to magnetically deflect and focus an electron beam along a stationary, reflection tungsten target in an ultra-high-vacuum stainless-steel chamber. Resulting X-rays were collimated into converging beamlets that span a large planar angle and converge at the system isocenter. In this study, radiochromic film dosimetry was done at 140 and 145 kVp for a designated planning treatment volume (PTV) of 4 cm diameter. An acrylic phantom was employed for dose distribution measurements of stationary and rotational delivery. Film dosimetry was evaluated in planes parallel to the source X-ray window at various depths, as well as in the plane of gantry rotation.

RESULTS

At 140 and 145 kVp and using a collimated 4 cm square field at depth, lesion-to-skin dose ratio was shown to improve with additional beams from different relative source positions, where the different beams are focused at the same isocenter and do not overlap at the phantom surface. It was only possible to achieve a 1:1 D_max -to-surface ratio with four delivery beams, but the ratio improved to 4:1 with 12 beams, focused at the same isocenter depth of 7.8 cm in an acrylic phantom. For the tests conducted, the following D_max -to-surface ratios were obtained: 0.4:1 lesion-to-skin ratio for stationary delivery from one entry beam, 0.71:1 lesion-to-skin ratio was obtained for two beams, 1.07:1 ratio for four beams, and 4:1 for 12 beams. Dose-depth profiles were evaluated for stationary and rotational dosimetry. Additionally, rotational dosimetry was measured for a case more analogous to a clinical scenario, where the isocenter was located at an off-center simulated lesion.

CONCLUSIONS

The results demonstrate potential dose-depth improvements with kV arc therapy by distributing the surface flux with a wide converging beam along with perpendicular mechanical source rotation of the LCRS. The system delivered tolerable dose to a large surface area when a threshold of multiple, separated beams was reached. The radiochromic film data support the feasibility of the construct of the LCRS kV radiotherapy system design.

fastCAT: Fast cone beam CT (CBCT) simulation

PURPOSE

To develop fastCAT, a fast cone-beam computed tomography (CBCT) simulator. fastCAT uses pre-calculated Monte Carlo (MC) CBCT phantom-specific scatter and detector response functions to reduce simulation time for megavoltage (MV) and kilovoltage (kV) CBCT imaging.

METHODS

Pre-calculated x-ray beam energy spectra, detector optical spread functions and energy deposition, and phantom scatter kernels are combined with GPU raytracing to produce CBCT volumes. MV x-ray beam spectra are simulated with EGSnrc for 2.5- and 6 MeV electron beams incident on a variety of target materials and kV x-ray beam spectra are calculated analytically for an x-ray tube with a tungsten anode. Detectors were modeled in Geant4 extended by Topas and included optical transport in the scintillators. Two MV detectors were modeled-a standard Varian AS1200 GOS detector and a novel CWO high detective quantum efficiency detector. A kV CsI detector was also modeled. Energy-dependent scatter kernels were created in Topas for two 16 cm diameter phantoms: A Catphan 515 contrast phantom and an anthropomorphic head phantom. The Catphan phantom contained inserts of 1-5 mm in diameter of six different tissue types: brain, deflated lung, compact and cortical bone, adipose, and B-100.

RESULTS

fastCAT simulations retain high fidelity to measurements and MC simulations: MTF curves were within 3.5% and 1.2% of measured values for the CWO and GOS detectors, respectively. HU values and CNR in a fastCAT Catphan 515 simulation were seen to be within 95% confidence intervals of an equivalent MC simulation for all of the tissues with root mean squared errors less than 16 HU and 1.6 in HU values and CNR comparisons, respectively. The anthropomorphic head phantom CWO detector CBCT image resulted in much higher tissue contrast and lower noise compared to the GOS detector CBCT image. A fastCAT simulation of the Catphan 515 module with an image size of 1024  ×  1024  ×  10 voxels took 61 s on a GPU while the equivalent Topas MC was estimated to take more than 0.3 CPU years.

CONCLUSIONS

We present an open source fast CBCT simulation with high fidelity to MC simulations. The fastCAT python package can be found at https://github.com/jerichooconnell/fastCAT.git.

High length-to-width aspect ratio lead bromide microwires via perovskite-induced local concentration gradient for X-ray detection

Well-oriented PbBr₂ microwires with a length-to-width ratio of up to 5000 were grown using a concentration gradient in co-crystallization with perovskite. Planar-integrated microwires showed a response to X-ray photons.

Physics and biology of ultrahigh dose-rate (FLASH) radiotherapy: a topical review

Ultrahigh dose-rate radiotherapy (RT), or 'FLASH' therapy, has gained significant momentum following various in vivo studies published since 2014 which have demonstrated a reduction in normal tissue toxicity and similar tumor control for FLASH-RT when compared with conventional dose-rate RT. Subsequent studies have sought to investigate the potential for FLASH normal tissue protection and the literature has been since been inundated with publications on FLASH therapies. Today, FLASH-RT is considered by some as having the potential to 'revolutionize radiotherapy'. FLASH-RT is considered by some as having the potential to 'revolutionize radiotherapy'. The goal of this review article is to present the current state of this intriguing RT technique and to review existing publications on FLASH-RT in terms of its physical and biological aspects. In the physics section, the current landscape of ultrahigh dose-rate radiation delivery and dosimetry is presented. Specifically, electron, photon and proton radiation sources capable of delivering ultrahigh dose-rates along with their beam delivery parameters are thoroughly discussed. Additionally, the benefits and drawbacks of radiation detectors suitable for dosimetry in FLASH-RT are presented. The biology section comprises a summary of pioneering in vitro ultrahigh dose-rate studies performed in the 1960s and early 1970s and continues with a summary of the recent literature investigating normal and tumor tissue responses in electron, photon and proton beams. The section is concluded with possible mechanistic explanations of the FLASH normal-tissue protection effect (FLASH effect). Finally, challenges associated with clinical translation of FLASH-RT and its future prospects are critically discussed; specifically, proposed treatment machines and publications on treatment planning for FLASH-RT are reviewed.

Monte Carlo simulations of EBT3 film dose deposition for percentage depth dose (PDD) curve evaluation

Purpose

To use Monte Carlo (MC) calculations to evaluate the effects of Gafchromic EBT3 film orientation on percentage depth dose (PDD) curves.

Methods

Dose deposition in films placed in a water phantom, and oriented either parallel or perpendicular with respect to beam axis, were simulated with MC and compared to PDDs scored in a homogenous water phantom. The effects of introducing 0.01–1.00 mm air gaps on each side of the film as well as a small 1°‐3° tilt for film placed in parallel orientation were studied. PDDs scored based on two published EBT3 film compositions were compared. Three photon beam energies of 120 kVp, 220 kVp, and 6 MV and three field sizes between 1 × 1 and 5 × 5 cm² were considered. Experimental PDDs for a 6‐MV 3 × 3 cm² beam were acquired.

Results

PDD curves for films in perpendicular orientation more closely agreed to water PDDs than films placed in parallel orientation. The maximum difference between film and water PDD for films in parallel orientation was −12.9% for the 220 kVp beam. For the perpendicular film orientation, the maximum difference decreased to 5.7% for the 120 kVp beam. The inclusion of an air gap had the largest effect on the 6‐MV 1 × 1 cm² beam, for which the dose in the buildup region was underestimated by 21.2% compared to the simulation with no air gap. A 2° film tilt decreased the difference between the parallel film and homogeneous water phantom PDDs from −5.0% to −0.5% for the 6 MV 3 × 3 cm² beam. The “newer” EBT3 film composition resulted in larger PDD discrepancies than the previous composition. Experimental film data qualitatively agreed with MC simulations.

Conclusions

PDD measurements with films should either be performed with film in perpendicular orientation to the beam axis or in parallel orientation with a ~ 2º tilt and no air gaps.

Technical Note: Synthesizing of lung tumors in computed tomography images

PURPOSE

When investigating new radiation therapy techniques in the treatment planning stage, it can be extremely time consuming to locate multiple patient scans that match the desired characteristics for the treatment. With the help of machine learning, we propose to bypass the difficulty in finding patient computed tomography (CT) scans that match the treatment requirements. Furthermore, we aim to provide the developed method as a tool that is easily accessible to interested researchers.

METHODS

We propose a generative adversarial network (GAN) to edit individual volumes of interest (VOIs) in pre-existing CT scans, translating features of the healthy VOIs into features of cancerous volumes. Training and testing was done using VOIs from a dataset of 460 diagnostic and lung cancer screening CT scans. Agreement between real tumors and those produced by the editor was tested by comparing the distributions of several histogram parameters and second-order statistics as well as using qualitative analysis.

RESULTS

After training, the network was successfully able to map healthy CT segments to realistic looking cancerous volumes. Based on visual inspection, tumors produced by the editor were found to be both realistic and visually consistent with the surrounding anatomy when placed back into the original CT scan. Furthermore, the network was found to be able to extrapolate well beyond the upper size limit of the training set. Lastly, a graphical user interface (GUI) was developed to easily interact with the resulting network.

CONCLUSION

The trained network and associated GUI can serve as a tool to develop an abundance of lung cancer patient data to be used in treatment planning. In addition, this method can be extended to a variety of cancer types if given an appropriate baseline dataset. The GUI and instructions on how to utilize the tool have been made publicly available at https://github.com/teaghan/CT_Editor.

Optimization of radiochromic film stacks to diagnose high-flux laser-accelerated proton beams

Here, we extend flatbed scanner calibrations of GafChromic EBT3, MD-V3, and HD-V2 radiochromic films using high-precision x-ray irradiation and monoenergetic proton bombardment. By computing a visibility parameter based on fractional errors, optimal dose ranges and transitions between film types are identified. The visibility analysis is used to design an ideal radiochromic film stack for the proton energy spectrum expected from the interaction of a petawatt laser with a cryogenic hydrogen jet target.

Photon-counting computed tomography of lanthanide contrast agents with a high-flux 330-μm-pitch cadmium zinc telluride detector in a table-top system

Purpose: We present photon-counting computed tomography (PCCT) imaging of contrast agent triplets similar in atomic number ( Z ) achieved with a high-flux cadmium zinc telluride (CZT) detector. Approach: The table-top PCCT imaging system included a 330 - μ m -pitch CZT detector of size 8    mm × 24    mm 2 capable of using six energy bins. Four 3D-printed 3-cm-diameter phantoms each contained seven 6-mm-diameter vials with water and low and high concentration solutions of various contrast agents. Lanthanum ( Z = 57 ), gadolinium (Gd) ( Z = 64 ), and lutetium ( Z = 71 ) were imaged together and so were iodine ( Z = 53 ), Gd, and holmium ( Z = 67 ). Each phantom was imaged with 1-mm aluminum-filtered 120-kVp cone beam x rays to produce six energy-binned computed tomography (CT) images. Results: K -edge images were reconstructed using a weighted sum of six CT images, which distinguished each contrast agent with a root-mean-square error (RMSE) of < 0.29 % and 0.51% for the 0.5% and 5% concentrations, respectively. Minimal cross-contamination in each K -edge image was seen, with RMSE values < 0.27 % in vials with no contrast. Conclusion: This is the first preliminary demonstration of simultaneously imaging three similar Z contrast agents with a difference in Z as low as 3.

A failure modes and effects analysis quality management framework for image-guided small animal irradiators: A change in paradigm for radiation biology

PURPOSE

Image-guided small animal irradiators (IGSAI) are increasingly being adopted in radiation biology research. These animal irradiators, designed to deliver radiation with submillimeter accuracy, exhibit complexity similar to that of clinical radiation delivery systems, including image guidance, robotic stage motion, and treatment planning systems. However, physics expertise and resources are scarcer in radiation biology, which makes implementation of conventional prescriptive QA infeasible. In this study, we apply the failure modes and effect analysis (FMEA) popularized by the AAPM task group 100 (TG-100) report to IGSAI and radiation biological research.

METHODS

Radiation biological research requires a change in paradigm where small errors to large populations of animals are more severe than grievous errors that only affect individuals. To this end, we created a new adverse effects severity table adapted to radiation biology research based on the original AAPM TG-100 severity table. We also produced a process tree which outlines the main components of radiation biology studies performed on an IGSAI, adapted from the original clinical IMRT process tree from TG-100. Using this process tree, we created and distributed a preliminary survey to eight expert IGSAI operators in four institutions. Operators rated proposed failure modes for occurrence, severity, and lack of detectability, and were invited to share their own experienced failure modes. Risk probability numbers (RPN) were calculated and used to identify the failure modes which most urgently require intervention.

RESULTS

Surveyed operators indicated a number of high (RPN >125) failure modes specific to small animal irradiators. Errors due to equipment breakdown, such as loss of anesthesia or thermal control, received relatively low RPN (12-48) while errors related to the delivery of radiation dose received relatively high RPN (72-360). Errors identified could either be improved by manufacturer intervention (e.g., electronic interlocks for filter/collimator) or physics oversight (errors related to tube calibration or treatment planning system commissioning). Operators identified a number of failure modes including collision between the collimator and the stage, misalignment between imaging and treatment isocenter, inaccurate robotic stage homing/translation, and incorrect SSD applied to hand calculations. These were all relatively highly rated (90-192), indicating a possible bias in operators towards reporting high RPN failure modes.

CONCLUSIONS

The first FMEA specific to radiation biology research was applied to image-guided small animal irradiators following the TG-100 methodology. A new adverse effects severity table and a process tree recognizing the need for a new paradigm were produced, which will be of great use to future investigators wishing to pursue FMEA in radiation biology research. Future work will focus on expanding scope of user surveys to users of all commercial IGSAI and collaborating with manufacturers to increase the breadth of surveyed expert operators.

X-Ray Fluorescence Computed Tomography Induced by Photon, Electron, and Proton Beams

X-ray fluorescence CT (XFCT) has shown promise for molecular imaging of gold nanoparticles. To date, XFCT has been induced by kilovoltage photon beams due to the high photoelectric interaction probability. We compare K-shell and L-shell XFCT induced by photon, electron, and proton beams for two phantom sizes. A 2.5 and 5.0-cm diameter phantom with four 5 mm and 10 mm vials, respectively, with gold-solutions of 0.1%-2% by weight were built in TOPAS, a GEANT4-based Monte Carlo simulation tool. The 2.5-cm phantom was imaged with XFCT induced by beams of 7.45×10⁴ 81 keV- and 5 MeV-photons, 220 kVp- and 6 MV-photons, 10 MeV- and 100 MeV-electrons, and 100 MeV- and 250 MeV-protons. The doses between each phantom size were equal. First-generation CT geometry with 0.5 mm × 0.5 mm pencil beams with 0.5 mm-translation and 2°-rotation steps over each phantom was modeled. The scattered x-rays were detected on an idealized spherical detector from which the K-shell and L-shell fluorescent x-rays were extracted in 0.5 keV and 0.2 keV bins. XFCT images were generated using iterative reconstruction algorithms. The highest gold sensitivity was seen in the 81 keV-photon K-shell and L-shell images (0.004% and 0.007%) of the 5.0 cm-phantom at 30 mGy. For the 2.5 cm-phantom, the detection limits were 0.006%, 0.62%, and 0.28% for 81 keV-photon K-shell, 100 MeV-electron K-shell, and 100 MeV-proton L-shell images, respectively. The mean imaging dose was approximately 2-3 orders of magnitude higher in electron- and proton-XFCT compared to 81keV-photon XFCT. Our MC study demonstrates that the small-object XFCT imaging achieves the best performance when induced with kilovoltage-photon beams. Due to high imaging doses, electron- and proton-induced XFCT might be feasible for guiding nanoparticle-enhanced charged-particle radiotherapy.

Preclinical dose verification using a 3D printed mouse phantom for radiobiology experiments

PURPOSE

Dose verification in preclinical radiotherapy is often challenged by a lack of standardization in the techniques and technologies commonly employed along with the inherent difficulty of dosimetry associated with small-field kilovoltage sources. As a consequence, the accuracy of dosimetry in radiobiological research has been called into question. Fortunately, the development and characterization of realistic small-animal phantoms has emerged as an effective and accessible means of improving dosimetric accuracy and precision in this context. The application of three-dimensional (3D) printing, in particular, has enabled substantial improvements in the conformity of representative phantoms with respect to the small animals they are modeled after. In this study, our goal was to evaluate a fully 3D printed mouse phantom for use in preclinical treatment verification of sophisticated therapies for various anatomical targets of therapeutic interest.

METHODS

An anatomically realistic mouse phantom was 3D printed based on segmented microCT data of a tumor-bearing mouse. The phantom was modified to accommodate both laser-cut EBT3 radiochromic film within the mouse thorax and a plastic scintillator dosimeter (PSD), which may be placed within the brain, abdomen, or 1-cm flank subcutaneous tumor. Various treatments were delivered on an image-guided small-animal irradiator in order to determine the doses to isocenter using a PSD and validate lateral- and depth-dose distributions using film dosimeters. On-board cone-beam CT imaging was used to localize isocenter to the film plane or PSD active element prior to irradiation. The PSD irradiations comprised a 3 × 3 mm² brain arc, 5 × 5 mm² parallel-opposed pair (POP), and 5-beam 10 × 10 mm² abdominal coplanar arrangement while two-dimensional (2D) film dose distributions were acquired using a 3 × 3 mm² arc and both 5 × 5 and 10 × 10 mm² 3-beam coplanar plans. A validated Monte Carlo (MC) model of the source was used as to verify the accuracy of the film and PSD dose measurements. computer-aided design (CAD) geometries for the mouse phantom and dosimeters were imported directly into the MC code to allow for highly accurate reproduction of the physical experiment conditions. Experimental and MC-derived film data were co-registered and film dose profiles were compared for points above 90% of the dose maximum. Point dose measurements obtained with the PSD were similarly compared for each of the candidate (brain, abdomen, and tumor) treatment sites.

RESULTS

For each treatment configuration and anatomical target, the MC-calculated and measured doses met the proposed 5% agreement goal for dose accuracy in radiobiology experiments. The 2D film and MC dose distributions were successfully registered and mean doses for lateral profiles were found to agree to within 2.3% in all cases. Isocentric point-dose measurements taken with the PSD were similarly consistent, with a maximum percentage deviation of 3.2%.

CONCLUSIONS

Our study confirms the utility of 3D printed phantom design in providing accurate dose estimates for a variety of preclinical treatment paradigms. As a tool for pretreatment dose verification, the phantom may be of particular interest to researchers for its ability to facilitate precise dosimetry while fostering a reduction in cost for radiobiology experiments.

A 3D printed modular phantom for quality assurance of image-guided small animal irradiators: Design, imaging experiments, and Monte Carlo simulations

PURPOSE

The goal of this work was to develop and test a cylindrical tissue-equivalent quality assurance (QA) phantom for micro computed tomography (microCT) image-guided small animal irradiators that overcomes deficiencies of existing phantoms due to its mouse-like dimensions and composition.

METHODS

The 8.6-cm-long and 2.4-cm-diameter phantom was three-dimensionally (3D) printed out of Somos NeXt plastic on a stereolithography (SLA) printer. The modular phantom consisted of four sections: (a) CT number evaluation section, (b) spatial resolution with slanted edge (for the assessment of longitudinal resolution) and targeting section, (c) spatial resolution with hole pattern (for the assessment of radial direction) section, and (d) uniformity and geometry section. A Python-based graphical user interface (GUI) was developed for automated analysis of microCT images and evaluated CT number consistency, longitudinal and radial modulation transfer function (MTF), image uniformity, noise, and geometric accuracy. The phantom was placed at the imaging isocenter and scanned with the small animal radiation research platform (SARRP) in the pancake geometry (long axis of the phantom perpendicular to the axis of rotation) with a variety of imaging protocols. Tube voltage was set to 60 and 70 kV, tube current was set to 0.5 and 1.2 mA, voxel size was set to 200 and 275 μm, imaging times of 1, 2, and 4 min were used, and frame rates of 6 and 12 frames per second (fps) were used. The phantom was also scanned in the standard (long axis of the phantom parallel to the axis of rotation) orientation. The quality of microCT images was analyzed and compared to recommendations presented in our previous work that was derived from a multi-institutional study. Additionally, a targeting accuracy test with a film placed in the phantom was performed. MicroCT imaging of the phantom was also simulated in a modified version of the EGSnrc/DOSXYZnrc code. Images of the resolution section with the hole pattern were acquired experimentally as well as simulated in both the pancake and the standard imaging geometries. The radial spatial resolution of the experimental and simulated images was evaluated and compared to experimental data.

RESULTS

For the centered phantom images acquired in the pancake geometry, all imaging protocols passed the spatial resolution criterion in the radial direction (>1.5 lp/mm @ 0.2 MTF), the geometric accuracy criterion (<200 μm), and the noise criterion (<55 HU). Only the imaging protocol with 200-μm voxel size passed the criterion for spatial resolution in the longitudinal direction (>1.5 lp/mm @ 0.2 MTF). The 70-kV tube voltage dataset failed the bone CT number consistency test (<55 HU). Due to cupping artifacts, none of the imaging protocols passed the uniformity test of <55 HU. When the phantom was scanned in the standard imaging geometry, image uniformity and longitudinal MTF were satisfactory; however, the CT number consistency failed the recommended limit. A targeting accuracy of 282 and 251 μm along the x- and z-direction was observed. Monte Carlo simulations confirmed that the radial spatial resolution for images acquired in the pancake geometry was higher than the one acquired in the standard geometry.

CONCLUSIONS

The new 3D-printed phantom presents a useful tool for microCT image analysis as it closely mimics a mouse. In order to image mouse-sized animals with acceptable image quality, the standard protocol with a 200-μm voxel size should be chosen and cupping artifacts need to be resolved.

[18F]-SuPAR: A Radiofluorinated Probe for Noninvasive Imaging of DNA Damage-Dependent Poly(ADP-ribose) Polymerase Activity

Poly(ADP ribose) polymerase (PARP) enzymes generate poly(ADP ribose) post-translational modifications on target proteins for an array of functions centering on DNA and cell stress. PARP isoforms 1 and 2 are critically charged with the surveillance of DNA integrity and are the first line guardians of the genome against DNA breaks. Here we present a novel probe ([¹⁸F]-SuPAR) for noninvasive imaging of PARP-1/2 activity using positron emission tomography (PET). [¹⁸F]-SuPAR is a radiofluorinated nicotinamide adenine dinucleotide (NAD) analog that can be recognized by PARP-1/2 and incorporated into the long branched polymers of poly(ADP ribose) (PAR). The measurement of PARP-1/2 activity was supported by a reduction of radiotracer uptake in vivo following PARP-1/2 inhibition with talazoparib treatment, a potent PARP inhibitor recently approved by FDA for treatment of breast cancer, as well as ex vivo colocalization of radiotracer analog and poly(ADP ribose). With [¹⁸F]-SuPAR, we were able to map the dose- and time-dependent activation of PARP-1/2 following radiation therapy in breast and cervical cancer xenograft mouse models. Tumor response to therapy was determined by [¹⁸F]-SuPAR PET within 8 h of administration of a single dose of radiation equivalent to one round of stereotactic ablative radiotherapy.

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/cm³ ) and bone (ρ=1.20 g/cm³ ). The excavated lungs were filled with polystyrene (ρ=0.32 g/cm³ ). 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.

Optimization of a table-top x-ray fluorescence computed tomography (XFCT) system

Pencil beam x-ray fluorescence computed tomography (XFCT) has typically used a single spectrometer and prohibitively long scan times. However, detecting backscattered fluorescent x-rays from gold nanoparticles (AuNPs) using multiple spectrometers greatly reduces image noise and scan time. The arrangement of eight spectrometers for combined K-shell and L-shell XFCT was investigated along with a variety of conditions. A 2.5 cm-diameter cylindrical water phantom containing 4 mm-diameter vials with 0.1%-2% AuNP concentrations by weight was modeled by TOPAS, a GEANT4-based Monte Carlo software. The phantom was irradiated to 30 mGy by a 0.5 mm Pb-filtered 120 kVp and 1 mm Al-filtered 30 kVp 1 mm² x-ray pencil beam to yield respective Au K-shell and L-shell fluorescent x-rays, with 50 0.5 mm translation and 2-degree rotation steps. Eight CdTe and silicon drift detector (SDD) spectrometers were placed 2.25 cm away from the isocentre. The respective energy resolution was applied to the detected energy spectra and the spectra were corrected for detector response before extracting the fluorescence signal. Three CdTe and SDD spectrometer configurations (isotropic/backscattered grid/backscattered row arrangements), two CdTe crystal sizes (9 mm²/25 mm²), two scanning techniques (moving/stationary spectrometers) and five vial-edge depths (0-4 mm) were considered in optimizing the contrast-to-noise ratio (CNR) for each XFCT image reconstructed with a maximum-likelihood expectation maximization (MLEM) algorithm. The isotropic spectrometer arrangement had AuNP detection sensitivities of 0.106% for K-shell and 0.132% for L-shell XFCT at 4 mm depth. Comparatively, the backscattered grid arrangement had the best AuNP sensitivity of 0.055% and 0.095%. The highest K-shell (0.044%) and L-shell (0.004%) AuNP sensitivities were found for vials at 0 mm depth. Using stationary spectrometers or the 9 mm² CdTe crystal compromised the CNR. For the best-case arrangement, L-shell XFCT is superior at vial-edge depths less than 3.0 mm. This work demonstrated the importance of spectrometer arrangement and vial depth for improving AuNP sensitivity and will guide the design for our table-top XFCT system.

Characterization of a plastic scintillating detector for the Small Animal Radiation Research Platform (SARRP)

PURPOSE

The purpose of this study was to characterize a small plastic scintillator developed for high resolution, real-time dosimetry of therapy and imaging x-ray beams delivered by an image-guided small animal irradiator.

MATERIALS AND METHODS

A 1 mm diameter, 1 mm long polystyrene BCF-60 scintillating fiber dosimeter was characterized with 220 kVp therapy and 40, 50, 60, 70, and 80 kVp imaging beams on the Small Animal Research Platform (SARRP). Scintillator output, sensitivity (charge per unit dose), linearity, and 0.2-mm resolution beam profile measurements were performed. A validated in-house Monte Carlo (MC) model of the SARRP was used to compute detailed energy spectra at locations of dosimetry, and validated scintillator measurement with MC simulations. Mass energy-absorption coefficients from the National Institute of Standards and Technology (NIST) tables convolved with MC-derived spectra were used in conjunction with Birks ionization quenching factors to correct scintillator output. An air kerma calibration method was employed to correct scintillator output for in-air beam profile measurements with open, 5 × 5, and 3 × 3 mm² square field sizes, and compared to MC simulations.

RESULTS

Scintillator dose response showed excellent linearity (R²  ≥ 0.999) for all sensitivity measurements, including output as a function of tube current. Detector sensitivity was 2.41 μC Gy^-1 for the 220 kVp therapy beam, and it ranged from 1.21 to 1.32 μC Gy^-1 for the 40-80 imaging beams. Percentage difference in sensitivity between the therapy and imaging beams before sensitivity correction and after using the Birks quenching factors were 52.3% and 10.2%, respectively. Percentage differences between the therapy and imaging beam sensitivities after using the air kerma calibration method for in-air measurements was excellent and below 0.3%. In-air beam profile measurements agreed to MC simulations within a mean difference of 2.4% for the 5 × 5 and 3 × 3 mm² field sizes, however, the scintillator showed signs of volume averaging at the penumbra edges.

CONCLUSIONS

A small plastic scintillator was characterized for therapy and imaging energies of a small animal irradiator, with output corrected for using an in-house MC model of the irradiator. The characterization of the scintillator detector system for small fields presents steps toward implementing real-time measurements for quality assurance and small animal treatment and imaging dose verification.

Inverse optimization of low-cost kilovoltage x-ray arc therapy plans

PURPOSE

The objective of this work was to investigate the benefits of using inverse optimization treatment planning for kilovoltage arc therapy (KVAT) and to assess the dosimetric limitations of KVAT.

METHODS

Monte Carlo (MC) calculated, inversely optimized KVAT plans of spherical, idealized breast, lung, and prostate lesions were calculated using the EGSnrc/BEAMnrc and DOSXYZnrc MC codes. The dose delivered with the KVAT system, which generates 200-225 kV photon beamlets, was calculated and inversely optimized using an optimization framework developed at McGill University. KVAT dose distributions were compared with inversely optimized and MC generated megavoltage (MV) volumetric modulated arc therapy (VMAT) plans as a reference. Prescription doses delivered to 95% of the planning target volume (PTV) were 38.5 (10 fractions), 60 (30 fractions) and 73.8 (41 fractions) Gy for the breast, lung and prostate patients, respectively. Dose distributions, dose volume histograms, and PTV homogeneity indices were used to evaluate KVAT and VMAT plans based on RTOG protocols.

RESULTS

All organ-at-risk (OAR) doses were within prescribed dose limits for KVAT and VMAT plans. Generally, KVAT plans delivered higher doses to OARs. For example, due to the lower energy of KVAT, 50% of the rib volume received 12.9 Gy from KVAT while only receiving 2.5 Gy from VMAT. OAR doses were especially high for the KVAT prostate plan due to the presence of large volumes of bony anatomy, which illustrates a limitation of the KVAT system. The KVAT treatment times per fraction for the breast, lung and prostate patients were 2.8, 2.6 and 5.5 min, respectively.

CONCLUSIONS

The inversely optimized KVAT plans presented in this work have demonstrated the ability of our novel low-cost, kilovoltage x-ray therapy system to safely treat deep-seated spherical lesions in breast and lung patients while meeting RTOG dose constraints on OARs.

MicroCT imaging dose to mouse organs using a validated Monte Carlo model of the small animal radiation research platform (SARRP)

The goal of this work was to establish imaging dose to mouse organs with a validated Monte Carlo (MC) model of the image-guided Small Animal Radiation Research Platform (SARRP) and to investigate the effect of scatter from the internal walls on animal therapy dose determination. A MC model of the SARRP was built in the BEAMnrc code and validated with a series of homogeneous and heterogeneous phantom measurements. A segmented microCT scan of a mouse was used in DOSXYZnrc to determine mouse organ microCT imaging doses to 15-35 g mice for the SARRP pancake (mouse lying on couch) and standard (mouse standing on couch) imaging geometries for 40-80 kVp tube voltages. Imaging dose for off-center positioning shifts and maintaining image noise across tube voltages were also calculated. Half-value layer (HVL) measurements for the 220 kVp therapy beam in the presence of the SARRP shielding cabinet were modeled in BEAMnrc and compared to the 100 cm source-to-detector distance (SDD) in the scatter free, narrow-beam geometry recommended by the American Association of Physicists in Medicine Task Group 61 (AAPM TG-61). For a 60 kVp, 0.8 mA, and 60 s scan protocol, maximum mean organ imaging doses to boney and non-boney structures were 10.5 cGy and 3.5 cGy, respectively, for an average size 20 g mouse. Current-exposure combinations above 323, 203, 147, 116, and 95 mAs for 40-80 kVp tube voltages, respectively, will increase body doses above 10 cGy. MicroCT mean body dose was 18% lower in pancake compared to standard imaging geometry. An 11% difference in measured HVL at a 50 cm SDD was found compared to MC simulated HVL for the AAPM TG-61 recommended scatter free geometry at a 100 cm SDD. This change in HVL resulted in a 0.5% change in absorbed dose to water calculations for the treatment beam.

Sheet beam x-ray fluorescence computed tomography (XFCT) imaging of gold nanoparticles

PURPOSE

X-ray fluorescence computed tomography (XFCT) experiments have typically used pencil beams for data acquisition, which yielded good quality images of gold nanoparticles (AuNP) but prolonged the imaging time. Here we propose three novel collimator geometries for use with faster sheet beam XFCT data acquisition. The feasibility of a multipinhole, parallel, and converging collimator was investigated in a Monte Carlo study.

METHODS

A cylindrical water phantom with 2 cm in diameter and 3 cm in height containing 0.5-2 mm diameter vials with 0.4%-1.6% AuNP concentrations was modelled by FLUKA. A 15 and 81 keV monoenergetic x-ray sheet beam of 0.4 mm in width was used to image the phantom with L-shell and K-shell XFCT, respectively, with a dose of 30 mGy. The collimator thickness for L-shell and K-shell data acquisition was 3.3 and 5.1 mm, respectively. The XFCT images resulting from three collimator geometries were generated using the maximum likelihood expectation maximization (MLEM) iterative reconstruction method. With a resolution of 0.4 mm they were corrected for x-ray attenuation. The sheet beam XFCT images were compared against pencil beam geometry images that were generated using 55 translations. To assess image quality, the contrast-to-noise ratio (CNR) was evaluated for each vial. The Rose criterion was used to determine the lowest AuNP concentration detectable for each image.

RESULTS

Among the three collimator geometry types, the sheet beam L-shell and K-shell parallel collimator XFCT images yielded AuNP sensitivity limits at 0.09% and 0.08%, respectively, for a 2 mm diameter vial. The AuNP sensitivity limits of the pencil beam XFCT images were 0.07% and 0.01% for L-shell and K-shell XFCT, respectively. The L-shell parallel collimator AuNP imaging sensitivity approached that of the pencil beam geometry with a 55-fold reduction in imaging time. The AuNP sensitivity limits for the 1 mm diameter vial for the L-shell and K-shell parallel collimator XFCT images were 0.19% and 0.16%, respectively, and those of the pencil beam XFCT images were 0.08% and 0.01% for L-shell and K-shell XFCT, respectively. The remaining two collimator geometries resulted in a lower CNR and poorer image quality. For a 2 mm diameter vial, the AuNP sensitivity limits for the L-shell and K-shell multipinhole collimator XFCT images were 0.23% and 0.52%, respectively, while for the L-shell and K-shell converging collimator XFCT images the AuNP sensitivity limits were 0.38% and 0.13%, respectively.

CONCLUSION

This work demonstrates the feasibility of sheet beam L-shell XFCT imaging for small animal studies using parallel-oriented lead collimators which can detect AuNP concentrations approaching the level of pencil beam images with reduced imaging time.

Brief Report: External Beam Radiation Therapy for the Treatment of Human Pluripotent Stem Cell-Derived Teratomas

Human pluripotent stem cells, including human embryonic stem cells (hESCs) and human induced PSCs (hiPSCs), have great potential as an unlimited donor source for cell-based therapeutics. The risk of teratoma formation from residual undifferentiated cells, however, remains a critical barrier to the clinical application of these cells. Herein, we describe external beam radiation therapy (EBRT) as an attractive option for the treatment of this iatrogenic growth. We present evidence that EBRT is effective in arresting growth of hESC-derived teratomas in vivo at day 28 post-implantation by using a microCT irradiator capable of targeted treatment in small animals. Within several days of irradiation, teratomas derived from injection of undifferentiated hESCs and hiPSCs demonstrated complete growth arrest lasting several months. In addition, EBRT reduced reseeding potential of teratoma cells during serial transplantation experiments, requiring irradiated teratomas to be seeded at 1 × 10³ higher doses to form new teratomas. We demonstrate that irradiation induces teratoma cell apoptosis, senescence, and growth arrest, similar to established radiobiology mechanisms. Taken together, these results provide proof of concept for the use of EBRT in the treatment of existing teratomas and highlight a strategy to increase the safety of stem cell-based therapies. Stem Cells 2017;35:1994-2000.

Multi-institutional MicroCT image comparison of image-guided small animal irradiators

To recommend imaging protocols and establish tolerance levels for microCT image quality assurance (QA) performed on conformal image-guided small animal irradiators. A fully automated QA software SAPA (small animal phantom analyzer) for image analysis of the commercial Shelley micro-CT MCTP 610 phantom was developed, in which quantitative analyses of CT number linearity, signal-to-noise ratio (SNR), uniformity and noise, geometric accuracy, spatial resolution by means of modulation transfer function (MTF), and CT contrast were performed. Phantom microCT scans from eleven institutions acquired with four image-guided small animal irradiator units (including the commercial PXi X-RAD SmART and Xstrahl SARRP systems) with varying parameters used for routine small animal imaging were analyzed. Multi-institutional data sets were compared using SAPA, based on which tolerance levels for each QA test were established and imaging protocols for QA were recommended. By analyzing microCT data from 11 institutions, we established image QA tolerance levels for all image quality tests. CT number linearity set to R ²  >  0.990 was acceptable in microCT data acquired at all but three institutions. Acceptable SNR  >  36 and noise levels  <55 HU were obtained at five of the eleven institutions, where failing scans were acquired with current-exposure time of less than 120 mAs. Acceptable spatial resolution (>1.5 lp mm^-1 for MTF  =  0.2) was obtained at all but four institutions due to their large image voxel size used (>0.275 mm). Ten of the eleven institutions passed the set QA tolerance for geometric accuracy (<1.5%) and nine of the eleven institutions passed the QA tolerance for contrast (>2000 HU for 30 mgI ml^-1). We recommend performing imaging QA with 70 kVp, 1.5 mA, 120 s imaging time, 0.20 mm voxel size, and a frame rate of 5 fps for the PXi X-RAD SmART. For the Xstrahl SARRP, we recommend using 60 kVp, 1.0 mA, 240 s imaging time, 0.20 mm voxel size, and 6 fps. These imaging protocols should result in high quality images that pass the set tolerance levels on all systems. Average SAPA computation time for complete QA analysis for a 0.20 mm voxel, 400 slice Shelley phantom microCT data set was less than 20 s. We present image quality assurance recommendations for image-guided small animal radiotherapy systems that can aid researchers in maintaining high image quality, allowing for spatially precise conformal dose delivery to small animals.