Difference between revisions of "Cosmological Parameters"
(49 intermediate revisions by 7 users not shown) | |||
Line 1: | Line 1: | ||
+ | {{DISPLAYTITLE: 2018 Cosmological parameters and MC chains}} | ||
== Description == | == Description == | ||
− | The 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 | + | 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 {{PlanckPapers|planck2013-p11}}, {{PlanckPapers|planck2014-a15||Planck-2015-A15}}, and the 2018 parameter paper {{PlanckPapers|planck2016-l06||Planck-2018-L06}}. |
==Production process== | ==Production process== | ||
+ | |||
+ | Parameter chains are produced using CosmoMC, a sampling package available [http://cosmologist.info/cosmomc 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 [http://getdist.readthedocs.io/en/latest/ 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 == | ||
+ | |||
+ | # 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. | ||
+ | # 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. | ||
+ | # 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. | ||
+ | # Quoted values for the helium abundance (<i>Y</i><sub>P</sub> and <i>Y</i><sub>P</sub><sup>BBN</sup>) and D/H, when <i>Y</i><sub>P</sub> is not varied as a parameter, are obtained assuming the Parthenope 1.1 BBN code, and do not include theoretical errors. | ||
+ | # Some nuisance parameters have slightly different definitions in the plikHM and CamSpecHM likelihoods, and cannot be directly compared. | ||
+ | |||
+ | == Related products == | ||
+ | |||
+ | Results of the baseline parameter exploration runs should be reproducible using CosmoMC with the Planck 2018 likelihood codes. | ||
+ | A new Python code called [https://github.com/JesusTorrado/cobaya Cobaya] also follows the same methodology and should produce identical results. | ||
+ | |||
+ | == Parameter tables == | ||
+ | |||
+ | These tables list parameter constraints for each considered model and data combination separately. For the baseline likelihood combinations see: | ||
+ | |||
+ | * PDF tables with 68% limits [[Media:baseline_params_table_2018_68pc_v2.pdf]] ; | ||
+ | * PDF tables with 95% limits [[Media:baseline_params_table_2018_95pc_v2.pdf]] . | ||
+ | |||
+ | There are also larger full files including results from the alternative CamSpec likelihood, lensing only, other data combinations and priors: | ||
+ | |||
+ | * PDF tables with 68% limits [[Media:result_table_2018_68pc_v2.pdf]] ; | ||
+ | * PDF tables with 95% limits [[Media:Result_table_2018_95pc_v2.pdf]] . | ||
+ | |||
+ | 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 {{PlanckPapers|planck2018-p06}} for full description and references). | ||
+ | |||
+ | |||
+ | {| class="wikitable" align="center" style="text-align:left" border="1" cellpadding="3" cellspacing="0" width=800px | ||
+ | |+ | ||
+ | |- bgcolor="ffdead" | ||
+ | ! 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 <i>z</i><sub>re</sub> > 6.5 | ||
+ | |- | ||
+ | | '''reion''' || A hard prior <i>z</i><sub>re</sub> > 6.5, combined with Gaussian prior <i>z</i><sub>re</sub> = 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), <i>H</i><sub>0</sub> = 73.45 ± 1.66 | ||
+ | |- | ||
+ | | '''BK15''' || Bicep-Keck (+Planck/WMAP) 2015 analysis (arXiv:1810.05216) | ||
+ | |- | ||
+ | | '''theta''' || θ<sub>MC</sub> = 1.0409 ± 0.0006 Gaussian prior | ||
+ | |- | ||
+ | | '''WMAP''' || The full WMAP (temperature and polarization) 9-year data | ||
+ | |- | ||
+ | | '''lenspriors''' || Standard base parameters with <i>n</i><sub>s</sub> = 0.96 ± 0.02, Ω<sub>b</sub><i>h</i><sup>2</sup> = 0.0222 ± 0.0005, 100><i>H</i><sub>0</sub>>40, τ=0.055 | ||
+ | |- | ||
+ | | '''DESpriors''' || DES cosmological parameter priors (flat on 0.1< Ω<sub>m</sub><0.9, 0.03<Ω<sub>b</sub><0.07, 55<<i>H</i><sub>0</sub><91, 0.5<10<sup>9</sup><i>A</i><sub>s</sub><5, <i>Y</i><sub>P</sub>=0.245341 and, if varied, 0.05 eV< Σ <i>m</i><sub>ν</sub> <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 == | ||
+ | |||
+ | 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: | ||
+ | * [https://pla.esac.esa.int/pla/aio/product-action?COSMOLOGY.FILE_ID=COM_CosmoParams_fullGrid_R3.01.zip COM_CosmoParams_fullGrid_R3.01.zip] . | ||
+ | You can also download smaller files containing key results: | ||
+ | * [https://pla.esac.esa.int/pla/aio/product-action?COSMOLOGY.FILE_ID=COM_CosmoParams_base-plikHM-TTTEEE-lowl-lowE_R3.00.zip COM_CosmoParams_base-plikHM-TTTEEE-lowl-lowE_R3.00.zip] – baseline LCDM chains with baseline likelihoods (based on plikHM_TTTEEE_lowl_lowE) ; | ||
+ | * [https://pla.esac.esa.int/pla/aio/product-action?COSMOLOGY.FILE_ID=COM_CosmoParams_base-plikHM_R3.01.zip COM_CosmoParams_base-plikHM_R3.01.zip] – baseline LCDM chains with all plikHM combinations and no <i>z</i><sub>re</sub> prior cut ; | ||
+ | * [https://pla.esac.esa.int/pla/aio/product-action?COSMOLOGY.FILE_ID=COM_CosmoParams_base-plikHM-zre6p5_R3.01.zip COM_CosmoParams_base-plikHM-zre6p5_R3.01.zip] – baseline LCDM chains with all plikHM combinations with <i>z</i><sub>re</sub>>6.5 prior cut ; | ||
+ | * [https://pla.esac.esa.int/pla/aio/product-action?COSMOLOGY.FILE_ID=COM_CosmoParams_lensonly_R3.01.zip COM_CosmoParams_lensonly_R3.01.zip] – CMB lensing + priors based constraints for LCDM and LCDM+mnu ; | ||
+ | * [https://pla.esac.esa.int/pla/aio/product-action?COSMOLOGY.FILE_ID=COM_CosmoParams_parameter-tables_R3.01.zip COM_CosmoParams_parameter-tables_R3.01.zip] – set of PDF parameter tables, including multipole cuts and a .tex file of σ<sub>8</sub> Ω<sub>m</sub><sup>x</sup> data constraints from the baseline likelihood combination. | ||
+ | |||
+ | 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. | ||
+ | |||
+ | |||
+ | {| class="wikitable" align="center" style="text-align:left" border="1" cellpadding="3" cellspacing="0" width=800px | ||
+ | |+ | ||
+ | |- bgcolor="ffdead" | ||
+ | ! 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. | ||
+ | |||
+ | |||
+ | {| class="wikitable" align="center" style="text-align:left" border="1" cellpadding="3" cellspacing="0" width=800px | ||
+ | |+ | ||
+ | |- bgcolor="ffdead" | ||
+ | ! 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 [http://cosmologist.info/cosmomc/readme_planck.html]. The config directory in the download includes information about the grid configuration used by the plotting and grid scripts. | ||
+ | |||
+ | == File formats == | ||
+ | |||
+ | The file formats are standard CosmoMC/GetDist outputs. [http://getdist.readthedocs.io/en/latest/ 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 ℓ, <i>D<sub>ℓ</sub><sup>TT</sup></i>, <i>D<sub>ℓ</sub><sup>TE</sup></i>, <i>D<sub>ℓ</sub><sup>EE</sup></i>, <i>D<sub>ℓ</sub><sup>BB</sup></i>, and <i>D<sub>ℓ</sub><sup>dd</sup></i>, where <i>D<sub>ℓ</sub></i> ≡ ℓ(ℓ+1)<i>C</i><sub>ℓ</sub>/(2π) in μK<sup>2</sup>. Also <i>D<sub>ℓ</sub><sup>dd</sup></i>= [ℓ(ℓ+1)]<sup>2</sup> <i>C</i><sub>ℓ</sub><sup>φφ</sup>/(2π) is the power spectrum of the lensing deflection angle, where <i>C</i><sub>ℓ</sub><sup>φφ</sup> 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 == | ||
+ | |||
+ | <div class="toccolours mw-collapsible mw-collapsed" style="background-color: #FFDAB9;width:80%"> | ||
+ | '''2015 Release of Cosmological Parameters and MC Chains''' | ||
+ | <div class="mw-collapsible-content"> | ||
+ | |||
+ | ''' Description ''' | ||
+ | |||
+ | The 2015 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 reference are contained in {{PlanckPapers|planck2013-p11}} {{PlanckPapers|planck2014-a15||Planck-2015-A15}}. | ||
+ | |||
+ | '''Production process''' | ||
Parameter chains are produced using CosmoMC, a sampling package available [http://cosmologist.info/cosmomc here]. This includes the sample analysis package (and GUI) GetDist, and the scripts for managing, analysing, and plotting results from the full grid or runs. Chain products provided here have had burn in removed. Some results with additional data are produced by importance sampling. | Parameter chains are produced using CosmoMC, a sampling package available [http://cosmologist.info/cosmomc here]. This includes the sample analysis package (and GUI) GetDist, and the scripts for managing, analysing, and plotting results from the full grid or runs. Chain products provided here have had burn in removed. Some results with additional data are produced by importance sampling. | ||
Line 9: | Line 201: | ||
Note that the baseline model includes one massive neutrino (0.06eV). Grid outputs include WMAP 9 results for consistent assumptions. | Note that the baseline model includes one massive neutrino (0.06eV). Grid outputs include WMAP 9 results for consistent assumptions. | ||
− | + | ''' Caveats and known issues ''' | |
# Confidence intervals are derived from the MCMC samples, and assume 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. | # Confidence intervals are derived from the MCMC samples, and assume 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. | ||
# Non-linear lensing modelling uses Halofit; for some extended models and CMB lensing only analyses, tails of the chains may be away from the domain of validity. | # Non-linear lensing modelling uses Halofit; for some extended models and CMB lensing only analyses, tails of the chains may be away from the domain of validity. | ||
# The CAMB version used for most results is Dec 2014; the Jan 2015 version is used for lensing-only models with neutrinos, and only differs in the neutrino corrections to the Halofit model. | # The CAMB version used for most results is Dec 2014; the Jan 2015 version is used for lensing-only models with neutrinos, and only differs in the neutrino corrections to the Halofit model. | ||
− | # There is evidence of temperature-polarization leakage that may affect results including high-L polarization | + | # There is evidence of temperature-polarization leakage that may affect results including high-L polarization; hence use caution in the interpretation of results including polarization |
# Alternative CamSpec likelihood results in the tables are generated using a slightly older CosmoMC version, with fewer derived parameters and a slightly different BBN predictions for the helium abundance. | # Alternative CamSpec likelihood results in the tables are generated using a slightly older CosmoMC version, with fewer derived parameters and a slightly different BBN predictions for the helium abundance. | ||
− | + | ''' Related products ''' | |
− | Results of the parameter exploration runs should be reproducible using CosmoMC with the Planck likelihood | + | Results of the baseline parameter exploration runs should be reproducible using CosmoMC with the Planck 2015 likelihood codes (when available). |
− | + | ''' Parameter Tables ''' | |
These list parameter constraints for each considered model and data combination separately. For the baseline likelihood see | These list parameter constraints for each considered model and data combination separately. For the baseline likelihood see | ||
− | * PDF tables with 68% limits [[ | + | * PDF tables with 68% limits [[Media:baseline_params_table_2015_limit68.pdf]] |
− | * PDF tables with 95% limits [[ | + | * PDF tables with 95% limits [[Media:baseline_params_table_2015_limit95.pdf]] |
There are also larger files including alternative CamSpec and DetSet likelihood results, along with shifts in parameters compared to baseline in units of the baseline error: | There are also larger files including alternative CamSpec and DetSet likelihood results, along with shifts in parameters compared to baseline in units of the baseline error: | ||
− | * PDF tables with 68% limits [[ | + | * PDF tables with 68% limits [[Media:params_table_2015_limit68.pdf]] |
− | * PDF tables with 95% limits [[ | + | * PDF tables with 95% limits [[Media:params_table_2015_limit95.pdf]] |
Line 61: | Line 253: | ||
|- | |- | ||
| '''lensingonly''' || Planck lensing power spectrum reconstruction only; T,E fixed to best-fit spectrum + priors | | '''lensingonly''' || Planck lensing power spectrum reconstruction only; T,E fixed to best-fit spectrum + priors | ||
+ | |- | ||
+ | | '''BKP''' || The Bicep-Keck-Planck fiducial B mode likelihood | ||
|- | |- | ||
| '''zre6p5''' || A hard prior z_re > 6.5 | | '''zre6p5''' || A hard prior z_re > 6.5 | ||
Line 84: | Line 278: | ||
− | Tags used to identify the model parameters that are varied are described in [[ | + | Tags used to identify the model parameters that are varied are described in [[Media:parameter_tag_definitions_2015.pdf]]. |
− | + | ''' Parameter Chains ''' | |
We provide the full chains and getdist outputs for our parameter results. The entire grid of results is available from as a 3.7GB compressed file: | We provide the full chains and getdist outputs for our parameter results. The entire grid of results is available from as a 3.7GB compressed file: | ||
− | * {{PLASingleFile|fileType=cosmo|name= | + | <!--- * {{PLASingleFile|fileType=cosmo|name= COM_CosmoParams_R2.00.tar.gz|link=Full Grid Download}} ----> |
+ | * ''COM_CosmoParams_R2.nn.tar.gz'' | ||
+ | Where ''nn'' is the most recent update. You can also download the Bicep2/Keck/Planck (BKP) joint constraints for +r models, and smaller files containing key results in the base model only: | ||
+ | <!---* {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_plikHM_TT_lowTEB_R2.00.tar.gz|link=Baseline LCDM chains with plikHM_TT_lowTEB}} | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_plikHM_TT_lowTEB_R2.00.tar.gz|link=Baseline LCDM chains with all plikHM combinations}} | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_lensonly_R2.00.tar.gz|link=CMB lensing only in LCDM}} ---> | ||
+ | * ''COM_CosmoParams_base_plikHM_TT_lowTEB_R2.nn.tar.gz'' | ||
+ | * ''COM_CosmoParams_base_plikHM_R2.nn.tar.gz'' | ||
+ | * ''COM_CosmoParams_base_lensonly_R2.nn.tar.gz'' | ||
+ | * ''COM_CosmoParams_base_r_plikHM_BKP_R2.nn.tar.gz'' | ||
+ | |||
+ | 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 | ||
+ | |||
+ | |||
+ | {| class="wikitable" align="center" style="text-align:left" border="1" cellpadding="3" cellspacing="0" width=800px | ||
+ | |+ | ||
+ | |- bgcolor="ffdead" | ||
+ | ! 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''' || The best-fit temperature and polarization power spectra and lensing potential (see below) | ||
+ | |- | ||
+ | | '''.minimum.plik_foregrounds''' || The 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'' | |
− | * {{PLASingleFile|fileType=cosmo|name= | + | |
− | * {{PLASingleFile|fileType=cosmo|name= | + | 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 |
− | * {{PLASingleFile|fileType=cosmo|name= | + | |
+ | |||
+ | {| class="wikitable" align="center" style="text-align:left" border="1" cellpadding="3" cellspacing="0" width=800px | ||
+ | |+ | ||
+ | |- bgcolor="ffdead" | ||
+ | ! 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''' || A 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 [http://cosmologist.info/cosmomc/readme_planck.html]. The config directory in the download includes information about the grid configuration used by the plotting and grid scripts. | ||
+ | |||
+ | ''' File formats ''' | ||
+ | |||
+ | The file formats are standard CosmoMC/GetDist outputs. [http://getdist.readthedocs.io/en/latest/ 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. Note that 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%, 3: 99%. The limit tags specify whether a given limit is one tail, two tail 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: <math>l</math>, <math>D^{TT}_l</math>, <math>D^{TE}_l</math>, <math>D^{EE}_l</math>, <math>D^{BB}_l</math>, and <math>D^{dd}_l</math>, were <math>D_l \equiv l(l+1) C_l / (2\pi)</math> in <math>\mu{\rm K}^2</math>. Also <math>D^{dd}_l= [l(l+1)]^2 C^{\phi\phi}_l/(2\pi)</math> is the power spectrum of the lensing deflection angle, where <math>C^{\phi\phi}_l</math> is the lensing potential power spectrum. Note that the lensing spectrum may not be accurate at L > 400 due to the maximum wavenumber and non-linear correction accuracy settings. | ||
+ | |||
+ | </div> | ||
+ | </div> | ||
+ | |||
+ | <div class="toccolours mw-collapsible mw-collapsed" style="background-color: #EEE8AA;width:80%"> | ||
+ | '''2013 Release of Cosmological parameters and MC Chains''' | ||
+ | <div class="mw-collapsible-content"> | ||
+ | |||
+ | ''' Description ''' | ||
+ | |||
+ | The 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 reference are contained in {{PlanckPapers|planck2013-p11}}. | ||
+ | |||
+ | '''Production process''' | ||
+ | |||
+ | Parameter chains are produced using CosmoMC, a sampling package available [http://cosmologist.info/cosmomc here]. This includes the sample analysis package GetDist, and the scripts for managing, analysing, and plotting results from the full grid or runs. 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 neutrino (0.06eV). Grid outputs include WMAP 9 results for consistent assumptions. | ||
+ | |||
+ | ''' Caveats and known issues ''' | ||
+ | |||
+ | # Confidence intervals are derived from the MCMC samples, and assume 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. We had some issues producing reliable results from the minimizer used to produce the best fits, so in some cases the quoted fits may be significantly improved. The chain outputs contain some parameters that are not used, for example the beam mode ranges for all but the first mode (the beam modes are marginalised over anlaytically internally to the likelihood). | ||
+ | # Where determined from BBN consistency, the <math>Y_P</math> parameter uses an interpolation table from{{BibCite|Hamann2007sb}} based on the 2008 version of the Parthenope BBN code. More recent updates to the neutron lifetime suggest that the <math>Y_P</math> values reported in the tables may be in error by around 0.0005. This has a negligible impact on the predicted CMB power spectrum or any of the parameter results reported in this series of papers. However, the difference should be taken into account when comparing with BBN results reported in Sect. 6.4. of {{PlanckPapers|planck2013-p11}}, which use an updated version for the neutron lifetime (and several other nuclear reaction rates that have negligible impact). Note also that the error on <math>Y_P</math> quoted in the tables here does not include theoretical errors in the BBN prediction. | ||
+ | |||
+ | ''' Related products ''' | ||
+ | |||
+ | Results of the parameter exploration runs should be reproducible using CosmoMC with the Planck likelihood code. | ||
+ | |||
+ | ''' Parameter Tables ''' | ||
+ | |||
+ | These list paramter constraints for each considered model and data combination separately | ||
+ | |||
+ | * PDF tables with 68% limits [[File:grid_limit68.pdf]] | ||
+ | * PDF tables with 95% limits [[File:grid_limit95.pdf]] | ||
+ | |||
+ | There are also summary comparison tables, showing how constraints for selected models vary with data used to constrain them: | ||
+ | |||
+ | * Comparison tables with 68% limits [[File:comparetables_limit68.pdf | bla]] | ||
+ | * Comparison tables with 95% limits [[File:comparetables_limit95.pdf]] | ||
+ | |||
+ | |||
+ | Data combination tags used to label results are as follows (see {{PlanckPapers|planck2013-p11}} for full description and references): | ||
+ | |||
+ | |||
+ | {| class="wikitable" align="center" style="text-align:left" border="1" cellpadding="3" cellspacing="0" width=800px | ||
+ | |+ | ||
+ | |- bgcolor="ffdead" | ||
+ | ! Tag|| Data | ||
+ | |- | ||
+ | | '''planck''' || high-L Planck temperature (CamSpec, 50 <= l <= 2500) | ||
+ | |- | ||
+ | | '''lowl''' || low-L: Planck temperature (2 <= l <= 49) | ||
+ | |- | ||
+ | | '''lensing''' || Planck lensing power spectrum reconstruction | ||
+ | |- | ||
+ | | '''lowLike''' || low-L WMAP 9 polarization (WP) | ||
+ | |- | ||
+ | | '''tauprior''' || A Gaussian prior on the optical depth, tau = 0.09 +- 0.013 | ||
+ | |- | ||
+ | | '''BAO''' || Baryon oscillation data from DR7, DR9 and and 6DF | ||
+ | |- | ||
+ | | '''SNLS''' || Supernova data from the Supernova Legacy Survey | ||
+ | |- | ||
+ | | '''Union2''' || Supernova data from the Union compilation | ||
+ | |- | ||
+ | | '''HST''' || Hubble parameter constraint from HST (Riess et al) | ||
+ | |- | ||
+ | | '''WMAP''' || The full WMAP (temperature and polarization) 9 year data | ||
+ | |} | ||
+ | |||
+ | |||
+ | Tags used to identify the model paramters that are varied are described in [[File:parameter_tag_definitions.pdf]]. Note that alpha1 results are not used in the parameter paper, and are separate from the isocurvature results in the inflation paper. | ||
+ | |||
+ | ''' Parameter Chains ''' | ||
+ | |||
+ | We provide the full chains and getdist outputs for our parameter results. The entire grid of results is available from as a 2.8GB compressed file: | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_FullGrid_R1.10.tar.gz|link=Full Grid Download}} | ||
+ | |||
+ | You can also download key chains for the baseline LCDM model here: | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_planck_lowl_post_lensing_R1.10.tar.gz|link=Planck+lensing}} | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_planck_lowl_lowLike_R1.10.tar.gz|link=Planck+WP}} | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_planck_lowl_lowLike_post_lensing_R1.10.tar.gz|link=Planck+WP+lensing}} | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_planck_lowl_lowLike_highL_R1.10.tar.gz|link=Planck+WP+highL}} | ||
+ | * {{PLASingleFile|fileType=cosmo|name=COM_CosmoParams_base_planck_lowl_lowLike_highL_post_lensing_R1.10.tar.gz|link=Planck+WP+highL+lensing}} | ||
The download contains a hierarchy of directories, with each separate chain in a separate directory. The structure for the directories is | The download contains a hierarchy of directories, with each separate chain in a separate directory. The structure for the directories is | ||
Line 128: | Line 490: | ||
− | In addition each directory contains any | + | In addition each directory contains any importanced sampled outputs with additional data. These have names of the form |
: ''base_AAA_BBB_XXX_YYY_post_ZZZ.ext'' | : ''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 | + | 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 convntions, with the following extensions |
Line 152: | Line 514: | ||
− | Python scripts for reading in chains and calculating new derived parameter constraints are available as part of CosmoMC, see the readme for details [http://cosmologist.info/cosmomc/readme_planck.html] | + | Python scripts for reading in chains and calculating new derived parameter constraints are available as part of CosmoMC, see the readme for details [http://cosmologist.info/cosmomc/readme_planck.html]. |
− | + | ''' File formats ''' | |
− | The file formats are standard | + | The file formats are standard March 2013 CosmoMC outputs. CosmoMC includes python scripts for generating tables, 1D, 2D and 3D plots using the provided data. The formats are summarised here: |
; Chain files | ; Chain files | ||
Line 163: | Line 525: | ||
: '' weight like param1 param2 param3 …'' | : '' 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. | + | : 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. Note that 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. | : 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. Note that 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. | ||
Line 175: | Line 537: | ||
;.bestfit_cl files | ;.bestfit_cl files | ||
− | : They contain the best-fit theoretical power spectra (without foregrounds) for each model. The columns are: <math>l</math>, <math>D^{TT}_l</math>, <math>D^{TE}_l</math>, <math>D^{EE}_l</math>, <math>D^{BB}_l</math>, and <math>D^{dd}_l</math>, were <math>D_l \equiv l(l+1) C_l / (2\pi)</math> in <math>\mu{\rm K}^2</math>. Also <math>D^{dd}_l= [l(l+1)]^2 C^{\phi\phi}_l/(2\pi)</math> is the power spectrum of the lensing deflection angle, where <math>C^{\phi\phi}_l</math> is the lensing potential power spectrum. | + | : They contain the best-fit theoretical power spectra (without foregrounds) for each model. The columns are: <math>l</math>, <math>D^{TT}_l</math>, <math>D^{TE}_l</math>, <math>D^{EE}_l</math>, <math>D^{BB}_l</math>, and <math>D^{dd}_l</math>, were <math>D_l \equiv l(l+1) C_l / (2\pi)</math> in <math>\mu{\rm K}^2</math>. Also <math>D^{dd}_l= [l(l+1)]^2 C^{\phi\phi}_l/(2\pi)</math> is the power spectrum of the lensing deflection angle, where <math>C^{\phi\phi}_l</math> is the lensing potential power spectrum. For results not including the lensing likelihood, this is the prediction from linear theory; for lensing outputs this includes corrections due to non-linear structure growth. The <math>D_l</math> are output to high <math>l</math>, but not actually computed above <math>l_{\rm max}=2500</math> (Planck), <math>l_{\rm max}=4500</math> (Planck+highL) or <math>l_{\rm max}=1500</math> (WMAP), and <math>l</math> values above these are fixed to a scaled fiducial template. |
+ | |||
+ | |||
+ | </div> | ||
+ | </div> | ||
+ | |||
− | |||
− | |||
== References == | == References == |
Latest revision as of 06:43, 31 July 2019
Contents
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]
- 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.
- 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.
- 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.
- 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.
- 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:
- PDF tables with 68% limits Media:baseline_params_table_2018_68pc_v2.pdf ;
- PDF tables with 95% limits Media:baseline_params_table_2018_95pc_v2.pdf .
There are also larger full files including results from the alternative CamSpec likelihood, lensing only, other data combinations and priors:
- PDF tables with 68% limits Media:result_table_2018_68pc_v2.pdf ;
- PDF tables with 95% limits Media:Result_table_2018_95pc_v2.pdf .
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 [4] 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:
- COM_CosmoParams_base-plikHM-TTTEEE-lowl-lowE_R3.00.zip – baseline LCDM chains with baseline likelihoods (based on plikHM_TTTEEE_lowl_lowE) ;
- COM_CosmoParams_base-plikHM_R3.01.zip – baseline LCDM chains with all plikHM combinations and no zre prior cut ;
- COM_CosmoParams_base-plikHM-zre6p5_R3.01.zip – baseline LCDM chains with all plikHM combinations with zre>6.5 prior cut ;
- COM_CosmoParams_lensonly_R3.01.zip – CMB lensing + priors based constraints for LCDM and LCDM+mnu ;
- COM_CosmoParams_parameter-tables_R3.01.zip – set of PDF parameter tables, including multipole cuts and a .tex file of σ8 Ωmx data constraints from the baseline likelihood combination.
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 ℓ, DℓTT, DℓTE, DℓEE, DℓBB, and Dℓdd, where Dℓ ≡ ℓ(ℓ+1)Cℓ/(2π) in μK2. Also Dℓdd= [ℓ(ℓ+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]
2015 Release of Cosmological Parameters and MC Chains
Description
The 2015 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 reference are contained in Planck-2013-XVI[1]Planck-2015-A15[2].
Production process
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 or runs. 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 neutrino (0.06eV). Grid outputs include WMAP 9 results for consistent assumptions.
Caveats and known issues
- Confidence intervals are derived from the MCMC samples, and assume 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.
- Non-linear lensing modelling uses Halofit; for some extended models and CMB lensing only analyses, tails of the chains may be away from the domain of validity.
- The CAMB version used for most results is Dec 2014; the Jan 2015 version is used for lensing-only models with neutrinos, and only differs in the neutrino corrections to the Halofit model.
- There is evidence of temperature-polarization leakage that may affect results including high-L polarization; hence use caution in the interpretation of results including polarization
- Alternative CamSpec likelihood results in the tables are generated using a slightly older CosmoMC version, with fewer derived parameters and a slightly different BBN predictions for the helium abundance.
Related products
Results of the baseline parameter exploration runs should be reproducible using CosmoMC with the Planck 2015 likelihood codes (when available).
Parameter Tables
These list parameter constraints for each considered model and data combination separately. For the baseline likelihood see
- PDF tables with 68% limits Media:baseline_params_table_2015_limit68.pdf
- PDF tables with 95% limits Media:baseline_params_table_2015_limit95.pdf
There are also larger files including alternative CamSpec and DetSet likelihood results, along with shifts in parameters compared to baseline in units of the baseline error:
- PDF tables with 68% limits Media:params_table_2015_limit68.pdf
- PDF tables with 95% limits Media:params_table_2015_limit95.pdf
Data combination tags used to label results are as follows (see Planck-2013-XVI[1] for full description and references):
Tag | Data |
---|---|
plikHM | baseline high-L Planck power spectra (plik cross half-mission, 30 <= l <= 2508) |
plikDS | high-L Planck power spectra (plik cross detsets, 30 <= l <= 2508) |
CamSpecHM | high-L Planck power spectra (CamSpec cross half-mission, 30 <= l <= 2500) |
CamSpecDS | high-L Planck power spectra (CamSpec cross detsets, 30 <= l <= 2500) |
lowl | low-L: Planck temperature only (2 <= l <= 29) |
lowTEB | low-L temperature and LFI polarization (2 <= l <= 29) |
lowEB | low-L LFI polarization only (2 <= l <= 29) |
WMAPTEB | low-L temperature, and LFI+WMAP polarization (2 <= l <= 29) |
lensing | Planck lensing power spectrum reconstruction |
lensingonly | Planck lensing power spectrum reconstruction only; T,E fixed to best-fit spectrum + priors |
BKP | The Bicep-Keck-Planck fiducial B mode likelihood |
zre6p5 | A hard prior z_re > 6.5 |
tau07 | A Gaussian prior on the optical depth, tau = 0.07 +- 0.02 |
reion | A hard prior z_re > 6.5, combined with Gaussian prior z_re = 7 +- 1 |
BAO | Baryon oscillation data from DR11LOWZ, DR11CMASS, MGS and 6DF |
JLA | Supernova data from the SDSS-II/SNLS3 Joint Light-curve Analysis |
H070p6 | Hubble parameter constraint, H_0 = 70.6 +- 3.3 |
theta | theta_MC fixed to 1.0408 |
WLonlyHeymans | Conservative cut of the CFHTLenS weak lensing data + priors |
WMAP | The full WMAP (temperature and polarization) 9 year data |
The high-L Planck likelihoods have TT, TE, EE variants from each spectrum alone, plus the TTTEEE joint constraint.
Tags used to identify the model parameters that are varied are described in Media:parameter_tag_definitions_2015.pdf.
Parameter Chains
We provide the full chains and getdist outputs for our parameter results. The entire grid of results is available from as a 3.7GB compressed file:
- COM_CosmoParams_R2.nn.tar.gz
Where nn is the most recent update. You can also download the Bicep2/Keck/Planck (BKP) joint constraints for +r models, and smaller files containing key results in the base model only:
- COM_CosmoParams_base_plikHM_TT_lowTEB_R2.nn.tar.gz
- COM_CosmoParams_base_plikHM_R2.nn.tar.gz
- COM_CosmoParams_base_lensonly_R2.nn.tar.gz
- COM_CosmoParams_base_r_plikHM_BKP_R2.nn.tar.gz
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
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 | The best-fit temperature and polarization power spectra and lensing potential (see below) |
.minimum.plik_foregrounds | The 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 | A 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 [2]. The config directory in the download includes information about the grid configuration used by the plotting and grid scripts.
File formats
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. Note that 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%, 3: 99%. The limit tags specify whether a given limit is one tail, two tail 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: , , , , , and , were in . Also is the power spectrum of the lensing deflection angle, where is the lensing potential power spectrum. Note that the lensing spectrum may not be accurate at L > 400 due to the maximum wavenumber and non-linear correction accuracy settings.
2013 Release of Cosmological parameters and MC Chains
Description
The 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 reference are contained in Planck-2013-XVI[1].
Production process
Parameter chains are produced using CosmoMC, a sampling package available here. This includes the sample analysis package GetDist, and the scripts for managing, analysing, and plotting results from the full grid or runs. 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 neutrino (0.06eV). Grid outputs include WMAP 9 results for consistent assumptions.
Caveats and known issues
- Confidence intervals are derived from the MCMC samples, and assume 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. We had some issues producing reliable results from the minimizer used to produce the best fits, so in some cases the quoted fits may be significantly improved. The chain outputs contain some parameters that are not used, for example the beam mode ranges for all but the first mode (the beam modes are marginalised over anlaytically internally to the likelihood).
- Where determined from BBN consistency, the [5] based on the 2008 version of the Parthenope BBN code. More recent updates to the neutron lifetime suggest that the values reported in the tables may be in error by around 0.0005. This has a negligible impact on the predicted CMB power spectrum or any of the parameter results reported in this series of papers. However, the difference should be taken into account when comparing with BBN results reported in Sect. 6.4. of Planck-2013-XVI[1], which use an updated version for the neutron lifetime (and several other nuclear reaction rates that have negligible impact). Note also that the error on quoted in the tables here does not include theoretical errors in the BBN prediction. parameter uses an interpolation table from
Related products
Results of the parameter exploration runs should be reproducible using CosmoMC with the Planck likelihood code.
Parameter Tables
These list paramter constraints for each considered model and data combination separately
- PDF tables with 68% limits File:Grid limit68.pdf
- PDF tables with 95% limits File:Grid limit95.pdf
There are also summary comparison tables, showing how constraints for selected models vary with data used to constrain them:
- Comparison tables with 68% limits File:Comparetables limit68.pdf
- Comparison tables with 95% limits File:Comparetables limit95.pdf
Data combination tags used to label results are as follows (see Planck-2013-XVI[1] for full description and references):
Tag | Data |
---|---|
planck | high-L Planck temperature (CamSpec, 50 <= l <= 2500) |
lowl | low-L: Planck temperature (2 <= l <= 49) |
lensing | Planck lensing power spectrum reconstruction |
lowLike | low-L WMAP 9 polarization (WP) |
tauprior | A Gaussian prior on the optical depth, tau = 0.09 +- 0.013 |
BAO | Baryon oscillation data from DR7, DR9 and and 6DF |
SNLS | Supernova data from the Supernova Legacy Survey |
Union2 | Supernova data from the Union compilation |
HST | Hubble parameter constraint from HST (Riess et al) |
WMAP | The full WMAP (temperature and polarization) 9 year data |
Tags used to identify the model paramters that are varied are described in File:Parameter tag definitions.pdf. Note that alpha1 results are not used in the parameter paper, and are separate from the isocurvature results in the inflation paper.
Parameter Chains
We provide the full chains and getdist outputs for our parameter results. The entire grid of results is available from as a 2.8GB compressed file:
You can also download key chains for the baseline LCDM model here:
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
Extension | Data |
---|---|
.txt | parameter chain file with burn in removed |
.paramnames | File that describes the parameters included in the chains |
.minimum | Best-fit parameter values, -log likelihoods and chi-square |
.bestfit_cl | The best-fit temperature and polarization power spectra and lensing potential (see below) |
.inputparams | Input parameters used when generating the chain |
.minimum.inputparams | Input parameters used when generating the best fit |
.ranges | prior ranges assumed for each parameter |
In addition each directory contains any importanced 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 convntions, 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 |
.corr | Correlation matrix for the parameters |
.converge | A 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 [3].
File formats
The file formats are standard March 2013 CosmoMC outputs. CosmoMC includes python scripts for generating tables, 1D, 2D and 3D plots using the provided data. 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. Note that 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%, 3: 99%. The limit tags specify whether a given limit is one tail, two tail or none (if no constraint within the assumed prior boundary).
- .bestfit_cl files
- They contain the best-fit theoretical power spectra (without foregrounds) for each model. The columns are: , , , , , and , were in . Also is the power spectrum of the lensing deflection angle, where is the lensing potential power spectrum. For results not including the lensing likelihood, this is the prediction from linear theory; for lensing outputs this includes corrections due to non-linear structure growth. The are output to high , but not actually computed above (Planck), (Planck+highL) or (WMAP), and values above these are fixed to a scaled fiducial template.
References[edit]
- ↑ 1.01.11.21.31.41.5 Planck 2013 results. XVI. Cosmological parameters, Planck Collaboration, 2014, A&A, 571, A16.
- ↑ 2.02.1 Planck 2015 results. XIII. Cosmological parameters, Planck Collaboration, 2016, A&A, 594, A13.
- ↑ Planck 2018 results. VI. Cosmological parameters, Planck Collaboration, 2020, A&A, 641, A6.
- ↑
- ↑ 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