A comprehensive Monte Carlo study of out-of-field secondary neutron spectra in a scanned-beam proton therapy gantry room.

PURPOSE: To simulate secondary neutron radiation fields that had been measured at different relative positions during phantom irradiation inside a scanning proton therapy gantry treatment room. Further, to identify origin, energy distribution, and angular emission of the secondary neutrons as a function of proton beam energy. METHODS: The FLUKA Monte Carlo code was used to model the relevant parts of the treatment room in a scanned pencil beam proton therapy gantry including shielding walls, floor, major metallic gantry-components, patient table, and a homogeneous PMMA target. The proton beams were modeled based on experimental beam ranges in water and spot shapes in air. Neutron energy spectra were simulated at 0°, 45°, 90° and 135° relative to the beam axis at 2m distance from isocenter for monoenergetic 11×11cm2 fields from 200MeV, 140MeV, 75MeV initial proton beams, as well as for 118MeV protons with a 5cm thick PMMA range shifter. The total neutron spectra were scored for these four positions and proton energies. FLUKA neutron spectra simulations were crosschecked with Geant4 simulations using initial proton beam properties from FLUKA-generated phase spaces. Additionally, the room-components generating secondary neutrons in the room and their contributions to the total spectrum were identified and quantified. RESULTS: FLUKA and Geant4 simulated neutron spectra showed good general agreement with published measurements in the whole simulated neutron energy range of 10-10 to 103MeV. As in previous studies, high-energy (E≥19.6MeV) neutrons from the phantom are most prevalent along 0°, while thermalized (1meV≤E<0.4eV) and fast (100keV≤E<19.4MeV) neutrons dominate the spectra in the lateral and backscatter direction. The iron of the large bending magnet and its counterweight mounted on the gantry were identified as the most determinant sources of secondary fast-neutrons, which have been lacking in simplified room simulations. CONCLUSIONS: The results helped disentangle the origin of secondary neutrons and their dominant contributions and were strengthened by the fact that a cross comparison was made using two independent Monte Carlo codes. The complexity of such room model can in future be limited using the result. They may further be generalized in that they can be used for an assessment of neutron fields, possibly even at facilities where detailed neutron measurements and simulations cannot be performed. They may also help to design future proton therapy facilities and to reduce unwanted radiation doses from secondary neutrons to patients.

Radiation protection modelling for 2.5 Petawatt-laser production of ultrashort x-ray, proton and ion bunches: Monte Carlo model of the Munich CALA facility.

The 'Centre for Advanced Laser Applications' (CALA) is a new research institute for laser-based acceleration of electron beams for brilliant x-ray generation, laser-driven sub-nanosecond bunches of protons and heavy ions for biomedical applications like imaging and tumour therapy as well as for nuclear and high-field physics.The radiation sources emerging from experiments using the up to 2.5 petawatt laser pulses with 25 femtosecond duration will be mixed particle-species of high intensity, high energy and pulsed, thus posing new challenges compared to conventional radiation protection. Such worldwide pioneering laser experiments result in source characteristics that require careful a-priori radiation safety simulations.The FLUKA Monte-Carlo code was used to model the five CALA experimental caves, including the corridors, halls and air spaces surrounding the caves. Beams of electrons (<5 GeV), protons (<200 MeV),12C (<400MeV/u) and197Au (<10MeV/u) ions were simulated using spectra, divergences and bunch-charges based on expectations from recent scientific progress.Simulated dose rates locally can exceed 1.5 kSv h-1inside beam dumps. Vacuum pipes in the cave walls for laser transport and extraction channels for the generated x-rays result in small dose leakage to neighboring areas. Secondary neutrons contribute to most of the prompt dose rate outside caves into which the beam is delivered. This secondary radiation component causes non-negligible dose rates to occur behind walls to which large fluences of secondary particles are directed.By employing adequate beam dumps matched to beam-divergence, magnets, passive shielding and laser pulse repetition limits, average dose rates in- and outside the experimental building stay below design specifications (<0.5µSv h-1) for unclassified areas,<2.5µSv h-1for supervised areas,<7.5µSv h-1maximum local dose rate) and regulatory limits (<1mSv a-1for unclassified areas).

A feasibility study of zebrafish embryo irradiation with laser-accelerated protons.

The development from single shot basic laser plasma interaction research toward experiments in which repetition rated laser-driven ion sources can be applied requires technological improvements. For example, in the case of radio-biological experiments, irradiation duration and reproducible controlled conditions are important for performing studies with a large number of samples. We present important technological advancements of recent years at the ATLAS 300 laser in Garching near Munich since our last radiation biology experiment. Improvements range from target positioning over proton transport and diagnostics to specimen handling. Exemplarily, we show the current capabilities by performing an application oriented experiment employing the zebrafish embryo model as a living vertebrate organism for laser-driven proton irradiation. The size, intensity, and energy of the laser-driven proton bunches resulted in evaluable partial body changes in the small (<1 mm) embryos, confirming the feasibility of the experimental system. The outcomes of this first study show both the appropriateness of the current capabilities and the required improvements of our laser-driven proton source for in vivo biological experiments, in particular the need for accurate, spatially resolved single bunch dosimetry and image guidance.

Optimization and performance study of a proton CT system for pre-clinical small animal imaging.

Proton computed tomography (pCT) promises to reduce or even eliminate range uncertainties inherent in the conversion of Hounsfield units into relative stopping power (RSP) for proton therapy treatment planning. This is of particular interest for proton irradiation studies in animal models due to the high precision required and uncertainties in tissue properties. We propose a dedicated single-particle tracking pCT system consisting of low material budget floating strip Micromegas detectors for tracking and a segmented time-projection-chamber with vertical Mylar absorbers, functioning as a range telescope. Based on Monte Carlo simulations of a realistic in silico beam and detector implementation, a geometrical optimization of the system components was conducted to safeguard an ideal operation close to intrinsic performance limits at 75 MeV. Moreover, the overall imaging capabilities relevant for pre-clinical proton therapy treatment planning were evaluated for a mouse model. In order to minimize extrinsic uncertainties in the estimated proton trajectories, a spacing of the two tracking planes of at least 7 cm is required in both tracking detectors. Additionally, novel in-house developed and produced aluminum-based readout electrodes promise superior performance with around 3 mm-1 spatial resolution due to the reduced material budget. Concerning the range telescope, an absorber thickness within 500 µm to 750 µm was found to yield the best compromise between water-equivalent path length resolution and complexity of the detector instrumentation, still providing sub-0.5% RSP accuracy. The optimized detector configuration enables better than 2% range accuracy for proton therapy treatment planning in pre-clinical data sets. This work outlines the potential of pCT for small animal imaging. The performance of the proposed and optimized system provides superior treatment planning accuracy compared to conventional x-ray CT. Thus, pCT can play an important role in translational and pre-clinical cancer research.

I-BEAT: Ultrasonic method for online measurement of the energy distribution of a single ion bunch.

The shape of a wave carries all information about the spatial and temporal structure of its source, given that the medium and its properties are known. Most modern imaging methods seek to utilize this nature of waves originating from Huygens' principle. We discuss the retrieval of the complete kinetic energy distribution from the acoustic trace that is recorded when a short ion bunch deposits its energy in water. This novel method, which we refer to as Ion-Bunch Energy Acoustic Tracing (I-BEAT), is a refinement of the ionoacoustic approach. With its capability of completely monitoring a single, focused proton bunch with prompt readout and high repetition rate, I-BEAT is a promising approach to meet future requirements of experiments and applications in the field of laser-based ion acceleration. We demonstrate its functionality at two laser-driven ion sources for quantitative online determination of the kinetic energy distribution in the focus of single proton bunches.

Time-of-flight spectrometry of ultra-short, polyenergetic proton bunches.

A common approach for spectrum determination of polyenergetic proton bunches from laser-ion acceleration experiments is based on the time-of-flight (TOF) method. However, spectra obtained using this method are typically given in relative units or are estimated based on some prior assumptions on the energy distribution of the accelerated ions. In this work, we present a new approach using the TOF method that allows for an absolute energy spectrum reconstruction from a current signal acquired with a sub-nanosecond fast and 10 µm thin silicon detector. The reconstruction is based on solving a linear least-squares problem, taking into account the response function of the detection system. The general principle of signal generation and spectrum reconstruction by setting up an appropriate system response matrix is presented. Proof-of-principle experiments at a 12 MV Tandem accelerator using different nanosecond-short (quasi-)monoenergetic and polyenergetic proton bunches at energies up to 20 MeV were successfully performed. Within the experimental uncertainties of 2.4% and 12.1% for energy and particle number, respectively, reconstructed energy distributions were found in excellent agreement with the spectra calculated using Monte Carlo simulations and measured by a magnetic spectrometer. This TOF method can hence be used for absolute online spectrometry of laser-accelerated particle bunches.