We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

We present a new estimate of the mass of the Milky Way, inferred via a Bayesian approach by making use of tracers of the circular velocity in the disk plane and stars in the stellar halo, as from the publicly available {\tt galkin} compilation. We use the rotation curve method to determine the dark matter distribution and total mass under different assumptions for the dark matter profile, while the total stellar mass is constrained by surface stellar density and microlensing measurements. We also include uncertainties on the baryonic morphology via Bayesian model averaging, thus converting a potential source of systematic error into a more manageable statistical uncertainty. We evaluate the robustness of our result against various possible systematics, including rotation curve data selection, uncertainty on the Sun's velocity $V_0$, dependence on the dark matter profile assumptions, and choice of priors. We find the Milky Way's dark matter virial mass to be $\log_{10}M_{200}^{\rm DM}/ {\rm M_\odot} = 11.92^{+0.06}_{-0.05}{\rm(stat)}\pm{0.28}\pm0.27{\rm(syst)}$ ($M_{200}^{\rm DM}=8.3^{+1.2}_{-0.9}{\rm(stat)}\times10^{11}\,{\rm M_\odot}$). We also apply our framework to Gaia DR2 rotation curve data and find good statistical agreement with the above results.

We present a new solution to the problem of classifying Type Ia supernovae from their light curves alone given a spectroscopically confirmed but biased training set, circumventing the need to obtain an observationally expensive unbiased training set. We use Gaussian processes (GPs) to model the supernovae's (SN) light curves, and demonstrate that the choice of covariance function has only a small influence on the GPs ability to accurately classify SNe. We extend and improve the approach of Richards et al (2012} -- a diffusion map combined with a random forest classifier -- to deal specifically with the case of biassed training sets. We propose a novel method, called STACCATO (SynThetically Augmented Light Curve ClassificATiOn') that synthetically augments a biased training set by generating additional training data from the fitted GPs. Key to the success of the method is the partitioning of the observations into subgroups based on their propensity score of being included in the training set. Using simulated light curve data, we show that STACCATO increases performance, as measured by the area under the Receiver Operating Characteristic curve (AUC), from 0.93 to 0.96, close to the AUC of 0.977 obtained using the 'gold standard' of an unbiased training set and significantly improving on the previous best result of 0.88. STACCATO also increases the true positive rate for SNIa classification by up to a factor of 50 for high-redshift/low brightness SNe.

Dwarf spheroidal galaxies that form in halo substructures provide stringent constraints on dark matter annihilation. Many ultrafaint dwarfs discovered with modern surveys contribute significantly to these constraints. At present, because of the lack of abundant stellar kinematic data for the ultrafaints, non-informative prior assumptions are usually made for the parameters of the density profiles. Based on semi-analytic models of dark matter subhalos and their connection to satellite galaxies, we present more informative and realistic satellite priors. We show that our satellite priors lead to constraints on the annihilation rate that are between a factor of 2 and a factor of 7 weaker than under non-informative priors. As a result, the thermal relic cross section can at best only be excluded (with 95% probability) for dark matter masses of $\lesssim 40$ GeV from dwarf spheroidal data, assuming annihilation into $b\bar{b}$.

We search for a dark matter signal in 11 years of Fermi-LAT gamma-ray data from 27 Milky Way dwarf spheroidal galaxies with spectroscopically measured $J$-factors. Our analysis includes uncertainties in $J$-factors and background normalisations and compares results from a Bayesian and a frequentist perspective. We revisit the dwarf spheroidal galaxy Reticulum II, confirming that the purported gamma-ray excess seen in Pass 7 data is much weaker in Pass 8, independently of the statistical approach adopted. We introduce for the first time posterior predictive distributions to quantify the probability of a dark matter detection from another dwarf galaxy given a tentative excess. A global analysis including all 27 dwarfs shows no indication for a signal in nine annihilation channels. We present stringent new Bayesian and frequentist upper limits on the dark matter cross section as a function of dark matter mass. The best-fit dark matter parameters associated with the Galactic Centre excess are excluded by at least 95% confidence level/posterior probability in the frequentist/Bayesian framework in all cases. However, from a Bayesian model comparison perspective, dark matter annihilation within the dwarfs is not strongly disfavoured compared to a background-only model. These results constitute the highest exposure analysis on the most complete sample of dwarfs to date. Posterior samples and likelihood maps from this study are publicly available.

We develop a novel Bayesian methodology aimed at reliably and precisely inferring the distribution of dark matter within the Milky Way using rotation curve data. We identify a subset of the available rotation curve tracers that are mutually consistent with each other, thus eliminating data sets that might suffer from systematic bias. We investigate three different models for the mass distribution of the luminous (baryonic) component that bracket the range of likely morphologies. We demonstrate the statistical performance of our method on simulated data in terms of coverage, fractional distance, and mean squared error. Applying it to Milky Way data we measure the local dark matter density at the solar circle $\rho_0$ to be $\rho_0 = 0.43\pm 0.02(\rm{stat})\pm0.01(\rm{sys})$ GeV/cm$^3$, with an accuracy $\sim$ 6%. This result is robust to the assumed baryonic morphology. The scale radius and inner slope of the dark matter profile are degenerate and cannot be individually determined with high accuracy. We show that these results are robust to several possible residual systematic errors in the rotation curve data.

We use multi-band imagery data from the Sloan Digital Sky Survey (SDSS) to measure projected distances of 302 supernova type Ia (SNIa) from the centre of their host galaxies, normalized to the galaxy's brightness scale length, with a Bayesian approach. We test the hypothesis that SNIas further away from the centre of their host galaxy are less subject to dust contamination (as the dust column density in their environment is smaller) and/or come from a more homogeneous environment. Using the Mann-Whitney U test, we find a statistically significant difference in the observed colour correction distribution between SNIas that are near and those that are far from the centre of their host. The local p-value is 3 x 10^{-3}, which is significant at the 5 per cent level after look-elsewhere effect correction. We estimate the residual scatter of the two subgroups to be 0.073 +/- 0.018 for the far SNIas, compared to 0.114 +/- 0.009 for the near SNIas -- an improvement of 30 per cent, albeit with a low statistical significance of 2sigma. This confirms the importance of host galaxy properties in correctly interpreting SNIa observations for cosmological inference.

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 non-trivial 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 point-like 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 WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b 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 WASP-12b. The superstellar carbon-to-oxygen 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 WASP-12b in a solar-composition 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 disk-free mechanism, and are inconsistent with both in-situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. (2017). We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.

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, research-level 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 p-value 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.

DARk matter WImp search with liquid xenoN (DARWIN) will be an experiment for the direct detection of dark matter using a multi-ton liquid xenon time projection chamber at its core. Its primary goal will be to explore the experimentally accessible parameter space for Weakly Interacting Massive Particles (WIMPs) in a wide mass-range, until neutrino interactions with the target become an irreducible background. The prompt scintillation light and the charge signals induced by particle interactions in the xenon will be observed by VUV sensitive, ultra-low background photosensors. Besides its excellent sensitivity to WIMPs above a mass of 5 GeV/c2, such a detector with its large mass, low-energy threshold and ultra-low background level will also be sensitive to other rare interactions. It will search for solar axions, galactic axion-like particles and the neutrinoless double-beta decay of 136-Xe, as well as measure the low-energy solar neutrino flux with <1% precision, observe coherent neutrino-nucleus interactions, and detect galactic supernovae. We present the concept of the DARWIN detector and discuss its physics reach, the main sources of backgrounds and the ongoing detector design and R&D efforts.

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 low-z SNIa, observed by the Carnegie Supernova Project (CSP) in the optical and near infra-red (NIR). We study the time of the second peak in the NIR band due to re-brightening, 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 results obtained by applying our BAyesian HierArchical Modeling for the Analysis of Supernova cosmology (BAHAMAS) software package to the 740 spectroscopically confirmed supernovae type Ia (SNIa) from the "Joint Light-curve Analysis" (JLA) dataset. We simultaneously determine cosmological parameters and standardization parameters, including host galaxy mass corrections, residual scatter and object-by-object intrinsic magnitudes. Combining JLA and Planck Cosmic Microwave Background data, we find significant discrepancies in cosmological parameter constraints with respect to the standard analysis: we find Omega_M = 0.399+/-0.027, 2.8\sigma\ higher than previously reported and w = -0.910+/-0.045, 1.6\sigma\ higher than the standard analysis. We determine the residual scatter to be sigma_res = 0.104+/-0.005. We confirm (at the 95% probability level) the existence of two sub-populations segregated by host galaxy mass, separated at log_{10}(M/M_solar) = 10, differing in mean intrinsic magnitude by 0.055+/-0.022 mag, lower than previously reported. Cosmological parameter constraints are however unaffected by inclusion of host galaxy mass corrections. We find ~4\sigma\ evidence for a sharp drop in the value of the color correction parameter, beta(z), at a redshift z_trans = 0.662+/-0.055. We rule out some possible explanations for this behaviour, which remains unexplained.

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 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 p-wave annihilation, and also do not contribute to the spin-independent 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 non-detection 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 2015-2025 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 B-modes 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 tensor-to-scalar ratio. Using only those four bandpowers, we carry out a complete analysis in the cosmological and inflationary slow-roll 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 slow-roll single-field 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 slow-roll 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 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.

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 self-consistent dark matter distribution function and the standard Maxwellian velocity distribution. We find that the systematic bias on the dark matter mass and cross-section 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 cross-section determination using the standard Maxwellian velocity distribution and the importance of accounting for astrophysical uncertainties in a statistically consistent fashion.

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 discuss irreducible statistical limitations of future ton-scale 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 ton-scale experiments. We find that approximate confidence intervals from a profile-likelihood analysis exactly cover or over-cover 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.

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.

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 ton-scale direct detectors and the IceCube neutrino telescope. We demonstrate that ton-scale 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.

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.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We present a new method based on a Bayesian hierarchical model to extract constraints on cosmological parameters from SNIa data obtained with the SALT-II lightcurve fitter. We demonstrate with simulated data sets that our method delivers tighter statistical constraints on the cosmological parameters over 90% of the time, that it reduces statistical bias typically by a factor ~ 2-3 and that it has better coverage properties than the usual chi-squared approach. As a further benefit, a full posterior probability distribution for the dispersion of the intrinsic magnitude of SNe is obtained. We apply this method to recent SNIa data, and by combining them with CMB and BAO data we obtain Omega_m=0.28 +/- 0.02, Omega_Lambda=0.73 +/- 0.01 (assuming w=-1) and Omega_m=0.28 +/- 0.01, w=-0.90 +/- 0.05 (assuming flatness; statistical uncertainties only). We constrain the intrinsic dispersion of the B-band magnitude of the SNIa population, obtaining sigma_mu^int = 0.13 +/- 0.01 [mag]. Applications to systematic uncertainties will be discussed in a forthcoming paper.

We investigate the reconstruction capabilities of Dark Matter mass and spin-independent cross-section from future ton-scale 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 self-calibration of some astrophysical parameters, and they will be able to constrain the WIMP mass with only very weak external astrophysical constraints.

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 well-defined 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 self-interacting 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.

We extend the Figure of Merit formalism usually adopted to quantify the statistical performance of future dark energy probes to assess the robustness of a future mission to plausible systematic bias. We introduce a new robustness Figure of Merit which can be computed in the Fisher Matrix formalism given arbitrary systematic biases in the observable quantities. We argue that robustness to systematics is an important new quantity that should be taken into account when optimizing future surveys. We illustrate our formalism with toy examples, and apply it to future type Ia supernova (SNIa) and baryonic acoustic oscillation (BAO) surveys. For the simplified systematic biases that we consider, we find that SNIa are a somewhat more robust probe of dark energy parameters than the BAO. We trace this back to a geometrical alignement of systematic bias direction with statistical degeneracy directions in the dark energy parameter space.

Bayesian model averaging is a procedure to obtain parameter constraints that account for the uncertainty about the correct cosmological model. We use recent cosmological observations and Bayesian model averaging to derive tight limits on the curvature parameter, as well as robust lower bounds on the curvature radius of the Universe and its minimum size, while allowing for the possibility of an evolving dark energy component. Because flat models are favoured by Bayesian model selection, we find that model-averaged constraints on the curvature and size of the Universe can be considerably stronger than non model-averaged ones. For the most conservative prior choice (based on inflationary considerations), our procedure improves on non model-averaged constraints on the curvature by a factor of ~ 2. The curvature scale of the Universe is conservatively constrained to be R_c > 42 Gpc (99%), corresponding to a lower limit to the number of Hubble spheres in the Universe N_U > 251 (99%).

We present a general Bayesian formalism for the definition of Figures of Merit (FoMs) quantifying the scientific return of a future experiment. We introduce two new FoMs for future experiments based on their model selection capabilities, called the decisiveness of the experiment and the expected strength of evidence. We illustrate these by considering dark energy probes, and compare the relative merits of stage II, III and IV dark energy probes. We find that probes based on supernovae and on weak lensing perform rather better on model selection tasks than is indicated by their Fisher matrix FoM as defined by the Dark Energy Task Force. We argue that our ability to optimize future experiments for dark energy model selection goals is limited by our current uncertainty over the models and their parameters, which is ignored in the usual Fisher matrix forecasts. Our approach gives a more realistic assessment of the capabilities of future probes and can be applied in a variety of situations.

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.

In theories with Universal Extra-Dimensions (UED), the gamma_1 particle, first excited state of the hypercharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in order to assess its detectability at accelerators and with DM experiments. We derive in particular the most probable range of mass and scattering cross sections off nucleons, keeping into account cosmological and electroweak precision constraints. The consequences for the detectability of the gamma_1 with direct and indirect experiments are dramatic. The spin-independent cross section probability distribution peaks at ~ 10^{-11} pb, i.e. below the sensitivity of ton-scale experiments. The spin-dependent cross-section drives the predicted neutrino flux from the center of the Sun below the reach of present and upcoming experiments. The only strategy that remains open appears to be direct detection with ton-scale experiments sensitive to spin-dependent cross-sections. On the other hand, the LHC with 1 1/fb of data should be able to probe the current best-fit UED parameters.

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 24-parameter supersymmetric model, with a neutralino LSP in the stau co-annihilation region. Our results show that future ton-scale 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.

While Bayesian model selection is a useful tool to discriminate between competing cosmological models, it only gives a relative rather than an absolute measure of how good a model is. Bayesian doubt introduces an unknown benchmark model against which the known models are compared, thereby obtaining an absolute measure of model performance in a Bayesian framework. We apply this new methodology to the problem of the dark energy equation of state, comparing an absolute upper bound on the Bayesian evidence for a presently unknown dark energy model against a collection of known models including a flat LambdaCDM scenario. We find a strong absolute upper bound to the Bayes factor B between the unknown model and LambdaCDM, giving B < 3. The posterior probability for doubt is found to be less than 6% (with a 1% prior doubt) while the probability for LambdaCDM rises from an initial 25% to just over 50% in light of the data. We conclude that LambdaCDM remains a sufficient phenomenological description of currently available observations and that there is little statistical room for model improvement.

In arXiv:0911.2150, Rutger van Haasteren seeks to criticize the nested sampling algorithm for Bayesian data analysis in general and its MultiNest implementation in particular. He introduces a new method for evidence evaluation based on the idea of Voronoi tessellation and requiring samples from the posterior distribution obtained through MCMC based methods. He compares its accuracy and efficiency with MultiNest, concluding that it outperforms MultiNest in several cases. This comparison is completely unfair since the proposed method can not perform the complete Bayesian data analysis including posterior exploration and evidence evaluation on its own while MultiNest allows one to perform Bayesian data analysis end to end. Furthermore, their criticism of nested sampling (and in turn MultiNest) is based on a few conceptual misunderstandings of the algorithm. Here we seek to set the record straight.

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 spin-independent cross section that is incorrect at the 15-sigma 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.

Cosmologists will soon be in a unique position. Observational noise will gradually be replaced by cosmic variance as the dominant source of uncertainty in an increasing number of observations. We reflect on the ramifications for the discovery and verification of new models. If there are features in the full data set that call for a new model, there will be no subsequent observations to test that model's predictions. We give specific examples of the problem by discussing the pitfalls of model discovery by prior adjustment in the context of dark energy models and inflationary theories. We show how the gradual release of data can mitigate this difficulty, allowing anomalies to be identified, and new models to be proposed and tested. We advocate that observers plan for the frugal release of data from future cosmic variance limited observations.

We study in detail the ability of the nominal configuration of the IceCube neutrino telescope (with 80 strings) to probe the parameter space of the Constrained MSSM (CMSSM) favoured by current collider and cosmological data. Adopting conservative assumptions about the galactic halo model and the expected experiment performance, we find that IceCube has a probability between 2% and 12% of achieving a 5sigma detection of dark matter annihilation in the Sun, depending on the choice of priors for the scalar and gaugino masses and on the astrophysical assumptions. We identify the most important annihilation channels in the CMSSM parameter space favoured by current constraints, and we demonstrate that assuming that the signal is dominated by a single annihilation channel canlead to large systematic errors in the inferred WIMP annihilation cross section. We demonstrate that ~ 66% of the CMSSM parameter space violates the equilibrium condition between capture and annihilation in the center of the Sun. By cross-correlating our predictions with direct detection methods, we conclude that if IceCube does detect a neutrino flux from the Sun at high significance while direct detection experiments do not find a signal above a spin-independent cross section sigma_SI^p larger than 5x10^{-9} pb, the CMSSM will be strongly disfavoured, given standard astrophysical assumptions for the WIMP distribution. This result is robust with respect to a change of priors. We argue that the proposed low-energy DeepCore extension of IceCube will be an ideal instrument to focus on relevant CMSSM areas of parameter space.

We present a general methodology for determining the gamma-ray 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 best-fitting 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 concentration-mass (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^-9-10^-6 solar masses. For the boost model where the observed power-law 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 power-law c(M) boost model and photon energies greater than a GeV, the Fermi space-telescope has about 20% chance of detecting a dark matter annihilation signal from Draco with signal-to-noise greater than 3 after about 5 years of observation.

The question of determining the spatial geometry of the Universe is of greater relevance than ever, as precision cosmology promises to verify inflationary predictions about the curvature of the Universe. We revisit the question of what can be learnt about the spatial geometry of the Universe from the perspective of a three-way Bayesian model comparison. We show that, given current data, the probability that the Universe is spatially infinite lies between 67% and 98%, depending on the choice of priors. For the strongest prior choice, we find odds of order 50:1 (200:1) in favour of a flat Universe when compared with a closed (open) model. We also report a robust, prior-independent lower limit to the number of Hubble spheres in the Universe, N_U > 5 (at 99% confidence). We forecast the accuracy with which future CMB and BAO observations will be able to constrain curvature, finding that a cosmic variance limited CMB experiment together with an SKA-like BAO observation will constrain curvature with a precision of about sigma ~ 4.5x10^{-4}. We demonstrate that the risk of 'model confusion' (i.e., wrongly favouring a flat Universe in the presence of curvature) is much larger than might be assumed from parameter errors forecasts for future probes. We argue that a 5-sigma detection threshold guarantees a confusion- and ambiguity-free model selection. Together with inflationary arguments, this implies that the geometry of the Universe is not knowable if the value of the curvature parameter is below |Omega_curvature| ~ 10^{-4}, a bound one order of magnitude larger than the size of curvature perturbations, ~ 10^{-5}. [abridged]

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.

In the framework of the Constrained MSSM we derive the most probable ranges of the diffuse gamma radiation flux from the direction of the Galactic center and of the positron flux from the Galactic halo due to neutralino dark matter annihilation. We find that, for a given halo model, and assuming flat priors, the 68% probability range of the integrated gamma-ray flux spans about one order of magnitude, while the 95% probability range can be much larger and extend over four orders of magnitude (even exceeding five for a tiny region at small neutralino mass). The detectability of the signal by GLAST depending primarily on the cuspiness of the halo profile. The positron flux, on the other hand, appears to be too small to be detectable by PAMELA, unless the boost factor is at least of order ten and/or the halo profile is extremely cuspy. We also briefly discuss the sensitivity of our results to the choice of priors.

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 re-examine the properties of the Constrained MSSM in light of updated constraints, paying particular attention to the impact of the recent substantial shift in the Standard Model prediction for BR(B to X_s gamma). With the help of a Markov Chain Monte Carlo scanning technique, we vary all relevant parameters simultaneously and derive Bayesian posterior probability maps. We find that the case of \mu>0 remains favored, and that for \mu<0 it is considerably more difficult to find a good global fit to current constraints. In both cases we find a strong preference for a focus point region. This leads to improved prospects for detecting neutralino dark matter in direct searches, while superpartner searches at the LHC become more problematic, especially when \mu<0. In contrast, prospects for exploring the whole mass range of the lightest Higgs boson at the Tevatron and the LHC remain very good, which should, along with dark matter searches, allow one to gain access to the otherwise experimentally challenging focus point region. An alternative measure of the mean quality-of-fit which we also employ implies that present data are not yet constraining enough to draw more definite conclusions. We also comment on the dependence of our results on the choice of priors and on some other assumptions.

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 supernova-driven outflows. In a monolithic model, our analysis points to the need for a mechanism that drives metal-poor winds, perhaps in the form of supermassive black hole-induced 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 quality-of-fit, it cannot yet be excluded.

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.

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 p-values onto upper bounds of the Bayes factor (or odds) for the new parameters. A p-value 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 p-value 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.

We present a Bayesian technique based on a maximum entropy method to reconstruct the dark energy equation of state $w(z)$ in a non--parametric way. This MaxEnt technique allows to incorporate relevant prior information while adjusting the degree of smoothing of the reconstruction in response to the structure present in the data. After demonstrating the method on synthetic data, we apply it to current cosmological data, separately analysing type Ia supernovae measurement from the HST/GOODS program and the first year Supernovae Legacy Survey (SNLS), complemented by cosmic microwave background and baryonic acoustic oscillations data. We find that the SNLS data are compatible with $w(z) = -1$ at all redshifts $0 \leq z \lsim 1100$, with errorbars of order 20% for the most constraining choice of priors. The HST/GOODS data exhibit a slight (about $1\sigma$ significance) preference for $w>-1$ at $z\sim 0.5$ and a drift towards $w>-1$ at larger redshifts, which however is not robust with respect to changes in our prior specifications. We employ both a constant equation of state prior model and a slowly varying $w(z)$ and find that our conclusions are only mildly dependent on this choice at high redshifts. Our method highlights the danger of employing parametric fits for the unknown equation of state, that can potentially miss or underestimate real structure in the data.

Bayesian model selection is a tool to decide whether the introduction of a new parameter is warranted by data. I argue that the usual sampling statistic significance tests for a null hypothesis can be misleading, since they do not take into account the information gained through the data, when updating the prior distribution to the posterior. On the contrary, Bayesian model selection offers a quantitative implementation of Occam's razor. I introduce the Savage-Dickey density ratio, a computationally quick method to determine the Bayes factor of two nested models and hence perform model selection. As an illustration, I consider three key parameters for our understanding of the cosmological concordance model. By using WMAP 3-year data complemented by other cosmological measurements, I show that a non-scale invariant spectral index of perturbations is favoured for any sensible choice of prior. It is also found that a flat Universe is favoured with odds of 29:1 over non--flat models, and that there is strong evidence against a CDM isocurvature component to the initial conditions which is totally (anti)correlated with the adiabatic mode (odds of about 2000:1), but that this is strongly dependent on the prior adopted. These results are contrasted with the analysis of WMAP 1-year data, which were not informative enough to allow a conclusion as to the status of the spectral index. In a companion paper, a new technique to forecast the Bayes factor of a future observation is presented.

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 scale-invariant spectrum. This result is robust with respect to a wide range of choices for the prior on n_S.

In astro-ph/0702542v2, Linder and Miquel seek to criticize the use of Bayesian model selection for data analysis and for survey forecasting and design. Their discussion is based on three serious misunderstandings of the conceptual underpinnings and application of model-level Bayesian inference, which invalidate all their main conclusions. Their paper includes numerous further inaccuracies, including an erroneous calculation of the Bayesian Information Criterion. Here we seek to set the record straight.

I revisit the question of the adiabaticity of initial conditions for cosmological perturbations in view of the 3-year WMAP data. I focus on the simplest alternative to pure adiabatic conditions, namely a superposition of the adiabatic mode and one of the 3 possible isocurvature modes, with the same spectral index as the adiabatic component. I discuss findings in terms of posterior bounds on the isocurvature fraction and Bayesian model selection. The Bayes factor (models likelihood ratio) and the effective Bayesian complexity are computed for several prior ranges for the isocurvature content. I find that the CDM isocurvature fraction is now constrained to be less than about 10%, while the fraction in either the neutrino entropy or velocity mode is below about 20%. Model comparison strongly disfavours mixed models that allow for isocurvature fractions larger than unity, while current data do not allow to distinguish between a purely adiabatic model and models with a moderate (ie, below about 10%) isocurvature contribution. The conclusion is that purely adiabatic conditions are strongly favoured from a model selection perspective. This is expected to apply in even stronger terms to more complicated superpositions of isocurvature contributions.

We revisit anthropic arguments purporting to explain the measured value of the cosmological constant. We argue that different ways of assigning probabilities to candidate universes lead to totally different anthropic predictions. As an explicit example, we show that weighting different universes by the total number of possible observations leads to an extremely small probability for observing a value of Lambda equal to or greater than what we now measure. We conclude that anthropic reasoning within the framework of probability as frequency is ill-defined and that in the absence of a fundamental motivation for selecting one weighting scheme over another the anthropic principle cannot be used to explain the value of Lambda, nor, likely, any other physical parameters.

Selection effects in cosmology are often invoked to "explain" why some of the fundamental constant of Nature, and in particular the cosmological constant, take on the value they do in our Universe. We briefly review this probabilistic "anthropic reasoning" and we argue that different (equally plausible) ways of assigning probabilities to candidate universes lead to totally different anthropic predictions, presenting an explicit example based on the total number of possible observations observers can carry out. We conclude that in absence of a fundamental motivation for selecting one weighting scheme over another the anthropic principle cannot be used to explain the value of Lambda. .

We perform a comprehensive exploration of the Constrained MSSM parameters using a Markov Chain Monte Carlo technique and a Bayesian analysis. We compute superpartner masses and other collider observables as well as a cold dark matter abundance, and compare them with experimental data. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant SM parameters. We delineate probability distributions of the CMSSM parameters, the collider and cosmological observables as well as a dark matter direct detection cross section. The 68% probability intervals of the CMSSM parameters are: 0.52 TeV < m_{1/2} < 1.26 TeV, m_0 <2.10 TeV, -0.34 TeV < A_0 < 2.41 TeV and 38.5< tan(beta) <54.6. Generally, large fractions of high probability ranges of the superpartner masses will be probed at the LHC. We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We further expose a number of correlations among the observables, in particular between BR(B_s \to \mu^+ \mu^-) and BR({\bar B}\to X_s\gamma) or sigma_p^{SI}. Once SUSY is discovered, this and other correlations may prove helpful in distinguishing the CMSSM from other supersymmetric models. The robustness of our results is investigated in terms of the assumed ranges of CMSSM parameters and the effect of the (g-2)_mu anomaly which shows some tension with the other observables. We find that the results for m_0, and the observables which strongly depend on it, are sensitive to our assumptions, while our conclusions for the other variables are robust.

We outline the WIMP dark matter parameter space in the Constrained MSSM by performing a comprehensive statistical analysis that compares with experimental data predicted superpartner masses and other collider observables as well as a cold dark matter abundance. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant Standard Model parameters. We present high-probability regions for neutralino dark matter direct detection cross section, and we find that 10^-10 pb < sigma_p^SI < 10^-8 for direct WIMP detection (with details slightly dependent on the assumptions made). We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We conclude that most of the 95% probability region for the cross section will be explored by future one-tonne detectors, that will therefore cover most of the currently favoured region of parameter space.

I review the observational prospects to constrain the equation of state parameter of dark energy and I discuss the potential of future imaging and redshift surveys. Bayesian model selection is used to address the question of the level of accuracy on the equation of state parameter that is required before explanations alternative to a cosmological constant become very implausible. I discuss results in the prediction space of dark energy models. If no significant departure from w=-1 is detected, a precision on w of order 1% will translate into strong evidence against fluid-like dark energy, while decisive evidence will require a precision of order 10^-3.

We examine the prospects for the next generation of surveys aimed at elucidating the nature of dark energy. We review the methods that can be used to determine the redshift evolution of the dark energy equation of state parameter w, highlighting their respective strengths and potential weaknesses. All of the attractive methods require surveys covering more than 5-10,000 square degrees of the sky. We examine the accuracy that each method is likely to deliver within a decade, and discuss the difficulties arising from systematic uncertainties associated with the techniques. We conclude that the proposed photometric and redshift surveys have the potential of delivering measurements of w with percent accuracy at several redshifts out to z ~ 3 . Of particular interest will be the combination of weak lensing and baryonic acoustic oscillations measurements. This exquisite precision is likely to have a fundamental impact on our understanding of the nature of dark energy, providing the necessary guidance for its theoretical explanation.

We introduce a statistical measure of the effective model complexity, called the Bayesian complexity. We demonstrate that the Bayesian complexity can be used to assess how many effective parameters a set of data can support and that it is a useful complement to the model likelihood (the evidence) in model selection questions. We apply this approach to recent measurements of cosmic microwave background anisotropies combined with the Hubble Space Telescope measurement of the Hubble parameter. Using mildly non-informative priors, we show how the 3-year WMAP data improves on the first-year data by being able to measure both the spectral index and the reionization epoch at the same time. We also find that a non-zero curvature is strongly disfavored. We conclude that although current data could constrain at least seven effective parameters, only six of them are required in a scheme based on the Lambda-CDM concordance cosmology.

We demonstrate that combining Cosmic Microwave Background anisotropy measurements from the 1st year WMAP observations with clustering data from the SLOAN galaxy redshift survey yields an indication for primordial anisotropies in the cosmological Neutrino Background.

In the first part of this work, I review the theoretical framework of cosmological perturbation theory necessary to understand the generation and evolution of cosmic microwave background (CMB) anisotropies. Using analytical and numerical techniques, in the second part I describe the impact on the CMB power spectra of the standard cosmological parameters (such as the matter-energy budget of the Universe, its curvature, the amplitude and spectral properties of the primordial fluctuations, etc.). I introduce the most general type of initial conditions for the primordial perturbations, deriving a new analytical approximation for the neutrino isocurvature modes. In the third part, I discuss the issue of extracting constraints on the parameters of interest from the recent, high-quality CMB measurements, presenting the relevant statistical tools and focusing on the Fisher matrix analysis as a technique to produce reliable forecasts for the performance of future observations. I then apply those tools to the study of several possible extensions of the (currently) standard Lambda CDM model: the presence of extra relativistic particles, possible time variations of the fine structure constant, and the value of the primordial Helium mass fraction. I also use the CMB as a tool to study the very early Universe, via its dependence on the type of initial conditions: I relax the assumption of purely adiabatic initial conditions and discuss the observational consequences and constraints on the presence of general isocurvature modes.

We consider a supersymmetric grand unified model which naturally solves the strong CP and mu problems via a Peccei-Quinn symmetry and leads to the standard realization of hybrid inflation. We show that the Peccei-Quinn field of this model can act as curvaton. In contrast to the standard curvaton hypothesis, both the inflaton and the curvaton contribute to the total curvature perturbation. The model predicts an isocurvature perturbation too which has mixed correlation with the adiabatic one. The cold dark matter of the universe is mostly constituted by axions plus a small amount of lightest sparticles. The predictions of the model are confronted with the Wilkinson microwave anisotropy probe and other cosmic microwave background radiation data. We analyze two representative choices of parameters and derive bounds on the curvaton contribution to the adiabatic perturbation. We find that, for the choice which provides the best fitting of the data, the curvaton contribution to the adiabatic amplitude must be smaller than about 67% (at 95% confidence level). The best-fit power spectra are dominated by the adiabatic part of the inflaton contribution. We use Bayesian model comparison to show that this choice of parameters is disfavored with respect to the pure inflaton scale-invariant case with odds of 50 to 1. For the second choice of parameters, the adiabatic mode is dominated by the curvaton, but this choice is strongly disfavored relative to the pure inflaton scale-invariant case (with odds of 10^7 to 1). We conclude that in the present framework the perturbations must be dominated by the adiabatic component from the inflaton.

We introduce the concepts of adiabatic (curvature) and isocurvature (entropy) cosmological perturbations and present their relevance for parameter estimation from cosmic microwave background anisotropies data. We emphasize that, while present-day data are in excellent agreement with pure adiabaticity, subdominant isocurvature contributions cannot be ruled out. We discuss model independent constraints on the isocurvature contribution. Finally, we argue that the Planck satellite will be able to do precision cosmology even if the assumption of adiabaticity is relaxed.

We present up-to-date constraints on the value of the fine-structure constant at the epoch of decoupling from the recent observations made by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. In the framework of models we considered, a positive (negative) variation of the value of $\alpha$ at decoupling with respect to the present-day value is now bounded to be smaller than 2% (6%) at 95% confidence level. We point out that the existence of an early reionization epoch as suggested by the above measurements will, when more accurate cosmic microwave background polarization data is available, lead to considerably tighter constraints. The so-called `reionization bump', in particular, will be extremely useful for this purpose. We find that constraints on $\alpha$ can in principle only be as small as $0.1 %$ using CMB data alone--tighter constraints will require further (non-CMB) priors.

We consider for the first time the ability of present-day cosmic microwave background (CMB) anisotropies data to determine the primordial helium mass fraction, Y_p. We find that CMB data alone gives the confidence interval 0.160 < Y_p < 0.501 (at 68% c.l.). We analyse the impact on the baryon abundance as measured by CMB and discuss the implications for big bang nucleosynthesis. We identify and discuss correlations between the helium mass fraction and both the redshift of reionization and the spectral index. We forecast the precision of future CMB observations, and find that Planck alone will measure Y_p with error-bars of 5%. We point out that the uncertainty in the determination of the helium fraction will have to be taken into account in order to correctly estimate the baryon density from Planck-quality CMB data.

We present a detailed analysis of present and future Cosmic Microwave Background constraints of the value of the fine-structure constant, $\alpha$. We carry out a more detailed analysis of the WMAP first-year data, deriving state-of-the-art constraints on $\alpha$ and discussing various other issues, such as the possible hints for the running of the spectral index. We find, at 95% C.L. that $0.95 < \alpha_{\text{dec}} / \alpha_0 < 1.02$. Setting $dn_S /dlnk=0$, yields $0.94< \alpha_{\text{dec}} / \alpha_0 < 1.01$ as previously reported. We find that a lower value of $\alpha / \alpha_0$ makes a value of $d n_S /dlnk = 0$ more compatible with the data. We also perform a thorough Fisher Matrix Analysis (including both temperature and polarization, as well as $\alpha$ and the optical depth $\tau$), in order to estimate how future CMB experiments will be able to constrain $\alpha$ and other cosmological parameters. We find that Planck data alone can constrain $\tau$ with a accuracy of the order 4% and that this constraint can be as small as 1.7% for an ideal cosmic variance limited experiment. Constraints on $\alpha$ are of the order 0.3% for Planck and can in principle be as small as 0.1% using CMB data alone - tighter constraints will require further (non-CMB) priors.

We present a summary of recent constraints on the value of the fine-structure constant at the epoch of decoupling from the recent observations made by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. Within the set of models considered, a variation of the value of $\alpha$ at decoupling with respect to the present-day value is now bounded to be smaller than 2% (6%) at 95% confidence level. We point out that the existence of an early reionization epoch as suggested by the above measurements will, when more accurate cosmic microwave background polarization data is available, lead to considerably tighter constraints. We find that the tightest possible constraint on $\alpha$ is about $0.1 %$ using CMB data alone--tighter constraints will require further (non-CMB) priors.

We discuss the value of the cosmological constant as recovered from CMB and LSS data and the robustness of the results when general isocurvature initial conditions are allowed for, as opposed to purely adiabatic perturbations. The Bayesian and frequentist statistical approaches are compared. It is shown that pre-WMAP CMB and LSS data tend to be incompatible with a non-zero cosmological constant, regardless of the type of initial conditions and of the statistical approach. The non-adiabatic contribution is constrained to be < 40% (2sigma c.l.).

We investigate in detail the question whether a non-vanishing cosmological constant is required by present-day cosmic microwave background and large scale structure data when general isocurvature initial conditions are allowed for. We also discuss differences between the usual Bayesian and the frequentist approaches in data analysis. We show that the COBE-normalized matter power spectrum is dominated by the adiabatic mode and therefore breaks the degeneracy between initial conditions which is present in the cosmic microwave background anisotropies. We find that in a flat universe the Bayesian analysis requires \Omega_\Lambda \neq 0 to more than 3 \sigma, while in the frequentist approach \Omega_\Lambda = 0 is still within 3 \sigma for a value of h < 0.48. Both conclusions hold regardless of initial conditions.

We extend our recent work on the effects of a time-varying fine-structure constant $\alpha$ in the cosmic microwave background, by providing a thorough analysis of the degeneracies between $\alpha$ and the other cosmological parameters, and discussing ways to break these with both existing and/or forthcoming data. In particular, we present the state-of-the-art CMB constraints on $\alpha$, through a combined analysis of the BOOMERanG, MAXIMA and DASI datasets. We also present a novel discussion of the constraints on $\alpha$ coming from large-scale structure observations, focusing in particular on the power spectrum from the 2dF survey. Our results are consistent with no variation in $\alpha$ from the epoch of recombination to the present day, and restrict any such (relative) variation to be less than about 4%. We show that the forthcoming MAP and (particularly) Planck experiments will be able to break most of the currently existing degeneracies between $\alpha$ and other parameters, and measure $\alpha$ to better than percent accuracy.

Recent estimates of cosmological parameters derived from Cosmic Microwave Background (CMB) anisotropies are based on the assumption that we know the precise amount of energy density in relativistic particles in the universe, $\omega_{rel}$, at all times. There are, however, many possible mechanisms that can undermine this assumption. In this paper we investigate the effect that removing this assumption has on the determination of the various cosmological parameters. We obtain fairly general bounds on the redshift of equality, $z_{eq}= \omega_{m}/\omega_{rel}=3100_{-400}^{+600}$. We show that $\omega_{rel}$ is nearly degenerate with the amount of energy in matter, $\omega_m$, and that its inclusion in CMB parameter estimation also affects the present constraints on other parameters such as the curvature or the scalar spectral index of primordial fluctuations. This degeneracy has the effect of limiting the precision of parameter estimation from the MAP satellite, but it can be broken by measurements on smaller scales such as those provided by the Planck satellite mission.

Recently high quality data of the cosmic microwave background anisotropies have been published. In this work we study to which extent the cosmological parameters determined by using this data depend on assumptions about the initial conditions. We show that for generic initial conditions, not only the best fit values are very different but, and this is our main result, the allowed parameter range enlarges dramatically.