Abstract: BRIKEN is a complex detection system to be installed at the RIB-facility of the RIKEN Nishina Center. It is aimed at the detection of heavy-ion implants, β-particles, γ-rays and β-delayed neutrons. The whole detection setup involves the Advanced Implantation Detection Array (AIDA), two HPGe Clover detectors and a large set of 166 counters of 3He embedded in a high-density polyethylene matrix. This article reports on a novel methodology developed for the conceptual design and optimisation of the 3He-tubes array, aiming at the best possible performance in terms of neutron detection. The algorithm is based on a geometric representation of two selected parameters of merit, namely, average neutron detection efficiency and efficiency flatness, as a function of a reduced number of geometric variables. The response of the detection system itself, for each configuration, is obtained from a systematic MC-simulation implemented realistically in Geant4. This approach has been found to be particularly useful. On the one hand, due to the different types and large number of 3He-tubes involved and, on the other hand, due to the additional constraints introduced by the ancillary detectors for charged particles and gamma-rays. Empowered by the robustness of the algorithm, we have been able to design a versatile detection system, which can be easily re-arranged into a compact mode in order to maximize the neutron detection performance, at the cost of the gamma-ray sensitivity. In summary, we have designed a system which shows, for neutron energies up to 1(5) MeV, a rather flat and high average efficiency of 68.6%(64%) and 75.7%(71%) for the hybrid and compact modes, respectively. The performance of the BRIKEN system has been also quantified realistically by means of MC-simulations made with different neutron energy distributions.

Abstract: It is well established that phenomenological two-zone diffusion models of the galactic halo can very well reproduce cosmic-ray nuclear data and the observed antiproton flux. Here, we consider lepton propagation in such models and compute the expected galactic population of electrons, as well as the diffuse synchrotron emission that results from their interaction with galactic magnetic fields. We find models in agreement not only with cosmic ray data but also with radio surveys at essentially all frequencies. Requiring such a globally consistent description strongly disfavors very large (L greater than or similar to 15 kpc) and, even stronger, small (L less than or similar to 1 kpc) effective diffusive halo sizes. This has profound implications for, e.g., in direct dark matter searches.

Abstract: A two-body fully relativistic three-dimensional scattering framework has been utilised to develop an alternative Compton scattering computational model to those adapted from Ribberfors' work for Monte Carlo modelling of Compton scattering. Using a theoretical foundation that ensures the conservation of energy and momentum in the relativistic impulse approximation, this new model, the Monash University Compton scattering model, develops energy and directional algorithms for both the scattered photon and ejected Compton electron from first principles. The Monash University Compton scattering model was developed to address the limitation of the Compton electron directionality algorithms of other computational models adapted from Ribberfors' work. Here the development of the Monash University Compton scattering model, including its implementation in a Geant4 low energy electromagnetic physics class, G4LowEPComptonModel, is outlined. Assessment of the performance of G4LowEPComptonModel was undertaken in two steps: (1) comparison with respect to the two standard Compton scattering classes of Geant4 version 9.5, G4LivermoreComptonModel and G4PenelopeComptonModel, and (2) experimental comparison with respect to Compton electron kinetic energy spectra obtained from the Compton scattering of 662 key photons off the K-shell of gold. Both studies illustrate that the Monash University Compton scattering model, and in turn G4LowEPComptonModel, is a viable replacement for the majority of computational models that have been adapted from Ribberfors' work. It was also shown that the Monash University Compton scattering model is able to reproduce the Compton scattering triply differential cross-section Compton electron kinetic energy spectra of 662 keV photons K-shell scattering off of gold to within experimental uncertainty.

Abstract: A high-resolution PET system, which incorporates a silicon detector probe into a conventional PET scanner, has been proposed to obtain increased image quality in a limited region of interest. Detailed simulation studies have previously shown that the additional probe information improves the spatial resolution of the reconstructed image and increases lesion detectability, with no cost to other image quality measures. The current study expands on the previous work by using a laboratory prototype of the silicon PET-probe system to examine the resolution improvement in an experimental setting. Two different versions of the probe prototype were assessed, both consisting of a back-to-back pair of 1-mm thick silicon pad detectors, one arranged in 32 x 16 arrays of 1.4mm x 1.4mm pixels and the other in 40 x 26 arrays of 1.0mm x 1.0mm pixels. Each detector was read out by a set of VATAGP7 ASICs and a custom-designed data acquisition board which allowed trigger and data interfacing with the PET scanner, itself consisting of BGO block detectors segmented into 8 x 6 arrays of 6mm x 12mm x 30mm crystals. Limited-angle probe data was acquired from a group of Na-22 point-like sources in order to observe the maximum resolution achievable using the probe system. Data from a Derenzo-like resolution phantom was acquired, then scaled to obtain similar statistical quality as that of previous simulation studies. In this case, images were reconstructed using measurements of the PET ring alone and with the inclusion of the probe data. Images of the Na-22 source demonstrated a resolution of 1.5mm FWHM in the probe data, the PET ring resolution being approximately 6 mm. Profiles taken through the image of the Derenzo-like phantom showed a clear increase in spatial resolution. Improvements in peak-to-valley ratios of 50% and 38%, in the 4.8mm and 4.0mm phantom features respectively, were observed, while previously unresolvable 3.2mm features were brought to light by the addition of the probe. These results support the possibility of improving the image resolution of a clinical PET scanner using the silicon PET-probe.

Abstract: Continuous crystals can potentially obtain better intrinsic detector spatial resolution compared to pixelated crystals, additionally providing depth of interaction (DoI) information from the light distribution. To achieve high performance sophisticated interaction position estimation algorithms are required. There are a number of algorithms in the literature applied to different crystal dimensions and different photodetectors. However, the different crystal properties and photodetector array geometries have an impact on the algorithm performance. In this work we analysed, through Monte Carlo simulations, different combinations of realistic crystals and photodetector parameters to better understand their influence on the interaction position estimation accuracy, with special emphasis on the DoI. We used an interaction position estimation based on an analytical model for the present work. Different photodetector granulation schemes were investigated. The impact of the number of crystal faces readout by photodetectors was studied by simulating scenarios with one and two photodetectors. In addition, crystals with different levels of reflection and aspect ratios (AR) were analysed. Results showed that the impact of photodetector granularity is mainly shown near the edges and specially in the corners of the crystal. The resulting intrinsic spatial resolution near the centre with a 12 x 12 x 10 mm(3) LYSO crystal was 0.7-0.9 mm, while the average spatial resolution calculated on the entire crystal was 0.77 +/- 0.18 mm for all the simulated geometries with one and two photodetectors. Having front and back photodetectors reduced the DoI bias (Euclidean distance between estimated DoI and real DoI) and improved the transversal resolution near the corners. In scenarios with one photodetector, small AR resulted in DoI inaccuracies for absorbed events at the entrance of the crystal. These inaccuracies were slightly reduced either by increasing the AR or reducing the amount of reflected light, and highly mitigated using two photodetectors. Using one photodetector, we obtained a piecewise DoI error model with a DoI resolution of 0.4-0.9 mm for a 1.2 AR crystal, and we observed that including a second photodetector or reducing the amount of reflections reduced the DoI bias but did not significantly improve the DoI resolution. Translating the piecewise DoI error model obtained in this study to image reconstruction we obtained a spatial resolution variability of 0.39 mm using 85% of the FoV, compared to 2.59 mm and 1.87 mm without DoI correction or with a dual layer system, respectively.