skip to main content

Explore Research Products in the NSF-PAR

Title: Design and validation of a synchrotron proton beam line for FLASH radiotherapy preclinical research experiments
Abstract Purpose

The main purpose of this work was to generate and validate the dosimetric accuracy of proton beams of dimensions that are appropriate for in vivo small animal and in vitro ultrahigh dose rate (FLASH) radiotherapy experiments using a synchrotron‐based treatment delivery system. This study was performed to enable future investigations of the relevance of a spread‐out Bragg peak (SOBP) under FLASH conditions.

 Methods

The spill characteristics of the small field fixed horizontal beam line were modified to deliver accelerated protons in times as short as 2 ms and to control the dose delivered. A Gaussian‐like transverse beam profile was transformed into a square uniform one at FLASH dose rates, while avoiding low‐dose regions, a crucial requirement to protect normal tissue during FLASH irradiation. Novel beam‐shaping devices were designed using Monte Carlo techniques to produce up to about 6 cm3of uniform dose in SOBPs while maximizing the dose rate. These included a scattering foil, a conical flattening filter to maximize the flux of protons into the region of interest, energy filters, range compensators, and collimators. The shapes, sizes, and positions of the components were varied to provide the required field sizes and SOBPs.

 Results

The designed and fabricated devices were used to produce 10‐, 15‐, and 20‐mm diameter, circular field sizes and 10‐, 15‐, and 20‐mm SOBP modulation widths at uniform physical dose rates of up to 375 Gy/s at the center of the SOBP and a minimum dose rate of about 255 Gy/s at the entrance, respectively, in cylindrical volumes. The flatness of lateral dose profiles at the center could be adjusted to within ±1.5% at the center of the SOBP. Assessment of systematic uncertainties, such as impact of misalignments and positioning uncertainties, was performed using simulations, and the results were used to provide appropriate adjustments to ensure high‐accuracy FLASH beam delivery for both in vitro and in vivo preclinical experiments.

 Conclusions

It is feasible to use synchrotron‐generated proton beams of sufficient dimensions for FLASH radiobiology experiments. We expect to use the system we developed to acquire in vitro and in vivo small animal FLASH radiobiology data as a function of dose, dose rate, oxygen content, and linear energy transfer to help us understand the underlying mechanisms of the FLASH phenomenon.

  more » « less
NSF-PAR ID:
10446137
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Medical Physics
Volume:
49
Issue:
1
ISSN:
0094-2405
Page Range / eLocation ID:
p. 497-509
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ( , Medical Physics)
    Abstract Purpose

    The radiobiological benefits afforded by spatially fractionated (GRID) radiation therapy pairs well with the dosimetric advantages of proton therapy. Inspired by the emergence of energy‐layer specific collimators in pencil beam scanning (PBS), this work investigates how the spot spacing and collimation can be optimized to maximize the therapeutic gains of a GRID treatment while demonstrating the integration of a dynamic collimation system (DCS) within a commercial beamline to deliver GRID treatments and experimentally benchmark Monte Carlo calculation methods.

     Methods

    GRID profiles were experimentally benchmarked using a clinical DCS prototype that was mounted to the nozzle of the IBA‐dedicated nozzle system. Integral depth dose (IDD) curves and lateral profiles were measured for uncollimated and GRID‐collimated beamlets. A library of collimated GRID dose distributions were simulated by placing beamlets within a specified uniform grid and weighting the beamlets to achieve a volume‐averaged tumor cell survival equivalent to an open field delivery. The healthy tissue sparing afforded by the GRID distribution was then estimated across a range of spot spacings and collimation widths, which were later optimized based on the radiosensitivity of the tumor cell line and the nominal spot size of the PBS system. This was accomplished by using validated models of the IBA universal and dedicated nozzles.

     Results

    Excellent agreement was observed between the measured and simulated profiles. The IDDs matched above 98.7% when analyzed using a 1%/1‐mm gamma criterion with some minor deviation observed near the Bragg peak for higher beamlet energies. Lateral profile distributions predicted using Monte Carlo methods agreed well with the measured profiles; a gamma passing rate of 95% or higher was observed for all in‐depth profiles examined using a 3%/2‐mm criteria. Additional collimation was shown to improve PBS GRID treatments by sharpening the lateral penumbra of the beamlets but creates a trade‐off between enhancing the valley‐to‐peak ratio of the GRID delivery and the dose‐volume effect. The optimal collimation width and spot spacing changed as a function of the tumor cell radiosensitivity, dose, and spot size. In general, a spot spacing below 2.0 cm with a collimation less than 1.0 cm provided a superior dose distribution among the specific cases studied.

     Conclusions

    The ability to customize a GRID dose distribution using different collimation sizes and spot spacings is a useful advantage, especially to maximize the overall therapeutic benefit. In this regard, the capabilities of the DCS, and perhaps alternative dynamic collimators, can be used to enhance GRID treatments. Physical dose models calculated using Monte Carlo methods were experimentally benchmarked in water and were found to accurately predict the respective dose distributions of uncollimated and DCS‐collimated GRID profiles.

     
    more » « less
  2. ; ; ( , Medical Physics)
    Abstract Purpose

    One standard method, proton resonance frequency shift, for measuring temperature using magnetic resonance imaging (MRI), in MRI‐guided surgeries, fails completely below the freezing point of water. Because of this, we have developed a new methodology for monitoring temperature with MRI below freezing. The purpose of this paper is to show that a strong temperature dependence of the nuclear relaxation timeT1in soft silicone polymers can lead to temperature‐dependent changes of MRI intensity acquired withT1weighting. We propose the use of silicone filaments inserted in tissue for measuring temperature during MRI‐guided cryoablations.

     Methods

    The temperature dependence ofT1in bio‐compatible soft silicone polymers was measured using nuclear magnetic resonance spectroscopy and MRI. Phantoms, made of bulk silicone materials and put in an MRI‐compatible thermal container with dry ice, allowed temperature measurements ranging from –60°C to + 20°C.T1‐weighted gradient echo images of the phantoms were acquired at spatially uniform temperatures and with a gradient in temperature to determine the efficacy of using these materials as temperature indicators in MRI. Ex vivo experiments on silicone rods, 4 mm in diameter, inserted in animal tissue were conducted to assess the practical feasibility of the method.

     Results

    Measurements of nuclear relaxation times of protons in soft silicone polymers show a monotonic, nearly linear, change with temperature (R2 > 0.98) and have a significant correlation with temperature (Pearson'sr > 0.99,p < 0.01). Similarly, the intensity of the MR images in these materials, taken with a gradient echo sequence, are also temperature dependent. There is again a monotonic change in MRI intensity that correlates well with the measured temperature (Pearson'sr < ‐0.98 andp < 0.01). The MRI experiments show that a temperature change of 3°C can be resolved in a distance of about 2.5 mm. Based on MRI images and external sensor calibrations for a sample with a gradient in temperature, temperature maps with 3°C isotherms are created for a bulk phantom. Experiments demonstrate that these changes in MRI intensity with temperature can also be seen in 4 mm silicone rods embedded in ex vivo animal tissue.

     Conclusions

    We have developed a new method for measuring temperature in MRI that potentially could be used during MRI‐guided cryoablation operations, reducing both procedure time and cost, and making these surgeries safer.

     
    more » « less
  3. ; ; ; ; ; ; ( , Medical Physics)
    Abstract Purpose

    Most commercially available treatment planning systems (TPSs) approximate the continuous delivery of volumetric modulated arc therapy (VMAT) plans with a series of discretized static beams for treatment planning, which can make VMAT dose computation extremely inefficient. In this study, we developed a polar‐coordinate‐based pencil beam (PB) algorithm for efficient VMAT dose computation with high‐resolution gantry angle sampling that can improve the computational efficiency and reduce the dose discrepancy due to the angular under‐sampling effect.

     Methods and Materials

    6 MV pencil beams were simulated on a uniform cylindrical phantom under an EGSnrc Monte Carlo (MC) environment. The MC‐generated PB kernels were collected in the polar coordinate system for each bixel on a fluence map and subsequently fitted via a series of Gaussians. The fluence was calculated using a detectors’ eye view with off‐axis and MLC transmission factors corrected. Doses of VMAT arc on the phantom were computed by summing the convolution results between the corresponding PB kernels and fluence for each bixel in the polar coordinate system. The convolution was performed using fast Fourier transform to expedite the computing speed. The calculated doses were converted to the Cartesian coordinate system and compared with the reference dose computed by a collapsed cone convolution (CCC) algorithm of the TPS. A heterogeneous phantom was created to study the heterogeneity corrections using the proposed algorithm. Ten VMAT arcs were included to evaluate the algorithm performance. Gamma analysis and computation complexity theory were used to measure the dosimetric accuracy and computational efficiency, respectively.

     Results

    The dosimetric comparisons on the homogeneous phantom between the proposed PB algorithm and the CCC algorithm for 10 VMAT arcs demonstrate that the proposed algorithm can achieve a dosimetric accuracy comparable to that of the CCC algorithm with average gamma passing rates of 96% (2%/2mm) and 98% (3%/3mm). In addition, the proposed algorithm can provide better computational efficiency for VMAT dose computation using a PC equipped with a 4‐core processor, compared to the CCC algorithm utilizing a dual 10‐core server. Moreover, the computation complexity theory reveals that the proposed algorithm has a great advantage with regard to computational efficiency for VMAT dose computation on homogeneous medium, especially when a fine angular sampling rate is applied. This can support a reduction in dose errors from the angular under‐sampling effect by using a finer angular sampling rate, while still preserving a practical computing speed. For dose calculation on the heterogeneous phantom, the proposed algorithm with heterogeneity corrections can still offer a reasonable dosimetric accuracy with comparable computational efficiency to that of the CCC algorithm.

     Conclusions

    We proposed a novel polar‐coordinate‐based pencil beam algorithm for VMAT dose computation that enables a better computational efficiency while maintaining clinically acceptable dosimetric accuracy and reducing dose error caused by the angular under‐sampling effect. It also provides a flexible VMAT dose computation structure that allows adjustable sampling rates and direct dose computation in regions of interest, which makes the algorithm potentially useful for clinical applications such as independent dose verification for VMAT patient‐specific QA.

     
    more » « less
  4. ; ; ; ; ( , Alcoholism: Clinical and Experimental Research)
    Abstract Background

    Fetal alcohol syndrome (FAS) due to gestational alcohol exposure represents one of the most common causes of nonheritable lifelong disability worldwide. In vitro and in vivo models have successfully recapitulated multiple facets of the disorder, including morphological and behavioral deficits, but far less is understood regarding the molecular and genetic mechanisms underlying FAS.

     Methods

    In this study, we utilized an in vitro human pluripotent stem cell‐based (hPSC) model of corticogenesis to probe the effects of early, chronic intermittent alcohol exposure on the transcriptome of first trimester‐equivalent cortical neurons.

     Results

    We used RNA sequencing of developing hPSC‐derived neurons treated for 50 days with 50 mM ethanol and identified a relatively small number of biological pathways significantly altered by alcohol exposure. These included cell‐type specification, axon guidance, synaptic function, and regional patterning, with a notable upregulation of WNT signaling‐associated transcripts observed in alcohol‐exposed cultures relative to alcohol‐naïve controls. Importantly, this effect paralleled a shift in gene expression of transcripts associated with regional patterning, such that caudal forebrain‐related transcripts were upregulated at the expense of more anterior ones. Results from H9 embryonic stem cells were largely replicated in an induced pluripotent stem cell line (IMR90‐4), indicating that these patterning alterations are not cell line‐specific.

     Conclusions

    We found that a major effect of chronic intermittent alcohol on the developing cerebral cortex is an overall imbalance in regionalization, with enrichment of gene expression related to the production of posterodorsal progenitors and a diminution of anteroventral progenitors. This finding parallels behavioral and morphological phenotypes observed in animal models of high‐dose prenatal alcohol exposure, as well as patients with FAS.

     
    more » « less
  5. ; ; ; ; ; ( , Medical Physics)
    Abstract Purpose

    Pre‐calculation of accurate dose deposition kernels for treatment planning of spot‐based radiotherapies, such as Gamma Knife (GK) and Gamma Pod (GP), can be very time‐consuming and may require large data storage with an enormous number of possible spots. We proposed a novel kernel decomposition (KD) model to address accurate and fast (real‐time) dose calculation with reduced data storage requirements for spot‐based treatment planning. The application of the KD model was demonstrated for clinical GK and GP radiotherapy platforms.

     Methods

    The dose deposition kernel at each spot (shot position) is modeled as the product of a shift‐invariant kernel based on a reference kernel and spatially variant scale factor. The reference kernel, one for each collimator, is defined at the center of the commissioning phantom for GK and at the center of the treatment target for GP and calculated using the Monte Carlo (MC) method. The spatially variant scale factor is defined as the ratio of the mean tissue maximum ratio (TMR) at the candidate shot position to that at the reference kernel position, and the mean TMR map is calculated within the entire volume through parallel beam ray tracing on the density image followed by averaging over all source directions. The proposed KD dose calculations were compared with the MC method and with the GK and GP treatment planning system (TPS) computations for various shot positions and collimator sizes utilizing a phantom and 14 and 12 clinical plans for GK and GP, respectively.

     Results

    For the phantom study, the KD Gamma index (3%/1 mm) passing rates were greater than 99% (median 100%) relative to the MC doses, except for the shots close to the boundary. The passing rates dropped below 90% for 8 mm (16 mm) shots positioned within ∼1 cm (∼2 cm) of the boundary. For the clinical GK plans, the KD Gamma passing rates were greater than 99% (median 100%) compared to the MC and greater than 92% (median 99%) compared to the TPS. For the clinical GP plans, the KD Gamma passing rates were greater than 95% (median 98%) compared to the MC and greater than 91% (median 97%) compared to the TPS. The scale factors were calculated in sub‐seconds with GPU implementation and only need to be calculated once before treatment plan optimization. The calculation of the dose kernel was also within sub‐seconds without requiring beam‐by‐beam calculation commonly done in the TPS.

     Conclusion

    The proposed model can provide an accurate dose and enables real‐time dose and derivative calculations by kernel shifting and scaling without pre‐calculating or requiring large data storage for GK and GP dose deposition kernels during treatment planning. This model could be useful for spot‐based radiotherapy treatment planning by allowing an efficient global fine search for optimal spots.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email