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

Towards the application of the Maximum Entropy Method to finite temperature Upsilon Spectroscopy

According to the Narnhofer Thirring Theorem interacting systems at finite temperature cannot be described by particles with a sharp dispersion law. It is therefore mandatory to develop new methods to extract particle masses at finite temperature. The Maximum Entropy method offers a path to obtain the spectral function of a particle correlation function directly. We have implemented the method and tested it with zero temperature Upsilon correlation functions obtained from an NRQCD simulation. Results for different smearing functions are discussed.Comment: Lattice 2000 (Finite Temperature

Using Relative Entropy to Find Optimal Approximations: an Application to Simple Fluids

We develop a maximum relative entropy formalism to generate optimal approximations to probability distributions. The central results consist in (a) justifying the use of relative entropy as the uniquely natural criterion to select a preferred approximation from within a family of trial parameterized distributions, and (b) to obtain the optimal approximation by marginalizing over parameters using the method of maximum entropy and information geometry. As an illustration we apply our method to simple fluids. The "exact" canonical distribution is approximated by that of a fluid of hard spheres. The proposed method first determines the preferred value of the hard-sphere diameter, and then obtains an optimal hard-sphere approximation by a suitably weighed average over different hard-sphere diameters. This leads to a considerable improvement in accounting for the soft-core nature of the interatomic potential. As a numerical demonstration, the radial distribution function and the equation of state for a Lennard-Jones fluid (argon) are compared with results from molecular dynamics simulations.Comment: 5 figures, accepted for publication in Physica A, 200

Analytical view of diffusive and convective cosmic ray transport in elliptical galaxies

Context: An analytical solution of the generalized diffusive and convective transport equation is derived to explain the transport of cosmic ray protons within elliptical galaxies. Aims: Cosmic ray transport within elliptical galaxies is an interesting element in understanding the origin of high energetic particles measured on Earth. As probable sources of those high energetic particles, elliptical galaxies show a dense interstellar medium as a consequence of activity in the galactic nucleus or merging events between galaxies. Thus it is necessary for an appropriate description of cosmic ray transport to take the diffusive and convective processes in a dense interstellar environment into account. Here we show that the transport equations can be solved analytically with respect to the given geometry and boundary conditions in position space, as well as in momentum space. Results: The spatial solution is shown using a generalized source of cosmic rays. Additionally, the special case of a jet-like source is illustrated. We present the solution in momentum space with respect to an escape term for cosmic ray protons depending on the spatial shape of the galaxy. For a delta-shape injection function, the momentum solution is obtained analytically. We find that the spectral index measured on Earth can be obtained by appropriately choosing of the strength of Fermi I and Fermi II processes. From these results we calculate the gamma-ray flux from pion decay due to proton-proton interaction to give connection to observations. Additionally we determine the escape-spectrum of cosmic rays. The results show that both spectra are harder than the intrinsic power-law spectrum for cosmic rays in elliptical galaxies.Comment: 23 pages, 7 figures, accepted for publication in A&

Cosmic Rays X. The cosmic ray knee and beyond: Diffusive acceleration at oblique shocks

Our purpose is to evaluate the rate of the maximum energy and the acceleration rate that cosmic rays acquire in the non-relativistic diffusive shock acceleration as it could apply during their lifetime in various astrophysical sites, where highly oblique shocks exist. We examine numerically (using Monte Carlo simulations) the effect of the diffusion coefficients on the energy gain and the acceleration rate, by testing the role between the obliquity of the magnetic field to the shock normal, and the significance of both perpendicular cross-field diffusion and parallel diffusion coefficients to the acceleration rate. We find (and justify previous analytical work - Jokipii 1987) that in highly oblique shocks the smaller the perpendicular diffusion gets compared to the parallel diffusion coefficient values, the greater the energy gain of the cosmic rays to be obtained. An explanation of the cosmic ray spectrum in high energies, between $10^{15}$eV and about $10^{18}$eV is claimed, as we estimate the upper limit of energy that cosmic rays could gain in plausible astrophysical regimes; interpreted by the scenario of cosmic rays which are injected by three different kind of sources, (a) supernovae which explode into the interstellar medium, (b) Red Supergiants, and (c) Wolf-Rayet stars, where the two latter explode into their pre-supernovae winds.Comment: Accepted in Astronomy and Astrophysics, 9 pages, 8 figures (for the 'Cosmic Rays' series papers

A Complete Spectroscopic Survey of the Milky Way satellite Segue 1: Dark matter content, stellar membership and binary properties from a Bayesian analysis

We introduce a comprehensive analysis of multi-epoch stellar line-of-sight velocities to determine the intrinsic velocity dispersion of the ultrafaint satellites of the Milky Way. Our method includes a simultaneous Bayesian analysis of both membership probabilities and the contribution of binary orbital motion to the observed velocity dispersion within a 14-parameter likelihood. We apply our method to the Segue 1 dwarf galaxy and conclude that Segue 1 is a dark-matter-dominated galaxy at high probability with an intrinsic velocity dispersion of 3.7^{+1.4}_{-1.1} km/sec. The dark matter halo required to produce this dispersion must have an average density of 2.5^{+4.1}_{-1.9} solar mass/pc^3 within a sphere that encloses half the galaxy's stellar luminosity. This is the highest measured density of dark matter in the Local Group. Our results show that a significant fraction of the stars in Segue 1 may be binaries with the most probable mean period close to 10 years, but also consistent with the 180 year mean period seen in the solar vicinity at about 1 sigma. Despite this binary population, the possibility that Segue 1 is a bound star cluster with the observed velocity dispersion arising from the orbital motion of binary stars is disfavored by the multi-epoch stellar velocity data at greater than 99% C.L. Finally, our treatment yields a projected (two-dimensional) half-light radius for the stellar profile of Segue 1 of 28^{+5}_{-4} pc, in excellent agreement with photometric measurements.Comment: 15 pages, 19 figure

Approximate inference of the bandwidth in multivariate kernel density estimation

Kernel density estimation is a popular and widely used non-parametric method for data-driven density estimation. Its appeal lies in its simplicity and ease of implementation, as well as its strong asymptotic results regarding its convergence to the true data distribution. However, a major difficulty is the setting of the bandwidth, particularly in high dimensions and with limited amount of data. An approximate Bayesian method is proposed, based on the Expectation–Propagation algorithm with a likelihood obtained from a leave-one-out cross validation approach. The proposed method yields an iterative procedure to approximate the posterior distribution of the inverse bandwidth. The approximate posterior can be used to estimate the model evidence for selecting the structure of the bandwidth and approach online learning. Extensive experimental validation shows that the proposed method is competitive in terms of performance with state-of-the-art plug-in methods

Constraints on cosmic-ray propagation models from a global Bayesian analysis

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.Comment: 19 figures, 3 tables, emulateapj.sty. A typo is fixed. To be published in the Astrophysical Journal v.728 (February 10, 2011 issue). Supplementary material can be found at http://www.g-vo.org/pub/GALPROP/GalpropBayesPaper

The Interaction of Cosmic Rays with Diffuse Clouds

We study the change in cosmic-ray pressure, the change in cosmic-ray density, and the level of cosmic-ray induced heating via Alfven-wave damping when cosmic rays move from a hot ionized plasma to a cool cloud embedded in that plasma. The general analysis method outlined here can apply to diffuse clouds in either the ionized interstellar medium or in galactic winds. We introduce a general-purpose model of cosmic-ray diffusion building upon the hydrodynamic approximation for cosmic rays (from McKenzie & Voelk and Breitschwerdt and collaborators). Our improved method self-consistently derives the cosmic-ray flux and diffusivity under the assumption that the streaming instability is the dominant mechanism for setting the cosmic-ray flux and diffusion. We find that, as expected, cosmic rays do not couple to gas within cool clouds (cosmic rays exert no forces inside of cool clouds), that the cosmic-ray density does not increase within clouds (it may slightly decrease in general, and decrease by an order of magnitude in some cases), and that cosmic-ray heating (via Alfven-wave damping and not collisional effects as for ~10 MeV cosmic rays) is only important under the conditions of relatively strong (10 micro-Gauss) magnetic fields or high cosmic-ray pressure (~10^{-11} ergs cm^{-3}).Comment: 15 pages, 9 figures, 1 table; Accepted to Ap

Hadronic Spectral Functions in Lattice QCD

QCD spectral functions of hadrons in the pseudo-scalar and vector channels are extracted from lattice Monte Carlo data of the imaginary time Green's functions. The maximum entropy method works well for this purpose, and the resonance and continuum structures in the spectra are obtained in addition to the ground state peaks.Comment: 4 pages, 3 eps-figures, revtex (minor modifications in the text and an added reference). To appear in Physical Review D Rapid Communication