We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics.

Two years ago, the AMS collaboration released the most precise measurement of the cosmic ray positron flux. It confirms that pure secondary predictions fall below the data above 10 GeV, suggesting the presence of a primary component, e.g. annihilations of WIMPs dark matter. Most analyses have focused on the high-energy part of the spectrum, disregarding the GeV energy region where cosmic ray transport is harder to model and solar modulation comes into play. Given the high quality of AMS measurements, it is timely to re-examine the positron anomaly over the entire energy range, taking into account transport processes so far neglected, e.g. convection or diffusive re-acceleration. We devise a new semi-analytical method to take into account transport processes so far neglected, but important below a few GeV. It is based on the pinching of inverse Compton and synchrotron energy losses inside the Galactic disc. It allows to carry out extensive scans over the cosmic ray propagation parameters, which we strongly constrain by requiring that the secondary component does not overshoot the AMS measurements. Only models with large diffusion coefficients survive this test. The positron flux is a powerful and independent probe of cosmic ray propagation, complementary to the boron-to-carbon ratio. We then scan over WIMP mass to fit the annihilation cross section and branching ratios, exploring both direct annihilations into standard model particles or through light mediators. In the former case, the best fit yields a $p$-value of 0.4% for a mass of 264 GeV, a value that does not allow to reproduce the highest energy data points. Worse quality fits are found in the latter case. The interpretation of the positron excess in terms of single DM species annihilations is strongly disfavored. This conclusion is based solely on the positron data, and no other observation needs to be invoked.

Primary cosmic-ray elemental spectra have been measured with the balloon-borne Cosmic Ray Energetics And Mass (CREAM) experiment since 2004. The third CREAM payload (CREAM-III) flew for 29 days during the 2007-2008 Antarctic season. Energies of incident particles above 1 TeV are measured with a calorimeter. Individual elements are clearly separated with a charge resolution of ~0.12 e (in charge units) and ~0.14 e for protons and helium nuclei, respectively, using two layers of silicon charge detectors. The measured proton and helium energy spectra at the top of the atmosphere are harder than other existing measurements at a few tens of GeV. The relative abundance of protons to helium nuclei is 9.53+-0.03 for the range of 1 TeV/n to 63 TeV/n. The ratio is considerably smaller than other measurements at a few tens of GeV/n. The spectra become softer above ~20 TeV. However, our statistical uncertainties are large at these energies and more data are needed.

Recent years have seen increased theoretical and experimental effort towards the first-ever detection of cosmic-ray antideuterons, in particular as an indirect signature of dark matter annihilation or decay. In contrast to indirect dark matter searches using positrons, antiprotons, or gamma-rays, which suffer from relatively high and uncertain astrophysical backgrounds, searches with antideuterons benefit from very suppressed conventional backgrounds, offering a potential breakthrough in unexplored phase space for dark matter. This article is based on the first dedicated cosmic-ray antideuteron workshop, which was held at UCLA in June 2014. It reviews broad classes of dark matter candidates that result in detectable cosmic-ray antideuteron fluxes, as well as the status and prospects of current experimental searches. The coalescence model of antideuteron production and the influence of antideuteron measurements at particle colliders are discussed. This is followed by a review of the modeling of antideuteron propagation through the magnetic fields, plasma currents, and molecular material of our Galaxy, the solar system, the Earth's geomagnetic field, and the atmosphere. Finally, the three ongoing or planned experiments that are sensitive to cosmic-ray antideuterons, BESS, AMS-02, and GAPS, are detailed. As cosmic-ray antideuteron detection is a rare event search, multiple experiments with orthogonal techniques and backgrounds are essential. Many theoretical and experimental groups have contributed to these studies over the last decade, this review aims to provide the first coherent discussion of the relevant dark matter theories that antideuterons probe, the challenges to predictions and interpretations of antideuteron signals, and the experimental efforts toward cosmic antideuteron detection.

PAMELA and, more recently, AMS-02, are ushering us into a new era of greatly reduced statistical uncertainties in experimental measurements of cosmic-ray fluxes. In particular, new determinations of traditional diagnostic tools such as the boron-to-carbon ratio (B/C) are expected to significantly reduce errors on cosmic-ray diffusion parameters, with important implications for astroparticle physics, ranging from inferring primary source spectra to indirect dark matter searches. It is timely to stress, however, that the conclusions obtained crucially depend on the framework in which the data are interpreted as well as from some nuclear input parameters. We aim at assessing the theoretical uncertainties affecting the outcome, with models as simple as possible while still retaining the key dependencies. We compare different semi-analytical, two-zone model descriptions of cosmic-ray transport in the Galaxy. We test for the effect of a primary source contamination in the boron flux by parametrically altering its flux, as well as for nuclear cross section uncertainties. Our study on preliminary results from AMS-02 suggests that, differently for instance from the leptonic case, realistic modelling of the geometry of the Galaxy and of the source distribution are of minor importance to correctly reproduce B/C data at high energies and thus, to a large extent, for the extraction of diffusion parameters. The Ansatz on the lack of primary injection of boron represents the most serious bias, and requires multi-messenger studies to be addressed. If this uncertainty could be lifted, nuclear uncertainties would still represent a serious concern, which degrade the systematic error on the inferred parameters to the 20% level, or three times the estimated experimental sensitivity. In order to reduce this, a new nuclear cross section measurement campaign is probably required.

The precise measurement of cosmic-ray antiparticles serves as important means for identifying the nature of dark matter. Recent years showed that identifying the nature of dark matter with cosmic-ray positrons and higher energy antiprotons is difficult, and has lead to a significantly increased interest in cosmic-ray antideuteron searches. Antideuterons may also be generated in dark matter annihilations or decays, offering a potential breakthrough in unexplored phase space for dark matter. Low-energy antideuterons are an important approach because the flux from dark matter interactions exceeds the background flux by more than two orders of magnitude in the low-energy range for a wide variety of models. This review is based on the "dbar14 - dedicated cosmic-ray antideuteron workshop", which brought together theorists and experimentalists in the field to discuss the current status, perspectives, and challenges for cosmic-ray antideuteron searches and discusses the motivation for antideuteron searches, the theoretical and experimental uncertainties of antideuteron production and propagation in our Galaxy, as well as give an experimental cosmic-ray antideuteron search status update. This report is a condensed summary of the article "Review of the theoretical and experimental status of dark matter identification with cosmic-ray antideuteron" (arXiv:1505.07785).

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.

The positron fraction in cosmic rays was found to be a steadily increasing in function of energy, above $\sim$ 10 GeV. This behaviour contradicts standard astrophysical mechanisms, in which positrons are secondary particles, produced in the interactions of primary cosmic rays during the propagation in the interstellar medium. The observed anomaly in the positron fraction triggered a lot of excitement, as it could be interpreted as an indirect signature of the presence of dark matter species in the Galaxy. Alternatively, it could be produced by nearby astrophysical sources, such as pulsars. Both hypotheses are probed in this work in light of the latest AMS-02 positron fraction measurements. The transport of the primary and secondary positrons in the Galaxy is described using a semi-analytic two-zone model. MicrOMEGAs is used to model the positron flux generated by dark matter species. The description of the positron fraction from astrophysical sources is based on the pulsar observations included in the ATNF catalogue. We find that the mass of the favoured dark matter candidates is always larger than 500 GeV. The only dark matter species that fulfils the numerous gamma ray and cosmic microwave background bounds is a particle annihilating into four leptons through a light scalar or vector mediator, with a mixture of tau (75%) and electron (25%) channels, and a mass between 0.5 and 1 TeV. The positron anomaly can also be explained by a single astrophysical source and a list of five pulsars from the ATNF catalogue is given. Those results are obtained with the cosmic ray transport parameters that best fit the B/C ratio. Uncertainties in the propagation parameters turn out to be very significant. In the WIMP annihilation cross section to mass plane for instance, they overshadow the error contours derived from the positron data.

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.

A proximity focusing Cherenkov imager called CHERCAM, has been built for the charge measurement of nuclear cosmic rays with the CREAM instrument. It consists of a silica aerogel radiator plane across from a detector plane equipped with 1,600 1" diameter photomultipliers. The two planes are separated by a ring expansion gap. The Cherenkov light yield is proportional to the charge squared of the incident particle. The expected relative light collection accuracy is in the few percents range. It leads to an expected single element separation over the range of nuclear charge Z of main interest 1 < Z < 26. CHERCAM is designed to fly with the CREAM balloon experiment. The design of the instrument and the implemented technical solutions allowing its safe operation in high altitude conditions (radiations, low pressure, cold) are presented.

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.

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.

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.

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.

The Alpha Magnetic Spectrometer (AMS), whose final version AMS-02 is to be installed on the International Space Station (ISS) for at least 3 years, is a detector designed to measure charged cosmic ray spectra with energies up to the TeV region and with high energy photon detection capability up to a few hundred GeV. It is equipped with several subsystems, one of which is a proximity focusing RICH detector with a dual radiator (aerogel+NaF) that provides reliable measurements for particle velocity and charge. The assembly and testing of the AMS RICH is currently being finished and the full AMS detector is expected to be ready by the end of 2008. The RICH detector of AMS-02 is presented. Physics prospects are briefly discussed.

The AMS detector, to be installed on the International Space Station, includes a Ring Imaging Cerenkov detector with two different radiators, silica aerogel (n=1.05) and sodium fluoride (n=1.334). This detector is designed to provide very precise measurements of velocity and electric charge in a wide range of cosmic nuclei energies and atomic numbers. The detector geometry, in particular the presence of a reflector for acceptance purposes, leads to complex Cerenkov patterns detected in a pixelized photomultiplier matrix. The results of different reconstruction methods applied to test beam data as well as to simulated samples are presented. To ensure nominal performances throughout the flight, several detector parameters have to be carefully monitored. The algorithms developed to fulfill these requirements are presented. The velocity and charge measurements provided by the RICH detector endow the AMS spectrometer with precise particle identification capabilities in a wide energy range. The expected performances on light isotope separation are discussed.

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.