Advanced aspects of radiation protection in the use of particle accelerators in the medical field

In this work, the well-known MC code FLUKA was used to simulate the GE PETrace cyclotron (16.5 MeV) installed at “S. Orsola-Malpighi” University Hospital (Bologna, IT) and routinely used in the production of positron emitting radionuclides. Simulations yielded estimates of various quantities of interest, including: the effective dose distribution around the equipment; the effective number of neutron produced per incident proton and their spectral distribution; the activation of the structure of the cyclotron and the vault walls; the activation of the ambient air, in particular the production of 41Ar, the assessment of the saturation yield of radionuclides used in nuclear medicine. The simulations were validated against experimental measurements in terms of physical and transport parameters to be used at the energy range of interest in the medical field. The validated model was also extensively used in several practical applications uncluding the direct cyclotron production of non-standard radionuclides such as 99mTc, the production of medical radionuclides at TRIUMF (Vancouver, CA) TR13 cyclotron (13 MeV), the complete design of the new PET facility of “Sacro Cuore – Don Calabria” Hospital (Negrar, IT), including the ACSI TR19 (19 MeV) cyclotron, the dose field around the energy selection system (degrader) of a proton therapy cyclotron, the design of plug-doors for a new cyclotron facility, in which a 70 MeV cyclotron will be installed, and the partial decommissioning of a PET facility, including the replacement of a Scanditronix MC17 cyclotron with a new TR19 cyclotron.

Improvement of photon transport model by including coupled photon-electron transport and kernel refinement

The first part of this work deals with the inverse problem solution in the X-ray spectroscopy field. An original strategy to solve the inverse problem by using the maximum entropy principle is illustrated. It is built the code UMESTRAT, to apply the described strategy in a semiautomatic way. The application of UMESTRAT is shown with a computational example. The second part of this work deals with the improvement of the X-ray Boltzmann model, by studying two radiative interactions neglected in the current photon models. Firstly it is studied the characteristic line emission due to Compton ionization. It is developed a strategy that allows the evaluation of this contribution for the shells K, L and M of all elements with Z from 11 to 92. It is evaluated the single shell Compton/photoelectric ratio as a function of the primary photon energy. It is derived the energy values at which the Compton interaction becomes the prevailing process to produce ionization for the considered shells. Finally it is introduced a new kernel for the XRF from Compton ionization. In a second place it is characterized the bremsstrahlung radiative contribution due the secondary electrons. The bremsstrahlung radiation is characterized in terms of space, angle and energy, for all elements whit Z=1-92 in the energy range 1–150 keV by using the Monte Carlo code PENELOPE. It is demonstrated that bremsstrahlung radiative contribution can be well approximated with an isotropic point photon source. It is created a data library comprising the energetic distributions of bremsstrahlung. It is developed a new bremsstrahlung kernel which allows the introduction of this contribution in the modified Boltzmann equation. An example of application to the simulation of a synchrotron experiment is shown.

E-ToF-Messungen mit kalorimetrischen Tieftemperatur-Detektoren zur Bestimmung spezischer Energieverluste und zur Massenidentifikation niederenergetischer Schwerionen

Kalorimetrische Tieftemperatur-Detektoren (Calorimetric Low Temperature Detectors, CLTDs) wurden erstmals in Messungen zur Bestimmung des spezifischen Energieverlustes (dE/dx) niederenergetischer Schwerionen beim Durchgang durch Materie eingesetzt. Die Messungen wurden im Energiebereich unterhalb des Bragg-Peaks, mit 0.1 - 1.4 MeV/u 238U-Ionen in Kohlenstoff und Gold sowie mit 0.05 - 1.0 MeV/u 131Xe-Ionen in Kohlenstoff, Nickel und Gold, durchgeführt. Die Kombination der CLTDs mit einem Flugzeitdetektor ermöglichte dabei, kontinuierliche dE/dx-Kurven über größere Energiebereiche hinweg simultan zu bestimmen. Im Vergleich zu herkömmlichen Meßsystemen, die Ionisationsdetektoren zur Energiemessung verwenden, erlaubten die höhere Energieauflösung und -linearität der CLTDs eine Verringerung der Kalibrierungsfehler sowie eine Erweiterung des zugänglichen Energiebereiches der dE/dx-Messungen in Richtung niedriger Energien. Die gewonnen Daten können zur Anpassung theoretischer und semi-empirischer Modelle und somit zu einer Erhöhung der Präzision bei der Vorhersage spezifischer Energieverluste schwerer Ionen beitragen. Neben der experimentellen Bestimmung neuer Daten wurden das alternative Detektionsprinzip der CLTDs, die Vorteile dieser Detektoren bezüglich Energieauflösung und -linearität sowie der modulare Aufbau des CLTD-Arrays aus mehreren Einzeldetektoren genutzt, um diese Art von Messung auf potentielle systematische Unsicherheiten zu untersuchen. Unter anderem wurden hierbei unerwartete Channeling-Effekte beim Durchgang der Ionen durch dünne polykristalline Absorberfolien beobachtet. Die koinzidenten Energie- und Flugzeitmessungen (E-ToF) wurden weiterhin genutzt, um das Auflösungsvermögen des Detektor-Systems bei der direkten in-flight Massenbestimmung langsamer und sehr schwerer Ionen zu bestimmen. Durch die exzellente Energieauflösung der CLTDs konnten hierbei Massenauflösungen von Delta-m(FWHM) = 1.3 - 2.5 u für 0.1 - 0.6 MeV/u 238U-Ionen erreicht werden. In einer E-ToF-Messung mit Ionisationsdetektoren sind solche Werte in diesem Energie- und Massenbereich aufgrund der Limitierung der Energieauflösung durch statistische Schwankungen von Verlustprozessen beim Teilchennachweis nicht erreichbar.

Probing the decay characteristics of the Pygmy Dipole Resonance in the semi-magic nucleus 140 Ce with gamma-gamma coincidence measurements

Im Rahmen dieser Arbeit wurde ein neuartiger Experimentaufbau -- das γ3 Experiment -- zur Messung von photoneninduzierten Kern-Dipolanregungen in stabilen Isotopen konzipiert und an der High Intensity γ-Ray Source (HIγS) an der Duke University installiert.rnDie hohe Energieauflösung und die hohe Nachweiseffizienz des Detektoraufbaus, welcher aus einer Kombination von LaBr Szintillatoren und hochreinen Germanium-Detektoren besteht, erlaubt erstmals die effiziente Messung von γ-γ-Koinzidenzen in Verbindung mit der Methode der Kernresonanzfluoreszenz.rnDiese Methode eröffnet den Zugang zum Zerfallsverhalten der angeregten Dipolzustände als zusätzlicher Observablen, die ein detaillierteres Verständnis der zugrunde liegenden Struktur dieser Anregungen ermöglicht.rnDer Detektoraufbau wurde bereits erfolgreich im Rahmen von zwei Experimentkampagnen in 2012 und 2013 für die Untersuchung von 13 verschiedenen Isotopen verwendet. Im Fokus dieser Arbeit stand die Analyse der Pygmy-Dipolresonanz (PDR) im Kern 140Ce im Energiebereich von 5,2 MeV bis 8,3 MeV basierend auf den mit dem γ3 Experimentaufbau gemessenen Daten. Insbesondere das Zerfallsverhalten der Zustände, die an der PDR beteiligt sind, wurde untersucht. Der Experimentaufbau, die Details der Analyse sowie die Resultate werden in der vorliegenden Arbeit präsentiert. Desweiteren erlaubt ein Vergleich der Ergebnisse mit theoretischen Rechnungen im quasi-particle phonon model (QPM) eine Interpretation des beobachteten Zerfallsverhaltens.

Calibration of a portal imaging device for high-precision dosimetry: a Monte Carlo study

Today electronic portal imaging devices (EPID's) are used primarily to verify patient positioning. They have, however, also the potential as 2D-dosimeters and could be used as such for transit dosimetry or dose reconstruction. It has been proven that such devices, especially liquid filled ionization chambers, have a stable dose response relationship which can be described in terms of the physical properties of the EPID and the pulsed linac radiation. For absolute dosimetry however, an accurate method of calibration to an absolute dose is needed. In this work, we concentrate on calibration against dose in a homogeneous water phantom. Using a Monte Carlo model of the detector we calculated dose spread kernels in units of absolute dose per incident energy fluence and compared them to calculated dose spread kernels in water at different depths. The energy of the incident pencil beams varied between 0.5 and 18 MeV. At the depth of dose maximum in water for a 6 MV beam (1.5 cm) and for a 18 MV beam (3.0 cm) we observed large absolute differences between water and detector dose above an incident energy of 4 MeV but only small relative differences in the most frequent energy range of the beam energy spectra. It is shown that for a 6 MV beam the absolute reference dose measured at 1.5 cm water depth differs from the absolute detector dose by 3.8%. At depth 1.2 cm in water, however, the relative dose differences are almost constant between 2 and 6 MeV. The effects of changes in the energy spectrum of the beam on the dose responses in water and in the detector are also investigated. We show that differences larger than 2% can occur for different beam qualities of the incident photon beam behind water slabs of different thicknesses. It is therefore concluded that for high-precision dosimetry such effects have to be taken into account. Nevertheless, the precise information about the dose response of the detector provided in this Monte Carlo study forms the basis of extracting directly the basic radiometric quantities photon fluence and photon energy fluence from the detector's signal using a deconvolution algorithm. The results are therefore promising for future application in absolute transit dosimetry and absolute dose reconstruction.

HELIOSPHERIC ENERGETIC NEUTRAL HYDROGEN MEASURED WITH ASPERA-3 AND ASPERA-4

We re-analyze the signal of non-planetary energetic neutral atoms (ENAs) in the 0.4-5.0 keV range measured with the Neutral Particle Detector (NPD) of the ASPERA-3 and ASPERA-4 experiments on board the Mars and Venus Express satellites. Due to improved knowledge of sensor characteristics and exclusion of data sets affected by instrument effects, the typical intensity of the ENA signal obtained by ASPERA-3 is an order of magnitude lower than in earlier reports. The ENA intensities measured with ASPERA-3 and ASPERA-4 now agree with each other. In the present analysis, we also correct the ENA signal for Compton-Getting and for ionization loss processes under the assumption of a heliospheric origin. We find spectral shapes and intensities consistent with those measured by the Interstellar Boundary Explorer (IBEX). The principal advantage of ASPERA with respect to the IBEX sensors is the two times better spectral resolution. In this study, we discuss the physical significance of the spectral shapes and their potential variation across the sky. At present, these observations are the only independent test of the heliospheric ENA signal measured with IBEX in this energy range. The ASPERA measurements also allow us to check for a temporal variation of the heliospheric signal as they were obtained between 2003 and 2007, whereas IBEX has been operational since the end of 2008.

Reflection of solar wind hydrogen from the lunar surface

The solar wind continuously flows out from the Sun and directly interacts with the surfaces of dust and airless planetary bodies throughout the solar system. A significant fraction of solar wind ions reflect from an object's surface as energetic neutral atoms (ENAs). ENA emission from the Moon was first observed during commissioning of the Interstellar Boundary Explorer (IBEX) mission on 3 December 2008. We present the analysis of 10 additional IBEX observations of the Moon while it was illuminated by the solar wind. For the viewing geometry and energy range (> 250 eV) of the IBEX-Hi ENA imager, we find that the spectral shape of the ENA emission from the Moon is well-represented by a linearly decreasing flux with increasing energy. The fraction of the incident solar wind ions reflected as ENAs, which is the ENA albedo and defined quantitatively as the ENA reflection coefficient RN, depends on the incident solar wind speed, ranging from ~0.2 for slow solar wind to ~0.08 for fast solar wind. The average energy per incident solar wind ion that is reflected to space is 30 eV for slow solar wind and 45 eV for fast solar wind. Once ionized, these ENAs can become pickup ions in the solar wind with a unique spectral signature that reaches 3vSW. These results apply beyond the solar system; the reflection process heats plasmas that have significant bulk flow relative to interstellar dust and cools plasmas having no net bulk flow relative to the dust.

Solar wind reflection from the lunar surface: The view from far and near

The Moon appears bright in the sky as a source of energetic neutral atoms (ENAs). These ENAs have recently been imaged over a broad energy range both from near the lunar surface, by India's Chandrayaan-1 mission (CH-1), and from a much more distant Earth orbit by NASA's Interstellar Boundary Explorer (IBEX) satellite. Both sets of observations have indicated that a relatively large fraction of the solar wind is reflected from the Moon as energetic neutral hydrogen. CH-1's angular resolution over different viewing angles of the lunar surface has enabled measurement of the emission as a function of angle. IBEX in contrast views not just a swath but a whole quadrant of the Moon as effectively a single pixel, as it subtends even at the closest approach no more than a few degrees on the sky. Here we use the scattering function measured by CH-1 to model global lunar ENA emission and combine these with IBEX observations. The deduced global reflection is modestly larger (by a factor of 1.25) when the angular scattering function is included. This provides a slightly updated IBEX estimate of AH=0.11±0.06 for the global neutralized albedo, which is ˜25% larger than the previous values of 0.09±0.05, based on an assumed uniform scattering distribution.

Thermal neutron capture effects in radioactive and stable nuclide systems

Neutron capture effects in meteorites and lunar surface samples have been successfully used in the past to study exposure histories and shielding conditions. In recent years, however, it turned out that neutron capture effects produce a nuisance for some of the short-lived radionuclide systems. The most prominent example is the 182Hf-182W system in iron meteorites, for which neutron capture effects lower the 182W/184W ratio, thereby producing too old apparent ages. Here, we present a thorough study of neutron capture effects in iron meteorites, ordinary chondrites, and carbonaceous chondrites, whereas the focus is on iron meteorites. We study in detail the effects responsible for neutron production, neutron transport, and neutron slowing down and find that neutron capture in all studied meteorite types is not, as usually expected, exclusively via thermal neutrons. In contrast, most of the neutron capture in iron meteorites is in the epithermal energy range and there is a significant contribution from epithermal neutron capture even in stony meteorites. Using sophisticated particle spectra and evaluated cross section data files for neutron capture reactions we calculate the neutron capture effects for Sm, Gd, Cd, Pd, Pt, and Os isotopes, which all can serve as neutron-dose proxies, either in stony or in iron meteorites. In addition, we model neutron capture effects in W and Ag isotopes. For W isotopes, the GCR-induced shifts perfectly correlate with Os and Pt isotope shifts, which therefore can be used as neutron-dose proxies and permit a reliable correction. We also found that GCR-induced effects for the 107Pd-107Ag system can be significant and need to be corrected, a result that is in contrast to earlier studies.

First direct observation of sputtered lunar oxygen

We present the first direct measurement of neutral oxygen in the lunar exosphere, detected by the Chandrayaan-1 Energetic Neutral Analyzer (CENA). With the lunar surface consisting of about 60% of oxygen in number, the neutral oxygen detected in CENA's energy range (11 eV−3.3 keV) is attributed to have originated from the lunar surface, where it was released through solar wind ion sputtering. Fitting of CENA's mass spectra with calibration spectra from ground and in-flight data resulted in the detection of a robust oxygen signal, with a flux of 0.2 to 0.4 times the flux of backscattered hydrogen, depending on the solar wind helium content and particle velocity. For the two solar wind types observed, we derive subsolar surface oxygen atom densities of N0= (1.1 ± 0.3) · 107m−3 and (1.4 ± 0.4) · 107m−3, respectively, which agree well with earlier model predictions and measured upper limits. From these surface densities, we derive column densities of NC= (1.5 ± 0.5) · 1013 m−2and (1.6 ± 0.5) · 1013 m−2. In addition, we identified for the first time a helium component. This helium is attributed to backscattering of solar wind helium (alpha particles) from the lunar surface as neutral energetic helium atoms, which has also been observed for the first time. This identification is supported by the characteristic energy of the measured helium atoms, which is roughly 4 times the energy of reflected solar wind hydrogen, and the correlation with solar wind helium content.

Measurement of the TeV atmospheric muon charge ratio with the complete OPERA data set

The OPERA detector, designed to search for νμ → ντ oscillations in the CNGS beam, is located in the underground Gran Sasso laboratory, a privileged location to study TeV-scale cosmic rays. For the analysis here presented, the detector was used to measure the atmospheric muon charge ratio in the TeV region. OPERA collected chargeseparated cosmic ray data between 2008 and 2012. More than 3 million atmospheric muon events were detected and reconstructed, among which about 110000 multiple muon bundles. The charge ratio Rμ ≡ Nμ+/Nμ− was measured separately for single and for multiple muon events. The analysis exploited the inversion of the magnet polarity which was performed on purpose during the 2012 Run. The combination of the two data sets with opposite magnet polarities allowedminimizing systematic uncertainties and reaching an accurate determination of the muon charge ratio. Data were fitted to obtain relevant parameters on the composition of primary cosmic rays and the associated kaon production in the forward fragmentation region. In the surface energy range 1–20 TeV investigated by OPERA, Rμ is well described by a parametric model including only pion and kaon contributions to themuon flux, showing no significant contribution of the prompt component. The energy independence supports the validity of Feynman scaling in the fragmentation region up to 200 TeV/nucleon primary energy.

Exclusion of leptophilic dark matter models using XENON100 electronic recoil data

We have searched for periodic variations of the electronic recoil event rate in the (2-6) keV energy range recorded between February 2011 and March 2012 with the XENON100 detector, adding up to 224.6 live days in total. Following a detailed study to establish the stability of the detector and its background contributions during this run, we performed an un-binned profile likelihood analysis to identify any periodicity up to 500 days. We find a global significance of less than 1 sigma for all periods suggesting no statistically significant modulation in the data. While the local significance for an annual modulation is 2.8 sigma, the analysis of a multiple-scatter control sample and the phase of the modulation disfavor a dark matter interpretation. The DAMA/LIBRA annual modulation interpreted as a dark matter signature with axial-vector coupling of WIMPs to electrons is excluded at 4.8 sigma.

Search for Event Rate Modulation in XENON100 Electronic Recoil Data

We have searched for periodic variations of the electronic recoil event rate in the (2-6) keV energy range recorded between February 2011 and March 2012 with the XENON100 detector, adding up to 224.6 live days in total. Following a detailed study to establish the stability of the detector and its background contributions during this run, we performed an un-binned profile likelihood analysis to identify any periodicity up to 500 days. We find a global significance of less than 1 sigma for all periods suggesting no statistically significant modulation in the data. While the local significance for an annual modulation is 2.8 sigma, the analysis of a multiple-scatter control sample and the phase of the modulation disfavor a dark matter interpretation. The DAMA/LIBRA annual modulation interpreted as a dark matter signature with axial-vector coupling of WIMPs to electrons is excluded at 4.8 sigma.

First results obtained using the CENBG nanobeam line: performances and applications

A high resolution focused beam line has been recently installed on the AIFIRA (“Applications Interdisciplinaires des Faisceaux d’Ions en Région Aquitaine”) facility at CENBG. This nanobeam line, based on a doublet–triplet configuration of Oxford Microbeam Ltd. OM-50™ quadrupoles, offers the opportunity to focus protons, deuterons and alpha particles in the MeV energy range to a sub-micrometer beam spot. The beam optics design has been studied in detail and optimized using detailed ray-tracing simulations and the full mechanical design of the beam line was reported in the Debrecen ICNMTA conference in 2008. During the last two years, the lenses have been carefully aligned and the target chamber has been fully equipped with particle and X-ray detectors, microscopes and precise positioning stages. The beam line is now operational and has been used for its firstapplications to ion beam analysis. Interestingly, this set-up turned out to be a very versatile tool for a wide range of applications. Indeed, even if it was not intended during the design phase, the ion optics configuration offers the opportunity to work either with a high current microbeam (using the triplet only) or with a lower current beam presenting a sub-micrometer resolution (using the doublet–triplet configuration). The performances of the CENBGnanobeam line are presented for both configurations. Quantitative data concerning the beam lateral resolutions at different beam currents are provided. Finally, the firstresults obtained for different types of application are shown, including nuclear reaction analysis at the micrometer scale and the firstresults on biological samples

Self-consistent, nonlocal electron heat flux at arbitrary ion charge number

A single, nonlocal expression for the electron heat flux, which closely reproduces known results at high and low ion charge number 2, and “exact” results for the local limit at all 2, is derived by solving the kinetic equation in a narrow, tail-energy range. The solution involves asymptotic expansions of Bessel functions of large argument, and (Z-dependent)order above or below it, corresponding to the possible parabolic or hyperbolic character of the kinetic equation; velocity space diffusion in self-scattering is treated similarly to isotropic thermalization of tail energies in large Z analyses. The scale length H characterizing nonlocal effects varies with Z, suggesting an equal dependence of any ad hoc flux limiter. The model is valid for all H above the mean-free path for thermal electrons.