Investigation of dosimetric effect of beam current fluctuations in synchrotron-based proton PBS continuous scanning.

Objective.In proton pencil beam scanning (PBS) continuous delivery, the beam is continuously delivered without interruptions between spots. For synchrotron-based systems, the extracted beam current exhibits a spill structure, and recent publications on beam current measurements have demonstrated significant fluctuations around the nominal values. These fluctuations potentially lead to dose deviations from those calculated assuming a stable beam current. This study investigated the dosimetric implications of such beam current fluctuations during proton PBS continuous scanning.Approach.Using representative clinical proton PBS plans, we performed simulations to mimic a worst-case clinical delivery environment with beam current varies from 50% to 250% of the nominal values. The simulations used the beam delivery parameters optimized for the best beam delivery efficiency of the upcoming particle therapy system at Mayo Clinic Florida. We reconstructed the simulated delivered dose distributions and evaluated the dosimetric impact of beam current fluctuations.Main results.Despite significant beam current fluctuations resulting in deviations at each spot level, the overall dose distributions were nearly identical to those assuming a stable beam current. The 1 mm/1% Gamma passing rate was 100% for all plans. Less than 0.2% root mean square error was observed in the planning target volume dose-volume histogram. Minimal differences were observed in all dosimetric evaluation metrics.Significance.Our findings demonstrate that with our beam delivery system and clinical planning practice, while significant beam current fluctuations may result in large local move monitor unit deviations at each spot level, the overall impact on the dose distribution is minimal.

Assessment of fluence- and dose-averaged linear energy transfer with passive luminescence detectors in clinical proton beams.

Objective.This work investigates the use of passive luminescence detectors to determine different types of averaged linear energy transfer (LET-) for the energies relevant to proton therapy. The experimental results are compared to reference values obtained from Monte Carlo simulations.Approach.Optically stimulated luminescence detectors (OSLDs), fluorescent nuclear track detectors (FNTDs), and two different groups of thermoluminescence detectors (TLDs) were irradiated at four different radiation qualities. For each irradiation, the fluence- (LET-f) and dose-averaged LET (LET-d) were determined. For both quantities, two sub-types of averages were calculated, either considering the contributions from primary and secondary protons or from all protons and heavier, charged particles. Both simulated and experimental data were used in combination with a phenomenological model to estimate the relative biological effectiveness (RBE).Main results.All types ofLET-could be assessed with the luminescence detectors. The experimental determination ofLET-fis in agreement with reference data obtained from simulations across all measurement techniques and types of averaging. On the other hand,LET-dcan present challenges as a radiation quality metric to describe the detector response in mixed particle fields. However, excluding secondaries heavier than protons from theLET-dcalculation, as their contribution to the luminescence is suppressed by ionization quenching, leads to equal accuracy betweenLET-fandLET-d. Assessment of RBE through the experimentally determinedLET-dvalues agrees with independently acquired reference values, indicating that the investigated detectors can determineLET-with sufficient accuracy for proton therapy.Significance.OSLDs, TLDs, and FNTDs can be used to determineLET-and RBE in proton therapy. With the capability to determine dose through ionization quenching corrections derived fromLET-, OSLDs and TLDs can simultaneously ascertain dose,LET-, and RBE. This makes passive detectors appealing for measurements in phantoms to facilitate validation of clinical treatment plans or experiments related to proton therapy.

Development and characterization of the first proton minibeam system for single-gantry proton facility.

BACKGROUND: Minibeam represents a preclinical spatially fractionated radiotherapy modality with great translational potential. The advantage lies in its high therapeutic index (compared to GRID and LATTICE) and ability to treat at greater depth (compared to microbeam). Proton minibeam radiotherapy (pMBRT) is a synergy of proton and minibeam. While the single-gantry proton facility has gained popularity due to its affordability and compact design, it often has limited beam time available for research purposes. Conversely, given the current requirement of pMBRT on specific minibeam hardware collimators, necessitates a reproducible and fast setup to minimize pMBRT treatment time and streamline the switching time between pMBRT and conventional treatment for clinically translation. PURPOSE: The contribution of this work is the development and characterization of the first pMBRT system tailored for single-gantry proton facility. The system allows for efficient and reproducible plug-and-play setup, achievable within minutes. METHODS: The single room pMBRT system is constructed based on IBA ProteusONE proton machine. The end of nozzle is attached with beam modifying accessories though an accessory drawer. A small snout is attached to the accessory drawer and used to hold apertures and range shifters. The minibeam aperture consists of two components: a fitting ring and an aperture body. Three minibeam apertures were manufactured. The first-generation apertures underwent qualitatively analysis with film, and the second generation aperture underwent more comprehensive quantitative measurement. The reproducibility of the setup is accessed, and the film measurements are performed to characterize the pMBRT system in cross validation with Monte Carlo (MC) simulations. RESULTS: We presented initial results of large field pMBRT aperture and the film measurements indicates the effect of source-to-isocenter distance = 930 cm in Y proton scanning direction. Consistent with TOPAS MC simulation, the dose uniformity of pMBRT field <2 cm is demonstrated to be better than 2%, rendering its suitability for pre-clinical studies. Subsequently, we developed the second generation of aperture with five slits and characterized the aperture with film dosimetry studies and compared the results to the benchmark MC. Comprehensive film measurements were also performed to evaluate the effect of divergence, air gap and gantry-angle dependency and repeatability and revealing a consistent performance within 5%. Furthermore, the 2D gamma analysis indicated a passing rate exceeding 99% using 3% dose difference and 0.2 mm distance agreement criteria. We also establish the peak valley dose ratio and the depth dose profile measurements, and the results are within 10% from MC simulation. CONCLUSIONS: We have developed the first pMBRT system tailored for a single-gantry proton facility, which has demonstrated accuracy in benchmark with MC simulations, and allows for efficient plug-and-play setup, emphasizing efficiency.

Experimental comparison of cylindrical and plane parallel ionization chambers for reference dosimetry in continuous and pulsed scanned proton beams.

Objective. In this experimental work we compared the determination of absorbed dose to water using four ionization chambers (ICs), a PTW-34045 Advanced Markus, a PTW-34001 Roos, an IBA-PPC05 and a PTW-30012 Farmer, irradiated under the same conditions in one continuous- and in two pulsed-scanned proton beams.Approach. The ICs were positioned at 2 cm depth in a water phantom in four square-field single-energy scanned-proton beams with nominal energies between 80 and 220 MeV and in the middle of 10 × 10 × 10 cm3dose cubes centered at 10 cm or 12.5 cm depth in water. The water-equivalent thickness (WET) of the entrance window and the effective point of measurement was considered when positioning the plane parallel (PP) ICs and the cylindrical ICs, respectively. To reduce uncertainties, all ICs were calibrated at the same primary standards laboratory. We used the beam quality (kQ) correction factors for the ICs under investigation from IAEA TRS-398, the newly calculated Monte Carlo (MC) values and the anticipated IAEA TRS-398 updated recommendations.Main results. Dose differences among the four ICs ranged between 1.5% and 3.7% using both the TRS-398 and the newly recommendedkQvalues. The spread among the chambers is reduced with the newlykQvalues. The largest differences were observed between the rest of the ICs and the IBA-PPC05 IC, obtaining lower dose with the IBA-PPC05.Significance. We provide experimental data comparing different types of chambers in different proton beam qualities. The observed dose differences between the ICs appear to be related to inconsistencies in the determination of thekQvalues. For PP ICs, MC studies account for the physical thickness of the entrance window rather than the WET. The additional energy loss that the wall material invokes is not negligible for the IBA-PPC05 and might partially explain the lowkQvalues determined for this IC. To resolve this inconsistency and to benchmark MC values,kQvalues measured using calorimetry are needed.

Future technological developments in proton therapy - A predicted technological breakthrough.

In the current spectrum of cancer treatments, despite high costs, a lack of robust evidence based on clinical outcomes or technical and radiobiological uncertainties, particle therapy and in particular proton therapy (PT) is rapidly growing. Despite proton therapy being more than fifty years old (first proposed by Wilson in 1946) and more than 220,000 patients having been treated with in 2020, many technological challenges remain and numerous new technical developments that must be integrated into existing systems. This article presents an overview of on-going technical developments and innovations that we felt were most important today, as well as those that have the potential to significantly shape the future of proton therapy. Indeed, efforts have been done continuously to improve the efficiency of a PT system, in terms of cost, technology and delivery technics, and a number of different developments pursued in the accelerator field will first be presented. Significant developments are also underway in terms of transport and spatial resolution achievable with pencil beam scanning, or conformation of the dose to the target: we will therefore discuss beam focusing and collimation issues which are important parameters for the development of these techniques, as well as proton arc therapy. State of the art and alternative approaches to adaptive PT and the future of adaptive PT will finally be reviewed. Through these overviews, we will finally see how advances in these different areas will allow the potential for robust dose shaping in proton therapy to be maximised, probably foreshadowing a future era of maturity for the PT technique.

Accounting for prompt gamma emission and detection for range verification in proton therapy treatment planning.

Prompt gamma (PG) imaging is widely investigated as one of the most promising methods for proton range verification in proton therapy. The performance of this technique is affected by several factors like tissue heterogeneity, number of protons in the considered pencil beam and the detection device. Our previous work proposed a new treatment planning concept which boosts the number of protons of a few PG monitoring-friendly pencil beams (PBs), selected on the basis of two proposed indicators quantifying the conformity between the dose and PG at the emission level, above the desired detectability threshold. To further explore this method at the detection level, in this work we investigated the response of a knife-edge slit PG camera which was deployed in the first clinical application of PG to proton therapy monitoring. The REGistration Graphical User Interface (REGGUI) is employed to simulate the PG emission, PG detection as well as the corresponding dose distribution. As the PG signal detected by this kind of PG camera is sensitive to the relative position of the camera and PG signal falloff, we optimized our PB selection method for this camera by introducing a new camera position indicator identifying whether the expected falloff of the PG signal is centered in the field of view of the camera or not. Our camera-adapted PB selection method is investigated using computed tomography (CT) scans at two different treatment time points of a head and neck, and a prostate cancer patient under scenarios considering different statistics level. The results show that a precision of 0.8 mm for PG falloff identification can be achieved when a PB has more than 2 × 108 primary protons. Except for one case due to unpredictable and comparably large anatomical changes, the PG signals of most of the PBs recommended by all our indicators are observed to be reliable for proton range verification with deviations between the inter-fractional shift of proton range (as deduced from the PB dose distribution) and the detected PG signal within 2.0 mm. In contrast, a shift difference up to 9.6 mm has been observed for the rejected PBs. The magnitude of the proton range shift due to the inter-fractional anatomical changes is observed to be up to 23 mm. The proposed indicators are shown to be valuable for identifying and recommending reliable PBs to create new PG monitoring-friendly TPs. Comparison between our PB boosting method and the alternative PB aggregation, which combines the signal of nearby PBs to reach the desired counting statistics, is also discussed.

A novel proton counting detector and method for the validation of tissue and implant material maps for Monte Carlo dose calculation.

The presence of artificial implants complicates the delivery of proton therapy due to inaccurate characterization of both the implant and the surrounding tissues. In this work, we describe a method to characterize implant and human tissue mimicking materials in terms of relative stopping power (RSP) using a novel proton counting detector. Each proton is tracked by directly measuring the deposited energy along the proton track using a fast, pixelated spectral detector AdvaPIX-TPX3 (TPX3). We considered three scenarios to characterize the RSPs. First, in-air measurements were made in the presence of metal rods (Al, Ti and CoCr) and bone. Then, measurements of energy perturbations in the presence of metal implants and bone in an anthropomorphic phantom were performed. Finally, sampling of cumulative stopping power (CSP) of the phantom were made at different locations of the anthropomorphic phantom. CSP and RSP information were extracted from energy spectra at each beam path. To quantify the RSP of metal rods we used the shift in the most probable energy (MPE) of CSP from the reference CSP without a rod. Overall, the RSPs were determined as 1.48, 2.06, 3.08, and 5.53 from in-air measurements; 1.44, 1.97, 2.98, and 5.44 from in-phantom measurements, for bone, Al, Ti and CoCr, respectively. Additionally, we sampled CSP for multiple paths of the anthropomorphic phantom ranging from 18.63 to 25.23 cm deriving RSP of soft tissues and bones in agreement within 1.6% of TOPAS simulations. Using minimum error of these multiple CSP, optimal mass densities were derived for soft tissue and bone and they are within 1% of vendor-provided nominal densities. The preliminary data obtained indicates the proposed novel method can be used for the validation of material and density maps, required by proton Monte Carlo Dose calculation, provided by competing multi-energy computed tomography and metal artifact reduction techniques.

Out-of-field doses for scanning proton radiotherapy of shallowly located paediatric tumours-a comparison of range shifter and 3D printed compensator.

The lowest possible energy of proton scanning beam in cyclotron proton therapy facilities is typically between 60 and 100 MeV. Treatment of superficial lesions requires a pre-absorber to deliver doses to shallower volumes. In most of the cases a range shifter (RS) is used, but as an alternative solution, a patient-specific 3D printed proton beam compensator (BC) can be applied. A BC enables further reduction of the air gap and consequently reduction of beam scattering. Such pre-absorbers are additional sources of secondary radiation. The aim of this work was the comparison of RS and BC with respect to out-of-field doses for a simulated treatment of superficial paediatric brain tumours. EURADOS WG9 performed comparative measurements of scattered radiation in the Proteus C-235 IBA facility (Cyclotron Centre Bronowice at the Institute of Nuclear Physics, CCB IFJ PAN, Kraków, Poland) using two anthropomorphic phantoms-5 and 10 yr old-for a superficial target in the brain. Both active detectors located inside the therapy room, and passive detectors placed inside the phantoms were used. Measurements were supplemented by Monte Carlo simulation of the radiation transport. For the applied 3D printed pre-absorbers, out-of-field doses from both secondary photons and neutrons were lower than for RS. Measurements with active environmental dosimeters at five positions inside the therapy room indicated that the RS/BC ratio of the out-of-field dose was also higher than one, with a maximum of 1.7. Photon dose inside phantoms leads to higher out-of-field doses for RS than BC to almost all organs with the highest RS/BC ratio 12.5 and 13.2 for breasts for 5 and 10 yr old phantoms, respectively. For organs closest to the isocentre such as the thyroid, neutron doses were lower for BC than RS due to neutrons moderation in the target volume, but for more distant organs like bladder-conversely-lower doses for RS than BC were observed. The use of 3D printed BC as the pre-absorber placed in the near vicinity of patient in the treatment of superficial tumours does not result in the increase of secondary radiation compared to the treatment with RS, placed far from the patient.

On-line range verification for proton beam therapy using spherical ionoacoustic waves with resonant frequency.

In contrast to conventional X-ray therapy, proton beam therapy (PBT) can confine radiation doses to tumours because of the presence of the Bragg peak. However, the precision of the treatment is currently limited by the uncertainty in the beam range. Recently, a unique range verification methodology has been proposed based on simulation studies that exploit spherical ionoacoustic waves with resonant frequency (SPIREs). SPIREs are emitted from spherical gold markers in tumours initially introduced for accurate patient positioning when the proton beam is injected. These waves have a remarkable property: their amplitude is linearly correlated with the residual beam range at the marker position. Here, we present proof-of-principle experiments using short-pulsed proton beams at the clinical dose to demonstrate the feasibility of using SPIREs for beam-range verification with submillimetre accuracy. These results should substantially contribute to reducing the range uncertainty in future PBT applications.

A new detector for the beam energy measurement in proton therapy: a feasibility study.

The proof of concept of a new device, capable of determining in a few seconds the energy of clinical proton beams by measuring the time of flight (ToF) of protons, is presented. The prototype consists of two thin ultra fast silicon detector (UFSD) pads, aligned along the beam direction in a telescope configuration and readout by a digitizer. The method developed for extracting the energy at the isocenter from the measured ToF, validated by Monte Carlo simulations, and the procedure used to calibrate the system are also presented and discussed in detail. The prototype was tested at the Centro Nazionale di Adroterapia Oncologica (CNAO, Pavia, Italy), at several beam energies, covering the entire clinical range, and using different distances between the sensors. The measured beam energies were benchmarked against the nominal CNAO energy values, obtained during the commissioning of the centre from the measured ranges in water. Deviations of few hundreds of keV have been achieved for all considered proton beam energies for distances between the two sensors larger than 60 cm, indicating a sensitivity to the corresponding beam range in water smaller than the clinical tolerance of 1 mm. Moreover, few seconds of irradiation were necessary to collect the required statistics. These preliminary results indicate that a telescope of UFSDs could achieve in a short time the accuracy required for the clinical application and therefore encourage further investigations towards the improvement and the optimization of the present prototype.

Preparation of a radiobiology beam line at the 18 MeV proton cyclotron facility at CNA.

Proton therapy has gained interest in recent years due to its excellent clinical outcomes. However, the lack of accurate biological data, especially in the Bragg peak region of clinical beams, makes it difficult to implement biophysically optimized treatment plans in clinical practice. In this context, low energy proton accelerator facilities provide the perfect environment to collect good radiobiological data, as they can produce high LET beams with narrow energy distributions. This study presents the radiobiology beam line that has been designed at the 18 MeV proton cyclotron facility at the National Centre of Accelerators (CNA, Seville, Spain), to perform irradiations of mono-layer cell cultures. To ensure that all the cells receive the same dose with a suitable dose rate, low beam intensities and broad and homogeneous beam profiles are necessary. To do so, at the CNA an unfocused beam has been used, broadened with a 500 µm thick aluminium scattering foil. Homogeneous dose profiles, with deviations lower than 10% have been obtained over a circular surface of 35 mm diameter for an incident average energy of 12.8 MeV. Further, a Monte Carlo simulation of the beam line has been developed with Geant4, and benchmarked towards experimental measurements, with differences generally below 1%. Once validated, the code has been used, together with an ionization chamber, for dosimetry studies, to characterize the beam and monitor the dose. Finally, cultures of Human Bone Osteosarcoma cells (U2OS) have been successfully irradiated at the radiobiology beam line, investigating the effects of radiation in terms of DNA damage induction.

Prompt gamma spectroscopy for absolute range verification of 12C ions at synchrotron-based facilities.

The physical range uncertainty limits the exploitation of the full potential of charged particle therapy. In this work, we face this issue aiming to measure the absolute Bragg peak position in the target. We investigate p, 4He, 12C and 16O beams accelerated at the Heidelberg Ion-Beam Therapy Center. The residual range of the primary 12C ions is correlated to the energy spectrum of the prompt gamma radiation. The prompt gamma spectroscopy method was demonstrated for proton beams accelerated by cyclotrons and is developed here for the first time for heavier ions accelerated by a synchrotron. We develop a detector system that includes (i) a spectroscopic unit based on cerium(III) bromide and bismuth germanium oxide scintillating crystals, (ii) a beam trigger based on an array of scintillating fibers and (iii) a data acquisition system based on a FlashADC. We test the system in two different scenarios. In the first series of experiments, we detect and identify 19 independent spectral lines over a wide gamma energy spectrum in the presence of the four ion species for different targets, including a water target with a titanium insert. In the second series of experiments, we introduce a collimator aiming to relate the spectral information to the range of the primary particles. We perform extensive measurements for a 12C beam and demonstrate submillimetric precision for the measurement of its Bragg peak position in the experimental setup. The features of the energy and time spectra for gamma radiation induced by p, 4He and 16O are investigated upstream and downstream from the Bragg peak position. We conclude the analysis by extrapolating the required future developments, which would be needed to achieve range verification with a 2 mm accuracy during a single fraction delivery of [Formula: see text] physical dose.

Sensitivity analysis of Monte Carlo model of a gantry-mounted passively scattered proton system.

PURPOSE: This study aimed to present guidance on the correlation between treatment nozzle and proton source parameters, and dose distribution of a passive double scattering compact proton therapy unit, known as Mevion S250. METHODS: All 24 beam options were modeled using the MCNPX MC code. The calculated physical dose for pristine peak, profiles, and spread out Bragg peak (SOBP) were benchmarked with the measured data. Track-averaged LET (LETt ) and dose-averaged LET (LETd ) distributions were also calculated. For the sensitivity investigations, proton beam line parameters including Average Energy (AE), Energy Spread (ES), Spot Size (SS), Beam Angle (BA), Beam Offset (OA), and Second scatter Offset (SO) from central Axis, and also First Scatter (FS) thickness were simulated in different stages to obtain the uncertainty of the derived results on the physical dose and LET distribution in a water phantom. RESULTS: For the physical dose distribution, the MCNPX MC model matched measurements data for all the options to within 2 mm and 2% criterion. The Mevion S250 was found to have a LETt between 0.46 and 8.76 keV.µm-1 and a corresponding LETd between 0.84 and 15.91 keV.µm-1 . For all the options, the AE and ES had the greatest effect on the resulting depth of pristine peak and peak-to-plateau ratio respectively. BA, OA, and SO significantly decreased the flatness and symmetry of the profiles. The LETs were found to be sensitive to the AE, ES, and SS, especially in the peak region. CONCLUSIONS: This study revealed the importance of considering detailed beam parameters, and identifying those that resulted in large effects on the physical dose distribution and LETs for a compact proton therapy machine.

Design and verification of an external radiobiological beam port on a 16.5 MeV GE PETtrace proton cyclotron.

PURPOSE: Protons and heavy ions are considered to be ideal particles for use in external beam radiotherapy due to the superior properties of the dose distribution. While a photon (x-ray) beam delivers considerable dose to healthy tissues around the tumor, a proton beam that is delivered with sufficient energies has: a low entrance dose (the dose in front of the tumor); a high-dose region within the tumor, known as the Bragg peak; and, no exit dose beyond the tumor. Proton therapy is the next major step in advancing radiotherapy treatment. The purpose of this project was to adapt an existing radioisotope production cyclotron, a General Electric (GE) PETtrace, to enable radiobiological studies using proton beams. During routine use the PETtrace delivers 16.5 MeV protons to target with beam currents in the range of 10-100 µA resulting in dose rates in the order of kGy/s. To achieve the aim of the project the dose rate had to be reduced to the Gy/min range, without attenuating the proton energy below 5 MeV. This paper covers the design, construction and validation of the beam port. METHODS: Monte Carlo simulations were performed, using GEANT4, SRIM and PACE4 to design the beam port and optimize its components. Once the beam port was fabricated, validation experiments were performed using EBT3 and HD-V2 Gafchromic™ films, and a Keithley 6485 picoampere meter. RESULTS AND CONCLUSION: The external beam port was successfully modeled, designed and fabricated. By using a 0.25 mm thick gold foil and a brass pin-hole collimator the beam was spread from a narrow full beam diameter of 10 mm to a wide beam with a 5% flatness area in the center of the beam that had a diameter of ~20 mm. In using this system the dose rate was reduced from kGy/s to ~30 Gy/min.

Verification of a Monte Carlo dose calculation engine in proton minibeam radiotherapy in a passive scattering beamline for preclinical trials.

OBJECTIVES: Proton minibeam radiation therapy (pMBRT) is a novel therapeutic strategy that combines the benefits of proton therapy with the remarkable normal tissue preservation observed with the use of submillimetric spatially fractionated beams. This promising technique has been implemented at the Institut Curie-Proton therapy centre (ICPO) using a first prototype of a multislit collimator. The purpose of this work was to develop a Monte Carlo-based dose calculation engine to reliably guide preclinical studies at ICPO. METHODS: The whole "Y1"-passive beamline at the ICPO, including pMBRT implementation, was modelled using the Monte Carlo GATE v. 7.0 code. A clinically relevant proton energy (100 MeV) was used as starting point. Minibeam generation by means of the brass collimator used in the first experiments was modelled. A virtual source was modelled at the exit of the beamline nozzle and outcomes were compared with dosimetric measurements performed with EBT3 gafchromic films and a diamond detector in water. Dose distributions were recorded in a water phantom and in rat CT images (7-week-old male Fischer rats). RESULTS: The dose calculation engine was benchmarked against experimental data and was then used to assess dose distributions in CT images of a rat, resulting from different irradiation configurations used in several experiments. It reduced computational time by an order of magnitude. This allows us to speed up simulations for in vivo trials, where we obtained peak-to-valley dose ratios of 1.20 ± 0.05 and 6.1 ± 0.2 for proton minibeam irradiations targeting the tumour and crossing the rat head. Tumour eradication was observed in the 67 and 22% of the animals treated respectively. CONCLUSION: A Monte Carlo dose calculation engine for pMBRT implementation with mechanical collimation has been developed. This tool can be used to guide and interpret the results of in vivo trials. ADVANCES IN KNOWLEDGE: This is the first Monte Carlo dose engine for pMBRT that is being used to guide preclinical trials in a clinical proton therapy centre.

Transitioning from measurement-based to combined patient-specific quality assurance for intensity-modulated proton therapy.

OBJECTIVE: This study is part of ongoing efforts aiming to transit from measurement-based to combined patient-specific quality assurance (PSQA) in intensity-modulated proton therapy (IMPT). A Monte Carlo (MC) dose-calculation algorithm is used to improve the independent dose calculation and to reveal the beam modeling deficiency of the analytical pencil beam (PB) algorithm. METHODS: A set of representative clinical IMPT plans with suboptimal PSQA results were reviewed. Verification plans were recalculated using an MC algorithm developed in-house. Agreements of PB and MC calculations with measurements that quantified by the γ passing rate were compared. RESULTS: The percentage of dose planes that met the clinical criteria for PSQA (>90% γ passing rate using 3%/3 mm criteria) increased from 71.40% in the original PB calculation to 95.14% in the MC recalculation. For fields without beam modifiers, nearly 100% of the dose planes exceeded the 95% γ passing rate threshold using the MC algorithm. The model deficiencies of the PB algorithm were found in the proximal and distal regions of the SOBP, where MC recalculation improved the γ passing rate by 11.27% (p < 0.001) and 16.80% (p < 0.001), respectively. CONCLUSIONS: The MC algorithm substantially improved the γ passing rate for IMPT PSQA. Improved modeling of beam modifiers would enable the use of the MC algorithm for independent dose calculation, completely replacing additional depth measurements in IMPT PSQA program. For current users of the PB algorithm, further improving the long-tail modeling or using MC simulation to generate the dose correction factor is necessary. ADVANCES IN KNOWLEDGE: We justified a change in clinical practice to achieve efficient combined PSQA in IMPT by using the MC algorithm that was experimentally validated in almost all the clinical scenarios in our center. Deficiencies in beam modeling of the current PB algorithm were identified and solutions to improve its dose-calculation accuracy were provided.

Simulation of spread-out bragg peaks in proton beams using Geant4/TOPAS.

The simulation of proton Spread-Out Bragg Peaks (SOBPs) was implemented using the Geant4-based TOPAS Monte Carlo software. Dynamic proton energy switching was implemented using TOPAS time features, while beam weights were calculated using an empirical power law formalism with Bragg peaks spaced by 0.5 mm. To find power parameters yielding flat SOBPs we sampled power parameters for maximum kinetic energies of 50 MeV to 250 MeV and SOBP widths of 15% to 40% of the depth of the distal SOBP end. Simulations were run in a 50 cm cubic water phantom using a uniform squared proton beam. Depth dose was scored along the central axis in a binned cylinder with 1 cm diameter in 2.5 mm increments. Power parameters yielding a flat SOBPs were found to vary with, both energy and SOBP width and differed significantly from previously reported values based on simulations with MCNPX.

Shaping of a laser-accelerated proton beam for radiobiology applications via genetic algorithm.

Laser-accelerated protons have a great potential for innovative experiments in radiation biology due to the sub-picosecond pulse duration and high dose rate achievable. However, the broad angular divergence makes them not optimal for applications with stringent requirements on dose homogeneity and total flux at the irradiated target. The strategy otherwise adopted to increase the homogeneity is to increase the distance between the source and the irradiation plane or to spread the beam with flat scattering systems or through the transport system itself. Such methods considerably reduce the proton flux and are not optimal for laser-accelerated protons. In this paper we demonstrate the use of a Genetic Algorithm (GA) to design an optimal non-flat scattering system to shape the beam and efficiently flatten the transversal dose distribution at the irradiated target. The system is placed in the magnetic transport system to take advantage of the presence of chromatic focusing elements to further mix the proton trajectories. The effect of a flat scattering system placed after the transport system is also presented for comparison. The general structure of the GA and its application to the shaping of a laser-accelerated proton beam are presented, as well as its application to the optimisation of dose distribution in a water target in air.

[Proton therapy in France in 2019]. / État des lieux de la protonthérapie en France en 2019.

Among over 100 proton therapy centres worldwide in operation or under construction, French proton therapy is coming to full maturity with the recent opening of the Nice (1991, upgrade in 2016) and Caen (2018) facilities next to the Orsay (1991, upgrade in 2010) centre. Proton therapy is a national priority for children and young adults in all three centres. The patient-related activity of the three French centres is coordinated via the Protonshare portal to optimise referral by type of indication and available expertise in coordination with the French society of radiation oncology SFRO and French radiotherapy centres. The centres are recognised by the French Health Care excellence initiative, promoted by the ministry of Foreign Affairs. The three centres collaborate structurally in terms of clinical research and are engaged at the international level in the participation to European databases and research initiatives. Concerted actions are now also promoted in preclinical research via the Radiotransnet network. Ongoing French developments in proton therapy are well presented in international hadron therapy meetings, including European Proton Therapy Network and Particle Therapy Cooperative Oncology Group. Proton therapy teaching in France is offered at several levels and is open to colleagues from all radiation oncology centres, so that they are fully informed, involved and trained to facility recognition of possible indications and thereby to contribute to appropriate patient referral. This close collaboration between all actors in French radiation oncology facilitates the work to demonstrate the required level of medical and scientific evidence for current and emerging indications for particle therapy. Based on that, the future might entail a possible creation of more proton therapy facilities in France.

Polynomial modelling of proton trajectories in homogeneous media for fast most likely path estimation and trajectory simulation.

Protons undergo many small angle deflections when traversing a medium, such as a patient. This effect, known as multiple Coulomb scattering (MCS), leads to degraded image resolution in proton radiography and computed tomography (CT) and to lateral spreading of the dose distribution in proton therapy. To optimally account for MCS in proton imaging, the most likely path (MLP) of a proton is estimated based on its position and propagation angle measured in front of and behind the object. In this work, we propose a functional which quantifies the likelihood of a proton trajectory and study how it can be used to model proton trajectories in a homogeneous medium. We focus on two aspects: first, we present an analytical method to quickly generate proton trajectories in a homogeneous medium based on the likelihood functional and validate it through Monte Carlo simulations. It could be used for fast generation of proton CT images without a full Monte Carlo simulation, or potentially to complement the components in a treatment planning Monte Carlo which simulate MCS. Second, by maximising the likelihood functional, we derive an expression for the MLP which is equivalent to the conventional ones reported in the literature yet computationally more convenient. Moreover, we show that the MLP is strictly a polynomial function if the protons' energy loss in the medium is approximated as a polynomial and that the orders of both are linked. We validate our MLP through Monte Carlo simulations and compare proton CT images reconstructed with our expression and with the conventional one. We find that an MLP polynomial of orders larger than five do not lead to increased spatial resolution compared to lower order expressions.