Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant 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, which uses astrophysical information, and nuclear and particle data as inputs 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.

Abstract: 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 similar to 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 fb(-1) of data should be able to probe the current best-fit UED parameters.

Abstract: We give a straightforward procedure to scan the seesaw parameter-space, using the common “R-parametrization”, in a complete way. This includes a very simple rule to incorporate the perturbativity requirement as a condition for the entries of the R-matrix. As a relevant application, we show that the somewhat propagated belief that BR(mu -> e, gamma) in supersymmetric seesaw models depends strongly on the value of theta(13) is an “optical effect” produced by incomplete scans, and does not hold after a careful analytical and numerical study. When the complete scan is done, BR(mu -> e, gamma) gets very insensitive to theta(13). This holds even if the right-handed neutrino masses are kept constant or under control (as is required for succesful leptogenesis). In most cases the values of BR(mu -> e, gamma) are larger than the experimental upper bound. Including (unflavoured) leptogenesis does not introduce any further dependence on theta(13), although decreases the typical value of BR(mu -> e, gamma).

Abstract: Supersymmetry has often been invoked 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 constrained minimal supersymmetric standard model (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 standard model 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 the current data (similar to 3.2 sigma) is already achieved. Requiring less than 3 sigma discrepancy, implies m(h) less than or similar to 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 an inconsistency.

Abstract: 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.