Intra-fractional per-beam adaptive workflow to mitigate the need for a rotating gantry during MRI-guided proton therapy.

The integration of real-time magnetic resonance imaging (MRI) guidance and proton therapy would potentially improve the proton dose steering capability by reducing daily uncertainties due to anatomical variations. The use of a fixed beamline coupled with an axial patient couch rotation would greatly simplify the proton delivery with MRI guidance. Nonetheless, it is mandatory to assure that the plan quality is not deteriorated by the anatomical deformations due to patient rotation. In this work, an in-house tool allowing for intra-fractional per-beam adaptation of intensity-modulated proton plans (BeamAdapt) was implemented through features available in RayStation. A set of three MRIs was acquired for two healthy volunteers (V1,V2): (1) no rotation/static, (2) rotation to the right and (3) left.V1was rotated by 15°, to simulate a clinical pediatric abdominal case andV2by 45°, to simulate an extreme patient rotation case. For each volunteer, a total of four intensity-modulated pencil beam scanning plans were optimized on the static MRI using virtual abdominal targets and two-three posterior-oblique beams. Beam angles were defined according to the angulations on the rotated MRIs. With BeamAdapt, each original plan was initially converted into separate plans with one beam per plan. In an iterative order, individual beam doses were non-rigidly deformed to the rotated anatomies and re-optimized accounting for the consequent deformations and the beam doses delivered so far. For evaluation, the final accumulated dose distribution was propagated back to the static MRI. Planned and adapted dose distributions were compared by computing relative differences between dose-volume histogram metrics. Absolute target dose differences were on average below 1% and organs-at-risk mean dose differences were below 3%. With BeamAdapt, not only intra-fractional per-beam proton plan adaptation coupled with axial patient rotation is possible but also the need for a rotating gantry during MRI guidance might be mitigated.

A filtering approach for PET and PG predictions in a proton treatment planning system.

Positron emission tomography (PET) and prompt gamma (PG) detection are promising proton therapy monitoring modalities. Fast calculation of the expected distributions is desirable for comparison to measurements and to develop/train algorithms for automatic treatment error detection. A filtering formalism was used for positron-emitter predictions and adapted to allow for its use for the beamline of any proton therapy centre. A novel approach based on a filtering formalism was developed for the prediction of energy-resolved PG distributions for arbitrary tissues. The method estimates PG yields and their energy spectra in the entire treatment field. Both approaches were implemented in a research version of the RayStation treatment planning system. The method was validated against PET monitoring data and Monte Carlo simulations for four patients treated with scanned proton beams. Longitudinal shifts between profiles from analytical and Monte Carlo calculations were within -1.7 and 0.9 mm, with maximum standard deviation of 0.9 mm and 1.1 mm, for positron-emitters and PG shifts, respectively. Normalized mean absolute errors were within 1.2 and 5.3%. When comparing measured and predicted PET data, the same more complex case yielded an average shift of 3 mm, while all other cases were below absolute average shifts of 1.1 mm. Normalized mean absolute errors were below 7.2% for all cases. A novel solution to predict positron-emitter and PG distributions in a treatment planning system is proposed, enabling calculation times of only a few seconds to minutes for entire patient cases, which is suitable for integration in daily clinical routine.

Validation of the analytical irradiator model and Monte Carlo dose engine in the small animal irradiation treatment planning system µ-RayStation 8B.

Dose calculation in preclinical context with a clinical level of accuracy is a challenge due to the small animal scale and the medium photon energy range. In this work, we evaluate the effectiveness and accuracy of an analytical irradiator model combined with Monte Carlo (MC) calculations in the irradiated volume to calculate the dose delivered by a modern small animal irradiator. A model of the XRAD225Cx was created in µ-RayStation 8B, a preclinical treatment planning system, allowing arc and static beams for seven cylindrical collimators. Calculations with the µ-RayStation MC dose engine were compared with EBT3 measurements in water for all static beams and with a validated GATE model in water, heterogeneous media and a mouse CT. The GATE model is a complete MC representation of the XRAD225Cx. In water, µ-RayStation calculations, compared to GATE calculations and EBT3 measurements, agreed within a maximal error of 3.2% (mean absolute error of 0.6% and 0.8% respectively) and maximal distance-to-agreement (DTA) was 0.2 mm at 50% of the central dose. For a 5 mm static beam in heterogeneous media, the maximal absolute error between µ-RayStation and GATE calculations was below 1.3% in each medium and DTA was 0.1 mm at interfaces. For calculations on a mouse CT, µ-RayStation and GATE calculations agreed well for both static and arc beams. The 2D local gamma passing rate was >98.9% for 1%/0.3 mm criteria and >92.9% for 1%/0.2 mm criteria. Moreover, µ-RayStation reduces calculation time significantly comparing with GATE (speed-up factor between 120 and 680). These findings show that the analytical irradiator model presented in this work combined with the µ-RayStation MC dose engine accurately computes dose for the XRAD225Cx irradiator. The improvements in calculation time and availability of functionality and tools for managing, planning and evaluating the irradiation makes this platform very useful for pre-clinical irradiation research.

End-to-end tests using alanine dosimetry in scanned proton beams.

This paper describes end-to-end test procedures as the last fundamental step of medical commissioning before starting clinical operation of the MedAustron synchrotron-based pencil beam scanning (PBS) therapy facility with protons. One in-house homogeneous phantom and two anthropomorphic heterogeneous (head and pelvis) phantoms were used for end-to-end tests at MedAustron. The phantoms were equipped with alanine detectors, radiochromic films and ionization chambers. The correction for the 'quenching' effect of alanine pellets was implemented in the Monte Carlo platform of the evaluation version of RayStation TPS. During the end-to-end tests, the phantoms were transferred through the workflow like real patients to simulate the entire clinical workflow: immobilization, imaging, treatment planning and dose delivery. Different clinical scenarios of increasing complexity were simulated: delivery of a single beam, two oblique beams without and with range shifter. In addition to the dose comparison in the plastic phantoms the dose obtained from alanine pellet readings was compared with the dose determined with the Farmer ionization chamber in water. A consistent systematic deviation of about 2% was found between alanine dosimetry and the ionization chamber dosimetry in water and plastic materials. Acceptable agreement of planned and delivered doses was observed together with consistent and reproducible results of the end-to-end testing performed with different dosimetric techniques (alanine detectors, ionization chambers and EBT3 radiochromic films). The results confirmed the adequate implementation and integration of the new PBS technology at MedAustron. This work demonstrates that alanine pellets are suitable detectors for end-to-end tests in proton beam therapy and the developed procedures with customized anthropomorphic phantoms can be used to support implementation of PBS technology in clinical practice.

Experimental assessment of proton dose calculation accuracy in inhomogeneous media.

PURPOSE: Proton therapy with Pencil Beam Scanning (PBS) has the potential to improve radiotherapy treatments. Unfortunately, its promises are jeopardized by the sensitivity of the dose distributions to uncertainties, including dose calculation accuracy in inhomogeneous media. Monte Carlo dose engines (MC) are expected to handle heterogeneities better than analytical algorithms like the pencil-beam convolution algorithm (PBA). In this study, an experimental phantom has been devised to maximize the effect of heterogeneities and to quantify the capability of several dose engines (MC and PBA) to handle these. METHODS: An inhomogeneous phantom made of water surrounding a long insert of bone tissue substitute (1×10×10 cm3) was irradiated with a mono-energetic PBS field (10×10 cm2). A 2D ion chamber array (MatriXX, IBA Dosimetry GmbH) lied right behind the bone. The beam energy was such that the expected range of the protons exceeded the detector position in water and did not attain it in bone. The measurement was compared to the following engines: Geant4.9.5, PENH, MCsquare, as well as the MC and PBA algorithms of RayStation (RaySearch Laboratories AB). RESULTS: For a γ-index criteria of 2%/2mm, the passing rates are 93.8% for Geant4.9.5, 97.4% for PENH, 93.4% for MCsquare, 95.9% for RayStation MC, and 44.7% for PBA. The differences in γ-index passing rates between MC and RayStation PBA calculations can exceed 50%. CONCLUSION: The performance of dose calculation algorithms in highly inhomogeneous media was evaluated in a dedicated experiment. MC dose engines performed overall satisfactorily while large deviations were observed with PBA as expected.

The effect of electron collimator leaf shape on the build-up dose in narrow electron MLC fields.

Previously, we have found that the build-up dose from abutting narrow electron beams formed with unfocussed electron multi-leaf collimator (eMLC) steal leaves was higher than with the respective open field. To investigate more closely the effect of leaf material and shape on dose in the build-up region, straight, round (radius 1.5 cm) and leaf ends with a different front face angle of alpha (leaf front face pointing towards the beam axis at an angle of 90 - alpha) made of steel, brass and tungsten were modelled using the BEAMnrc code. Based on a treatment head simulation of a Varian 2100 C/D linac, depth-dose curves and profiles in water were calculated for narrow 6, 12 and 20 MeV eMLC beams (width 1.0 cm, length 10 cm) at source-to-surface distances (SSD) of 102 and 105 cm. The effects of leaf material and front face angle were evaluated based on electron fluence, angle and energy spectra. With a leaf front face angle of 15 degrees, the dose in the build-up region of the 6 MeV field varied between 91 and 100%, while for straight and round leaf shapes the dose varied between 89 and 100%. The variation was between 94 and 100% for 12 and 20 MeV. For abutting narrow 6 MeV fields with total field size 5 x 10 cm(2), the build-up doses at 5 mm depth for the face angle 15 degrees and straight and round leaf shapes were 96% and 86% (SSD 102 cm) and 89% and 85% (SSD 105 cm). With higher energies, the effect of eMLC leaf shape on dose at 5 mm was slight (3-4% units with 12 MeV) and marginal with 20 MeV. The fluence, energy and angle spectra for total and leaf scattered electrons were practically the same for different leaf materials with 6 MeV. With high energies, the spectra for tungsten were more peaked due to lower leaf transmission. Compared with straight leaf ends, the face angle of 15 degrees and round leaf ends led to a 1 mm (for 6 MeV) and between 1 and 5 mm (12 and 20 MeV at a SSD of 105 cm) decrease of therapeutic range and increase of the field size, respectively. However, profile flatness was better for abutting 6 MeV beams with round (2.5%) and face angle 15 degrees leaves (3.0%) compared to straight leaf shape (5.2%). The eMLC leaves with a face angle of 15 degrees resulted in a marked increase in the build-up dose for the single narrow eMLC beam and thus in the dose in the build-up region from matched abutting fields.

Enhancement of electron-beam surface dose with an electron multi-leaf collimator (eMLC): a feasibility study.

Use of a water-equivalent bolus in electron-beam radiotherapy is sometimes impractical and non-hygienic. Therefore, the feasibility of applying adjacent narrow beams for producing high surface dose electron beams without a bolus was investigated. Depth dose curves and profiles in water were calculated and measured for 6 and 9 MeV electron-beam segments (width 0.3-1.5 cm, length 10 cm) for source-to-surface distances (SSD) 102 and 105 cm. Segment shaping was performed with an add-on electron multi-leaf collimator prototype attached to the Varian 2100 C/D linac. Dose calculations were performed with the Voxel Monte Carlo++ algorithm. Resulting dose distributions in typical clinical cases were compared with the bolus technique. With a composite segmental field with 1.0 cm wide segments the surface dose was over 90% of the depth dose maximum for both energies. The build-up area practically disappeared with a 0.5 cm wide single beam. This led to decrease in the therapeutic range for composite fields with segment widths smaller than 1.0 cm. The new technique yielded similar surface doses as the bolus technique. The photon contamination was 4% with a 9 x 10 cm(2) field (1.0 cm wide segments) compared to 1% for the respective open field with 9 MeV with a bolus. The calculated dose agreed within 2 mm and 3% of the measured dose in 93.7% and 85.2% of the voxels. Adjacent narrow eMLC beams with a 1.0 cm width are suitable to produce electron fields with high surface dose. Despite a slight nonuniformity in the surface profiles in the lateral part of the field at SSD 102 cm, surface dose and target coverage are comparable with the bolus technique.

Dosimetric verification of a Monte Carlo electron beam model for an add-on eMLC.

Dosimetric verification of a new Monte Carlo beam model for multi-leaf collimated electrons was performed using experimental data from an add-on electron multi-leaf collimator (eMLC) prototype. The measurements were compared against calculations using an electron phase space sampled from a parameterized electron beam model and the voxel Monte Carlo++ (VMC++) code for in-phantom energy deposition. Verification of the calculations was performed in a water phantom with the developed eMLC attached to a Varian 2100 C/D radiotherapy accelerator with nominal energies 6 MeV, 9 MeV, 12 MeV, 16 MeV and 20 MeV. The eMLC prototype consisting of 2 cm thick and 5 mm wide steel leaves is fixed under the 20 x 20 cm(2) electron applicator with a source-to-leaf distance 97.2 cm. The eMLC prototype has non-motorized leaves with straight leaf edges and a maximum field size of 20 x 20 cm(2) at SSD 100 cm. The beam model is a coupled multi-source model with parameters derived from detailed beam characterization measurements and a kernel model for the indirect leaf-scattered electrons. Typical calculation times with a 2% mean statistical uncertainty was under 5 min. In extensive set of in-water measurements 88% of the voxels were within 2% /2 mm acceptance criterion. Although at SSD 100 cm the dose near the phantom surface is slightly pronounced due to the short collimator-to-surface distance, the new beam model was suitable for dose calculation of the add-on type eMLC.

Observation of the alpha particle "Knock-On" neutron emission from magnetically confined DT fusion plasmas

Suprathermal fuel ions from alpha-particle knock-on collisions in fusion DT plasmas are predicted to cause a weak feature in the neutron spectrum of d+t-->alpha+n. The knock-on feature has been searched for in the neutron emission of high ( >1 MW) fusion-power plasmas produced at JET and was found using a magnetic proton recoil type neutron spectrometer of high performance. Measurement and predictions agree both in absolute amplitude and in plasma-parameter dependence, supporting the interpretation and model. Moreover, the results provide input to projecting alpha-particle diagnostics for future self-heated fusion plasmas.