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.
To set Priors on meffstetile in cosmomc 2014
-
- Posts: 1
- Joined: April 08 2014
- Affiliation: University of Rome, Tor Vergata
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: To set Priors on meffstetile in cosmomc 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.