In the new precision era for cosmic ray astrophysics, theoretical predictions cannot content themselves with average trends, but need to correctly take into account intrinsic uncertainties. The space-time discreteness of the cosmic ray sources, together with a substantial ignorance of their precise epochs and locations (with the possible exception of the most recent and close ones) play an important role in this sense. We elaborate a statistical theory to deal with this problem, relating the composite probability P({\Psi}) to obtain a flux {\Psi} at the Earth and the single-source probability p({\psi}) to contribute with a flux {\psi}. The main difficulty arises from the fact that p({\psi}) is a fat tail distribution, characterized by power-law or broken power-law behavior up to very large fluxes for which central limit theorem does not hold, and leading to well-known stable laws as opposed Gaussian. We find that relatively simple recipes provide a satisfactory description of the probability P({\Psi}). We also find that a naive Gaussian fit to simulation results would underestimate the probability of very large fluxes, that is, several times above the average, while overestimating the probability of relatively milder excursions. At large energies, large flux fluctuations are prevented by causal considerations, while at low energies a partial knowledge on the recent and nearby population of sources plays an important role. A few proposal have been discussed in the literature to account for spectral breaks recently reported in cosmic ray data in terms of local contributions. We apply our newly developed theory to assess their probabilities, finding that they are relatively small.

A hardening of the proton and helium fluxes is observed above a few hundreds of GeV/nuc. The distribution of local sources of primary cosmic rays has been suggested as a potential solution to this puzzling behavior. Some authors even claim that a single source is responsible for the observed anomalies. But how probable these explanations are? To answer that question, our current description of cosmic ray Galactic propagation needs to be replaced by the Myriad model. In the former approach, sources of protons and helium nuclei are treated as a jelly continuously spread over space and time. A more accurate description is provided by the Myriad model where sources are considered as point-like events. This leads to a probabilistic derivation of the fluxes of primary species, and opens the possibility that larger-than-average values may be observed at the Earth. For a long time though, a major obstacle has been the infinite variance associated to the probability distribution function which the fluxes follow. Several suggestions have been made to cure this problem but none is entirely satisfactory. We go a step further here and solve the infinite variance problem of the Myriad model by making use of the generalized central limit theorem. We find that primary fluxes are distributed according to a stable law with heavy tail, well-known to financial analysts. The probability that the proton and helium anomalies are sourced by local SNR can then be calculated. The p-values associated to the CREAM measurements turn out to be small, unless somewhat unrealistic propagation parameters are assumed.

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).

Recent measurements of cosmic ray proton and helium spectra show a hardening above a few hundreds of GeV. This excess is hard to understand in the framework of the conventional models of Galactic cosmic ray production and propagation. We propose here to explain this anomaly by the presence of local sources (myriad model). Cosmic ray propagation is described as a diffusion process taking place inside a two-zone magnetic halo. We calculate the proton and helium fluxes at the Earth between 50 GeV and 100 TeV. Improving over a similar analysis, we consistently derive these fluxes by taking into account both local and remote sources for which a unique injection rate is assumed. We find cosmic ray propagation parameters compatible with B/C measurements and for which the proton and helium spectra remarkably agree with the PAMELA and CREAM measurements over four decades in energy.

Measurements of cosmic ray fluxes by the PAMELA and CREAM experiments show unexpected spectral features between 200 GeV and 100 TeV. They could be due to the presence of nearby and young cosmic ray sources. This can be studied in the myriad model, in which cosmic rays diffuse from point-like instantaneous sources located randomly throughout the Galaxy. To test this hypothesis, one must compute the flux due to a catalog of local sources, but also the error bars associated to this quantity. This turns out not to be as straightforward as it seems, as the standard deviation is infinite when computed for the most general statistical ensemble. The goals of this paper are to provide a method to associate error bars to the flux measurements which has a clear statistical meaning, and to explore the relation between the myriad model and the more usual source model based on a continuous distribution. To this end, we show that the quantiles of the flux distribution are well-defined, even though the standard deviation is infinite. They can be used to compute 68% confidence levels, for instance. We also use the fact that local sources have known positions and ages to reduce the statistical ensemble from which random sources are drawn in the myriad model. We present a method to evaluate meaningful error bars for the flux obtained in the myriad model. In this context, we also discuss the status of the spectral features observed in the proton flux by CREAM and PAMELA.

Secondary positrons are produced by spallation of cosmic rays within the interstellar gas. Measurements have been typically expressed in terms of the positron fraction, which exhibits an increase above 10 GeV. Many scenarios have been proposed to explain this feature, among them some additional primary positrons originating from dark matter annihilation in the Galaxy. The PAMELA satellite has provided high quality data that has enabled high accuracy statistical analyses to be made, showing that the increase in the positron fraction extends up to about 100 GeV. It is therefore of paramount importance to constrain theoretically the expected secondary positron flux to interpret the observations in an accurate way. We find the secondary positron flux to be reproduced well by the available observations, and to have theoretical uncertainties that we quantify to be as large as about one order of magnitude. We also discuss the positron fraction issue and find that our predictions may be consistent with the data taken before PAMELA. For PAMELA data, we find that an excess is probably present after considering uncertainties in the positron flux, although its amplitude depends strongly on the assumptions made in relation to the electron flux. By fitting the current electron data, we show that when considering a soft electron spectrum, the amplitude of the excess might be far lower than usually claimed. We provide fresh insights that may help to explain the positron data with or without new physical model ingredients. PAMELA observations and the forthcoming AMS-02 mission will allow stronger constraints to be aplaced on the cosmic--ray transport parameters, and are likely to reduce drastically the theoretical uncertainties.

Recent PAMELA and ATIC data seem to indicate an excess in positron cosmic rays above approximately 10 GeV which might be due to galactic Dark Matter particle annihilation. However the background of this signal suffers many uncertainties that make our task difficult in constraining Dark Matter or any other astrophysical explanation for these recent surprising data.

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.

The annihilation of dark matter (DM) in the Galaxy could produce specific imprints on the spectra of antimatter species in Galactic cosmic rays, which could be detected by upcoming experiments such as PAMELA and AMS02. Recent studies show that the presence of substructures can enhance the annihilation signal by a "boost factor" that not only depends on energy, but that is intrinsically a statistical property of the distribution of DM substructures inside the Milky Way. We investigate a scenario in which substructures consist of $\sim 100$ "mini-spikes" around intermediate-mass black holes. Focusing on primary positrons and antiprotons, we find large boost factors, up to a few thousand, that exhibit a large variance at high energy in the case of positrons and at low energy in the case of antiprotons. As a consequence, an estimate of the DM particle mass based on the observed cut-off in the positron spectrum could lead to a substantial underestimate of its actual value.

We present a measurement of the distant Type Ia supernova rate derived from the first two years of the Canada -- France -- Hawaii Telescope Supernova Legacy Survey. We observed four one-square degree fields with a typical temporal frequency of <Delta t> ~ 4 observer-frame days over time spans of from 158 to 211 days per season for each field, with breaks during full moon. We used 8-10 meter-class telescopes for spectroscopic followup to confirm our candidates and determine their redshifts. Our starting sample consists of 73 spectroscopically verified Type Ia supernovae in the redshift range 0.2 < z < 0.6. We derive a volumetric SN Ia rate of r_V(<z>=0.47) = 0.42^{+0.13}_{-0.09} (systematic) +- 0.06 (statistical) X 10^-4 yr^-1 Mpc^3, assuming h = 0.7, Omega_m = 0.3 and a flat cosmology. Using recently published galaxy luminosity functions derived in our redshift range, we derive a SN Ia rate per unit luminosity of r_L(<z>=0.47) = 0.154^{+0.048}_{-0.033} (systematic) ^{+0.039}_{-0.031} (statistical) SNu. Using our rate alone, we place an upper limit on the component of SN Ia production that tracks the cosmic star formation history of 1 SN Ia per 10^3 M_sun of stars formed. Our rate and other rates from surveys using spectroscopic sample confirmation display only a modest evolution out to z=0.55.

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.

We present a progress report on a project to derive the evolution of the volumetric supernova Type Ia rate from the Supernova Legacy Survey. Our preliminary estimate of the rate evolution divides the sample from Neill et al. (2006) into two redshift bins: 0.2 < z < 0.4, and 0.4 < z < 0.6. We extend this by adding a bin from the sample analyzed in Sullivan et al. (2006) in the range 0.6 < z < 0.75 from the same time period. We compare the derived trend with previously published rates and a supernova Type Ia production model having two components: one component associated closely with star formation and an additional component associated with host galaxy mass. Our observed trend is consistent with this model, which predicts a rising SN Ia rate out to at least z=2.

We compare the rise times of nearby and distant Type Ia supernovae (SNe Ia) as a test for evolution using 73 high-redshift spectroscopically-confirmed SNe Ia from the first two years of the five year Supernova Legacy Survey (SNLS) and published observations of nearby SN. Because of the ``rolling'' search nature of the SNLS, our measurement is approximately 6 times more precise than previous studies, allowing for a more sensitive test of evolution between nearby and distant supernovae. Adopting a simple $t^2$ early-time model (as in previous studies), we find that the rest-frame $B$ rise times for a fiducial SN Ia at high and low redshift are consistent, with values $19.10^{+0.18}_{-0.17}({stat}) \pm 0.2 ({syst})$ and $19.58^{+0.22}_{-0.19}$ days, respectively; the statistical significance of this difference is only 1.4 \sg . The errors represent the uncertainty in the mean rather than any variation between individual SN. We also compare subsets of our high-redshift data set based on decline rate, host galaxy star formation rate, and redshift, finding no substantive evidence for any subsample dependence.

The small-scale distribution of dark matter in Galactic halos is poorly known. Several studies suggest that it could be very clumpy, which turns out to be of paramount importance when investigating the annihilation signal from exotic particles (e.g. supersymmetric or Kaluza-Klein). In this paper we focus on the annihilation signal in positrons. We estimate the associated uncertainty, due to the fact that we do not know exactly how the clumps are distributed in the Galactic halo. To this aim, we perform a statistical study based on analytical computations, as well as numerical simulations. In particular, we study the average and variance of the annihilation signal over many Galactic halos having the same statistical properties. We find that the so-called boost factor used by many authors should be handled with care, as i) it depends on energy and ii) it may be different for positrons, antiprotons and gamma rays, a fact which has not received any attention before. As an illustration, we use our results to discuss the positron spectrum measurements by the HEAT experiment.

(ABRIDGED) We show that Type Ia supernovae (SNe Ia) are formed within both very young and old stellar populations, with observed rates that depend on the stellar mass and mean star-formation rates (SFRs) of their host galaxies. Models where the SN Ia rate depends solely on host galaxy stellar mass are ruled out with 99% confidence. Our analysis is based on 100 spectroscopically-confirmed SNe Ia, plus 24 photometrically-classified events, all from the Supernova Legacy Survey (SNLS) and distributed over 0.2<z<0.75. Using multi-band photometry, we estimate stellar masses and SFRs for the SN Ia host galaxies by fitting their broad-band spectral energy distributions with the galaxy spectral synthesis code, PEGASE.2. We show that the SN Ia rate per unit mass is proportional to the specific SFR of the parent galaxies -- more vigorously star-forming galaxies host more SNe Ia per unit stellar mass, broadly equivalent to the trend of increasing SN Ia rate in later-type galaxies seen in the local universe. Following earlier suggestions for a simple "two-component" model approximating the SN Ia rate, we find bivariate linear dependencies of the SN Ia rate on both the stellar masses and the mean SFRs of the host systems. We also demonstrate a dependence of distant SN Ia light-curve shapes on star-formation in the host galaxy, similar to trends observed locally. Passive galaxies, with no star-formation, preferentially host faster-declining/dimmer SNe Ia, while slower-declining/brighter events are only found in systems with ongoing star-formation. We model the light-curve width distribution in star-forming galaxies as the sum of a young component, and an old component taken from the distribution in non-star-forming galaxies.

We present a method for selecting high-redshift type Ia supernovae (SNe Ia) located via rolling SN searches. The technique, using both color and magnitude information of events from only 2-3 epochs of multi-band real-time photometry, is able to discriminate between SNe Ia and core collapse SNe. Furthermore, for the SNe Ia, the method accurately predicts the redshift, phase and light-curve parameterization of these events based only on pre-maximum-light data. We demonstrate the effectiveness of the technique on a simulated survey of SNe Ia and core-collapse SNe, where the selection method effectively rejects most core-collapse SNe while retaining SNe Ia. We also apply the selection code to real-time data acquired as part of the Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS). During the period May 2004 to January 2005 in the SNLS, 440 SN candidates were discovered of which 70 were confirmed spectroscopically as SNe Ia and 15 as core-collapse events. For this test dataset, the selection technique correctly identifies 100% of the identified SNe II as non-SNe Ia with only a 1-2% false rejection rate. The predicted parameterization of the SNe Ia has a precision of |delta_z|/(1+z_spec)<0.09 in redshift, and +/- 2-3 rest-frame days in phase, providing invaluable information for planning spectroscopic follow-up observations. We also investigate any bias introduced by this selection method on the ability of surveys such as SNLS to measure cosmological parameters (e.g., w and omega matter), and find any effect to be negligible.

We present distance measurements to 71 high redshift type Ia supernovae discovered during the first year of the 5-year Supernova Legacy Survey (SNLS). These events were detected and their multi-color light-curves measured using the MegaPrime/MegaCam instrument at the Canada-France-Hawaii Telescope (CFHT), by repeatedly imaging four one-square degree fields in four bands. Follow-up spectroscopy was performed at the VLT, Gemini and Keck telescopes to confirm the nature of the supernovae and to measure their redshift. With this data set, we have built a Hubble diagram extending to z=1, with all distance measurements involving at least two bands. Systematic uncertainties are evaluated making use of the multi-band photometry obtained at CFHT. Cosmological fits to this first year SNLS Hubble diagram give the following results : Omega_M = 0.263 +/- 0.042(stat) +/- 0.032(sys) for a flat LambdaCDM model; and w = -1.023 +/- 0.090(stat) +/- 0.054(sys) for a flat cosmology with constant equation of state w when combined with the constraint from the recent Sloan Digital Sky Survey measurement of baryon acoustic oscillations.

We present new techiques for improving the efficiency of supernova (SN) classification at high redshift using 64 candidates observed at Gemini North and South during the first year of the Supernova Legacy Survey (SNLS). The SNLS is an ongoing five year project with the goal of measuring the equation of state of Dark Energy by discovering and following over 700 high-redshift SNe Ia using data from the Canada-France-Hawaii Telescope Legacy Survey. We achieve an improvement in the SN Ia spectroscopic confirmation rate: at Gemini 71% of candidates are now confirmed as SNe Ia, compared to 54% using the methods of previous surveys. This is despite the comparatively high redshift of this sample, where the median SN Ia redshift is z=0.81 (0.155 <= z <= 1.01). These improvements were realized because we use the unprecedented color coverage and lightcurve sampling of the SNLS to predict whether a candidate is an SN Ia and estimate its redshift, before obtaining a spectrum, using a new technique called the "SN photo-z." In addition, we have improved techniques for galaxy subtraction and SN template chi^2 fitting, allowing us to identify candidates even when they are only 15% as bright as the host galaxy. The largest impediment to SN identification is found to be host galaxy contamination of the spectrum -- when the SN was at least as bright as the underlying host galaxy the target was identified more than 90% of the time. However, even SNe on bright host galaxies can be easily identified in good seeing conditions. When the image quality was better than 0.55 arcsec the candidate was identified 88% of the time. Over the five-year course of the survey, using the selection techniques presented here we will be able to add approximately 170 more confirmed SNe Ia than would be possible using previous methods.

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.

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 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.

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.

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.

It is widely believed that dark matter exists within galaxies and clusters of galaxies. Under the assumption that this dark matter is composed of the lightest, stable supersymmetric particle, assumed to be the neutralino, the feasibility of its indirect detection via observations of a diffuse gamma-ray signal due to neutralino annihilation within M31 is examined.

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.

It is widely believed that dark matter exists within galaxies and clusters of galaxies. Under the assumption that this dark matter is composed of the lightest, stable supersymmetric particle, assumed to be the neutralino, the feasibility of its indirect detection via observations of a diffuse gamma-ray signal due to neutralino annihilations within M31 is examined. To this end, first the dark matter halo of the close spiral galaxy M31 is modeled from observations, then the resultant gamma-ray flux is estimated within supersymmetric model configurations. We conclude that under favorable conditions such as the rapid accretion of neutralinos on the central black hole in M31 and/or the presence of many clumps inside its halo with $r^{-3/2}$ inner profiles, a neutralino annihilation gamma-ray signal is marginally detectable by the ongoing collaboration CELESTE.

The distant globular cluster Palomar 13 has been found to have a very high mass-to-light ratio and its profile can be well fitted either by a King model with a tail, or with a NFW model. This cluster may be the first case of the many clumps predicted by CDM simulations that would not be disrupted by the galactic halo potential. We make the hypothesis that Pal 13 is made of neutralinos and run the DarkSuspect code to estimate the high-energy photon flux due to the annihilation of neutralinos through various channels in some benchmark scenarios. These low fluxes may be used as targets to be reached in proposals for future ground-based high altitude Cerenkov telescopes.

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.

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$).

Recent observations give some clues that the lenses discovered by the microlensing experiments in the direction of the Magellanic Clouds may be located in these satellite galaxies. We re-examine the possibility that self-lensing alone may account for the optical depth measured towards the Large Magellanic Cloud (LMC). We present a self-consistent multi-component model of the LMC consisting of distinct stellar populations, each associated to a vertical velocity dispersion ranging from 10 to 60 km/s. The present work focuses on showing that such dispersions comply with current 20 - 30 km/s limits set by observation on specific LMC populations. We also show that this model reproduces both the 1-2 10^-7 observed optical depth and the event duration distribution.

We consider neutralino annihilation in dense extragalactic systems known to be dominated by dark matter, in particular M87 and several local dwarf spheroidal galaxies. These annihilations can produce energetic gamma rays which may be visible to atmospheric Cerenkov telescopes. We explore the supersymmetric parameter space, and compute the expected flux of gamma--rays coming from these objects. It is shown that some parts of the parameter space lead to a signal observable with the next generation of Cerenkov telescopes, provided the supersymmetric dark matter has a clumpy structure, as may be expected in a hierarchical scenario for structure formation.

The observed microlensing events towards the LMC do not have yet a coherent explanation. If they are due to Galactic Halo objects, the nature of these objects is puzzling --- half the halo in dark 0.5 Msol objects. On the other hand, traditional models of the LMC predict a self-lensing optical depth about an order of magnitude too low, although characteristics of some of the observed events favor a self-lensing explanation. We present here two models of the LMC taking into account the correlation between the mass of the stars and their velocity dispersion: a thin Mestel disk, and an ellipsoidal model. Both yield optical depths, event rates, and event duration distributions compatible with the observations. The grounds for such models are discussed, as well as their observational consequences.

Measurements of the surface brightness distribution and of the velocity dispersion profile have been so far used to infer the inner dynamics of globular clusters. We show that those observations do not trace back the dark matter potentially concealed in these systems in the form of low-mass compact objects. We have built Michie models of globular clusters which contain both massive and low-mass stars. An analytic expression for the stellar mass densities has been explicitely derived in terms of the usual error function and Dawson's integral. While the heavy population is kept fixed, the abundance of the light species of our models is varied. When stellar velocities are anisotropic, both the surface brightness and the velocity dispersion profiles of the cluster become insensitive to the abundance of low-mass stars. This suggests that the actual stellar mass function of many globular clusters is still to be discovered.

The leaky box model is now ruled out by measurements of a cosmic ray gradient throughout the galactic disk. It needs to be replaced by a more refined treatment which takes into account the diffusion of cosmic rays in the magnetic fields of the Galaxy. We have estimated the flux of antiprotons on the Earth in the framework of a two-zone diffusion model. Those species are created by the spallation reactions of high-energy nuclei with the interstellar gas. Another potential source of antiprotons is the annihilation of supersymmetric particles in the dark halo that surrounds our Galaxy. In this letter, we investigate both processes. Special emphasis is given to the antiproton signature of supersymmetric dark matter. The corresponding signal exceeds the conventional spallation flux below 300 MeV, a domain that will be thoroughly explored by the Antimatter Spectrometer experiment. The propagation of the antiprotons produced in the remote regions of the halo back to the Earth plays a crucial role. Depending on the energy, the leaky box estimates are wrong by a factor varying from 0.5 up to 3.