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: 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: 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: In this paper, we present steady-state relativistic magnetohydrodynamic simulations that include a mass-load term to study the process of jet deceleration. The mass load mimics the injection of a proton-electron plasma from stellar winds within the host galaxy into initially pair plasma jets, with mean stellar mass-losses ranging from 10(-14) to 10(-9) M-circle dot yr(-1). The spatial jet evolution covers similar to 500 pc from jet injection in the grid at 10 pc from the jet nozzle. Our simulations use a relativistic gas equation of state and a pressure profile for the ambient medium. We compare these simulations with previous dynamical simulations of relativistic, non-magnetized jets. Our results show that toroidal magnetic fields can prevent fast jet expansion and the subsequent embedding of further stars via magnetic tension. In this sense, magnetic fields avoid a runaway deceleration process. Furthermore, when the mass load is large enough to increase the jet density and produce fast, differential jet expansion, the conversion of magnetic energy flux into kinetic energy flux (i.e. magnetic acceleration), helps to delay the deceleration process with respect to non-magnetized jets. We conclude that the typical stellar population in elliptical galaxies cannot explain jet deceleration in classical Fanaroff-Riley type I radio galaxies. However, we observe a significant change in the jet composition, thermodynamical parameters, and energy dissipation along its evolution, even for moderate values of the mass load.

Abstract: Addressing the origin of the astrophysical neutrino flux observed by IceCube is of paramount importance. Gamma-Ray Bursts (GRBs) are among the few astrophysical sources capable of achieving the required energy to contribute to such neutrino flux through p gamma interactions. In this work, ANTARFS data have been used to search for upward going muon neutrinos in spatial and temporal coincidence with 784 GRBs occurred from 2007 to 2017. For each GRB, the expected neutrino flux has been calculated in the framework of the internal shock model and the impact of the lack of knowledge on the majority of source redshifts and on other intrinsic parameters of the emission mechanism has been quantified. It is found that the model parameters that set the radial distance where shock collisions occur have the largest impact on neutrino flux expectations. In particular, the bulk Lorentz factor of the source ejecta and the minimum variability time-scale are found to contribute significantly to the GRB-neutrino flux uncertainty. For the selected sources, ANTARES data have been analysed by maximizing the discovery probability of the stacking sample through an extended maximum-likelihood strategy. Since no neutrino event passed the quality cuts set by the optimization procedure, 90 per cent confidence level upper limits (with their uncertainty) on the total expected diffuse neutrino flux have been derived, according to the model. The GRB contribution to the observed diffuse astrophysical neutrino flux around 100 TeV is constrained to be less than 10 percent.