Proton Range Research Articles

Technical Note: Monte Carlo methods to comprehensively evaluate the robustness of 4D treatments in proton therapy.

Proton therapy is very sensitive to treatment uncertainties. These uncertainties can induce proton range variations and may lead to severe dose distortions. However, most commercial tools only offer a limited integration of these uncertainties during treatment planning. In order to verify the robustness of a treatment plan, this study aims at developing a comprehensive Monte Carlo simulation of the treatment delivery, including the simulation of setup and range errors, variation of the breathing motion, and interplay effect. Most clinically relevant uncertainties have been modeled and implemented in the fast Monte Carlo dose engine MCsquare. Especially, variation of the breathing motion is taken into account by deforming the initial Four-dimensional computed tomography (4DCT) series and generating multiple new 4DCT series with scaled motion. Systematic and random errors are randomly sampled, following a Monte Carlo approach, to generate individual erroneous treatment scenarios. The robustness of treatment plans is analyzed and reported with dose-volume histogram(DVH) bands. The statistical uncertainty coming from the Monte Carlo scenario sampling is studied. A validation demonstrated the ability of the motion model to generate new 4DCT series with scaled motion amplitude and improved image quality in comparison to the initial 4DCT. The robustness analysis is applied to a lung tumor treatment. Considering the proposed uncertainty model, the simulation of 300 treatment scenarios was necessary to reach an acceptable level of statistical uncertainty on the DVH band. A comprehensive and statistically sound method of treatment plan robustness verification is proposed. The uncertainty model presented in this paper is not specific to protons and can also be applied to photon treatments. Moreover, the generated 4DCT series, with scaled motion, can be imported in commercial TPSs.

First Cerenkov charge-induction (CCI) TlBr detector for TOF-PET and proton range verification

Thallium bromide (TlBr) is a semiconductor material and, simultaneously, a good Cerenkov radiator. The performance of a TlBr detector that integrates two different readouts, the charge induction readout and the detection of Cerenkov light, was evaluated. A TlBr detector with dimensions of 4 × 4 × 5 mm3, with a monolithic cathode and an anode segmented into strips, was manufactured. One of the bare and polished 4 × 4 mm2 faces of the detector was coupled to a silicon photomultiplier (SiPM) to read out the Cerenkov light. Simultaneous timing and energy resolutions of <400 ps full width at half maximum (FWHM) and ~8.5% at 511 keV were measured using the Cerenkov detection and charge induction readouts, respectively. A coincidence time resolution of 330 ps was obtained when selecting Cerenkov events with amplitudes above 70 mV. The combination of both readouts showed the potential to resolve the depth-of-interaction (DOI) positioning, based on the improvement of energy resolution when selecting events with similar electron drift times. This manuscript sets the stage for a new family of semiconductor detectors that combine charge induction readout with the Cerenkov light detection. Such detectors can provide, simultaneously, outstanding timing, energy, and spatial resolution, and will be an excellent fit for applications that require the detection of high-energy gamma photons with high timing accuracy, such as time-of-flight positron emission tomography (TOF-PET) and prompt gamma imaging (PGI) to assess the particle range in hadron therapy.

Compact Method for Proton Range Verification Based on Coaxial Prompt Gamma-Ray Monitoring: a Theoretical Study.

Range uncertainties in proton therapy hamper treatment precision. Prompt gamma-rays were suggested 16 years ago for real-time range verification, and have already shown promising results in clinical studies with collimated cameras. Simultaneously, alternative imaging concepts without collimation are investigated to reduce the footprint and price of current prototypes. In this manuscript, a compact range verification method is presented. It monitors prompt gamma-rays with a single scintillation detector positioned coaxially to the beam and behind the patient. Thanks to the solid angle effect, proton range deviations can be derived from changes in the number of gamma-rays detected per proton, provided that the number of incident protons is well known. A theoretical background is formulated and the requirements for a future proof-of-principle experiment are identified. The potential benefits and disadvantages of the method are discussed, and the prospects and potential obstacles for its use during patient treatments are assessed. The final milestone is to monitor proton range differences in clinical cases with a statistical precision of 1 mm, a material cost of 25000 USD and a weight below 10 kg. This technique could facilitate the widespread application of in vivo range verification in proton therapy and eventually the improvement of treatment quality.

Statistical evaluation of worst-case robust optimization intensity-modulated proton therapy plans using an exhaustive sampling approach

PurposeTo assess the worst-case robust optimization IMPT plans with setup and range uncertainties and to test the hypothesis that the worst-case robust optimization strategies could cover most possible setup and range uncertainties in the real scenarios.MethodsWe analyzed the nominal and worst-case robust optimization IMPT plans of seven patients with head and neck cancer patients. To take uncertainties into account for the dose calculation, we performed a comprehensive simulation in which the dose was recalculated 625 times per given plan using Gaussian systematic setup and proton range uncertainties. Subsequently, based on the simulation results, we calculated the target coverage in all perturbation scenarios, as well as the ratios of target coverage located within the threshold of eight worst-case scenarios. We set the criteria for the optimized plan to be the ratios of 1) the dose delivered to 95% (D95%) of clinical target volumes 1 and 2 (CTV1 and CTV2) above 95% of the prescribed dose, and 2) the D95% of clinical target volume 3 (CTV3) above 90% of the prescribed dose in worst-case situations.ResultsThe probability that the perturbed-dose indices of the CTVs in each scenario were within the worst-case scenario limits ranged from 89.51 to 91.22% for both the nominal and worst-case robust optimization IMPT plans. A quartile analysis showed that the selective robust optimization IMPT plans all had higher D95% values for CTV1, CTV2, and CTV3 than did the nominal IMPT plans.ConclusionsThe worst-case strategy for robust optimization is adequately models and covers most of the setup and range uncertainties for the IMPT treatment of head and neck patients in our center.

Towards proton arc therapy: physical and biologically equivalent doses with increasing number of beams in pediatric brain irradiation

Background: Proton arc therapy may improve physical dose conformity and reduce concerns of elevated linear energy transfer (LET) and relative biological effectiveness (RBE) at the end of the proton range, while offering more degrees of freedom for normal tissue sparing. To explore the potential of proton arc therapy, we studied the effect of increasing the number of beams on physical and biologically equivalent dose conformity in the setting of pediatric brain tumors.Material and methods: A cylindrical phantom (Ø = 150 mm) with central cylindrical targets (Ø = 25 and 30 mm) was planned with increasing number of equiangular coplanar proton beams (from 3 to 36). For four anonymized pediatric brain tumor patients, two ‘surrogate’ proton arc plans (18 equiangular coplanar or sagittal beams) and a reference plan with 3 non-coplanar beams were constructed. Biologically equivalent doses were calculated using two RBE scenarios: RBE1.1; and RBELET, the physical dose weighted by the LET. For both RBE scenarios, dose gradients were assessed, and doses to cognitive brain structures were reported.Results: Increasing the number of beams resulted in an improved dose gradient and reduced volume exposed to intermediate LET levels, at the expense of increased low-dose and low-LET volumes. Most of the differences between the two RBE scenarios were seen around the prescription dose level, where the isodose volumes increased with the RBELET plans, e.g. up to 63% in the 3-beam plan for the smallest phantom target. Overall, the temporal lobes were better spared with the sagittal proton arc surrogate plans, e.g. a mean dose of 3.9 Gy compared to 6 Gy in the reference 3-beam plan (median value, RBE1.1).Conclusion: Proton arc therapy has the potential to improve dose gradients to better spare cognitive brain structures. However, this is at the expense of increased low-dose/low-LET volumes, with possible implications for secondary cancer risks.

Evaluating the benefit of PBS vs. VMAT dose distributions in terms of dosimetric sparing and robustness against inter-fraction anatomical changes for pediatric abdominal tumors

Background and purposeTo evaluate the dosimetric sparing and robustness against inter-fraction anatomical changes between photon and proton dose distributions for children with abdominal tumors. Material and methodsVolumetric modulated arc therapy (VMAT) and intensity-modulated pencil beam scanning (PBS) proton dose distributions were calculated for 20 abdominal pediatric cases (average 3, range 1–8 years). VMAT plans were based on a full-arc while PBS plans on 2–3 posterior-oblique irradiation fields. Plans were robustly optimized on a patient-specific internal target volume (ITV) using a uniform 5 mm set-up uncertainty. Additionally, for the PBS plans a ± 3% proton range uncertainty was accounted for. Fractional dose re-calculations were performed using the planning computed tomography (CT) deformably registered to the daily cone-beam CT (CBCT) images. Fractional doses were accumulated rigidly. Planned and CBCT accumulated VMAT and PBS dose distributions were compared using dose–volume histogram (DVH) parameters. ResultsSignificant better sparing of the organs at risk with a maximum reduction in the mean dose of 40% was achieved with PBS. Mean ITV DVH parameters differences between planned and CBCT accumulated dose distributions were smaller than 0.5% for both VMAT and PBS. However, the ITV coverage (V95% > 99%) was not reached for one patient for the accumulated VMAT dose distribution. ConclusionsFor pediatric patients with abdominal tumors, improved dosimetric sparing was obtained with PBS compared to VMAT. In addition, PBS delivered by posterior-oblique irradiation fields demonstrated to be robust against anatomical inter-fraction changes. Compared to PBS, daily anatomical changes proved to affect the target coverage of VMAT dose distributions to a higher extent.

Experimental evaluation of depth-encoding absorber designs for prompt-gamma Compton imaging in proton therapy

Prompt-gamma Compton imaging (PGCI) has been presented as a promising in-vivo method for proton range verification. An accurate estimation of the Bragg-peak position can potentially be achieved by imaging high-energy prompt-gamma rays (in the range of several MeV). Therefore, scintillation detectors with thick and high Z crystals are mostly used as Compton absorbers for high-energy gamma-ray detection. However, an absorber using thick crystals degrades the angular resolution unless the absorber can provide continuous depth-of-interaction (DOI) measurement. The study investigated various detector configurations using thick crystals for developing a Compton absorber with DOI resolving capability without compromising energy performance. Two groups of lutetium yttrium oxyorthosilicate (LYSO) arrays with dimensions of 1.8 × 1.8 × 50 mm3 and 1.8 × 1.8 × 20 mm3, and each with four different surface treatments (combining crystal surface finishing and type of reflector coverage), were constructed for the study. The DOI detector utilized the dual-ended readout of pixelated scintillator arrays for depth encoding. The results revealed that the influences of type of the reflector coverage and crystal surface roughness on the performance of the DOI detectors for the 50- and 20-mm-thick LYSO differed greatly. Moreover, the combined use of the proposed partial inter-crystal reflector coverage and the unpolished surface finishing on the 50-mm-thick crystals were shown to improve both the DOI resolution and the energy performance but degraded the flood map quality. These results provide useful guidance for selecting an applicable Compton DOI absorber design for developing a high-performance PGCI system.

Design and commissioning of an image-guided small animal radiation platform and quality assurance protocol for integrated proton and x-ray radiobiology research

Small animal x-ray irradiation platforms are expanding the capabilities and future pathways for radiobiology research. Meanwhile, proton radiotherapy is transitioning to a standard treatment modality in the clinician’s precision radiotherapy toolbox, highlighting a gap between state-of-the-art clinical radiotherapy and small animal radiobiology research. Comparative research of the biological differences between proton and x-ray beams could benefit from an integrated small animal irradiation system for in vivo experiments and corresponding quality assurance (QA) protocols to ensure rigor and reproducibility. The objective of this study is to incorporate a proton beam into a small animal radiotherapy platform while implementing QA modelled after clinical protocols.A 225 kV x-ray small animal radiation research platform (SARRP) was installed on rails to align with a modified proton experimental beamline from a 230 MeV cyclotron-based clinical system. Collimated spread out Bragg peaks (SOBP) were produced with beam parameters compatible with small animal irradiation. Proton beam characteristics were measured and alignment reproducibility with the x-ray system isocenter was evaluated. A QA protocol was designed to ensure consistent proton beam quality and alignment. As a preliminary study, cellular damage via γ-H2AX immunofluorescence staining in an irradiated mouse tumor model was used to verify the beam range in vivo.The beam line was commissioned to deliver Bragg peaks with range 4–30 mm in water at 2 Gy min−1. SOBPs were delivered with width up to 25 mm. Proton beam alignment with the x-ray system agreed within 0.5 mm. A QA phantom was created to ensure reproducible alignment of the platform and verify beam delivery. γ-H2AX staining verified expected proton range in vivo.An image-guided small animal proton/x-ray research system was developed to enable in vivo investigations of radiobiological effects of proton beams, comparative studies between proton and x-ray beams, and investigations into novel proton treatment methods.

Fast In Situ Image Reconstruction for Proton Radiography.

Proton beam therapy is an emerging modality for cancer treatment that, compared to X-ray radiation therapy, promises to provide better dose delivery to clinical targets with lower doses to normal tissues. Crucial to accurate treatment planning and dose delivery is knowledge of the water equivalent path length (WEPL) of each ray, or pencil beam, from the skin to every point in the target. For protons, this length is estimated from relative stopping power based on X-ray Hounsfield units. Unfortunately, such estimates lead to 3 to 4% uncertainties in the proton range prediction. Therefore, protons in the Bragg peak may overshoot (or undershoot) the desired stopping depth in the target causing tissue damage beyond the target volume. Recent studies indicate that tomographic imaging using protons has the potential to provide directly more accurate measurement of RSPs with significantly lower radiation dose than X-rays. We are currently working on a proton radiography system that promises to provide accurate two-dimensional (2D) images of WEPL values for protons that pass through the body. These will be suitable for positioning and range verification in daily treatments. In this study, we demonstrate that this system is capable of rapidly achieving such accurate images in clinically meaningful times. We have developed a software platform to characterize the potential performance of the prototype proton radiography system. We use Geant4 to simulate raw data detected by the device. An especially-written software - pRad - was written to process these data as they are received and uses iterative methods to generate radiographs. The software has been designed to generate a radiograph from a few million protons in under a minute after receiving the first proton from the device. We used a head phantom with known chemical compositions that could be modelled quite accurately in Geant4 simulations of proton radiographs. The radiographs are displayed as pixelated WEPL values displayed on a 2D gray scale image of WEPL values. Rapid radiograph reconstruction of 3D phantoms using simulated proton pencil beams have been achieved with our software platform. On a modest desktop computer with a single central processing unit (CPU) and a single graphics processing unit (GPU), it takes about 11 seconds to reconstruct images using iterative linear algorithms to reconstruct a radiograph from 7.6 million protons. For the radiographic reconstructions of the head phantom described here, the mean WEPL errors, in the proton radiograph using a large majority of the pixels in the complete image were less than 1 mm when compared to images obtained without proton scattering and without detector resolution included. We have demonstrated, through computer simulations of proton irradiation of a pediatric head phantom using the newly built pRad detector and image reconstruction software, that high quality proton radiographs can be generated for patient alignment and verification of water equivalent thickness of the patient before each treatment.

Processing of prompt gamma-ray timing data for proton range measurements at a clinical beam delivery

In proton therapy, patients benefit from the precise deposition of the dose in the tumor volume due to the interaction of charged particles with matter. Currently, the determination of the beam range in the patient’s body during the treatment is not a clinical standard. This lack causes broad safety margins around the tumor, which limits the potential of proton therapy. To overcome this obstacle, different methods are under investigation aiming at the verification of the proton range in real time during the irradiation. One approach is the prompt gamma-ray timing (PGT) method, where the range of the primary protons is derived from time-resolved profiles (PGT spectra) of promptly emitted gamma rays, which are produced along the particle track in tissue. After verifying this novel technique in an experimental environment but far away from treatment conditions, the translation of PGT into clinical practice is intended. Therefore, new hardware was extensively tested and characterized using short irradiation times of 70 ms and clinical beam currents of 2 nA. Experiments were carried out in the treatment room of the University Proton Therapy Dresden. A pencil beam scanning plan was delivered to a target without and with cylindrical air cavities of down to 5 mm thickness. The range shifts of the proton beam induced due to the material variation could be identified from the corresponding PGT spectra, comprising events collected during the delivery of a whole energy layer. Additionally, an assignment of the PGT data to the individual pencil beam spots allowed a spot-wise analysis of the variation of the PGT distribution mean and width, corresponding to range shifts produced by the different air cavities. Furthermore, the paper presents a comprehensive software framework which standardizes future PGT analysis methods and correction algorithms for technical limitations that have been encountered in the presented experiments.

OC-0187 Comparison of proton range predictions between Single- and Dual-Energy CT using prompt gamma imaging

OC-0187 Comparison of proton range predictions between Single- and Dual-Energy CT using prompt gamma imaging

An analytical approximation of proton Bragg curves in lithium fluoride for beam energy distribution analysis

An analytical model that approximates Bragg curves due to protons in water is extended to lithium fluoride (LiF), both in bulk and lower-density homogeneous thin film form. To that purpose, Monte Carlo simulations of proton range for energies from 5 to 65 MeV are best fitted so that optimal values of the model parameters are deduced for LiF. Tests are performed by comparing analytically computed Bragg curves with corresponding ones simulated by Monte Carlo software in LiF, bulk and lower density, at two example proton energies and energy-straggling widths. The model is then applied to the reconstruction of the energy spectra of proton beams by best fitting experimental Bragg profiles due to visible photoluminescence of optically-active F2 and F3+ color centers created in LiF crystals by the ionizing interaction of protons with the material.

Iterative optimization of relative stopping power by single detector based multi-projection proton radiography

In proton therapy, range uncertainties induced by the conversion from x-ray CT (xCT) Hounsfield units (HU) to relative stopping power (RSP) compromise the precision of dose delivery. To reduce range uncertainties induced by HU-converted RSPs, this study investigates optimizing the RSP of individual voxels in xCT iteratively based on multi-projection proton radiography (pRG) acquired using a single amorphous silicon flat panel imager. Time-resolved dose rate functions (DRF) were measured by the imager placed downstream of a test phantom consisting of tissue substitute materials. Water equivalent path lengths (WEPL) in the pRG were derived. By rotating the phantom, multiple pRG projections were acquired at angles from 0 to 358° with an increment of 2°. X-ray CT of the phantom was acquired and co-registered with the pRG acquisition coordinates. RSPs of individual xCT voxels were optimized iteratively by minimizing the difference between the measured WEPLs and the calculated WEPLs by ray tracing with HU-converted RSPs. Pixels in pRGs that exhibited severe proton range mixing were rejected for the optimization. Tikhonov regularization was applied under the assumption that HU-converted RSPs are inaccurate, but the inaccuracy is within a few percent. While ~50% of WEPL pixels were rejected due to severe range mixing in pRG, RSPs of >90% CT voxels could still be optimized if multiple pRG projections, e.g. ⩾12, around the phantom are utilized. For tissue substitute materials in a cylindrical phantom, percentage errors of RSPs were reduced from a range of −8% to +4% to be within ±2%. Further optimization, achieved by implementing a material-specific regularization parameter, reduced percent errors to be within ±0.5%. This study demonstrates the concept of optimizing RSPs of individual CT voxels with multi-projection pRGs acquired by a single flat panel imager, which could be further explored and implemented in proton therapy to reduce range uncertainties.

A Monte Carlo feasibility study for neutron based real-time range verification in proton therapy

Uncertainties in the proton range in tissue during proton therapy limit the precision in treatment delivery. These uncertainties result in expanded treatment margins, thereby increasing radiation dose to healthy tissue. Real-time range verification techniques aim to reduce these uncertainties in order to take full advantage of the finite range of the primary protons. In this paper, we propose a novel concept for real-time range verification based on detection of secondary neutrons produced in nuclear interactions during proton therapy. The proposed detector concept is simple; consisting of a hydrogen-rich converter material followed by two charged particle tracking detectors, mimicking a proton recoil telescopic arrangement. Neutrons incident on the converter material are converted into protons through elastic and inelastic (n,p) interactions. The protons are subsequently detected in the tracking detectors. The information on the direction and position of these protons is then utilized in a new reconstruction algorithm to estimate the depth distribution of neutron production by the proton beam, which in turn is correlated with the primary proton range. In this paper, we present the results of a Monte Carlo feasibility study and show that the proposed concept could be used for real-time range verification with millimetric precision in proton therapy.

Compton-based prompt gamma imaging using ordered origin ensemble algorithm with resolution recovery in proton therapy

Prompt gamma ray (PG) imaging based on Compton camera (CC) is promising to realize in vivo verification during the proton therapy. However, the finite spatial and energy resolution of current CC, as well as the Doppler broaden effect, degrade the quality and resolution of PG images. In addition, due to the inherent geometrical complexity of Compton camera data, PG imaging can be time-consuming and difficult to reconstruct in real-time, while using standard techniques such as filtered back-projection or maximum likelihood-expectation maximization. In this paper, we propose three modifications of origin ensembles with resolution recovery (OE-RR) algorithm based on Markov chains to accelerate the convergence to equilibrium of OE-RR algorithm and improve the image quality. For evaluation, we performed a Monte Carlo simulation of a three-stage CZT Compton camera with resolution loss to detect the PG produced by a proton beam in a water phantom, and evaluate image quality of the gamma rays emitted during proton irradiation. The results show that our ordered OE-RR algorithm realized a good resolution recovery and accurate estimation of the position, including the peak and the distal falloff of the PG emission with remarkably faster reconstruction, thus demonstrating the feasibility of this new method in non-idealized PG-based proton range verification.

Measurement of the Beam Energy Distribution of a Medical Cyclotron with a Multi-Leaf Faraday Cup

Accurate knowledge of the beam energy distribution is crucial for particle accelerators, compact medical cyclotrons for the production of radioisotopes in particular. For this purpose, a compact instrument was developed, based on a multi-leaf Faraday cup made of thin aluminum foils interleaved with plastic absorbers. The protons stopping in the aluminum foils produce a measurable current that is used to determine the range distribution of the proton beam. On the basis of the proton range distribution, the beam energy distribution is assessed by means of stopping-power Monte Carlo simulations. In this paper, we report on the design, construction, and testing of this apparatus, as well as on the first measurements performed with the IBA Cyclone 18-MeV medical cyclotron in operation at the Bern University Hospital.

Robust optimization to reduce the impact of biological effect variation from physical uncertainties in intensity-modulated proton therapy

Robust optimization (RO) methods are applied to intensity-modulated proton therapy (IMPT) treatment plans to ensure their robustness in the face of treatment delivery uncertainties, such as proton range and patient setup errors. However, the impact of those uncertainties on the biological effect of protons has not been specifically considered. In this study, we added biological effect-based objectives into a conventional RO cost function for IMPT optimization to minimize the variation in biological effect.One brain tumor case, one prostate tumor case and one head & neck tumor case were selected for this study. Three plans were generated for each case using three different optimization approaches: planning target volume (PTV)-based optimization, conventional RO, and RO incorporating biological effect (BioRO). In BioRO, the variation in biological effect caused by IMPT delivery uncertainties was minimized for voxels in both target volumes and critical structures, in addition to a conventional voxel-based worst-case RO objective function. The biological effect was approximated by the product of dose-averaged linear energy transfer (LET) and physical dose. All plans were normalized to give the same target dose coverage, assuming a constant relative biological effectiveness (RBE) of 1.1. Dose, biological effect, and their uncertainties were evaluated and compared among the three optimization approaches for each patient case.Compared with PTV-based plans, RO plans achieved more robust target dose coverage and reduced biological effect hot spots in critical structures near the target. Moreover, with their sustained robust dose distributions, BioRO plans not only reduced variations in biological effect in target and normal tissues but also further reduced biological effect hot spots in critical structures compared with RO plans.Our findings indicate that IMPT could benefit from the use of conventional RO, which would reduce the biological effect in normal tissues and produce more robust dose distributions than those of PTV-based optimization. More importantly, this study provides a proof of concept that incorporating biological effect uncertainty gap into conventional RO would not only control the IMPT plan robustness in terms of physical dose and biological effect but also achieve further reduction of biological effect in normal tissues.

Modeling the Acoustic Field Generated by a Pulsed Beam for Experimental Proton Range Verification

Proton range verification by ionoacoustic wave sensing is a technique under development for applications in adron therapy as an alternative to nuclear imaging. It provides an acoustic imaging of the proton energy deposition vs. depth using the acoustic wave Time of Flight (ToF). State-of-the-art (based on simulations and experimental results) points out that this detection technique achieves better spatial resolution (&lt; 1 mm) of the proton range comparing with Positron-Emission-Tomography (PET) and prompt gamma ray techniques. This work presents a complete Geant4/k-Wave model that allows to understand several physical phenomena and to evaluate the key parameters that affect the acoustic field generated by the incident proton radiation.

A single detector energy-resolved proton radiography system: a proof of principle study by Monte Carlo simulations

Taking advantage of Bragg peak and small spot size, pencil beam scanning proton therapy can deliver a highly conformal dose distribution to target while sparing normal tissues. However, such dose distributions can be highly sensitive to the proton range uncertainty which can reach 5% or higher in lung tissue. One proposed method for reducing range uncertainty is to measure the water equivalent path length (WEPL) by proton radiography. In this study, we followed a newly proposed proton beam radiography technique based on energy resolved dose functions (ERDF) to construct a Monte Carlo model for a single detector energy-resolved proton radiography system (SDPRS). This SDPRS model was constructed in the Monte Carlo software package TOPAS (TOol for PArticle Simulation) and it includes the Mevion HYPERSCAN™ pencil beam scanning treatment head and a 2D dose detector positioned downstream as the imager. A calibration phantom containing a number of tissue equivalent materials was simulated to evaluate the accuracy in WEPL measurement by SDPRS. The mean deviation of the obtained relative stopping power (RSP) from the reference values was 0.31%. Proton radiographs of an anthropomorphic head phantom were also generated to demonstrate the clinical relevance of the technique. Effects of different energy layer spacing and measurement noise were also studied.

A quality assurance for respiratory gated proton irradiation with range modulation wheel.

The purpose of this study was to provide periodic quality assurance (QA) methods for respiratory‐gated proton beam with a range modulation wheel (RMW) and to clarify the characteristics and long‐term stability of the respiratory‐gated proton beam. A two‐dimensional detector array and a solid water phantom were used to measure absolute dose, spread‐out Bragg peak (SOBP) width and proton range for monthly QA. SOBP width and proton range were measured using an oblique incidence beam to the lateral side of a solid water phantom and compared between with and without a gating proton beam. To measure the delay time of beam‐on/off for annual QA, we collected the beam‐on/off signals and the dose monitor‐detected pulse. We analyzed the results of monthly QA over a 15‐month period and investigated the delay time by machine signal analysis. The dose deviations at proximal, SOBP center and distal points were −0.083 ± 0.25%, 0.026 ± 0.20%, and −0.083 ± 0.35%, respectively. The maximum dose deviation between with and without respiratory gating was −0.95% at the distal point and other deviations were within ±0.5%. Proximal and SOBP center doses showed the same trend over a 15‐month period. Delay times of beam‐on/off for 200 MeV/SOBP 16 cm were 140.5 ± 0.8 ms and 22.3 ± 13.0 ms, respectively. Delay times for 160 MeV/SOBP 10 cm were 167.5 ± 15.1 ms and 19.1 ± 9.8 ms. Our beam delivery system with the RMW showed sufficient stability for respiratory‐gated proton therapy and the system did not show dependency on the energy and the respiratory wave form. The delay times of beam‐on/off were within expectations. The proposed QA methods will be useful for managing the quality of respiratory‐gated proton beams and other beam delivery systems.