COSMOMC: Prior on \Omega_{DE}>0.1

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Ewan Tarrant
Posts: 1
Joined: November 22 2010
Affiliation: University of Sussex

COSMOMC: Prior on \\Omega_{DE}>0.1

Post by Ewan Tarrant » November 22 2010

Dear All,

I am studying a set of quintessence cosmologies where I struggle to find solutions to the background equations if I demand \Omega_{DE}<~0.1 today. (For particular parameter values in the quintessence potentials I can find no initial conditions for the quintessence field and initial CDM density if I demand \Omega_{DE}<~0.1 today.)

I therefore cannot call CAMB from COSMOMC with \Omega_{DE}~<0.1

To get around the issue I have modified subroutine ParamsToCBMParams (params_CMB.f90):

Code: Select all

try_b = 50
call SetForH&#40;Params,CMB,try_b, .true.&#41;
if &#40;CMB%omv >= 0.1&#41; then
   D_b_OK = .true.
   D_b = CMBToTheta&#40;CMB&#41;
end if

try_t = 100
call SetForH&#40;Params,CMB,try_t&#41;
if &#40;CMB%omv >=0.1 &#41; then
   D_t_OK = .true.
   D_t = CMBToTheta&#40;CMB&#41;
end if

if &#40;D_b_OK==.false. .or. D_t_OK==.false.&#41; then
   CMB%omv=0 ! Reject it

else if &#40;DA < D_b .or. DA > D_t&#41; then
   cmb%H0=0 ! Reject it


and have hardcoded the prior in function GetLogPrior (calclike.f90):

Code: Select all

if &#40;CMB%Omv < 0.1&#41; return
I have also tightened the prior on H_0: 50<H_0<100

Is there anything wrong with this approach?

( I have run this modified code with a standard LCDM cosmology, keeping the default parameter limits in params.ini for COSMOMC, apart from the SZ template normalization which I do not vary. Comparing the results obtained from this run with results obtained from an unmodified COSMOMC, all is as I would expect.)


Ewan Tarrant

Post Reply