The modelling of positron emitter production and PET imaging during carbon ion therapy

Evaluating an analytical prediction algorithm of positron emitter distributions in patient data for PET monitoring of carbon ion therapy: A simulation study

In vivo treatment monitoring in ion therapy is one of the key issues for improving the treatment quality assurance procedures. Range verification is one of the most relevant and yet complex task used for in vivo treatment monitoring. In carbon ion therapy, positron emission tomography is the most widely used method. This technique exploits the β+-activity of positron emitters created by nuclear interactions between the incoming beam and the irradiated tissue. Currently, high computational efforts and time-consuming Monte Carlo simulation platforms are typically used to predict positron emitter distributions. In order to avoid time-consuming simulations, an extended filtering approach was suggested to analytically predict positron emitter profiles from depth dose distributions in carbon ion therapy. The purpose of this work is to investigate such an analytical prediction model in patient anatomies of varying complexity, highlighting its potential and the need of further improvements, especially in highly heterogeneous anatomies where many air cavities are present in the beam path. The accuracy of range verification showed a mean relative error of ∼3% and a deviation between the simulation and the prediction below 2mm for the three patient cases analysed: a brain case and two head and neck cases. Additional investigations demonstrated the region of applicability of the method for cases of patient data. The analytical method enables range verification in carbon ion therapy by replacing computing-intensive Monte Carlo simulations and thus minimize the PET monitoring burden on the clinical workflow.

Emerging technologies for cancer therapy using accelerated particles

Cancer therapy with accelerated charged particles is one of the most valuable biomedical applications of nuclear physics. The technology has vastly evolved in the past 50 years, the number of clinical centers is exponentially growing, and recent clinical results support the physics and radiobiology rationale that particles should be less toxic and more effective than conventional X-rays for many cancer patients. Charged particles are also the most mature technology for clinical translation of ultra-high dose rate (FLASH) radiotherapy. However, the fraction of patients treated with accelerated particles is still very small and the therapy is only applied to a few solid cancer indications. The growth of particle therapy strongly depends on technological innovations aiming to make the therapy cheaper, more conformal and faster. The most promising solutions to reach these goals are superconductive magnets to build compact accelerators; gantryless beam delivery; online image-guidance and adaptive therapy with the support of machine learning algorithms; and high-intensity accelerators coupled to online imaging. Large international collaborations are needed to hasten the clinical translation of the research results.

Treatment-integrated imaging, radiomics, and personalised radiotherapy: the future is at hand

Since the introduction of computed tomography for planning purposes in the 1970s, we have been observing a continuous development of different imaging methods in radiotherapy. The current achievements of imaging technologies in radiotherapy enable more than just improvement of accuracy on the planning stage. Through integrating imaging with treatment machines, they allow advanced control methods of dose delivery during the treatment. This article reviews how the integration of existing and novel forms of imaging changes radiotherapy and how these advances can allow a more individualised approach to cancer therapy. We believe that the significant challenge for the next decade is the continued integration of a range of different imaging devices into linear accelerators. These imaging modalities should show intra-fraction changes in body morphology and inter-fraction metabolic changes. As the use of these more advanced, integrated machines grows, radiotherapy delivery will become more accurate, thus resulting in better clinical outcomes: higher cure rates with fewer side effects.

Quantitative assessment of radionuclide production yields in in-beam and offline PET measurements at different proton irradiation facilities

Objective. Reliable radionuclide production yield data are a prerequisite for positron-emission-tomography (PET) based in vivo proton treatment verification. In this context, activation data acquired at two different treatment facilities with different imaging systems were analyzed to provide experimentally determined radionuclide yields in thick targets and were compared with each other to investigate the impact of the respective imaging technique. Approach. Homogeneous thick targets (PMMA, gelatine, and graphite) were irradiated with mono-energetic proton pencil-beams at two distinct energies. Material activation was measured (i) in-beam during and after beam delivery with a double-head prototype PET camera and (ii) offline shortly after beam delivery with a commercial full-ring PET/CT scanner. Integral as well as depth-resolved β +-emitter yields were determined for the dominant positron-emitting radionuclides 11C, 15O, 13N and (in-beam only) 10C. In-beam data were used to investigate the qualitative impact of different monitoring time schemes on activity depth profiles and their quantitative impact on count rates and total activity. Main results. Production yields measured with the in-beam camera were comparable to or higher compared to respective offline results. Depth profiles of radionuclide-specific yields obtained from the double-head camera showed qualitative differences to data acquired with the full-ring camera with a more convex profile shape. Considerable impact of the imaging timing scheme on the activity profile was observed for gelatine only with a range variation of up to 3.5 mm. Evaluation of the coincidence rate and the total number of observed events in the considered workflows confirmed a strongly decreasing rate in targets with a large oxygen fraction. Significance. The observed quantitative and qualitative differences between the datasets underline the importance of a thorough system commissioning. Due to the lack of reliable cross-section data, in-house phantom measurements are still considered a gold standard for careful characterization of the system response and to ensure a reliable beam range verification.

The history of ion beam therapy in Germany

The advantageous depth dose profile of ion beams together with state of the art beam delivery and treatment planning systems allow for highly conformal tumor treatments in patients. First treatments date back to 1954 at the Lawrence Berkeley Laboratory (LBL) and in Europe, ion beam therapy started in the mid-1990s at the Paul-Scherrer Institute (PSI) with protons and at the Helmholtz Center for Heavy Ion Research (GSI) with carbon ions, followed by the Heidelberg Ion Therapy Center (HIT) in Heidelberg. This review describes the historical development of ion beam therapy in Germany based on the pioneering work at LBL and in the context of simultaneous developments in other countries as well as recent developments.

Physics and biomedical challenges of cancer therapy with accelerated heavy ions

Radiotherapy should have low toxicity in the entrance channel (normal tissue) and be very effective in cell killing in the target region (tumour). In this regard, ions heavier than protons have both physical and radiobiological advantages over conventional X-rays. Carbon ions represent an excellent combination of physical and biological advantages. There are a dozen carbon-ion clinical centres in Europe and Asia, and more under construction or at the planning stage, including the first in the USA. Clinical results from Japan and Germany are promising, but a heated debate on the cost-effectiveness is ongoing in the clinical community, owing to the larger footprint and greater expense of heavy ion facilities compared with proton therapy centres. We review here the physical basis and the clinical data with carbon ions and the use of different ions, such as helium and oxygen. Research towards smaller and cheaper machines with more effective beam delivery is necessary to make particle therapy affordable. The potential of heavy ions has not been fully exploited in clinics and, rather than there being a single 'silver bullet', different particles and their combination can provide a breakthrough in radiotherapy treatments in specific cases.

Offline imaging of positron emitters induced by therapeutic helium, carbon and oxygen ion beams with a full-ring PET/CT scanner: experiments in reference targets

In vivo verification of light ion therapy based on positron-emission tomography (PET) imaging of irradiation induced patient activation relies on activity predictions from Monte-Carlo (MC) or analytical computational engines for comparison to the measurements. In order to achieve the necessary accuracy, experimental data are indispensable for the validation of the calculation models. For this we irradiated thick reference targets with mono-energetic helium, carbon and oxygen ion beams and measured the resulting material activation offline with a commercial full-ring PET/CT scanner located nearby the treatment room. Acquired PET data were analysed over time to separate the activity contribution of different radionuclides. Determined production yields were compared to published findings obtained from in-beam activation measurements with a limited-angle double-head PET camera. In addition, we investigated the time-dependence of the measured radionuclide-specific contributions and of the distal activity range, as well as the lateral spread of the activity signal as a function of beam penetration depth. We present radionuclide-specific depth-resolved activity distributions and production yields for the radionuclides [Formula: see text], [Formula: see text] and [Formula: see text], dominating irradiation-induced patient activation. We observe systematically lower production yields with a ratio between the dual-head and our full-ring PET measurements of, on average, 1.7 and 1.3 for the oxygen and carbon beam irradiations, and 1.7 (2.1) for the high (low) energy helium beam irradiations. Findings on the temporal development of the activity range confirm the expectation, with the oxygen beam induced signal being the most sensitive scenario. The experimental data reported in this work, acquired with a state-of-the-art full ring PET scanner, provide a comprehensive and consistent basis for the benchmarking of PET signal calculation engines. In particular, they can support a fine-tuning of the underlying physics models used by the respective implementation and therefore improve the accuracy of PET-based therapy verifications at current and future treatment facilities.

Measurement of PET isotope production cross sections for protons and carbon ions on carbon and oxygen targets for applications in particle therapy range verification

Measured cross sections for the production of the PET isotopes [Formula: see text], [Formula: see text] and [Formula: see text] from carbon and oxygen targets induced by protons (40-220 [Formula: see text]) and carbon ions (65-430 [Formula: see text]) are presented. These data were obtained via activation measurements of irradiated graphite and beryllium oxide targets using a set of three scintillators coupled by a coincidence logic. The measured cross sections are relevant for the PET particle range verification method where accurate predictions of the [Formula: see text] emitter distribution produced by therapeutic beams in the patient tissue are required. The presented dataset is useful for validation and optimization of the nuclear reaction models within Monte Carlo transport codes. For protons the agreement of a radiation transport calculation using the measured cross sections with a thick target PET measurement is demonstrated.

Comparative study of alternative Geant4 hadronic ion inelastic physics models for prediction of positron-emitting radionuclide production in carbon and oxygen ion therapy

The distribution of fragmentation products predicted by Monte Carlo simulations of heavy ion therapy depend on the hadronic physics model chosen in the simulation. This work aims to evaluate three alternative hadronic inelastic fragmentation physics options available in the Geant4 Monte Carlo radiation physics simulation framework to determine which model most accurately predicts the production of positron-emitting fragmentation products observable using in-beam PET imaging. Fragment distributions obtained with the BIC, QMD, and INCL + + physics models in Geant4 version 10.2.p03 are compared to experimental data obtained at the HIMAC heavy-ion treatment facility at NIRS in Chiba, Japan. For both simulations and experiments, monoenergetic beams are applied to three different block phantoms composed of gelatin, poly(methyl methacrylate) and polyethylene. The yields of the positron-emitting nuclei ¹¹C, ¹⁰C and ¹⁵O obtained from simulations conducted with each model are compared to the experimental yields estimated by fitting a multi-exponential radioactive decay model to dynamic PET images using the normalised mean square error metric in the entrance, build up/Bragg peak and tail regions. Significant differences in positron-emitting fragment yield are observed among the three physics models with the best overall fit to experimental ¹²C and ¹⁶O beam measurements obtained with the BIC physics model.

Monte Carlo investigation of the characteristics of radioactive beams for heavy ion therapy

This work presents a simulation study evaluating relative biological effectiveness at 10% survival fraction (RBE10) of several different positron-emitting radionuclides in heavy ion treatment systems, and comparing these to the RBE10s of their non-radioactive counterparts. RBE10 is evaluated as a function of depth for three positron-emitting radioactive ion beams (¹⁰C, ¹¹C and ¹⁵O) and two stable ion beams (¹²C and ¹⁶O) using the modified microdosimetric kinetic model (MKM) in a heterogeneous skull phantom subject to a rectangular 50 mm × 50 mm × 60 mm spread out Bragg peak. We demonstrate that the RBE10 of the positron-emitting radioactive beams is almost identical to the corresponding stable isotopes. The potential improvement in PET quality assurance image quality which is obtained when using radioactive beams is evaluated by comparing the signal to background ratios of positron annihilations at different intra- and post-irradiation time points. Finally, the incidental dose to the patient resulting from the use of radioactive beams is also quantified and shown to be negligible.

In vivo range verification in particle therapy

Exploitation of the full potential offered by ion beams in clinical practice is still hampered by several sources of treatment uncertainties, particularly related to the limitations of our ability to locate the position of the Bragg peak in the tumor. To this end, several efforts are ongoing to improve the characterization of patient position, anatomy, and tissue stopping power properties prior to treatment as well as to enable in vivo verification of the actual dose delivery, or at least beam range, during or shortly after treatment. This contribution critically reviews methods under development or clinical testing for verification of ion therapy, based on pretreatment range and tissue probing as well as the detection of secondary emissions or physiological changes during and after treatment, trying to disentangle approaches of general applicability from those more specific to certain anatomical locations. Moreover, it discusses future directions, which could benefit from an integration of multiple modalities or address novel exploitation of the measurable signals for biologically adapted therapy.

Experimental validation of the FLUKA Monte Carlo code for dose and [Formula: see text]-emitter predictions of radioactive ion beams

In the context of hadrontherapy, whilst ions are capable of effectively destroying radio resistant, deep seated tumors, their treatment localization must be well assessed to ensure the sparing of surrounding healthy tissue and treatment effectiveness. Thus, range verification techniques, such as online positron-emission-tomography (PET) imaging, hold great potential in clinical practice, providing information on the in vivo beam range and consequent tumor targeting. Furthermore, [Formula: see text] emitting radioactive ions can be an asset in online PET imaging, depending on their half-life, compared to their stable counterparts. It is expected that using these radioactive ions the signal obtained by a PET apparatus during beam delivery will be greatly increased, and exhibit a better correlation to the Bragg Peak. To this end, FLUKA Monte Carlo particle transport and interaction code was used to evaluate, in terms of annihilation events at rest and dose, the figure of merit in using [Formula: see text] emitter, radioactive ion beams (RI [Formula: see text]). For this purpose, the simulation results were compared with experimental data obtained with an openPET prototype in various online PET acquisitions at the Heavy Ion Medical Accelerator in Chiba (HIMAC), in collaboration with colleagues from the National Institute of Radiological Sciences' (NIRS) Imaging Physics Team. The dosimetry performance evaluation with FLUKA benefits from its recent developments in fragmentation production models. The present work estimated that irradiations with RI [Formula: see text], produced via projectile fragmentation and their signal acquisition with state-of-the-art PET scanner, lead to nearly a factor of two more accurate definition of the signals' peak position. In addition to its more advantageous distribution shape, it was observed at least an order magnitude higher signal acquired from ¹¹C and ¹⁵O irradiations, with respect to their stable counterparts.

An overview of recent developments in FLUKA PET tools

The new developments of the FLUKA Positron-Emission-Tomography (PET) tools are detailed. FLUKA is a fully integrated Monte Carlo (MC) particle transport code, used for an extended range of applications, including Medical Physics. Recently, it provided the medical community with dedicated simulation tools for clinical applications, including the PET simulation package. PET is a well-established imaging technique in nuclear medicine, and a promising method for clinical in vivo treatment verification in hadrontherapy. The application of clinically established PET scanners to new irradiation environments such as hadrontherapy requires further experimental and theoretical research to which MC simulations could be applied. The FLUKA PET tools, besides featuring PET scanner models in its library, allow the configuration of new PET prototypes via the FLUKA Graphical User Interface (GUI) Flair. Both the beam time structure and scan time can be specified by the user, reproducing PET acquisitions in time, in a particle therapy scenario. Furthermore, different scoring routines allow the analysis of single and coincident events, and identification of parent isotopes generating annihilation events. Two reconstruction codes are currently supported: the Filtered Back-Projection (FBP) and Maximum-Likelihood Expectation Maximization (MLEM), the latter embedded in the tools. Compatibility with other reconstruction frameworks is also possible. The FLUKA PET tools package has been successfully tested for different detectors and scenarios, including conventional functional PET applications and in beam PET, either using radioactive sources, or simulating hadron beam irradiations. The results obtained so far confirm the FLUKA PET tools suitability to perform PET simulations in R&D environment.

Towards offline PET monitoring of proton therapy at MedAustron

The characteristic depth-dose profile of protons traveling through material is the main advantage of proton therapy over conventional radiotherapy with photons or electrons. However, uncertainties regarding the range of the protons in human tissue prevent to exploit the full potential of proton therapy. Therefore, a non-invasive in-vivo dose monitoring is desired. At the ion beam center MedAustron in Wiener Neustadt/Austria, patient treatment with proton beams started in December 2016. A PET/CT is available in close vicinity of the treatment rooms, exclusively dedicated to offline PET monitoring directly after the therapeutic irradiation. Preparations for a patient study include workflow tests under realistic clinical conditions using two different phantoms, irradiated with protons prior to the scan in the PET/CT. GATE simulations of the C-11 production are used as basis for the prediction of the PET measurement. We present results from the workflow tests in comparison with simulation results, and by this, we demonstrate the applicability of the PET monitoring at the MedAustron facility.

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.

Short-lived positron emitters in beam-on PET imaging during proton therapy

The only method for in vivo dose delivery verification in proton beam radiotherapy in clinical use today is positron emission tomography (PET) of the positron emitters produced in the patient during irradiation. PET imaging while the beam is on (so called beam-on PET) is an attractive option, providing the largest number of counts, the least biological washout and the fastest feedback. In this implementation, all nuclides, independent of their half-life, will contribute. As a first step towards assessing the relevance of short-lived nuclides (half-life shorter than that of (10)C, T1/2  =  19 s) for in vivo dose delivery verification using beam-on PET, we measured their production in the stopping of 55 MeV protons in water, carbon, phosphorus and calcium The most copiously produced short-lived nuclides and their production rates relative to the relevant long-lived nuclides are: (12)N (T1/2  =  11 ms) on carbon (9% of (11)C), (29)P (T1/2  =  4.1 s) on phosphorus (20% of (30)P) and (38m)K (T1/2  =  0.92 s) on calcium (113% of (38g)K). No short-lived nuclides are produced on oxygen. The number of decays integrated from the start of an irradiation as a function of time during the irradiation of PMMA and 4 tissue materials has been determined. For (carbon-rich) adipose tissue, (12)N dominates up to 70 s. On bone tissue, (12)N dominates over (15)O during the first 8-15 s (depending on carbon-to-oxygen ratio). The short-lived nuclides created on phosphorus and calcium provide 2.5 times more beam-on PET counts than the long-lived ones produced on these elements during a 70 s irradiation. From the estimated number of (12)N PET counts, we conclude that, for any tissue, (12)N PET imaging potentially provides equal to superior proton range information compared to prompt gamma imaging with an optimized knife-edge slit camera. The practical implementation of (12)N PET imaging is discussed.

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.

Comparing the biological washout of β+-activity induced in mice brain after 12C-ion and proton irradiation

In clinical ion beam therapy, protons as well as heavier ions such as carbon are used for treatment. For protons, β(+)-emitters are only induced by fragmentation reactions in the target (target fragmentation), whereas for heavy ions, they are additionally induced by fragmentations of the projectile (further referred to as autoactivation). An approach utilizing these processes for treatment verification, by comparing measured Positron Emission Tomography (PET) data to predictions from Monte Carlo simulations, has already been clinically implemented. For an accurate simulation, it is important to consider the biological washout of β(+)-emitters due to vital functions. To date, mathematical expressions for washout have mainly been determined by using radioactive beams of (10)C- and (11)C-ions, both β(+)-emitters, to enhance the counting statistics in the irradiated area. Still, the question of how the choice of projectile (autoactivating or non-autoactivating) influences the washout coefficients, has not been addressed. In this context, an experiment was carried out at the Heidelberg Ion Beam Therapy Center with the purpose of directly comparing irradiation-induced biological washout coefficients in mice for protons and (12)C-ions. To this aim, mice were irradiated in the brain region with protons and (12)C-ions and measured after irradiation with a PET/CT scanner (Siemens Biograph mCT). After an appropriate waiting time, the mice were sacrificed, then irradiated and measured again under similar conditions. The resulting data were processed and fitted numerically to deduce the main washout parameters. Despite the very low PET counting statistics, a consistent difference could be identified between (12)C-ion and proton irradiated mice, with the (12)C data being described best by a two component fit with a combined medium and slow washout fraction of 0.50 ± 0.05 and the proton mice data being described best by a one component fit with only one (slow) washout fraction of 0.73 ± 0.06.

Systematic analysis on the achievable accuracy of PT-PET through automated evaluation techniques

INTRODUCTION

Particle Therapy Positron Emission Tomography (PT-PET) is currently the only clinically applied method for in vivo verification of ion-beam radiotherapy during or close in time to the treatment. Since a direct deduction of the delivered dose from the measured activity is not feasible, images are compared to a reference distribution. The achievable accuracy of two image analysis approaches was investigated by means of reproducible phantom benchmark tests. This is an objective method that excludes patient related factors of influence.

MATERIAL AND METHODS

Two types of phantoms were designed to produce well defined deviations in the activity distributions. Pure range differences were simulated using the first phantom type while the other emulated cavity structures. The phantoms were irradiated with (12)C-ions. PT-PET measurements were performed by means of a camera system installed at the beamline. Different measurement time scenarios were investigated, assuming a PET scanner directly at the irradiation site or placed within the treatment room. The images were analyzed by means of the Pearson Correlation Coefficient (PCC) and a range calculation algorithm combined with a dedicated cavity filling detection method.

RESULTS

Range differences could be measured with an error of less than 2 mm. The range comparison algorithm yielded slightly better results than the PCC method. The filling of a cavity structure could be safely detected if its inner diameter was at least 5 mm.

CONCLUSION

Both approaches evaluate the PT-PET data in an objective way and deliver promising results for in-beam and in-room PET for clinical realistic dose rates.

Automated evaluation of setup errors in carbon ion therapy using PET: feasibility study

PURPOSE

To investigate the possibility of detecting patient mispositioning in carbon-ion therapy with particle therapy positron emission tomography (PET) in an automated image registration based manner.

METHODS

Tumors in the head and neck (H&N), pelvic, lung, and brain region were investigated. Biologically optimized carbon ion treatment plans were created with TRiP98. From these treatment plans, the reference β(+)-activity distributions were calculated using a Monte Carlo simulation. Setup errors were simulated by shifting or rotating the computed tomography (CT). The expected β(+) activity was calculated for each plan with shifts. Finally, the reference particle therapy PET images were compared to the "shifted" β(+)-activity distribution simulations using the Pearson's correlation coefficient (PCC). To account for different PET monitoring options the inbeam PET was compared to three different inroom scenarios. Additionally, the dosimetric effects of the CT misalignments were investigated.

RESULTS

The automated PCC detection of patient mispositioning was possible in the investigated indications for cranio-caudal shifts of 4 mm and more, except for prostate tumors. In the rather homogeneous pelvic region, the generated β(+)-activity distribution of the reference and compared PET image were too much alike. Thus, setup errors in this region could not be detected. Regarding lung lesions the detection strongly depended on the exact tumor location: in the center of the lung tumor misalignments could be detected down to 2 mm shifts while resolving shifts of tumors close to the thoracic wall was more challenging. Rotational shifts in the H&N and lung region of +6° and more could be detected using inroom PET and partly using inbeam PET. Comparing inroom PET to inbeam PET no obvious trend was found. However, among the inroom scenarios a longer measurement time was found to be advantageous.

CONCLUSIONS

This study scopes the use of various particle therapy PET verification techniques in four indications. The automated detection of patients' setup errors was investigated in a broad accumulation of data sets. The evaluation of introduced setup errors is performed automatically, which is of utmost importance to introduce highly required particle therapy monitoring devices into the clinical routine.

Future radiation therapy: photons, protons and particles

Radiation therapy plays a critical role in the current management of cancer patients. The most common linear accelerator-based treatment device delivers photons of radiation. In an ever more precise fashion, state-of-the-art technology has recently allowed for both modulation of the radiation beam and imaging for this treatment delivery. This has resulted in better patient outcome with far fewer side effects than were achieved even a decade ago. Recently, a push has begun for proton therapy, which may have clinical advantage in select indications, although significant limitations for these devices have become apparent. In addition, currently, heavy particle therapy has been touted as a potential means to improve cancer patient outcomes. This article will highlight current benefits and drawbacks to modern radiation therapy and speculate on future tools that will likely dramatically improve radiation oncology.

Dependence of simulated positron emitter yields in ion beam cancer therapy on modeling nuclear fragmentation

In ion beam cancer therapy, range verification in patients using positron emission tomography (PET) requires the comparison of measured with simulated positron emitter yields. We found that (1) changes in modeling nuclear interactions strongly affected the positron emitter yields and that (2) Monte Carlo simulations with SHIELD-HIT10Areasonably matched the most abundant PET isotopes (11)C and (15)O. We observed an ion-energy (i.e., depth) dependence of the agreement between SHIELD-HIT10Aand measurement. Improved modeling requires more accurate measurements of cross-section values.

Radiation dose detection by imaging response in biological targets

Imaging was one of the earliest techniques to quantify radiation dose. While films and active fluorescent detectors are still commonly used in physical dosimetry, biological imaging is emerging as a new method to visualize and quantify radiation dose in biological targets. Methods for biological imaging are normally based on molecular fluorescent probes, labeling chromatin-conjugated molecules or specific repair proteins. Examples are chromatin-binding coumarin compounds, which become fluorescent under irradiation, or the H2AX histone, which is rapidly phosphorylated at sites of DNA double-strand breaks and can be visualized by immunostaining. Many other DNA repair proteins can be expressed with fluorescent targets, such as green fluorescent protein, thus becoming visible for dose estimation in vivo. The possibility to visualize radiation damage in living biological targets is particularly important for repair kinetic studies, for estimating individual radiation response, and for remote control of living samples exposed to radiation, for instance in robotic space missions. In vivo dose monitoring in particle therapy exploits the production of positron emitters by nuclear interaction of the incident beam in the patient's body. Positron emission tomography (PET) can then be used to visualize and quantify the particle dose in the patient, and it can in principle also be used for radiotherapy with high-energy X rays. Alternatively, prompt γ rays or scattered secondary particles are under study for in vivo dosimetry of ion beams in therapy.

Development of activity pencil beam algorithm using measured distribution data of positron emitter nuclei generated by proton irradiation of targets containing (12)C, (16)O, and (40)Ca nuclei in preparation of clinical application

PURPOSE

The purpose of this study is to develop a new calculation algorithm that is satisfactory in terms of the requirements for both accuracy and calculation time for a simulation of imaging of the proton-irradiated volume in a patient body in clinical proton therapy.

METHODS

The activity pencil beam algorithm (APB algorithm), which is a new technique to apply the pencil beam algorithm generally used for proton dose calculations in proton therapy to the calculation of activity distributions, was developed as a calculation algorithm of the activity distributions formed by positron emitter nuclei generated from target nuclear fragment reactions. In the APB algorithm, activity distributions are calculated using an activity pencil beam kernel. In addition, the activity pencil beam kernel is constructed using measured activity distributions in the depth direction and calculations in the lateral direction. (12)C, (16)O, and (40)Ca nuclei were determined as the major target nuclei that constitute a human body that are of relevance for calculation of activity distributions. In this study, "virtual positron emitter nuclei" was defined as the integral yield of various positron emitter nuclei generated from each target nucleus by target nuclear fragment reactions with irradiated proton beam. Compounds, namely, polyethylene, water (including some gelatin) and calcium oxide, which contain plenty of the target nuclei, were irradiated using a proton beam. In addition, depth activity distributions of virtual positron emitter nuclei generated in each compound from target nuclear fragment reactions were measured using a beam ON-LINE PET system mounted a rotating gantry port (BOLPs-RGp). The measured activity distributions depend on depth or, in other words, energy. The irradiated proton beam energies were 138, 179, and 223 MeV, and measurement time was about 5 h until the measured activity reached the background level. Furthermore, the activity pencil beam data were made using the activity pencil beam kernel, which was composed of the measured depth data and the lateral data including multiple Coulomb scattering approximated by the Gaussian function, and were used for calculating activity distributions.

RESULTS

The data of measured depth activity distributions for every target nucleus by proton beam energy were obtained using BOLPs-RGp. The form of the depth activity distribution was verified, and the data were made in consideration of the time-dependent change of the form. Time dependence of an activity distribution form could be represented by two half-lives. Gaussian form of the lateral distribution of the activity pencil beam kernel was decided by the effect of multiple Coulomb scattering. Thus, the data of activity pencil beam involving time dependence could be obtained in this study.

CONCLUSIONS

The simulation of imaging of the proton-irradiated volume in a patient body using target nuclear fragment reactions was feasible with the developed APB algorithm taking time dependence into account. With the use of the APB algorithm, it was suggested that a system of simulation of activity distributions that has levels of both accuracy and calculation time appropriate for clinical use can be constructed.

Basics of particle therapy I: physics

With the advance of modern radiation therapy technique, radiation dose conformation and dose distribution have improved dramatically. However, the progress does not completely fulfill the goal of cancer treatment such as improved local control or survival. The discordances with the clinical results are from the biophysical nature of photon, which is the main source of radiation therapy in current field, with the lower linear energy transfer to the target. As part of a natural progression, there recently has been a resurgence of interest in particle therapy, specifically using heavy charged particles, because these kinds of radiations serve theoretical advantages in both biological and physical aspects. The Korean government is to set up a heavy charged particle facility in Korea Institute of Radiological & Medical Sciences. This review introduces some of the elementary physics of the various particles for the sake of Korean radiation oncologists' interest.

The reliability of proton-nuclear interaction cross-section data to predict proton-induced PET images in proton therapy

In vivo PET range verification relies on the comparison of measured and simulated activity distributions. The accuracy of the simulated distribution depends on the accuracy of the Monte Carlo code, which is in turn dependent on the accuracy of the available cross-section data for β(+) isotope production. We have explored different cross-section data available in the literature for the main reaction channels ((16)O(p,pn)(15)O, (12)C(p,pn)(11)C and (16)O(p,3p3n)(11)C) contributing to the production of β(+) isotopes by proton beams in patients. Available experimental and theoretical values were implemented in the simulation and compared with measured PET images obtained with a high-resolution PET scanner. Each reaction channel was studied independently. A phantom with three different materials was built, two of them with high carbon or oxygen concentration and a third one with average soft tissue composition. Monoenergetic and SOBP field irradiations of the phantom were accomplished and measured PET images were compared with simulation results. Different cross-section values for the tissue-equivalent material lead to range differences below 1 mm when a 5 min scan time was employed and close to 5 mm differences for a 30 min scan time with 15 min delay between irradiation and scan (a typical off-line protocol). The results presented here emphasize the need of more accurate measurement of the cross-section values of the reaction channels contributing to the production of PET isotopes by proton beams before this in vivo range verification method can achieve mm accuracy.

Heavy charged particle radiobiology: using enhanced biological effectiveness and improved beam focusing to advance cancer therapy

Ionizing radiation causes many types of DNA damage, including base damage and single- and double-strand breaks. Photons, including X-rays and γ-rays, are the most widely used type of ionizing radiation in radiobiology experiments, and in radiation cancer therapy. Charged particles, including protons and carbon ions, are seeing increased use as an alternative therapeutic modality. Although the facilities needed to produce high energy charged particle beams are more costly than photon facilities, particle therapy has shown improved cancer survival rates, reflecting more highly focused dose distributions and more severe DNA damage to tumor cells. Despite early successes of charged particle radiotherapy, there is room for further improvement, and much remains to be learned about normal and cancer cell responses to charged particle radiation.

The impact of uncertainties in the CT conversion algorithm when predicting proton beam ranges in patients from dose and PET-activity distributions

The advantages of a finite range of proton beams can only be partly exploited in radiation therapy unless the range can be predicted in patient anatomy with <2 mm accuracy (for non-moving targets). Monte Carlo dose calculation aims at 1-2 mm accuracy in dose prediction, and proton-induced PET imaging aims at ∼2 mm accuracy in range verification. The latter is done using Monte Carlo predicted PET images. Monte Carlo methods are based on CT images to describe patient anatomy. The dose calculation algorithm and the CT resolution/artifacts might affect dose calculation accuracy. Additionally, when using Monte Carlo for PET range verification, the biological decay model and the cross sections for positron emitter production affect predicted PET images. The goal of this work is to study the effect of uncertainties in the CT conversion on the proton beam range predicted by Monte Carlo dose calculations and proton-induced PET signals. Conversion schemes to assign density and elemental composition based on a CT image of the patient define a unique Hounsfield unit (HU) to tissue parameters relationship. Uncertainties are introduced because there is no unique relationship between HU and tissue parameters. In this work, different conversion schemes based on a stoichiometric calibration method as well as different numbers of tissue bins were considered in three head and neck patients. For Monte Carlo dose calculation, the results show close to zero (<0.5 mm) differences in range using different conversion schemes. Further, a reduction of the number of bins used to define individual tissues down to 13 did not affect the accuracy. In the case of simulated PET images we found a more pronounced sensitivity on the CT conversion scheme with a mean fall-off position variation of about 1 mm. We conclude that proton dose distributions based on Monte Carlo calculation are only slightly affected by the uncertainty on density and elemental composition introduced by unique assignment to each HU if a stoichiometric calibration is used. Calculated PET images used for range verification are more sensitive to conversion uncertainties causing an intrinsic limitation due to CT conversion alone of at least 1 mm.

On the effectiveness of ion range determination from in-beam PET data

At present, in-beam positron emission tomography (PET) is the only method for in vivo and in situ range verification in ion therapy. At the GSI Helmholtzzentrum für Schwerionenforschung GmbH (GSI) Darmstadt, Germany, a unique in-beam PET installation has been operated from 1997 until the shut down of the carbon ion therapy facility in 2008. Therapeutic irradiation by means of (12)C ion beams of more than 400 patients have been monitored. In this paper a first quantitative study on the accuracy of the in-beam PET method to detect range deviations between planned and applied treatment in clinically relevant situations using simulations based on clinical data is presented. Patient treatment plans were used for performing simulations of positron emitter distributions. For each patient a range difference of + or - 6 mm in water was applied and compared to simulations without any changes. The comparisons were performed manually by six experienced evaluators for data of 81 patients. The number of patients required for the study was calculated using the outcome of a pilot study. The results indicate a sensitivity of (91 + or - 3)% and a specificity of (96 + or - 2)% for detecting an overrange, a reduced range is recognized with a sensitivity of (92 + or - 3)% and a specificity of (96 + or - 2)%. The positive and the negative predictive value of this method are 94% and 87%, respectively. The interobserver coefficient of variation is between 3 and 8%. The in-beam PET method demonstrated a high sensitivity and specificity for the detection of range deviations. As the range is a most indicative factor of deviations in the dose delivery, the promising results shown in this paper confirm the in-beam PET method as an appropriate tool for monitoring ion therapy.