Comparison of GEANT4 simulations with experimental data for thick al absorbers

Proton beams in medical applications deal with relatively thick targets like the human head or trunk. Therefore, relatively small differences in the total proton stopping power given, for example, by the different models provided by GEANT4 can lead to significant disagreements in the final proton energy spectra when integrated along lengthy proton trajectories. This work presents proton energy spectra obtained by GEANT4.8.2 simulations using ICRU49, Ziegler1985 and Ziegler2000 models for 19.68MeV protons passing through a number of Al absorbers with various thicknesses. The spectra were compared with the experimental data, with TRIM/SRIM2008 and MCNPX2.4.0 simulations, and with the Payne analytical solution for the transport equation in the Fokker-Plank approximation. It is shown that the MCNPX simulations reasonably reproduce well all experimental spectra. For the relatively thin targets all the methods give practically identical results but this is not the same for the thick absorbers. It should be noted that all the spectra were measured at the proton energies significantly above 2MeV, i.e., in the so-called Bethe-Bloch region. Therefore the observed disagreements in GEANT4 results, simulated with different models, are somewhat unexpected. Further studies are necessary for better understanding and definitive conclusions. © 2009 American Institute of Physics.

Comparison of the GEANT4 releases 8.2 and 9.2 in terms of a pCT reduced calibration curve

The GEANT4 simulations are essential for the development of medical tomography with proton beams pCT. In the case of thin absorbers the latest releases of GEANT4 generate very similar final spectra which agree well with the results of other popular Monte Carlo codes like TRIM/SRIM, or MCNPX. For thick absorbers, however, the disagreements became evident. In a part, these disagreements are due to the known contradictions in the NIST PSTAR and SRIM reference data. Therefore, it is interesting to compare the GEANT4 results with each other, with experiment, and with diverse code results in a reduced form, which is free from this kind of doubts. In this work such comparison is done within the Reduced Calibration Curve concept elaborated for the proton beam tomography. © 2010 IEEE.

GEANT4 tuning for pCT development

Proton beams in medical applications deal with relatively thick targets like the human head or trunk. Thus, the fidelity of proton computed tomography (pCT) simulations as a tool for proton therapy planning depends in the general case on the accuracy of results obtained for the proton interaction with thick absorbers. GEANT4 simulations of proton energy spectra after passing thick absorbers do not agree well with existing experimental data, as showed previously. Moreover, the spectra simulated for the Bethe-Bloch domain showed an unexpected sensitivity to the choice of low-energy electromagnetic models during the code execution. These observations were done with the GEANT4 version 8.2 during our simulations for pCT. This work describes in more details the simulations of the proton passage through aluminum absorbers with varied thickness. The simulations were done by modifying only the geometry in the Hadrontherapy Example, and for all available choices of the Electromagnetic Physics Models. As the most probable reasons for these effects is some specific feature in the code, or some specific implicit parameters in the GEANT4 manual, we continued our study with version 9.2 of the code. Some improvements in comparison with our previous results were obtained. The simulations were performed considering further applications for pCT development. © 2011 American Institute of Physics.

Estudo teórico da relação entre número de CT de raios-X e Stopping Power de prótons

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo da interação de prótons com alvos não homogêneos, aplicados a tomografia com feixes de prótons

Pós-graduação em Biologia Geral e Aplicada - IBB

Estudo teórico da contribuição à dose de partículas secundárias geradas por prótons de energias entre 100 e 200 MeV num Phantom de substituto de tecido MS20

The goal of this project is the reproduction, through the simulation code based on the MCNPX (Monte Carlo N-Particle eXtended) v2.50 method, of the proton beam interaction with the material, since, in proton therapy, only the particle ionization and excitation are analyzed and the occurence of nuclear interactive inelastic process are not considered. This work will help the development of studies concerning the contribution to the total dosis of secondary particles generated by nuclear interaction in proton therapy. They are: alpha particles ( ), deuterium(2H), tritium (3H), neutron (n) and helium (3He). A MS20 tissue substitute phantom was used as the target and the energy of the proton beams was within an interest range of 100 to 200MeV. With the results obtained, it was possible to generate graphics which allows the analysis of the dosis deposition relation with and without nuclear interaction, the percentage of secondary particles deposited dosis, the radial dispersion of neutrons in the material, the secondary particles multiplicity, as well as the relation between the secondary particles spectrum with the próton generated spectrum

Estudo teórico da contribuição à dose de partículas secundárias geradas por Prótons de energias entre 100 e 200 MeV num Phantom de osso cortical

In proton therapy, the deposition of secondary particles energy originated by nuclear inelastic process (n, 2H, 3H, 3He and α) has a contribution in the total dose that deserves to be discussed. In calculations of plans implemented for routine treatment, the paid dose is calculated whereas the proton loses energy by ionization and or coulomb excitement. The contribution of inelastic processes associated with nuclear reactions is not considered. There are only estimates for pure materials or simple composition (water, for example), because of the difficulty of processing targets consisting of different materials. For this project, we use the Monte Carlo method employing the code MCNPX v2.50 (Monte Carlo N-Particle eXtended) to present results of the contribution to the total dose of secondary particles. In this work, it was implemented a cylindrical phantom composed by cortical bone, for proton beams between 100 and 200 MeV. With the results obtained, it was possible to generate graphics to analyze: the dose deposition relation with and without nuclear interaction, the multiplicity and percentage of deposited dose for each secondary particle and a radial dispersion of neutrons in the material

Análise da viabilidade de simulações com feixes de prótons de 100, 150, 200 e 250 MeV em alvos heterogêneos

In the recent years, the use of proton beams in radiotherapy has been an outstanding progress (SMITH, 2006). Up to now, computed tomography (CT) is a prerequisite for treatment planning in this kind of therapy because it provides the electron density distribution required for calculation of dose and the interval of doses. However, the use of CT images for proton treatment planning ignores fundamental differences in physical interaction processes between photons and protons and is, therefore, potentially inaccurate (SADROZINSKI, 2004). Proton CT (pCT) can in principle directly measure the density distribution needed in a patient for the dose distribution (SCHULTE, et al, 2004). One important problem that should be solved is the implementation of image reconstruction algorithms. In this sense, it is necessary to know how the presence of materials with different density and composition interfere in the energy deposition by ionization and coulomb excitation, during its trajectory. The study was conducted in two stages, was used in both the program SRIM (The Stopping and Range of Ions in Matter) to perform simulations of the interaction of proton beams with pencil beam type. In the first step we used the energies in the range of 100-250 MeV (ZIEGLER, 1999). The targets were set to 50 mm in length for the beam of 100 MeV, due to its interaction with the target, and short-range, and 70 mm for 150, 200 and 250 MeV The target was composed of liquid water and a layer of 6 mm cortical bone (ICRP). It were made 9 simulations varying the position of the heterogeneity of 5 mm. In the second step the energy of 250 MeV was taken out from the simulations, due to its greater energy and less interaction. The targets were diminished to 50 mm thick to standardize the simulations. The layer of bone was divided into two equal parts and both were put in the ends of the target... (Complete abstract click electronic access below)

Electronic stopping cross sections for protons in Al2O3: an experimental and theoretical study

The electronic stopping cross section (SCS) of Al2O3 for proton beams is studied both experimentally and theoretically. The measurements are made for proton energies from 40 keV up to 1 MeV, which cover the maximum stopping region, using two experimental methods, the transmission technique at low energies (similar to 40-175 keV) and the Rutherford backscattering at high energies (approximate to 190-1000 keV). These new data reveal an increment of 16% in the SCS around the maximum stopping with respect to older measurements. The theoretical study includes electronic stopping power calculations based on the dielectric formalism and on the transport cross section (TCS) model to describe the electron excitations of Al2O3. The non-linear TCS calculations of the SCS for valence electrons together with the generalized oscillator strengths (GOS) model for the core electrons compare well with the experimental data in the whole range of energies considered.

Tracking, b-tagging and measurement of the b-jet production cross section with the ATLAS detector

The Standard Model of elementary particle physics was developed to describe the fundamental particles which constitute matter and the interactions between them. The Large Hadron Collider (LHC) at CERN in Geneva was built to solve some of the remaining open questions in the Standard Model and to explore physics beyond it, by colliding two proton beams at world-record centre-of-mass energies. The ATLAS experiment is designed to reconstruct particles and their decay products originating from these collisions. The precise reconstruction of particle trajectories plays an important role in the identification of particle jets which originate from bottom quarks (b-tagging). This thesis describes the step-wise commissioning of the ATLAS track reconstruction and b-tagging software and one of the first measurements of the b-jet production cross section in pp collisions at sqrt(s)=7 TeV with the ATLAS detector. The performance of the track reconstruction software was studied in great detail, first using data from cosmic ray showers and then collisions at sqrt(s)=900 GeV and 7 TeV. The good understanding of the track reconstruction software allowed a very early deployment of the b-tagging algorithms. First studies of these algorithms and the measurement of the b-tagging efficiency in the data are presented. They agree well with predictions from Monte Carlo simulations. The b-jet production cross section was measured with the 2010 dataset recorded by the ATLAS detector, employing muons in jets to estimate the fraction of b-jets. The measurement is in good agreement with the Standard Model predictions.

Improving Cervical Cancer Nodal Boost Radiation Therapy by Quantifying Uncertainties and Exploring Advanced Radiation Therapy Modalities

Radiation therapy for patients with intact cervical cancer is frequently delivered using primary external beam radiation therapy (EBRT) followed by two fractions of intracavitary brachytherapy (ICBT). Although the tumor is the primary radiation target, controlling microscopic disease in the lymph nodes is just as critical to patient treatment outcome. In patients where gross lymphadenopathy is discovered, an extra EBRT boost course is delivered between the two ICBT fractions. Since the nodal boost is an addendum to primary EBRT and ICBT, the prescription and delivery must be performed considering previously delivered dose. This project aims to address the major issues of this complex process for the purpose of improving treatment accuracy while increasing dose sparing to the surrounding normal tissues. Because external beam boosts to involved lymph nodes are given prior to the completion of ICBT, assumptions must be made about dose to positive lymph nodes from future implants. The first aim of this project was to quantify differences in nodal dose contribution between independent ICBT fractions. We retrospectively evaluated differences in the ICBT dose contribution to positive pelvic nodes for ten patients who had previously received external beam nodal boost. Our results indicate that the mean dose to the pelvic nodes differed by up to 1.9 Gy between independent ICBT fractions. The second aim is to develop and validate a volumetric method for summing dose of the normal tissues during prescription of nodal boost. The traditional method of dose summation uses the maximum point dose from each modality, which often only represents the worst case scenario. However, the worst case is often an exaggeration when highly conformal therapy methods such as intensity modulated radiation therapy (IMRT) are used. We used deformable image registration algorithms to volumetrically sum dose for the bladder and rectum and created a voxel-by-voxel validation method. The mean error in deformable image registration results of all voxels within the bladder and rectum were 5 and 6 mm, respectively. Finally, the third aim explored the potential use of proton therapy to reduce normal tissue dose. A major physical advantage of protons over photons is that protons stop after delivering dose in the tumor. Although theoretically superior to photons, proton beams are more sensitive to uncertainties caused by interfractional anatomical variations, and must be accounted for during treatment planning to ensure complete target coverage. We have demonstrated a systematic approach to determine population-based anatomical margin requirements for proton therapy. The observed optimal treatment angles for common iliac nodes were 90° (left lateral) and 180° (posterior-anterior [PA]) with additional 0.8 cm and 0.9 cm margins, respectively. For external iliac nodes, lateral and PA beams required additional 0.4 cm and 0.9 cm margins, respectively. Through this project, we have provided radiation oncologists with additional information about potential differences in nodal dose between independent ICBT insertions and volumetric total dose distribution in the bladder and rectum. We have also determined the margins needed for safe delivery of proton therapy when delivering nodal boosts to patients with cervical cancer.

Neutronic design for ESS-Bilbao neutron source

The European Spallation Source-Bilbao (ESS-Bilbao) project plans to build an accelerator facility compliant with the ESS-AB requirements which will be able to drive several experimental stations for research purposes involving intense proton beams with currents up to 75 mA, 50 MeV of final energy, 1.5 ms of pulse length and up to 50 Hz repetition rate. The accelerator will also drive a compact neutron source which will provide useful neutron beams to carry out experiments on moderator optimization, neutron optics devices and general neutron instrumentation as well as preparation work for experiments to be carried out by neutron beam users at the large facilities.

Comparing gold nano-particle enhanced radiotherapy with protons, megavoltage photons and kilovoltage photons: a Monte Carlo simulation

Gold nanoparticles (GNPs) have shown potential to be used as a radiosensitizer for radiation therapy. Despite extensive research activity to study GNP radiosensitization using photon beams, only a few studies have been carried out using proton beams. In this work Monte Carlo simulations were used to assess the dose enhancement of GNPs for proton therapy. The enhancement effect was compared between a clinical proton spectrum, a clinical 6 MV photon spectrum, and a kilovoltage photon source similar to those used in many radiobiology lab settings. We showed that the mechanism by which GNPs can lead to dose enhancements in radiation therapy differs when comparing photon and proton radiation. The GNP dose enhancement using protons can be up to 14 and is independent of proton energy, while the dose enhancement is highly dependent on the photon energy used. For the same amount of energy absorbed in the GNP, interactions with protons, kVp photons and MV photons produce similar doses within several nanometers of the GNP surface, and differences are below 15% for the first 10 nm. However, secondary electrons produced by kilovoltage photons have the longest range in water as compared to protons and MV photons, e.g. they cause a dose enhancement 20 times higher than the one caused by protons 10 μm away from the GNP surface. We conclude that GNPs have the potential to enhance radiation therapy depending on the type of radiation source. Proton therapy can be enhanced significantly only if the GNPs are in close proximity to the biological target.

LHC main dipole magnet circuits: sustaining near-nominal beam energies

Crossing the Franco-Swiss border, the Large Hadron Collider (LHC), designed to collide 7 TeV proton beams, is the world's largest and most powerful particle accelerator the operation of which was originally intended to commence in 2008. Unfortunately, due to an interconnect discontinuity in one of the main dipole circuit's 13 kA superconducting busbars, a catastrophic quench event occurred during initial magnet training, causing significant physical system damage. Furthermore, investigation into the cause found that such discontinuities were not only present in the circuit in question, but throughout the entire LHC. This prevented further magnet training and ultimately resulted in the maximum sustainable beam energy being limited to approximately half that of the design nominal, 3.5-4 TeV, for the first three years of operation (Run 1, 2009-2012) and a major consolidation campaign being scheduled for the first long shutdown (LS 1, 2012-2014). Throughout Run 1, a series of studies attempted to predict the amount of post-installation training quenches still required to qualify each circuit to nominal-energy current levels. With predictions in excess of 80 quenches (each having a recovery time of 8-12+ hours) just to achieve 6.5 TeV and close to 1000 quenches for 7 TeV, it was decided that for Run 2, all systems be at least qualified for 6.5 TeV operation. However, even with all interconnect discontinuities scheduled to be repaired during LS 1, numerous other concerns regarding circuit stability arose. In particular, observations of an erratic behaviour of magnet bypass diodes and the degradation of other potentially weak busbar sections, as well as observations of seemingly random millisecond spikes in beam losses, known as unidentified falling object (UFO) events, which, if persist at 6.5 TeV, may eventually deposit sufficient energy to quench adjacent magnets. In light of the above, the thesis hypothesis states that, even with the observed issues, the LHC main dipole circuits can safely support and sustain near-nominal proton beam energies of at least 6.5 TeV. Research into minimising the risk of magnet training led to the development and implementation of a new qualification method, capable of providing conclusive evidence that all aspects of all circuits, other than the magnets and their internal joints, can safely withstand a quench event at near-nominal current levels, allowing for magnet training to be carried out both systematically and without risk. This method has become known as the Copper Stabiliser Continuity Measurement (CSCM). Results were a success, with all circuits eventually being subject to a full current decay from 6.5 TeV equivalent current levels, with no measurable damage occurring. Research into UFO events led to the development of a numerical model capable of simulating typical UFO events, reproducing entire Run 1 measured event data sets and extrapolating to 6.5 TeV, predicting the likelihood of UFO-induced magnet quenches. Results provided interesting insights into the involved phenomena as well as confirming the possibility of UFO-induced magnet quenches. The model was also capable of predicting that such events, if left unaccounted for, are likely to be commonplace or not, resulting in significant long-term issues for 6.5+ TeV operation. Addressing the thesis hypothesis, the following written works detail the development and results of all CSCM qualification tests and subsequent magnet training as well as the development and simulation results of both 4 TeV and 6.5 TeV UFO event modelling. The thesis concludes, post-LS 1, with the LHC successfully sustaining 6.5 TeV proton beams, but with UFO events, as predicted, resulting in otherwise uninitiated magnet quenches and being at the forefront of system availability issues.

Instrumentation for diagnostics and control of laser-accelerated proton (ion) beams

Suitable instrumentation for laser-accelerated proton (ion) beams is critical for development of integrated, laser-driven ion accelerator systems. Instrumentation aimed at beam diagnostics and control must be applied to the driving laser pulse, the laser-plasma that forms at the target and the emergent proton (ion) bunch in a correlated way to develop these novel accelerators. This report is a brief overview of established diagnostic techniques and new developments based on material presented at the first workshop on 'Instrumentation for Diagnostics and Control of Laser-accelerated Proton (Ion) Beams' in Abingdon, UK. It includes radiochromic film (RCF), image plates (IP), micro-channel plates (MCP), Thomson spectrometers, prompt inline scintillators, time and space-resolved interferometry (TASRI) and nuclear activation schemes. Repetition-rated instrumentation requirements for target metrology are also addressed. (C) 2013 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.