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.
|
Stoppa, F., Bhattacharyya, S., Ruiz de Austri, R., Vreeswijk, P., Caron, S., Zaharijas, G., et al. (2023). AutoSourceID-Classifier Star-galaxy classification using a convolutional neural network with spatial information. Astron. Astrophys., 680, A109–16pp.
Abstract: Aims. Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification's reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images.Methods. The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results.Results. We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C's direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.
|
Beenakker, W., Caron, S., Kip, J., Ruiz de Austri, R., & Zhang, Z. (2023). New energy spectra in neutrino and photon detectors to reveal hidden dark matter signals. J. High Energy Phys., 11(11), 028–13pp.
Abstract: Neutral particles capable of travelling cosmic distances from a source to detectors on Earth are limited to photons and neutrinos. Examination of the Dark Matter annihilation/decay spectra for these particles reveals the presence of continuum spectra (e.g. due to fragmentation and W or Z decay) and peaks (due to direct annihilations/decays). However, when one explores extensions of the Standard Model (BSM), unexplored spectra emerge that differ significantly from those of the Standard Model (SM) for both neutrinos and photons. In this paper, we argue for the inclusion of important spectra that include peaks as well as previously largely unexplored entities such as boxes and combinations of box, peak and continuum decay spectra.
|
Stoppa, F., Ruiz de Austri, R., Vreeswijk, P., Bhattacharyya, S., Caron, S., Bloemen, S., et al. (2023). AutoSourceID-FeatureExtractor Optical image analysis using a two-step mean variance estimation network for feature estimation and uncertainty characterisation. Astron. Astrophys., 680, A108–14pp.
Abstract: Aims. In astronomy, machine learning has been successful in various tasks such as source localisation, classification, anomaly detection, and segmentation. However, feature regression remains an area with room for improvement. We aim to design a network that can accurately estimate sources' features and their uncertainties from single-band image cutouts, given the approximated locations of the sources provided by the previously developed code AutoSourceID-Light (ASID-L) or other external catalogues. This work serves as a proof of concept, showing the potential of machine learning in estimating astronomical features when trained on meticulously crafted synthetic images and subsequently applied to real astronomical data.Methods. The algorithm presented here, AutoSourceID-FeatureExtractor (ASID-FE), uses single-band cutouts of 32x32 pixels around the localised sources to estimate flux, sub-pixel centre coordinates, and their uncertainties. ASID-FE employs a two-step mean variance estimation (TS-MVE) approach to first estimate the features and then their uncertainties without the need for additional information, for example the point spread function (PSF). For this proof of concept, we generated a synthetic dataset comprising only point sources directly derived from real images, ensuring a controlled yet authentic testing environment.Results. We show that ASID-FE, trained on synthetic images derived from the MeerLICHT telescope, can predict more accurate features with respect to similar codes such as SourceExtractor and that the two-step method can estimate well-calibrated uncertainties that are better behaved compared to similar methods that use deep ensembles of simple MVE networks. Finally, we evaluate the model on real images from the MeerLICHT telescope and the Zwicky Transient Facility (ZTF) to test its transfer learning abilities.
|
Ferrer-Sanchez, A., Martin-Guerrero, J., Ruiz de Austri, R., Torres-Forne, A., & Font, J. A. (2024). Gradient-annihilated PINNs for solving Riemann problems: Application to relativistic hydrodynamics. Comput. Meth. Appl. Mech. Eng., 424, 116906–18pp.
Abstract: We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that forces the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting function. The method relies on a set of hyperparameters that control how gradients are treated in the physical loss. The performance of our methodology is demonstrated by solving Riemann problems in special relativistic hydrodynamics, extending earlier studies with PINNs in the context of the classical Euler equations. The solutions obtained with the GA-PINN model correctly describe the propagation speeds of discontinuities and sharply capture the associated jumps. We use the relative l(2) error to compare our results with the exact solution of special relativistic Riemann problems, used as the reference ''ground truth'', and with the corresponding error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by the GA-PINN model is comparable to that obtained with a shock-capturing scheme, achieving a performance superior to that of the baseline PINN algorithm in general. An additional benefit worth stressing is that our PINN-based approach sidesteps the costly recovery of the primitive variables from the state vector of conserved variables, a well-known drawback of grid-based solutions of the relativistic hydrodynamics equations. Due to its inherent generality and its ability to handle steep gradients, the GA-PINN methodology discussed in this paper could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions.
|
Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2024). The Strong Force meets the Dark Sector: a robust estimate of QCD uncertainties for anti-matter dark matter searches. J. High Energy Phys., 02(2), 119–48pp.
Abstract: In dark-matter annihilation channels to hadronic final states, stable particles – such as positrons, photons, antiprotons, and antineutrinos – are produced via complex sequences of phenomena including QED/QCD radiation, hadronisation, and hadron decays. These processes are normally modelled by Monte Carlo (MC) event generators whose limited accuracy imply intrinsic QCD uncertainties on the predictions for indirect-detection experiments like Fermi-LAT, Pamela, IceCube or Ams-02. In this article, we perform a comprehensive analysis of QCD uncertainties, meaning both perturbative and nonperturbative sources of uncertainty are included – estimated via variations of MC renormalization-scale and fragmentation-function parameters, respectively – in antimatter spectra from dark-matter annihilation, based on parametric variations of the Pythia 8 event generator. After performing several retunings of light-quark fragmentation functions, we define a set of variations that span a conservative estimate of the QCD uncertainties. We estimate the effects on antimatter spectra for various annihilation channels and final-state particle species, and discuss their impact on fitted values for the dark-matter mass and thermally-averaged annihilation cross section. We find dramatic impacts which can go up to O(10%) for the annihilation cross section. We provide the spectra in tabulated form including QCD uncertainties and code snippets to perform fast dark-matter fits, in this github repository.
|
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.
|
Caron, S., Ruiz de Austri, R., & Zhang, Z. Y. (2023). Mixture-of-Theories training: can we find new physics and anomalies better by mixing physical theories? J. High Energy Phys., 03(3), 004–37pp.
Abstract: Model-independent search strategies have been increasingly proposed in recent years because on the one hand there has been no clear signal for new physics and on the other hand there is a lack of a highly probable and parameter-free extension of the standard model. For these reasons, there is no simple search target so far. In this work, we try to take a new direction and ask the question: bearing in mind that we have a large number of new physics theories that go beyond the Standard Model and may contain a grain of truth, can we improve our search strategy for unknown signals by using them “in combination”? In particular, we show that a signal hypothesis based on a large, intermingled set of many different theoretical signal models can be a superior approach to find an unknown BSM signal. Applied to a recent data challenge, we show that “mixture-of-theories training” outperforms strategies that optimize signal regions with a single BSM model as well as most unsupervised strategies. Applications of this work include anomaly detection and the definition of signal regions in the search for signals of new physics.
|
van Beekveld, M., Beenakker, W., Caron, S., Kip, J., Ruiz de Austri, R., & Zhang, Z. Y. (2023). Non-standard neutrino spectra from annihilating neutralino dark matter. SciPost Phys. Core, 6(1), 006–23pp.
Abstract: Neutrino telescope experiments are rapidly becoming more competitive in indirect de-tection searches for dark matter. Neutrino signals arising from dark matter annihilations are typically assumed to originate from the hadronisation and decay of Standard Model particles. Here we showcase a supersymmetric model, the BLSSMIS, that can simulta-neously obey current experimental limits while still providing a potentially observable non-standard neutrino spectrum from dark matter annihilation.
|
Kim, J. S., Lopez-Fogliani, D. E., Perez, A. D., & Ruiz de Austri, R. (2023). Right-handed sneutrino and gravitino multicomponent dark matter in light of neutrino detectors. J. Cosmol. Astropart. Phys., 04(4), 050–32pp.
Abstract: We investigate the possibility that right-handed (RH) sneutrinos and gravitinos can coexist and explain the dark matter (DM) problem. We compare extensions of the minimal supersymmetric standard model (MSSM) and the next-to-MSSM (NMSSM) adding RH neutrinos superfields, with special emphasis on the latter. If the gravitino is the lightest supersymmetric particle (LSP) and the RH sneutrino the next-to-LSP (NLSP), the heavier particle decays to the former plus left-handed (LH) neutrinos through the mixing between the scalar partners of the LH and RH neutrinos. However, the interaction is suppressed by the Planck mass, and if the LH-RH sneutrino mixing parameter is small, << O(10-2), a long-lived RH sneutrino NLSP is possible even surpassing the age of the Universe. As a byproduct, the NLSP to LSP decay produces monochromatic neutrinos in the ballpark of current and planned neutrino telescopes like Super-Kamiokande, IceCube and Antares that we use to set constraints and show prospects of detection. In the NMSSM+RHN, assuming a gluino mass parameter M3 = 3 TeV we found the following lower limits for the gravitino mass m3/2 >= 1-600 GeV and the reheating temperature TR >= 105-3 x 107 GeV, for m nu similar to R similar to 10-800 GeV. If we take M3 = 10 TeV, then the limits on TR are relaxed by one order of magnitude.
|
