Johannesson, G., Ruiz de Austri, R., Vincent, A. C., Moskalenko, I. V., Orlando, E., Porter, T. A., et al. (2016). Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion. Astrophys. J., 824(1), 16–19pp.
Abstract: We present the results of the most complete scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine-learning package. This is the first study to separate out low-mass isotopes (p, (p) over bar and He) from the usual light elements (Be, B, C, N, and O). We find that the propagation parameters that best-fit p, (p) over bar, and He data are significantly different from those that fit light elements, including the B/C and Be-10/Be-9 secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests that each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best-fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.
|
Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2023). Impact of QCD uncertainties on antiproton spectra from dark-matter annihilation. J. Cosmol. Astropart. Phys., 04(4), 068–15pp.
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.
|
Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2024). The Strong Force meets the Dark Sector: a robust estimate of QCD uncertainties for anti-matter dark matter searches. J. High Energy Phys., 02(2), 119–48pp.
Abstract: In dark-matter annihilation channels to hadronic final states, stable particles – such as positrons, photons, antiprotons, and antineutrinos – are produced via complex sequences of phenomena including QED/QCD radiation, hadronisation, and hadron decays. These processes are normally modelled by Monte Carlo (MC) event generators whose limited accuracy imply intrinsic QCD uncertainties on the predictions for indirect-detection experiments like Fermi-LAT, Pamela, IceCube or Ams-02. In this article, we perform a comprehensive analysis of QCD uncertainties, meaning both perturbative and nonperturbative sources of uncertainty are included – estimated via variations of MC renormalization-scale and fragmentation-function parameters, respectively – in antimatter spectra from dark-matter annihilation, based on parametric variations of the Pythia 8 event generator. After performing several retunings of light-quark fragmentation functions, we define a set of variations that span a conservative estimate of the QCD uncertainties. We estimate the effects on antimatter spectra for various annihilation channels and final-state particle species, and discuss their impact on fitted values for the dark-matter mass and thermally-averaged annihilation cross section. We find dramatic impacts which can go up to O(10%) for the annihilation cross section. We provide the spectra in tabulated form including QCD uncertainties and code snippets to perform fast dark-matter fits, in this github repository.
|
Menjo, H. et al, Faus-Golfe, A., & Velasco, J. (2011). Monte Carlo study of forward pi(0) production spectra to be measured by the LHCf experiment for the purpose of benchmarking hadron interaction models at 10(17) eV. Astropart Phys., 34(7), 513–520.
Abstract: The LHCf experiment aims to improve knowledge of forward neutral particle production spectra at the LHC energy which is relevant for the interpretation of air shower development of high energy cosmic rays. Two detectors, each composed of a pair of sampling and imaging calorimeters, have been installed at the forward region of IP1 to measure pi(0) energy spectra above 600 GeV. In this paper, we present a Monte Carlo study of the pi(0) measurements to be performed with one of the LHCf detectors for proton-proton collisions at root s = 14 TeV. In approximately 40 min of operation at luminosity 0.8 x 10(29) cm(-2) s(-1) during the beam commissioning phase of LHC, about 1.5 x 10(4) pi(0) events are expected to be obtained at two transverse positions of the detector. The backgrounds from interactions of secondary particles with beam pipes and interactions of beam particles with residual gas in the beam pipes are expected to be less than 0.1% of the signal from pi(0)s. We also discuss the capability of LHCf measurements to discriminate between the various hadron interaction models that are used for simulation of high energy air showers, such as DPMJET3.03, QGSJETII-03, SIBYLL2.1 and EPOS1.99.
|
Pierre Auger Collaboration(Aab, A. et al), & Pastor, S. (2014). A search for point sources of EeV photons. Astrophys. J., 789(2), 160–12pp.
Abstract: Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85 degrees to +20 degrees, in an energy range from 10(17.3) eV to 10(18.5) eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm(-2) s(-1), and no celestial direction exceeds 0.25 eV cm(-2) s(-1). These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.
|
Pierre Auger Collaboration(Aab, A. et al), & Pastor, S. (2014). A targeted search for point sources of EeV neutrons. Astrophys. J. Lett., 789(2), L34–7pp.
Abstract: A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine “target sets,” in addition to the search for a neutron flux from the Galactic center or from the Galactic plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.
|
Pierre Auger Collaboration(Aab, A. et al), & Pastor, S. (2014). Origin of atmospheric aerosols at the Pierre Auger Observatory using studies of air mass trajectories in South America. Atmos. Res., 149, 120–135.
Abstract: The Pierre Auger Observatory is making significant contributions towards understanding the nature and origin of ultra-high energy cosmic rays. One of its main challenges is the monitoring of the atmosphere, both in terms of its state variables and its optical properties. The aim of this work is to analyse aerosol optical depth tau(a)(z) values measured from 2004 to 2012 at the observatory, which is located in a remote and relatively unstudied area of Pampa Amarilla, Argentina. The aerosol optical depth is in average quite low – annual mean tau(a)(3.5 km) similar to 0.04 – and shows a seasonal trend with a winter minimum – tau(a)(3.5 km) – 0.03 -, and a summer maximum – tau(a)(3.5 km) similar to 0.06 -, and an unexpected increase from August to September tau(a)(35 km) similar to 0.055. We computed backward trajectories for the years 2005 to 2012 to interpret the air mass origin. Winter nights with low aerosol concentrations show air masses originating from the Pacific Ocean. Average concentrations are affected by continental sources (wind-blown dust and urban pollution), whilst the peak observed in September and October could be linked to biomass burning in the northern part of Argentina or air pollution coming from surrounding urban areas.
|
Pierre Auger Collaboration(Aab, A. et al), & Pastor, S. (2014). Reconstruction of inclined air showers detected with the pierre Auger Observatory. J. Cosmol. Astropart. Phys., 08(8), 019–32pp.
Abstract: We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than 60 degrees detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.
|
Pierre Auger Collaboration(Abraham, J. et al), & Pastor, S. (2010). The fluorescence detector of the Pierre Auger Observatory. Nucl. Instrum. Methods Phys. Res. A, 620(2-3), 227–251.
Abstract: The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.
|
Pierre Auger Collaboration(Abraham, J. et al), & Pastor, S. (2010). Measurement of the energy spectrum of cosmic rays above 10(18) eV using the Pierre Auger Observatory. Phys. Lett. B, 685(4-5), 239–246.
Abstract: We report a measurement of the flux of cosmic rays with unprecedented precision and Statistics using the Pierre Auger Observatory Based on fluorescence observations in coincidence with at least one Surface detector we derive a spectrum for energies above 10(18) eV We also update the previously published energy spectrum obtained with the surface detector array The two spectra are combined addressing the systematic uncertainties and, in particular. the influence of the energy resolution on the spectral shape The spectrum can be described by a broken power law E-gamma with index gamma = 3 3 below the ankle which is measured at log(10)(E-ankle/eV) = 18 6 Above the ankle the spectrum is described by a power law with index 2 6 followed by a flux suppression, above about log(10)(E/eV) = 19 5, detected with high statistical significance.
|
