Design study of a novel geometrical arrangement for an in-beam small animal positron emission tomography scanner.

Objective.We designed a geometrical solution for a small animal in-beam positron emission tomography (PET) scanner to be used in the project SIRMIO (Small animal proton irradiator for research in molecular image-guided radiation-oncology). The system is based on 56 scintillator blocks of pixelated LYSO crystals. The crystals are arranged providing a pyramidal-step shape to optimize the geometrical coverage in a spherical configuration.Approach.Different arrangements have been simulated and compared in terms of spatial resolution and sensitivity. The chosen setup enables us to reach a good trade-off between a solid angle coverage and sufficient available space for the integration of additional components of the first design prototype of the SIRMIO platform. The possibility of moving the mouse holder inside the PET scanner furthermore allows for achieving the optimum placement of the irradiation area for all the possible tumor positions in the body of the mouse. The work also includes a study of the scintillator material where LYSO and GAGG are compared with a focus on the random coincidence noise due to the natural radioactivity of Lutetium in LYSO, justifying the choice of LYSO for the development of the final system.Main results.The best imaging performance can be achieved with a sub-millimeter spatial resolution and sensitivity of 10% in the center of the scanner, as verified in thorough simulations of point sources. The simulation of realistic irradiation scenarios of proton beams in PMMA targets with/without air gaps indicates the ability of the proposed PET system to detect range shifts down to 0.2 mm.Significance.The presented results support the choice of the identified optimal design for a novel spherical in-beam PET scanner which is currently under commissioning for application to small animal proton and light ion irradiation, and which might find also application, e.g. for biological image-guidance in x-ray irradiation.

Radioactive Beams for Image-Guided Particle Therapy: The BARB Experiment at GSI.

Several techniques are under development for image-guidance in particle therapy. Positron (ß+) emission tomography (PET) is in use since many years, because accelerated ions generate positron-emitting isotopes by nuclear fragmentation in the human body. In heavy ion therapy, a major part of the PET signals is produced by ß+-emitters generated via projectile fragmentation. A much higher intensity for the PET signal can be obtained using ß+-radioactive beams directly for treatment. This idea has always been hampered by the low intensity of the secondary beams, produced by fragmentation of the primary, stable beams. With the intensity upgrade of the SIS-18 synchrotron and the isotopic separation with the fragment separator FRS in the FAIR-phase-0 in Darmstadt, it is now possible to reach radioactive ion beams with sufficient intensity to treat a tumor in small animals. This was the motivation of the BARB (Biomedical Applications of Radioactive ion Beams) experiment that is ongoing at GSI in Darmstadt. This paper will present the plans and instruments developed by the BARB collaboration for testing the use of radioactive beams in cancer therapy.

Toward a new treatment planning approach accounting for in vivo proton range verification.

Protons with modern pencil-beam scanning delivery are widely used in state-of-the-art radiotherapy. To reduce the unwanted effect of proton range uncertainties, prompt gamma (PG) monitoring is investigated and considered one of the most promising methods for real-time, in vivo range verification. Despite good correlation between the penetration depth of the PG signal and proton range in most cases, mismatch can occur especially because of tissue heterogeneities. Moreover, detectability and reproducibility of the prompt gamma signal critically depends on counting statistics. Nowadays, conventional treatment planning systems do not account for the degree of correlation between dose and PG signal nor the expected PG signal counting statistics, which considerably influences the possibility of a reliable verification of the intended beam range. Hence, in this project, we investigate a new treatment planning approach, in which the spot-by-spot conformities between PG and dose profiles (PG-dose correlation) as well as PG signal detectability and precision are taken into account based on a TPS optimizer. To investigate the feasibility of this idea, a research computational platform, combining Monte Carlo (MC, Geant4) pre-calculated pencil beams with the analytical Matlab-based TPS engine CERR, is used for treatment planning. Geant4 is employed for realistic simulation of the dose delivery and PG generation of all spots in the heterogeneous patient anatomy given by CT images. First of all, a treatment plan is created using a charged particle extension of CERR. Secondly, the PG fall-off positions of all individual pencil beams are evaluated and compared to the 80% distal dose fall-off positions. Thirdly, the PG-dose correlations of all spots are quantified. A new plan, in which a few spots with the best PG-dose correlation are boosted to ensure PG detectability with good precision, is then made. Finally, the optimized plan is fully recalculated on the same patient CT using Geant4, and the result is evaluated considering both plan quality and beam range monitorability. The evaluation shows that the re-optimized treatment plan is comparable to the initial plan in terms of dose distribution, dose averaged LET distribution and robustness, while fulfilling the set statistical conditions for reliable PG monitoring of the few automatically or manually selected spots. The method could thus complement, and for the selected pencil beams even overcome limitations of, alternative suggested approach such as pencil beam aggregation to provide sufficient counting statistics for precise PG range retrieval with good correlation to the treatment dose.

Full Monte Carlo-Based Biologic Treatment Plan Optimization System for Intensity Modulated Carbon Ion Therapy on Graphics Processing Unit.

PURPOSE: One of the major benefits of carbon ion therapy is enhanced biological effectiveness at the Bragg peak region. For intensity modulated carbon ion therapy (IMCT), it is desirable to use Monte Carlo (MC) methods to compute the properties of each pencil beam spot for treatment planning, because of their accuracy in modeling physics processes and estimating biological effects. We previously developed goCMC, a graphics processing unit (GPU)-oriented MC engine for carbon ion therapy. The purpose of the present study was to build a biological treatment plan optimization system using goCMC. METHODS AND MATERIALS: The repair-misrepair-fixation model was implemented to compute the spatial distribution of linear-quadratic model parameters for each spot. A treatment plan optimization module was developed to minimize the difference between the prescribed and actual biological effect. We used a gradient-based algorithm to solve the optimization problem. The system was embedded in the Varian Eclipse treatment planning system under a client-server architecture to achieve a user-friendly planning environment. We tested the system with a 1-dimensional homogeneous water case and 3 3-dimensional patient cases. RESULTS: Our system generated treatment plans with biological spread-out Bragg peaks covering the targeted regions and sparing critical structures. Using 4 NVidia GTX 1080 GPUs, the total computation time, including spot simulation, optimization, and final dose calculation, was 0.6 hour for the prostate case (8282 spots), 0.2 hour for the pancreas case (3795 spots), and 0.3 hour for the brain case (6724 spots). The computation time was dominated by MC spot simulation. CONCLUSIONS: We built a biological treatment plan optimization system for IMCT that performs simulations using a fast MC engine, goCMC. To the best of our knowledge, this is the first time that full MC-based IMCT inverse planning has been achieved in a clinically viable time frame.

Initial development of goCMC: a GPU-oriented fast cross-platform Monte Carlo engine for carbon ion therapy.

Monte Carlo (MC) simulation is considered as the most accurate method for calculation of absorbed dose and fundamental physics quantities related to biological effects in carbon ion therapy. To improve its computational efficiency, we have developed a GPU-oriented fast MC package named goCMC, for carbon therapy. goCMC simulates particle transport in voxelized geometry with kinetic energy up to 450 MeV u-1. Class II condensed history simulation scheme with a continuous slowing down approximation was employed. Energy straggling and multiple scattering were modeled. δ-electrons were terminated with their energy locally deposited. Four types of nuclear interactions were implemented in goCMC, i.e. carbon-hydrogen, carbon-carbon, carbon-oxygen and carbon-calcium inelastic collisions. Total cross section data from Geant4 were used. Secondary particles produced in these interactions were sampled according to particle yield with energy and directional distribution data derived from Geant4 simulation results. Secondary charged particles were transported following the condensed history scheme, whereas secondary neutral particles were ignored. goCMC was developed under OpenCL framework and is executable on different platforms, e.g. GPU and multi-core CPU. We have validated goCMC with Geant4 in cases with different beam energy and phantoms including four homogeneous phantoms, one heterogeneous half-slab phantom, and one patient case. For each case [Formula: see text] carbon ions were simulated, such that in the region with dose greater than 10% of maximum dose, the mean relative statistical uncertainty was less than 1%. Good agreements for dose distributions and range estimations between goCMC and Geant4 were observed. 3D gamma passing rates with 1%/1 mm criterion were over 90% within 10% isodose line except in two extreme cases, and those with 2%/1 mm criterion were all over 96%. Efficiency and code portability were tested with different GPUs and CPUs. Depending on the beam energy and voxel size, the computation time to simulate [Formula: see text] carbons was 9.9-125 s, 2.5-50 s and 60-612 s on an AMD Radeon GPU card, an NVidia GeForce GTX 1080 GPU card and an Intel Xeon E5-2640 CPU, respectively. The combined accuracy, efficiency and portability make goCMC attractive for research and clinical applications in carbon ion therapy.