We assess Galactic Dark Matter (DM) sensitivities to photons from annihilation and decay using the spatial and kinematic information determined by state-of-the-art simulations in the Latte suite of Feedback In Realistic Environments (FIRE-2). For kinematic information, we study the energy shift pattern of DM narrow emission lines predicted in FIRE-2 and discuss its potential as DM-signal diagnosis, showing for the first time the power of symmetric observations around $l=0^{\circ}$. We find that the exposures needed to resolve the line separation of DM to gas by XRISM at $5\sigma$ to be large, $\gtrsim 4$ Ms, while exposures are smaller for Athena ($\lesssim 50$ ks) and Lynx ($\lesssim 100$ ks). We find that large field-of-view exposures remain the most sensitive methods for detection of DM annihilation or decay by the luminosity of signals in the field of view dominating velocity information. The $\sim$4 sr view of the Galactic Center region by the Wide Field Monitor (WFM) aboard the eXTP mission will be highly sensitive to DM signals, with a prospect of $\sim 10^5$ to $10^6$ events from the 3.5 keV line in a 100 ks exposure, with the range dependent on photon acceptance in WFM's field of view. We also investigate detailed all-sky luminosity maps for both DM annihilation and decay signals - evaluating the signal-to-noise for a DM detection with realistic X-ray and gamma-ray backgrounds - as a guideline for what could be a forthcoming era of DM astronomy.

We point out an anti-correlation between the central dark matter (DM) densities of the bright Milky Way dwarf spheroidal galaxies (dSphs) and their orbital pericenter distances inferred from Gaia data. The dSphs that have not come close to the Milky Way center (like Fornax, Carina and Sextans) are less dense in DM than those that have come closer (like Draco and Ursa Minor). The same anti-correlation cannot be inferred for the ultra-faint dSphs due to large scatter. Including ultra-faints, a trend that dSphs with more extended stellar distributions tend to have lower DM densities emerges. A fresh look at solutions to the too-big-to-fail problem is warranted in light of these observations.

Milky Way dwarf spheroidal satellites are a prime target for Dark Matter (DM) indirect searches. There have been recent reassessments of the expected DM gamma-ray signals in case of long-range interactions, commonly known as Sommerfeld enhancement. Since details of the underlying DM phase-space distribution function become critical, there are potentially large uncertainties in the final result. We provide here a first attempt towards a comprehensive investigation of these systematics, addressing the impact on the expected DM flux from Milky Way dwarfs via Bayesian inference on the available stellar kinematic datasets. After reconsidering the study case of ergodic systems, we investigate for the first time scenarios where DM particle orbits may have a radial or tangential bias. We consider both cuspy and cored parametric DM density profiles, together with the case of a non-parametric halo modelling directly connected to line-of-sight observable quantities. The main findings of our work highlight the relevance of the assumed phase-space distribution: Referring to a generalized J-factor, namely the line-of-sight convolution of the spatial part in case of velocity-dependent annihilation rate, an enhancement (suppression) with respect to the limit of isotropic phase-space distributions is obtained for the case of tangentially (radially) biased DM particle orbits. We provide new estimates for J-factors for the eight brightest Milky Way dwarfs also in the limit of velocity-independent DM annihilation, in good agreement with previous results in literature, and derive data-driven lower-bounds based on the non-parametric modelling of the halo density. The outcome of our broad study stands out as a representative of the state-of-the-art in the field, and falls within the interest of current and future experimental collaborations involved in DM indirect detection programs.

Dwarf spheroidal galaxies provide well-known challenges to the standard cold and collisionless dark matter scenario: The too-big-to-fail problem, namely the mismatch between the observed mass enclosed within the half-light radius of dwarf spheroidals and cold dark matter N-body predictions; The hints for inner constant-density cores. While these controversies may be alleviated by baryonic physics and environmental effects, revisiting the standard lore of cold and collisionless dark matter remains an intriguing possibility. Self-interacting dark matter may be the successful proposal to such a small-scale crisis. Self-interactions correlate dark matter and baryon distributions, allowing for constant-density cores in low surface brightness galaxies. Here we report the first data-driven study of the too-big-to-fail of Milky Way dwarf spheroidals within the self-interacting dark matter paradigm. We find good description of stellar kinematics and compatibility with the concentration-mass relation from recent pure cold dark matter simulations. Within this concentration-mass relation, a subset of Milky Way dwarfs are well fitted by cross sections of 0.5-3 cm$^2$/g, while others point to values greater than 10 cm$^2$/g.

The quest for the elusive dark matter (DM) that permeates the Universe (and in general the search for signatures of Physics beyond the Standard Model at astronomical scales) provides a unique opportunity and a tough challenge to the high-energy astrophysics community. In particular, the so-called DM indirect searches - mostly focused on a class of theoretically well-motivated DM candidates such as the weakly-interacting massive particles - are affected by a complex astrophysical background of cosmic radiation. The understanding and modeling of such background requires a deep comprehension of an intricate classical plasma physics problem, i.e. the interaction between high-energy charged particles, accelerated in peculiar astrophysical environments, and magneto-hydrodynamic turbulence in the interstellar medium of our Galaxy. In this review we highlight several aspects of this exciting interplay between the most recent claims of DM annihilation/decay signatures from the sky and the Galactic cosmic-ray research field. Our purpose is to further stimulate the debate about viable astrophysical explanations, discussing possible directions that would help breaking degeneracy patterns in the interpretation of current data. We eventually aim to emphasize how a deep knowledge on the physics of CR transport is therefore required to tackle the DM indirect search program at present and in the forthcoming years.

Within the standard approach of effective field theory of weak interactions for $\Delta B = 1$ transitions, we look for possibly unexpected subtle New Physics effects, here dubbed "flavourful Easter eggs". We perform a Bayesian global fit using the publicly available HEPfit package, taking into account state-of-the-art experimental information concerning these processes, including the suggestive measurements from LHCb of $R_{K}$ and $R_{K^{*}}$, the latter available only very recently. We parametrize New Physics contributions to $b \to s$ transitions in terms of shifts of Wilson coefficients of the electromagnetic dipole and semi-leptonic operators, assuming CP-conserving effects, but allowing in general for violation of lepton flavour universality. We show how optimistic/conservative hadronic estimates can impact quantitatively the size of New Physics extracted from the fit. With a conservative approach to hadronic uncertainties we find nonzero New Physics contributions to Wilson coefficients at the level of $\sim 3\sigma$, depending on the model chosen. Furthermore, given the interplay between hadronic contributions and New Physics effects in the leptonic vector current, a scenario with nonstandard leptonic axial currents is comparable to the more widely advocated one with New Physics in the leptonic vector current.

Dwarf satellite galaxies are ideal laboratories for identifying particle Dark Matter signals. When setting limits on particle Dark Matter properties from null searches, it becomes however crucial the level at which the Dark Matter density profile within these systems is constrained by observations. In the limit in which the spherical Jeans equation is assumed to be valid for a given tracer stellar population, we study the solution of this equation having the Dark Matter mass profile as an output rather than as a trial parametric input. Within our new formulation, we address to what level dwarf spheroidal galaxies feature a reliable mass estimator. We assess then possible extrapolation of the density profiles in the inner regions and -- keeping explicit the dependence on the orbital anisotropy profile of the tracer population -- we derive general trends on the line-of-sight integral of the density profile squared, a quantity commonly dubbed $J$-factor and crucial to estimate fluxes from prompt Dark Matter pair annihilations. Taking Ursa Minor as a study case among Milky Way satellites, we perform Bayesian inference using the available kinematical data for this galaxy. Contrary to all previous studies, we avoid marginalization over quantities poorly constrained by observations or by theoretical arguments. We find minimal $J$-factors to be about 2 to 4 times smaller than commonly quoted estimates, approximately relaxing by the same amount the limit on Dark Matter pair annihilation cross section from gamma-ray surveys of Ursa Minor. At the same time, if one goes back to a fixed trial parametric form for the density, e.g. using a NFW or Burkert profile, we show that the minimal $J$ can hardly be reduced by more than a factor of 1.5.

The last IceCube catalog of High Energy Starting Events (HESE) obtained with a livetime of 1347 days comprises 54 neutrino events equally-distributed between the three families with energies between 25 TeV and few PeVs. Considering the homogeneous flavors distribution (1:1:1) and the spectral features of these neutrinos the IceCube collaboration claims the astrophysical origin of these events with more than $5\sigma$. The spatial distribution of cited events does not show a clear correlation with known astrophysical accelerators leaving opened both the Galactic and the extra-Galactic origin interpretations. Here, we compute the neutrino diffuse emission of our Galaxy on the basis of a recently proposed phenomenological model characterized by radially-dependent cosmic-ray (CR) transport properties. We show that the astrophysical spectrum measured by IceCube experiment can be well explained adding to the diffuse Galactic neutrino flux (obtained with this new model) a extra-Galactic component derived from the astrophysical muonic neutrinos reconstructed in the Northern hemisphere. A good agreement between the expected astrophysical neutrino flux and the IceCube data is found for the full sky as well as for the Galactic plane region.

A spherical-symmetric gamma-ray emission from the central region of the Galaxy has been recently identified in Fermi-LAT data, and initially associated to dark matter particle annihilations. Guided by the evidence for a high gas density in the inner kpc of the Galaxy correlated with a very large Supernova rate, and hence with ongoing cosmic-ray acceleration, we investigate instead the possibility of addressing this excess in terms of ordinary cosmic-ray sources and standard steady-state diffusion. We alter the source term, and consistently the correlated gamma-ray emissions, in the context of a template-fitting analysis. We focus on a region of interest (ROI) defined as: $|l|<20^{\circ}; \hspace{5mm} 2^{\circ}<|b|<20^{\circ},$ with $l$ and $b$ the Galactic longitude and latitude coordinates. We analyze in detail the overall goodness of the fit of our framework, and perform a detailed direct comparison against data examining profiles in different directions. Remarkably, the test statistic of the fit related to our scenario turns out to be as good as the Dark Matter one in the ROI here considered.

We compute the gamma-ray and neutrino diffuse emission of the Galaxy on the basis of a recently proposed phenomenological model characterized by radially dependent cosmic-ray (CR) transport properties. We show how this model, designed to reproduce both Fermi-LAT gamma-ray data and local CR observables, naturally reproduces the anomalous TeV diffuse emission observed by Milagro in the inner Galactic plane. Above 100 TeV our picture predicts a neutrino flux that is about five (two) times larger than the neutrino flux computed with conventional models in the Galactic Center region (full-sky). Explaining in that way up to $\sim 25 \%$ of the flux measured by IceCube, we reproduce the full-sky IceCube spectrum adding an extra-Galactic component derived from the muonic neutrino flux in the northern hemisphere. We also present precise predictions for the Galactic plane region where the flux is dominated by the Galactic emission.

Several observations from Fermi-LAT, up to few hundred GeV, and from H.E.S.S., up to $\sim$ 10 TeV, reported an intense $\gamma$-ray emission from the inner part of the Galactic plane. After the subtraction of point-like contributions, the remaining $\gamma$-ray spectrum can provide important hints about the cosmic-ray (CR) population in that region. In particular, the diffuse spectrum measured by both Fermi-LAT and H.E.S.S. in the Galactic Ridge is significantly harder with respect to the rest of the Galaxy. These results were recently interpreted in terms of a comprehensive CR transport model which, adopting a spatial dependent diffusion coefficient and convective velocity, reproduces Fermi-LAT results on the whole sky as well as local CR spectra. We showed as that model predicts a significantly harder neutrino diffuse emission compared to conventional scenarios: The predicted signal is able to account for a significant fraction of the astrophysical flux measured by IceCube. In this contribution, we use the same setup to calculate the expected neutrino flux from several windows in the inner Galactic plane and compare the results with IceCube observations and the sensitivities of Mediterranean neutrino telescopes. In particular, for the ANTARES experiment, we compare the model expectations with the upper limits obtained from a recent unblinded data-analysis focused on the galactic ridge region. Moreover, we also show the expectations from the galactic ridge for the future KM3NeT observatory, whose position is optimal to observe this portion of the sky.

As recently shown, Fermi-LAT measurements of the diffuse gamma-ray emission from the Galaxy favor the presence of a smooth softening in the primary cosmic-ray spectrum with increasing Galactocentric distance. This result can be interpreted in terms of a spatial-dependent rigidity scaling of the diffusion coefficient. The DRAGON code was used to build a model based on such feature. That scenario correctly reproduces the latest Fermi-LAT results as well as local cosmic-ray measurements from PAMELA, AMS-02 and CREAM. Here we show that the model, if extrapolated at larger energies, grasps both the gamma-ray flux measured by MILAGRO at 15 TeV and the H.E.S.S. data from the Galactic ridge, assuming that the cosmic-ray spectral hardening found by those experiments at about 250 GeV/n is present in the whole inner Galactic plane region. Moreover, we show as that model also predicts a neutrino emission which may account for a significant fraction, as well as for the correct spectral shape, of the astrophysical flux measured by IceCube above 25 TeV.

The standard approach to cosmic-ray (CR) propagation in the Galaxy is based on the assumption that local transport properties can be extrapolated to the whole CR confining volume. Such models tend to underestimate the gamma-ray flux above few GeV measured by the Fermi Large Area Telescope towards the inner Galactic plane. We consider here for the first time a phenomenological scenario allowing for both the rigidity scaling of the diffusion coefficient and convective effects to be position-dependent. We show that within this approach we can reproduce the observed gamma-ray spectra at both low and mid Galactic latitudes - including the Galactic center - without spoiling any local CR observable.

The calculation of transport properties of Fermi liquids, based on the formalism developed by Abrikosov and Khalatnikov, requires the knowledge of the probability of collisions between quasiparticles in the vicinity of the Fermi surface. We have carried out a numerical study of the shear viscosity of pure neutron matter, whose value plays a pivotal role in determining the stability of rotating neutron stars, in which these processes are described using a state-of-the-art nucleon-nucleon potential model. Within our approach medium modifications of the scattering cross section are consistently taken into account, through an effective interaction obtained from the matrix elements of the bare interaction between correlated states. Inclusion of medium effects lead to a large increase of the viscosity at densities larger than $\sim 0.1$ fm^{-3}.