Particle Accelerators and Detectors for medical Diagnostics and Therapy

This Habilitationsschrift (Habilitation thesis) is focused on my research activities on medical applications of particle physics and was written in 2013 to obtain the Venia Docendi (Habilitation) in experimental physics at the University of Bern. It is based on selected publications, which represented at that time my major scientific contributions as an experimental physicist to the field of particle accelerators and detectors applied to medical diagnostics and therapy. The thesis is structured in two parts. In Part I, Chapter 1 presents an introduction to accelerators and detectors applied to medicine, with particular focus on cancer hadrontherapy and on the production of radioactive isotopes. In Chapter 2, my publications on medical particle accelerators are introduced and put into their perspective. In particular, high frequency linear accelerators for hadrontherapy are discussed together with the new Bern cyclotron laboratory. Chapter 3 is dedicated to particle detectors with particular emphasis on three instruments I contributed to propose and develop: segmented ionization chambers for hadrontherapy, a proton radiography apparatus with nuclear emulsion films, and a beam monitor detector for ion beams based on doped silica fibres. Selected research and review papers are contained in Part II. For copyright reasons, they are only listed and not reprinted in this on-line version. They are available on the websites of the journals.

Dark matter sensitivity of multi-ton liquid xenon detectors

We study the sensitivity of multi ton-scale time projection chambers using a liquid xenon target, e.g., the proposed DARWIN instrument, to spin-independent and spin-dependent WIMP-nucleon scattering interactions. Taking into account realistic backgrounds from the detector itself as well as from neutrinos, we examine the impact of exposure, energy threshold, background rejection efficiency and energy resolution on the dark matter sensitivity. With an exposure of 200 t x y and assuming detector parameters which have been already demonstrated experimentally, spin-independent cross sections as low as 2.5×10−49 cm2 can be probed for WIMP masses around 40 GeV/c2. Additional improvements in terms of background rejection and exposure will further increase the sensitivity, while the ultimate WIMP science reach will be limited by neutrinos scattering coherently off the xenon nuclei.

Measurement of isotopic uranium in water for compliance monitoring by liquid scintillation counting with alpha/beta discrimination

A simple and inexpensive method is described for analysis of uranium (U) activity and mass in water by liquid scintillation counting using $\alpha$/$\beta$ discrimination. This method appears to offer a solution to the need for an inexpensive protocol for monitoring U activity and mass simultaneously and an alternative to the potential inaccuracy involved when depending on the mass-to-activity conversion factor or activity screen.^ U is extracted virtually quantitatively into 20 ml extractive scintillator from a 1-$\ell$ aliquot of water acidified to less than pH 2. After phase separation, the sample is counted for a 20-minute screening count with a minimum detection level of 0.27 pCi $\ell\sp{-1}$. $\alpha$-particle emissions from the extracted U are counted with close to 100% efficiency with a Beckman LS6000 LL liquid scintillation counter equipped with pulse-shape discrimination electronics. Samples with activities higher than 10 pCi $\ell\sp-1$ are recounted for 500-1000 minutes for isotopic analysis. Isotopic analysis uses events that are automatically stored in spectral files and transferred to a computer during assay. The data can be transferred to a commercially available spreadsheet and retrieved for examination or data manipulation. Values for three readily observable spectral features can be rapidly identified by data examination and substituted into a simple formula to obtain $\sp{234}$U/$\sp{238}$U ratio for most samples. U mass is calculated by substituting the isotopic ratio value into a simple equation.^ The utility of this method for the proposed compliance monitoring of U in public drinking water supplies was field tested with a survey of drinking water from Texas supplies that had previously been known to contain elevated levels of gross $\alpha$ activity. U concentrations in 32 samples from 27 drinking water supplies ranged from 0.26 to 65.5 pCi $\ell\sp{-1}$, with seven samples exceeding the proposed Maximum Contaminant Level of 20 $\mu$g $\ell\sp{-1}$. Four exceeded the proposed activity screening level of 30 pCi $\ell\sp{-1}$. Isotopic ratios ranged from 0.87 to 41.8, while one sample contained $\sp{234}$U activity of 34.6 pCi $\ell\sp{-1}$ in the complete absence of its parent, $\sp{238}$U. U mass in the samples with elevated activity ranged from 0.0 to 103 $\mu$g $\ell\sp{-1}$. A limited test of screening surface and groundwaters for contamination by U from waste sites and natural processes was also successful. ^

Development and Implementation of the Use of Optically Stimulated Luminescent Detectors in the Radiological Physics Center Anthropomorphic Quality Assurance Phantoms

The Radiological Physics Center (RPC) uses both on-site and remote reviews to credential institutions for participation in clinical trials. Anthropomorphic quality assurance (QA) phantoms are one tool the RPC uses to remotely audit institutions, which include thermoluminescent dosimeters (TLDs) and radiochromic film. The RPC desires to switch from TLD as the absolute dosimeter in the phantoms, to optically stimulated luminescent dosimeters (OSLDs), but a problem lies in the angular dependence exhibited by the OSLD. The purpose of this study was to characterize the angular dependence of OSLD and establish a correction factor if necessary, to provide accurate dosimetric measurements as a replacement for TLD in the QA phantoms. A 10 cm diameter high-impact polystyrene spherical phantom was designed and constructed to hold an OSLD to study the angular response of the dosimeter under the simplest of circumstances for both coplanar and non-coplanar treatment deliveries. OSLD were irradiated in the spherical phantom, and the responses of the dosimeter from edge-on angles were normalized to the response when irradiated with the beam incident normally on the surface of the dosimeter. The average normalized response was used to establish an angular correction factor for 6 MV and 18 coplanar treatments, and for 6 MV non-coplanar treatments specific to CyberKnife. The RPC pelvic phantom dosimetry insert was modified to hold OSLD, in addition to the TLD, adjacent to the planes of film. Treatment plans of increasing angular beam delivery were developed, three in Pinnacle v9.0 (4-field box, IMRT, and VMAT) and one in Accuray’s MultiPlan v3.5.3 (CyberKnife). The plans were delivered to the pelvic phantom containing both TLD and OSLD in the target volume. The pelvic phantom was also sent to two institutions to be irradiated as trials, one delivering IMRT, and the other a CyberKnife treatment. For the IMRT deliveries and the two institution trials, the phantom also included film in the sagittal and coronal planes. The doses measured from the TLD and OSLD were calculated for each irradiation, and the angular correction factors established from the spherical phantom irradiations were applied to the OSLD dose. The ratio of the TLD dose to the angular corrected OSLD dose was calculated for each irradiation. The corrected OSLD dose was found to be within 1% of the TLD measured dose for all irradiations, with the exception of the in-house CyberKnife deliveries. The films were normalized to both TLD measured dose and the corrected OSLD dose. Dose profiles were obtained and gamma analysis was performed using a 7%/4 mm criteria, to compare the ability of the OSLD, when corrected for the angular dependence, to provide equivalent results to TLD. The results of this study indicate that the OSLD can effectively be used as a replacement for TLD in the RPC’s anthropomorphic QA phantoms for coplanar treatment deliveries when a correction is applied for the dosimeter’s angular dependence.

Characterization of the dose distribution in the halo region of a clinical proton pencil beam using emulsion film detectors

Proton therapy is a high precision technique in cancer radiation therapy which allows irradiating the tumor with minimal damage to the surrounding healthy tissues. Pencil beam scanning is the most advanced dose distribution technique and it is based on a variable energy beam of a few millimeters FWHM which is moved to cover the target volume. Due to spurious effects of the accelerator, of dose distribution system and to the unavoidable scattering inside the patient's body, the pencil beam is surrounded by a halo that produces a peripheral dose. To assess this issue, nuclear emulsion films interleaved with tissue equivalent material were used for the first time to characterize the beam in the halo region and to experimentally evaluate the corresponding dose. The high-precision tracking performance of the emulsion films allowed studying the angular distribution of the protons in the halo. Measurements with this technique were performed on the clinical beam of the Gantry1 at the Paul Scherrer Institute. Proton tracks were identified in the emulsion films and the track density was studied at several depths. The corresponding dose was assessed by Monte Carlo simulations and the dose profile was obtained as a function of the distance from the center of the beam spot.

Emulsion detectors for the antihydrogen detection in AEgIS

The AEgIS experiment at CERN aims to perform the first direct measurement of gravitational interaction between matter and antimatter by measuring the deviation of a cold antihydrogen beam in the Earth gravitational field. The design of the experiment has been recently updated to include emulsion films as position sensitive detector. The submicrometric position accuracy of emulsions leads indeed to a significant improvement of the experimental sensitivity. We present results of preliminary tests and discuss perspectives for the final measurement.

Scintillation efficiency of liquid argon in low energy neutron-argon scattering

Experiments searching for weak interacting massive particles with noble gases such as liquid argon require very low detection thresholds for nuclear recoils. A determination of the scintillation efficiency is crucial to quantify the response of the detector at low energy. We report the results obtained with a small liquid argon cell using a monoenergetic neutron beam produced by a deuterium-deuterium fusion source. The light yield relative to electrons was measured for six argon recoil energies between 11 and 120 keV at zero electric drift field.

Room temperature absorption in laterally biased Quantum Infrared Detectors fabricated by MBE regrowth

In this paper, we show room temperature operation of a quantum well infrared photodetector (QWIP) using lateral conduction through ohmic contacts deposited at both sides of two n-doped quantum wells. To reduce the dark current due to direct conduction in the wells, we apply an electric field between the quantum wells and two pinch-off Schottky gates, in a fashion similar to a field effect device. Since the normal incidence absorption is strongly reduced in intersubband transitions in quantum wells, we first analyze the response of a detector based on quantum dots (QD). This QD device shows photocurrent signal up to 150 K when it is processed in conventional vertical detector. However, it is possible to observe room temperature signal when it is processed in a lateral structure. Finally, the room temperature photoresponse of the QWIP is demonstrated, and compared with theory. An excellent agreement between the estimated and measured characteristics of the device is found

Projector Model for Efficient List-Mode Reconstruction in PET Scanners with Parallel Planar Detectors

We have developed a new projector model specifically tailored for fast list-mode tomographic reconstructions in Positron emission tomography (PET) scanners with parallel planar detectors. The model provides an accurate estimation of the probability distribution of coincidence events defined by pairs of scintillating crystals. This distribution is parameterized with 2D elliptical Gaussian functions defined in planes perpendicular to the main axis of the tube of response (TOR). The parameters of these Gaussian functions have been obtained by fitting Monte Carlo simulations that include positron range, acolinearity of gamma rays, as well as detector attenuation and scatter effects. The proposed model has been applied efficiently to list-mode reconstruction algorithms. Evaluation with Monte Carlo simulations over a rotating high resolution PET scanner indicates that this model allows to obtain better recovery to noise ratio in OSEM (ordered-subsets, expectation-maximization) reconstruction, if compared to list-mode reconstruction with symmetric circular Gaussian TOR model, and histogram-based OSEM with precalculated system matrix using Monte Carlo simulated models and symmetries.

Tropospheric scintillation with concurrent rain attenuation at 50 GHz in Madrid

Tropospheric scintillation can become a significant impairment in satellite communication systems, especially in those with low fade-margin. Moreover, fast amplitude fluctuations due to scintillation are even larger when rain is present on the propagation path. Few studies of scintillation during rain have been reported and the statistical characterization is still not totally clear. This paper presents experimental results on the relationship between scintillation and rain attenuation obtained from slant-path attenuation measurements at 50 GHz. The study is focused on the probability density function (PDF) of various scintillation parameters. It is shown that scintillation intensity, measured as the standard deviation of the amplitude fluctuations, increases with rain attenuation; in the range 1-10 dB this relationship can be expressed by power-law or linear equations. The PDFs of scintillation intensity conditioned to a given rain attenuation level are lognormal, while the overall long-term PDF is well fltted by a generalized extreme valué (GEV) distribution. The short-term PDFs of amplitude conditioned to a given intensity are normal, although skewness effects are observed for the strongest intensities. A procedure is given to derive numerically the overall PDF of scintillation amplitude using a combination of conditional PDFs and local statistics of rain attenuation.

Failure Detectors in Homonymous Distributed Systems (with an Application to Consensus).

This paper is on homonymous distributed systems where processes are prone to crash failures and have no initial knowledge of the system membership (?homonymous? means that several processes may have the same identi?er). New classes of failure detectors suited to these systems are ?rst de?ned. Among them, the classes H? and H? are introduced that are the homonymous counterparts of the classes ? and ?, respectively. (Recall that the pair h?,?i de?nes the weakest failure detector to solve consensus.) Then, the paper shows how H? and H? can be implemented in homonymous systems without membership knowledge (under different synchrony requirements). Finally, two algorithms are presented that use these failure detectors to solve consensus in homonymous asynchronous systems where there is no initial knowledge ofthe membership. One algorithm solves consensus with hH?, H?i, while the other uses only H?, but needs a majority of correct processes. Observe that the systems with unique identi?ers and anonymous systems are extreme cases of homonymous systems from which follows that all these results also apply to these systems. Interestingly, the new failure detector class H? can be implemented with partial synchrony, while the analogous class A? de?ned for anonymous systems can not be implemented (even in synchronous systems). Hence, the paper provides us with the ?rst proof showing that consensus can be solved in anonymous systems with only partial synchrony (and a majority of correct processes).

Design of a training tool for improving the use of hand-held detectors in humanitarian demining

Purpose – The purpose of this paper is to introduce the design of a training tool intended to improve deminers' technique during close-in detection tasks. Design/methodology/approach – Following an introduction that highlights the impact of mines and improvised explosive devices (IEDs), and the importance of training for enhancing the safety and the efficiency of the deminers, this paper considers the utilization of a sensory tracking system to study the skill of the hand-held detector expert operators. With the compiled information, some critical performance variables can be extracted, assessed, and quantified, so that they can be used afterwards as reference values for the training task. In a second stage, the sensory tracking system is used for analysing the trainee skills. The experimentation phase aims to test the effectiveness of the elements that compose the sensory system to track the hand-held detector during the training sessions. Findings – The proposed training tool will be able to evaluate the deminers' efficiency during the scanning tasks and will provide important information for improving their competences. Originality/value – This paper highlights the need of introducing emerging technologies for enhancing the current training techniques for deminers and proposes a sensory tracking system that can be successfully utilised for evaluating trainees' performance with hand-held detectors.

A function using cubic splines for the analysis of alpha--particles spectra from silicon detectors

A function based on the characteristics of the alpha-particle lines obtained with silicon semiconductor detectors and modi"ed by using cubic splines is proposed to parametrize the shape of the peaks. A reduction in the number of parameters initially considered in other proposals was carried out in order to improve the stability of the optimization process. It was imposed by the boundary conditions for the cubic splines term. This function was then able to describe peaks with highly anomalous shapes with respect to those expected from this type of detector. Some criteria were implemented to correctly determine the area of the peaks and their errors. Comparisons with other well-established functions revealed excellent agreement in the "nal values obtained from both "ts. Detailed studies on reliability of the "tting results were carried out and the application of the function is proposed. Although the aim was to correct anomalies in peak shapes, the peaks showing the expected shapes were also well "tted. Accordingly, the validity of the proposal is quite general in the analysis of alpha-particle spectrometry with semiconductor detectors.

Real-time Surveillance Application by Multiple Detectors and Compressive Trackers

A real-time surveillance system for IP network cameras is presented. Motion, part-body, and whole-body detectors are efficiently combined to generate robust and fast detections, which feed multiple compressive trackers. The generated trajectories are then improved using a reidentification strategy for long term operation.

Detection of Biological CO2 and 1,3-Pentadiene Using Non-refrigerated Low-Cost MWIR Detectors

The early detection of spoiling metabolic products in contaminated food is a very important tool to control quality. Some volatile compounds produce unpleasant odours at very low concentrations, making their early detection very challenging. This is the case of 1,3-pentadiene produced by microorganisms through decarboxylation of the preservative sorbate. In this work, we have developed a methodology to use the data produced by a low-cost, compact MWIR (Mid-Wave IR) spectrometry device without moving parts, which is based on a linear array of 128 elements of VPD PbSe coupled to a linear variable filter (LVF) working in the spectral range between 3 and 4.6 ?m. This device is able to analyze food headspace gases through dedicated sample presentation setup. This methodology enables the detection of CO2 and the volatile compound 1,3-pentadiene, as compared to synthetic patrons. Data analysis is based on an automated multidimensional dynamic processing of the MWIR spectra. Principal component and discriminant analysis allow segregating between four yeast strains including producers and no producers. The segregation power is accounted as a measure of the discrimination quality.