Abstract: Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron de-cays, before particles that are stable on astrophysical time scales are produced. Antiprotons produced in this way may leave footprints in experiments such as AMS-02. Several groups have reported an excess of events in the antiproton flux in the rigidity range of 10-20 GV. However, the theoretical modeling of baryon production is not straightforward and relies in part on phenomenological models in Monte Carlo event generators. In this work, we assess the impact of QCD uncertainties on the spectra of antiprotons from dark-matter annihila-tion. As a proof-of-principle, we show that for a two-parameter model that depends only on the thermally-averaged annihilation cross section ((o -v)) and the dark-matter mass (Mx), QCD uncertainties can affect the best-fit mass by up to ti 14% (with large uncertainties for large DM masses), depending on the choice of Mx and the annihilation channel (bb over bar or W+W-), and (o -v) by up to ti 10%. For comparison, changes to the underlying diffusion parameters are found to be within 1%-5%, and the results are also quite resilient to the choice of cosmic-ray propagation model. These findings indicate that QCD uncertainties need to be included in future DM analyses. To facilitate full-fledged analyses, we provide the spectra in tabulated form including QCD uncertainties and code snippets to perform mass interpolations and quick DM fits. The code can be found in this GitHub [1] repository.

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.

Abstract: We consider a scenario in which dark matter particles are accelerated to semirelativistic velocities through their scattering with the Diffuse Supernova Neutrino Background. Such a subdominant, but more energetic dark matter component can be then detected via its scattering on the electrons and nucleons inside direct detection experiments. This opens up the possibility to probe the sub -GeV mass range, a region of parameter space that is usually not accessible at such facilities. We analyze current data from the XENONnT and LUX-ZEPLIN experiments and we obtain novel constraints on the scattering cross sections of sub -GeV boosted dark matter with both nucleons and electrons. We also highlight the importance of carefully taking into account Earth's attenuation effects as well as the finite nuclear size into the analysis. By comparing our results to other existing constraints, we show that these effects lead to improved and more robust constraints.

Abstract: We examine the cosmology of warm dark matter (WDM), both stable and decaying, from the point of view of structure formation. We compare the matter power spectrum associated to WDM masses of 1.5 keV and 0.158 keV, with that expected for the stable cold dark matter ACDM Xi SCDM paradigm, taken as our reference model. We scrutinize the effects associated to the warm nature of dark matter, as well as the fact that it decays. The decaying warm dark matter (DWDM) scenario is well-motivated, emerging in a broad class of particle physics theories where neutrino masses arise from the spontaneous breaking of a continuous global lepton number symmetry. The majoron arises as a Nambu-Goldstone boson, and picks up a mass from gravitational effects, that explicitly violate global symmetries. The majoron necessarily decays to neutrinos, with an amplitude proportional to their tiny mass, which typically gives it cosmologically long lifetimes. Using N-body simulations we show that our DWDM picture leads to a viable alternative to the ACDM scenario, with predictions that can differ substantially on small scales.

Abstract: Accurate power spectrum (or correlation function) covariance matrices are a crucial requirement for cosmological parameter estimation from large scale structure surveys. In order to minimize reliance on computationally expensive mock catalogs, it is important to have a solid analytic understanding of the different components that make up a covariance matrix. Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This beat coupling effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated average density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this local average effect and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only similar to 10% of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations, using both full simulation boxes and subboxes thereof to study cases without beat coupling, with beat coupling and with both beat coupling and the local average effect. For the variances, we find excellent agreement with the analytic expressions for k < 0.2 hMpc(-1) at z = 0.5, while the correlation coefficients agree to beyond k = 0.4 hMpc(-1). As expected, the range of agreement increases towards higher redshift and decreases slightly towards z = 0. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix.