Uncertainties in ocular proton planning and their impact on required margins.

PURPOSE: To review required margins in ocular proton therapy (OPT) based on an uncertainty estimation and to compare them with widely used values. Further, uncertainties when using registered funduscopy images in the 3D model is investigated. METHODS: An uncertainty budget in planning and delivery was defined to determine required aperture and range margins. Setup uncertainties were considered for a cohort of treated patients and tested in a worst-case estimation. Other uncertainties were based on a best-guess and knowledge of institutional specifics, e.g. range reproducibility. Margins for funduscopy registration were defined resulting from scaling, rotation and translation of the image. Image formation for a wide-field fundus camera was reviewed and compared to the projection employed in treatment planning systems. RESULTS: Values for aperture and range with margins of 2.5 mm as reported in literature could be determined. Aperture margins appear appropriate for setup uncertainties below 0.5 mm, but depend on lateral penumbra. Range margins depend on depth and associated density uncertainty in tissue. Registration of funduscopy images may require margins of >2 mm, increasing towards the equator. Difference in the projection may lead to discrepancies of several mm. CONCLUSIONS: The commonly used 2.5 mm aperture margin was validated as an appropriate choice, while range margins could be reduced for lower ranges. Margins may however not include uncertainties in contouring and possible microscopic spread. If a target base is contoured on registered funduscopy images care must be taken as they are subject to larger uncertainties. Multimodal imaging approach in OPT remains advisable.

Consistency of Faraday cup and ionization chamber dosimetry of proton fields and the role of nuclear interactions.

BACKGROUND: A Faraday cup (FC) facilitates a quite clean measurement of the proton fluence emerging from clinical spot-scanning nozzles with narrow pencil-beams. The utilization of FCs appears to be an attractive option for high dose rate delivery modes and the source models of Monte-Carlo (MC) dose engines. However, previous studies revealed discrepancies of 3%-6% between reference dosimetry with ionization chambers (ICs) and FC-based dosimetry. This has prevented the widespread use of FCs for dosimetry in proton therapy. PURPOSE: The current study aims at bridging the gap between FC dosimetry and IC dosimetry of proton fields delivered with spot-scanning treatment heads. Particularly, a novel method to evaluate FC measurements is introduced. METHODS: A consistency check is formulated, which makes use of the energy balance and the reciprocity theorem. The measurement data comprise central-axis depth distributions of the absorbed dose of quasi-monochromatic fields with a width of about 28.5 cm and FC measurements of the reciprocal fields with a single spot. These data are complemented by a look-up of energy-range tables, the average Q-value of transmutations, and the escape energy carried away by neutrons and photons. The latter data are computed by MC simulations, which in turn are validated with measurements of the distal dose tail and neutron out-of-field doses. For comparison, the conventional approach of FC evaluation is performed, which computes absorbed dose from the product of fluence and stopping power. The results from the FC measurements are compared with the standard dosimetry protocols and improved reference dosimetry methods. RESULTS: The deviation between the conventional FC-based dosimetry and the IC-based one according to standard dosimetry protocols was -4.7 ( ± $\pm$ 3.3)% for a 100 MeV field and -3.6 ( ± $\pm$ 3.5)% for 200 MeV, thereby agreeing within the reported uncertainties. The deviations could be reduced to -4.0 ( ± $\pm$ 2.9)% and -3.0 ( ± $\pm$ 3.1)% by adopting state-of-the-art reference dosimetry methods. The alternative approach using the energy balance gave deviations of only -1.9% (100 MeV) and -2.6% (200 MeV) using state-of-the-art dosimetry. The standard uncertainty of this novel approach was estimated to be about 2%. CONCLUSIONS: An alternative concept has been established to determine the absorbed dose of monoenergetic proton fields with an FC. It eliminates the strong dependence of the conventional FC-based approach on the MC simulation of the stopping-power and of the secondary ions, which according to the study at hand is the major contributor to the underestimation of the absorbed dose. Some contributions to the uncertainty of the novel approach could potentially be reduced in future studies. This would allow for accurate consistency tests of conventional dosimetry procedures.

The dirty and clean dose concept: Towards creating proton therapy treatment plans with a photon-like dose response.

BACKGROUND: Applying tolerance doses for organs at risk (OAR) from photon therapy introduces uncertainties in proton therapy when assuming a constant relative biological effectiveness (RBE) of 1.1. PURPOSE: This work introduces the novel dirty and clean dose concept, which allows for creating treatment plans with a more photon-like dose response for OAR and, thus, less uncertainties when applying photon-based tolerance doses. METHODS: The concept divides the 1.1-weighted dose distribution into two parts: the clean and the dirty dose. The clean and dirty dose are deposited by protons with a linear energy transfer (LET) below and above a set LET threshold, respectively. For the former, a photon-like dose response is assumed, while for the latter, the RBE might exceed 1.1. To reduce the dirty dose in OAR, a MaxDirtyDose objective was added in treatment plan optimization. It requires setting two parameters: LET threshold and max dirty dose level. A simple geometry consisting of one target volume and one OAR in water was used to study the reduction in dirty dose in the OAR depending on the choice of the two MaxDirtyDose objective parameters during plan optimization. The best performing parameter combinations were used to create multiple dirty dose optimized (DDopt) treatment plans for two cranial patient cases. For each DDopt plan, 1.1-weighted dose, variable RBE-weighted dose using the Wedenberg RBE model and dose-average LETd distributions as well as resulting normal tissue complication probability (NTCP) values were calculated and compared to the reference plan (RefPlan) without MaxDirtyDose objectives. RESULTS: In the water phantom studies, LET thresholds between 1.5 and 2.5 keV/µm yielded the best plans and were subsequently used. For the patient cases, nearly all DDopt plans led to a reduced Wedenberg dose in critical OAR. This reduction resulted from an LET reduction and translated into an NTCP reduction of up to 19 percentage points compared to the RefPlan. The 1.1-weighted dose in the OARs was slightly increased (patient 1: 0.45 Gy(RBE), patient 2: 0.08 Gy(RBE)), but never exceeded clinical tolerance doses. Additionally, slightly increased 1.1-weighted dose in healthy brain tissue was observed (patient 1: 0.81 Gy(RBE), patient 2: 0.53 Gy(RBE)). The variation of NTCP values due to variation of α/ß from 2 to 3 Gy was much smaller for DDopt (2 percentage points (pp)) than for RefPlans (5 pp). CONCLUSIONS: The novel dirty and clean dose concept allows for creating biologically more robust proton treatment plans with a more photon-like dose response. The reduced uncertainties in RBE can, therefore, mitigate uncertainties introduced by using photon-based tolerance doses for OAR.

Complete patient exposure during paediatric brain cancer treatment for photon and proton therapy techniques including imaging procedures.

Background: In radiotherapy, especially when treating children, minimising exposure of healthy tissue can prevent the development of adverse outcomes, including second cancers. In this study we propose a validated Monte Carlo framework to evaluate the complete patient exposure during paediatric brain cancer treatment. Materials and methods: Organ doses were calculated for treatment of a diffuse midline glioma (50.4 Gy with 1.8 Gy per fraction) on a 5-year-old anthropomorphic phantom with 3D-conformal radiotherapy, intensity modulated radiotherapy (IMRT), volumetric modulated arc therapy (VMAT) and intensity modulated pencil beam scanning (PBS) proton therapy. Doses from computed tomography (CT) for planning and on-board imaging for positioning (kV-cone beam CT and X-ray imaging) accounted for the estimate of the exposure of the patient including imaging therapeutic dose. For dose calculations we used validated Monte Carlo-based tools (PRIMO, TOPAS, PENELOPE), while lifetime attributable risk (LAR) was estimated from dose-response relationships for cancer induction, proposed by Schneider et al. Results: Out-of-field organ dose equivalent data of proton therapy are lower, with doses between 0.6 mSv (testes) and 120 mSv (thyroid), when compared to photon therapy revealing the highest out-of-field doses for IMRT ranging between 43 mSv (testes) and 575 mSv (thyroid). Dose delivered by CT ranged between 0.01 mSv (testes) and 72 mSv (scapula) while a single imaging positioning ranged between 2 µSv (testes) and 1.3 mSv (thyroid) for CBCT and 0.03 µSv (testes) and 48 µSv (scapula) for X-ray. Adding imaging dose from CT and daily CBCT to the therapeutic demonstrated an important contribution of imaging to the overall radiation burden in the course of treatment, which is subsequently used to predict the LAR, for selected organs. Conclusion: The complete patient exposure during paediatric brain cancer treatment was estimated by combining the results from different Monte Carlo-based dosimetry tools, showing that proton therapy allows significant reduction of the out-of-field doses and secondary cancer risk in selected organs.

Monte Carlo calculated ionization chamber correction factors in clinical proton beams - deriving uncertainties from published data.

For the update of the IAEA TRS-398 Code of Practice (CoP), global ionization chamber factors (fQ) and beam quality correction factors (kQ) for air-filled ionization chambers in clinical proton beams have been calculated with different Monte Carlo codes. In this study, average Monte Carlo calculated fQ and kQ factors are provided and the uncertainty of these factors is estimated. Average fQ factors in monoenergetic proton beams with energies between 60 MeV and 250 MeV were derived from Monte Carlo calculated fQ factors published in the literature. Altogether, 195 fQ factors for six plane-parallel and three cylindrical ionization chambers calculated with penh, fluka and geant4 were incorporated. Additionally, a weighted standard deviation of fQ factors was calculated, where the same weight was assigned to each Monte Carlo code. From average fQ factors, kQ factors were derived and compared to the values from the IAEA TRS-398 CoP published in 2000 as well as to the values of the upcoming version. Average Monte Carlo calculated fQ factors are constant within 0.6% over the energy range investigated. In general, the different Monte Carlo codes agree within 1% for low energies and show larger differences up to 2% for high energies. As a result, the standard deviation of fQ factors increases with energy and is ∼0.3% for low energies and ∼0.8% for high energies. kQ factors derived from average Monte Carlo calculated fQ factors differ from the values presented in the IAEA TRS-398 CoP by up to 2.4%. The overall estimated uncertainty of Monte Carlo calculated kQ factors is ∼0.5%-1% smaller than the uncertainties estimated in IAEA TRS-398 CoP since the individual ionization chamber characteristics (e.g. fluence perturbations) are considered in detail in Monte Carlo calculations. The agreement between Monte Carlo calculated kQ factors and the values of the upcoming version of IAEA TRS-398 CoP is better with deviations smaller than 1%.

Comprehensive investigation of lateral dose profile and output factor measurements in small proton fields from different delivery techniques.

BACKGROUND AND PURPOSE: As a part of the commissioning and quality assurance in proton beam therapy, lateral dose profiles and output factors have to be acquired. Such measurements can be performed with point detectors and are especially challenging in small fields or steep lateral penumbra regions as the detector's volume effect may lead to perturbations. To address this issue, this work aims to quantify and correct for such perturbations of six point detectors in small proton fields created via three different delivery techniques. METHODS: Lateral dose profile and output measurements of three proton beam delivery techniques (pencil beam scanning, pencil beam scanning combined with collimators, passive scattering with collimators) were performed using high-resolution EBT3 films, a PinPoint 3D 31022 ionization chamber, a microSilicon diode 60023 and a microDiamond detector 60019 (all PTW Freiburg, Germany). Detector specific lateral dose response functions K(x,y) acting as the convolution kernel transforming the undisturbed dose distribution D(x,y) into the measured signal profiles M(x,y) were applied to quantify perturbations of the six investigated detectors in the proton fields and correct the measurements. A signal theoretical analysis in Fourier space of the dose distributions and detector's K(x,y) was performed to aid the understanding of the measurement process with regard to the combination of detector choice and delivery technique. RESULTS: Quantification of the lateral penumbra broadening and signal reduction at the fields center revealed that measurements in the pencil beam scanning fields are only compromised slightly even by large volume ionization chambers with maximum differences in the lateral penumbra of 0.25 mm and 4% signal reduction at the field center. In contrast, radiation techniques with collimation are not accurately represented by the investigated detectors as indicated by a penumbra broadening up to 1.6 mm for passive scattering with collimators and 2.2 mm for pencil beam scanning with collimators. For a 3 mm diameter collimator field, a signal reduction at field center between 7.6% and 60.7% was asserted. Lateral dose profile measurements have been corrected via deconvolution with the corresponding K(x,y) to obtain the undisturbed D(x,y). Corrected output ratios of the passively scattered collimated fields obtained for the microDiamond, microSilicon and PinPoint 3D show agreement better than 0.9% (one standard deviation) for the smallest field size of 3 mm. CONCLUSION: Point detector perturbations in small proton fields created with three delivery techniques were quantified and found to be especially pronounced for collimated small proton fields with steep dose gradients. Among all investigated detectors, the microSilicon diode showed the smallest perturbations. The correction strategies based on detector's K(x,y) were found suitable for obtaining unperturbed lateral dose profiles and output factors. Approximation of K(x,y) by considering only the geometrical averaging effect has been shown to provide reasonable prediction of the detector's volume effect. The findings of this work may be used to guide the choice of point detectors in various proton fields and to contribute toward the development of a code of practice for small field proton dosimetry.

Is an Endorectal Balloon Beneficial for Rectal Sparing after Spacer Implantation in Prostate Cancer Patients Treated with Hypofractionated Intensity-Modulated Proton Beam Therapy? A Dosimetric and Radiobiological Comparison Study.

BACKGROUND: The aim of this study is to examine the dosimetric influence of endorectal balloons (ERB) on rectal sparing in prostate cancer patients with implanted hydrogel rectum spacers treated with dose-escalated or hypofractionated intensity-modulated proton beam therapy (IMPT). METHODS: Ten patients with localized prostate cancer included in the ProRegPros study and treated at our center were investigated. All patients underwent placement of hydrogel rectum spacers before planning. Two planning CTs (with and without 120 cm3 fluid-filled ERB) were applied for each patient. Dose prescription was set according to the h strategy, with 72 Gray (Gy)/2.4 Gy/5× weekly to prostate + 1 cm of the seminal vesicle, and 60 Gy/2 Gy/5× weekly to prostate + 2 cm of the seminal vesicle. Planning with two laterally opposed IMPT beams was performed in both CTs. Rectal dosimetry values including dose-volume statistics and normal tissue complication probability (NTCP) were compared for both plans (non-ERB plans vs. ERB plans). RESULTS: For ERB plans compared with non-ERB, the reductions were 8.51 ± 5.25 Gy (RBE) (p = 0.000) and 15.76 ± 11.11 Gy (p = 0.001) for the mean and the median rectal doses, respectively. No significant reductions in rectal volumes were found after high dose levels. The use of ERB resulted in significant reduction in rectal volume after receiving 50 Gy (RBE), 40 Gy (RBE), 30 Gy (RBE), 20 Gy (RBE), and 10 Gy (RBE) with p values of 0.034, 0.008, 0.003, 0.001, and 0.001, respectively. No differences between ERB and non-ERB plans for the anterior rectum were observed. ERB reduced posterior rectal volumes in patients who received 30 Gy (RBE), 20 Gy (RBE), or 10 Gy (RBE), with p values of 0.019, 0.003, and 0.001, respectively. According to the NTCP models, no significant reductions were observed in mean or median rectal toxicity (late rectal bleeding ≥ 2, necrosis or stenosis, and late rectal toxicity ≥ 3) when using the ERB. CONCLUSION: ERB reduced rectal volumes exposed to intermediate or low dose levels. However, no significant reduction in rectal volume was observed in patients receiving high or intermediate doses. There was no benefit and also no disadvantage associated with the use of ERB for late rectal toxicity, according to available NTCP models.

Optimization of proton pencil beam positioning in collimated fields.

BACKGROUND: The addition of static or dynamic collimator systems to the pencil beam scanning delivery technique increases the number of options for lateral field shaping. The collimator shape needs to be optimized together with the intensity modulation of spots. PURPOSE: To minimize the proton field's lateral penumbra by investigating the fundamental relations between spot and collimating aperture edge position. METHODS: Analytical approaches describing the effect of spot position on the resulting spot profile are presented. The theoretical description is then compared with Monte Carlo simulations in TOPAS and in the RayStation treatment planning system, as well as with radiochromic film measurements at a clinical proton therapy facility. In the model, one single spot profile is analyzed for various spot positions in air. Further, irradiation setups in water with different energies, the combination with a range shifter, and two-dimensional proton fields were investigated in silico. RESULTS: The further the single spot is placed beyond the collimating aperture edge ('overscanning'), the sharper the relative lateral dose fall-off and thus the lateral penumbra. Overscanning up to 5 mm $5\,\text{mm}$ reduced the lateral penumbra by about 20% on average after a propagation of 13 cm $13\,\text{cm}$ in air. This benefit from overscanning is first predicted by the analytical proofs and later verified by simulations and measurements. Corresponding analyses in water confirm the benefit in lateral penumbra with spot position optimization as observed theoretically and in air. The combination of spot overscanning with fluence modulation facilitated an additional improvement. CONCLUSIONS: The lateral penumbra of single spots in collimated scanned proton fields can be improved by the method of spot overscanning. This suggests a better sparing of proximal organs at risk in smaller water depths at higher energies, especially in the plateau of the depth dose distribution. All in all, spot overscanning in collimated scanned proton fields offers particular potential in combination with techniques such as fluence modulation or dynamic collimation for optimizing the lateral penumbra to spare normal tissue.

Overcoming inter-observer planning variability in target volume contouring and dose planning for high-risk neuroblastoma - a European multicenter effort of the SIOPEN radiotherapy committee.

BACKGROUND AND PURPOSE: To establish an international quality standard for contouring and planning for high-risk neuroblastoma within the prospective High-Risk Neuroblastoma Study 2 of SIOP-Europe-Neuroblastoma (SIOPEN HR-NBL2), which includes a randomized question on dose escalation for residual disease. MATERIALS AND METHODS: Data on four patients with high-risk neuroblastoma were selected and distributed to the radiotherapy committee of the HR-NBL2 study for independent contouring and planning. Differences in contouring were analyzed using apparent and kappa-corrected agreement. Plans were analyzed regarding the dose-volume histogram metrics. Results were discussed among experts and agreement was obtained. RESULTS: Substantial agreement was found for contouring of the heart (0.64), liver (0.70), left lung (0.74), and right lung (0.74). For contouring of the gastrointestinal tract (0.54), left kidney (0.60), and right kidney (0.59) moderate agreement was obtained. For target volume delineation, agreement for preoperative tumour extent was moderate (0.42), for CTV fair (0.35) and only low (0.06) for residual tumour, respectively. The dose planning strategies appeared to be relatively homogeneous among all experts. CONCLUSION: Considerable variability was found for the delineation of target volumes, particularly the boost volume, whereas the contouring of the organs at risk and the planning strategy were reasonably consistent. In order to obtain reliable results from the randomized HR-NBL2 trial, standardization of target volume delineation based on adequate imaging is crucial.

Validating a double Gaussian source model for small proton fields in a commercial Monte-Carlo dose calculation engine.

PURPOSE: The primary fluence of a proton pencil beam exiting the accelerator is enveloped by a region of secondaries, commonly called "spray". Although small in magnitude, this spray may affect dose distributions in pencil beam scanning mode e.g., in the calculation of the small field output, if not modelled properly in a treatment planning system (TPS). The purpose of this study was to dosimetrically benchmark the Monte Carlo (MC) dose engine of the RayStation TPS (v.10A) in small proton fields and systematically compare single Gaussian (SG) and double Gaussian (DG) modeling of initial proton fluence providing a more accurate representation of the nozzle spray. METHODS: The initial proton fluence distribution for SG/DG beam modeling was deduced from two-dimensional measurements in air with a scintillation screen with electronic readout. The DG model was either based on direct fits of the two Gaussians to the measured profiles, or by an iterative optimization procedure, which uses the measured profiles to mimic in-air scan-field factor (SF) measurements. To validate the DG beam models SFs, i.e. relative doses to a 10 × 10 cm2 field, were measured in water for three different initial proton energies (100MeV, 160MeV, 226.7MeV) and square field sizes from 1×1cm2 to 10×10cm2 using a small field ionization chamber (IBA CC01) and an IBA ProteusPlus system (universal nozzle). Furthermore, the dose to the center of spherical target volumes (diameters: 1cm to 10cm) was determined using the same small volume ionization chamber (IC). A comprehensive uncertainty analysis was performed, including estimates of influence factors typical for small field dosimetry deduced from a simple two-dimensional analytical model of the relative fluence distribution. Measurements were compared to the predictions of the RayStation TPS. RESULTS: SFs deviated by more than 2% from TPS predictions in all fields <4×4cm2 with a maximum deviation of 5.8% for SG modeling. In contrast, deviations were smaller than 2% for all field-sizes and proton energies when using the directly fitted DG model. The optimized DG model performed similarly except for slightly larger deviations in the 1×1cm2 scan-fields. The uncertainty estimates showed a significant impact of pencil beam size variations (±5%) resulting in up to 5.0% standard uncertainty. The point doses within spherical irradiation volumes deviated from calculations by up to 3.3% for the SG model and 2.0% for the DG model. CONCLUSION: Properly representing nozzle spray in RayStation's MC-based dose engine using a DG beam model was found to reduce the deviation to measurements in small spherical targets to below 2%. A thorough uncertainty analysis shows a similar magnitude for the combined standard uncertainty of such measurements.

Commissioning and validation of a novel commercial TPS for ocular proton therapy.

BACKGROUND: Until today, the majority of ocular proton treatments worldwide were planned with the EYEPLAN treatment planning system (TPS). Recently, the commercial, computed tomography (CT)-based TPS for ocular proton therapy RayOcular was released, which follows the general concepts of model-based treatment planning approach in conjunction with a pencil-beam-type dose algorithm (PBA). PURPOSE: To validate RayOcular with respect to two main features: accurate geometrical representation of the eye model and accuracy of its dose calculation algorithm in combination with an Ion Beam Applications (IBA) eye treatment delivery system. METHODS: Different 3D-printed eye-ball-phantoms were fabricated to test the geometrical representation of the corresponding CT-based model, both in orthogonal 2D images for X-ray image overlay and in fundus view overlaid with a funduscopy. For the latter, the phantom was equipped with a lens matching refraction of the human eye. Funduscopy was acquired in a Zeiss Claus 500 camera. Tantalum clips and fiducials attached to the phantoms were localized in the TPS model, and residual deviations to the actual position in X-ray images for various orientations of the phantom were determined, after the nominal eye orientation was corrected in RayOcular to obtain a best overall fit. In the fundus view, deviations between known and displayed distances were measured. Dose calculation accuracy of the PBA on a 0.2 mm grid was investigated by comparing between measured lateral and depth-dose profiles in water for various combinations of range, modulation, and field-size. Ultimately, the modeling of dose distributions behind wedges was tested. A 1D gamma-test was applied, and the lateral and distal penumbra were further compared. RESULTS: Average residuals between model clips and visible clips/fiducials in orthogonal X-ray images were within 0.3 mm, including different orientations of the phantom. The differences between measured distances on the registered funduscopy image in the RayOcular fundus view and the known ground-truth were within 1 mm up to 10.5 mm distance from the posterior pole. No clear benefit projection of either polar mode or camera mode could be identified, the latter mimicking camera properties. Measured dose distributions were reproduced with gamma-test pass-rates of >95% with 2%/0.3 mm for depth and lateral profiles in the middle of spread-out Bragg-peaks. Distal falloff and lateral penumbra were within 0.2 mm for fields without a wedge. For shallow depths, the agreement was worse, reaching pass-rates down to 80% with 5%/0.3 mm when comparing lateral profiles in air. This is caused by low-energy protons from a scatter source in the IBA system not modeled by RayOcular. Dose distributions modified by wedges were reproduced, matching the wedge-induced broadening of the lateral penumbra to within 0.4 mm for the investigated cases and showing the excess dose within the field due to wedge scatter. CONCLUSION: RayOcular was validated for its use with an IBA single scattering delivery nozzle. Geometric modeling of the eye and representation of 2D projections fulfill clinical requirements. The PBA dose calculation reproduces measured distributions and allows explicit handling of wedges, overcoming approximations of simpler dose calculation algorithms used in other systems.

Comparing biological effectiveness guided plan optimization strategies for cranial proton therapy: potential and challenges.

BACKGROUND: To introduce and compare multiple biological effectiveness guided (BG) proton plan optimization strategies minimizing variable relative biological effectiveness (RBE) induced dose burden in organs at risk (OAR) while maintaining plan quality with a constant RBE. METHODS: Dose-optimized (DOSEopt) proton pencil beam scanning reference treatment plans were generated for ten cranial patients with prescription doses ≥ 54 Gy(RBE) and ≥ 1 OAR close to the clinical target volume (CTV). For each patient, four additional BG plans were created. BG objectives minimized either proton track-ends, dose-averaged linear energy transfer (LETd), energy depositions from high-LET protons or variable RBE-weighted dose (DRBE) in adjacent serially structured OARs. Plan quality (RBE = 1.1) was assessed by CTV dose coverage and robustness (2 mm setup, 3.5% density), dose homogeneity and conformity in the planning target volumes and adherence to OAR tolerance doses. LETd, DRBE (Wedenberg model, α/ßCTV = 10 Gy, α/ßOAR = 2 Gy) and resulting normal tissue complication probabilities (NTCPs) for blindness and brainstem necrosis were derived. Differences between DOSEopt and BG optimized plans were assessed and statistically tested (Wilcoxon signed rank, α = 0.05). RESULTS: All plans were clinically acceptable. DOSEopt and BG optimized plans were comparable in target volume coverage, homogeneity and conformity. For recalculated DRBE in all patients, all BG plans significantly reduced near-maximum DRBE to critical OARs with differences up to 8.2 Gy(RBE) (p < 0.05). Direct DRBE optimization primarily reduced absorbed dose in OARs (average ΔDmean = 2.0 Gy; average ΔLETd,mean = 0.1 keV/µm), while the other strategies reduced LETd (average ΔDmean < 0.3 Gy; average ΔLETd,mean = 0.5 keV/µm). LET-optimizing strategies were more robust against range and setup uncertaintes for high-dose CTVs than DRBE optimization. All BG strategies reduced NTCP for brainstem necrosis and blindness on average by 47% with average and maximum reductions of 5.4 and 18.4 percentage points, respectively. CONCLUSIONS: All BG strategies reduced variable RBE-induced NTCPs to OARs. Reducing LETd in high-dose voxels may be favourable due to its adherence to current dose reporting and maintenance of clinical plan quality and the availability of reported LETd and dose levels from clinical toxicity reports after cranial proton therapy. These optimization strategies beyond dose may be a first step towards safely translating variable RBE optimization in the clinics.

Validation of a Monte Carlo Framework for Out-of-Field Dose Calculations in Proton Therapy.

Proton therapy enables to deliver highly conformed dose distributions owing to the characteristic Bragg peak and the finite range of protons. However, during proton therapy, secondary neutrons are created, which can travel long distances and deposit dose in out-of-field volumes. This out-of-field absorbed dose needs to be considered for radiation-induced secondary cancers, which are particularly relevant in the case of pediatric treatments. Unfortunately, no method exists in clinics for the computation of the out-of-field dose distributions in proton therapy. To help overcome this limitation, a computational tool has been developed based on the Monte Carlo code TOPAS. The purpose of this work is to evaluate the accuracy of this tool in comparison to experimental data obtained from an anthropomorphic phantom irradiation. An anthropomorphic phantom of a 5-year-old child (ATOM, CIRS) was irradiated for a brain tumor treatment in an IBA Proteus Plus facility using a pencil beam dedicated nozzle. The treatment consisted of three pencil beam scanning fields employing a lucite range shifter. Proton energies ranged from 100 to 165 MeV. A median dose of 50.4 Gy(RBE) with 1.8 Gy(RBE) per fraction was prescribed to the initial planning target volume (PTV), which was located in the cerebellum. Thermoluminescent detectors (TLDs), namely, Li-7-enriched LiF : Mg, Ti (MTS-7) type, were used to detect gamma radiation, which is produced by nuclear reactions, and secondary as well as recoil protons created out-of-field by secondary neutrons. Li-6-enriched LiF : Mg,Cu,P (MCP-6) was combined with Li-7-enriched MCP-7 to measure thermal neutrons. TLDs were calibrated in Co-60 and reported on absorbed dose in water per target dose (µGy/Gy) as well as thermal neutron dose equivalent per target dose (µSv/Gy). Additionally, bubble detectors for personal neutron dosimetry (BD-PND) were used for measuring neutrons (>50 keV), which were calibrated in a Cf-252 neutron beam to report on neutron dose equivalent dose data. The Monte Carlo code TOPAS (version 3.6) was run using a phase-space file containing 1010 histories reaching an average standard statistical uncertainty of less than 0.2% (coverage factor k = 1) on all voxels scoring more than 50% of the maximum dose. The primary beam was modeled following a Fermi-Eyges description of the spot envelope fitted to measurements. For the Monte Carlo simulation, the chemical composition of the tissues represented in ATOM was employed. The dose was tallied as dose-to-water, and data were normalized to the target dose (physical dose) to report on absorbed doses per target dose (mSv/Gy) or neutron dose equivalent per target dose (µSv/Gy), while also an estimate of the total organ dose was provided for a target dose of 50.4 Gy(RBE). Out-of-field doses showed absorbed doses that were 5 to 6 orders of magnitude lower than the target dose. The discrepancy between TLD data and the corresponding scored values in the Monte Carlo calculations involving proton and gamma contributions was on average 18%. The comparison between the neutron equivalent doses between the Monte Carlo simulation and the measured neutron doses was on average 8%. Organ dose calculations revealed the highest dose for the thyroid, which was 120 mSv, while other organ doses ranged from 18 mSv in the lungs to 0.6 mSv in the testes. The proposed computational method for routine calculation of the out-of-the-field dose in proton therapy produces results that are compatible with the experimental data and allow to calculate out-of-field organ doses during proton therapy.

Dosimetric feasibility of moderately hypofractionated/dose escalated radiation therapy for localised prostate cancer with intensity-modulated proton beam therapy using simultaneous integrated boost (SIB-IMPT) and impact of hydrogel prostate-rectum spacer.

PURPOSE: To examine the dosimetric feasibility of hypofractionated/dose escalated radiation therapy in patients with localized prostate carcinoma using simultaneous integrated boost intensity-modulated proton beam therapy (SIB-IMPT) in absence or presence of prostate-rectum spacer. METHODS: IMPT technique was implemented in 23 patients with intermediate- and high-risk prostate cancer treated at West German Proton Therapy Centre from March 2016 till June 2018, using SIB technique prescribing 60 GyRBE and 72 GyRBE in 30 fractions to PTV1 (prostate and seminal vesicle) and PTV2 boost (prostate and proximal seminal vesicle), respectively. In 15 patients, a transperineal injection of hydrogel was applied prior to radiotherapy to increase the distance between prostate and rectum. Planning and all treatments were performed with a 120 ml fluid-filled endorectal balloon customised daily for each patient. For each patient, 2 lateral IMPT beams were implemented taking a field-specific range uncertainty (RU) into account. Dose volume histograms (DVH) were analyzed for PTV2, PTV2 with range uncertainty margin (PTV2RU), rectum, bladder, right/left femoral heads, and penile bulb. For late rectal toxicities, the normal tissue complication probabilities (NTCP) were calculated using different biological models. A DVH- and NTCP-based dosimetric comparison was carried out between non-spacer and spacer groups. RESULTS: For the 23 patients, high-quality plans could be achieved for target volume and for other organs at risk (OARs). For PTV2, the V107% was 0% and the Dmax did not exceed 106.2% of the prescribed dose. The volume PTV2RU covered by 95% of the dose ranged from 96.16 to 99.95%. The conformality index for PTV2RU was 1.12 ± 0.057 and the homogeneity index (HI) was 1.04 ± 0.014. Rectum Dmax and rectal volume receiving 73-50 Gy could be further reduced for the spacer-group. Significant reductions in mean and median rectal NTCPs (stenosis/necrosis, late rectal bleeding ≥ 2, and late rectal toxicities ≥ 3) were predicted for the spacer group in comparison to the non-spacer group. CONCLUSION: Hypofractionated/dose escalated radiotherapy with SIB-IMPT is dosimetrically feasible. Further reduction of the rectal volumes receiving high and medium dose levels (73-50 Gy) and rectal NTCP could be achieved through injection of spacers between rectum and prostate.

Technical note: Impact of beam properties for uveal melanoma proton therapy-An in silico planning study.

PURPOSE: To evaluate the impact of beam quality in terms of distal fall-off (DFO, 90%-10%) and lateral penumbra (LP, 80%-20%) of single beam ocular proton therapy (OPT) and to derive resulting ideal requirements for future systems. METHODS: Nine different beam models with DFO varying between 1 and 4 mm and LP between 1 and 4 mm were created. Beam models were incorporated into the RayStation with RayOcular treatment planning system version 10 B (RaySearch Laboratories, Stockholm, Sweden). Each beam model was applied for eight typical clinical cases, covering different sizes and locations of uveal melanoma. Plans with and without an additional wedge were created, resulting in 117 plans with a total prescribed median dose of 60 Gy(RBE) to the clinical target volume. Treatment plans were analyzed in terms of V20-V80 penumbra volume, D1 (dose to 1% of the volume) for optic disc and macula, optic nerve V30 (volume receiving 30 Gy(RBE), i.e., 50% of prescription), as well as average dose to lens and ciliary body. An LP-dependent aperture margin was based on estimated uncertainties, ranging from 1.7 to 4.0 mm. RESULTS: V20-V80 showed a strong influence by LP, while DFO was less relevant. The optic disc D1 reached an extra dose of up to 3000 cGy(RBE), comparing the defined technical limit of DFO = LP = 1 mm with DFO = 3 mm/LP = 4 mm. The latter may result from a pencil-beam scanning (PBS) system with static apertures. Plans employing a wedge showed an improvement for organs at risk sparing. CONCLUSION: Plan quality is strongly influenced by initial beam parameters. The impact of LP is more pronounced when compared to DFO. The latter becomes important in the treatment of posterior tumors near the macula, optic disc or optic nerve. The plan quality achieved by dedicated OPT nozzles in single- or double-scattering design might not be achievable with modified PBS systems.

Mitigation of motion effects in pencil-beam scanning - Impact of repainting on 4D robustly optimized proton treatment plans for hepatocellular carcinoma.

Proton fields delivered by the active scanning technique can be interfered with the intrafractional motion. This in-silico study seeks to mitigate the dosimetric impacts of motion artifacts, especially its interplay with the time-modulated dose delivery. Here four-dimensional (4d) robust optimization and dose repainting, which is the multiple application of the same field with reduced fluence, were combined. Two types of repainting were considered: layered and volumetric repainting. The time-resolved dose calculation, which is necessary to quantify the interplay effect, was integrated into the treatment planning system and validated. Nine clinical cases of hepatocellular carcinoma (HCC) showing motion in the range of 0.4-1.5cm were studied. It was found that the repainted delivery of 4D robustly optimized plans reduced the impact of interplay effect as quantified by the homogeneity index within the clinical target volume (CTV) to a tolerable level. Similarly, the fractional over- and underdosage was reduced sufficiently for some HCC cases to achieve the purpose of motion management. This holds true for both investigated types of repainting with small dosimetric advantages of volume repainting over layered repainting. Volume repainting, however, cannot be applied clinically in proton centers with slow energy changes. Thus, it served as a reference in the in-silico evaluation. It is recommended to perform the dynamic dose calculation for individual cases to judge if robust optimization in conjunction with repainting is sufficient to keep the interplay effect within bounds.

Technical note: Providing proton fields down to the few-MeV level at clinical pencil beam scanning facilities for radiobiological experiments.

PURPOSE: The adequate performance of radiobiological experiments using clinical proton beams typically requires substantial preparations to provide the appropriate setup for specific experiments. Providing radiobiologically interesting low-energy protons is a particular challenge, due to various physical effects that become more pronounced with larger absorber thickness and smaller proton energy. This work demonstrates the generation of decelerated low-energy protons from a clinical proton beam. METHODS: Monte Carlo simulations of proton energy spectra were performed for energy absorbers with varying thicknesses to reduce the energy of the clinical proton beam down to the few-MeV level corresponding to µ m-ranges. In this way, a setup with an optimum thickness of the absorber with a maximum efficiency of the proton fluence for the provisioning of low-energy protons is supposed to be found. For the specific applications of 2.5-3.3 MeV protons and α -particle range equivalent protons, the relative depth dose was measured and simulated together with the dose-averaged linear energy transfer (LETd) distribution. RESULTS: The resulting energy spectra from Monte Carlo simulations indicate an optimal absorber thickness for providing low-energy protons with maximum efficiency of proton fluence at an user-requested energy range for experiments. For instance, providing energies lower than 5 MeV, an energy spectrum with a relative total efficiency of 38.6 % to the initial spectrum was obtained with the optimal setup. The measurements of the depth dose, compared to the Monte Carlo simulations, showed that the dosimetry of low-energy protons works and protons with high LETd down to the range of α -particles can be produced. CONCLUSIONS: This work provides a method for generating all clinically and radiobiologically relevant energies - especially down to the few-MeV level - at one clinical facility with pencil beam scanning. Thereby, it enables radiobiological experiments under environmentally uniform conditions.

Estimating the modulating effect of lung tissue in particle therapy using a clinical CT voxel histogram analysis.

To treat lung tumours with particle therapy, different additional tasks and challenges in treatment planning and application have to be addressed thoroughly. One of these tasks is the quantification and consideration of the Bragg peak (BP) degradation due to lung tissue: as lung is an heterogeneous tissue, the BP is broadened when particles traverse the microscopic alveoli. These are not fully resolved in clinical CT images and thus, the effect is not considered in the dose calculation. In this work, a correlation between the CT histograms of heterogeneous material and the impact on the BP curve is presented. Different inorganic materials were scanned with a conventional CT scanner and additionally, the BP degradation was measured in a proton beam and was then quantified. A model is proposed that allows an estimation of the modulation power by performing a histogram analysis on the CT scan. To validate the model for organic samples, a second measurement series was performed with frozen porcine lunge samples. This allows to investigate the possible limits of the proposed model in a set-up closer to clinical conditions. For lung substitutes, the agreement between model and measurement is within ±0.05 mm and for the organic lung samples, within ±0.15 mm. This work presents a novel, simple and efficient method to estimate if and how much a material or a distinct region (within the lung) is degrading the BP on the basis of a common clinical CT image. Up until now, only a direct in-beam measurement of the region or material of interest could answer this question.