I am trying to run Cobaya for a LCDM + Neff + Yp model using some CMB power spectra. While adding the two extension parameters individually produces reasonable results, freeing them simultaneously yields the monstrosity below.
These chains have run for O(10h) on a multi-processor system and have >100k samples, so I don't think convergence is an issue here. Instead, I suspect this is due to some relation between Neff and Yp being imposed that I'm not aware of. In the .yaml file I simply add the two as free parameters ala
Code: Select all
(standard LCDM parameters)
YHe:
prior:
min: 0.05
max: 0.8
ref:
dist: norm
loc: 0.245
scale: 0.01
proposal: 0.001
latex: Y_\mathrm{P}
nnu:
prior:
min: 0
max: 10
ref:
dist: norm
loc: 3.046
scale: 0.5
proposal: 0.001
latex: N_\mathrm{eff}
Code: Select all
bbn_predictor: PArthENoPE_880.2_standard.dat
Code: Select all
bbn_consistency: False
Thanks in advance!