Balaudo, A., Calore, F., De Romeri, V., & Donato, F. (2024). NAJADS: a self-contained framework for the direct determination of astrophysical J-factors. J. Cosmol. Astropart. Phys., 02(2), 001–33pp.
Abstract: Cosmological simulations play a pivotal role in understanding the properties of the dark matter (DM) distribution in both galactic and galaxy -cluster environments. The characterization of DM structures is crucial for informing indirect DM searches, aiming at the detection of the annihilation (or decay) products of DM particles. A fundamental quantity in these analyses is the astrophysical J -factor. In the DM phenomenology community, J -factors are typically computed through the semi -analytical modelling of the DM mass distribution, which is affected by large uncertainties. With the scope of addressing and possibly reducing these uncertainties, we present NAJADS, a self-contained framework to derive the DM J -factor directly from the raw simulations data. We show how this framework can be used to compute all -sky maps of the J -factor, automatically accounting for the complex 3D structure of the simulated halos and for the boosting of the signal due to the density fluctuations along the line of sight. After validating our code, we present a proof -of -concept application of NAJADS to a realistic halo from the IllustrisTNG suite, and exploit it to make a thorough comparison between our numerical approach and traditional semi -analytical methods. JCAP02(2024)001
|
Bernal, N., Forero-Romero, J. E., Garani, R., & Palomares-Ruiz, S. (2014). Systematic uncertainties from halo asphericity in dark matter searches. J. Cosmol. Astropart. Phys., 09(9), 004–30pp.
Abstract: Although commonly assumed to be spherical, dark matter halos are predicted to be non-spherical by N-body simulations and their asphericity has a potential impact on the systematic uncertainties in dark matter searches. The evaluation of these uncertainties is the main aim of this work, where we study the impact of aspherical dark matter density distributions in Milky-Way-like halos on direct and indirect searches. Using data from the large N-body cosmological simulation Bolshoi, we perform a statistical analysis and quantify the systematic uncertainties on the determination of local dark matter density and the so-called J factors for dark matter annihilations and decays from the galactic center. We find that, due to our ignorance about the extent of the non-sphericity of the Milky Way dark matter halo, systematic uncertainties can be as large as 35%, within the 95% most probable region, for a spherically averaged value for the local density of 0.3-0.4 GeV/cm(3). Similarly, systematic uncertainties on the J factors evaluated around the galactic center can be as large as 10% and 15%, within the 95% most probable region, for dark matter annihilations and decays, respectively.
|
Moline, A., Ibarra, A., & Palomares-Ruiz, S. (2015). Future sensitivity of neutrino telescopes to dark matter annihilations from the cosmic diffuse neutrino signal. J. Cosmol. Astropart. Phys., 06(6), 005–34pp.
Abstract: Cosmological observations and cold dark matter N-body simulations indicate that our Universe is populated by numerous halos, where dark matter particles annihilate, potentially producing Standard Model particles. In this paper we calculate the contribution to the diffuse neutrino background from dark matter annihilations in halos at all redshifts and we estimate the future sensitivity to the annihilation cross section of neutrino telescopes such as IceCube or ANTARES. We consider various parametrizations to describe the internal halo properties and for the halo mass function in order to bracket the theoretical uncertainty in the limits from the modeling of the cosmological annihilation flux. We find that observations of the cosmic diffuse neutrino flux at large angular distances from the galactic center lead to constraints on the dark matter annihilation cross section which are complementary to ( and for some extrapolations of the astrophysical parameters, better than) those stemming from observations of the Milky Way halo, especially for neutrino telescopes not pointing directly to the Milky Way center, as is the case of IceCube.
|
Achterberg, A., Amoroso, S., Caron, S., Hendriks, L., Ruiz de Austri, R., & Weniger, C. (2015). A description of the Galactic Center excess in the Minimal Supersymmetric Standard Model. J. Cosmol. Astropart. Phys., 08(8), 006–27pp.
Abstract: Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used to search for solutions within the pMSSM. We found solutions that are consistent with astroparticle physics and collider experiments, and provide a fit to the energy spectrum of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with a mass in the range 84-92 GeV or 87-97 GeV annihilating into W bosons. A third solutions is found for a neutralino of mass 174-187 GeV annihilating into top quarks. The best solutions yield a Dark Matter relic density 0.06 < Omega h(2) < 0.13. These pMSSM solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If the pMSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be discovered soon.
|
Amoroso, S., Caron, S., Jueid, A., Ruiz de Austri, R., & Skands, P. (2019). Estimating QCD uncertainties in Monte Carlo event generators for gamma-ray dark matter searches. J. Cosmol. Astropart. Phys., 05(5), 007–44pp.
Abstract: Motivated by the recent galactic center gamma-ray excess identified in the Fermi-LAT data, we perform a detailed study of QCD fragmentation uncertainties in the modeling of the energy spectra of gamma-rays from Dark-Matter (DM) annihilation. When Dark-Matter particles annihilate to coloured final states, either directly or via decays such as W(*) -> qq-', photons are produced from a complex sequence of shower, hadronisation and hadron decays. In phenomenological studies their energy spectra are typically computed using Monte Carlo event generators. These results have however intrinsic uncertainties due to the specific model used and the choice of model parameters, which are difficult to asses and which are typically neglected. We derive a new set of hadronisation parameters (tunes) for the PYTHIA 8.2 Monte Carlo generator from a fit to LEP and SLD data at the Z peak. For the first time we also derive a conservative set of uncertainties on the shower and hadronisation model parameters. Their impact on the gamma-ray energy spectra is evaluated and discussed for a range of DM masses and annihilation channels. The spectra and their uncertainties are also provided in tabulated form for future use. The fragmentation-parameter uncertainties may be useful for collider studies as well.
|
n_TOF Collaboration(Alcayne, V. et al), Balibrea-Correa, J., Domingo-Pardo, C., Lerendegui-Marco, J., Babiano-Suarez, V., & Ladarescu, I. (2024). A Segmented Total Energy Detector (sTED) optimized for (n,γ) cross-section measurements at n_TOF EAR2. Radiat. Phys. Chem., 217, 11pp.
Abstract: The neutron time-of-flight facility nTOF at CERN is a spallation source dedicated to measurements of neutroninduced reaction cross-sections of interest in nuclear technologies, astrophysics, and other applications. Since 2014, Experimental ARea 2 (EAR2) is operational and delivers a neutron fluence of similar to 4 center dot 10(7) neutrons per nominal proton pulse, which is similar to 50 times higher than the one of Experimental ARea 1 (EAR1) of similar to 8 center dot 10(5) neutrons per pulse. The high neutron flux at EAR2 results in high counting rates in the detectors that challenged the previously existing capture detection systems. For this reason, a Segmented Total Energy Detector (sTED) has been developed to overcome the limitations in the detector's response, by reducing the active volume per module and by using a photo-multiplier (PMT) optimized for high counting rates. This paper presents the main characteristics of the sTED, including energy and time resolution, response to gamma-rays, and provides as well details of the use of the Pulse Height Weighting Technique (PHWT) with this detector. The sTED has been validated to perform neutron-capture cross-section measurements in EAR2 in the neutron energy range from thermal up to at least 400 keV. The detector has already been successfully used in several measurements at nTOF EAR2.
|
Caron, S., Eckner, C., Hendriks, L., Johannesson, G., Ruiz de Austri, R., & Zaharijas, G. (2023). Mind the gap: the discrepancy between simulation and reality drives interpretations of the Galactic Center Excess. J. Cosmol. Astropart. Phys., 06(6), 013–56pp.
Abstract: The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model (-yEM) used in Galactic center analyses affects the interpretation of the GCE's nature. To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of-yEMs with increasing parametric freedom. We calculate the fractional contribution (fsrc) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the ryEM. For the simplest ryEM, we obtain fsrc = 0.10 f 0.07, while the most complex model yields fsrc = 0.79 f 0.24. In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed ryEM. The quoted results for fsrc do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated ryEM iterations. We quantify the reality gap between our ryEMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed ryEMs have gaps to reality. Our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked “out of domain” uncertainties in future interpretations.
|
Mertsch, P., Parimbelli, G., de Salas, P. F., Gariazzo, S., Lesgourgues, J., & Pastor, S. (2020). Neutrino clustering in the Milky Way and beyond. J. Cosmol. Astropart. Phys., 01(1), 015–23pp.
Abstract: The standard cosmological model predicts the existence of a Cosmic Neutrino Background, which has not yet been observed directly. Some experiments aiming at its detection are currently under development, despite the tiny kinetic energy of the cosmological relic neutrinos, which makes this task incredibly challenging. Since massive neutrinos are attracted by the gravitational potential of our Galaxy, they can cluster locally. Neutrinos should be more abundant at the Earth position than at an average point in the Universe. This fact may enhance the expected event rate in any future experiment. Past calculations of the local neutrino clustering factor only considered a spherical distribution of matter in the Milky Way and neglected the influence of other nearby objects like the Virgo cluster, although recent N-body simulations suggest that the latter may actually be important. In this paper, we adopt a back-tracking technique, well established in the calculation of cosmic rays fluxes, to perform the first three-dimensional calculation of the number density of relic neutrinos at the Solar System, taking into account not only the matter composition of the Milky Way, but also the contribution of the Andromeda galaxy and the Virgo cluster. The effect of Virgo is indeed found to be relevant and to depend non-trivially on the value of the neutrino mass. Our results show that the local neutrino density is enhanced by 0.53% for a neutrino mass of 10 meV, 12% for 50 meV, 50% for 100 meV or 500% for 300 meV.
|
Tain, J. L., Agramunt, J., Algora, A., Aprahamian, A., Cano-Ott, D., Fraile, L. M., et al. (2015). The sensitivity of LaBr3:Ce scintillation detectors to low energy neutrons: Measurement and Monte Carlo simulation. Nucl. Instrum. Methods Phys. Res. A, 774, 17–24.
Abstract: The neutron sensitivity of a cylindrical circle minus 1.5 in x 1.5 in LaBr3:Ce scintillation detector was measured using quasi-monoenergetic neutron beams in the energy range from 40 keV to 2.5 MeV. In this energy range the detector is sensitive to gamma-rays generated in neutron inelastic and capture processes. The experimental energy response was compared with Monte Carlo simulations performed with the Geant4 simulation toolkit using the so-called High Precision Neutron Models. These models rely on relevant information stored in evaluated nuclear data libraries. The performance of the Geant4 Neutron Data Library as well as several standard nuclear data libraries was investigated. In the latter case this was made possible by the use of a conversion tool that allowed the direct use of the data from other libraries in Geant4. Overall it was found that there was good agreement with experiment for some of the neutron data bases like ENDF/B-VII.0 or JENDL-3.3 but not with the others such as ENDF/B-VI.8 or JEFF-3.1.
|
de Putter, R., Wagner, C., Mena, O., Verde, L., & Percival, W. J. (2012). Thinking outside the box: effects of modes larger than the survey on matter power spectrum covariance. J. Cosmol. Astropart. Phys., 04(4), 019–31pp.
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.
|
