We present an update on CRDB (https://lpsc.in2p3.fr/crdb), the cosmic-ray database for charged species. CRDB is based on MySQL, queried and sorted by jquery and tsorter libraries, and displayed via PHP web pages through AJAX protocol. We review the modifications made on the structure and outputs of the database since the first release (Maurin et al., 2014). For this update, the most important feature is the inclusion of ultra-heavy nuclei ($Z>30$), ultra-high energy nuclei (from $10^{15}$ to $10^{20}$ eV), and limits on anti-nuclei fluxes ($Z\leq -1$ for $A>1$); more than 100 experiments, 350 publications, and 40000 data points are now available in CRDB. We also revisited and simplified how users can retrieve data and submit new data. For questions and request, please contact crdb@lpsc.in2p3.fr.

The vertical diffusive halo size of the Galaxy, $L$, is a key parameter for dark matter indirect searches. It can be better determined thanks to recent AMS-02 data. We set constraints on $L$ from Be/B and $^{10}$Be/Be data, and we performed a consistency check with positron data. We detail the dependence of Be/B and $^{10}$Be/Be on $L$ and forecast on which energy range better data would be helpful for future $L$ improvements. We used USINE v3.5 for the propagation of nuclei, and $e^+$ were calculated with the pinching method of Boudaud et al. (2017). The current AMS-02 Be/B ($\sim3\%$ precision) and ACE-CRIS $^{10}$Be/Be ($\sim 10\%$ precision) data bring similar and consistent constraints on $L$. The AMS-02 Be/B data alone constrain $L=5^{+3}_{-2}$~kpc at a 68\% confidence level (spanning different benchmark transport configurations), a range for which most models do not overproduce positrons. Future experiments need to deliver percent-level accuracy on $^{10}$Be/$^9$Be anywhere below 10 GV to further constrain $L$. Forthcoming AMS-02, HELIX, and PAMELA $^{10}$Be/$^9$Be results will further test and possibly tighten the limits derived here. Elemental ratios involving radioactive species with different lifetimes (e.g. Al/Mg and Cl/Ar) are also awaited to provide complementary and robuster constraints.

The AMS-02 experiment has ushered cosmic-ray physics into precision era. In a companion paper, we designed an improved method to calibrate propagation models on B/C data. Here we provide a robust prediction of the $\bar{p}$ flux, accounting for several sources of uncertainties and their correlations. Combined with a correlation matrix for the $\bar{p}$ data, we show that the latter are consistent with a secondary origin. This paper presents key elements relevant to the dark matter search in this channel, notably by pointing out the inherent difficulties in achieving predictions at the percent-level precision.

The AMS-02 experiment measured several secondary-to-primary ratios enabling a detailed study of Galactic cosmic-ray transport. We constrain previously derived benchmark scenarios (based on AMS-02 B/C data only) using other secondary-to-primary ratios, to test the universality of transport and the presence of a low-rigidity diffusion break. We use the 1D thin disc/thick halo propagation model of USINE and a $\chi^2$ minimisation accounting for a covariance matrix of errors (AMS-02 systematics) and nuisance parameters (cross-sections and solar modulation uncertainties). The combined analysis of AMS-02 Li/C, Be/C, and B/C strengthens the case for a diffusion slope of $\delta=0.50\pm 0.03$ with a low-rigidity break or upturn of the diffusion coefficient at GV rigidities. Our simple model can successfully reproduce all considered data (Li/C, Be/C, B/C, N/O, and 3He/4He), although several issues remain: (i) the quantitative agreement depends on the assumptions made on the not well constrained correlation lengths of AMS-02 data systematics; (ii) combined analyses are very sensitive to production cross sections, and we find post-fit values differing by $\sim5-15\%$ from their most likely values (roughly within currently estimated nuclear uncertainties); (iii) two very distinct regions of the parameter space remain viable, either with reacceleration and convection, or with purely diffusive transport. To take full benefit of combined analyses of AMS-02 data, better nuclear data and a better handle on energy correlations in the data systematic are required. AMS-02 data on heavier species are eagerly awaited to further explore cosmic-ray propagation scenarios.

AMS-02 on the International Space Station has been releasing data of unprecedented accuracy. This poses new challenges for their interpretation. We refine the methodology to get a statistically sound determination of the cosmic-ray propagation parameters. We inspect the numerical precision of the model calculation, nuclear cross-section uncertainties, and energy correlations in data systematic errors. We used the 1D diffusion model in USINE. Our $\chi^2$ analysis includes a covariance matrix of errors for AMS-02 systematics and nuisance parameters to account for cross-section uncertainties. Mock data were used to validate some of our choices. We show that any mis-modelling of nuclear cross-section values or the energy correlation length of the covariance matrix of errors biases the analysis. It also makes good models ($\chi^2_{\rm min}/{\rm dof}\approx1$) appear as excluded ($\chi^2_{\rm min}/{\rm dof}\gg1$). We provide a framework to mitigate these effects (AMS-02 data are interpreted in a companion paper). New production cross-section data and the publication by the AMS-02 collaboration of a covariance matrix of errors for each data set would be an important step towards an unbiased view of cosmic-ray propagation in the Galaxy.

I present the first public releases (v3.4 and v3.5) of the USINE code for cosmic-ray propagation in the Galaxy (https://lpsc.in2p3.fr/usine). It contains several semi-analytical propagation models previously used in the literature (leaky-box model, 2-zone 1D and 2D diffusion models) for the calculation of nuclei ($Z=1-30$), anti-protons, and anti-deuterons. For minimisations, the geometry, transport, and source parameters of all models can be enabled as free parameters, whereas nuisance parameters are enabled on solar modulation levels, cross sections (inelastic and production), and systematics of the CR data. With a single ASCII initialisation file to configure runs, its many displays, and the speed associated to semi-analytical approaches, USINE should be a useful tool for beginners, but also for experts to perform statistical analyses of high-precision cosmic-ray data.

The distribution of dark matter (DM) subhalos in our Galaxy remains disputed, leading to varying $\gamma$-ray and $\nu$ flux predictions from their annihilation or decay. In this work, we study how, in the inner Galaxy, subhalo tidal disruption from the Galactic baryonic potential impacts these signals. Based on state-of-the art modeling of this effect from numerical simulations and semi-analytical results, updated subhalo spatial distributions are derived and included in the CLUMPY code. The latter is used to produce a thousand realizations of the $\gamma$-ray and $\nu$ sky. Compared to predictions based on DM only, we conclude a decrease of the flux of the brightest subhalo by a factor 2 to 7 for annihilating DM and no impact on decaying DM: the discovery prospects or limits subhalos can set on DM candidates are affected by the same factor. This study also provides probability density functions for the distance, mass, and angular distribution of the brightest subhalo, among which the mass may hint at its nature: it is most likely a dwarf spheroidal galaxy in the case of strong tidal effects from the baryonic potential, whereas it is lighter and possibly a dark halo for DM only or less pronounced tidal effects.

This article aims at establishing new benchmark scenarios for Galactic cosmic-ray propagation in the GV-TV rigidity range, based on fits to the AMS-02 B/C data with the USINE v3.5 propagation code. We employ a new fitting procedure, cautiously taking into account data systematic error correlations in different rigidity bins and considering Solar modulation potential and leading nuclear cross-section as nuisance parameters. We delineate specific low, intermediate, and high-rigidity ranges that can be related to both features in the data and peculiar microphysics mechanisms resulting in spectral breaks. We single out a scenario which yields excellent fits to the data and includes all the presumably relevant complexity, the BIG model. This model has two limiting regimes: (i) the SLIM model, a minimal diffusion-only setup, and (ii) the QUAINT model, a convection-reacceleration model where transport is tuned by non-relativistic effects. All models lead to robust predictions in the high-energy regime ($\gtrsim10$GV), i.e. independent of the propagation scenario: at $1\sigma$, the diffusion slope $\delta$ is $[0.43-0.53]$, whereas $K_{10}$, the diffusion coefficient at 10GV, is $[0.26-0.36]$kpc$^2$Myr$^{-1}$; we confirm the robustness of the high-energy break, with a typical value $\Delta_h\sim 0.2$. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak ($\sim 4$GV) which might be related to some effective damping scale in the magnetic turbulence.

We present the third release of the CLUMPY code for calculating $\gamma$-ray and $\nu$ signals from annihilations or decays in dark matter structures. This version includes the mean extragalactic signal with several pre-defined options and keywords related to cosmological parameters, mass functions for the dark matter structures, and $\gamma$-ray absorption up to high redshift. For more flexibility and consistency, dark matter halo masses and concentrations are now defined with respect to a user-defined overdensity $\Delta$. We have also made changes for the user's benefit: distribution and versioning of the code via \git{}, less dependencies and a simplified installation, better handling of options in run command lines, consistent naming of parameters, and a new Sphinx documentation at http://lpsc.in2p3.fr/clumpy/.

We present methods to assess whether gamma-ray excesses towards Milky Way dwarf galaxies can be attributed to astrophysical sources rather than to dark matter annihilation. As a case study we focus on Reticulum II, the dwarf which shows the strongest evidence for a gamma-ray signal in Fermi data. Dark matter models and those with curved energy spectra provide good fits to the data, while a simple power law is ruled out at 97.5% confidence. We compare RetII's spectrum to known classes of gamma-ray sources and find a useful representation in terms of spectral curvature and the energy at which the spectral energy distribution peaks. In this space the blazar classes appear segregated from the confidence region occupied by RetII. Pulsars have similar gamma-ray spectra to RetII but we show that RetII is unlikely to host a pulsar population detectable in gamma rays. Tensions with astrophysical explanations are stronger when analyzing 6.5 years of Pass 7 than with the same amount of Pass 8 data, where the excess is less significant. These methods are applicable to any dwarf galaxy which is a promising dark matter target and shows signs of gamma-ray emission along its line of sight.

The accuracy of the current generation of cosmic-ray (CR) experiments, such as AMS-02, PAMELA, CALET, and ISS-CREAM, is now reaching $\sim$1--3\% in a wide range in energy per nucleon from GeV/n to multi-TeV/n. Their correct interpretation could potentially lead to discoveries of new physics and subtle effects that were unthinkable just a decade ago. However, a major obstacle in doing so is the current uncertainty in the isotopic production cross sections that can be as high as 20--50\% or even larger in some cases. While there is a recently reached consensus in the astrophysics community that new measurements of cross sections are desirable, no attempt to evaluate the importance of particular reaction channels and their required accuracy has been made yet. It is, however, clear that it is a huge work that requires an incremental approach. The goal of this study is to provide the ranking of the isotopic cross sections contributing to the production of the most astrophysically important CR Li, Be, B, C, and N species. In this paper, we (i) rank the reaction channels by their importance for a production of a particular isotope, (ii) provide comparisons plots between the models and data used, and (iii) evaluate a generic beam time necessary to reach a 3\% precision in the production cross-sections pertinent to the AMS-02 experiment. This first roadmap may become a starting point in the planning of new measurement campaigns that could be carried out in several nuclear and/or particle physics facilities around the world. A comprehensive evaluation of other isotopes $Z\leq30$ will be a subject of follow-up studies.

Recent high-energy data from Fermi-LAT on the diffuse gamma-ray background have been used to set among the best constraints on annihilating TeV cold dark matter candidates. In order to assess the robustness of these limits, we revisit and update the calculation of the isotropic extragalactic gamma-ray intensity from dark matter annihilation. The emission from halos with masses $\geq10^{10}\,M_{\odot}$ provides a robust lower bound on the predicted intensity. The intensity including smaller halos whose properties are extrapolated from their higher mass counterparts is typically 5 times higher, and boost from subhalos yields an additional factor ~1.5. We also rank the uncertainties from all ingredients and provide a detailed error budget for them. Overall, our fiducial intensity is a factor 5 lower than the one derived by the Fermi-LAT collaboration in their latest analysis. This indicates that the limits set on extragalactic dark matter annihilations could be relaxed by the same factor. We also calculate the expected intensity for self-interacting dark matter in massive halos and find the emission reduced by a factor 3 compared to the collisionless counterpart. The next release of the CLUMPY code will provide all the tools necessary to reproduce and ease future improvements of this prediction.

Using cosmic-ray boron to carbon ratio (B/C) data recently released by the AMS-02 experiment, we find tantalizing indications ({\it decisive evidence}, in Bayesian terms) in favor of a diffusive origin for the broken power-law spectra found in protons ($p$) and helium nuclei (He). The result is robust with respect to currently estimated uncertainties in the cross sections, and in the presence of a small component of primary boron, expected because of spallation at the acceleration site. Reduced errors at high energy as well as further cosmic ray nuclei data (as absolute spectra of C, N, O, Li, Be) may definitively confirm this scenario.

Top-of-atmosphere (TOA) cosmic-ray (CR) fluxes from satellites and balloon-borne experiments are snapshots of the solar activity imprinted on the interstellar (IS) fluxes. Given a series of snapshots, the unknown IS flux shape and the level of modulation (for each snapshot) can be recovered. We wish (i) to provide the most accurate determination of the IS H and He fluxes from TOA data alone, (ii) to obtain the associated modulation levels (and uncertainties) while fully accounting for the correlations with the IS flux uncertainties, and (iii) to inspect whether the minimal force-field approximation is sufficient to explain all the data at hand. Using H and He TOA measurements, including the recent high-precision AMS, BESS-Polar, and PAMELA data, we performed a non-parametric fit of the IS fluxes $J^{\rm IS}_{\rm H,~He}$ and modulation level $\phi_i$ for each data-taking period. We relied on a Markov chain Monte Carlo (MCMC) engine to extract the probability density function and correlations (hence the credible intervals) of the sought parameters. Although H and He are the most abundant and best measured CR species, several datasets had to be excluded from the analysis because of inconsistencies with other measurements. From the subset of data passing our consistency cut, we provide ready-to-use best-fit and credible intervals for the H and He IS fluxes from MeV/n to PeV/n energy (with a relative precision in the range [2-10\%] at 1$\sigma$). Given the strong correlation between $J^{\rm IS}$ and $\phi_i$ parameters, the uncertainties on $J^{\rm IS}$ translate into $\Delta\phi\approx \pm 30$~MV (at 1$\sigma$) for all experiments. We also find that the presence of $^3$He in He data biases $\phi$ towards higher $\phi$ values by $\sim 30$~MV. The force-field approximation gives an excellent ($\chi^2/$dof$=1.02$) description of the recent high-precision TOA H and He fluxes.

Hierarchical structure formation leads to a clumpy distribution of dark matter in the Milky Way. These clumps are possible targets to search for dark matter annihilation with present and future $\gamma$-ray instruments. Many uncertainties exist on the clump distribution, leading to disputed conclusions about the expected number of detectable clumps and the ensuing limits that can be obtained from non-detection. In this paper, we use the CLUMPY code to simulate thousands of skymaps for several clump distributions. This allows us to statistically assess the typical properties (mass, distance, angular size, luminosity) of the detectable clumps. Varying parameters of the clump distributions allows us to identify the key quantities to which the number of detectable clumps is the most sensitive. Focusing our analysis on two extreme clump configurations, yet consistent with results from numerical simulations, we revisit and compare various calculations made for the Fermi-LAT instrument, in terms of number of dark clumps expected and the angular power spectrum for the Galactic signal. We then focus on the prospects of detecting dark clumps with the future CTA instrument, for which we make a detailed sensitivity analysis using open-source CTA software. Based on a realistic scenario for the foreseen CTA extragalactic survey, and accounting for a post-trial sensitivity in the survey, we show that we obtain competitive and complementary limits to those based on long observation of a single bright dwarf spheroidal galaxy.

The expected gamma-ray flux coming from dark matter annihilation in dwarf spheroidal (dSph) galaxies depends on the so-called `J-factor', the integral of the squared dark matter density along the line-of-sight. We examine the degree to which estimates of J are sensitive to contamination (by foreground Milky Way stars and stellar streams) of the stellar-kinematic samples that are used to infer dark matter densities in `ultrafaint' dSphs. Applying standard kinematic analyses to hundreds of mock data sets that include varying levels of contamination, we find that mis-classified contaminants can cause J-factors to be overestimated by orders of magnitude. Stellar-kinematic data sets for which we obtain such biased estimates tend 1) to include relatively large fractions of stars with ambiguous membership status, and 2) to give estimates for J that are sensitive to specific choices about how to weight and/or to exclude stars with ambiguous status. Comparing publicly-available stellar-kinematic samples for the nearby dSphs Reticulum~II and Segue~I, we find that only the latter displays both of these characteristics. Estimates of Segue~I's J-factor should therefore be regarded with a larger degree of caution when planning and interpreting gamma-ray observations. Moreover, robust interpretations regarding dark matter annihilation in dSph galaxies in general will require explicit examination of how interlopers might affect the inferred dark matter density profile.

The level of solar modulation at different times (related to the solar activity) is a central question of solar and galactic cosmic-ray physics. In the first paper of this series, we have established a correspondence between the uncertainties on ground-based detectors count rates and the parameter $\phi$ (modulation level in the force-field approximation) reconstructed from these count rates. In this second paper, we detail a procedure to obtain a reference $\phi$ time series from neutron monitor data. We show that we can have an unbiased and accurate $\phi$ reconstruction ($\Delta\phi/\phi\simeq 10\%$). We also discuss the potential of Bonner spheres spectrometers and muon detectors to provide $\phi$ time series. Two by-products of this calculation are updated $\phi$ values for the cosmic-ray database and a web interface to retrieve and plot $\phi$ from the 50's to today (\url{http://lpsc.in2p3.fr/crdb}).

We present an update of the CLUMPY code for the calculation of the astrophysical J-factors (from dark matter annihilation/decay) for any Galactic or extragalactic dark matter halo including substructures: halo-to-halo concentration scatter may now be enabled, boost factors can include several levels of substructures, and triaxiality is a new option for dark matter haloes. This new version takes advantage of the CFITSIO and HEALPix libraries to propose FITS output maps using the HEALPix pixelisation scheme. Skymaps for $\gamma$-ray and $\nu$ signals from generic annihilation/decay spectra are now direct outputs of CLUMPY. Making use of HEALPix routines, smoothing by a user-defined instrumental Gaussian beam and computing the angular power spectrum of the maps is now possible. In addition to these improvements, the main novelty is the implementation of a Jeans analysis module, to obtain dark matter density profiles from kinematic data in relaxed spherical systems (e.g., dwarf spheroidal galaxies). The code is also interfaced with the GreAT toolkit designed for Markov Chain Monte Carlo analyses, from which probability density functions and credible intervals can be obtained for velocity dispersions, dark matter profiles, and J-factors.

The dwarf spheroidal galaxies (dSph) of the Milky Way are among the most attractive targets for indirect searches of dark matter. In this work, we reconstruct the dark matter annihilation (J-factor) and decay profiles for the newly discovered dSph Reticulum II. Using an optimized spherical Jeans analysis of kinematic data obtained from the Michigan/Magellan Fiber System (M2FS), we find Reticulum II's J-factor to be among the largest of any Milky Way dSph. We have checked the robustness of this result against several ingredients of the analysis. Unless it suffers from tidal disruption or significant inflation of its velocity dispersion from binary stars, Reticulum II may provide a unique window on dark matter particle properties.

Dwarf spheroidal (dSph) galaxies are prime targets for present and future gamma-ray telescopes hunting for indirect signals of particle dark matter. The interpretation of the data requires careful assessment of their dark matter content in order to derive robust constraints on candidate relic particles. Here, we use an optimised spherical Jeans analysis to reconstruct the `astrophysical factor' for both annihilating and decaying dark matter in 21 known dSphs. Improvements with respect to previous works are: (i) the use of more flexible luminosity and anisotropy profiles to minimise biases, (ii) the use of weak priors tailored on extensive sets of contamination-free mock data to improve the confidence intervals, (iii) systematic cross-checks of binned and unbinned analyses on mock and real data, and (iv) the use of mock data including stellar contamination to test the impact on reconstructed signals. Our analysis provides updated values for the dark matter content of 8 `classical' and 13 `ultrafaint' dSphs, with the quoted uncertainties directly linked to the sample size; the more flexible parametrisation we use results in changes compared to previous calculations. This translates into our ranking of potentially-brightest and most robust targets---viz., Ursa Minor, Draco, Sculptor---, and of the more promising, but uncertain targets---viz., Ursa Major 2, Coma---for annihilating dark matter. Our analysis of Segue 1 is extremely sensitive to whether we include or exclude a few marginal member stars, making this target one of the most uncertain. Our analysis illustrates challenges that will need to be addressed when inferring the dark matter content of new `ultrafaint' satellites that are beginning to be discovered in southern sky surveys.

Galactic Cosmic-ray (CR) transport parameters are usually constrained by the boron-to-carbon ratio. This procedure is generically plagued with degeneracies between the diffusion coefficient and the vertical extent of the Galactic magnetic halo. The latter is of paramount importance for indirect dark matter (DM) searches, because it fixes the amount of DM annihilation or decay that contributes to the local antimatter CR flux. These degeneracies could be broken by using secondary radioactive species, but the current data still have large error bars, and this method is extremely sensitive to the very local interstellar medium (ISM) properties. Here, we propose to use the low-energy CR positrons in the GeV range as another direct constraint on diffusion models. We show that the PAMELA data disfavor small diffusion halo ($L\lesssim 3$ kpc) and large diffusion slope models, and exclude the minimal ({\em min}) configuration (Maurin et al. 2001, Donato et al. 2004) widely used in the literature to bracket the uncertainties in the DM signal predictions. This is complementary to indirect constraints (diffuse radio and gamma-ray emissions) and has strong impact on DM searches. Indeed this makes the antiproton constraints more robust while enhancing the discovery/exclusion potential of current and future experiments, like AMS-02 and GAPS, especially in the antiproton and antideuteron channels.

Dwarf spheroidal (dSph) galaxies are among the most promising targets for the indirect detection of dark matter (DM) from annihilation and/or decay products. Empirical estimates of their DM content - and hence the magnitudes of expected signals - rely on inferences from stellar-kinematic data. However, various kinematic analyses can give different results and it is not obvious which are most reliable. Using extensive sets of mock data of various sizes (mimicking 'ultra-faint' and 'classical' dSphs) and an MCMC engine, here we investigate biases, uncertainties, and limitations of analyses based on parametric solutions to the spherical Jeans equation. For a variety of functional forms for the tracer and DM density profiles, as well as the orbital anisotropy profile, we examine reliability of estimates for the astrophysical J- and D-factors for annihilation and decay, respectively. For large (N > 1000) stellar-kinematic samples typical of 'classical' dSphs, errors tend to be dominated by systematics, which can be reduced through the use of sufficiently general and flexible functional forms. For small (N < 100) samples typical of 'ultrafaints', statistical uncertainties tend to dominate systematic errors and flexible models are less necessary. We define an optimal strategy that would mitigate sensitivity to priors and other aspects of analyses based on the spherical Jeans equation. We also find that the assumption of spherical symmetry can bias estimates of J (with the 95% credibility intervals not encompassing the true J-factor) when the object is mildly triaxial (axis ratios b/a = 0.8, c/a = 0.6). A concluding table summarises the typical error budget and biases for the different sample sizes considered.

Particles count rates at given Earth location and altitude result from the convolution of (i) the interstellar (IS) cosmic-ray fluxes outside the solar cavity, (ii) the time-dependent modulation of IS into Top-of-Atmosphere (TOA) fluxes, (iii) the rigidity cut-off (or geomagnetic transmission function) and grammage at the counter location, (iv) the atmosphere response to incoming TOA cosmic rays (shower development), and (v) the counter response to the various particles/energies in the shower. Count rates from neutron monitors or muon counters are therefore a proxy to solar activity. In this paper, we review all ingredients, discuss how their uncertainties impact count rate calculations, and how they translate into variation/uncertainties on the level of solar modulation $\varphi$ (in the simple Force-Field approximation). The main uncertainty for neutron monitors is related to the yield function. However, many other effects have a significant impact, at the 5-10\% level on $\varphi$ values. We find no clear ranking of the dominant effects, as some depend on the station position and/or the weather and/or the season. An abacus to translate any variation of count rates (for neutron and $\mu$ detectors) to a variation of the solar modulation $\varphi$ is provided.

This paper gives a description of a new on-line database http://lpsc.in2p3.fr/crdb and associated on-line tools (data selection, data export, plots, etc.) for charged cosmic-ray measurements. The experimental setups (type, flight dates, techniques) from which the data originate are included in the database, along with the references to all relevant publications. The database relies on the MySQL5 engine. The web pages and queries are based on PHP, AJAX and the jquery, jquery.cluetip, jquery-ui, and table-sorter third-party libraries. In this first release, we restrict ourselves to Galactic cosmic rays with Z<=30 and a kinetic energy per nucleon up to a few tens of TeV/n. This corresponds to more than 200 different sub-experiments (i.e., different experiments, or data from the same experiment flying at different times) in as many publications. We set up a cosmic-ray database and provide tools to sort and visualise the data. New data can be submitted, providing the community with a collaborative tool to archive past and future cosmic-ray measurements. Any help/ideas to further expand and/or complement the database is welcome (please contact crdatabase@lpsc.in2p3.fr).

The Fornax galaxy cluster was observed with the High Energy Stereoscopic System (H.E.S.S.) for a total live time of 14.5 hours, searching for very-high-energy (VHE, E>100 GeV) gamma-rays from dark matter (DM) annihilation. No significant signal was found in searches for point-like and extended emissions. Using several models of the DM density distribution, upper limits on the DM velocity-weighted annihilation cross-section <sigma v> as a function of the DM particle mass are derived. Constraints are derived for different DM particle models, such as those arising from Kaluza-Klein and supersymmetric models. Various annihilation final states are considered. Possible enhancements of the DM annihilation gamma-ray flux, due to DM substructures of the DM host halo, or from the Sommerfeld effect, are studied. Additional gamma-ray contributions from internal bremsstrahlung and inverse Compton radiation are also discussed. For a DM particle mass of 1 TeV, the exclusion limits at 95% of confidence level reach values of <sigma v> ~ 10^-23cm^3s^-1, depending on the DM particle model and halo properties. Additional contribution from DM substructures can improve the upper limits on <\sigma v> by more than two orders of magnitude. At masses around 4.5 TeV, the enhancement by substructures and the Sommerfeld resonance effect results in a velocity-weighted annihilation cross-section upper limit at the level of <\sigma v> ~ 10^-26cm^3s^-1.

Galaxy clusters are among the best targets for indirect dark matter detection in gamma-rays, despite the large astrophysical background expected from these objects. Detection is now within reach of current observatories (Fermi-LAT or Cerenkov telescopes); however, assessing the origin of this signal might be difficult. We investigate whether the behaviour of the number of objects per `flux' bin (logN-logF) and that of the stacked signal could be used as a signature of the dominant process at stake.We use the CLUMPY code to integrate the signal from decaying or annihilating dark matter and cosmic rays along the line of sight. We assume the standard NFW profile for the dark matter density and rely on a parametrised emissivity for the cosmic-ray component. In this context, the consequences of stacking are explored using the MCXC meta-catalogue of galaxy clusters. We find the value of the slope of the logN-logF power law (or the increase of the signal with the number of stacked objects) to be a clear diagnosis to disentangle decaying dark matter from cosmic-ray induced gamma-rays. For dark matter annihilation, depending on the signal boost from the substructures, it is either similar to the cosmic-ray signal (no boost) or similar to the decay case (large boosts). The shift between the brightest object and its followers also depends on the signal origin. For annihilation, this shift and the stacked signal are poorly constrained because of the large uncertainty affecting the boost. We also underline that the angular dependence of the annihilation signal is not universal because of the substructure contribution.

Clusters of galaxies are potentially important targets for indirect searches for dark matter annihilation. Here we reassess the detection prospects for annihilation in massive halos, based on a statistical investigation of 1743 clusters in the new Meta-Catalog of X-ray Clusters. We derive a new limit for the extra-galactic dark matter annihilation background of at least 20% of that originating from the Galaxy for an integration angle of 0.1 deg. The number of clusters scales as a power law with their brightness, suggesting that stacking may provide a significant improvement over a single target analysis. The mean angle containing 80% of the dark-matter signal for the sample is ~0.15 deg, indicating that instruments with this angular resolution or better would be optimal for a cluster annihilation search based on stacking. A detailed study based on the Fermi-LAT performance and position-dependent background, suggests that stacking may result in a factor ~2 improvement in sensitivity, depending on the source selection criteria. Based on the expected performance of CTA, we find no improvement with stacking, due to the requirement for pointed observations. We note that several potentially important targets: Opiuchius, A2199, A3627 (Norma) and CIZAJ1324.7-5736 may be disfavoured due to a poor contrast with respect to the Galactic dark-matter signal. The use of the homogenised MCXC meta-catalogue provides a robust ranking of the targets, although the absolute value of their signal depends on the exact dark matter substructure content. For conservative assumptions, we find that galaxy clusters (with or without stacking) can probe <sigma v> down to 1e-25-1e-24 cm3/s for dark matter masses in the range 10 GeV-100 GeV. For more favourable substructure configurations, <sigma v>~1e-26 cm3/s may be reached.

The observational coverage with HESS of the Carina region in VHE gamma-rays benefits from deep exposure (40 h) of the neighboring open cluster Westerlund 2. The observations have revealed a new extended region of VHE gamma-ray emission. The new VHE source HESS J1018-589 shows a bright, point-like emission region positionally coincident with SNR G284.3-1.8 and 1FGL J1018.6 - 5856 and a diffuse extension towards the direction of PSR J1016-5857. A soft Gamma=2.7+-0.5 photon index, with a differential flux at 1TeV of N0=(4.2+-1.1)10^-13 TeV^-1 cm^-2 s^-1 is found for the point-like source, whereas the total emission region including the diffuse emission region is well fit by a power-law function with spectral index Gamma=2.9+-0.4 and differential flux at 1TeV of N0=(6.8+-1.6) 10^-13 TeV^-1 cm^-2 s^-1. This H.E.S.S. detection motivated follow-up X-ray observations with the XMM-Newton satellite to investigate the origin of the VHE emission. The analysis of the XMM-Newton data resulted in the discovery of a bright, non-thermal point-like source (XMMU J101855.4-58564) with a photon index of Gamma=1.65+-0.08 in the center of SNRG284.3-1.8, and a thermal, extended emission region coincident with its bright northern filament. The characteristics of this thermal emission are used to estimate the plasma density in the region as n~0.5 cm^-3(2.9kpc/d)^2. The position of XMMUJ101855.4-58564 is compatible with the position reported by the Fermi-LAT collaboration for the binary system 1FGL J1018.6-5856 and the variable Swift XRT source identified with it. The new X-ray data are used alongside archival multi-wavelength data to investigate the relationship between the VHE gamma-ray emission from HESSJ1018-589 and the various potential counterparts in the Carina arm region.

We show that a stacking approach to galaxy clusters can improve current limits on decaying dark matter by a factor $\gtrsim 5-100$, with respect to a single source analysis, for all-sky instruments such as Fermi-LAT. Based on the largest sample of X-ray-selected galaxy clusters available to date (the MCXC meta-catalogue), we provide all the astrophysical information, in particular the astrophysical term for decaying dark matter, required to perform an analysis with current instruments.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

We present the first public code for semi-analytical calculation of the gamma-ray flux astrophysical J-factor from dark matter annihilation/decay in the Galaxy, including dark matter substructures. The core of the code is the calculation of the line of sight integral of the dark matter density squared (for annihilations) or density (for decaying dark matter). The code can be used in three modes: i) to draw skymaps from the Galactic smooth component and/or the substructure contributions, ii) to calculate the flux from a specific halo (that is not the Galactic halo, e.g. dwarf spheroidal galaxies) or iii) to perform simple statistical operations from a list of allowed DM profiles for a given object. Extragalactic contributions and other tracers of DM annihilation (e.g. positrons, antiprotons) will be included in a second release.

Multiwavelength (MWL) observations of the blazar PKS 2155-304 during two weeks in July and August 2006, the period when two exceptional flares at very high energies (VHE, E>= 100 GeV) occurred, provide a detailed picture of the evolution of its emission. The complete data set from this campaign is presented, including observations in VHE gamma-rays (H.E.S.S.), X-rays (RXTE, CHANDRA, SWIFT XRT), optical (SWIFT UVOT, Bronberg, Watcher, ROTSE), and in the radio band (NRT, HartRAO, ATCA). Optical and radio light curves from 2004 to 2008 are compared to the available VHE data from this period, to put the 2006 campaign into the context of the long-term evolution of the source. The X-ray and VHE gamma-ray emission are correlated during the observed high state of the source, but show no direct connection with longer wavelengths. The long-term flux evolution in the optical and radio bands is found to be correlated and shows that the source reaches a high state at long wavelengths after the occurrence of the VHE flares. Spectral hardening is seen in the SWIFT XRT data. The nightly averaged high-energy spectra of the non-flaring nights can be reproduced by a stationary one-zone SSC model, with only small variations in the parameters. The spectral and flux evolution in the high-energy band during the night of the second VHE flare is modelled with multi-zone SSC models, which can provide relatively simple interpretations for the hour time-scale evolution of the high-energy emission, even for such a complex data set. For the first time in this type of source, a clear indication is found for a relation between high activity at high energies and a long-term increase in the low frequency fluxes.

1ES 0414+009 (z = 0.287) is a distant high-frequency-peaked BL Lac object, and has long been considered a likely emitter of very-high energy (VHE, E>100 GeV) gamma-rays due to its high X-ray and radio flux. Observations in the VHE gamma-ray band and across the electromagnetic spectrum can provide insights into the origin of highly energetic particles present in the source and the radiation processes at work. Because of the distance of the source, the gamma-ray spectrum might provide further limits on the level of the Extragalactic Background Light (EBL). We report observations made between October 2005 and December 2009 with H.E.S.S., an array of four imaging atmospheric Cherenkov telescopes. Observations at high energies (HE, 100 MeV - 100 GeV) with the Fermi-LAT instrument in the first 20 months of its operation are also reported. To complete the multi-wavelength picture, archival UV and X-ray observations with the Swift satellite and optical observations with the ATOM telescope are also used. Based on the observations with H.E.S.S., 1ES 0414+009 is detected for the first time in the VHE band. An excess of 224 events is measured, corresponding to a significance of 7.8 sigma. The photon spectrum of the source is well described by a power law, with photon index of 3.45 \pm 0.25stat \pm 0.20syst. The integral flux above 200 GeV is (1.88 \pm 0.20stat \pm 0.38syst) \times10-12 cm-2 s-1. Observations with the Fermi-LAT in the first 20 months of operation show a flux between 200 MeV and 100 GeV of (2.3 \pm 0.2stat) \times 10-9 erg cm-2 s-1, and a spectrum well described by a power-law function with a photon index 1.85 \pm 0.18. Swift/XRT observations show an X-ray flux between 2 and 10 keV of (0.8 - 1) \times 10-11 erg cm-2 s-1, and a steep spectrum (2.2 - 2.3). Combining X-ray with optical-UV data, a fit with a log-parabolic function locates the synchrotron peak around 0.1 keV. ...

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

Results obtained in very-high-energy (VHE; E > 100 GeV) \gamma-ray observations performed with the H.E.S.S. telescope array are used to investigate particle acceleration processes in the vicinity of the young massive stellar cluster Westerlund 1 (Wd 1). Imaging of Cherenkov light from \gamma-ray induced particle cascades in the Earth's atmosphere is used to search for VHE \gamma\ rays from the region around Wd 1. Possible catalogued counterparts are searched for and discussed in terms of morphology and energetics of the H.E.S.S. source. The detection of the degree-scale extended VHE \gamma-ray source HESS J1646-458 is reported based on 45 hours of H.E.S.S. observations performed between 2004 and 2008. The VHE \gamma-ray source is centred on the nominal position of Wd 1 and detected with a total statistical significance of ~20\sigma. The emission region clearly extends beyond the H.E.S.S. point-spread function (PSF). The differential energy spectrum follows a power law in energy with an index of \Gamma=2.19 \pm 0.08_{stat} \pm 0.20_{sys} and a flux normalisation at 1 TeV of \Phi_0 = (9.0 \pm 1.4_{stat} \pm 1.8_{sys}) x 10^{-12} TeV^{-1} cm^{-2} s^{-1}. The integral flux above 0.2 TeV amounts to (5.2 \pm 0.9) x 10^{-11} cm^{-2} s^{-1}. Four objects coincident with HESS J1646-458 are discussed in the search of a counterpart, namely the magnetar CXOU J164710.2-455216, the X-ray binary 4U 1642-45, the pulsar PSR J1648-4611 and the massive stellar cluster Wd 1. In a single-source scenario, Wd 1 is favoured as site of VHE particle acceleration. Here, a hadronic parent population would be accelerated within the stellar cluster. Beside this, there is evidence for a multi-source origin, where a scenario involving PSR J1648-4611 could be viable to explain parts of the VHE \gamma-ray emission of HESS J1646-458.

The secondary-to-primary B/C ratio is widely used to study Galactic cosmic-ray propagation processes. The 2H/4He and 3He/4He ratios probe a different Z/A regime, therefore testing the `universality' of propagation. We revisit the constraints on diffusion-model parameters set by the quartet (1H, 2H, 3He, 4He), using the most recent data as well as updated formulae for the inelastic and production cross-sections. The analysis relies on the USINE propagation package and a Markov Chain Monte Carlo technique to estimate the probability density functions of the parameters. Simulated data are also used to validate analysis strategies. The fragmentation of CNO cosmic rays (resp. NeMgSiFe) on the ISM during their propagation contributes to 20% (resp. 20%) of the 2H and 15% (resp. 10%) of the 3He flux at high energy. The C to Fe elements are also responsible for up to 10% of the 4He flux measured at 1 GeV/n. The analysis of 3He/4He (and to a less extent 2H/4He) data shows that the transport parameters are consistent with those from the B/C analysis: the diffusion model with delta~0.7 (diffusion slope), Vc~20 km/s (galactic wind), Va~40 km/s (reacceleration) is favoured, but the combination delta~0.2, Vc~0, and Va~80 km/s is a close second. The confidence intervals on the parameters show that the constraints set by the quartet data are competitive with those brought by the B/C data. These constraints are tighter when adding the 3He (or 2H) flux measurements, and the tightest when further adding the He flux. For the latter, the analysis of simulated and real data show an increased sensitivity to biases. Using secondary-to-primary ratio along with a loose prior on the source parameters is recommended to get the most robust constraints on the transport parameters.

Several newly discovered very-high-energy (VHE; E > 100 GeV) gamma-ray sources in the Galaxy are thought to be associated with energetic pulsars. Among them, middle-aged (> 1E+4 yr) systems exhibit large centre-filled VHE nebulae, offset from the pulsar position, which result from the complex relationship between the pulsar wind and the surrounding medium, and reflect the past evolution of the pulsar. Imaging Atmospheric Cherenkov Telescopes (IACTs) have been successful in revealing extended emission from these sources in the VHE regime. Together with radio and X-ray observations, this observational window allows one to probe the energetics and magnetic field inside these large-scale nebulae. H.E.S.S., with its large field of view, angular resolution of < 0.1deg and unprecedented sensitivity, has been used to discover a large population of such VHE sources. In this paper, the H.E.S.S. data from the continuation of the Galactic Plane Survey (-80deg < l < 60deg, |b| < 3deg), together with the existing multi-wavelength observations, are used. A new VHE gamma-ray source was discovered at R.A. (J2000) = 13h56m00s, Dec. (J2000) = -64d30m00s with a 2' statistical error in each coordinate, namely HESS J1356-645. The source is extended, with an intrinsic Gaussian width of (0.20 +/- 0.02)deg. Its integrated energy flux between 1 and 10 TeV of 8E-12 erg cm-2 s-1 represents ~ 11% of the Crab Nebula flux in the same energy band. The energy spectrum between 1 and 20 TeV is well described by a power law dN/dE ~ E-Gamma with photon index Gamma = 2.2 +/- 0.2stat +/- 0.2sys. The inspection of archival radio images at three frequencies and the analysis of X-ray data from ROSAT/PSPC and XMM-Newton/MOS reveal the presence of faint non-thermal diffuse emission coincident with HESS J1356-645. HESS J1356-645 is most likely associated with the young and energetic pulsar PSR J1357-6429 (Abridged)

Due to their large dynamical mass-to-light ratios, dwarf spheroidal galaxies (dSphs) are promising targets for the indirect detection of dark matter (DM) in gamma-rays. We examine their detectability by present and future gamma-ray observatories. The key innovative features of our analysis are: (i) We take into account the angular size of the dSphs; while nearby objects have higher gamma ray flux, their larger angular extent can make them less attractive targets for background-dominated instruments. (ii) We derive DM profiles and the astrophysical J-factor (which parameterises the expected gamma-ray flux, independently of the choice of DM particle model) for the classical dSphs directly from photometric and kinematic data. We assume very little about the DM profile, modelling this as a smooth split-power law distribution, with and without sub-clumps. (iii) We use a Markov Chain Monte Carlo (MCMC) technique to marginalise over unknown parameters and determine the sensitivity of our derived J-factors to both model and measurement uncertainties. (iv) We use simulated DM profiles to demonstrate that our J-factor determinations recover the correct solution within our quoted uncertainties. Our key findings are: (i) Sub-clumps in the dSphs do not usefully boost the signal; (ii) The sensitivity of atmospheric Cherenkov telescopes to dSphs within 20 kpc with cored halos can be up to ~50 times worse than when estimated assuming them to be point-like. Even for the satellite-borne Fermi-LAT the sensitivity is significantly degraded on the relevant angular scales for long exposures, hence it is vital to consider the angular extent of the dSphs when selecting targets; (iii) No DM profile has been ruled out by current data, but using a prior on the inner dark matter cusp slope 0<=gamma<=1 provides J-factor estimates accurate to a factor of a few if an appropriate angular scale [abridged]

The H.E.S.S. very-high-energy (VHE, E > 0.1 TeV) gamma-ray telescope system has discovered a new source, HESS J1747-248. The measured integral flux is (1.2 +/- 0.3) \times 10^-12 cm-2 s-1 above 440 GeV for a power-law photon spectral index of 2.5 +/- 0.3 stat +/- 0.2 sys. The VHE gamma-ray source is located in the close vicinity of the Galactic globular cluster Terzan 5 and extends beyond the H.E.S.S. point spread function (0.07 degree). The probability of a chance coincidence with Terzan 5 and an unrelated VHE source is quite low (~ 10^-4). With the largest population of identified millisecond pulsars (msPSRs), a very high core stellar density and the brightest GeV range flux as measured by Fermi-LAT, Terzan 5 stands out among Galactic globular clusters. The properties of the VHE source are briefly discussed in the context of potential emission mechanisms, notably in relation to msPSRs. Interpretation of the available data accommodates several possible origins for this VHE gamma-ray source, although none of them offers a satisfying explanation of its peculiar morphology.

The recent discovery of the radio shell-type supernova remnant (SNR), G353.6-0.7, in spatial coincidence with the unidentified TeV source HESS J1731-347 has motivated further observations of the source with the High Energy Stereoscopic System (H.E.S.S.) Cherenkov telescope array to test a possible association of the gamma-ray emission with the SNR. With a total of 59 hours of observation, representing about four times the initial exposure available in the discovery paper of HESS J1731-347, the gamma-ray morphology is investigated and compared with the radio morphology. An estimate of the distance is derived by comparing the interstellar absorption derived from X-rays and the one obtained from 12CO and HI observations. The deeper gamma-ray observation of the source has revealed a large shell-type structure with similar position and extension (r~0.25{\deg}) as the radio SNR, thus confirming their association. By accounting for the H.E.S.S. angular resolution and projection effects within a simple shell model, the radial profile is compatible with a thin, spatially unresolved, rim. Together with RX J1713.7-3946, RX J0852.0-4622 and SN 1006, HESS J1731-347 is now the fourth SNR with a significant shell morphology at TeV energies. The derived lower limit on the distance of the SNR of 3.2 kpc is used together with radio and X-ray data to discuss the possible origin of the gamma-ray emission, either via inverse Compton scattering of electrons or the decay of neutral pions resulting from proton-proton interaction.

Observations of the globular clusters NGC 6388 and M 15 were carried out by the H.E.S.S. array of Cherenkov telescopes for a live time of 27.2 and 15.2 hours respectively. No gamma-ray signal is found at the nominal target position of NGC 6388 and M 15. In the primordial formation scenario, globular clusters are formed in a dark matter halo and dark matter could still be present in the baryon-dominated environment of globular clusters. This opens the possibility of observing a dark matter self-annihilation signal. The dark matter content of the globular clusters NGC 6388 and M 15 is modelled taking into account the astrophysical processes that can be expected to influence the dark matter distribution during the evolution of the globular cluster: the adiabatic contraction of dark matter by baryons, the adiabatic growth of a black hole in the dark matter halo and the kinetic heating of dark matter by stars. 95% confidence level exclusion limits on the dark matter particle velocity-weighted annihilation cross section are derived for these dark matter haloes. In the TeV range, the limits on the velocity-weighted annihilation cross section are derived at the 10-25 cm3 s-1 level and a few 10-24 cm3 s-1 for NGC 6388 and M 15 respectively.

We present a new analysis of the relative detectability of dark matter annihilation in the Milky Way's eight classical dwarf spheroidal satellite galaxies. Ours is similar to previous analyses in that we use Markov-Chain Monte Carlo techniques to fit dark matter halo parameters to empirical velocity dispersion profiles via the spherical Jeans equation, but more general in the sense that we do not adopt priors derived from cosmological simulations. We show that even without strong constraints on the shapes of dSph dark matter density profiles (we require only that the inner profile satisfies -lim(r->0) [dlnrho/dln r] <=1), we obtain a robust and accurate constraint on the astrophysical component of a prospective dark matter annihilation signal, provided that the integration angle is approximately twice the projected half-light radius of the dSph divided by distance to the observer, alpha_int~2rh/d. Using this integration angle, which represents a compromise between maximizing prospective flux and minimizing uncertainty in the dSph's dark matter distribution, we calculate the relative detectability of the classical dSphs by ground- and space-based gamma-ray observatories.

A search for a very-high-energy (VHE; >= 100 GeV) gamma-ray signal from self-annihilating particle Dark Matter (DM) is performed towards a region of projected distance r ~ 45-150 pc from the Galactic Center. The background-subtracted gamma-ray spectrum measured with the High Energy Stereoscopic System (H.E.S.S.) gamma-ray instrument in the energy range between 300 GeV and 30 TeV shows no hint of a residual gamma-ray flux. Assuming conventional Navarro-Frenk-White (NFW) and Einasto density profiles, limits are derived on the velocity-weighted annihilation cross section < \sigma v> as a function of the DM particle mass. These are among the best reported so far for this energy range. In particular, for the DM particle mass of ~1 TeV, values for <\sigma v> above 3 * 10^(-25) cm^3 s^(-1) are excluded for the Einasto density profile. The limits derived here differ much less for the chosen density profile parametrizations, as opposed to limits from gamma-ray observations of dwarf galaxies or the very center of the Milky Way, where the discrepancy is significantly larger.

We report on a newly detected point-like source, HESS J1943+213 located in the Galactic plane. This source coincides with an unidentified hard X-ray source IGR J19443+2117, which was proposed to have radio and infrared counterparts. HESS J1943+213 is detected at the significance level of 7.9 \sigma (post-trials) at RA(J2000)=19h 43m 55s +- 1s (stat) +- 1s (sys), DEC(J2000) = +21deg 18' 8" +- 17" (stat) +- 20" (sys). The source has a soft spectrum with photon index Gamma = 3.1 +- 0.3 (stat) +- 0.2 (sys) and a flux above 470 GeV of 1.3 +- 0.2 (stat) +- 0.3 (sys) x 10^{-12} cm^{-2} s^{-1}. There is no Fermi/LAT counterpart down to a flux limit of 6 x 10^{-9} cm^{-2} s^{-1} in the 0.1-100 GeV energy range (95% confidence upper limit calculated for an assumed power-law model with a photon index Gamma=2.0). The data from radio to VHE gamma-rays do not show any significant variability. We combine new H.E.S.S., Fermi/LAT and Nancay Radio Telescope observations with pre-existing non-simultaneous multi-wavelength observations of IGR J19443+2117 and discuss the likely source associations as well as the interpretation as an active galactic nucleus, a gamma-ray binary or a pulsar wind nebula. The lack of a massive stellar counterpart disfavors the binary hypothesis, while the soft VHE spectrum would be very unusual in case of a pulsar wind nebula. In addition, the distance estimates for Galactic counterparts places them outside of the Milky Way. All available observations favor an interpretation as an extreme, high-frequency peaked BL Lac object with a redshift z>0.14. This would be the first time a blazar is detected serendipitously from ground-based VHE observations, and the first VHE AGN detected in the Galactic Plane.

The gamma-ray pulsar PSR B1706-44 and the adjacent supernova remnant (SNR) candidate G343.1-2.3 were observed by H.E.S.S. during a dedicated observation campaign in 2007. As a result of this observation campaign, a new source of very-high-energy (VHE; E > 100 GeV) gamma-ray emission, HESS J1708-443, was detected with a statistical significance of 7 sigma, although no significant point-like emission was detected at the position of the energetic pulsar itself. In this paper, the morphological and spectral analyses of the newly-discovered TeV source are presented. The centroid of HESS J1708-443 is considerably offset from the pulsar and located near the apparent center of the SNR, at RA(J2000) = 17h08m11s +/- 17s and Dec(J2000) = -44d20' +/- 4'. The source is found to be significantly more extended than the H.E.S.S. point spread function (~0.1 deg), with an intrinsic Gaussian width of 0.29 deg +/- 0.04 deg. Its integral flux between 1 and 10 TeV is ~ 3.8 x 10^-12 ph cm^-2 s^-1, equivalent to 17% of the Crab Nebula flux in the same energy range. The measured energy spectrum is well-fit by a power law with a relatively hard photon index Gamma = 2.0 +/- 0.1 (stat) +/- 0.2 (sys). Additional multi-wavelength data, including 330 MHz VLA observations, were used to investigate the VHE gamma-ray source's possible associations with the pulsar wind nebula of PSR B1706-44 and/or with the complex radio structure of the partial shell-type SNR G343.1-2.3.

Several models of Quantum Gravity predict Lorentz Symmetry breaking at energy scales approaching the Planck scale (10^{19} GeV). With present photon data from the observations of distant astrophysical sources, it is possible to constrain the Lorentz Symmetry breaking linear term in the standard photon dispersion relations. Gamma-ray Bursts (GRB) and flaring Active Galactic Nuclei (AGN) are complementary to each other for this purpose, since they are observed at different distances in different energy ranges and with different levels of variability. Following a previous publication of the High Energy Stereoscopic System (H.E.S.S.) collaboration, a more sensitive event-by-event method consisting of a likelihood fit is applied to PKS 2155-304 flare data of MJD 53944 (July 28, 2006) as used in the previous publication. The previous limit on the linear term is improved by a factor of ~3 up to M^{l}_{QG} > 2.1x10^{18} GeV and is currently the best result obtained with blazars. The sensitivity to the quadratic term is lower and provides a limit of M^{q}_{QG} > 6.4x10^10 GeV, which is the best value obtained so far with an AGN and similar to the best limits obtained with GRB.

The Sculptor and Carina Dwarf spheroidal galaxies were observed with the H.E.S.S. Cherenkov telescope array between January 2008 and December 2009. The data sets consist of a total of 11.8 and 14.8 hours of high quality data, respectively. No gamma-ray signal was detected at the nominal positions of these galaxies above 220 GeV and 320 GeV, respectively. Upper limits on the gamma-ray fluxes at 95% C.L. assuming two forms for the spectral energy distribution (a power law shape and one derived from dark matter annihilation) are obtained at the level of 10^-13 to 10^-12 cm^-2s^-1 in the TeV range. Constraints on the velocity weighted dark matter particle annihilation cross section for both Sculptor and Carina dwarf galaxies range from <\sigma v> ~ 10^-21 cm^3s^-1 down to <\sigma v> ~ 10^-22 cm^3s^-1 depending on the dark matter halo model used. Possible enhancements of the gamma-ray flux are studied: the Sommerfeld effect, which is found to exclude some dark matter particle masses, the internal Bremsstrahlung and clumps in the dark-matter halo distributions.

The propagated fluxes of proton, helium, and heavier primary cosmic-ray species (up to Fe) are a means to indirectly access the source spectrum of cosmic rays. We check the compatibility of the primary fluxes with the transport parameters derived from the B/C analysis, but also if they bring further constraints. Proton data are well described in the simplest model defined by a power-law source spectrum and plain diffusion. They can also be accommodated by models with, e.g., convection and/or reacceleration. There is no need for breaks in the source spectral indices below $\sim 1$ TeV/n. Fits on the primary fluxes alone do not provide physical constraints on the transport parameters. If we let free the source spectrum $dQ/dE = q \beta^{\eta_S} {\cal R}^{-\alpha}$ and fix the diffusion coefficient $K(R)= K_0\beta^{\eta_T} {\cal R}^{\delta}$ such as to reproduce the B/C ratio, the MCMC analysis constrains the source spectral index $\alpha$ to be in the range $2.2-2.5$ for all primary species up to Fe, regardless of the value of the diffusion slope $\delta$. The $\eta_S$ low-energy shape of the source spectrum is degenerate with the low-energy shape $\eta_T$ of the diffusion coefficient: we find $\eta_S-\eta_T\approx 0$ for p and He data, but $\eta_S-\eta_T\approx 1$ for C to Fe primary species. This is consistent with the toy-model calculation in which the shape of the p/He and C/O to Fe/O data is reproduced if $\eta_S-\eta_T\approx 0-1$ (no need for different slopes $\alpha$). When plotted as a function of the kinetic energy per nucleon, the low-energy p/He ratio is shaped mostly by the modulation effect, whereas primary/O ratios are mostly shaped by their destruction rate.

Context: The detection of gamma-rays in the very-high-energy (VHE) range (100 GeV-100 TeV) offers the possibility of studying the parent population of ultrarelativistic particles found in astrophysical sources, so it is useful for understanding the underlying astrophysical processes in nonthermal sources. Aim: The discovery of the VHE gamma-ray source HESS J1507-622 is reported and possibilities regarding its nature are investigated. Methods: The H.E.S.S. array of imaging atmospheric Cherenkov telescopes (IACTs) has a high sensitivity compared with previous instruments (~1% of the Crab flux in 25 hours observation time for a 5 sigma point-source detection) and has a large field of view (~5 deg in diameter). HESS J1507-622 was discovered within the ongoing H.E.S.S. survey of the inner Galaxy, and the source was also studied by means of dedicated multiwavelength observations. Results: A Galactic gamma-ray source, HESS J1507-622, located ~3.5 deg from the Galactic plane was detected with a statistical significance > 9 sigma. Its energy spectrum is well fitted by a power law with spectral index \Gamma = 2.24 +/- 0.16_{stat} +/- 0.20_{sys} and a flux above 1 TeV of (1.5 +/- 0.4_{stat} +/- 0.3_{sys}) X 10^{-12} cm^{-2} s^{-1}. Possible interpretations (considering both hadronic and leptonic models) of the VHE gamma-ray emission are discussed in the absence of an obvious counterpart.

Aims. Previous observations with the H.E.S.S. telescope array revealed the existence of extended very-high-energy (VHE; E>100 GeV) {\gamma}-ray emission, HESS J1023-575, coincident with the young stellar cluster Westerlund 2. At the time of discovery, the origin of the observed emission was not unambiguously identified, and follow-up observations have been performed to further investigate the nature of this {\gamma}-ray source. Methods. The Carina region towards the open cluster Westerlund 2 has been re-observed, increasing the total exposure to 45.9 h. The combined dataset includes 33 h of new data and now permits a search for energy-dependent morphology and detailed spectroscopy. Results. A new, hard spectrum VHE {\gamma}-ray source, HESSJ1026-582, was discovered with a statistical significance of 7{\sigma}. It is positionally coincident with the Fermi LAT pulsar PSR J1028-5819. The positional coincidence and radio/{\gamma}-ray characteristics of the LAT pulsar favors a scenario where the TeV emission originates from a pulsar wind nebula. The nature of HESS J1023-575 is discussed in light of the deep H.E.S.S. observations and recent multi-wavelength discoveries, including the Fermi LAT pulsar PSRJ1022-5746 and giant molecular clouds in the region. Despite the improved VHE dataset, a clear identification of the object responsible for the VHE emission from HESS J1023-575 is not yet possible, and contribution from the nearby high-energy pulsar and/or the open cluster remains a possibility.

On-going measurements of the cosmic radiation (nuclear, electronic, and gamma-ray) are shedding new light on cosmic-ray physics. A comprehensive picture of these data relies on an accurate determination of the transport and source parameters of propagation models. A Markov Chain Monte Carlo is used to obtain these parameters in a diffusion model. From the measurement of the B/C ratio and radioactive cosmic-ray clocks, we calculate their probability density functions, with a special emphasis on the halo size L of the Galaxy and the local underdense bubble of size r_h. The analysis relies on the USINE code for propagation and on a Markov Chain Monte Carlo technique (Putze et al. 2009, paper I of this series) for the parameter determination. As found in previous studies, the B/C best-fit model favours diffusion/convection/reacceleration (Model III) over diffusion/reacceleration (Model II). A combined fit on B/C and the isotopic ratios (10Be/9Be, 26Al/27Al, 36Cl/Cl) leads to L ~ 8 kpc and r_h ~ 120 pc for the best-fit Model III. This value for r_h is consistent with direct measurements of the local interstallar medium. For Model II, L ~ 4 kpc and r_h is consistent with zero. We showed the potential and usefulness of the Markov Chain Monte Carlo technique in the analysis of cosmic-ray measurements in diffusion models. The size of the diffusive halo depends crucially on the value of the diffusion slope delta, and also on the presence/absence of the local underdensity damping effect on radioactive nuclei. More precise data from on-going experiments are expected to clarify this issue.

The B/C ratio is used in cosmic-ray physics to constrain the transport parameters. However, from the same set of data, the various published values show a puzzling large scatter of these parameters. We investigate the impact of using different inputs (gas density and hydrogen fraction in the Galactic disc, source spectral shape, low-energy dependence of the diffusion coefficient, and nuclear fragmentation cross-sections) on the best-fit values of the transport parameters. We quantify the systematics produced when varying these inputs, and compare them to statistical uncertainties. We discuss the consequences for the slope of the diffusion coefficient delta. The analysis relies on the propagation code USINE interfaced with the Minuit minimisation routines. We find the typical systematic uncertainties to be larger than the statistical ones. The several published values of delta (from 0.3 to 0.8) can be recovered when varying the low-energy shape of the diffusion coefficient and the convective wind strength. Models including a convective wind are characterised by delta > 0.6, which cannot be reconcile with the expected theoretical values (1/3 and 1/2). However, from a statistical point of view (chi^2 analysis), models with both reacceleration and convection-hence large delta-are favoured. The next favoured models in line yield delta that can be accommodated with 1/3 and 1/2, but require a strong upturn of the diffusion coefficient at low energy (and no convection). To date, using the best statistical tools, the transport parameter determination is still plagued by many unknowns at low energy (~ GeV/n). To disentangle between all these configurations, measurements of the B/C ratio at TeV/n energies and/or combination with other secondary-to-primary ratios is necessary.

PKS 0548-322 (z=0.069) is a ``high-frequency-peaked'' BL Lac object and a candidate very high energy (VHE, E>100 GeV) gamma-ray emitter, due to its high X-ray and radio flux. Observations at the VHE band provide insights into the origin of very energetic particles present in this source and the radiation processes at work. We report observations made between October 2004 and January 2008 with the H.E.S.S. array, a four imaging atmospheric-Cherenkov telescopes. Contemporaneous UV and X-ray observations with the Swift satellite in November 2006 are also reported. PKS 0548-322 is detected for the first time in the VHE band with H.E.S.S. We measure an excess of 216 gamma-rays corresponding to a significance of 5.6 standard deviations. The photon spectrum of the source is described by a power-law, with a photon index of Gamma=2.86 +/- 0.34 (stat) +/- 0.10 (sys). The integral flux above 200 GeV is 1.3 % of the flux of the Crab Nebula, and is consistent with being constant in time. Contemporaneous Swift/XRT observations reveal an X-ray flux between 2 and 10 keV of F_{2-10 keV}=2.3 +/- 0.2 x 10^{-11} erg.cm^{-2}. s^{-1}, an intermediate intensity state with respect to previous observations. The spectral energy distribution can be reproduced using a simple one-zone synchrotron self Compton model, with parameters similar those observed for other sources of this type.

Observations of very high energy gamma-rays from blazars provide information about acceleration mechanisms occurring in their innermost regions. Studies of variability in these objects allow a better understanding of the mechanisms at play. To investigate the spectral and temporal variability of VHE (>100 GeV) gamma-rays of the well-known high-frequency-peaked BL Lac object PKS 2155-304 with the H.E.S.S. imaging atmospheric Cherenkov telescopes over a wide range of flux states. Data collected from 2005 to 2007 are analyzed. Spectra are derived on time scales ranging from 3 years to 4 minutes. Light curve variability is studied through doubling timescales and structure functions, and is compared with red noise process simulations. The source is found to be in a low state from 2005 to 2007, except for a set of exceptional flares which occurred in July 2006. The quiescent state of the source is characterized by an associated mean flux level of 4.32 +/-0.09 x 10^-11 cm^-2 s^-1 above 200 GeV, or approximately 15% of the Crab Nebula, and a power law photon index of 3.53 +/-0.06. During the flares of July 2006, doubling timescales of ~2 min are found. The spectral index variation is examined over two orders of magnitude in flux, yielding different behaviour at low and high fluxes,which is a new phenomenon in VHE gamma-ray emitting blazars. The variability amplitude characterized by the fractional r.m.s. is strongly energy-dependent and is proportional to E^(0.19 +/- 0.01). The light curve r.m.s. correlates with the flux. This is the signature of a multiplicative process which can be accounted for as a red noise with a Fourier index of ~2. This unique data set shows evidence for a low level gamma-ray emission state from PKS 2155-304, which possibly has a different origin than the outbursts. The discovery of the light curve lognormal behaviour might be an indicator ..

AIMS: The properties of the broad-band emission from the high-frequency peaked BL Lac H 2356-309 (z=0.165) are investigated. METHODS: Very High Energy (VHE; E > 100 GeV) observations of H 2356-309 were performed with the High Energy Stereoscopic System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM optical monitoring observations were also made in 2007. RESULTS: A strong VHE signal, ~13 sigma total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The integral flux above 240 GeV is I(>240 GeV) = (3.06 +- 0.26 {stat} +- 0.61 {syst}) x 10^{-12} cm^{-2} s^{-1}, corresponding to ~1.6% of the flux observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a power law (photon index of Gamma = 3.06 +- 0.15 {stat} +- 0.10 {syst}). Significant small-amplitude variations in the VHE flux from H 2356-309 are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by a hard, broken-power-law spectrum on both nights. CONCLUSIONS: The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed X-ray flux is at a historically low level.

Recent theoretical predictions of the lowest very high energy (VHE) luminosity of SN 1006 are only a factor 5 below the previously published H.E.S.S. upper limit, thus motivating further in-depth observations of this source. Deep observations at VHE energies (above 100 GeV) were carried out with the High Energy Stereoscopic System (H.E.S.S.) of Cherenkov Telescopes from 2003 to 2008. More than 100 hours of data have been collected and subjected to an improved analysis procedure. Observations resulted in the detection of VHE gamma-rays from SN 1006. The measured gamma-ray spectrum is compatible with a power-law, the flux is of the order of 1% of that detected from the Crab Nebula, and is thus consistent with the previously established H.E.S.S. upper limit. The source exhibits a bipolar morphology, which is strongly correlated with non-thermal X-rays. Because the thickness of the VHE-shell is compatible with emission from a thin rim, particle acceleration in shock waves is likely to be the origin of the gamma-ray signal. The measured flux level can be accounted for by inverse Compton emission, but a mixed scenario that includes leptonic and hadronic components and takes into account the ambient matter density inferred from observations also leads to a satisfactory description of the multi-wavelength spectrum.

New observations from the Hubble ultra deep field suggest that the star formation rate at z>7 drops off faster than previously thought. Using a newly determined star formation rate for the normal mode of Population II/I stars (PopII/I), including this new constraint, we compute the Thomson scattering optical depth and find a result that is marginally consistent with WMAP5 results. We also reconsider the role of Population III stars (PopIII) in light of cosmological and stellar evolution constraints. While this input may be needed for reionization, we show that it is essential in order to account for cosmic chemical evolution in the early Universe. We investigate the consequences of PopIII stars on the local metallicity distribution function of the Galactic halo (from the recent Hamburg/ESO survey of metal-poor stars) and on the evolution of abundances with metallicity (based on the ESO large program on very metal-poor stars), with special emphasis on carbon-enhanced metal-poor stars. Our most important results show that the nucleosynthetic yields of PopIII stars lead to abundance patterns in agreement with those observed in extremely metal-poor stars. In this chemical approach to cosmic evolution, PopIII stars prove to be a compulsory ingredient, and extremely metal-poor stars are inevitably born at high redshift. (Abridged)

A new calculation of the $\bar{p}/p$ ratio in cosmic rays is compared to the recent PAMELA data. The good match up to 100 GeV allows to set constraints on exotic contributions from thermal WIMP dark matter candidates. We derive stringent limits on possible enhancements of the WIMP \pbar flux: a $m_{\rm WIMP}$=100 GeV (1 TeV) signal cannot be increased by more than a factor 6 (40) without overrunning PAMELA data. Annihilation through the $W^+W^-$ channel is also inspected and cross-checked with $e^+/(e^-+e^+)$ data. This scenario is strongly disfavored as it fails to simultaneously reproduce positron and antiproton measurements.

Propagation of charged cosmic-rays in the Galaxy depends on the transport parameters, whose number can be large depending on the propagation model under scrutiny. A standard approach for determining these parameters is a manual scan, leading to an inefficient and incomplete coverage of the parameter space. We implement a Markov Chain Monte Carlo (MCMC), which is well suited to multi-parameter determination. Its specificities (burn-in length, acceptance, and correlation length) are discussed in the phenomenologically well-understood Leaky-Box Model. From a technical point of view, a trial function based on binary-space partitioning is found to be extremely efficient, allowing a simultaneous determination of up to nine parameters, including transport and source parameters, such as slope and abundances. Our best-fit model includes both a low energy cut-off and reacceleration, whose values are consistent with those found in diffusion models. A Kolmogorov spectrum for the diffusion slope (delta=1/3) is excluded. The marginalised probability-density function for delta and alpha (the slope of the source spectra) are delta~0.55-0.60 and alpha~2.14-2.17, depending on the dataset used and the number of free parameters in the fit. All source-spectrum parameters (slope and abundances) are positively correlated among themselves and with the reacceleration strength, but are negatively correlated with the other propagation parameters. A forthcoming study will extend our analysis to more physical diffusion models.

We use a 1D propagation model to study the dependence of the pbar and dbar exotic fluxes on the transport parameters. The simple analytical solutions allow us i) to clarify the origin of the astrophysical uncertainties, and ii) to compare two models used for {\em signal} predictions, namely the constant and the linear Galactic wind models. We also study how these uncertainties should be reduced using forthcoming nuclear cosmic ray data. We confirm that the degeneracy of the transport parameters for a given propagation model leads to very different fluxes for primary antinuclei (~10^2). However, we show that with forthcoming data, these uncertainties could be greatly reduced (~2). As the precision will increase, the astrophysical uncertainty could then be dominated by our ignorance of the correct spatial dependence for some of the transport parameters: for instance, the constant and the linear wind models do not predict the same amount of exotic pbar at low energy.

Antideuterons are among the most promising galactic cosmic ray-related targets for dark matter indirect detection. Currently only upper limits exist on the flux, but the development of new experiments, such as GAPS and AMS-02, provides exciting perspectives for a positive measurement in the near future. In this Paper, we present a novel and updated calculation of both the secondary and primary antideuteron fluxes. We employ a two-zone diffusion model which successfully reproduces cosmic-ray nuclear data and the observed antiproton flux. We review the nuclear and astrophysical uncertainties and provide an up to date secondary (i.e. background) antideuteron flux. The primary (i.e. signal) contribution is calculated for generic WIMPs annihilating in the galactic halo: we explicitly consider and quantify the various sources of uncertainty in the theoretical evaluations. Propagation uncertainties, as is the case of antiprotons, are sizeable. Nevertheless, antideuterons offer an exciting target for indirect dark matter detection for low and intermediate mass WIMP dark matter. We then show the reaching capabilities of the future experiments for neutralino dark matter in a variety of supersymmetric models.

Anti-proton and positron Galactic cosmic ray spectra are among the key targets for indirect detection of dark matter. The boost factors, corresponding to an enhancement of the signal and linked to the clumpiness properties of the dark matter distribution, have been taken as high as thousands in the past. The dramatic impact of these boost factors for indirect detection of antiparticles, for instance with the PAMELA satellite or the coming AMS-02 experiment, asks for their detailed calculation. We take into account the results of high resolution N-body dark matter simulations to calculate the most likely energy dependent boost factors linked to the cosmic ray propagation properties, for anti-protons and positrons. Starting from the mass and space distributions of sub-halos, the anti-proton and positron propagators are used to calculate the mean value and the variance of the boost factor for the primary fluxes. We take advantage of the statistical method introduced in Lavalle et al. (2007) and cross-check the results with Monte Carlo computations. By spanning some extreme configurations of sub-halo and propagation properties, we find that the average contribution of the clumps is negligible compared to that of the smooth dark matter component. Sub-halos do not lead to enhancement of the signals, unless they are taken with some extreme (unexpected) properties. This result is independent of the nature of the self-annihilating dark matter candidate considered, and provides precise estimates of the theoretical and the statistical uncertainties of the antimatter flux from dark matter substructures. Spectral distortions can still be expected in antimatter flux measurements, but scenarios invoking large and even mild clumpiness boost factors are strongly disfavoured by our analysis.

We calculate the antideuteron flux expected from dark matter annihilation in the galactic halo. The propagation is treated in a full 2-D propagation model consistent with the results obtained from the propagation of B/C and other galactic species. We discuss the potentials of this indirect dark matter detection means, with special emphasis on the possible sources of uncertainties affecting future measurements

Recently, new observations of Li6 in Pop II stars of the galactic halo have shown a surprisingly high abundance of this isotope, about a thousand times higher than its predicted primordial value. In previous papers, a cosmological model for the cosmic ray-induced production of this isotope in the IGM has been developed to explain the observed abundance at low metallicity. In this paper, given this constraint on the Li6, we calculate the non-thermal evolution with redshift of D, Be, and B in the IGM. In addition to cosmological cosmic ray interactions in the IGM, we include additional processes driven by SN explosions: neutrino spallation and a low energy component in the structures ejected by outflows to the IGM. We take into account CNO CRs impinging on the intergalactic gas. Although subdominant in the galactic disk, this process is shown to produce the bulk of Be and B in the IGM, due to the differential metal enrichment between structures (where CRs originate) and the IGM. We also consider the resulting extragalactic gamma-ray background which we find to be well below existing data. The computation is performed in the framework of hierarchical structure formation considering several star formation histories including Pop III stars. We find that D production is negligible and that a potentially detectable Be and B plateau is produced by these processes at the time of the formation of the Galaxy (z ~ 3).

The transport of the nuclear cosmic ray flux in the atmosphere is studied and the atmospheric corrections to be applied to the measurements are calculated. The contribution of the calculated corrections to the accuracy of the experimental results are discussed and evaluated over the kinetic energy range 10-10$^{3}$ GeV/n. The Boron (B) and Carbon (C) elements system is used as a test case. It is shown that the required corrections become largely dominant at the highest energies investigated. The results are discussed.

The Galactic secondary cosmic ray anti-proton flux calculated with different propagation models is fairly consistent with data, and the associated propagation uncertainty is small. This is not the case for any anti-proton exotic component of the dark matter halo. Detailed propagation models are mandatory if the ultimate goal is to explain an excess. However, simpler and faster approximate formulae for anti-protons are an attractive alternative to quickly check that a given dark matter model is not inconsistent with the anti-proton observed flux. This paper provides such formulae. In addition, they could be used to put constraints on new physics in this channel, where an extensive scan of a large parameter space could otherwise be quite expensive in computer ressources.

The Vela region is a promising target for the detection of 1.8 MeV gamma-rays emitted by the decays of radioactive 26Al isotopes produced in hydrostatic or explosive stellar nucleosynthesis processes. COMPTEL has claimed 1.8 MeV gamma-ray detection from Vela at a 3sigma level with a flux of 3.6 10^-5 ph/cm^2/s. In this paper, we present first results of our search for 1.8 MeV gamma-rays from Vela with the spectrometer SPI aboard INTEGRAL. Using the data set acquired during 1.7 Ms at the end of 2005 in the frame of our AO-3 open-time observation, we determine a flux of (6.5 \pm 1.9(stat) \pm 2.4(syst)) 10^-5 ph/cm^2/s from 26Al decays in the Vela region.

Extra dimensions offer new ways to address long-standing problems in beyond the standard model particle physics. In some classes of extra-dimensional models, the lightest Kaluza-Klein particle is a viable dark matter candidate. In this work, we study indirect detection of Kaluza-Klein dark matter via its annihilation into antiprotons. We use a sophisticated galactic cosmic ray diffusion model whose parameters are fully constrained by an extensive set of experimental data. We discuss how fluxes of cosmic antiprotons can be used to exclude low Kaluza-Klein masses.

The fluxes of light antinuclei A<=4 induced near earth by Cosmic Ray intera ctions with the interstellar matter in the Galaxy and with the Earth atmosphere, are calculated in a phenomenological framework. The hadronic production cross section for antinucleons is based on a recent parametrization of a wide set of accelerator data. The production of light nuclei is calculated using coalescence models. For the standard coalescence model, the coalescence radius is fitted to the available experimental data. The non annihilating inelastic scattering process for the antideuterons is discussed and taken into account for the first time via a more realistic procedure than used so far for antiprotons.

Measurements of ultra heavy nuclei at GeV/n energies in the galactic cosmic radiation address the question of the sources (nucleosynthetic s- and r-processes). As such, the determination of CR source abundances is a promising way to discriminate between existing nucleosynthesis models. For primary species (nuclei present and accelerated at sources), it is generally assumed that the relative propagated abundances, if they are close in mass, are not too different from their relative source abundances. Besides, the range of the correction factor associated to propagation has been estimated in weighted slab models only. Heavy CRs that are detected near Earth were accelerated from regions that are closer to us than were the light nuclei. Hence, the geometry of sources in the Solar neighbourhood, and as equally important, the geometry of gas in the same region, must be taken into account. In this paper, a two zone diffusion model is used, and as was previously investigated for radioactive species, we report here on the impact of the local interstellar medium (LISM) feature (under-dense medium over a scale ~100 pc) on primary and secondary stable nuclei propagated abundances. Going down to Fe nuclei, the connection between heavy and light abundances is also inspected. A general trend is found that decreases the UHCR source abundances relative to the HCR ones. This could have an impact on the level of r-process required to reproduce the data.

Most studies of GeV Galactic Cosmic Rays (GCR) nuclei assume a steady state/continuous distribution for the sources of cosmic rays, but this distribution is actually discrete in time and in space. The current progress in our understanding of cosmic ray physics (acceleration, propagation), the required consistency in explaining several GCRs manifestation (nuclei, $\gamma$,...) as well as the precision of present and future space missions (e.g. INTEGRAL, AMS, AGILE, GLAST) point towards the necessity to go beyond this approximation. A steady state semi-analytical model that describes well many nuclei data has been developed in the past years based on this approximation, as well as others. We wish to extend it to a time dependent version, including discrete sources. As a first step, the validity of several approximations of the model we use are checked to validate the approach: i) the effect of the radial variation of the interstellar gas density is inspected and ii) the effect of a specific modeling for the galactic wind (linear vs constant) is discussed. In a second step, the approximation of using continuous sources in space is considered. This is completed by a study of time discreteness through the time-dependent version of the propagation equation. A new analytical solution of this equation for instantaneous point-like sources, including the effect of escape, galactic wind and spallation, is presented. Application of time and space discretness to definite propagation conditions and realistic distributions of sources will be presented in a future paper.

We calculate the antiproton flux due to relic neutralino annihilations, in a two-dimensional diffusion model compatible with stable and radioactive cosmic ray nuclei. We find that the uncertainty in the primary flux induced by the propagation parameters alone is about two orders of magnitude at low energies, and it is mainly determined by the lack of knowledge on the thickness of the diffusive halo. On the contrary, different dark matter density profiles do not significantly alter the flux: a NFW distribution produces fluxes which are at most 20% higher than an isothermal sphere. The most conservative choice for propagation parameters and dark matter distribution normalization, together with current data on cosmic antiprotons, cannot lead to any definitive constraint on the supersymmetric parameter space, neither in a low-energy effective MSSM, or in a minimal SUGRA scheme. However, if the best choice for propagation parameters - corresponding to a diffusive halo of L=4 kpc - is adopted, some supersymmetric configurations with the neutralino mass of about 100 GeV should be considered as excluded. An enhancement flux factor - due for instance to a clumpy dark halo or to a higher local dark matter density - would imply a more severe cut on the supersymmetric parameters.

We use the constraints on the diffusion parameters as obtained with stable nuclei to calculate the cosmic antiproton flux from annihilating relic neutralinos. We discuss the relevance of each characteristic parameter, describing our two dimension diffusion model, on the flux of antiprotons produced in the dark halo of our Galaxy. We estimate a two orders of magnitude uncertainty on the flux due to the unknowledge of the propagation parameters. A conservative and systematic evaluation of the flux in the supersymmetric parameter space is done in order to exclude configurations providing a total (secondary plus primary) flux in excess with observations. We also study the effect on the flux induced by modifications in the distribution of cold dark matter in the Galaxy.

We present the results on the source spectrum function for primary nuclei in galactic cosmic rays, where two distinct energy dependences are used for the source spectra. We discuss the evolution of the goodness of fit to B/C data with the propagation parameters and also show that the results are not much affected by a different choice for the diffusion scheme. We apply the constraints on the diffusion scheme as derived from stable nuclei to calculate the propagation of beta-radioactive isotopes. The diffusion model is refined to properly take into account the effect of a local bubble in the interstellar medium. Our calculations are compared to existing data, which prefer the local bubble description instead of the homogeneus scheme.

In a companion paper, we investigated the question of the spatial origin of the cosmic rays detected in the Solar neighborhood, in the case of standard sources located in the Galactic disk. There are some reasons to believe that there may also be a large number of sources located in the halo, for example if the Galactic dark matter is made of supersymmetric particles or if Primordial Black Holes are present. These exotic sources could enhance the $\bar{p}$, $\bar{d}$ or positrons above the standard background, indicating the existence of new physics. The spatial distribution of these hypothetical sources, though an important ingredient to evaluate these exotic signals, is poorly known. The aim of this paper is to point out that this discussion should not be disconnected from that of the propagation properties in the Galaxy. More precisely, we determine the regions of the halo from which a significant fraction f of cosmic rays antiprotons and antideuterons detected in the Solar neighborhood were emitted (we refer to these regions as f-volumes), for different sets of propagation parameters consistent with B/C data, as derived in Maurin et al (2002). It is found that some of them lead to rather small f-volumes, indicating that the exotic cosmic rays could have a local origin (in particular for a small halo or a large Galactic convective wind), coming from the solar neighborhood or the Galactic center region. It is also found that the dark matter density enhancement (spike) due to the accretion around the central supermassive black hole gives a negligible contribution to the exotic charged particle signal on Earth. The case of electrons and positrons is also discussed.

The propagation of Galactic Cosmic Ray nuclei having energies between 100 MeV/nuc and several PeV/nuc is strongly believed to be of diffusive nature. The particles emitted by a source located in the disk do not pervade the whole Galaxy, but are rather confined to a smaller region whose spatial extension is related to the height of the diffusive halo, the Galactic wind and the spallation rate. Following the pioneering work of Jones (1978), this paper presents a general study on the spatial origin of cosmic rays, with a particular attention to the role of spallations and Galactic wind. This question is different, and to a certain extent disconnected, from that of the origin of cosmic rays. We find the regions of the disk from which a given fraction of cosmic rays detected in the Solar neighborhood were emitted (f-surfaces). After a general study, we apply the results to a realistic source distribution, with the propagation parameters obtained in our previous systematic analysis of the observed secondary-to-primary ratios Maurin et al (2002). The shape and size of these f-surfaces depend on the species as well as on the values of the propagation parameters. For some of the models preferred by our previous analysis (i.e. large diffusion slope $\delta$), these f-surfaces are small and in some extreme cases only a fraction of a percent of the whole Galactic sources actually contribute to the Solar neighborhood Cosmic Ray flux. Moreover, a very small number of sources may be responsible for more than 15 % of the flux detected in the Solar neighborhood. This may point towards the necessity to go beyond the approximations of both homogeneity and stationarity. Finally, the observed primary composition is dominated by sources within a few kpc.

In most cosmological models, primordial black holes (PBHs) should have formed in the early Universe. Their Hawking evaporation into particles could eventually lead to the formation of antideuterium nuclei. This paper is devoted to a first computation of this antideuteron flux. The production of these antinuclei is studied with a simple coalescence scheme, and their propagation in the Galaxy is treated with a well-constrained diffusion model. We compare the resulting primary flux to the secondary background, due to the spallation of protons on the interstellar matter. Antideuterons are shown to be a very sensitive probe for primordial black holes in our Galaxy. The next generation of experiments should allow investigators to significantly improve the current upper limit, nor even provide the first evidence of the existence of evaporating black holes.

Cosmic Ray nuclei in the energy range 100 MeV/nuc - 100 GeV/nuc provide crucial information about the physical properties of the Galaxy. They can also be used to answer questions related to astroparticle physics. This paper reviews the results obtained in this direction, with a strong bias towards the work done by the authors at {\sc lapth}, {\sc isn} and {\sc iap}. The propagation of these nuclei is studied quantitatively in the framework of a semi-analytical two-zone diffusion model taking into account the effect of galactic wind, diffuse reacceleration and energy losses. The parameters of this model are severely constrained by an analysis of the observed B/C ratio. These constraints are then used to study other species such as radioactive species and light antinuclei. Finally, we focus on the astroparticle subject and we study the flux of antiprotons and antideuterons that might be due to neutralino annihilations or primordial black hole evaporation. The question of the spatial origin of all these species is also addressed.

In a previous study (Maurin et al., 2001), we explored the set of parameters describing diffusive propagation of cosmic rays (galactic convection, reacceleration, halo thickness, spectral index and normalization of the diffusion coefficient), and we identified those giving a good fit to the measured B/C ratio. This study is now extended to take into account a sixth free parameter, namely the spectral index of sources. We use an updated version of our code where the reacceleration term comes from standard minimal reacceleration models. The goal of this paper is to present a general view of the evolution of the goodness of fit to B/C data with the propagation parameters. In particular, we find that, unlike the well accepted picture, and in accordance with our previous study, a Kolmogorov-like power spectrum for diffusion is strongly disfavored. Rather, the $\chi^2$ analysis points towards $\delta\gtrsim 0.7$ along with source spectra index $\lesssim 2.0$. Two distinct energy dependences are used for the source spectra: the usual power-law in rigidity and a law modified at low energy, the second choice being only slightly preferred. We also show that the results are not much affected by a different choice for the diffusion scheme. Finally, we compare our findings to recent works, using other propagation models. This study will be further refined in a companion paper, focusing on the fluxes of cosmic ray nuclei.

A substantial increase of the mean logarithmic mass $<\ln A>$ of galactic cosmic rays {\em vs} energy has been observed . We study three effects that could explain this trend i) different source spectra for protons and heavy nuclei ii) a selective nuclear destruction in flight of heavies iii) a gradient of the source number and chemical composition in the galactic disk. We take advantage of the diffusive cosmic ray propagation model developed at LAPTH to study specifically the geometrical aspects of the propagation and extend it to high energy. Using a simple modeling of the spectral knee around $10^{15}$ eV, a bump in $<\ln A>$ appears. This feature is smoother when the spectral index of protons is steeper than Fe's. We analyze the effects of the rigidity dependence of the diffusion coefficient and the scale height of the confinement halo and we show that $<\ln A>$ is most sensitive to the first parameter. Pure geometrical effects are less determining than the diffusion coefficient spectral index. Subsequently, we conclude that the physics of cosmic ray confinement is the essential cause of the heavy nuclei enrichment until $\sim 10^{15}$ eV.

Primordial black holes (PBHs) have motivated many studies since it was shown that they should evaporate and produce all kinds of particles. Recent experimental measurements of cosmic rays with great accuracy, theoretical investigations on the possible formation mechanisms and detailed evaporation processes have revived the interest in such astrophysical objects. This article aims at using the latest developments on antiproton propagation models together with new data from BESS, CAPRICE and AMS experiments to constrain the local amount of PBH dark matter. Depending on the diffusion halo parameters and on the details of emission mechanism, we derive an average upper limit of the order of rho < 1.7E-33 g cm^-3.

In the present paper, we extend the analysis of Maurin et al. (2001) and Donato et al. (2001) to the $\beta$-radioactive nuclei $^{10}$Be, $^{26}$Al and $^{36}$Cl. These species are be shown to be particularly sensitive to the properties of the local interstellar medium (LISM). As studies of the LISM suggest that we live in an underdense bubble of extent $r_{hole} \sim 50 - 200 \unit{pc}$, this local feature must be taken into account. We present a modified version of our diffusion model which describes the underdensity as a hole in the galactic disc. It is found that the presence of the bubble leads to a decrease in the radioactive fluxes which can be approximated by a simple factor $\exp(-r_{hole}/l_{rad})$ where $l_{rad}=\sqrt{K \gamma \tau_0}$ is the typical distance travelled by a radioactive nucleus before it decays. We find that each of the radioactive nuclei independently point towards a bubble radius $\lesssim 100 \unit{pc}$. If these nuclei are considered simultaneously, only models with a bubble radius $r_{hole} \sim 60 - 100 \unit{pc}$ are marginally consistent with data. In particular, the standard case $r_{hole}=0 \unit{pc}$ is disfavoured. Our main concern is about the consistency of the currently available data, especially $^{26}$Al/$^{27}$Al.

Cosmic ray antiprotons provide an important probe for the study of the galactic Dark Matter, as they could be produced by exotic sources. On the other hand, antiprotons are anyway produced by standard nuclear reactions of cosmic ray nuclei on interstellar matter. This process is responsible for a background flux that must be carefully determined to estimate the detectability of an hypothetical exotic signal. Estimates of this background suffer from potential uncertainties of various origins. The propagation of cosmic antiprotons depends on several physical characteristics of the Galaxy which are poorly known. Antiprotons are created from cosmic protons and helium nuclei whose fluxes were not measured with great accuracy until very recently. Calculations of antiproton fluxes make use of nuclear physics models with their own shortcomings and uncertainties. The goal of this paper is to give a new evaluation of the cosmic antiproton flux along with the associated uncertainties. The propagation parameters were tightly constrained in Maurin et al. 2001 by an analysis of cosmic ray nuclei data in the framework of a two-zone diffusion model and we apply these parameters to the propagation of antiprotons. We use the recently published data on proton and helion fluxes, and we find that this particular source of uncertainty has become negligible. The Monte Carlo program DTUNUC was used to carefully examine nuclear reactions. We find that all the cosmic antiproton fluxes naturally coming out of the calculation are fully compatible with experimental data. Uncertainties in this flux have been strongly reduced. Those related to propagation are less than 25%. All other possible sources of uncertainty have also been studied.

Cosmic ray nuclei fluxes are expected to be measured with high precision in the near future. For instance, high quality data on the antiproton component could give important clues about the nature of the astronomical dark matter. A very good understanding of the different aspects of cosmic ray propagation is therefore necessary. In this paper, we use cosmic ray nuclei data to give constraints on the diffusion parameters. Propagation is studied with semi-analytical solutions of a diffusion model, and we give new analytical solutions for radioactively produced species. Our model includes convection and reacceleration as well as the standard energy losses. We perform a $\chi^2$ analysis over B/C data for a large number of configurations obtained by varying the relevant parameters of the diffusion model. A very good agreement with B/C data arises for a number of configurations, all of which are compatible with sub-Fe/Fe data. Different source spectra $Q(E)$ and diffusion coefficients $K(E)$ have been tried, but for both parameters only one form gives a good fit. Another important result is that models without convection or without reacceleration are excluded. We find that the various parameters, i.e. the diffusion coefficient normalisation $K_0$ and spectral index $\delta$, the halo thickness $L$, the Alfv\'en velocity $V_a$, and the convection velocity $V_c$ are strongly correlated. We obtain limits on the spectral index $\delta$ of thediffusion coefficient, and in particular we exclude a Kolmogorov spectrum ($\delta = 1/3$).