skip to main content

Title: The dosimetric enhancement of GRID profiles using an external collimator in pencil beam scanning proton therapy
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.


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.


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.


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
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Medical Physics
Page Range / eLocation ID:
p. 2684-2698
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.


    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.


    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.


    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
  2. 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.


    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.


    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.


    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
  3. Purpose

    The combination of nonhuman primates (NHPs) with the state‐of‐the‐art molecular imaging technologies allows for within‐subject longitudinal research aiming at gaining new insights into human normal and disease conditions and provides an ideal foundation for future translational studies of new diagnostic tools, medical interventions, and therapies. However, radiation dose estimations for nonhuman primates from molecular imaging probes are lacking and are difficult to perform experimentally. The aim of this work is to construct age‐dependent NHP computational model series to estimate the absorbed dose to NHP specimens in common molecular imaging procedures.

    Materials and methods

    A series of NHP models from baby to adult were constructed based on nonuniform rational B‐spline surface (NURBS) representations. Particle transport was simulated using Monte Carlo calculations to estimate S‐values from nine positron‐emitting radionuclides and absorbed doses from PET radiotracers.


    Realistic age‐dependent NHP computational model series were developed. For most source‐target pairs in computational NHP models, differences between C‐11 S‐values were between −13.4% and −8.8%/kg difference in body weight while differences between F‐18 S‐values were between −12.9% and −8.0%/kg difference in body weight. The absorbed doses of11C‐labeled brain receptor substances,18F‐labeled brain receptor substances, and18F‐FDG in the brain ranged within 0.047–0.32 mGy/MBq, 0.25–1.63 mGy/MBq, and 0.32–2.12 mGy/MBq, respectively.


    The absorbed doses to organs are significantly higher in the baby NHP model than in the adult model. These results can be used in translational longitudinal studies to estimate the cumulated absorbed organ doses in NHPs at various ages.

    more » « less
  4. Abstract

    While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement.

    Significance Statement

    In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.

    more » « less
  5. Abstract

    Recent work demonstrates that high‐resolution global models forced simultaneously by atmospheric fields and the astronomical tidal potential contain a partial internal (gravity) wave (IW) spectral continuum. Regional simulations of the MITgcm forced at the horizontal boundaries by a global run that carries a partial IW continuum spectrum are performed at the same grid spacing as the global run and at finer grid spacings in an attempt to fill out more of the IW spectral continuum. Decreasing only the horizontal grid spacing from 2 to 0.25 km greatly improves the frequency spectra and slightly improves the vertical wavenumber spectra of the horizontal velocity. Decreasing only the vertical grid spacing by a factor of 3 does not yield any significant improvements. Decreasing both horizontal and vertical grid spacings yields the greatest degree of improvement, filling the frequency spectrum out to 72 cpd. Our results suggest that improved IW spectra in regional models are possible if they are run at finer grid spacings and are forced at their lateral boundaries by remotely generated IWs. Additionally, consistency relations demonstrate that improvements in the spectra are indeed due to the existence of IWs at higher frequencies and vertical wavenumbers when remote IW forcing is included and model grid spacings decrease. By being able to simulate an IW spectral continuum to 0.25 km scales, these simulations demonstrate that one may be able to track the energy pathways of IWs from generation to dissipation and improve the understanding of processes such as IW‐driven mixing.

    more » « less