Fast and reliable MCMC for cosmological parameter estimation

Abstract

Markov Chain Monte Carlo (MCMC) techniques are now widely used for cosmological parameter estimation. Chains are generated to sample the posterior probability distribution obtained following the Bayesian approach. An important issue is how to optimize the efficiency of such sampling and how to diagnose whether a finite-length chain has adequately sampled the underlying posterior probability distribution. We show how the power spectrum of a single such finite chain may be used as a convergence diagnostic by means of a fitting function, and discuss strategies for optimizing the distribution for the proposed steps. The methods developed are applied to current CMB and LSS data interpreted using both a pure adiabatic cosmological model and a mixed adiabatic/isocurvature cosmological model including possible correlations between modes. For the latter application, because of the increased dimensionality and the presence of degeneracies, the need for tuning MCMC methods for maximum efficiency becomes particularly acute.

12 pages, 17 figures. Submitted to MNRAS

References in corpus (13)

• Cosmological parameters from SDSS and WMAP
• The 3D power spectrum of galaxies from the SDSS
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Implications for Inflation
• Wilkinson Microwave Anisotropy Probe (WMAP) First Year Observations: TE Polarization
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Angular Power Spectrum
• The 2dF Galaxy Redshift Survey: The power spectrum and the matter content of the universe
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Parameter Estimation Methodology
• Efficient Cosmological Parameter Estimation from Microwave Background Anisotropies
• II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements
• The Age of the Universe and the Cosmological Constant Determined from Cosmic Microwave Background Anisotropy Measurements
• Analyze This! A Cosmological Constraint Package for CMBEASY
• The initial conditions of the universe: how much isocurvature is allowed?
• Fast Cosmological Parameter Estimation from Microwave Background Temperature and Polarization Power Spectra

Cited by in corpus (186)

• emcee: The MCMC Hammer
• Wilkinson Microwave Anisotropy Probe (WMAP) Three Year Results: Implications for Cosmology
• The Average Star Formation Histories of Galaxies in Dark Matter Halos from z=0-8
• Axion Cosmology
• The propagation of uncertainties in stellar population synthesis modeling I: The relevance of uncertain aspects of stellar evolution and the IMF to the derived physical properties of galaxies
• New Hubble Space Telescope Discoveries of Type Ia Supernovae at z > 1: Narrowing Constraints on the Early Behavior of Dark Energy
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP data
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• A Comprehensive Analysis of Uncertainties Affecting the Stellar Mass - Halo Mass Relation for 0<z<4
• The Atlas3D project - XV. Benchmark for early-type galaxies scaling relations from 260 dynamical models: mass-to-light ratio, dark matter, Fundamental Plane and Mass Plane
• Efficient sampling of fast and slow cosmological parameters
• The Stellar Initial Mass Function in Early-Type Galaxies From Absorption Line Spectroscopy. II. Results
• Cosmological Constraints from the SDSS maxBCG Cluster Catalog
• The Atacama Cosmology Telescope: Cosmological Parameters from the 2008 Power Spectra
• Two accurate time-delay distances from strong lensing: Implications for cosmology
• The Flat Transmission Spectrum of the Super-Earth GJ1214b from Wide Field Camera 3 on the Hubble Space Telescope
• Evolution of the Galaxy - Dark Matter Connection and the Assembly of Galaxies in Dark Matter Halos
• Correlation of CMB with large-scale structure: I. ISW Tomography and Cosmological Implications
• A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old
• Cosmology with the lights off: Standard sirens in the Einstein Telescope era
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Galactic Foreground Emission
• The Halos of Satellite Galaxies: the Companion of the Massive Elliptical Lens SL2S J08544-0121
• CLASH-VLT: Insights on the mass substructures in the Frontier Fields Cluster MACS J0416.1-2403 through accurate strong lens modeling
• The Atacama Cosmology Telescope: likelihood for small-scale CMB data
• Using the Full Power of the Cosmic Microwave Background to Probe Axion Dark Matter
• Modeling the evolution of infrared galaxies: A Parametric backwards evolution model
• Cluster Mass Profiles from a Bayesian Analysis of Weak Lensing Distortion and Magnification Measurements: Applications to Subaru Data
• Modelling clusters of galaxies by f(R)-gravity
• A New Photometric Model of the Galactic Bar using Red Clump Giants
• Markov Chain Monte Carlo Exploration of Minimal Supergravity with Implications for Dark Matter
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Bayesian Estimation of CMB Polarization Maps
• Cosmological constraints on Brans-Dicke theory
• A precise and accurate determination of the cosmic microwave background temperature at z=0.89
• The Atacama Cosmology Telescope: a measurement of the primordial power spectrum
• HerMES: Deep Galaxy Number Counts from a P(D) Fluctuation Analysis of SPIRE Science Demonstration Phase Observations
• Fast Hamiltonian sampling for large scale structure inference
• Stellar Populations in the Central 0.5 pc of the Galaxy I: A New Method for Constructing Luminosity Functions and Surface-Density Profiles
• The Impact of Modeling Assumptions in Galactic Chemical Evolution Models
• Constraining Isocurvature Initial Conditions with WMAP 3-year data
• Unbiased constraints on ultralight axion mass from dwarf spheroidal galaxies
• Comprehensive cosmographic analysis by Markov Chain Method
• Redshift-Space Distortions with the Halo Occupation Distribution II: Analytic Model
• Bayesian Analysis of Inflation: Parameter Estimation for Single Field Models
• Bayesian non-linear large scale structure inference of the Sloan Digital Sky Survey data release 7
• Model-independent constraints on reionization from large-scale CMB polarization
• Multiple inflation and the WMAP 'glitches' II. Data analysis and cosmological parameter extraction
• Lensing Constraints on the Mass Profile Shape and the Splashback Radius of Galaxy Clusters
• SED fitting with Markov Chain Monte Carlo: Methodology and Application to z=3.1 Lyman Alpha Emitting Galaxies
• Estimation of cosmological parameters using adaptive importance sampling
• The Megamaser Cosmology Project. VI. Observations of NGC 6323
• Direct reconstruction of the quintessence potential
• Galactic Archaeology with asteroseismic ages. II. Confirmation of a delayed gas infall using Bayesian analysis based on MCMC methods
• Comparison of sampling techniques for Bayesian parameter estimation
• Cosmological constraints from baryon acoustic oscillations and clustering of large-scale structure
• Falsifying Paradigms for Cosmic Acceleration
• CMB constraints on primordial black hole dark matter
• Breaking the spell of Gaussianity: forecasting with higher order Fisher matrices
• An improved test of the strong equivalence principle with the pulsar in a triple star system
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Geometrical observational bounds on a fractal horizon holographic dark energy
• Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module
• On the Luminosity Dependence of the Galaxy Pairwise Velocity Dispersion
• Rigorous luminosity function determination in presence of a background: theory and application to two intermediate redshift clusters
• Constraining Lorentz violation with cosmology
• A Markov Chain Monte Carlo for Galactic Cosmic Ray physics: I. Method and results for the Leaky Box Model
• Crossing SNe Ia and BAO observational constraints with local ones in hybrid metric-Palatini gravity
• The Extreme Small Scales: Do Satellite Galaxies Trace Dark Matter?
• Fuzzy Dark Matter and the Dark Energy Survey Year 1 Data
• SHARP - III: First Use Of Adaptive Optics Imaging To Constrain Cosmology With Gravitational Lens Time Delays
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• Cosmography from two-image lens systems: overcoming the lens profile slope degeneracy
• Spin Dependence of Dark Matter Scattering
• The Galaxy-Halo/Subhalo connection: Mass relations and implications for some satellite occupational distributions
• The XMM Cluster Survey: Forecasting cosmological and cluster scaling-relation parameter constraints
• Parameters Estimation for the Cosmic Microwave Background with Bayesian Neural Networks
• The Spectrum of the Axion Dark Sector
• The mass and anisotropy profiles of galaxy clusters from the projected phase space density: testing the method on simulated data
• Oscillations in the dark energy EoS: new MCMC lessons
• Constraints on primordial isocurvature perturbations and spatial curvature by Bayesian model selection
• Measuring the X-ray luminosities of SDSS DR7 clusters from RASS
• Interpreting the strongly lensed supernova iPTF16geu: time delay predictions, microlensing, and lensing rates
• The Strikingly Similar Relation between Satellite and Central Galaxies and Their Dark Matter Halos Since z=2
• MACS J0416.1-2403: Impact of line-of-sight structures on strong gravitational lensing modelling of galaxy clusters
• MAXIPOL: Data Analysis and Results
• Testable dark energy predictions from current data
• Robust sampling for weak lensing and clustering analyses with the Dark Energy Survey
• The Cosmic Far-Infrared Background Buildup Since Redshift 2 at 70 and 160 microns in the COSMOS and GOODS fields
• The Hubble Space Telescope Cluster Supernova Survey: III. Correlated Properties of Type Ia Supernovae and Their Hosts at 0.9 < z < 1.46
• Sampling using a `bank' of clues
• Observational constraints on cosmological future singularities
• Power-Law Template for IR Point Source Clustering
• Modeling the Very Small-Scale Clustering of Luminous Red Galaxies
• Cold Dark Matter Isocurvature Perturbations: Constraints and Model Selection
• Cosmological constraints on the Generalized Uncertainty Principle from modified Friedmann equations
• Anisotropic expansion and SNIa: an open issue
• Figures of merit for present and future dark energy probes
• Constraining thermal dust emission in distant galaxies with number counts and angular power spectra
• Dark Matter and Dark Energy from a single scalar field and CMB data
• TDCOSMO X. Automated Modeling of 9 Strongly Lensed Quasars and Comparison Between Lens Modeling Software
• Galaxy halo truncation and Giant Arc Surface Brightness Reconstruction in the Cluster MACSJ1206.2-0847
• Bayesian Methods in Cosmology
• Bayesian Limits on Primordial Isotropy Breaking
• Annealing a Follow-up Program: Improvement of the Dark Energy Figure of Merit for Optical Galaxy Cluster Surveys
• A spatial correlation analysis for a toroidal universe
• Constraints on scalar and tensor perturbations in phenomenological and two-field inflation models: Bayesian evidences for primordial isocurvature and tensor modes
• A comparison of optimisation algorithms for high-dimensional particle and astrophysics applications
• candl: Cosmic Microwave Background Analysis with a Differentiable Likelihood
• Measuring the geometry of the Universe in the presence of isocurvature modes
• Breaking the Vainshtein screening in clusters of galaxies
• Impact of reionization on CMB polarization tests of slow-roll inflation
• The Galaxy Cross-Correlation Function as a Probe of the Spatial Distribution of Galactic Satellites
• Testing flatness of the universe with probes of cosmic distances and growth
• Determining neutrino properties using future galaxy redshift surveys
• The global mass - to - light ratio of SLACS lenses
• Physics of the cosmic microwave background anisotropy
• Physical component analysis of galaxy cluster weak gravitational lensing data
• Cosmography and large scale structure by f(R) gravity: new results
• Two-component scalar dark matter in $Z_{2n}$ scenarios
• J0454-0309: Evidence for a strong lensing fossil group falling into a poor galaxy cluster
• Likelihood reconstruction method of real-space density and velocity power spectra from a redshift galaxy survey
• Ultra high energy neutrino-nucleon cross section from cosmic ray experiments and neutrino telescopes
• Testing non-local gravity by clusters of galaxies
• The 2L1S/1L2S Degeneracy for Two Microlensing Planet Candidates Discovered by the KMTNet Survey in 2017
• New method for Earth neutral atmospheric density retrieval based on energy spectrum fitting during occultation with LE/\emph{Insight}-HXMT
• HOLISMOKES -- X. Comparison between neural network and semi-automated traditional modeling of strong lenses
• A Study of the Point Spread Function in SDSS Images
• On the Redshift Evolution of MgII Absorption Systems
• Validating the Fisher approach for stage IV spectroscopic surveys
• OGLE-2016-BLG-1045: A Test of Cheap Space-Based Microlens Parallaxes
• Improving Dark Energy Constraints with High Redshift Type Ia Supernovae from CANDELS and CLASH
• A generalized mass-to-horizon relation: a new global approach to entropic cosmologies and its connection to $\Lambda$CDM
• HST/STIS capability for Love number measurement of WASP-121b
• Model independent approaches to reionization in the analysis of upcoming CMB data
• Multi-component DHOST analysis in galaxy clusters
• Small scale contributions to CMB: A coherent analysis
• Forecast and analysis of the cosmological redshift drift
• Two Jupiter-Mass Planets Discovered by the KMTNet Survey in 2017
• Euclid preparation. LIII. LensMC, weak lensing cosmic shear measurement with forward modelling and Markov Chain Monte Carlo sampling
• Cosmological constraints on a unified dark matter-energy scalar field model with fast transition
• Cosmological Constraints on Entropic Cosmology with Matter Creation
• First cosmological constraints on the Superfluid Chaplygin gas model
• AGEL: Is the Conflict Real? Investigating Galaxy Evolution Models using Strong Lensing at 0.3 < z < 0.9
• Observational constraints from SNe Ia and Gamma-Ray Bursts on a clumpy universe
• Observational constraints on electromagnetic Born-Infeld cosmology
• Interacting quintessence cosmology from Noether symmetries: comparing theoretical predictions with observational data
• On the concordance of cosmological data in the case of the generalized Chaplygin gas
• Energy and temperature dependencies for electron-induced sputtering from H${}_2$O-ice: Implications for the icy Galilean moons
• The XMM Cluster Survey: The Halo Occupation Number of BOSS galaxies in X-ray clusters
• Parameter estimation of a nonlinear magnetic universe from observations
• Cosmological Parameters from CMB Maps without Likelihood Approximation
• Testing feasibility of scalar-tensor gravity by scale dependent mass and coupling to matter
• DHOST gravity in Ultra-diffuse galaxies -- Part II: NGC 1052-DF4 and Dragonfly 44
• Large-Scale Polarized Foreground Component Separation for Planck
• Unifying static analysis of gravitational structures with a scale-dependent scalar field gravity as an alternative to dark matter
• Tests of Standard Cosmology in Horava Gravity and Bayesian Evidence for a Closed Universe, and the Hubble Tension
• Measurement of the vertical atmospheric density profile from the X-ray Earth occultation of the Crab Nebula with Insight-HXMT
• BEAMS: separating the wheat from the chaff in supernova analysis
• Deep Learning and genetic algorithms for cosmological Bayesian inference speed-up
• The flat-sky approximation to galaxy number counts -- redshift space correlation function
• Old Dark Energy
• Cosmological constrains on fast transition Unified Dark Energy/Matter models
• Deep Involutive Generative Models for Neural MCMC
• First measurement of the total gravitational quadrupole moment of a black widow companion
• Observational constraints on the unified dark matter and dark energy model based on the quark bag model
• Mapping the Cosmological Confidence Ball Surface
• Numerical search for a potential planet sculpting the young disc of HD 115600
• A Bayesian chemical evolution model of the DustPedia Galaxy M74
• Constraining coronal heating: employing Bayesian analysis techniques to improve the determination of solar atmospheric plasma parameters
• Constraining Ultralight Axions with Galaxy Surveys
• Quantifying the Impact of Starspot-Crossing Events on Retrieved Parameters from Transit Lightcurves
• HOLISMOKES -- XIV. Time-delay and differential dust extinction determination with lensed type II supernova color curves
• Accelerating MCMC algorithms through Bayesian Deep Networks
• Using Source Proper Motion to Validate Terrestrial Parallax: OGLE-2019-BLG-1058
• Sunyaev-Zel'dovich contribution in CMB analysis
• A possible geometrical origin of the accelerated expansion of the universe
• Dark energy constraints and correlations with systematics from CFHTLS weak lensing, SNLS supernovae Ia and WMAP5
• A Panchromatic JWST Spectrum of a Giant Starspot on the Fully Convective M-dwarf TOI-3884
• Explanation of the exceptionally strong timing noise of PSR J0337+1715 by a circum-ternary planet and consequences for gravity tests
• DHOST gravity in Ultra-diffuse galaxies -- Part I: the case of NGC1052-DF2
• Bayesian Inversion of Stokes Profiles
• Determination of Dark Matter Properties at High-Energy Colliders
• BICEP/Keck XXI: Constraints on Early-Universe Parity Violation from Multipole-Dependent Birefringence
• Strong Lensing Model and Dust Extinction Maps of the Host Galaxy of Type Ia Supernova H0pe
• Optimisation of the Population Monte Carlo algorithm: Application to constraining isocurvature models with cosmic microwave background data
• On the energy flow of $\lambda$ in Hořava-Lifshitz cosmology
• APES: Approximate Posterior Ensemble Sampler