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.
|
Dorigo, T. et al, Ramos, A., & Ruiz de Austri, R. (2023). Toward the end-to-end optimization of particle physics instruments with differentiable programming. Rev. Phys., 10, 100085– pp.
Abstract: The full optimization of the design and operation of instruments whose functioning relies on the interaction of radiation with matter is a super-human task, due to the large dimensionality of the space of possible choices for geometry, detection technology, materials, data-acquisition, and information-extraction techniques, and the interdependence of the related parameters. On the other hand, massive potential gains in performance over standard, “experience-driven” layouts are in principle within our reach if an objective function fully aligned with the final goals of the instrument is maximized through a systematic search of the configuration space. The stochastic nature of the involved quantum processes make the modeling of these systems an intractable problem from a classical statistics point of view, yet the construction of a fully differentiable pipeline and the use of deep learning techniques may allow the simultaneous optimization of all design parameters.
|
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.
|
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.
|
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.
|
