Helium ions at the heidelberg ion beam therapy center: comparisons between FLUKA Monte Carlo code predictions and dosimetric measurements.

In the field of particle therapy helium ion beams could offer an alternative for radiotherapy treatments, owing to their interesting physical and biological properties intermediate between protons and carbon ions. We present in this work the comparisons and validations of the Monte Carlo FLUKA code against in-depth dosimetric measurements acquired at the Heidelberg Ion Beam Therapy Center (HIT). Depth dose distributions in water with and without ripple filter, lateral profiles at different depths in water and a spread-out Bragg peak were investigated. After experimentally-driven tuning of the less known initial beam characteristics in vacuum (beam lateral size and momentum spread) and simulation parameters (water ionization potential), comparisons of depth dose distributions were performed between simulations and measurements, which showed overall good agreement with range differences below 0.1 mm and dose-weighted average dose-differences below 2.3% throughout the entire energy range. Comparisons of lateral dose profiles showed differences in full-width-half-maximum lower than 0.7 mm. Measurements of the spread-out Bragg peak indicated differences with simulations below 1% in the high dose regions and 3% in all other regions, with a range difference less than 0.5 mm. Despite the promising results, some discrepancies between simulations and measurements were observed, particularly at high energies. These differences were attributed to an underestimation of dose contributions from secondary particles at large angles, as seen in a triple Gaussian parametrization of the lateral profiles along the depth. However, the results allowed us to validate FLUKA simulations against measurements, confirming its suitability for 4He ion beam modeling in preparation of clinical establishment at HIT. Future activities building on this work will include treatment plan comparisons using validated biological models between proton and helium ions, either within a Monte Carlo treatment planning engine based on the same FLUKA code, or an independent analytical planning system fed with a validated database of inputs calculated with FLUKA.

Noise evaluation of Compton camera imaging for proton therapy.

Compton Cameras emerged as an alternative for real-time dose monitoring techniques for Particle Therapy (PT), based on the detection of prompt-gammas. As a consequence of the Compton scattering process, the gamma origin point can be restricted onto the surface of a cone (Compton cone). Through image reconstruction techniques, the distribution of the gamma emitters can be estimated, using cone-surfaces backprojections of the Compton cones through the image space, along with more sophisticated statistical methods to improve the image quality. To calculate the Compton cone required for image reconstruction, either two interactions, the last being photoelectric absorption, or three scatter interactions are needed. Because of the high energy of the photons in PT the first option might not be adequate, as the photon is not absorbed in general. However, the second option is less efficient. That is the reason to resort to spectral reconstructions, where the incoming γ energy is considered as a variable in the reconstruction inverse problem. Jointly with prompt gamma, secondary neutrons and scattered photons, not strongly correlated with the dose map, can also reach the imaging detector and produce false events. These events deteriorate the image quality. Also, high intensity beams can produce particle accumulation in the camera, which lead to an increase of random coincidences, meaning events which gather measurements from different incoming particles. The noise scenario is expected to be different if double or triple events are used, and consequently, the reconstructed images can be affected differently by spurious data. The aim of the present work is to study the effect of false events in the reconstructed image, evaluating their impact in the determination of the beam particle ranges. A simulation study that includes misidentified events (neutrons and random coincidences) in the final image of a Compton Telescope for PT monitoring is presented. The complete chain of detection, from the beam particle entering a phantom to the event classification, is simulated using FLUKA. The range determination is later estimated from the reconstructed image obtained from a two and three-event algorithm based on Maximum Likelihood Expectation Maximization. The neutron background and random coincidences due to a therapeutic-like time structure are analyzed for mono-energetic proton beams. The time structure of the beam is included in the simulations, which will affect the rate of particles entering the detector.

Distributions of secondary particles in proton and carbon-ion therapy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes.

Monte Carlo simulations play a crucial role for in-vivo treatment monitoring based on PET and prompt gamma imaging in proton and carbon-ion therapies. The accuracy of the nuclear fragmentation models implemented in these codes might affect the quality of the treatment verification. In this paper, we investigate the nuclear models implemented in GATE/Geant4 and FLUKA by comparing the angular and energy distributions of secondary particles exiting a homogeneous target of PMMA. Comparison results were restricted to fragmentation of (16)O and (12)C. Despite the very simple target and set-up, substantial discrepancies were observed between the two codes. For instance, the number of high energy (>1 MeV) prompt gammas exiting the target was about twice as large with GATE/Geant4 than with FLUKA both for proton and carbon ion beams. Such differences were not observed for the predicted annihilation photon production yields, for which ratios of 1.09 and 1.20 were obtained between GATE and FLUKA for the proton beam and the carbon ion beam, respectively. For neutrons and protons, discrepancies from 14% (exiting protons-carbon ion beam) to 57% (exiting neutrons-proton beam) have been identified in production yields as well as in the energy spectra for neutrons.

Diabetes-specific variables associated with quality of life changes in young diabetic people: the type 1 diabetes Registry of Turin (Italy).

BACKGROUND AND AIMS: Type 1 diabetes (T1DM) affects young people during the most active years of their life. Our aim was to assess quality of life (QoL) and associated variables in a large cohort of adults with childhood-onset and adult-onset T1DM. METHODS: A cohort of adult patients (18 years and older) from the T1DM Registry of Turin, Italy, was recruited. Clinical characteristics and Diabetes QoL (DQOL) questionnaire were assessed by standardized procedures. RESULTS: 310 adults completed the questionnaire. Age and diabetes duration at assessment (mean ± SD) were 32.8 ± 7.3 years and 17.3 ± 6.3 years, respectively. DQOL and its subscores were in the lower quartiles of their distributions, indicating a good level of QoL. However, scores were significantly higher in females than in males, particularly for the subscale of diabetes-related worries. In multivariate analysis, lower QoL was independently associated with female sex (ß = 1.07, 95% CI 1.03-1.11, p = 0.003), higher age at onset (ß = 1.03, 1.00-1.05, p = 0.009), lower schooling (ß = 1.05, 1.00-1.09, p = 0.02), higher fasting plasma glucose (ß = 1.03, 1.01-1.05, p = 0.008), daily SMBG >4 (ß = 1.06, 1.01-1.10, p = 0.01), severe hypoglycemia over the last year (ß = 1.06, 1.01-1.11, p = 0.02), lower numbers of diabetologic visits (ß = 1.07, 1.01-1.13, p = 0.02) and hypertension (ß = 1.06, 1.02-1.10, p = 0.005). Autonomic neuropathy was associated with diabetes impact. Female sex (ß = 4.36, 2.43-7.83) and daily SMBG >4 (ß = 3.77, 1.72-8.30) were independently associated with worst level and CSII with better level (ß = 0.22, 0.07-0.68) of diabetes-related worries. CONCLUSIONS: The impact of T1DM on QoL may depend on demographic, metabolic control-related variables, presence of complications and insulin delivery modality.

The ALBA project: an evaluation of needs, management, fears of Italian young patients with type 1 diabetes in a school setting and an evaluation of parents' and teachers' perceptions.

OBJECTIVE: To determine how Italian parents and school personnel of 6-13-year-old children with type 1 diabetes (T1D) manage during school hours, including insulin administration, management of hypoglycemia, and glucagon use. A further aim was an investigation into the responsibilities and training of school personnel regarding diabetes. RESEARCH DESIGN AND METHODS: After an initial qualitative phase, semi-structured questionnaires were completed by a sample of parents and teachers. RESULTS: 220 parent and 52 teacher questionnaires were completed. 43.6% of parents said diabetes had negatively influenced school activities. Children either self-administer insulin, or have help from a parent, since there is very rarely a nurse present (3.6%) or a teacher who will take responsibility for the treatment (2.9%). Most parents (55.9%) stated either that the school had no refrigerator to store glucagon or that they did not know if the school was so equipped. A small percentage of teachers considered their schools to be equipped to manage an emergency (23%) and said they would use glucagon directly in an emergency (14.9%). Only 40.4% of teachers said that they had received any specific training. CONCLUSIONS: The study shows that people who are not directly involved have superficial knowledge of the different aspects of diabetes, even though no parents reported episodes of neglect/incorrect management. There is no legislation which clearly defines the role of the school in the care of children with T1D, and teachers are not trained to help them. Training sessions for school personnel and greater legislative clarity about the 'insulin and glucagon question' are key factors that may improve the full integration of the child with diabetes.

Benchmarking nuclear models of FLUKA and GEANT4 for carbon ion therapy.

As carbon ions, at therapeutic energies, penetrate tissue, they undergo inelastic nuclear reactions and give rise to significant yields of secondary fragment fluences. Therefore, an accurate prediction of these fluences resulting from the primary carbon interactions is necessary in the patient's body in order to precisely simulate the spatial dose distribution and the resulting biological effect. In this paper, the performance of nuclear fragmentation models of the Monte Carlo transport codes, FLUKA and GEANT4, in tissue-like media and for an energy regime relevant for therapeutic carbon ions is investigated. The ability of these Monte Carlo codes to reproduce experimental data of charge-changing cross sections and integral and differential yields of secondary charged fragments is evaluated. For the fragment yields, the main focus is on the consideration of experimental approximations and uncertainties such as the energy measurement by time-of-flight. For GEANT4, the hadronic models G4BinaryLightIonReaction and G4QMD are benchmarked together with some recently enhanced de-excitation models. For non-differential quantities, discrepancies of some tens of percent are found for both codes. For differential quantities, even larger deviations are found. Implications of these findings for the therapeutic use of carbon ions are discussed.

The FLUKA Monte Carlo code coupled with the local effect model for biological calculations in carbon ion therapy.

Clinical Monte Carlo (MC) calculations for carbon ion therapy have to provide absorbed and RBE-weighted dose. The latter is defined as the product of the dose and the relative biological effectiveness (RBE). At the GSI Helmholtzzentrum für Schwerionenforschung as well as at the Heidelberg Ion Therapy Center (HIT), the RBE values are calculated according to the local effect model (LEM). In this paper, we describe the approach followed for coupling the FLUKA MC code with the LEM and its application to dose and RBE-weighted dose calculations for a superimposition of two opposed (12)C ion fields as applied in therapeutic irradiations. The obtained results are compared with the available experimental data of CHO (Chinese hamster ovary) cell survival and the outcomes of the GSI analytical treatment planning code TRiP98. Some discrepancies have been observed between the analytical and MC calculations of absorbed physical dose profiles, which can be explained by the differences between the laterally integrated depth-dose distributions in water used as input basic data in TRiP98 and the FLUKA recalculated ones. On the other hand, taking into account the differences in the physical beam modeling, the FLUKA-based biological calculations of the CHO cell survival profiles are found in good agreement with the experimental data as well with the TRiP98 predictions. The developed approach that combines the MC transport/interaction capability with the same biological model as in the treatment planning system (TPS) will be used at HIT to support validation/improvement of both dose and RBE-weighted dose calculations performed by the analytical TPS.

In-beam PET monitoring of mono-energetic (16)O and (12)C beams: experiments and FLUKA simulations for homogeneous targets.

(16)O and (12)C ion beams will be used-besides lighter ions-for cancer treatment at the Heidelberg Ion Therapy Center (HIT), Germany. It is planned to monitor the treatment by means of in-beam positron emission tomography (PET) as it is done for therapy with (12)C beams at the experimental facility at the Gesellschaft für Schwerionenforschung (GSI), Darmstadt, Germany. To enable PET also for (16)O beams, experimental data of the beta(+)-activity created by these beams are needed. Therefore, in-beam PET measurements of the activity created by (16)O beams of various energies on targets of PMMA, water and graphite were performed at GSI for the first time. Additionally reference measurements of (12)C beams on the same target materials were done. The results of the measurements are presented. The deduction of clinically relevant results from in-beam PET data requires reliable simulations of the beta(+)-activity production, which is done presently by a dedicated code limited to (12)C beams. Because this code is not extendable to other ions in an easy way, a new code, capable of simulating the production of the beta(+)-activity by all ions of interest, is needed. Our choice is the general purpose Monte Carlo code FLUKA which was used to simulate the ion transport, the beta(+)-active isotope production, the decay, the positron annihilation and the transport of the annihilation photons. The detector response was simulated with an established software that gives the output in the same list-mode data format as in the experiment. This allows us to use the same software to reconstruct measured and simulated data, which makes comparisons easier and more reliable. The calculated activity distribution shows general good agreement with the measurements.

The FLUKA code: new developments and application to 1 GeV/n iron beams.

The modeling of ion transport and interactions in matter is a subject of growing interest, driven by the continuous increase of possible application fields. These include hadron therapy, dosimetry, and space missions, but there are also several issues involving fundamental research, accelerator physics, and cosmic ray physics, where a reliable description of heavy ion induced cascades is important. In the present work, the capabilities of the FLUKA code for ion beams will be briefly recalled and some recent developments presented. Applications of the code to the simulation of therapeutic carbon, nitrogen and oxygen ion beams, and of iron beams, which are of direct interest for space mission related experiments, will be also presented together with interesting consideration relative to the evaluation of dosimetric quantities. Both applications involve ion beams in the AGeV range.

The application of FLUKA to dosimetry and radiation therapy.

The FLUKA Monte Carlo code has been evolving over the last several decades and is now widely used for radiation shielding calculations. In order to facilitate the use of FLUKA in dosimetry and therapy applications, supporting software has been developed to allow the direct conversion of the output files from standard CT-scans directly into a voxel geometry for transport within FLUKA. Since the CT-scan information essentially contains only the electron density information over the scanned volume, one needs the specific compositions for each voxel individually. We present here the results of a simple algorithm to assign tissues in the human body to one of four categories: soft-tissue, hard-bone, trabecular-bone and porous-lung. In addition, we explore the problem of the pathlength distributions in porous media such as trabecular bone. A mechanism will be implemented within FLUKA to allow for variable multipal fixed density materials to accommodate the pathlength distributions discovered.

The FLUKA code for space applications: recent developments.

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions.