Source author record

Roberto Trotta

Roberto Trotta appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2026arXiv

Simulation-based population inference of LISA's Galactic binaries: Bypassing the global fit

The Laser Interferometer Space Antenna (LISA) is expected to detect thousands of individually resolved gravitational wave sources, overlapping in time and frequency, on top of unresolved astrophysical and/or primordial backgrounds. Disentangling resolved sources from backgrounds and extracting their parameters in a computationally intensive "global fit" is normally regarded as a necessary step toward reconstructing the properties of the underlying astrophysical populations. Here, we show that it is in principle feasible to infer the population properties of the most numerous of LISA sources -- Galactic double white dwarfs -- directly from the frequency (or, equivalently, time) strain series by adopting a simulation-based approach, without extracting and estimating the parameters of each single source. By training a normalizing flow on a custom-designed compression of simulated LISA frequency series from the Galactic double white dwarf population, we demonstrate how to infer the posterior distribution of population parameters (e.g., mass function, frequency, and spatial distributions). This allows for extracting information on the population parameters from both resolved and unresolved sources simultaneously and in a computationally efficient manner. This approach can be extended to other source classes (e.g., massive and stellar-mass black holes, extreme mass ratio inspirals) and to scenarios involving non-Gaussian or non-stationary noise (e.g., data gaps), provided that fast and accurate simulations are available.

preprint2022arXiv

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.

preprint2020arXiv

A robust estimate of the Milky Way mass from rotation curve data

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.

preprint2020arXiv

STACCATO: A Novel Solution to Supernova Photometric Classification with Biased Training Sets

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.

preprint2020arXiv

Structure Formation Models Weaken Limits on WIMP Dark Matter from Dwarf Spheroidal Galaxies

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}$.

preprint2016arXiv

Effective Field Theory of Dark Matter: a Global Analysis

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.

preprint2016arXiv

Global analysis of the pMSSM in light of the Fermi GeV excess: prospects for the LHC Run-II and astroparticle experiments

We present a new global fit of the 19-dimensional phenomenological Minimal Supersymmetric Standard Model (pMSSM-19) that comply with all the latest experimental results from dark matter indirect, direct and accelerator dark matter searches. We show that the model provides a satisfactory explanation of the excess of gamma-rays from the Galactic centre observed by the Fermi~Large Area Telescope, assuming that it is produced by the annihilation of neutralinos in the Milky Way halo. We identify two regions that pass all the constraints: the first corresponds to neutralinos with a mass ~80-100 GeV annihilating into WW with a branching ratio of 95% ; the second to heavier neutralinos, with mass ~180-200 GeV annihilating into t tbar with a branching ratio of 87%. We show that neutralinos compatible with the Galactic centre GeV excess will soon be within the reach of LHC run-II -- notably through searches for charginos and neutralinos, squarks and light smuons -- and of Xenon1T, thanks to its unprecedented sensitivity to spin-dependent cross-section off neutrons.

preprint2016arXiv

Standardizing Type Ia supernovae using Near Infrared rebrightening time

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.

preprint2015arXiv

Cosmology and fundamental physics with the Euclid satellite

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.

preprint2014arXiv

Compatibility of Planck and BICEP2 in the Light of Inflation

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.

preprint2014arXiv

The Best Inflationary Models After Planck

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.

preprint2013arXiv

Taming astrophysical bias in direct dark matter searches

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σ. 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.

preprint2012arXiv

Dark Matter Searches: The Nightmare Scenario

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.

preprint2012arXiv

Evidence for dark matter modulation in CoGeNT?

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.

preprint2012arXiv

Fundamental statistical limitations of future dark matter direct detection experiments

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.

preprint2011arXiv

Complementarity of Dark Matter Direct Detection Targets

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.

preprint2011arXiv

Hunting Down the Best Model of Inflation with Bayesian Evidence

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.

preprint2011arXiv

Quantifying the tension between the Higgs mass and (g-2)_mu in the CMSSM

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.

preprint2010arXiv

Global fits of the Minimal Universal Extra Dimensions scenario

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.

preprint2010arXiv

Identification of Dark Matter particles with LHC and direct detection data

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.

preprint2009arXiv

How flat can you get? A model comparison perspective on the curvature of the Universe

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]

preprint2009arXiv

The Virtues of Frugality - Why cosmological observers should release their data slowly

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.

preprint2007arXiv

On the detectability of the CMSSM light Higgs boson at the Tevatron

We examine the prospects of detecting the light Higgs h^0 of the Constrained MSSM at the Tevatron. To this end we explore the CMSSM parameter space with μ>0, using a Markov Chain Monte Carlo technique, and apply all relevant collider and cosmological constraints including their uncertainties, as well as those of the Standard Model parameters. Taking 50 GeV < m_{1/2}, m_0 < 4 TeV, |A_0| < 7 TeV and 2 < tan(beta) < 62 as flat priors and using the formalism of Bayesian statistics we find that the 68% posterior probability region for the h^0 mass lies between 115.4 GeV and 120.4 GeV. Otherwise, h^0 is very similar to the Standard Model Higgs boson. Nevertheless, we point out some enhancements in its couplings to bottom and tau pairs, ranging from a few per cent in most of the CMSSM parameter space, up to several per cent in the favored region of tan(beta)\sim 50 and the pseudoscalar Higgs mass of m_A\lsim 1 TeV. We also find that the other Higgs bosons are typically heavier, although not necessarily much heavier. For values of the h^0 mass within the 95% probability range as determined by our analysis, a 95% CL exclusion limit can be set with about 2/fb of integrated luminosity per experiment, or else with 4/fb (12/fb) a 3 sigma evidence (5 sigma discovery) will be guaranteed. We also emphasize that the alternative statistical measure of the mean quality-of-fit favors a somewhat lower Higgs mass range; this implies even more optimistic prospects for the CMSSM light Higgs search than the more conservative Bayesian approach. In conclusion, for the above CMSSM parameter ranges, especially m_0, either some evidence will be found at the Tevatron for the light Higgs boson or, at a high confidence level, the CMSSM will be ruled out.

preprint2007arXiv

The isocurvature fraction after WMAP 3-year data

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.

preprint2006arXiv

A Markov Chain Monte Carlo Analysis of the CMSSM

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 μ^+ μ^-) and BR({\bar B}\to X_sγ) 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.

preprint2003arXiv

The cosmological constant and the paradigm of adiabaticity

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.).