Development and validation of an automatic commissioning tool for the Monte Carlo dose engine in myQA iON.

Independent dose verification with Monte Carlo (MC) simulations is an important feature of proton therapy quality assurance (QA). However, clinical integration of such tools often generates an additional and complex workload for medical physicists. The preparation of the necessary clinical inputs, such as the machine beam model, should therefore be automated. In this work, a methodology for automatic MC commissioning has been devised, validated, and developed into a MATLAB tool for the users of myQA iON, the recent QA platform of IBA Dosimetry. With this workflow, all necessary parameters can easily be tuned using dedicated optimization methods. For the geometrical beam parameters (phase space), the assumption of a single or double Gaussian is made. To model the energy spectrum, a Gaussian function is assumed and parameters are optimized using either MC simulations or a library of pre-computed Bragg peaks. For the absolute dose calibration, commissioning fields can be reproduced with the dose engine to retrieve the necessary parameters. We discuss in a first time the tool efficiency and show that one can optimize all parameters in less than 4 min per energy with excellent accuracy. We then validate a beam model obtained with the tool by simulating homogeneous spread-out Bragg peaks (SOBPs) and patient QA plans previously measured in water. An average range agreement of 0.29 ± 0.34 mm is achieved for the SOBPs while 3%/3 mm local gamma passing rates reach 99.3% on average over all 62 measured patient QA planes, which is well within clinical tolerances.

Locally tuned deformation fields combination for 2D cine-MRI-based driving of 3D motion models.

PURPOSE: To target mobile tumors in radiotherapy with the recent MR-Linac hardware solutions, research is being conducted to drive a 3D motion model with 2D cine-MRI to reproduce the breathing motion in 4D. This work presents a method to combine several deformation fields using local measures to better drive 3D motion models. METHODS: The method uses weight maps, each representing the proximity with a specific area of interest. The breathing state is evaluated on cine-MRI frames in these areas and a different deformation field is estimated for each using a 2D to 3D motion model. The different deformation fields are multiplied by their respective weight maps and combined to form the final field to apply to a reference image. A global motion model is adjusted locally on the selected areas and creates a 3DCT for each cine-MRI frame. RESULTS: The 13 patients on which it was tested showed on average an improvement of the accuracy of our model of 0.71 mm for areas selected to drive the model and 0.5 mm for other areas compared to our previous method without local adjustment. The additional computation time for each region was around 40 ms on a modern laptop. CONCLUSION: The method improves the accuracy of the2D-based driving of 3D motion models. It can be used on top of existing methods relying on deformation fields. It does add some computation time but, depending on the area to deform and the number of regions of interests, offers the potential of online use.

Development of robustness evaluation strategies for enabling statistically consistent reporting.

Robustness evaluation of proton therapy treatment plans is essential for ensuring safe treatment delivery. However, available evaluation procedures feature a limited exploration of the actual robustness of the plan and generally do not provide confidence levels. This study compared established and more sophisticated robustness evaluation procedures, with quantified confidence levels. We have evaluated several robustness evaluation methods for 5 bilateral head-and-neck patients optimized considering spot scanning delivery and with a conventional CTV-to-PTV margin of 4 mm. Method (1) good practice scenario selection (GPSS) (e.g. +/- 4 mm setup error 3% range uncertainty); (2) statistically sound scenario selection (SSSS) either only on or both on and inside isoprobability hypersurface encompassing 90% of the possible errors; (3) statistically sound dosimetric selection (SSDS). In the last method, the 90% best plans were selected according to either target coverage quantified by D 95 (SSDS_D 95) or to an approximation of the final objective function (OF) used during treatment optimization (SSDS_OF). For all methods, we have considered systematic setup and systematic range errors. A mix of systematic and random setup errors were also simulated for SSDS, but keeping the same conventional margin of 4 mm. All robustness evaluations have been performed using the fast Monte Carlo dose engine MCsquare. Both SSSS strategies yielded on average very similar results. SSSS and GPSS yield comparable values for target coverage (within 0.5 Gy). The most noticeable differences were found for the CTV between GPSS, on the one hand, and SSDS_D 95 and SSDS_OF, on the other hand (average worst-case D 98 were 2.8 and 2.0 Gy larger than for GPSS, respectively). Simulating explicitly random errors in SSDS improved almost all DVH metrics. We have observed that the width of DVH-bands and the confidence levels depend on the method chosen to sample the scenarios. Statistically sound estimation of the robustness of the plan in the dosimetric space may provide an improved insight on the actual robustness of the plan for a given confidence level.

A kernel-based algorithm for the spectral fluence of clinical proton beams to calculate dose-averaged LET and other dosimetric quantities of interest.

PURPOSE: To introduce a new analytical methodology to calculate quantities of interest in particle radiotherapy inside the treatment planning system. Models are proposed to calculate dose-averaged LET (LETd) in proton radiotherapy. MATERIAL AND METHODS: A kernel-based approach for the spectral fluence of particles is developed by means of analytical functions depending on depth and lateral position. These functions are obtained by fitting them to data calculated with Monte Carlo (MC) simulations using Geant4 in liquid water for energies from 50 to 250 MeV. Contributions of primary, secondary protons and alpha particles are modeled separately. Lateral profiles and spectra are modeled as Gaussian functions to be convolved with the fluence coming from the nozzle. LETd is obtained by integrating the stopping power curves from the PSTAR and ASTAR databases weighted by the spectrum at each position. The fast MC code MCsquare is employed to benchmark the results. RESULTS: Considering the nine energies simulated, fits for the functions modeling the fluence in-depth provide an average R 2 equal to 0.998, 0.995 and 0.986 for each one of the particles considered. Fits for the Gaussian lateral functions yield average R 2 of 0.997, 0.982 and 0.993, respectively. Similarly, the Gaussian functions fitted to the computed spectra lead to average R 2 of 0.995, 0.938 and 0.902. LETd calculation in water shows mean differences of -0.007 ± 0.008 keV/µm with respect to MCsquare if only protons are considered and 0.022 ± 0.007 keV/µm including alpha particles. In a prostate case, mean difference for all voxels with dose >5% of prescribed dose is 0.28 ± 0.23 keV/µm. CONCLUSION: This new spectral fluence-based methodology allows for simultaneous calculations of quantities of interest in proton radiotherapy such as dose, LETd or microdosimetric quantities. The method also enables the inclusion of more particles by following an analogous process.

Platform for automatic patient quality assurance via Monte Carlo simulations in proton therapy.

For radiation therapy, it is crucial to ensure that the delivered dose matches the planned dose. Errors in the dose calculations done in the treatment planning system (TPS), treatment delivery errors, other software bugs or data corruption during transfer might lead to significant differences between predicted and delivered doses. As such, patient specific quality assurance (QA) of dose distributions, through experimental validation of individual fields, is necessary. These measurement based approaches, however, are performed with 2D detectors, with limited resolution and in a water phantom. Moreover, they are work intensive and often impose a bottleneck to treatment efficiency. In this work, we investigated the potential to replace measurement-based approach with a simulation-based patient specific QA using a Monte Carlo (MC) code as independent dose calculation engine in combination with treatment log files. Our developed QA platform is composed of a web interface, servers and computation scripts, and is capable to autonomously launch simulations, identify and report dosimetric inconsistencies. To validate the beam model of independent MC engine, in-water simulations of mono-energetic layers and 30 SOBP-type dose distributions were performed. Average Gamma passing ratio 99 ± 0.5% for criteria 2%/2 mm was observed. To demonstrate feasibility of the proposed approach, 10 clinical cases such as head and neck, intracranial indications and craniospinal axis, were retrospectively evaluated via the QA platform. The results obtained via QA platform were compared to QA results obtained by measurement-based approach. This comparison demonstrated consistency between the methods, while the proposed approach significantly reduced in-room time required for QA procedures.

Impact of machine log-files uncertainties on the quality assurance of proton pencil beam scanning treatment delivery.

Irradiation log-files store useful information about the plan delivery, and together with independent Monte Carlo dose engine calculations can be used to reduce the time needed for patient-specific quality assurance (PSQA). Nonetheless, machine log-files carry an uncertainty associated to the measurement of the spot position and intensity that can influence the correct evaluation of the quality of the treatment delivery. This work addresses the problem of the inclusion of these uncertainties for the final verification of the treatment delivery. Dedicated measurements performed in an IBA Proteus Plus gantry with a pencil beam scanning (PBS) dedicated nozzle have been carried out to build a 'room-dependent' model of the spot position uncertainties. The model has been obtained through interpolation of the look-up tables describing the systematic and random uncertainties, and it has been tested for a clinical case of a brain cancer patient irradiated in a dry-run. The delivered dose has been compared with the planned dose with the inclusion of the errors obtained applying the model. Our results suggest that the accuracy of the treatment delivery is higher than the spot position uncertainties obtained from the log-file records. The comparison in terms of DVHs shows that the log-reconstructed dose is compatible with the planned dose within the 95% confidence interval obtained applying our model. The initial mean dose difference between the calculated dose to the patient based on the plan and recorded data is around 1%. The difference is essentially due to the log-file uncertainties and it can be removed with a correct treatment of these errors. In conclusion our new PSQA protocol allows for a fast verification of the dose delivered after every treatment fraction through the use of machine log-files and an independent Monte Carlo dose engine. Moreover, the inclusion of log-file uncertainties in the dose calculation allows for a correct evaluation of the quality of the treatment plan delivery.

Combined influence of CT random noise and HU-RSP calibration curve nonlinearities on proton range systematic errors.

Proton range random and systematic uncertainties are the major factors undermining the advantages of proton therapy, namely, a sharp dose falloff and a better dose conformality for lower doses in normal tissues. The influence of CT artifacts such as beam hardening or scatter can easily be understood and estimated due to their large-scale effects on the CT image, like cupping and streaks. In comparison, the effects of weakly-correlated stochastic noise are more insidious and less attention is drawn on them partly due to the common belief that they only contribute to proton range uncertainties and not to systematic errors thanks to some averaging effects. A new source of systematic errors on the range and relative stopping powers (RSP) has been highlighted and proved not to be negligible compared to the 3.5% uncertainty reference value used for safety margin design. Hence, we demonstrate that the angular points in the HU-to-RSP calibration curve are an intrinsic source of proton range systematic error for typical levels of zero-mean stochastic CT noise. Systematic errors on RSP of up to 1% have been computed for these levels. We also show that the range uncertainty does not generally vary linearly with the noise standard deviation. We define a noise-dependent effective calibration curve that better describes, for a given material, the RSP value that is actually used. The statistics of the RSP and the range continuous slowing down approximation (CSDA) have been analytically derived for the general case of a calibration curve obtained by the stoichiometric calibration procedure. These models have been validated against actual CSDA simulations for homogeneous and heterogeneous synthetical objects as well as on actual patient CTs for prostate and head-and-neck treatment planning situations.

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.

[Robust treatment planning in proton therapy]. / Planification de traitement robuste en protonthérapie.

The concentration of the dose delivered by protons at the end of their path, the Bragg peak, has the potential to improve external radiotherapy treatments. Unfortunately, the main strength of the protons, their finite range, is also their greatest weakness. Any uncertainty on the range may lead to inadequate target coverage or excessive toxicity. The uncertainties have multiple origins and include, among others, ballistic errors, morphological modifications or inaccurate estimations of the physical quantities necessary to predict the proton range. Uncertainties have been part of daily practice in conventional radiotherapy with X-rays for a long time. However, dose distributions delivered with X-rays are much less sensitive to uncertainties than the ones delivered with protons. This relative insensitivity enabled the management of uncertainties through safety margins using a simple formalism. The conditions of validity of this formalism are much more restrictive for proton therapy, leading to the need of developing new tools and adapted strategies to manage accurately these uncertainties. The objective of this paper is to present a vision for the management of uncertainties in proton therapy in the continuity of formalisms established for X-rays. The latter are first summarized before discussing the necessary developments in order to consistently apply them to protons.

LET dependence of the response of a PTW-60019 microDiamond detector in a 62MeV proton beam.

This study was initiated following conclusions from earlier experimental work, performed in a low-energy carbon ion beam, indicating a significant LET dependence of the response of a PTW-60019 microDiamond detector. The purpose of this paper is to present a comparison between the response of the same PTW-60019 microDiamond detector and an IBA Roos-type ionization chamber as a function of depth in a 62MeV proton beam. Even though proton beams are considered as low linear energy transfer (LET) beams, the LET value increases slightly in the Bragg peak region. Contrary to the observations made in the carbon ion beam, in the 62MeV proton beam good agreement is found between both detectors in both the plateau and the distal edge region. No significant LET dependent response of the PTW-60019 microDiamond detector is observed consistent with other findings for proton beams in the literature, despite this particular detector exhibiting a substantial LET dependence in a carbon ion beam.