
An excess of gamma rays has been identified at the centre of the Milky Way,
and annihilation of dark matter has been posited as a potential source. This
hypothesis faces significant challenges: difficulty characterizing
astrophysical backgrounds, the need for a nontrivial adiabatic contraction of
the inner part of the Milky Way's dark matter halo, and recent observations of
photon clustering, which suggest that the majority of the excess is due to
unresolved point sources. Here we point out that the apparent pointlike nature
of the emission rules out the dark matter interpretation of the excess
entirely. Attempting to model the emission with dark matter point sources
either worsens the problem with the inner slope, requires an unrealistically
large minihalo fraction toward the Galactic Centre, or overproduces the
observed emission at higher latitudes.

We analyze the emission spectrum of the hot Jupiter WASP12b using our
HELIOSR retrieval code and HELIOSK opacity calculator. When interpreting
Hubble and Spitzer data, the retrieval outcomes are found to be
priordominated. When the prior distributions of the molecular abundances are
assumed to be loguniform, the volume mixing ratio of HCN is found to be
implausibly high. A VULCAN chemical kinetics model of WASP12b suggests that
chemical equilibrium is a reasonable assumption even when atmospheric mixing is
implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian
priors on the elemental abundances of carbon, oxygen and nitrogen with the
Gaussian peaks being centered on the measured C/H, O/H and N/H values of the
star. By enforcing chemical equilibrium, we find substellar O/H and stellar to
slightly superstellar C/H for the dayside atmosphere of WASP12b. The
superstellar carbontooxygen ratio is just above unity, regardless of whether
clouds are included in the retrieval analysis, consistent with Madhusudhan et
al. (2011). Furthermore, whether a temperature inversion exists in the
atmosphere depends on one's assumption for the Gaussian width of the priors.
Our retrieved posterior distributions are consistent with the formation of
WASP12b in a solarcomposition protoplanetary disk, beyond the water iceline,
via gravitational instability or pebble accretion (without core erosion) and
migration inwards to its present orbital location via a diskfree mechanism,
and are inconsistent with both insitu formation and core accretion with disk
migration, as predicted by Madhusudhan et al. (2017). We predict that the
interpretation of James Webb Space Telescope WASP12b data will not be
priordominated.

These notes aim at presenting an overview of Bayesian statistics, the
underlying concepts and application methodology that will be useful to
astronomers seeking to analyse and interpret a wide variety of data about the
Universe. The level starts from elementary notions, without assuming any
previous knowledge of statistical methods, and then progresses to more
advanced, researchlevel topics. After an introduction to the importance of
statistical inference for the physical sciences, elementary notions of
probability theory and inference are introduced and explained. Bayesian methods
are then presented, starting from the meaning of Bayes Theorem and its use as
inferential engine, including a discussion on priors and posterior
distributions. Numerical methods for generating samples from arbitrary
posteriors (including Markov Chain Monte Carlo and Nested Sampling) are then
covered. The last section deals with the topic of Bayesian model selection and
how it is used to assess the performance of models, and contrasts it with the
classical pvalue approach. A series of exercises of various levels of
difficulty are designed to further the understanding of the theoretical
material, including fully worked out solutions for most of them.

Euclid is a European Space Agency medium class mission selected for launch in
2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to
understand the origin of the accelerated expansion of the universe. Euclid will
explore the expansion history of the universe and the evolution of cosmic
structures by measuring shapes and redshifts of galaxies as well as the
distribution of clusters of galaxies over a large fraction of the sky. Although
the main driver for Euclid is the nature of dark energy, Euclid science covers
a vast range of topics, from cosmology to galaxy evolution to planetary
research. In this review we focus on cosmology and fundamental physics, with a
strong emphasis on science beyond the current standard models. We discuss five
broad topics: dark energy and modified gravity, dark matter, initial
conditions, basic assumptions and questions of methodology in the data
analysis. This review has been planned and carried out within Euclid's Theory
Working Group and is meant to provide a guide to the scientific themes that
will underlie the activity of the group during the preparation of the Euclid
mission.

Accurate standardisation of Type Ia supernovae (SNIa) is instrumental to the
usage of SNIa as distance indicators. We analyse a homogeneous sample of 22
lowz SNIa, observed by the Carnegie Supernova Project (CSP) in the optical and
near infrared (NIR). We study the time of the second peak in the NIR band due
to rebrightening, t2, as an alternative standardisation parameter of SNIa peak
brightness. We use BAHAMAS, a Bayesian hierarchical model for SNIa cosmology,
to determine the residual scatter in the Hubble diagram. We find that in the
absence of a colour correction, t2 is a better standardisation parameter
compared to stretch: t2 has a 1 sigma posterior interval for the Hubble
residual scatter of [0.250, 0.257] , compared to [0.280, 0.287] when stretch
(x1) alone is used. We demonstrate that when employed together with a colour
correction, t2 and stretch lead to similar residual scatter. Using colour,
stretch and t2 jointly as standardisation parameters does not result in any
further reduction in scatter, suggesting that t2 carries redundant information
with respect to stretch and colour. With a much larger SNIa NIR sample at
higher redshift in the future, t2 could be a useful quantity to perform
robustness checks of the standardisation procedure.

We present global fits of an effective field theory description of real, and
complex scalar dark matter candidates. We simultaneously take into account all
possible dimension 6 operators consisting of dark matter bilinears and gauge
invariant combinations of quark and gluon fields. We derive constraints on the
free model parameters for both the real (five parameters) and complex (seven)
scalar dark matter models obtained by combining Planck data on the cosmic
microwave background, direct detection limits from LUX, and indirect detection
limits from the Fermi Large Area Telescope. We find that for real scalars
indirect dark matter searches disfavour a dark matter particle mass below 100
GeV. For the complex scalar dark matter particle current data have a limited
impact due to the presence of operators that lead to pwave annihilation, and
also do not contribute to the spinindependent scattering cross section.
Although current data are not informative enough to strongly constrain the
theory parameter space, we demonstrate the power of our formalism to
reconstruct the theoretical parameters compatible with an actual dark matter
detection, by assuming that the excess of gamma rays observed by the Fermi
Large Area Telescope towards the Galactic centre is entirely due to dark matter
annihilations. Please note that the excess can very well be due to
astrophysical sources such as millisecond pulsars. We find that scalar dark
matter interacting via effective field theory operators can in principle
explain the Galactic centre excess, but that such interpretation is in strong
tension with the nondetection of gamma rays from dwarf galaxies in the real
scalar case. In the complex scalar case there is enough freedom to relieve the
tension.

Euclid is a European Space Agency medium class mission selected for launch in
2019 within the Cosmic Vision 20152025 programme. The main goal of Euclid is
to understand the origin of the accelerated expansion of the Universe. Euclid
will explore the expansion history of the Universe and the evolution of cosmic
structures by measuring shapes and redshifts of galaxies as well as the
distribution of clusters of galaxies over a large fraction of the sky. Although
the main driver for Euclid is the nature of dark energy, Euclid science covers
a vast range of topics, from cosmology to galaxy evolution to planetary
research. In this review we focus on cosmology and fundamental physics, with a
strong emphasis on science beyond the current standard models. We discuss five
broad topics: dark energy and modified gravity, dark matter, initial
conditions, basic assumptions and questions of methodology in the data
analysis. This review has been planned and carried out within Euclid's Theory
Working Group and is meant to provide a guide to the scientific themes that
will underlie the activity of the group during the preparation of the Euclid
mission.

We investigate the implications for inflation of the detection of Bmodes
polarization in the Cosmic Microwave Background (CMB) by BICEP2. We show that
the hypothesis of primordial origin of the measurement is only favored by the
first four bandpowers, while the others would prefer unreasonably large values
of the tensortoscalar ratio. Using only those four bandpowers, we carry out a
complete analysis in the cosmological and inflationary slowroll parameter
space using the BICEP2 polarization measurements alone and extract the Bayesian
evidences and complexities for all the Encyclopaedia Inflationaris models. This
allows us to determine the most probable and simplest BICEP2 inflationary
scenarios. Although this list contains the simplest monomial potentials, it
also includes many other scenarios, suggesting that focusing model building
efforts on large field models only is unjustified at this stage. We demonstrate
that the sets of inflationary models preferred by Planck alone and BICEP2 alone
are almost disjoint, indicating a clear tension between the two data sets. We
address this tension with a Bayesian measure of compatibility between BICEP2
and Planck. We find that for models favored by Planck the two data sets tend to
be incompatible, whereas there is a moderate evidence of compatibility for the
BICEP2 preferred models. As a result, it would be premature to draw any
conclusion on the best Planck models, such as Starobinsky and/or Kahler moduli
inflation. For the subset of scenarios not exhibiting data sets
incompatibility, we update the evidences and complexities using both data sets
together.

We compute the Bayesian evidence and complexity of 193 slowroll singlefield
models of inflation using the Planck 2013 Cosmic Microwave Background data,
with the aim of establishing which models are favoured from a Bayesian
perspective. Our calculations employ a new numerical pipeline interfacing an
inflationary effective likelihood with the slowroll library ASPIC and the
nested sampling algorithm MULTINEST. The models considered represent a complete
and systematic scan of the entire landscape of inflationary scenarios proposed
so far. Our analysis singles out the most probable models (from an Occam's
razor point of view) that are compatible with Planck data, while ruling out
with very strong evidence 34% of the models considered. We identify 26% of the
models that are favoured by the Bayesian evidence, corresponding to 15
different potential shapes. If the Bayesian complexity is included in the
analysis, only 9% of the models are preferred, corresponding to only 9
different potential shapes. These shapes are all of the plateau type.

We explore systematic biases in the identification of dark matter in future
direct detection experiments and compare the reconstructed dark matter
properties when assuming a selfconsistent dark matter distribution function
and the standard Maxwellian velocity distribution. We find that the systematic
bias on the dark matter mass and crosssection determination arising from wrong
assumptions for its distribution function is of order ~1\sigma. A much larger
systematic bias can arise if wrong assumptions are made on the underlying Milky
Way mass model. However, in both cases the bias is substantially mitigated by
marginalizing over galactic model parameters. We additionally show that the
velocity distribution can be reconstructed in an unbiased manner for typical
dark matter parameters. Our results highlight both the robustness of the dark
matter mass and crosssection determination using the standard Maxwellian
velocity distribution and the importance of accounting for astrophysical
uncertainties in a statistically consistent fashion.

We discuss irreducible statistical limitations of future tonscale dark
matter direct detection experiments. We focus in particular on the coverage of
confidence intervals, which quantifies the reliability of the statistical
method used to reconstruct the dark matter parameters, and the bias of the
reconstructed parameters. We study 36 benchmark dark matter models within the
reach of upcoming tonscale experiments. We find that approximate confidence
intervals from a profilelikelihood analysis exactly cover or overcover the
true values of the WIMP parameters, and are hence conservative. We evaluate the
probability that unavoidable statistical fluctuations in the data might lead to
a biased reconstruction of the dark matter parameters, or large uncertainties
on the reconstructed parameter values. We show that this probability can be
surprisingly large, even for benchmark models leading to a large event rate of
order a hundred counts. We find that combining data sets from two different
targets leads to improved coverage properties, as well as a substantial
reduction of statistical bias and uncertainty on the dark matter parameters.

We investigate the question of whether the recent modulation signal claimed
by CoGeNT is best explained by the dark matter (DM) hypothesis from a Bayesian
model comparison perspective. We consider five phenomenological explanations
for the data: no modulation signal, modulation due to DM, modulation due to DM
compatible with the total CoGeNT rate, and a signal coming from other physics
with a free phase but annual period, or with a free phase and a free period. In
each scenario, we assign to the free parameters physically motivated priors. We
find that while the DM models are weakly preferred to the no modulation model,
but when compared to models where the modulation is due to other physics, the
DM hypothesis is favoured with odds ranging from 185:1 to 560:1. This result is
robust even when astrophysical uncertainties are taken into account and the
impact of priors assessed. Interestingly, the odds for the DM model in which
the modulation signal is compatible with the total rate against a DM model in
which this prior is not implemented is only 5:8, in spite of the former's
prediction of a modulation amplitude in the energy range 0.9 to 3.0 keVee that
is significantly smaller than the value observed by CoGeNT. Classical
hypothesis testing also rules out the null hypothesis of no modulation at the
1.6 sigma to 2.3 sigma level, depending on the details of the alternative.
Lastly, we investigate whether anisotropic velocity distributions can help to
mitigate the tension between the CoGeNT total and modulated rates, and find
encouraging results.

The unfortunate case where the Large Hadron Collider (LHC) fails to discover
physics Beyond the Standard Model (BSM) is sometimes referred to as the
"Nightmare scenario" of particle physics. We study the consequences of this
hypothetical scenario for Dark Matter (DM), in the framework of the constrained
Minimal Supersymmetric Standard Model (cMSSM). We evaluate the surviving
regions of the cMSSM parameter space after null searches at the LHC, using
several different LHC configurations, and study the consequences for DM
searches with tonscale direct detectors and the IceCube neutrino telescope. We
demonstrate that tonscale direct detection experiments will be able to
conclusively probe the cMSSM parameter space that would survive null searches
at the LHC with 100fb$^{1}$ of integrated luminosity at 14TeV. We also
demonstrate that IceCube (80 strings plus DeepCore) will be able to probe as
much as 17% of the currently favoured parameter space after 5 years of
observation.

Supersymmetry has been often invoqued as the new physics that might reconcile
the experimental muon magnetic anomaly, a_mu, with the theoretical prediction
(basing the computation of the hadronic contribution on e^+ e^ data). However,
in the context of the CMSSM, the required supersymmetric contributions (which
grow with decreasing supersymmetric masses) are in potential tension with a
possibly large Higgs mass (which requires large stop masses). In the limit of
very large m_h supersymmetry gets decoupled, and the CMSSM must show the same
discrepancy as the SM with a_mu . But it is much less clear for which size of
m_h does the tension start to be unbearable. In this paper, we quantify this
tension with the help of Bayesian techniques. We find that for m_h > 125 GeV
the maximum level of discrepancy given current data (~ 3.3 sigma) is already
achieved. Requiring less than 3 sigma discrepancy, implies m_h < 120 GeV. For a
larger Higgs mass we should give up either the CMSSM model or the computation
of a_mu based on e^+ e^; or accept living with such inconsistency.

We investigate the reconstruction capabilities of Dark Matter mass and
spinindependent crosssection from future tonscale direct detection
experiments using germanium, xenon or argon as targets. Adopting realistic
values for the exposure, energy threshold and resolution of Dark Matter
experiments which will come online within 5 to 10 years, the degree of
complementarity between different targets is quantified. We investigate how the
uncertainty in the astrophysical parameters controlling the local Dark Matter
density and velocity distribution affects the reconstruction. For a 50 GeV
WIMP, astrophysical uncertainties degrade the accuracy in the mass
reconstruction by up to a factor of $\sim 4$ for xenon and germanium, compared
to the case when astrophysical quantities are fixed. However, combination of
argon, germanium and xenon data increases the constraining power by a factor of
$\sim 2$ compared to germanium or xenon alone. We show that future direct
detection experiments can achieve selfcalibration of some astrophysical
parameters, and they will be able to constrain the WIMP mass with only very
weak external astrophysical constraints.

We present recent results aiming at assessing the coverage properties of
Bayesian and frequentist inference methods, as applied to the reconstruction of
supersymmetric parameters from simulated LHC data. We discuss the statistical
challenges of the reconstruction procedure, and highlight the algorithmic
difficulties of obtaining accurate profile likelihood estimates.

We present the first calculation of the Bayesian evidence for different
prototypical single field inflationary scenarios, including representative
classes of small field and large field models. This approach allows us to
compare inflationary models in a welldefined statistical way and to determine
the current "best model of inflation". The calculation is performed numerically
by interfacing the inflationary code FieldInf with MultiNest. We find that
small field models are currently preferred, while large field models having a
selfinteracting potential of power p>4 are strongly disfavoured. The class of
small field models as a whole has posterior odds of approximately 3:1 when
compared with the large field class. The methodology and results presented in
this article are an additional step toward the construction of a full numerical
pipeline to constrain the physics of the early Universe with astrophysical
observations. More accurate data (such as the Planck data) and the techniques
introduced here should allow us to identify conclusively the best inflationary
model.

Dark matter (DM) is currently searched for with a variety of detection
strategies. Accelerator searches are particularly promising, but even if Weakly
Interacting Massive Particles (WIMPs) are found at the Large Hadron Collider
(LHC), it will be difficult to prove that they constitute the bulk of the DM in
the Universe. We show that a significantly better reconstruction of the DM
properties can be obtained with a combined analysis of LHC and direct detection
(DD) data, by making a simple Ansatz on the WIMP local density, i.e. by
assuming that the local density scales with the cosmological relic abundance.
We demonstrate this method in an explicit example in the context of a
24parameter supersymmetric model, with a neutralino LSP in the stau
coannihilation region. Our results show that future tonscale DD experiments
will allow to break degeneracies in the SUSY parameter space and achieve a
significantly better reconstruction of the neutralino composition and its relic
density than with LHC data alone.

We present a new method for determining Weakly Interacting Massive Particle
(WIMP) properties in future tonne scale direct detection experiments which
accounts for uncertainties in the Milky Way (MW) smooth dark matter
distribution. Using synthetic data on the kinematics of MW halo stars matching
present samples from the Sloan Digital Sky Survey, complemented by local escape
velocity constraints, we demonstrate that the local dark matter density can be
constrained to approximately 20% accuracy. For low mass WIMPs, we find that a
factor of two error in the assumed local dark matter density leads to a
severely biased reconstruction of the WIMP spinindependent cross section that
is incorrect at the 15sigma level. We show that this bias may be overcome by
marginalizing over parameters that describe the MW potential, and use this
formalism to project the accuracy attainable on WIMP properties in future 1
tonne Xenon detectors. Our method can be readily applied to different detector
technologies and extended to more detailed MW halo models.

We present a general methodology for determining the gammaray flux from
annihilation of dark matter particles in Milky Way satellite galaxies, focusing
on two promising satellites as examples: Segue 1 and Draco. We use the
SuperBayeS code to explore the bestfitting regions of the Constrained Minimal
Supersymmetric Standard Model (CMSSM) parameter space, and an independent MCMC
analysis of the dark matter halo properties of the satellites using published
radial velocities. We present a formalism for determining the boost from halo
substructure in these galaxies and show that its value depends strongly on the
extrapolation of the concentrationmass (c(M)) relation for CDM subhalos down
to the minimum possible mass. We show that the preferred region for this
minimum halo mass within the CMSSM with neutralino dark matter is ~10^910^6
solar masses. For the boost model where the observed powerlaw c(M) relation is
extrapolated down to the minimum halo mass we find average boosts of about 20,
while the Bullock et al (2001) c(M) model results in boosts of order unity. We
estimate that for the powerlaw c(M) boost model and photon energies greater
than a GeV, the Fermi spacetelescope has about 20% chance of detecting a dark
matter annihilation signal from Draco with signaltonoise greater than 3 after
about 5 years of observation.

There are things we know, things we know we don't know, and then there are
things we don't know we don't know. In this paper we address the latter two
issues in a Bayesian framework, introducing the notion of doubt to quantify the
degree of (dis)belief in a model given observational data in the absence of
explicit alternative models. We demonstrate how a properly calibrated doubt can
lead to model discovery when the true model is unknown.

The application of Bayesian methods in cosmology and astrophysics has
flourished over the past decade, spurred by data sets of increasing size and
complexity. In many respects, Bayesian methods have proven to be vastly
superior to more traditional statistical tools, offering the advantage of
higher efficiency and of a consistent conceptual basis for dealing with the
problem of induction in the presence of uncertainty. This trend is likely to
continue in the future, when the way we collect, manipulate and analyse
observations and compare them with theoretical models will assume an even more
central role in cosmology.
This review is an introduction to Bayesian methods in cosmology and
astrophysics and recent results in the field. I first present Bayesian
probability theory and its conceptual underpinnings, Bayes' Theorem and the
role of priors. I discuss the problem of parameter inference and its general
solution, along with numerical techniques such as Monte Carlo Markov Chain
methods. I then review the theory and application of Bayesian model comparison,
discussing the notions of Bayesian evidence and effective model complexity, and
how to compute and interpret those quantities. Recent developments in
cosmological parameter extraction and Bayesian cosmological model building are
summarized, highlighting the challenges that lie ahead.

We consider an analytic model of cosmic star formation which incorporates
supernova feedback, gas accretion and enriched outflows, reproducing the
history of cosmic star formation, metallicity, supernovae type II rates and the
fraction of baryons allocated to structures. We present a new statistical
treatment of the available observational data on the star formation rate and
metallicity that accounts for the presence of possible systematics. We then
employ a Bayesian Markov Chain Monte Carlo method to compare the predictions of
our model with observations and derive constraints on the 7 free parameters of
the model. We find that the dust correction scheme one chooses to adopt for the
star formation data is critical in determining which scenario is favoured
between a hierarchical star formation model, where star formation is prolonged
by accretion, infall and merging, and a monolithic scenario, where star
formation is rapid and efficient. We distinguish between these modes by
defining a characteristic minimum mass, M > 10^{11} M_solar, in our fiducial
model, for early type galaxies where star formation occurs efficiently. Our
results indicate that the hierarchical star formation model can achieve better
agreement with the data, but that this requires a high efficiency of
supernovadriven outflows. In a monolithic model, our analysis points to the
need for a mechanism that drives metalpoor winds, perhaps in the form of
supermassive black holeinduced outflows. Furthermore, the relative absence of
star formation beyond z ~ 5 in the monolithic scenario requires an alternative
mechanism to dwarf galaxies for reionizing the universe at z ~ 11, as required
by observations of the microwave background. While the monolithic scenario is
less favoured in terms of its qualityoffit, it cannot yet be excluded.

Bayesian model selection provides a formal method of determining the level of
support for new parameters in a model. However, if there is not a specific
enough underlying physical motivation for the new parameters it can be hard to
assign them meaningful priors, an essential ingredient of Bayesian model
selection. Here we look at methods maximizing the prior so as to work out what
is the maximum support the data could give for the new parameters. If the
maximum support is not high enough then one can confidently conclude that the
new parameters are unnecessary without needing to worry that some other prior
may make them significant. We discuss a computationally efficient means of
doing this which involves mapping pvalues onto upper bounds of the Bayes
factor (or odds) for the new parameters. A pvalue of 0.05 ($1.96\sigma$)
corresponds to odds less than or equal to 5:2 which is below the `weak' support
at best threshold. A pvalue of 0.0003 ($3.6\sigma$) corresponds to odds of
less than or equal to 150:1 which is the `strong' support at best threshold.
Applying this method we find that the odds on the scalar spectral index being
different from one are 49:1 at best. We also find that the odds that there is
primordial hemispherical asymmetry in the cosmic microwave background are 9:1
at best.

I present a new procedure to forecast the Bayes factor of a future
observation by computing the Predictive Posterior Odds Distribution (PPOD).
This can assess the power of future experiments to answer model selection
questions and the probability of the outcome, and can be helpful in the context
of experiment design.
As an illustration, I consider a central quantity for our understanding of
the cosmological concordance model, namely the scalar spectral index of
primordial perturbations, n_S. I show that the Planck satellite has over 90%
probability of gathering strong evidence against n_S = 1, thus conclusively
disproving a scaleinvariant spectrum. This result is robust with respect to a
wide range of choices for the prior on n_S.