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.

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.

Abstract: We present updated LHC limits on the minimal universal extra dimensions (MUEDs) model from the Run 2 searches. We scan the parameter space against a number of searches implemented in the public code CheckMATE and derive up-to-date limits on the MUED parameter space from 13TeV searches. The strongest constraints come from a search dedicated to squarks and gluinos with one isolated lepton, jets and missing transverse energy. In the procedure, we take into account initial state radiation and stress its importance in the MUED searches, which is not always appreciated.

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.

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.