Big-Bang nucleosynthesis (BBN) represents one of the earliest phenomena that can lead to observational constraints on the early Universe properties. It is well-known that many important mechanisms and phase transitions occurred before BBN. We discuss the possibility of gaining insight into the primordial Universe through studies of dark matter in cosmology, astroparticle physics and colliders. For this purpose, we assume that dark matter is a thermal relic, and show that combining collider searches with dark matter observables can lead to strong constraints on the period of freeze-out before BBN.

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 show how knowledge of the cold dark matter (CDM) density can be used, in conjunction with measurements of the parameters of a scenario for beyond the Standard Model (BSM) physics, to provide information about the evolution of the Universe before Big Bang Nucleosynthesis (BBN). As examples of non-standard evolution, we consider models with a scalar field that may decay into BSM particles, and quintessence models. We illustrate our calculations using various supersymmetric models as representatives of classes of BSM scenarios in which the CDM density is either larger or smaller than the observed density when the early Universe is assumed to be radiation-dominated. In the case of a decaying scalar field, we show how the CDM density can constrain the initial scalar density and the reheating temperature after it decays in BSM scenarios that would yield overdense dark matter in standard radiation-dominated cosmology, and how the decays of the scalar field into BSM particles can be constrained in scenarios that would otherwise yield underdense CDM. We also show how the early evolution of the quintessence field can be constrained in BSM scenarios.

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 describe SuperIso Relic, a public program for the calculation of dark matter relic density and direct and indirect detection rates, which includes in addition the SuperIso routines for the calculation of flavour physics observables. SuperIso Relic v4 incorporates many new features, namely the possibility of multiprocessor calculation of the relic density, new cosmological models, and the implementations of the calculation of the observables related to direct and indirect detection experiments. Furthermore, the new version includes an implementation of the nuclear and astrophysical uncertainties, from namely nuclear form factors, dark matter density and velocity, as well as cosmic ray propagation through the galactic medium.

We study the implications of dark matter searches, together with collider constraints, on the phenomenological MSSM with neutralino dark matter and focus on the consequences of the related uncertainties in some detail. We consider, inter alia, the latest results from AMS-02, Fermi-LAT and XENON1T. In particular, we examine the impact of the choice of the dark matter halo profile, as well as the propagation model for cosmic rays, for dark matter indirect detection and show that the constraints on the MSSM differ by one to two orders of magnitude depending on the astrophysical hypotheses. On the other hand, our limited knowledge of the local relic density in the vicinity of the Earth and the velocity of Earth in the dark matter halo leads to a factor 3 in the exclusion limits obtained by direct detection experiments. We identified the astrophysical models leading to the most conservative and the most stringent constraints and for each case studied the complementarities with the latest LHC measurements and limits from Higgs, SUSY and monojet searches. We show that combining all data from dark matter searches and colliders, a large fraction of our supersymmetric sample could be probed. Whereas the direct detection constraints are rather robust under the astrophysical assumptions, the uncertainties related to indirect detection can have an important impact on the number of the excluded points.

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.

It is well known that dark matter density measurements, indirect and direct detection experiments, importantly complement the LHC in setting strong constraints on new physics scenarios. Yet, dark matter searches are subject to limitations which need to be considered for realistic analyses. For illustration, we explore the parameter space of the phenomenological MSSM and discuss the interplay of the constraints from dark matter searches and the LHC, and analyse the impact of the astrophysical uncertainties in some detail.

We revisit MSSM scenarios with light neutralino as a dark matter candidate in view of the latest LHC and dark matter direct and indirect detection experiments. We show that scenarios with a very light neutralino (~ 10 GeV) and a scalar bottom quark close in mass, can satisfy all the available constraints from LEP, Tevatron, LHC, flavour and low energy experiments and provide solutions in agreement with the bulk of dark matter direct detection experiments, and in particular with the recent CDMS results.

We investigate the viability of light neutralino scenarios as promoted by dark matter direct detection experiments. Using high statistics scans in the pMSSM we have identified several scenarios which give rise to very light neutralinos with large direct detection scattering cross sections. Our results are challenged with constraints from dark matter relic density, direct detection, indirect detection, as well as flavour physics, electroweak precision tests, LEP and Tevatron limits, LHC limits on SUSY, Higgs and monojet searches. In particular we study the effect of a Higgs boson in the range 122.5 < Mh < 127.5 GeV. We show that several scenarios emerge in agreement with all the constraints, and we study their characteristics and the LHC sensitivity.

The sensitivity of the lightest supersymmetric particle relic density calculation to different cosmological scenarios is discussed. In particular, we investigate the effects of modifications of the expansion rate and of the entropy content in the Early Universe. These effects, even with no observational consequences, can still drastically modify the relic density constraints on the SUSY parameter space. We suggest general parametrizations to evaluate such effects, and derive also constraints from Big-Bang nucleosynthesis. We show that using the relic density in the context of supersymmetric constraints requires a clear statement of the underlying cosmological model assumptions to avoid misinterpretations. On the other hand, we note that combining the relic density calculation with the eventual future discoveries at the LHC will hopefully shed light on the Very Early Universe properties.

Relic density calculations are often used to constrain particle physics models, and in particular supersymmetry. We will show that the presence of additional energy or entropy before the Big-Bang nucleosynthesis can however completely change the relic density constraints on the SUSY parameter space. Therefore one should be extremely careful when using the relic density to constrain supersymmetry as it could give misleading results, especially if combined with the future collider data. Alternatively, we will also show that combining the discoveries of the future colliders with relic density calculations can shed light on the inaccessible pre-BBN dark time physics. Finally we will present SuperIso Relic, a new relic density calculator code in Supersymmetry, which incorporates alternative cosmological models, and is publicly available.

We describe SuperIso Relic, a public program for evaluation of relic density and flavor physics observables in the minimal supersymmetric extension of the Standard Model (MSSM). SuperIso Relic is an extension of the SuperIso program which adds to the flavor observables of SuperIso the computation of all possible annihiliation and coannihilation processes of the LSP which are required for the relic density calculation. All amplitudes have been generated at the tree level with FeynArts/FormCalc, and widths of the Higgs bosons are computed with FeynHiggs at the two-loop level. SuperIso Relic also provides the possibility to modify the assumptions of the cosmological model, and to study their consequences on the relic density.

The sensitivity of the lightest supersymmetric particle relic density calculation to the variation of the cosmological expansion rate before nucleosynthesis is discussed. We show that such a modification, even extremely modest and with no consequence on the cosmological observations, can greatly enhance the calculated relic density, and therefore change the constraints on the SUSY parameter space drastically. We illustrate this variation in two examples of SUSY models, and show that it is unsafe to use the lower bound of the WMAP limits in order to constrain supersymmetry. We therefore suggest to use only the upper value Omega_DM h^2 < 0.135.

We study the loop corrections to potentials of complex or coupled real scalar fields used in cosmology to account for dark energy, dark matter or dark fluid. We show that the SUGRA quintessence and dark matter scalar field potentials are stable against the quantum fluctuations, and we propose solutions to the instability of the potentials of coupled quintessence and dark fluid scalar fields. We also find that a coupling to fermions is very restricted, unless this coupling has a structure which already exists in the scalar field potential or which can be compensated by higher order corrections. Finally, we study the influence of the curvature and kinetic term corrections.