Inelastic Nuclear Interactions Research Articles

A GPU-based fast Monte Carlo code that supports proton transport in magnetic field for radiation therapy.

This paper presents the effort to extend a previously reported code ARCHER, a GPU-based Monte Carlo (MC) code for coupled photon and electron transport, into protons including the consideration of magnetic fields. The proton transport is modeled using a Class-II condensed-history algorithm with continuous slowing-down approximation. The model includes ionization, multiple scattering, energy straggling, elastic and inelastic nuclear interactions, as well as deflection due to the Lorentz force in magnetic fields. An additional direction change is added for protons at the end of each step in the presence of the magnetic field. Secondary charge particles, except for protons, are terminated depositing kinetic energies locally, whereas secondary neutral particles are ignored. Each proton is transported step by step until its energy drops to below 0.5 MeV or when the proton leaves the phantom. The code is implemented using the compute unified device architecture (CUDA) platform for optimized GPU thread-level parallelism and efficiency. The code is validated by comparing it against TOPAS. Comparisons of dose distributions between our code and TOPAS for several exposure scenarios, ranging from single square beams in water to patient plan with magnetic fields, show good agreement. The 3D-gamma pass rate with a 2mm/2% criterion in the region with dose greater than 10% of the maximum dose is computed to be over 99% for all tested cases. Using a single NVIDIA TITAN V GPU card, the computational time of ARCHER is found to range from 0.82 to 4.54 seconds for 1×107 proton histories. Compared to a few hours running on TOPAS, this speed improvement is significant. This work presents, for the first time, the performance of a GPU-based MC code to simulate proton transportation magnetic fields, demonstrating the feasibility of accurate and efficient dose calculations in potential magnetic resonance imaging (MRI)-guided proton therapy.

Nuclear reaction model for deuterons below 150 MeV/n in FLUKA

The lack of a dedicated model for nuclear reactions of deuterons at energies below 150 MeV/n has been a long-standing shortcoming of the general-purpose radiation-transport code FLUKA, preventing its use in applications either driven by deuteron beams or where reinteracting secondary deuterons play a significant role. A dedicated modelling effort has been undertaken to overcome this limitation. In this contribution, the essential ingredients of FLUKA’s novel model for low-energy deuteron nuclear inelastic interactions are briefly described, including nucleon stripping to the continuum, nucleon transfer to discrete states, and elastic break-up. The resulting model has been extensively benchmarked against experimental data prior to its inclusion in FLUKA v4-2.0.

Development and benchmarking of the first fast Monte Carlo engine for helium ion beam dose calculation: MonteRay.

Monte Carlo (MC) simulations are considered the gold-standard for accuracy in radiotherapy dose calculation; however, general purpose MC engines are computationally demanding and require long runtimes. For this reason, several groups have recently developed fast MC systems dedicated mainly to photon and proton external beam therapy, affording both speed and accuracy. To support research and clinical activities at the Heidelberg Ion-beam Therapy Center (HIT) with actively scanned helium ion beams, this work presents MonteRay, the first fast MC dose calculation engine for helium ion therapy. MonteRay is a CPU MC dose calculation engine written in C++, capable of simulating therapeutic proton and helium ion beams. In this work, development steps taken to include helium ion beams in MonteRay are presented. A detailed description of the newly implemented physics models for helium ions, for example, for multiple coulomb scattering and inelastic nuclear interactions, is provided. MonteRay dose computations of helium ion beams are evaluated using a comprehensive validation dataset, including measurements of spread-out Bragg peaks (SOBPs) with varying penetration depths/field sizes, measurements with an anthropomorphic phantom and FLUKA simulations of a patient plan. Improvement in computational speed is demonstrated in comparison against reference FLUKA simulations. Dosimetric comparisons between MonteRay and measurements demonstrated good agreement. Comparing SOBPs at 5, 12.5, and 20cm depth, mean absolute percent dose differences were 0.7%, 0.7%, and 1.4%, respectively. Comparison against measurements behind an anthropomorphic head phantom revealed mean absolute dose differences of about 1.2% (FLUKA: 1.5%) with per voxel errors ranging from -4.5% to 4.1% (FLUKA: -6% to 3%). Computed global 3%/3mm 3D-gamma passing rates of ∼99% were achieved, exceeding those previously reported for an analytical dose engine. Comparisons against FLUKA simulations for a patient plan revealed local 2%/2mm 3D-gamma passing rates of 98%. Compared to FLUKA in voxelized geometries, MonteRay saw run-time reductions ranging from 20× to 60×, depending on the beam's energy. MonteRay, the first fast MC engine dedicated to helium ion therapy, has been successfully developed with a focus on both speed and accuracy. Validations against dosimetric measurements in homogeneous and heterogeneous scenarios and FLUKA MC calculations have proven the validity of the physical models implemented. Timing comparisons have shown significant speedups between 20 and 60 when compared to FLUKA, making MonteRay viable for clinical routine. MonteRay will support research and clinical practice at HIT, for example, TPS development, validation and treatment design for upcoming clinical trials for raster-scanned helium ion therapy.

Feasibility of crystal-assisted collimation in the CERN accelerator complex

Bent silicon crystals mounted on high-accuracy angular actuators were installed in the CERN Super Proton Synchrotron (SPS) and extensively tested to assess the feasibility of crystal-assisted collimation in circular hadron colliders. The adopted layout was exploited and regularly upgraded for about a decade by the UA9 Collaboration. The investigations provided the compelling evidence of a strong reduction of beam losses induced by nuclear inelastic interactions in the aligned crystals in comparison with amorphous orientation. A conceptually similar device, installed in the betatron cleaning insertion of CERN Large Hadron Collider (LHC), was operated through the complete acceleration and storage cycle and demonstrated a large reduction of the background leaking from the collimation region and radiated into the cold sections of the accelerator and the experimental detectors. The implemented layout and the relevant results of the beam tests performed in the SPS and in the LHC with stored proton and ion beams are extensively discussed.

An Analysis Scheme for 3D Visualization of Positron Emitting Radioisotopes Using Positron Emission Mammography System

Proton range monitoring and verification is important to enhance the effectiveness of treatment by ensuring that the correct dose is delivered to the correct location. Upon proton irradiation, different positron emitting radioisotopes are produced by the inelastic nuclear interactions of protons with the target elements. Recently, it was reported that the 16O(p,2p2n)13N reaction has a relatively low threshold energy, and it could be potentially used for proton range verification. In the present work, we have proposed an analysis scheme (i.e., algorithm) for the extraction and three-dimensional visualization of positron emitting radioisotopes. The proposed step-by-step analysis scheme was tested using our own experimentally obtained dynamic data from a positron emission mammography (PEM) system (our developed PEMGRAPH system). The experimental irradiation was performed using an azimuthally varying field (AVF) cyclotron with a 80 MeV monoenergetic pencil-like beam. The 3D visualization showed promising results for proton-induced radioisotope distribution. The proposed scheme and developed tools would be useful for the extraction and 3D visualization of positron emitting radioisotopes and in turn for proton range monitoring and verification.

Proton path reconstruction for proton computed tomography using neural networks

The most likely path formalism (MLP) is widely established as the most statistically precise method for proton path reconstruction in proton computed tomography. However, while this method accounts for small-angle multiple coulomb scattering (MCS) and energy loss, inelastic nuclear interactions play an influential role in a significant number of proton paths. By applying cuts based on energy and direction, tracks influenced by nuclear interactions are largely discarded from the MLP analysis. In this work we propose a new method to estimate the proton paths based on a deep neural network (DNN). Through this approach, estimates of proton paths equivalent to MLP predictions have been achieved in the case where only MCS occurs, together with an increased accuracy when nuclear interactions are present. Moreover, our tests indicate that the DNN algorithm can be considerably faster than the MLP algorithm.

Helium radiography with a digital tracking calorimeter—a Monte Carlo study for secondary track rejection

Radiation therapy using protons and heavier ions is a fast-growing therapeutic option for cancer patients. A clinical system for particle imaging in particle therapy would enable online patient position verification, estimation of the dose deposition through range monitoring and a reduction of uncertainties in the calculation of the relative stopping power of the patient.Several prototype imaging modalities offer radiography and computed tomography using protons and heavy ions. A Digital Tracking Calorimeter (DTC), currently under development, has been proposed as one such detector. In the DTC 43 longitudinal layers of laterally stacked ALPIDE CMOS monolithic active pixel sensor chips are able to reconstruct a large number of simultaneously recorded proton tracks.In this study, we explored the capability of the DTC for helium imaging which offers favorable spatial resolution over proton imaging. Helium ions exhibit a larger cross section for inelastic nuclear interactions, increasing the number of produced secondaries in the imaged object and in the detector itself. To that end, a filtering process able to remove a large fraction of the secondaries was identified, and the track reconstruction process was adapted for helium ions.By filtering on the energy loss along the tracks, on the incoming angle and on the particle ranges, 97.5% of the secondaries were removed. After passing through 16 cm water, 50.0% of the primary helium ions survived; after the proposed filtering 42.4% of the primaries remained; finally after subsequent image reconstruction 31% of the primaries remained. Helium track reconstruction leads to more track matching errors compared to protons due to the increased available focus strength of the helium beam. In a head phantom radiograph, the Water Equivalent Path Length error envelope was 1.0 mm for helium and 1.1 mm for protons. This accuracy is expected to be sufficient for helium imaging for pre-treatment verification purposes.

Angular asymmetry of the nuclear interaction probability of high energy particles in short bent crystals

The rate of inelastic nuclear interactions in a short bent silicon crystal was precisely measured for the first time using a 180 GeV/c positive hadron beam produced in the North Experimental Area of the CERN SPS. An angular asymmetry dependence on the crystal orientation in the vicinity of the planar channeling minimum has been observed. For the inspected crystal, this probability is about sim 20% larger than in the amorphous case because of the atomic density increase along the particle trajectories in the angular range of volume reflection, whose dimension is determined by the crystal bending angle. Instead, for the opposite angular orientation with respect to the planar channeling, there is a smaller probability excess of sim 4%.

Enhancement of the Inelastic Nuclear Interaction Rate in Crystals via Antichanneling.

The interaction rate of a charged particle beam with the atomic nuclei of a target varies significantly if the target has a crystalline structure. In particular, under specific orientations of the target with respect to the incident beam, the probability of inelastic interaction with nuclei can be enhanced with respect to the unaligned case. This effect, which can be named antichanneling, can be advantageously used in the cases where the interaction between beam and target has to be maximized. Here we propose to use antichanneling to increase the radioisotope production yield via cyclotron. A dedicated set of experimental measurements was carried out at the INFN Legnaro Laboratories with the AN2000 and CN accelerators to prove the existence of the antichanneling effect. The variation of the interaction yield at hundreds of keV to MeV energies was observed by means of sapphire and indium phosphide crystals, achieving an enhancement of the interaction rate up to 73% and 25%, respectively. Such a result may pave the way to the development of a novel type of nozzle for the existing cyclotrons, which can exploit crystalline materials as targets for radioisotope production, especially to enhance the production rate for expensive prime materials with minor upgrades of the current instrumentation.

Recent improvements for the lepton propagator PROPOSAL

The lepton propagator PROPOSAL is a Monte-Carlo Simulation library written in C++, propagating high energy muons and other charged particles over large distances in media. In this article, a restructuring of the code is described, which yields a performance improvement of up to 30%. For an improved accuracy of the propagation processes, more exact calculations of the leptonic and hadronic decay process and more precise parametrizations for the interaction cross sections are now available. The new modular structure allows a more flexible and customized usage, which is further facilitated with a Python interface. Program summaryProgram Title: PROPOSALProgram Files doi:http://dx.doi.org/10.17632/g478pjdcxy.1Licensing provisions: LGPLProgramming language: C++Nature of problem: Propagation of charged particles over large distances in three dimensions through different kinds of media. These particles lose their energy stochastically via the processes of ionization, pair production, bremsstrahlung and inelastic nuclear interaction and eventually decay, producing secondary particles along the trajectory.Solution method: Monte-Carlo simulation. The program samples the next stochastic interaction point, the type of interaction and the amount of energy lost in the interaction until either the particle decays, its energy is below a certain threshold or it reaches a given distance. To improve the performance and to deal with the bremsstrahlung divergence at small energy losses, an adaptable energy cut is used below which all losses are treated continuously. The use of interpolation tables further reduces computation time.The sampled energy till the next stochastic loss is smeared out with a Gaussian randomization inside the physically allowed limits of the continuous losses using the second moment of the summed processes to avoid artifacts introduced by the energy cut. The deviation from a straight trajectory is evaluated using the multiple scattering calculation by Molière or the Highland parametrization, a Gaussian approximation to Molière’s theory. Multiple kinds of parametrizations are also available for bremsstrahlung, pair production and inelastic nuclear interaction to study the effects of the uncertainty of the cross sections on the propagation and further simulation steps.

Theoretical Research and Development of a Clinical Setup Prototype for Online Monitoring of the Bragg Peak Position for the Prometheus Proton Therapy Complex

In this paper we report a result of theoretical studies of the method for real-time monitoring of the Bragg peak position in a water phantom during scanned proton pencil-beam irradiation. This method is based on the detection of the prompt gamma rays emitted orthogonally to the beam direction produced as a result of inelastic nuclear interactions of primary particles. The principal parameters of the clinical setup prototype and the accuracy of determining the longitudinal coordinates of the Bragg peak position were found on the basis of statistical simulation using the RTS&T Monte Carlo multiparticle transport code.

Study of inelastic nuclear interactions of 400 GeV/c protons in bent silicon crystals for beam steering purposes

Inelastic nuclear interaction probability of 400 GeV/c protons interacting with bent silicon crystals was investigated, in particular for both types of crystals installed at the CERN Large Hadron Collider for beam collimation purposes. In comparison to amorphous scattering interaction, in planar channeling this probability is sim 36% for the quasi-mosaic type (planes (111)), and sim 27% for the strip type (planes (110)). Moreover, the absolute inelastic nuclear interaction probability in the axial channeling orientation, along the langle 110rangle axis, was estimated for the first time, finding a value of 0.6% for a crystal 2 mm long along the beam direction, with a bending angle of 55 upmu rad. This value is more than two times lower with respect to the planar channeling orientation of the same crystal, and increases with the vertical angular misalignment. Finally, the correlation between the inelastic nuclear interaction probability in the planar channeling and the silicon crystal curvature is reported.

Simulation of orientational coherent effects via Geant4

Simulation of orientational coherent effects via Geant4 beam manipulation of high-and very-high-energy particle beams is a hot topic in accelerator physics. Coherent effects of ultra-relativistic particles in bent crystals allow the steering of particle trajectories thanks to the strong electrical field generated between atomic planes. Recently, a collimation experiment with bent crystals was carried out at the CERN-LHC, paving the way to the usage of such technology in current and future accelerators. Geant4 is a widely used object-oriented tool-kit for the Monte Carlo simulation of the interaction of particles with matter in high-energy physics. Moreover, its areas of application include also nuclear and accelerator physics, as well as studies in medical and space science. We present the first Geant4 extension for the simulation of orientational effects in straight and bent crystals for high energy charged particles. The model allows the manipulation of particle trajectories by means of straight and bent crystals and the scaling of the cross sections of hadronic and electromagnetic processes for channeled particles. Based on such a model, an extension of the Geant4 toolkit has been developed. The code and the model have been validated by comparison with published experimental data regarding the deflection efficiency via channeling and the variation of the rate of inelastic nuclear interactions.

Impact of the Border Crossing Effects on the DCNU for Pixel Arrays Irradiated With High Energy Protons

A 3D Monte Carlo model based on GEANT4, dedicated to the simulation of the Dark Current Non Uniformity (DCNU) in image sensors, is used to study the impact of the geometrical dimensions of pixels, on the shape of the DCNU. The study focuses on the inelastic nuclear interaction which is the dominant process for high degradation levels (tail of the DCNU). Simulations made on different pixel geometries demonstrate that border crossing effect impacts the shape of the DCNU only for very small pixel dimensions (depleted volume <; 8 μm <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> ).

High-efficiency deflection of high energy protons due to channeling along the 〈110〉 axis of a bent silicon crystal

Abstract A deflection efficiency of about 61 % was observed for 400 GeV/c protons due to channeling, most strongly along the 〈110〉 axis of a bent silicon crystal. It is comparable with the deflection efficiency in planar channeling and considerably larger than in the case of the 〈111〉 axis. The measured probability of inelastic nuclear interactions of protons in channeling along the 〈110〉 axis is only about 10 % of its amorphous level whereas in channeling along the (110) planes it is about 25 % . High efficiency deflection and small beam losses make this axial orientation of a silicon crystal a useful tool for the beam steering of high energy charged particles.

Observation of channeling for 6500 GeV/ c protons in the crystal assisted collimation setup for LHC

Abstract Two high-accuracy goniometers equipped with two bent silicon crystals were installed in the betatron cleaning insertion of the CERN Large Hadron Collider (LHC) during its long shutdown. First beam tests were recently performed at the LHC with 450 GeV/ c and 6500 GeV/ c stored proton beams to investigate the feasibility of beam halo collimation assisted by bent crystals. For the first time channeling of 6500 GeV/ c protons was observed in a particle accelerator. A strong reduction of beam losses due to nuclear inelastic interactions in the aligned crystal in comparison with its amorphous orientation was detected. The loss reduction value was about 24. Thus, the results show that deflection of particles by a bent crystal due to channeling is effective for this record particle energy.

Fast multipurpose Monte Carlo simulation for proton therapy using multi- and many-core CPU architectures.

Accuracy in proton therapy treatment planning can be improved using Monte Carlo (MC) simulations. However the long computation time of such methods hinders their use in clinical routine. This work aims to develop a fast multipurpose Monte Carlo simulation tool for proton therapy using massively parallel central processing unit (CPU) architectures. A new Monte Carlo, called MCsquare (many-core Monte Carlo), has been designed and optimized for the last generation of Intel Xeon processors and Intel Xeon Phi coprocessors. These massively parallel architectures offer the flexibility and the computational power suitable to MC methods. The class-II condensed history algorithm of MCsquare provides a fast and yet accurate method of simulating heavy charged particles such as protons, deuterons, and alphas inside voxelized geometries. Hard ionizations, with energy losses above a user-specified threshold, are simulated individually while soft events are regrouped in a multiple scattering theory. Elastic and inelastic nuclear interactions are sampled from ICRU 63 differential cross sections, thereby allowing for the computation of prompt gamma emission profiles. MCsquare has been benchmarked with the gate/geant4 Monte Carlo application for homogeneous and heterogeneous geometries. Comparisons with gate/geant4 for various geometries show deviations within 2%-1 mm. In spite of the limited memory bandwidth of the coprocessor simulation time is below 25 s for 10(7) primary 200 MeV protons in average soft tissues using all Xeon Phi and CPU resources embedded in a single desktop unit. MCsquare exploits the flexibility of CPU architectures to provide a multipurpose MC simulation tool. Optimized code enables the use of accurate MC calculation within a reasonable computation time, adequate for clinical practice. MCsquare also simulates prompt gamma emission and can thus be used also for in vivo range verification.

Influence of nuclear interactions in body tissues on tumor dose in carbon-ion radiotherapy.

In carbon-ion radiotherapy treatment planning, the planar integrated dose (PID) measured in water is applied to the patient dose calculation with density scaling using the stopping power ratio. Since body tissues are chemically different from water, this dose calculation can be subject to errors, particularly due to differences in inelastic nuclear interactions. In recent studies, the authors proposed and validated a PID correction method for these errors. In the present study, the authors used this correction method to assess the influence of these nuclear interactions in body tissues on tumor dose in various clinical cases. Using 10-20 cases each of prostate, head and neck (HN), bone and soft tissue (BS), lung, liver, pancreas, and uterine neoplasms, the authors first used treatment plans for carbon-ion radiotherapy without nuclear interaction correction to derive uncorrected dose distributions. The authors then compared these distributions with recalculated distributions using the nuclear interaction correction (corrected dose distributions). Median (25%/75% quartiles) differences between the target mean uncorrected doses and corrected doses were 0.2% (0.1%/0.2%), 0.0% (0.0%/0.0%), -0.3% (-0.4%/-0.2%), -0.1% (-0.2%/-0.1%), -0.1% (-0.2%/0.0%), -0.4% (-0.5%/-0.1%), and -0.3% (-0.4%/0.0%) for the prostate, HN, BS, lung, liver, pancreas, and uterine cases, respectively. The largest difference of -1.6% in target mean and -2.5% at maximum were observed in a uterine case. For most clinical cases, dose calculation errors due to the water nonequivalence of the tissues in nuclear interactions would be marginal compared to intrinsic uncertainties in treatment planning, patient setup, beam delivery, and clinical response. In some extreme cases, however, these errors can be substantial. Accordingly, this correction method should be routinely applied to treatment planning in clinical practice.

SU‐E‐T‐663: Radiation Properties of a Water‐Equivalent Material Formulated Using the Stoichiometric Analysis Method in Heavy Charged Particle Therapy

Purpose:A material has been designed to be employed as water‐equivalent in particle therapy using a previously established stoichiometric analysis method (SAM). After manufacturing, experimental verification of the material's water‐equivalent path length (WEPL) and analysis of its total inelastic nuclear interaction cross sections for proton beams were performed.Methods:Using the SAM, we optimized the material composed of three base materials, i.e., polyurethane, calcium carbonate and microspheres. From the elemental composition of the compound, electron density, linear attenuation coefficients, particle stopping powers and inelastic nuclear cross sections for protons using data from ICRU 63 were calculated. The calculations were then compared to Hounsfield units (HUs) measured with 350 mAs at 80, 100, 120 and 140 kV and the WEPLs measured with three different ions: proton (106.8 MeV/u), helium (107.93 MeV/u) and carbon (200.3 MeV/u).Results:The material's measured HUs (0.7±3.0 to 2.6±6.2 HU) as well as its calculated relative electron density (1.0001) are in close agreement with water as reference. The WEPLs measured on a 20.00 mm thick target were 20.16±0.12, 20.29±0.12 and 20.38±0.12 mmH2O for proton, helium and carbon ions, respectively. Within measurement uncertainties, these values verified the calculated WEPLs of 20.28 mmH2O (proton), 20.28 mmH2O (helium) and 20.26 mmH2O (carbon). Moreover, the calculated proton inelastic cross sections of the material differed only by 0.89% (100 MeV/u) and 0.01% (200 MeV/u) when compared to water.Conclusion:The SAM is capable of optimizing material with defined properties, e.g., HU, electron density, WEPL and inelastic nuclear interaction cross section for particle therapy. Such material will have a wide range of applications amongst others absolute dosimetry.This work was supported by grant ZIM KF2137107AK4 from the German Federal Ministry for Economic Affairs and Energy.

Influence of the patient size on secondary doses delivered in light ion beams

*Secondary radiation exposure of patients in light ion therapy is of concern due to possible normal tissue damage and risk of induction of secondary cancers. Neutrons, protons and heavier ions are generated by nuclear inelastic interactions of primary ions both in the beam line and in the patient. The patient is exposed to a complex radiation field and secondary doses can be deposited in normal tissues both close to and relatively far from the treated volume. The energy distribution of secondary particles and secondary doses delivered to different organs were studied by the MC code SHIELD-HIT10 using the anthropomorphic phantoms representing a 10-year-old child (CHILD-HIT) and an adult male (ADAM-HIT). Brain tumor irradiations were simulated with approximated scanned beams of 1 H, 7 Li and 12 C ions. The influence of patient size on the secondary dose distributions were studied using the same target definition and irradiation geometry with a lateral beam. For the scanned beams, the secondary organ absorbed doses normalized per absorbed dose to the treated volume (brain tumor) were in the range 1 nGy/Gy - 0.1 mGy/Gy and the absorbed doses in the CHILD-HIT phantom were higher than in ADAM-HIT by up to a factor of 5, depending on organ.