PET monitoring of cancer therapy with3He and12C beams: a study with the GEANT4 toolkit

The microdosimetric extension in TOPAS: development and comparison with published data

Microdosimetric energy depositions have been suggested as a key variable for the modeling of the relative biological effectiveness (RBE) in proton and ion radiation therapy. However, microdosimetry has been underutilized in radiation therapy. Recent advances in detector technology allow the design of new mico- and nano-dosimeters. At the same time Monte Carlo (MC) simulations have become more widely used in radiation therapy. In order to address the growing interest in the field, a microdosimetric extension was developed in TOPAS. The extension provides users with the functionality to simulate microdosimetric spectra as well as the contribution of secondary particles to the spectra, calculate microdosimetric parameters, and determine RBE with a biological weighting function approach or with the microdosimetric kinetic (MK) model. Simulations were conducted with the extension and the results were compared with published experimental data and other simulation results for three types of microdosimeters, a spherical tissue equivalent proportional counter (TEPC), a cylindrical TEPC and a solid state microdosimeter. The corresponding microdosimetric spectra obtained with TOPAS from the plateau region to the distal tail of the Bragg curve generally show good agreement with the published data.

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.

Experimental Approach to Evaluate the 11C Perfusion and Diffusion in Small Animal Tissues for HadronPET Applications

The development of a reliable dose monitoring system in hadron therapy is essential in order to control the treatment plan delivery. Positron Emission Tomography (PET) is the only method used in clinics nowadays for quality assurance. However, the accuracy of this method is limited by the loss of signal due to the biological washout processes. Up to the moment, very few studies measured the washout processes and there is no database of washout data as a function of the tissue and radioisotope. One of the main difficulties is related to the complexity of such measurements, along with the limited time slots available in hadron therapy facilities. Thus, in this work, we proposed an alternative in vivo methodology for the measurement and modeling of the biological washout parameters without any radiative devices. It consists in the implementation of a point-like radioisotope source by direct injection on the tissues of interest and its measurement by means of high-resolution preclinical PET systems. In particular, the washout of 11C carbonate radioisotopes was assessed, considering that 11C is is the most abundant β+ emitter produced by carbon beams. 11C washout measurements were performed in several tissues of interest (brain, muscle and 9L tumor xenograf) in rodents (Wistar rat). Results show that the methodology presented is sensitive to the washout variations depending on the selected tissue. Finally, a first qualitative correlation between 11C tumor washout properties and tumor metabolism (via 18F-FDG tracer uptake) was found.

PETSTEP: Generation of synthetic PET lesions for fast evaluation of segmentation methods

Purpose

This work describes PETSTEP (PET Simulator of Tracers via Emission Projection): a faster and more accessible alternative to Monte Carlo (MC) simulation generating realistic PET images, for studies assessing image features and segmentation techniques.

Methods

PETSTEP was implemented within Matlab as open source software. It allows generating three-dimensional PET images from PET/CT data or synthetic CT and PET maps, with user-drawn lesions and user-set acquisition and reconstruction parameters. PETSTEP was used to reproduce images of the NEMA body phantom acquired on a GE Discovery 690 PET/CT scanner, and simulated with MC for the GE Discovery LS scanner, and to generate realistic Head and Neck scans. Finally the sensitivity (S) and Positive Predictive Value (PPV) of three automatic segmentation methods were compared when applied to the scanner-acquired and PETSTEP-simulated NEMA images.

Results

PETSTEP produced 3D phantom and clinical images within 4 and 6 min respectively on a single core 2.7 GHz computer. PETSTEP images of the NEMA phantom had mean intensities within 2% of the scanner-acquired image for both background and largest insert, and 16% larger background Full Width at Half Maximum. Similar results were obtained when comparing PETSTEP images to MC simulated data. The S and PPV obtained with simulated phantom images were statistically significantly lower than for the original images, but led to the same conclusions with respect to the evaluated segmentation methods.

Conclusions

PETSTEP allows fast simulation of synthetic images reproducing scanner-acquired PET data and shows great promise for the evaluation of PET segmentation methods.

Range Verification Methods in Particle Therapy: Underlying Physics and Monte Carlo Modeling

Hadron therapy allows for highly conformal dose distributions and better sparing of organs-at-risk, thanks to the characteristic dose deposition as function of depth. However, the quality of hadron therapy treatments is closely connected with the ability to predict and achieve a given beam range in the patient. Currently, uncertainties in particle range lead to the employment of safety margins, at the expense of treatment quality. Much research in particle therapy is therefore aimed at developing methods to verify the particle range in patients. Non-invasive in vivo monitoring of the particle range can be performed by detecting secondary radiation, emitted from the patient as a result of nuclear interactions of charged hadrons with tissue, including β (+) emitters, prompt photons, and charged fragments. The correctness of the dose delivery can be verified by comparing measured and pre-calculated distributions of the secondary particles. The reliability of Monte Carlo (MC) predictions is a key issue. Correctly modeling the production of secondaries is a non-trivial task, because it involves nuclear physics interactions at energies, where no rigorous theories exist to describe them. The goal of this review is to provide a comprehensive overview of various aspects in modeling the physics processes for range verification with secondary particles produced in proton, carbon, and heavier ion irradiation. We discuss electromagnetic and nuclear interactions of charged hadrons in matter, which is followed by a summary of some widely used MC codes in hadron therapy. Then, we describe selected examples of how these codes have been validated and used in three range verification techniques: PET, prompt gamma, and charged particle detection. We include research studies and clinically applied methods. For each of the techniques, we point out advantages and disadvantages, as well as clinical challenges still to be addressed, focusing on MC simulation aspects.

An introduction to molecular imaging in radiation oncology: a report by the AAPM Working Group on Molecular Imaging in Radiation Oncology (WGMIR)

Molecular imaging is the direct or indirect noninvasive monitoring and recording of the spatial and temporal distribution of in vivo molecular, genetic, and/or cellular processes for biochemical, biological, diagnostic, or therapeutic applications. Molecular images that indicate the presence of malignancy can be acquired using optical, ultrasonic, radiologic, radionuclide, and magnetic resonance techniques. For the radiation oncology physicist in particular, these methods and their roles in molecular imaging of oncologic processes are reviewed with respect to their physical bases and imaging characteristics, including signal intensity, spatial scale, and spatial resolution. Relevant molecular terminology is defined as an educational assist. Current and future clinical applications in oncologic diagnosis and treatment are discussed. National initiatives for the development of basic science and clinical molecular imaging techniques and expertise are reviewed, illustrating research opportunities in as well as the importance of this growing field.

In vivo proton range verification: a review

Protons are an interesting modality for radiotherapy because of their well defined range and favourable depth dose characteristics. On the other hand, these same characteristics lead to added uncertainties in their delivery. This is particularly the case at the distal end of proton dose distributions, where the dose gradient can be extremely steep. In practice however, this gradient is rarely used to spare critical normal tissues due to such worries about its exact position in the patient. Reasons for this uncertainty are inaccuracies and non-uniqueness of the calibration from CT Hounsfield units to proton stopping powers, imaging artefacts (e.g. due to metal implants) and anatomical changes of the patient during treatment. In order to improve the precision of proton therapy therefore, it would be extremely desirable to verify proton range in vivo, either prior to, during, or after therapy. In this review, we describe and compare state-of-the art in vivo proton range verification methods currently being proposed, developed or clinically implemented.

In vivo photoactivation without "light": use of Cherenkov radiation to overcome the penetration limit of light

PURPOSE

The poor tissue penetration of visible light has been a major barrier for optical imaging, photoactivatable conversions, and photodynamic therapy for in vivo targets with depths beyond 10 mm. In this report, as a proof-of-concept, we demonstrated that a positron emission tomography (PET) radiotracer, 2-deoxy-2-[(18)F]fluoro-D-glucose ((18)FDG), could be used as an alternative light source for photoactivation.

PROCEDURES

We utilized (18)FDG, which is a metabolic activity-based PET probe, as a source of light to photoactivate caged luciferin in a breast cancer animal model expressing luciferase.

RESULTS

Bioluminescence produced from luciferin allowed for the real-time monitoring of Cherenkov radiation-promoted uncaging of the substrate.

CONCLUSION

The proposed method may provide a very important option for in vivo photoactivation, in particular for activation of photosensitizers for photodynamic therapy and eventually for combining radioisotope therapy and photodynamic therapy.

Comparison of basic features of proton and helium ion pencil beams in water using GATE

PURPOSE

The aim of this study was to investigate the basic features of helium ions for their possible application in advanced radiotherapy and to benchmark them against protons, the current particle of choice in the low linear energy transfer (LET) range.

MATERIAL AND METHODS

Geant4 Application for Emission Tomography (GATE) simulations were performed with beams of 1x10(7) monodirectional particles traversing a water phantom. Initial energies ranged from 50 to 250 MeV per nucleon (MeV/A). The following parameters were evaluated: particle range at the distal 80% of maximum energy deposition (E(max)), width of the Bragg peak (BP) at 60% of E(max), and dose fall-off width between 80% and 20% of E(max) for longitudinal spectra. In addition the fragmentation tail was quantified in terms of length, percental energy deposition, and contributing particles. For each energy lateral profiles were registered along the beam axis and the FWHM at four different depths was extracted. Besides the comparison of parameters between the two particle types, results were also compared to data in the literature.

RESULTS

As expected, the position of the BP as a function of initial kinetic energy showed similar values for protons and helium ions, with deviations smaller than 1.3%. The quantitative results of the Monte Carlo (MC) study showed less range straggling effects and smaller lateral deflections for helium ions compared to protons for the investigated energy range. On average, an about 56% reduction of the width of the BP and a 48% reduction of the dose fall-off was observed for helium ions compared to protons. Both the width of the BP and the dose fall-off width as a function of particle range or energy showed an almost linear increase with increasing energy. The tail length increased from 55.9 mm to 592.7 mm and the deposited energy increased from 0.5% to 7.3% for energies between 90 and 250 MeV/A. Lateral profiles of helium ions were about 52% narrower than those of protons.

CONCLUSIONS

Due to their mass and charge helium ions distinguish themselves from protons in reduced range straggling effects, smaller lateral deflections, and a fragmentation tail. The MC based comprehensive data set for 21 clinically relevant energies can be used to create look-up tables for semi-analytical pencil beam model for helium ions.