2018 Cosmological parameters and MC chains

From Planck Legacy Archive Wiki
Jump to: navigation, search

Description[edit]

The 2018 cosmological parameter results explore a variety of cosmological models with combinations of Planck and other data. We provide results from MCMC exploration chains, as well as best fits, and sets of parameter tables. Definitions, conventions, and references are contained in Planck-2013-XVI[1], Planck-2015-A15[2], and the 2018 parameter paper Planck-2018-L06[3].

Production process[edit]

Parameter chains are produced using CosmoMC, a sampling package available here. This includes the sample analysis package (and GUI) GetDist, and the scripts for managing, analysing, and plotting results from the full grid of runs. The Python GetDist sample analysis package is also available separately.

Chain products provided here have had burn-in removed. Some results with additional data are produced by importance sampling. Note that the baseline model includes one massive (0.06eV) neutrino.

Caveats and known issues[edit]

  1. Confidence intervals are derived from the MCMC samples, and assume that the input likelihoods are exactly correct, so there is no quantification for systematic errors other than via the covariance, foreground, and beam-error models assumed in the likelihood codes.
  2. Nonlinear lensing modelling uses Halofit (HMCode); for some extended models and CMB lensing-only analyses, tails of the chains may lie away from the domain of validity.
  3. Polarization results are sensitive to details of the polarization modelling; systematic uncertainties is not accounted for and may change results by up to one standard deviation in some cases (but is difficult to fully quantify). Alternative CamSpecHM likelihood results are also provided and give some idea of differences that can be obtained with different analysis choices.
  4. Quoted values for the helium abundance (YP and YPBBN) and D/H, when YP is not varied as a parameter, are obtained assuming the Parthenope 1.1 BBN code, and do not include theoretical errors.
  5. Some nuisance parameters have slightly different definitions in the plikHM and CamSpecHM likelihoods, and cannot be directly compared.

Related products[edit]

Results of the baseline parameter exploration runs should be reproducible using CosmoMC with the Planck 2018 likelihood codes. A new Python code called Cobaya also follows the same methodology and should produce identical results.

Parameter tables[edit]

These tables list parameter constraints for each considered model and data combination separately. For the baseline likelihood combinations see:

There are also larger full files including results from the alternative CamSpec likelihood, lensing only, other data combinations and priors:

A full set of parameter tables, including various ways of comparing data combinations, are available in one bundle from the PLA as shown below.

Data combination tags used to label results are as follows (see Planck-2020-A6[3] for full description and references).


Tag Data
plikHM Baseline high-ℓ Planck power spectra (plik cross-half-mission, 30 ≤ ℓ ≤ 2508)
CamSpecHM High-ℓ Planck power spectra (CamSpec cross-half-mission, 30 ≤ ℓ ≤ 2500)
CleanedCamSpecHM Foreground-cleaned high-L Planck CamSpec spectra (cross-half-mission, 30 ≤ ℓ ≤ 2500)
lowl Low-ℓ: Planck temperature only (2 ≤ ℓ ≤ 29)
lowE Low-ℓ HFI polarization EE likelihood (2 ≤ ℓ ≤ 29)
lensing Planck lensing conservative power spectrum reconstruction likelihood
zre6p5 A hard prior zre > 6.5
reion A hard prior zre > 6.5, combined with Gaussian prior zre = 7 ± 1
BAO Baryon oscillation data from DR12, MGS, and 6DF
Pantheon18 Supernova data from the Pantheon sample, with updated main distance file with heliocentric redshifts
JLA Supernova data from the SDSS-II/SNLS3 Joint Light-curve Analysis
Riess18 Local Hubble parameter measurement from Riess et al.(a), H0 = 73.45 ± 1.66
BK15 Bicep-Keck (+Planck/WMAP) 2015 analysis (arXiv:1810.05216)
theta θMC = 1.0409 ± 0.0006 Gaussian prior
WMAP The full WMAP (temperature and polarization) 9-year data
lenspriors Standard base parameters with ns = 0.96 ± 0.02, Ωbh2 = 0.0222 ± 0.0005, 100>H0>40, τ=0.055
DESpriors DES cosmological parameter priors (flat on 0.1< Ωm<0.9, 0.03<Ωb<0.07, 55<H0<91, 0.5<109As<5, YP=0.245341 and, if varied, 0.05 eV< Σ mν <1 eV)
DES DES 1yr, cosmic shear+galaxy auto+cross
DESlens DES 1yr, cosmic shear only
DESwt DES 1yr, galaxy auto+cross only

The high-ℓ Planck likelihoods have TT, TE, EE variants from each spectrum alone, plus the TTTEEE joint constraint. When the lensing likelihood is used with DESpriors or lenspriors, it is marginalized over the theoretical CMB power spectra (as described in the 2018 lensing paper).


Tags used to identify the model parameters that are varied are described in the introduction to the PDF table files.

Parameter chains[edit]

We provide the full chains and getdist outputs for our parameter results. The entire grid of results, including likelihood variations, various external data combinations and CMB lensing only results, is available as a 9GB compressed file:

You can also download smaller files containing key results:

The download contains a hierarchy of directories, with each separate chain in a separate directory. The structure for the directories is

base_AAA_BBB/XXX_YYY_.../ ,

where AAA and BBB are any additional parameters that are varied in addition to the six parameters of the baseline model. XXX, YYY, etc encode the data combinations used. These follow the naming conventions described above under Parameter Tables. Each directory contains the main chains, 4-8 text files with one chain in each, and various other files all with names of the form

base_AAA_BBB_XXX_YYY.ext ,

where ext describes the type of file, and the possible values or ext are as follows.


Extension Data
.txt Parameter chain file with burn-in removed
.paramnames File that describes the parameters included in the chains
.inputparams Input parameters used when generating the chain
.minimum Best-fit parameter values, −log likelihoods and chi-square
.minimum.theory_cl Best-fit temperature and polarization power spectra and lensing potential (see below)
.minimum.plik_foregrounds Best-fit foreground model (additive component) for each data power spectrum used
.minimum.inputparams Input parameters used when generating the best fit
.ranges Prior ranges assumed for each parameter


In addition each directory contains any importance sampled outputs with additional data. These have names of the form

base_AAA_BBB_XXX_YYY_post_ZZZ.ext ,

where ZZZ is the data likelihood that is added by importance sampling. Finally, each directory contains a dist subdirectory, containing results of chain analysis. File names follow the above conventions, with the following extensions.


Extension Data
.margestats Mean, variance and 68, 95 and 99% limits for each parameter (see below)
.likestats Parameters of best-fitting sample in the chain (generally different from the .minmum global best fit)
.covmat Covariance matrix for the MCMC parameters
.converge Summary of various convergence diagnostics


Python scripts for reading in chains and calculating new derived parameter constraints are available as part of CosmoMC, see the readme for details [1]. The config directory in the download includes information about the grid configuration used by the plotting and grid scripts.

File formats[edit]

The file formats are standard CosmoMC/GetDist outputs. GetDist includes python scripts for generating tables, 1D, 2D, and 3D plots using the provided data, as well as a GUI for conveniently making plots from grid downloads. The formats are summarised here.

Chain files
Each chain file is ASCII and contains one sample on each line. Each line is of the format
weight like param1 param2 param3 …
Here weight is the importance weight or multiplicity count, and like is the total −log Likelihood. param1,param2, etc are the parameter values for the sample, where the numbering is defined by the position in the accompanying .paramnames files.
Note that burn-in has been removed from the cosmomc outputs, so full chains provided can be used for analysis. Importance sampled results (with _post) in the name have been thinned by a factor of 10 compared to the original chains, so the files are smaller, but this does not significantly affect the effective number of samples. Due to the way MCMC works, the samples in the chain outputs are not independent, but it is safe to use all the samples for estimating posterior averages.
.margestats files
Each row contains the marginalized constraint on individual parameters. The format is fairly self explanatory given the text description in the file, with each line of the form
parameter mean sddev lower1 upper1 limit1 lower2 upper2 limit2 lower3 upper3 limit3
where sddev is the standard deviation, and the limits are: 1, 68%; 2, 95%; and 3, 99%. The limit tags specify whether a given limit is one tailed, two tailed or none (if no constraint within the assumed prior boundary).
.minimum.theory_cl files
They contain the best-fit theoretical power spectra (without foregrounds) for each model. The columns are ℓ, DTT, DTE, DEE, DBB, and Ddd, where D ≡ ℓ(ℓ+1)C/(2π) in μK2. Also Ddd= [ℓ(ℓ+1)]2Cφφ/(2π) is the power spectrum of the lensing deflection angle, where Cφφ is the lensing potential power spectrum. Note that the lensing spectrum may not be accurate at ℓ > 400 due to the maximum wavenumber and nonlinear correction accuracy settings.


Previous Releases: (2015) and (2013) Cosmological Parameters and MC Chains[edit]

Expand

2015 Release of Cosmological Parameters and MC Chains

Expand

2013 Release of Cosmological parameters and MC Chains


References[edit]

  1. Jump up to: 1.01.11.21.31.41.5 Planck 2013 results. XVI. Cosmological parameters, Planck Collaboration, 2014, A&A, 571, A16.
  2. Jump up to: 2.02.1 Planck 2015 results. XIII. Cosmological parameters, Planck Collaboration, 2016, A&A, 594, A13.
  3. Jump up to: 3.03.1 Planck 2018 results. VI. Cosmological parameters, Planck Collaboration, 2020, A&A, 641, A6.
  4. Jump up Using BBN in cosmological parameter extraction fromCMB: A Forecast for PLANCK, J. Hamann, J. Lesgourgues, G. Mangano, J. Cosmology Astropart. Phys., 0803, 004, (2008).

Cosmic Microwave background

Planck Legacy Archive

(Planck) High Frequency Instrument

(Planck) Low Frequency Instrument