To set Priors on meffstetile in cosmomc 2014

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Balakrishna sandeep
Posts: 1
Joined: April 08 2014
Affiliation: University of Rome, Tor Vergata

To set Priors on meffstetile in cosmomc 2014

Post by Balakrishna sandeep » July 11 2014

Hi,
I am trying to work on Warm dark matter scenario with sterile neutrinos as WDM candidates using COSMOMC 2014, sterile neutrino mass range in KeV.

Priors used :

param[omegach2] = 0.01 0.001 0.99 0.001 0.002
param[nnu] = 3.4 3.046 6 0.0001 0.0002
param[omeganuh2] = 0.11 0.001 0.99 0.001 0.002
param[meffsterile] = 100 0 5000 0.1 0.2
param[mnu] = 0.06

I have changed the "sterile_mphys_max" to 5000 in file CosmologyParameterizations.f90, since the units are defined in eV.


I get an error "Couldn't start after 1000 tries - check starting ranges" when trying to run with the above priors.

I believe it is not straight forward to set priors on meffsterile as it is not being passed to CAMB.
If I can, what is the right way to set priors on mass of the sterile neutrino ?

Thank you in advance.

Antony Lewis
Posts: 1941
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: To set Priors on meffstetile in cosmomc 2014

Post by Antony Lewis » July 12 2014

If you want to set a prior on the physical mass (for thermal or DW case), set it along with sterile_mphys_max in CosmologyParameterizations. The default code is only tested for masses less than about 10eV.

Post Reply