We update our previous search for trapped magnetic monopoles in LHC Run 2 using nearly six times more integrated luminosity and including additional models for the interpretation of the data. The MoEDAL forward trapping detector, comprising 222~kg of aluminium samples, was exposed to 2.11~fb$^{-1}$ of 13 TeV proton-proton collisions near the LHCb interaction point and analysed by searching for induced persistent currents after passage through a superconducting magnetometer. Magnetic charges equal to the Dirac charge or above are excluded in all samples. The results are interpreted in Drell-Yan production models for monopoles with spins 0, 1/2 and 1: in addition to standard point-like couplings, we also consider couplings with momentum-dependent form factors. The search provides the best current laboratory constraints for monopoles with magnetic charges ranging from two to five times the Dirac charge.

MoEDAL is designed to identify new physics in the form of long-lived highly-ionising particles produced in high-energy LHC collisions. Its arrays of plastic nuclear-track detectors and aluminium trapping volumes provide two independent passive detection techniques. We present here the results of a first search for magnetic monopole production in 13 TeV proton-proton collisions using the trapping technique, extending a previous publication with 8 TeV data during LHC run-1. A total of 222 kg of MoEDAL trapping detector samples was exposed in the forward region and analysed by searching for induced persistent currents after passage through a superconducting magnetometer. Magnetic charges exceeding half the Dirac charge are excluded in all samples and limits are placed for the first time on the production of magnetic monopoles in 13 TeV $pp$ collisions. The search probes mass ranges previously inaccessible to collider experiments for up to five times the Dirac charge.

We present the results of the most complete ever scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine learning package. This is the first such study to separate out low-mass isotopes ($p$, $\bar p$ and He) from the usual light elements (Be, B, C, N, O). We find that the propagation parameters that best fit $p$, $\bar p$, He data are significantly different from those that fit light elements, including the B/C and $^{10}$Be/$^9$Be secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.

We study SO(10), SU(5) and flipped SU(5) GUT models with non-universal soft supersymmetry-breaking scalar masses, exploring how they are constrained by LHC supersymmetry searches and cold dark matter experiments, and how they can be probed and distinguished in future experiments. We find characteristic differences between the various GUT scenarios, particularly in the coannihilation region, which is very sensitive to changes of parameters. For example, the flipped SU(5) GUT predict the possibility of $\tilde{t}_1-\chi$ coannihilation, which is absent in the regions of the SO(10) and SU(5) GUT parameter spaces that we study. We use the relic density predictions in different models to determine upper bounds for the neutralino masses, and we find large differences between different GUT models in the sparticle spectra for the same LSP mass, leading to direct connections of distinctive possible experimental measurements with the structure of the GUT group. We find that future LHC searches for generic missing $E_T$, charginos and stops will be able to constrain the different GUT models in complementary ways, as will the Xenon 1 ton and Darwin dark matter scattering experiments and future FERMI or CTA $\gamma$-ray searches.

Motivated by the recent diphoton excesses reported by both ATLAS and CMS collaborations, we suggest that a new heavy spinless particle is produced in gluon fusion at the LHC and decays to a couple of lighter pseudoscalars which then decay to photons. The new resonances could arise from a new strongly interacting sector and couple to Standard Model gauge bosons only via the corresponding Wess-Zumino-Witten anomaly. We present a detailed recast of the newest 13 TeV data from ATLAS and CMS together with the 8 TeV data to scan the consistency of the parameter space for those resonances.

We explore the effects of three-loop minimal supersymmetric standard model renormalisation group equation terms and some leading two-loopthreshold corrections on gauge and Yukawa unification: each being one loop higher order than current public spectrum calculators. We also explore the effect of the higher order terms (often 2-3 GeV) on the lightest CP even Higgs mass prediction. We illustrate our results in the constrained minimal supersymmetric standard model. Neglecting threshold corrections at the grand unified scale, the discrepancy between the unification scale $\alpha_s$ and the other two unified gauge couplings changes by 0.1$\%$ due to the higher order corrections and the difference between unification scale bottom-tau Yukawa couplings neglecting unification scale threshold corrections changes by up to 1$\%$. The difference between unification scale bottom and top Yukawa couplings changes by a few percent. Differences due to the higher order corrections also give an estimate of the size of theoretical uncertainties in the minimal supersymmetric standard model spectrum. We use these to provide estimates of theoretical uncertainties in predictions of the dark matter relic density (which can be of order one due to its strong dependence on sparticle masses) and the LHC sparticle production cross-section (often around 30$\%$). The additional higher order corrections have been incorporated into SOFTSUSY, and we provide details on how to compile and use the program. We also provide a summary of the approximations used in the higher order corrections.

We present statistically convergent profile likelihood maps obtained via global fits of a phenomenological Minimal Supersymmetric Standard Model with 15 free parameters (the MSSM-15), based on over 250M points. We derive constraints on the model parameters from direct detection limits on dark matter, the Planck relic density measurement and data from accelerator searches. We provide a detailed analysis of the rich phenomenology of this model, and determine the SUSY mass spectrum and dark matter properties that are preferred by current experimental constraints. We evaluate the impact of the measurement of the anomalous magnetic moment of the muon ($g-2$) on our results, and provide an analysis of scenarios in which the lightest neutralino is a subdominant component of the dark matter. The MSSM-15 parameters are relatively weakly constrained by current data sets, with the exception of the parameters related to dark matter phenomenology ($M_1$, $M_2$, $\mu$), which are restricted to the sub-TeV regime, mainly due to the relic density constraint. The mass of the lightest neutralino is found to be < 1.5 TeV at 99% C.L., but can extend up to 3 TeV when excluding the $g - 2$ constraint from the analysis. Low-mass bino-like neutralinos are strongly favoured, with spin-independent scattering cross-sections extending to very small values, $\sim 10^{-20}$ pb. ATLAS SUSY null searches strongly impact on this mass range, and thus rule out a region of parameter space that is outside the reach of any current or future direct detection experiment. The best-fit point obtained after inclusion of all data corresponds to a squark mass of 2.3 TeV, a gluino mass of 2.1 TeV and a 130 GeV neutralino with a spin-independent cross-section of $2.4 \times 10^{-10}$ pb, which is within the reach of future multi-ton scale direct detection experiments and of the upcoming LHC run at increased centre-of-mass energy.

The discovery by the ATLAS and CMS experiments of a new boson with mass around 125 GeV and with measured properties compatible with those of a Standard-Model Higgs boson, coupled with the absence of discoveries of phenomena beyond the Standard Model at the TeV scale, has triggered interest in ideas for future Higgs factories. A new circular e+e- collider hosted in a 80 to 100 km tunnel, TLEP, is among the most attractive solutions proposed so far. It has a clean experimental environment, produces high luminosity for top-quark, Higgs boson, W and Z studies, accommodates multiple detectors, and can reach energies up to the t-tbar threshold and beyond. It will enable measurements of the Higgs boson properties and of Electroweak Symmetry-Breaking (EWSB) parameters with unequalled precision, offering exploration of physics beyond the Standard Model in the multi-TeV range. Moreover, being the natural precursor of the VHE-LHC, a 100 TeV hadron machine in the same tunnel, it builds up a long-term vision for particle physics. Altogether, the combination of TLEP and the VHE-LHC offers, for a great cost effectiveness, the best precision and the best search reach of all options presently on the market. This paper presents a first appraisal of the salient features of the TLEP physics potential, to serve as a baseline for a more extensive design study.

We study in detail the impact of the current uncertainty in nucleon matrix elements on the sensitivity of direct and indirect experimental techniques for dark matter detection. We perform two scans in the framework of the cMSSM: one using recent values of the pion-sigma term obtained from Lattice QCD, and the other using values derived from experimental measurements. The two choices correspond to extreme values quoted in the literature and reflect the current tension between different ways of obtaining information about the structure of the nucleon. All other inputs in the scans, astrophysical and from particle physics, are kept unchanged. We use two experiments, XENON100 and IceCube, as benchmark cases to illustrate our case. We find that the interpretation of dark matter search results from direct detection experiments is more sensitive to the choice of the central values of the hadronic inputs than the results of indirect search experiments. The allowed regions of cMSSM parameter space after including XENON100 constrains strongly differ depending on the assumptions on the hadronic matrix elements used. On the other hand, the constraining potential of IceCube is almost independent of the choice of these values.

We present global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM) and the Non-Universal Higgs Model (NUHM), including the most recent CMS constraint on the Higgs boson mass, 5.8/fb integrated luminosity null Supersymmetry searches by ATLAS, the new LHCb measurement of the Bs to mu+mu- branching ratio and the 7-year WMAP dark matter relic abundance determination. We include the latest dark matter constraints from the XENON100 experiment, marginalising over astrophysical and particle physics uncertainties. We present Bayesian posterior and profile likelihood maps of the highest resolution available today, obtained from up to 350M points. We find that the new constraint on the Higgs boson mass has a dramatic impact, ruling out large regions of previously favoured cMSSM and NUHM parameter space. In the cMSSM, light sparticles and predominantly gaugino-like dark matter with a mass of a few hundred GeV are favoured. The NUHM exhibits a strong preference for heavier sparticle masses and a Higgsino-like neutralino with a mass of 1 TeV. The future ton-scale XENON1T direct detection experiment will probe large portions of the currently favoured cMSSM and NUHM parameter space. The LHC operating at 14 TeV collision energy will explore the favoured regions in the cMSSM, while most of the regions favoured in the NUHM will remain inaccessible. Our best-fit points achieve a satisfactory quality-of-fit, with p-values ranging from 0.21 to 0.35, so that none of the two models studied can be presently excluded at any meaningful significance level.

We present new global fits of the cMSSM, including LHC 1/fb integrated luminosity SUSY exclusion limits, recent LHC 5/fb constraints on the mass of the Higgs boson and XENON100 direct detection data. Our analysis fully takes into account astrophysical and hadronic uncertainties that enter the analysis when translating direct detection limits into constraints on the cMSSM parameter space. We provide results for both a Bayesian and a Frequentist statistical analysis. We find that LHC 2011 constraints in combination with XENON100 data can rule out a significant portion of the cMSSM parameter space. Our results further emphasise the complementarity of collider experiments and direct detection searches in constraining extensions of Standard Model physics. The LHC 2011 exclusion limit strongly impacts on low-mass regions of cMSSM parameter space, such as the stau co-annihilation region, while direct detection data can rule out regions of high SUSY masses, such as the Focus-Point region, which is unreachable for the LHC in the near future. We show that, in addition to XENON100 data, the experimental constraint on the anomalous magnetic moment of the muon plays a dominant role in disfavouring large scalar and gaugino masses. We find that, should the LHC 2011 excess hinting towards a Higgs boson at 126 GeV be confirmed, currently favoured regions of the cMSSM parameter space will be robustly ruled out from both a Bayesian and a profile likelihood statistical perspective.

We present updated global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM), including the most recent constraints from the ATLAS and CMS detectors at the LHC, as well as the most recent results of the XENON100 experiment. Our robust analysis takes into account both astrophysical and hadronic uncertainties that enter in the calculation of the rate of WIMP-induced recoils in direct detection experiment. We study the consequences for neutralino Dark Matter, and show that current direct detection data already allow to robustly rule out the so-called Focus Point region, therefore demonstrating the importance of particle astrophysics experiments in constraining extensions of the Standard Model of Particle Physics. We also observe an increased compatibility between results obtained from a Bayesian and a Frequentist statistical perspective. We find that upcoming ton-scale direct detection experiments will probe essentially the entire currently favoured region (at the 99% level), almost independently of the statistical approach used. Prospects for indirect detection of the cMSSM are further reduced.

Even if Supersymmetric particles are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the Dark Matter (DM) in the Universe using LHC data alone. We study the complementarity of LHC and DM indirect searches, working out explicitly the reconstruction of the DM properties for a specific benchmark model in the coannihilation region of a 24-parameters supersymmetric model. Combining mock high-luminosity LHC data with present-day null searches for gamma-rays from dwarf galaxies with the Fermi LAT, we show that current Fermi LAT limits already have the capability of ruling out a spurious Wino-like solution that would survive using LHC data only, thus leading to the correct identification of the cosmological solution. We also demonstrate that upcoming Planck constraints on the reionization history will have a similar constraining power, and discuss the impact of a possible detection of gamma-rays from DM annihilation in Draco with a CTA-like experiment. Our results indicate that indirect searches can be strongly complementary to the LHC in identifying the DM particles, even when astrophysical uncertainties are taken into account.

Statistical inference of the fundamental parameters of supersymmetric theories is a challenging and active endeavor. Several sophisticated algorithms have been employed to this end. While Markov-Chain Monte Carlo (MCMC) and nested sampling techniques are geared towards Bayesian inference, they have also been used to estimate frequentist confidence intervals based on the profile likelihood ratio. We investigate the performance and appropriate configuration of MultiNest, a nested sampling based algorithm, when used for profile likelihood-based analyses both on toy models and on the parameter space of the Constrained MSSM. We find that while the standard configuration is appropriate for an accurate reconstruction of the Bayesian posterior, the profile likelihood is poorly approximated. We identify a more appropriate MultiNest configuration for profile likelihood analyses, which gives an excellent exploration of the profile likelihood (albeit at a larger computational cost), including the identification of the global maximum likelihood value. We conclude that with the appropriate configuration MultiNest is a suitable tool for profile likelihood studies, indicating previous claims to the contrary are not well founded.

Research in many areas of modern physics such as, e.g., indirect searches for dark matter and particle acceleration in SNR shocks, rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions, The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, that uses astrophysical information, nuclear and particle data as input to self-consistently predict CRs, gamma rays, synchrotron and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.

The start of LHC has motivated an effort to determine the relative probability of the different regions of the MSSM parameter space, taking into account the present, theoretical and experimental, wisdom about the model. Since the present experimental data are not powerful enough to select a small region of the MSSM parameter space, the choice of a judicious prior probability for the parameters becomes most relevant. Previous studies have proposed theoretical priors that incorporate some (conventional) measure of the fine-tuning, to penalize unnatural possibilities. However, we show that such penalization arises from the Bayesian analysis itself (with no ad hoc assumptions), upon the marginalization of the mu-parameter. Furthermore the resulting effective prior contains precisely the Barbieri-Giudice measure, which is very satisfactory. On the other hand we carry on a rigorous treatment of the Yukawa couplings, showing in particular that the usual practice of taking the Yukawas "as required", approximately corresponds to taking logarithmically flat priors in the Yukawa couplings. Finally, we use an efficient set of variables to scan the MSSM parameter space, trading in particular B by tan beta, giving the effective prior in the new parameters. Beside the numerical results, we give accurate analytic expressions for the effective priors in all cases. Whatever experimental information one may use in the future, it is to be weighted by the Bayesian factors worked out here.

We use cosmic microwave background and large scale structure data to test a broad and physically well-motivated class of inflationary models: those with flat tree-level potentials (typical in supersymmetry). The non-trivial features of the potential arise from radiative corrections which give a simple logarithmic dependence on the inflaton field, making the models very predictive. We also consider a modified scenario with new physics beyond a certain high-energy cut-off showing up as non-renormalizable operators (NRO) in the inflaton field. We find that both kinds of models fit remarkably well CMB and LSS data, with very few free parameters. Besides, a large part of these models naturally predict a reasonable number of e-folds. A robust feature of these scenarios is the smallness of tensor perturbations (r < 10^{-3}). The NRO case can give a sizeable running of the spectral index while achieving a sufficient number of e-folds. We use Bayesian model comparison tools to assess the relative performance of the models. We believe that these scenarios can be considered as a standard physical class of inflationary models, on a similar footing with monomial potentials.

We update and extend to larger masses our previous analysis of the MSSM with minimal $SO_{10}$ [MSO$_{10}$SM] soft SUSY breaking boundary conditions. We find a well--defined, narrow region of parameter space which provides the observed relic density of dark matter, in a domain selected to fit precision electroweak data, including top, bottom and tau masses. The model is highly constrained which allows us to make several predictions. We find the light Higgs mass $m_h \leq 121 \pm 3$ GeV and also upper bounds on the mass of the gluino $\mgluino\lsim3.1$ TeV and lightest neutralino $\mchi\lsim450$ GeV. As the CP odd Higgs mass $m_A$ increases, the region of parameter space consistent with WMAP data is forced to larger values of $M_{1/2}$ and smaller values of $m_h$. Hence, we find an upper bound $m_A \lsim 1.3$ TeV. This in turn leads to lower bounds on ${\rm BR}(B_s\to \mu^+ \mu^-) > 10^{-8}$ (assuming minimal flavor violation) and on the dark matter spin independent detection cross section $\sigsip > 10^{-9}$ pb. Finally, we extend our previous analysis to include WIMP signals in indirect detection and find prospects for WIMP detection generally much less promising than in direct WIMP searches.

CMSSM boundary conditions are usually used when calculating cosmological dark matter densities. In this paper we calculate the cosmological density of dark matter in the MSSM using minimal $SO_{10}$ soft SUSY breaking boundary conditions. These boundary conditions incorporate several attractive features: they are consistent with $SO_{10}$ Yukawa unification, they result in a "natural" inverted scalar mass hierarchy and they reduce the dimension 5 operator contribution to the proton decay rate. With regards to dark matter, on the other hand, this is to a large extent an unexplored territory with large squark and slepton masses $m_{16}$, large $A_0$ and small $ \{\mu, M_{1/2} \} $. We find that in most regions of parameter space the cosmological density of dark matter is considerably less than required by the data. However there is a well--defined, narrow region of parameter space which provides the observed relic density of dark matter, as well as a good fit to precision electroweak data, including top, bottom and tau masses, and acceptable bounds on the branching fraction of $B_s \to \mu^+ \mu^-$. We present predictions for Higgs and SUSY spectra, the dark matter detection cross section and the branching ratio ${\rm BR}(B_s\to \mu^+ \mu^-)$ in this region of parameter space.