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 Monte Carlo feasibility study on quantitative laser-driven proton radiography.

Laser-accelerated proton bunches with kinetic energies up to several tens of MeV and at repetition rates in the order of Hz are nowadays achievable at several research centres housing high-power laser system. The unique features of such ultra-short bunches are also arousing interest in the field of radiological and biomedical applications. For many of these applications, accurate positioning of the biological target is crucial, raising the need for on-site imaging. One convenient option is proton radiography, which can exploit the polyenergetic spectrum of laser-accelerated proton bunches. We present a Monte Carlo (MC) feasibility study to assess the applicability and potential of laser-driven proton radiography of millimetre to centimetre sized objects. Our radiography setup consists of a thin time-of-flight spectrometer operated in transmission prior to the object and a pixelated silicon detector for imaging. Proton bunches with kinetic energies up to 20MeV and up to 100MeV were investigated. The water equivalent thickness (WET) of the traversed material is calculated from the energy deposition inside an imaging detector, using an online generated calibration curve that is based on a MC generated look-up table and the reconstructed proton energy distribution. With a dose of 43mGy for a 1mm thin object imaged with protons up to 20MeV, the reconstructed WET of defined regions-of-interest was within 1.5% of the ground truth values. The spatial resolution, which strongly depends on the gap between object and imaging detector, was 2.5lpmm-1 for a realistic distance of 5mm. Due to this relatively high imaging dose, our proposed setup for laser-driven proton radiography is currently limited to objects with low radio-sensitivity, but possibilities for further dose reduction are presented and discussed.

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.

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.

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.

Gel dosimetry for three dimensional proton range measurements in anthropomorphic geometries.

Proton beams used for radiotherapy have potential for superior sparing of normal tissue, although range uncertainties are among the main limiting factors in the accuracy of dose delivery. The aim of this study was to benchmark an N-vinylpyrrolidone based polymer gel to perform three-dimensional measurement of geometric proton beam characteristics and especially to test its suitability as a range probe in combination with an anthropomorphic phantom. For single proton pencil beams as well as for 3×3cm2 mono-energy layers depth dose profiles, lateral dose distribution at different depths and proton range were evaluated in simple cubic gel phantoms at different energies from 75 to 115MeV and different dose levels. In addition, a 90MeV mono-energetic beam was delivered to an anthropomorphic 3D printed head phantom, which was filled with gel. Subsequently, all phantoms underwent magnetic resonance imaging using an axial pixel size of 0.68-0.98mm and with slice thicknesses of 2 or 3mm to derive a 3-dimensional distribution of the T2 relaxation time, which correlates with radiation dose. Indices describing lateral dose distribution and proton range were compared against predictions from a treatment planning system (TPS, for cubic and head phantoms) and Monte Carlo simulations (MC, for the head phantom) after manual rigid co-registration with the T2 relaxation time datasets. For all pencil beams, the FWHM agreement with TPS was better than 1mm or 7%. For the mono-energetic layer, the agreement with TPS in this respect was even better than 0.3mm in each case. With respect to range, results from gel measurements differed no more than 0.9mm (1.6%) from values predicted by TPS. In case of the anthropomorphic phantom, deviations with respect to a nominal range of about 61mm as well as in FWHM were slightly higher, namely within 1.0mm and 1.1mm respectively. Average deviations between gel and TPS/MC were similar (-0.3mm±0.4mm/-0.2±0.5mm). In conclusion, polymer gel dosimetry was found to be a valuable tool to determine geometric proton beam properties three-dimensionally and with high spatial resolution in simple cubic as well as in a more complex anthropomorphic phantom. Post registration range errors of the order of 1mm could be achieved. The additional registration uncertainty (95%) was 1mm.

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.

AN ONLINE, RADIATION HARD PROTON ENERGY-RESOLVING SCINTILLATOR STACK FOR LASER-DRIVEN PROTON BUNCHES.

We report on a scintillator-based online detection system for the spectral characterization of polychromatic proton bunches. Using up to nine stacked layers of radiation hard polysiloxane scintillators, coupled to and readout edge-on by a large area pixelated CMOS detector, impinging polychromatic proton bunches were characterized. The energy spectra were reconstructed using calibration data and simulated using Monte-Carlo simulations. Despite the scintillator stack showed some problems like thickness inhomogeneities and unequal layer coupling, the prototype allows to obtain a first estimate of the energy spectrum of proton beams.

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.

Systematic out-of-field secondary neutron spectrometry and dosimetry in pencil beam scanning proton therapy.

BACKGROUND AND PURPOSE: Systematic investigation of the energy and angular dependence of secondary neutron fluence energy distributions and ambient dose equivalents values (H*(10)) inside a pencil beam scanning proton therapy treatment room using a gantry. MATERIALS AND METHODS: Neutron fluence energy distributions were measured with an extended-range Bonner sphere spectrometer featuring ³He proportional counters, at four positions at 0°, 45°, 90°, and 135° with respect to beam direction and at a distance of 2 m from the isocenter. The energy distribution of secondary neutrons was investigated for initial proton beam energies of 75 MeV, 140 MeV, and 200 MeV, respectively, using a 2D scanned irradiation field of 11 × 11 cm² delivered to a 30 × 30 × 30 cm³ PMMA phantom. Additional measurements were performed at a proton energy of 118 MeV including a 5 cm range-shifter (PMMA), yielding a Bragg peak position similar to that of 75 MeV protons. RESULTS: Ambient dose equivalent values from 0.3 µSv/Gy (75 MeV; 90°) to 24 µSv/Gy (200 MeV; 0°) were measured inside the treatment room at a distance of 2 m from the isocenter. H*(10) values were lower (by factors of up to 7.2 (at 45°)) at 75 MeV compared to those at 118 MeV with the 5 cm range-shifter. At 0° and 45°, an evaporation peak was found in the measured neutron fluence energy distributions, at neutron energies around MeV, which contributes about 50% to total H*(10) values, for all investigated proton beam energies. CONCLUSIONS: This study showed a pronounced increase of secondary neutron H*(10) values inside the proton treatment room with increasing proton energy without beam modifiers. For example, in beam direction this increase was about a factor of 50 when protons of 75 MeV and 200 MeV were compared. The existence of a peak of secondary neutrons in the MeV region was demonstrated in beam direction (0°). This peak is due to evaporation neutrons produced in the existing surrounding materials such as those used for the gantry. Therefore, any simulation of the secondary neutrons within a proton treatment room must take these materials into account. In addition, the results obtained here show that the use of a range-shifter increases the production of secondary neutrons inside the treatment room. Using a range-shifter, the higher neutron doses observed mainly result from the higher incident proton energy (118 MeV instead of 75 MeV when no range-shifter was used), due to higher neutron production cross-sections.