Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2023). Impact of QCD uncertainties on antiproton spectra from dark-matter annihilation. J. Cosmol. Astropart. Phys., 04(4), 068–15pp.
Abstract: Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron de-cays, before particles that are stable on astrophysical time scales are produced. Antiprotons produced in this way may leave footprints in experiments such as AMS-02. Several groups have reported an excess of events in the antiproton flux in the rigidity range of 10-20 GV. However, the theoretical modeling of baryon production is not straightforward and relies in part on phenomenological models in Monte Carlo event generators. In this work, we assess the impact of QCD uncertainties on the spectra of antiprotons from dark-matter annihila-tion. As a proof-of-principle, we show that for a two-parameter model that depends only on the thermally-averaged annihilation cross section ((o -v)) and the dark-matter mass (Mx), QCD uncertainties can affect the best-fit mass by up to ti 14% (with large uncertainties for large DM masses), depending on the choice of Mx and the annihilation channel (bb over bar or W+W-), and (o -v) by up to ti 10%. For comparison, changes to the underlying diffusion parameters are found to be within 1%-5%, and the results are also quite resilient to the choice of cosmic-ray propagation model. These findings indicate that QCD uncertainties need to be included in future DM analyses. To facilitate full-fledged analyses, we provide the spectra in tabulated form including QCD uncertainties and code snippets to perform mass interpolations and quick DM fits. The code can be found in this GitHub [1] repository.
|
DUNE Collaboration(Abud, A. A. et al), Amedo, P., Antonova, M., Barenboim, G., Cervera-Villanueva, A., De Romeri, V., et al. (2023). Highly-parallelized simulation of a pixelated LArTPC on a GPU. J. Instrum., 18(4), P04034–35pp.
Abstract: The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we present the first implementation of a full microphysical simulator of a liquid argon time projection chamber (LArTPC) equipped with light readout and pixelated charge readout, developed for the DUNE Near Detector. The software is implemented with an end-to-end set of GPU-optimized algorithms. The algorithms have been written in Python and translated into CUDA kernels using Numba, a just-in-time compiler for a subset of Python and NumPy instructions. The GPU implementation achieves a speed up of four orders of magnitude compared with the equivalent CPU version. The simulation of the current induced on 103 pixels takes around 1 ms on the GPU, compared with approximately 10 s on the CPU. The results of the simulation are compared against data from a pixel-readout LArTPC prototype.
|
Borja-Lloret, M., Barrientos, L., Bernabeu, J., Lacasta, C., Muñoz, E., Ros, A., et al. (2023). Influence of the background in Compton camera images for proton therapy treatment monitoring. Phys. Med. Biol., 68(14), 144001–16pp.
Abstract: Objective. Background events are one of the most relevant contributions to image degradation in Compton camera imaging for hadron therapy treatment monitoring. A study of the background and its contribution to image degradation is important to define future strategies to reduce the background in the system. Approach. In this simulation study, the percentage of different kinds of events and their contribution to the reconstructed image in a two-layer Compton camera have been evaluated. To this end, GATE v8.2 simulations of a proton beam impinging on a PMMA phantom have been carried out, for different proton beam energies and at different beam intensities. Main results. For a simulated Compton camera made of Lanthanum (III) Bromide monolithic crystals, coincidences caused by neutrons arriving from the phantom are the most common type of background produced by secondary radiations in the Compton camera, causing between 13% and 33% of the detected coincidences, depending on the beam energy. Results also show that random coincidences are a significant cause of image degradation at high beam intensities, and their influence in the reconstructed images is studied for values of the time coincidence windows from 500 ps to 100 ns. Significance. Results indicate the timing capabilities required to retrieve the fall-off position with good precision. Still, the noise observed in the image when no randoms are considered make us consider further background rejection methods.
|
Caron, S., Eckner, C., Hendriks, L., Johannesson, G., Ruiz de Austri, R., & Zaharijas, G. (2023). Mind the gap: the discrepancy between simulation and reality drives interpretations of the Galactic Center Excess. J. Cosmol. Astropart. Phys., 06(6), 013–56pp.
Abstract: The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model (-yEM) used in Galactic center analyses affects the interpretation of the GCE's nature. To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of-yEMs with increasing parametric freedom. We calculate the fractional contribution (fsrc) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the ryEM. For the simplest ryEM, we obtain fsrc = 0.10 f 0.07, while the most complex model yields fsrc = 0.79 f 0.24. In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed ryEM. The quoted results for fsrc do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated ryEM iterations. We quantify the reality gap between our ryEMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed ryEMs have gaps to reality. Our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked “out of domain” uncertainties in future interpretations.
|
Balaudo, A., Calore, F., De Romeri, V., & Donato, F. (2024). NAJADS: a self-contained framework for the direct determination of astrophysical J-factors. J. Cosmol. Astropart. Phys., 02(2), 001–33pp.
Abstract: Cosmological simulations play a pivotal role in understanding the properties of the dark matter (DM) distribution in both galactic and galaxy -cluster environments. The characterization of DM structures is crucial for informing indirect DM searches, aiming at the detection of the annihilation (or decay) products of DM particles. A fundamental quantity in these analyses is the astrophysical J -factor. In the DM phenomenology community, J -factors are typically computed through the semi -analytical modelling of the DM mass distribution, which is affected by large uncertainties. With the scope of addressing and possibly reducing these uncertainties, we present NAJADS, a self-contained framework to derive the DM J -factor directly from the raw simulations data. We show how this framework can be used to compute all -sky maps of the J -factor, automatically accounting for the complex 3D structure of the simulated halos and for the boosting of the signal due to the density fluctuations along the line of sight. After validating our code, we present a proof -of -concept application of NAJADS to a realistic halo from the IllustrisTNG suite, and exploit it to make a thorough comparison between our numerical approach and traditional semi -analytical methods. JCAP02(2024)001
|
Arina, C., Di Mauro, M., Fornengo, N., Heisig, J., Jueid, A., & Ruiz de Austri, R. (2024). CosmiXs: cosmic messenger spectra for indirect dark matter searches. J. Cosmol. Astropart. Phys., 03(3), 035–41pp.
Abstract: The energy spectra of particles produced from dark matter (DM) annihilation or decay are one of the fundamental ingredients to calculate the predicted fluxes of cosmic rays and radiation searched for in indirect DM detection. We revisit the calculation of the source spectra for annihilating and decaying DM using the VINCIA shower algorithm in PYTHIA to include QED and QCD final state radiation and diagrams for the EW corrections with massive bosons, not present in the default PYTHIA shower model. We take into account the spin information of the particles during the entire EW shower and the off -shell contributions from massive gauge bosons. Furthermore, we perform a dedicated tuning of the VINCIA and PYTHIA parameters to LEP data on the production of pions, photons, and hyperons at the Z resonance and discuss the underlying uncertainties. To enable the use of our results in DM studies, we provide the tabulated source spectra for the most relevant cosmic messenger particles, namely antiprotons, positrons, gamma rays and the three neutrino flavors, for all the fermionic and bosonic channels and DM masses between 5 GeV and 100 TeV, on github.
|
De Romeri, V., Majumdar, A., Papoulias, D. K., & Srivastava, R. (2024). XENONnT and LUX-ZEPLIN constraints on DSNB-boosted dark matter. J. Cosmol. Astropart. Phys., 03(3), 028–34pp.
Abstract: We consider a scenario in which dark matter particles are accelerated to semirelativistic velocities through their scattering with the Diffuse Supernova Neutrino Background. Such a subdominant, but more energetic dark matter component can be then detected via its scattering on the electrons and nucleons inside direct detection experiments. This opens up the possibility to probe the sub -GeV mass range, a region of parameter space that is usually not accessible at such facilities. We analyze current data from the XENONnT and LUX-ZEPLIN experiments and we obtain novel constraints on the scattering cross sections of sub -GeV boosted dark matter with both nucleons and electrons. We also highlight the importance of carefully taking into account Earth's attenuation effects as well as the finite nuclear size into the analysis. By comparing our results to other existing constraints, we show that these effects lead to improved and more robust constraints.
|
Guadilla, V., Algora, A., Estienne, M., Fallot, M., Gelletly, W., Porta, A., et al. (2024). First measurements with a new fl-electron detector for spectral shape studies. J. Instrum., 19(2), P02027–21pp.
Abstract: The shape of the electron spectrum emitted in /3 decay carries a wealth of information about nuclear structure and fundamental physics. In spite of that, few dedicated measurements have been made of /3 -spectrum shapes. In this work we present a newly developed detector for /3 electrons based on a telescope concept. A thick plastic scintillator is employed in coincidence with a thin silicon detector. The first measurements employing this detector have been carried out with mono -energetic electrons from the high-energy resolution electron -beam spectrometer at Bordeaux. Here we report on the good reproduction of the experimental spectra of mono -energetic electrons using Monte Carlo simulations. This is a crucial step for future experiments, where a detailed Monte Carlo characterization of the detector is needed to determine the shape of the /3 -electron spectra by deconvolution of the measured spectra with the response function of the detector. A chamber to contain two telescope assemblies has been designed for future /3 -decay experiments at the Ion Guide Isotope Separator On -Line facility in Jyvaskyla, aimed at improving our understanding of reactor antineutrino spectra.
|
Hornillos, M. B. G., Gorlychev, V., Caballero, R., Cortes, G., Poch, A., Pretel, C., et al. (2011). Monte Carlo Simulations for the Study of a Moderated Neutron Detector. J. Korean Phys. Soc., 59(2), 1573–1576.
Abstract: This work presents the Monte Carlo simulations performed with the MCNPX and GEANT4 codes for the design of a BEta deLayEd Neutron detector, BELEN-20. This detector will be used for the study of beta delayed neutron emission and consists of a block of polyethylene with dimensions 90 x 90 x 80 cm(3) and 20 cylindrical (3)He gas counters. The results of these simulations have been validated experimentally with a (252)Cf source in the laboratory at UPC, Barcelona. Also the first experiment with this detector has been carried out in November 2009 in JYFL, Finland. In this experiment the neutron emission probability after beta decay of the fission products (88)Br, (94,95)Rb, and (138)I has been measured; this data is still under analysis. Simulations with MCNPX and GEANT4 have been performed in order to obtain the efficiency of the BELEN-20 detector for each of the above nuclei using the neutron energy distribution corresponding to each nucleus.
|
n_TOF Collaboration(Mendoza, E. et al), Giubrone, G., & Tain, J. L. (2011). Improved Neutron Capture Cross Section Measurements with the n_TOF Total Absorption Calorimeter. J. Korean Phys. Soc., 59(2), 1813–1816.
Abstract: The n_TOF collaboration operates a Total Absorption Calorimeter (TAC) [1] for measuring neutron capture cross-sections of low-mass and/or radioactive samples. The results obtained with the TAC have led to a substantial improvement of the capture cross sections of (237)Np and (240)Pu [2]. The experience acquired during the first measurements has allowed us to optimize the performance of the TAC and to improve the capture signal to background ratio, thus opening the way to more complex and demanding measurements on rare radioactive materials. The new design has been reached by a series of detailed Monte Carlo simulations of complete experiments and dedicated test measurements. The new capture setup will be presented and the main achievements highlighted.
|
