Abstract: Axions, if realized in nature, can be copiously produced in the early universe via thermal processes, contributing to the mass-energy density of thermal hot relics. In light of the most recent cosmological observations, we analyse two different thermal processes within a realistic mixed hot dark matter scenario which includes also massive neutrinos. Considering the axion-gluon thermalization channel, we derive our most constraining bounds on the hot relic masses m(a) < 7.46 eV and Sigma m(nu) < 0.114 eV both at 95 percent CL; while studying the axion-pion scattering, without assuming any specific model for the axion-pion interactions, and remaining in the range of validity of the chiral perturbation theory, our most constraining bounds are improved to m(a) < 0.91 eV and Sigma m(nu) < 0.105 eV, both at 95 percent CL. Interestingly, in both cases, the total neutrino mass lies very close to the inverted neutrino mass ordering prediction. If future terrestrial double beta decay and/or long-baseline neutrino experiments find that the nature mass ordering is the inverted one, this could rule out a wide region in the currently allowed thermal axion window. Our results therefore, strongly support multi messenger searches of axions and neutrino properties, together with joint analyses of their expected sensitivities.

Abstract: Around one-third of the point-like sources in the Fermi-LAT catalogues remain as unidentified sources (unIDs) today. Indeed, these unIDs lack a clear, univocal association with a known astrophysical source. If dark matter (DM) is composed of weakly interacting massive particles (WIMPs), there is the exciting possibility that some of these unIDs may actually be DM sources, emitting gamma-rays from WIMPs annihilation. We propose a new approach to solve the standard, machine learning (ML) binary classification problem of disentangling prospective DM sources (simulated data) from astrophysical sources (observed data) among the unIDs of the 4FGL Fermi-LAT catalogue. We artificially build two systematic features for the DM data which are originally inherent to observed data: the detection significance and the uncertainty on the spectral curvature. We do it by sampling from the observed population of unIDs, assuming that the DM distributions would, if any, follow the latter. We consider different ML models: Logistic Regression, Neural Network (NN), Naive Bayes, and Gaussian Process, out of which the best, in terms of classification accuracy, is the NN, achieving around 93 . 3 per cent +/- 0 . 7 per cent performance. Other ML evaluation parameters, such as the True Ne gativ e and True Positive rates, are discussed in our work. Applying the NN to the unIDs sample, we find that the de generac y between some astrophysical and DM sources can be partially solved within this methodology. None the less, we conclude that there are no DM source candidates among the pool of 4FGL Fermi-LAT unIDs.

Abstract: We analyse the f(R) gravity in the so-called Jordan frame, as implemented to the isotropic Universe dynamics. The goal of the present study is to show that according to recent data analyses of the supernovae Ia Pantheon sample, it is possible to account for an effective redshift dependence of the Hubble constant. This is achieved via the dynamics of a non-minimally coupled scalar field, as it emerges in the f(R) gravity. We face the question both from an analytical and purely numerical point of view, following the same technical paradigm. We arrive to establish that the expected decay of the Hubble constant with the redshift z is ensured by a form of the scalar field potential, which remains essentially constant for z less than or similar to 0.3, independently if this request is made a priori, as in the analytical approach, or obtained a posteriori, when the numerical procedure is addressed. Thus, we demonstrate that an f(R) dark energy model is able to account for an apparent variation of the Hubble constant due to the rescaling of the Einstein constant by the f(R) scalar mode.

Abstract: Developing sustainable software for the scientific community requires expertise in software engineering and domain science. This can be challenging due to the unique needs of scientific software, the insufficient resources for software engineering practices in the scientific community, and the complexity of developing for evolving scientific contexts. While open-source software can partially address these concerns, it can introduce complicating dependencies and delay development. These issues can be reduced if scientists and software developers collaborate. We present a case study wherein scientists from the SuperNova Early Warning System collaborated with software developers from the Scalable Cyberinfrastructure for Multi-Messenger Astrophysics project. The collaboration addressed the difficulties of open-source software development, but presented additional risks to each team. For the scientists, there was a concern of relying on external systems and lacking control in the development process. For the developers, there was a risk in supporting a user-group while maintaining core development. These issues were mitigated by creating a second Agile Scrum framework in parallel with the developers' ongoing Agile Scrum process. This Agile collaboration promoted communication, ensured that the scientists had an active role in development, and allowed the developers to evaluate and implement the scientists' software requirements. The collaboration provided benefits for each group: the scientists actuated their development by using an existing platform, and the developers utilized the scientists' use-case to improve their systems. This case study suggests that scientists and software developers can avoid scientific computing issues by collaborating and that Agile Scrum methods can address emergent concerns.

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.