Perspectives in linear accelerator for FLASH VHEE: Study of a compact C-band system

PURPOSE

In order to translate the FLASH effect in clinical use and to treat deep tumors, Very High Electron Energy irradiations could represent a valid technique. Here, we address the main issues in the design of a VHEE FLASH machine. We present preliminary results for a compact C-band system aiming to reach a high accelerating gradient and high current necessary to deliver a Ultra High Dose Rate with a beam pulse duration of 3μs.

METHODS

The proposed system is composed by low energy high current injector linac followed by a high acceleration gradient structure able to reach 60-160 MeV energy range. To obtain the maximum energy, an energy pulse compressor options is considered. CST code was used to define the specifications RF parameters of the linac. To optimize the accelerated current and therefore the delivered dose, beam dynamics simulations was performed using TSTEP and ASTRA codes.

RESULTS

The VHEE parameters Linac suitable to satisfy FLASH criteria were simulated. Preliminary results allow to obtain a maximum energy of 160 MeV, with a peak current of 200 mA, which corresponds to a charge of 600 nC.

CONCLUSIONS

A promising preliminary design of VHEE linac for FLASH RT has been performed. Supplementary studies are on going to complete the characterization of the machine and to manufacture and test the RF prototypes.

4He dose- and track-averaged linear energy transfer: Monte Carlo algorithms and experimental verification

Objective. In the present hadrontherapy scenario, there is a growing interest in exploring the capabilities of different ion species other than protons and carbons. The possibility of using different ions paves the way for new radiotherapy approaches, such as the multi-ions treatment, where radiation could vary according to target volume, shape, depth and histologic characteristics of the tumor. For these reasons, in this paper, the study and understanding of biological-relevant quantities was extended for the case of 4He ion. Approach. Geant4 Monte Carlo based algorithms for dose- and track-averaged LET (Linear Energy Transfer) calculations, were validated for 4He ions and for the case of a mixed field characterised by the presence of secondary ions from both target and projectile fragmentation. The simulated dose and track averaged LETs were compared with the corresponding dose and frequency mean values of the lineal energy, y D ¯ and y ¯ F , derived from experimental microdosimetric spectra. Two microdosimetric experimental campaigns were carried out at the Italian eye proton therapy facility of the Laboratori Nazionali del Sud of Istituto Nazionale di Fisica Nucleare (INFN-LNS, Catania, I) using two different microdosimeters: the MicroPlus probe and the nano-TEPC (Tissue Equivalent Proportional Counter). Main results. A good agreement of L ¯ d Total and L ¯ t Total with y ¯ D and y ¯ T experimentally measured with both microdosimetric detectors MicroPlus and nano-TEPC in two configurations: full energy and modulated 4He ion beam, was found. Significance. The results of this study certify the use of a very effective tool for the precise calculation of LET, given by a Monte Carlo approach which has the advantage of allowing detailed simulation and tracking of nuclear interactions, even in complex clinical scenarios.

Experimental investigation at CATANA facility of n-10B and p-11B reactions for the enhancement of proton therapy

The aim of the NEPTUNE (Nuclear process-driven Enhancement of Proton Therapy UNravEled) project is to investigate in detail both the physical and radiobiological phenomena that could justify an increase of the proton-induced cytogenetic effects in cells irradiated in presence of an agent containing natural boron. In this work, a double-stage silicon telescope coupled to different boron converters was irradiated at the CATANA proton therapy facility (INFN-LNS) for studying the proton boron fusion and the neutron boron capture reactions by discriminating secondary particles from primary protons. Different boron targets were developed by depositing boric acid, enriched with a higher than 99% content of ¹⁰B or ¹¹B, on a 50 µm thick PolyMethilMetacrylate (PMMA) substrate. The ¹⁰B target allows to evaluate the contribution of lithium and alpha particles produced by the boron neutron capture reaction triggered by secondary thermal neutrons, while the ¹¹B target is exploited for studying the effect of the p + ¹¹B → 3α nuclear reaction directly triggered by primary protons. Experimental results clearly show the presence of alpha particles from both the reactions. The silicon telescope is capable of discriminating, by means of the so-called "scatter plots", the contribution of alpha particles originated by thermal neutrons on ¹⁰B with respect to the ones produced by protons impinging on ¹¹B. Although a reliable quantitative study of the alpha production rate has not been achieved yet, this work demonstrates that low energy and, therefore, high-LET particles from both the reactions can be measured.

Report on G4-Med, a Geant4 benchmarking system for medical physics applications developed by the Geant4 Medical Simulation Benchmarking Group

BACKGROUND

Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro- and nanodosimetry, imaging, radiation protection, and nuclear medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing.

AIMS

To respond to these needs, we developed G4-Med, a benchmarking and regression testing system of Geant4 for medical physics.

MATERIALS AND METHODS

G4-Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4-Med are executed on the CERN computing infrastructure via the use of the geant-val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes.

RESULTS

This paper describes the tests included in G4-Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data.

DISCUSSION

Our results indicate that the Geant4 electromagnetic physics constructor G4EmStandardPhysics_option4 gives a good agreement with the reference data for all the tests. The QGSP_BIC_HP physics list provided an overall adequate description of the physics involved in hadron therapy, including proton and carbon ion therapy. New tests should be included in the next stage of the project to extend the benchmarking to other physical quantities and application scenarios of interest for medical physics.

CONCLUSION

The results presented and discussed in this paper will aid users in tailoring physics lists to their particular application.

Microdosimetry of a therapeutic proton beam with a mini-TEPC and a MicroPlus-Bridge detector for RBE assessment

Proton beams are widely used worldwide to treat localized tumours, the lower entrance dose and no exit dose, thus sparing surrounding normal tissues, being the main advantage of this treatment modality compared to conventional photon techniques. Clinical proton beam therapy treatment planning is based on the use of a general relative biological effectiveness (RBE) of 1.1 along the whole beam penetration depth, without taking into account the documented increase in RBE at the end of the depth dose profile, in the Bragg peak and beyond. However, an inaccurate estimation of the RBE can cause both underdose or overdose, in particular it can cause the unfavourable situation of underdosing the tumour and overdosing the normal tissue just beyond the tumour, which limits the treatment success and increases the risk of complications. In view of a more precise dose delivery that takes into account the variation of RBE, experimental microdosimetry offers valuable tools for the quality assurance of LET or RBE-based treatment planning systems. The purpose of this work is to compare the response of two different microdosimetry systems: the mini-TEPC and the MicroPlus-Bridge detector. Microdosimetric spectra were measured across the 62 MeV spread out Bragg peak of CATANA with the mini-TEPC and with the Bridge microdosimeter. The frequency and dose distributions of lineal energy were compared and the different contributions to the spectra were analysed, discussing the effects of different site sizes and chord length distributions. The shape of the lineal energy distributions measured with the two detectors are markedly different, due to the different water-equivalent sizes of the sensitive volumes: 0.85 μm for the TEPC and 17.3 μm for the silicon detector. When the Loncol's biological weighting function is applied to calculate the microdosimetric assessment of the RBE, both detectors lead to results that are consistent with biological survival data for glioma U87 cells. Both the mini-TEPC and the MicroPlus-Bridge detector can be used to assess the RBE variation of a 62 MeV modulated proton beam along its penetration depth. The microdosimetric assessment of the RBE based on the Loncol's weighting function is in good agreement with radiobiological results when the 10% biological uncertainty is taken into account.

A new energy spectrum reconstruction method for time-of-flight diagnostics of high-energy laser-driven protons

The Time-of-Flight (TOF) technique coupled with semiconductorlike detectors, as silicon carbide and diamond, is one of the most promising diagnostic methods for high-energy, high repetition rate, laser-accelerated ions allowing a full on-line beam spectral characterization. A new analysis method for reconstructing the energy spectrum of high-energy laser-driven ion beams from TOF signals is hereby presented and discussed. The proposed method takes into account the detector's working principle, through the accurate calculation of the energy loss in the detector active layer, using Monte Carlo simulations. The analysis method was validated against well-established diagnostics, such as the Thomson parabola spectrometer, during an experimental campaign carried out at the Rutherford Appleton Laboratory (UK) with the high-energy laser-driven protons accelerated by the VULCAN Petawatt laser.

Radiobiological quantities in proton-therapy: Estimation and validation using Geant4-based Monte Carlo simulations

PURPOSE

The Geant4 Monte Carlo simulation toolkit was used to reproduce radiobiological parameters measured by irradiating three different cancerous cell lines with monochromatic and clinical proton beams.

METHODS

The experimental set-up adopted for irradiations was fully simulated with a dedicated open-source Geant4 application. Cells survival fractions was calculated coupling the Geant4 simulations with two analytical radiobiological models: one based on the LEM (Local Effect Model) approach and the other on a semi-empirical parameterisation. Results was evaluated and compared with experimental data.

RESULTS AND CONCLUSIONS

The results demonstrated the Geant4 ability to reproduce radiobiological quantities for different cell lines.

Transversal dose distribution optimization for laser-accelerated proton beam medical applications by means of Geant4

The main purpose of this paper is to quantitatively study the possibility of delivering dose distributions of clinical relevance with laser-driven proton beams. A Monte Carlo application has been developed with the Geant4 toolkit, simulating the ELIMED (MEDical and multidisciplinary application at ELI-Beamlines) transport and dosimetry beam line which is being currently installed at the ELI-Beamlines in Prague (CZ). The beam line will be used to perform irradiations for multidisciplinary studies, with the purpose of demonstrating the possible use of optically accelerated ion beams for therapeutic purposes. The ELIMED Geant4-based application, already validated against reference transport codes, accurately simulates each single element of the beam line, necessary to collect the accelerated beams and to select them in energy. Transversal dose distributions at the irradiation point have been studied and optimized to try to quantitatively answer the question if such kind of beam lines, and specifically the systems developed for ELIMED in Prague, will be actually able to transport ion beams not only for multidisciplinary applications, such as pitcher-catcher nuclear reactions (e.g. neutrons), PIXE analysis for cultural heritage and space radiation, but also for delivering dose patterns of clinical relevance in a future perspective of possible medical applications.

First experimental proof of Proton Boron Capture Therapy (PBCT) to enhance protontherapy effectiveness

Protontherapy is hadrontherapy's fastest-growing modality and a pillar in the battle against cancer. Hadrontherapy's superiority lies in its inverted depth-dose profile, hence tumour-confined irradiation. Protons, however, lack distinct radiobiological advantages over photons or electrons. Higher LET (Linear Energy Transfer) 12C-ions can overcome cancer radioresistance: DNA lesion complexity increases with LET, resulting in efficient cell killing, i.e. higher Relative Biological Effectiveness (RBE). However, economic and radiobiological issues hamper 12C-ion clinical amenability. Thus, enhancing proton RBE is desirable. To this end, we exploited the p + 11B → 3α reaction to generate high-LET alpha particles with a clinical proton beam. To maximize the reaction rate, we used sodium borocaptate (BSH) with natural boron content. Boron-Neutron Capture Therapy (BNCT) uses 10B-enriched BSH for neutron irradiation-triggered alpha particles. We recorded significantly increased cellular lethality and chromosome aberration complexity. A strategy combining protontherapy's ballistic precision with the higher RBE promised by BNCT and 12C-ion therapy is thus demonstrated.

Continuous monitoring of noise levels in the Gulf of Catania (Ionian Sea). Study of correlation with ship traffic

Acoustic noise levels were measured in the Gulf of Catania (Ionian Sea) from July 2012 to May 2013 by a low frequency (<1000Hz) hydrophone, installed on board the NEMO-SN1 multidisciplinary observatory. NEMO-SN1 is a cabled node of EMSO-ERIC, which was deployed at a water depth of 2100m, 25km off Catania. The study area is characterized by the proximity of mid-size harbors and shipping lanes. Measured noise levels were correlated with the passage of ships tracked with a dedicated AIS antenna. Noise power was measured in the frequency range between 10Hz and 1000Hz. Experimental data were compared with the results of a fast numerical model based on AIS data to evaluate the contribution of shipping noise in six consecutive 1/3 octave frequency bands, including the 1/3 octave frequency bands centered at 63Hz and 125Hz, indicated by the Marine Strategy Framework Directive (2008/56/EC).

Nuclear reaction measurements on tissue-equivalent materials and GEANT4 Monte Carlo simulations for hadrontherapy

When a carbon beam interacts with human tissues, many secondary fragments are produced into the tumor region and the surrounding healthy tissues. Therefore, in hadrontherapy precise dose calculations require Monte Carlo tools equipped with complex nuclear reaction models. To get realistic predictions, however, simulation codes must be validated against experimental results; the wider the dataset is, the more the models are finely tuned.Since no fragmentation data for tissue-equivalent materials at Fermi energies are available in literature, we measured secondary fragments produced by the interaction of a 55.6 MeV u(-1) (12)C beam with thick muscle and cortical bone targets. Three reaction models used by the Geant4 Monte Carlo code, the Binary Light Ions Cascade, the Quantum Molecular Dynamic and the Liege Intranuclear Cascade, have been benchmarked against the collected data. In this work we present the experimental results and we discuss the predictive power of the above mentioned models.

Reference dosimetry for light-ion beams based on graphite calorimetry

Developments in hadron therapy require efforts to improve the accuracy of the dose delivered to a target volume. Here, the determination of the absorbed dose under reference conditions was analysed. Based on the International Atomic Energy Agency TRS-398 code of practice, for hadron beams, the combined standard uncertainty on absorbed dose to water under reference conditions, derived from ionisation chambers, is too large. This uncertainty is dominated by the beam quality correction factors, [Formula: see text], mainly due to the mean energy to produce one ion pair in air, wair. A method to reduce this uncertainty is to carry out primary dosimetry, using calorimetry. A [Formula: see text]-value can be derived from a direct comparison between calorimetry and ionometry. Here, this comparison is performed using a graphite calorimeter in an 80-MeV A(-1) carbon ion beam. Assuming recommended TRS-398 values of water-to-graphite stopping power ratio and the perturbation factor for an ionisation chamber, preliminary results indicate a wair-value of 35.5 ± 0.9 J C(-1).

Carbon fragmentation measurements and validation of the Geant4 nuclear reaction models for hadrontherapy

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 (12)C 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.

Charged particle's flux measurement from PMMA irradiated by 80 MeV/u carbon ion beam

Hadrontherapy is an emerging technique in cancer therapy that uses beams of charged particles. To meet the improved capability of hadrontherapy in matching the dose release with the cancer position, new dose-monitoring techniques need to be developed and introduced into clinical use. The measurement of the fluxes of the secondary particles produced by the hadron beam is of fundamental importance in the design of any dose-monitoring device and is eagerly needed to tune Monte Carlo simulations. We report the measurements carried out with charged secondary particles produced from the interaction of a 80 MeV/u fully stripped carbon ion beam at the INFN Laboratori Nazionali del Sud, Catania, with a poly-methyl methacrylate target. Charged secondary particles, produced at 90° with respect to the beam axis, have been tracked with a drift chamber, while their energy and time of flight have been measured by means of a LYSO scintillator. Secondary protons have been identified exploiting the energy and time-of-flight information, and their emission region has been reconstructed backtracking from the drift chamber to the target. Moreover, a position scan of the target indicates that the reconstructed emission region follows the movement of the expected Bragg peak position. Exploiting the reconstruction of the emission region, an accuracy on the Bragg peak determination in the submillimeter range has been obtained. The measured differential production rate for protons produced with E(Prod)(kin) > 83 MeV and emitted at 90° with respect to the beam line is dN(P)/(dN(C)dΩ) (E(Prod)(kin) > 83 MeV, θ = 90°) = (2.69 ± 0.08(stat) ± 0.12(sys)) × 10⁻⁴ sr⁻¹.

Response of human HTB140 melanoma cells to conventional radiation and hadrons

Conventional radiotherapy with X- and gamma-rays is one of the common and effective treatments of cancer. High energy hadrons, i.e., charged particles like protons and (12)C ions, due to their specific physics and radiobiological advantages are increasingly used. In this study, effectiveness of different radiation types is evaluated on the radio-resistant human HTB140 melanoma cells. The cells were irradiated with gamma-rays, the 62 MeV protons at the Bragg peak and in the middle of the spread-out Bragg peak (SOBP), as well as with the 62 MeV/u (12)C ions. The doses ranged from 2 to 24 Gy. Cell survival and proliferation were assessed 7 days after irradiation, whereas apoptosis was evaluated after 48 h. The acquired results confirmed the high radio-resistance of cells, showing better effectiveness of protons than gamma-rays. The best efficiency was obtained with (12)C ions due to higher linear energy transfer. All analyzed radiation qualities reduced cell proliferation. The highest proliferation was detected for (12)C ions because of their large killing capacity followed by small induction of reparable lesions. This enabled unharmed cells to preserve proliferative activity. Irradiations with protons and (12)C ions revealed similar moderate pro-apoptotic ability that is in agreement with the level of cellular radio-resistance.

Cellular and molecular effects of protons: apoptosis induction and potential implications for cancer therapy

Due to their ballistic precision, apoptosis induction by protons could be a strategy to specifically eliminate neoplastic cells. To characterize the cellular and molecular effects of these hadrons, we performed dose-response and time-course experiments by exposing different cell lines (PC3, Ca301D, MCF7) to increasing doses of protons and examining them with FACS, RT-PCR, and electron spin resonance (ESR). Irradiation with a dose of 10 Gy of a 26,7 Mev proton beam altered cell structures such as membranes, caused DNA double strand breaks, and significantly increased intracellular levels of hydroxyl ions, are active oxygen species (ROS). This modified the transcriptome of irradiated cells, activated the mitochondrial (intrinsic) pathway of apoptosis, and resulted in cycle arrest at the G2/M boundary. The number of necrotic cells within the irradiated cell population did not significantly increase with respect to the controls. The effects of irradiation with 20 Gy were qualitatively as well as quantitatively similar, but exposure to 40 Gy caused massive necrosis. Similar experiments with photons demonstrated that they induce apoptosis in a significantly lower number of cells and in a temporally delayed manner. These data advance our knowledge on the cellular and molecular effects of proton irradiation and could be useful for improving current hadrontherapy protocols.

Two-dimensional dosimetry of radiotherapeutical proton beams using thermoluminescence foils

In modern radiation therapy such as intensity modulated radiation therapy or proton therapy, one is able to cover the target volume with improved dose conformation and to spare surrounding tissue with help of modern measurement techniques. Novel thermoluminescence dosimetry (TLD) foils, developed from the hot-pressed mixture of LiF:Mg,Cu,P (MCP TL) powder and ethylene-tetrafluoroethylene (ETFE) copolymer, have been applied for 2-D dosimetry of radiotherapeutical proton beams at INFN Catania and IFJ Krakow. A TLD reader with 70 mm heating plate and CCD camera was used to read the 2-D emission pattern of irradiated foils. The absorbed dose profiles were evaluated, taking into account correction factors specific for TLD such as dose and energy response. TLD foils were applied for measuring of dose distributions within an eye phantom and compared with predictions obtained from the MCNPX code and Eclipse Ocular Proton Planning (Varian Medical Systems) clinical radiotherapy planning system. We demonstrate the possibility of measuring 2-D dose distributions with point resolution of about 0.5 x 0.5 mm(2).

DNA fragments induction in human fibroblasts by radiations of different qualities

Experimental data on DNA double strand break (DSB) induction in human fibroblasts (AG1522), following irradiation with several radiation qualities, namely gamma rays, 0.84 MeV protons, 58.9 MeV u(-1) carbon ions, iron ions of 115 MeV u(-1), 414 MeV u(-1), 1 GeV u(-1), and 5 GeV u(-1), are presented. DSB yields were measured by calibrated Pulsed Field Gel Electrophoresis in the DNA fragment size range 0.023-5.7 Mbp. The DSB yields show little LET dependence, in spite of the large variation of the latter among the beams, and are slightly higher than that obtained using gamma rays. The highest yield was found for the 5 GeV u(-1) iron beam, that gave a value 30% higher than the 1 GeV u(-1) iron beam. A phenomenological method is used to parametrise deviation from randomness in fragment size spectra.