Proton transmission imaging has been proposed and investigated as imaging modality complementary to x-ray based techniques in proton beam therapy. In particular, it addresses the issue of range ...uncertainties due to the conversion of an x-ray patient computed tomography (CT) image expressed in Hounsfield Units (HU) to relative stopping power (RSP) needed as input to the treatment planning system. One approach to exploit a single proton radiographic projection is to perform a patient-specific calibration of the CT to RSP conversion curve by optimising the match between a measured and a numerically integrated proton radiography. In this work, we develop the mathematical tools needed to perform such an optimisation in an efficient and robust way. Our main focus lies on set-ups which combine pencil beam scanning with a range telescope detector, although most of our methods can be employed in combination with other set-ups as well. Proton radiographies are simulated in Monte Carlo using an idealised detector and applying the same data processing chain used with experimental data. This approach allows us to have a ground truth CT-RSP curve to compare the optimisation results with. Our results show that the parameters of the CT-RSP curve are strongly correlated when using a pencil beam based set-up, which leads to unrealistic variation in the optimised CT-RSP curves. To address this issue, we introduce a regularisation procedure which guarantees a plausible degree of smoothness in the optimised CT-RSP curves. We investigate three different methods to perform the numerical projection operation needed to generate a proton digitally reconstructed radiography. We find that the approximate and computationally faster method performs as well as the more accurate but more demanding method. We perform a Monte Carlo experiment based on a head and neck patient to evaluate the range accuracy achievable with the optimised CT-RSP curves and find an agreement with the ground truth expectation of better than . Our results further indicate that the region in the patient in which the proton radiography is acquired does not necessarily have to correspond to the treatment volume to achieve this accuracy. This is important as the imaged region could be freely chosen, e.g. in order to spare organs at risk.
In the field of radiotherapy, Monte Carlo (MC) particle transport calculations are recognized for their superior accuracy in predicting dose and fluence distributions in patient geometries compared ...to analytical algorithms which are generally used for treatment planning due to their shorter execution times. In this work, a newly developed MC-based treatment planning (MCTP) tool for proton therapy is proposed to support treatment planning studies and research applications. It allows for single-field and simultaneous multiple-field optimization in realistic treatment scenarios and is based on the MC code FLUKA. Relative biological effectiveness (RBE)-weighted dose is optimized either with the common approach using a constant RBE of 1.1 or using a variable RBE according to radiobiological input tables. A validated reimplementation of the local effect model was used in this work to generate radiobiological input tables. Examples of treatment plans in water phantoms and in patient-CT geometries together with an experimental dosimetric validation of the plans are presented for clinical treatment parameters as used at the Italian National Center for Oncological Hadron Therapy. To conclude, a versatile MCTP tool for proton therapy was developed and validated for realistic patient treatment scenarios against dosimetric measurements and commercial analytical TP calculations. It is aimed to be used in future for research and to support treatment planning at state-of-the-art ion beam therapy facilities.
Ion beam therapy is a rapidly growing technique for tumor radiation therapy. Ions allow for a high dose deposition in the tumor region, while sparing the surrounding healthy tissue. For this reason, ...the highest possible accuracy in the calculation of dose and its spatial distribution is required in treatment planning. On one hand, commonly used treatment planning software solutions adopt a simplified beam-body interaction model by remapping pre-calculated dose distributions into a 3D water-equivalent representation of the patient morphology. On the other hand, Monte Carlo (MC) simulations, which explicitly take into account all the details in the interaction of particles with human tissues, are considered to be the most reliable tool to address the complexity of mixed field irradiation in a heterogeneous environment. However, full MC calculations are not routinely used in clinical practice because they typically demand substantial computational resources. Therefore MC simulations are usually only used to check treatment plans for a restricted number of difficult cases. The advent of general-purpose programming GPU cards prompted the development of trimmed-down MC-based dose engines which can significantly reduce the time needed to recalculate a treatment plan with respect to standard MC codes in CPU hardware. In this work, we report on the development of fred, a new MC simulation platform for treatment planning in ion beam therapy. The code can transport particles through a 3D voxel grid using a class II MC algorithm. Both primary and secondary particles are tracked and their energy deposition is scored along the trajectory. Effective models for particle-medium interaction have been implemented, balancing accuracy in dose deposition with computational cost. Currently, the most refined module is the transport of proton beams in water: single pencil beam dose-depth distributions obtained with fred agree with those produced by standard MC codes within 1-2% of the Bragg peak in the therapeutic energy range. A comparison with measurements taken at the CNAO treatment center shows that the lateral dose tails are reproduced within 2% in the field size factor test up to 20 cm. The tracing kernel can run on GPU hardware, achieving 10 million primary s−1 on a single card. This performance allows one to recalculate a proton treatment plan at 1% of the total particles in just a few minutes.
. The Monte Carlo simulation software is a valuable tool in radiation therapy, in particular to achieve the needed accuracy in the dose evaluation for the treatment plans optimisation. The current ...challenge in this field is the time reduction to open the way to many clinical applications for which the computational time is an issue. In this manuscript we present an innovative GPU-accelerated Monte Carlo software for dose valuation in electron and photon based radiotherapy, developed as an update of the FRED (Fast paRticle thErapy Dose evaluator) software.
. The code transports particles through a 3D voxel grid, while scoring their energy deposition along their trajectory. The models of electromagnetic interactions in the energy region between 1 MeV-1 GeV available in literature have been implemented to efficiently run on GPUs, allowing to combine a fast tracking while keeping high accuracy in dose assessment. The FRED software has been bench-marked against state-of-art full MC (FLUKA, GEANT4) in the realm of two different radiotherapy applications: Intra-Operative Radio Therapy and Very High Electron Energy radiotherapy applications.
. The single pencil beam dose-depth profiles in water as well as the dose map computed on non-homogeneous phantom agree with full-MCs at 2% level, observing a gain in processing time from 200 to 5000.
. Such performance allows for computing a plan with electron beams in few minutes with an accuracy of ∼%, demonstrating the FRED potential to be adopted for fast plan re-calculation in photon or electron radiotherapy applications.
Nuclear fragmentation measurements are necessary when using heavy-ion beams in hadrontherapy to predict the effects of the ion nuclear interactions within the human body. Moreover, they are also ...fundamental to validate and improve the Monte Carlo codes for their use in planning tumor treatments. Nowadays, a very limited set of carbon fragmentation cross sections are being measured, and in particular, to our knowledge, no double-differential fragmentation cross sections at intermediate energies are available in the literature. In this work, we have measured the double-differential cross sections and the angular distributions of the secondary fragments produced in the 12C fragmentation at 62 A MeV on a thin carbon target. The experimental data have been used to benchmark the prediction capability of the Geant4 Monte Carlo code at intermediate energies, where it was never tested before. In particular, we have compared the experimental data with the predictions of two Geant4 nuclear reaction models: the Binary Light Ions Cascade and the Quantum Molecular Dynamic. From the comparison, it has been observed that the Binary Light Ions Cascade approximates the angular distributions of the fragment production cross sections better than the Quantum Molecular Dynamic model. However, the discrepancies observed between the experimental data and the Monte Carlo simulations lead to the conclusion that the prediction capability of both models needs to be improved at intermediate energies.
Tumour control is performed in particle therapy using particles and ions, whose high irradiation precision enhances the effectiveness of the treatment, while sparing the healthy tissue surrounding ...the target volume. Dose range monitoring devices using photons and charged particles produced by the beam interacting with the patient's body have already been proposed, but no attempt has been made yet to exploit the detection of the abundant neutron component. Since neutrons can release a significant dose far away from the tumour region, precise measurements of their flux, production energy and angle distributions are eagerly sought in order to improve the treatment planning system (TPS) software. It will thus be possible to predict not only the normal tissue toxicity in the target region, but also the risk of late complications in the whole body. The aforementioned issues underline the importance of an experimental effort devoted to the precise characterisation of neutron production, aimed at the measurement of their abundance, emission point and production energy. The technical challenges posed by a neutron detector aimed at high detection efficiency and good backtracking precision are addressed within the MONDO (monitor for neutron dose in hadrontherapy) project, whose main goal is to develop a tracking detector that can target fast and ultrafast neutrons. A full reconstruction of two consecutive elastic scattering interactions undergone by the neutrons inside the detector material will be used to measure their energy and direction. The preliminary results of an MC simulation performed using the FLUKA software are presented here, together with the DSiPM (digital SiPM) readout implementation. New detector readout implementations specifically tailored to the MONDO tracker are also discussed, and the neutron detection efficiency attainable with the proposed neutron tracking strategy are reported.
The calculation algorithm of a modern treatment planning system for ion-beam radiotherapy should ideally be able to deal with different ion species (e.g. protons and carbon ions), to provide relative ...biological effectiveness (RBE) evaluations and to describe different beam lines. In this work we propose a new approach for ion irradiation outcomes computations, the beamlet superposition (BS) model, which satisfies these requirements. This model applies and extends the concepts of previous fluence-weighted pencil-beam algorithms to quantities of radiobiological interest other than dose, i.e. RBE- and LET-related quantities. It describes an ion beam through a beam-line specific, weighted superposition of universal beamlets. The universal physical and radiobiological irradiation effect of the beamlets on a representative set of water-like tissues is evaluated once, coupling the per-track information derived from FLUKA Monte Carlo simulations with the radiobiological effectiveness provided by the microdosimetric kinetic model and the local effect model. Thanks to an extension of the superposition concept, the beamlet irradiation action superposition is applicable for the evaluation of dose, RBE and LET distributions. The weight function for the beamlets superposition is derived from the beam phase space density at the patient entrance. A general beam model commissioning procedure is proposed, which has successfully been tested on the CNAO beam line. The BS model provides the evaluation of different irradiation quantities for different ions, the adaptability permitted by weight functions and the evaluation speed of analitical approaches. Benchmarking plans in simple geometries and clinical plans are shown to demonstrate the model capabilities.
Uncertainties in determining clinically used relative biological effectiveness (RBE) values for ion beam therapy carry the risk of absolute and relative misestimations of RBE-weighted doses for ...clinical scenarios. This study assesses the consequences of hypothetical misestimations of input parameters to the RBE modelling for carbon ion treatment plans by a variational approach. The impact of the variations on resulting cell survival and RBE values is evaluated as a function of the remaining ion range. In addition, the sensitivity to misestimations in RBE modelling is compared for single fields and two opposed fields using differing optimization criteria. It is demonstrated for single treatment fields that moderate variations (up to ±50%) of representative nominal input parameters for four tumours result mainly in a misestimation of the RBE-weighted dose in the planning target volume (PTV) by a constant factor and only smaller RBE-weighted dose gradients. Ensuring a more uniform radiation quality in the PTV eases the clinical importance of uncertainties in the radiobiological treatment parameters, as for such a condition uncertainties tend to result only in a systematic misestimation of RBE-weighted dose in the PTV by a constant factor. Two opposed carbon ion fields with a constant RBE in the PTV are found to result in rather robust conditions. Treatments using two ion species may be used to achieve a constant RBE in the PTV irrespective of the size and depth of the spread-out Bragg peak.
Radiation therapy is currently the most utilized technique for the treatment of tumors by means of ionizing radiation, such as electrons, protons and x/gamma rays, depending on the type, size and ...depth of the cancer mass. Radiation therapy has in general fulfilled the main requirement of targeting thus damaging the malignant cells and sparing the healthy tissues as best as possible. In this scenario, electron linear accelerators have been operated as viable tools for the delivery of both high-energetic electrons and x-ray beams, which are obtained via the bremsstrahlung process of the electrons hitting on a high-Z material. Recently, it has been experimentally demonstrated that ultrahigh dose-rate bursts of electrons and x-ray beams increase the differential response between healthy and tumor tissues. This beneficial response is referred to as the FLASH effect. For this purpose, we have developed the first dedicated compactS-band linear accelerator for FLASH radiotherapy. This linac is optimized for a nominal energy of 7 MeV and a pulsed electron beam current of 100 mA and above. The accelerator is mounted on a remote-controlled system for preclinical research studies in the FLASH regime. We will show the rf and beam dynamics design of theS-band linac as well as the commissioning and high-power rf tests. Furthermore, the results of the dosimetric measurements will be illustrated.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UL, UM, UPUK