Definition of dose rate for FLASH pencil-beam scanning proton therapy: A comparative study

Purpose: Highlight the distinctions, both in terms of concept and numerical values, of the various definitions that can be established for the dose rate in PBS proton therapy. Methods: In an in silico study, five definitions of the dose rate, namely the PBS dose rate, the percentile dose rate, the maximum percentile dose rate, the average dose rate, and the dose averaged dose rate (DADR) were analyzed first through theoretical comparison, and then applied to a head and neck case. To carry out this study, a treatment plan utilizing a single energy level and requiring the use of a patient-specific range modulator was employed. The dose rate values were compared both locally and by means of dose rate volume histograms (DRVHs). Results: The PBS dose rate, the percentile dose rate, and the maximum percentile dose are definitions that are specifically designed to take into account the time structure of the delivery of a PBS treatment plan. Although they may appear similar, our study shows that they can vary locally by up to 10%. On the other hand, the DADR values were approximately twice as high as those of the PBS, percentile and maximum percentile dose rates, since the DADR disregards the periods when a voxel does not receive any dose. Finally, the average dose rate can be defined in various ways, as discussed in this paper. The average dose rate is found to be lower by a factor of approximately 1/2 than the PBS, percentile and maximum percentile dose rates. Conclusions: We have shown that using different definitions for the dose rate in FLASH proton therapy can lead to variations in calculated values ranging from a few percent to a factor of two. Since the dose rate is a critical parameter in FLASH radiation therapy, it is essential to carefully consider the choice of definition. However, to make an informed decision, additional biological data and models are needed

Integrate range shifter in immobilization for proton therapy: 3D printed materials characterization

3D printing is investigated for application in patient immobilization during proton therapy (PT). It potentially enables a merge of immobilization, range shifting and other functionality into one patient-specific structure. Beside minimizing the lateral beam spread due to the removal of air gap it could also reduce the collision risk and the treatment time compared to movable nozzle snouts. In this first study, 9 different 3D printed materials were characterized in detail. The resulting data (Table 1) will serve as input for the design of a printed immobilization structure. The printed test objects showed reduced geometric printing accuracy for 3 materials. Compression testing yielded Young’s moduli from 0.6 MPa to 3445 MPa, without deterioration after exposure to 100 Gy in a MV photon beam. Dual-energy CT methods were used to estimate the effective atomic number Zeff, the relative electron density e and the stopping power ratio SPR. Zeff ranged from 5.91 to 10.43. The SPR and e both ranged from 0.6 to 1.22. The measured photon attenuation coefficients at therapeutic energies scaled linearly with e. In a 62 MeV proton beam, good agreement was seen between the DECT estimated SPR and the measured range shift, except for the higher Zeff. As opposed to the photon attenuation, the proton range shifting was printing orientation dependent for certain materials. In conclusion printed materials exhibit a wide variation in structural and radiological properties. The quantification of these characteristics enables optimal material selection for the design of a multifunctional 3D printed immobilization structure for PT

Three-Dimensional Dose Prediction for Lung IMRT Patients with Deep Neural Networks: Robust Learning from Heterogeneous Beam Configurations

The use of neural networks to directly predict three-dimensional dose distributions for automatic planning is becoming popular. However, the existing methods only use patient anatomy as input and assume consistent beam configuration for all patients in the training database. The purpose of this work is to develop a more general model that, in addition to patient anatomy, also considers variable beam configurations, to achieve a more comprehensive automatic planning with a potentially easier clinical implementation, without the need of training specific models for different beam settings

Towards 3D printed multifunctional immobilization for proton therapy: initial materials characterization

Purpose: 3D printing technology is investigated for the purpose of patient immobilization during proton therapy. It potentially enables a merge of patient immobilization, bolus range shifting, and other functions into one single patient-speci c structure. In this rst step, a set of 3D printed materials is characterized in detail, in terms of structural and radiological properties, elemental composition, directional dependence, and structural changes induced by radiation damage. These data will serve as inputs for the design of 3D printed immobilization structure prototypes. Methods: Using four di erent 3D printing techniques, in total eight materials were subjected to testing. Samples with a nominal dimension of 20×20×80 mm3 were 3D printed. The geometrical printing accuracy of each test sample was measured with a dial gage. To assess the mechanical response of the samples, standardized compression tests were performed to determine the Young’s modulus. To investigate the e ect of radiation on the mechanical response, the mechanical tests were performed both prior and after the administration of clinically relevant dose levels (70 Gy), multiplied with a safety factor of 1.4. Dual energy computed tomography (DECT) methods were used to calculate the relative electron density to water ρe, the e ective atomic number Ze , and the proton stopping power ratio (SPR) to water SPR. In order to validate the DECT based calculation of radiological properties, beam measurements were performed on the 3D printed samples as well. Photon irradiations were performed to measure the photon linear attenuation coe cients, while proton irradiations were performed to measure the proton range shift of the samples. The direc- tional dependence of these properties was investigated by performing the irradiations for di erent orientations of the samples. Results: The printed test objects showed reduced geometric printing accuracy for 2 materials (deviation > 0.25 mm). Compression tests yielded Young’s moduli ranging from 0.6 to 2940 MPa. No deterioration in the mechanical response was observed after exposure of the samples to 100 Gy in a therapeutic MV photon beam. The DECT-based characterization yielded Ze ranging from 5.91 to 10.43. The SPR and ρe both ranged from 0.6 to 1.22. The measured photon attenuation coe cients at clinical energies scaled linearly with ρe. Good agreement was seen between the DECT estimated SPR and the measured range shift, except for the higher Ze . As opposed to the photon attenuation, the proton range shifting appeared to be printing orientation dependent for certain materials. Conclusions: In this study, the rst step toward 3D printed, multifunctional immobilization was performed, by going through a candidate clinical work ow for the rst time: from the material printing to DECT characterization with a veri cation through beam measurements. Besides a proof of concept for beam modi cation, the mechanical response of printed materials was also investigated to assess their capabilities for positioning functionality. For the studied set of printing techniques and materials, a wide variety of mechanical and radiological properties can be selected from for the intended purpose. Moreover the elaborated hybrid DECT methods aid in performing in-house quality assurance of 3D printed components, as these methods enable the estimation of the radiological properties relevant for use in radiation therapy

A hybrid multi-particle approach to range assessment-based treatment verification in particle therapy

Particle therapy (PT) used for cancer treatment can spare healthy tissue and reduce treatment toxicity. However, full exploitation of the dosimetric advantages of PT is not yet possible due to range uncertainties, warranting development of range-monitoring techniques. This study proposes a novel range-monitoring technique introducing the yet unexplored concept of simultaneous detection and imaging of fast neutrons and prompt-gamma rays produced in beam-tissue interactions. A quasimonolithic organic detector array is proposed, and its feasibility for detecting range shifts in the context of proton therapy is explored through Monte Carlo simulations of realistic patient models and detector resolution efects. The results indicate that range shifts of 1 mm can be detected at relatively low proton intensities (22.30(13) × 107 protons/spot) when spatial information obtained through imaging of both particle species are used simultaneously. This study lays the foundation for multiparticle detection and imaging systems in the context of range verifcation in PTpublishedVersio

Application of GAFchromic films for tomotherapy

n

Potential pitfalls of the PTV concept in dose-to-medium planning optimization

In typical treatment planning of 3D IMRT, the incident energy fluence is optimized to achieve a homoge- neous dose distribution to the PTV. The PTV includes the tumour but also healthy tissues that may have a different dose response for the same incident energy fluence, like bony structures included in the PTV (mandibles in head and neck tumours or femoral bones in sarcomas). Dose to medium optimization com- pensates for this heterogeneous response, leading to a non-homogeneous energy fluence in the PTV and a non-homogeneous dose in the CTV in the presence of geometric errors. We illustrate qualitatively this statement in a cylindrical geometry where the PTV includes a CTV (7 cm diameter) made of water sur- rounded by ICRU compact bone (1.2 cm thickness); such configuration was chosen to exaggerate the aforementioned effect. Optimization was performed assuming dose equals photon energy fluence times mass energy absorption coefficient. Bone has a 4% lower dose response in a 6 MV flattening filter free spectrum. After optimization either in medium or assuming everything as water composition, the geom- etry was shifted by 1.2 cm and dose recomputed. As expected, compensating for the under-response of the bone material during optimization in medium leads to an overdosage of the CTV when patient geo- metric errors are taken into account. Optimization in dose assuming everything as water composition leads to a uniform coverage. Robust optimization or forcing a uniform atomic composition in the PTV margin may resolve this incompatibility between the PTV concept and dose to medium optimization

Monte Carlo evaluation of the dose calculation algorithm of TomoTherapy for clinical cases in dynamic jaws mode

Purpose: For the TomoTherapy® system, longitudinal conformation can be improved by selecting a smaller field width but at the expense of longer treatment time. Recently, the TomoEdge® feature has been released with the possibility to move dynamically the jaws at the edges of the target volume, improving longitudinal penumbra and enabling faster treatments. Such delivery scheme requires additional modeling of treatment delivery. Using a previously validated Monte Carlo model (TomoPen), we evaluated the accuracy of the implementation of TomoEdge in the new dose engine of TomoTherapy for 15 clinical cases. Methods: TomoPen is based on PENELOPE. Particle tracking in the treatment head is performed almost instantaneously by 1) reading a particle from a phase-space file corresponding to the largest field and 2) correcting the weight of the particle depending on the actual jaw and MLC configurations using Monte Carlo pre-generated data. 15 clinical plans (5 head-and-neck, 5 lung and 5 prostate tumors) planned with TomoEdge and with the last release of the treatment planning system (VoLO®) were re-computed with TomoPen. The resulting dose-volume histograms were compared. Results: Good agreement was achieved overall, with deviations for the target volumes typically within 2% (D95), excepted for small lung tumors (17 cm3) where a maximum deviation of 4.4% was observed for D95. The results were consistent with previously reported values for static field widths. Conclusions: For the clinical cases considered in the present study, the introduction of TomoEdge did not impact significantly the accuracy of the computed dose distributions

Monte Carlo-based simulation of dynamic jaws tomotherapy

Purpose: Original TomoTherapy systems may involve a trade-off between conformity and treatment speed, the user being limited to three slice widths (1.0, 2.5, and 5.0 cm). This could be overcome by allowing the jaws to define arbitrary fields, including very small slice widths (<1 cm), which are challenging for a beam model. The aim of this work was to incorporate the dynamic jaws feature into a Monte Carlo (MC) model called TomoPen, based on the MC code PENELOPE, previously validated for the original TomoTherapy system. Methods: To keep the general structure of TomoPen and its efficiency, the simulation strategy introduces several techniques: (1) weight modifiers to account for any jaw settings using only the 5 cm phase-space file; (2) a simplified MC based model called FastStatic to compute the modifiers faster than pure MC; (3) actual simulation of dynamic jaws. Weight modifiers computed with both FastStatic and pure MC were compared. Dynamic jaws simulations were compared with the convolution/superposition (C/S) of TomoTherapy in the “cheese” phantom for a plan with two targets longitudinally separated by a gap of 3 cm. Optimization was performed in two modes: asymmetric jaws-constant couch speed (“running start stop,” RSS) and symmetric jaws-variable couch speed (“symmetric running start stop,” SRSS). Measurements with EDR2 films were also performed for RSS for the formal validation of TomoPen with dynamic jaws. Results: Weight modifiers computed with FastStatic were equivalent to pure MC within statistical uncertainties (0.5% for three standard deviations). Excellent agreement was achieved between TomoPen and C/S for both asymmetric jaw opening/constant couch speed and symmetric jaw opening/variable couch speed, with deviations well within 2%/2 mm. For RSS procedure, agreement between C/S and measurements was within 2%/2 mm for 95% of the points and 3%/3 mm for 98% of the points, where dose is greater than 30% of the prescription dose (gamma analysis). Dose profiles acquired in transverse and longitudinal directions through the center of the phantom were also compared with excellent agreement (2%/2 mm) between all modalities. Conclusions: The combination of weights modifiers and interpolation allowed implementing efficiently dynamic jaws and dynamic couch features into TomoPen at a minimal cost in terms of efficiency (simulation around 8 h on a single CPU)