
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 protonproton 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 DrellYan
production models for monopoles with spins 0, 1/2 and 1: in addition to
standard pointlike couplings, we also consider couplings with
momentumdependent 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 longlived
highlyionising particles produced in highenergy LHC collisions. Its arrays of
plastic nucleartrack 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 protonproton
collisions using the trapping technique, extending a previous publication with
8 TeV data during LHC run1. 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 lowmass 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
secondarytoprimary 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 nonuniversal soft
supersymmetrybreaking 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 WessZuminoWitten 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 threeloop minimal supersymmetric standard model
renormalisation group equation terms and some leading twoloopthreshold
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 23 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 bottomtau 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 crosssection (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 MSSM15), 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 ($g2$) on our results, and provide an
analysis of scenarios in which the lightest neutralino is a subdominant
component of the dark matter. The MSSM15 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
subTeV 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. Lowmass
binolike neutralinos are strongly favoured, with spinindependent scattering
crosssections 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 bestfit 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 spinindependent crosssection of $2.4 \times 10^{10}$ pb,
which is within the reach of future multiton scale direct detection
experiments and of the upcoming LHC run at increased centreofmass 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
StandardModel 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
topquark, Higgs boson, W and Z studies, accommodates multiple detectors, and
can reach energies up to the ttbar threshold and beyond. It will enable
measurements of the Higgs boson properties and of Electroweak SymmetryBreaking
(EWSB) parameters with unequalled precision, offering exploration of physics
beyond the Standard Model in the multiTeV range. Moreover, being the natural
precursor of the VHELHC, a 100 TeV hadron machine in the same tunnel, it
builds up a longterm vision for particle physics. Altogether, the combination
of TLEP and the VHELHC 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 pionsigma 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 NonUniversal 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 7year 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 gauginolike dark matter with a mass
of a few hundred GeV are favoured. The NUHM exhibits a strong preference for
heavier sparticle masses and a Higgsinolike neutralino with a mass of 1 TeV.
The future tonscale 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 bestfit points achieve a satisfactory qualityoffit, with pvalues
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 lowmass regions of cMSSM
parameter space, such as the stau coannihilation region, while direct
detection data can rule out regions of high SUSY masses, such as the
FocusPoint 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
WIMPinduced 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 socalled 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 tonscale 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 24parameters supersymmetric model. Combining mock highluminosity LHC
data with presentday null searches for gammarays from dwarf galaxies with the
Fermi LAT, we show that current Fermi LAT limits already have the capability of
ruling out a spurious Winolike 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 gammarays from DM annihilation in Draco with a CTAlike
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 MarkovChain 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 likelihoodbased 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, Xrays,
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 semianalytical 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
selfconsistently 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
bestfit 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
finetuning, 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 muparameter. Furthermore the
resulting effective prior contains precisely the BarbieriGiudice 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 wellmotivated class of inflationary models: those with
flat treelevel potentials (typical in supersymmetry). The nontrivial 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 highenergy cutoff showing up as nonrenormalizable 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 efolds. 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 efolds. 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 welldefined, 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
welldefined, 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.