minimally coupled dark energy solver introduced in camb for later use in cosmomc

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Nishanth Sasankan
Posts: 27
Joined: September 07 2016
Affiliation: university of notre dame

minimally coupled dark energy solver introduced in camb for

Post by Nishanth Sasankan » November 21 2017

I wanted to introduce a dynamic dark energy into CAMB, I am aware of the ppf module for camb, but i wanted to make a (w,a) array inside camb for future use with cosmomc, so I made a few modifications to camb.

Firstly I made a quintessence [V(phi)] solver which makes a table of a, wde and Grho(phi).
Then I set w not equal to −1 in params.ini (set it to w=−0.9)
Secondly i changed the subroutine DarkEnergy_ReadParams(Ini) to assign appf= a
wppf=wde
from the array generated by my solver
also nwmax and nw_ppf was changed to the length of the array made by the quintessence solver.

Thirdly subroutine wde(a) and grho_de(a) was changed to interpolate using the wde and Grho(phi) made by the solver.

after all this change I compiled camb and ran with params.ini

Now I have the following problem

Warning: Rombint failed to converge;
integral, error, tol: −3.0072019202949190E−003 1.6916759104756984E−005 1.0000000000000001E−005
Warning: Rombint failed to converge;
integral, error, tol: 7.8744507688045242E−004 1.2436078607158674E−004 1.0000000000000001E−005
Warning: Rombint failed to converge;
integral, error, tol: −1.2034253426512022E−004 4.5098700359258448E−004 1.0000000000000001E−005

after the first 2 warning lines, the code takes so much time before the 3rd warning comes up and then the following results show up

Reion redshift = 10.533
Om_b h2 = 0.022600
Om_c h2 = 0.112000
Om_nu h2 = 0.000640
Om_Lambda = 0.724000
Om_K = 0.000000
Om_m (1-Om_K-Om_L) = 0.276000
100 theta (CosmoMC) = 0.929762
Neff (total) = 3.046000
1 nu, g= 1.0153 mnu*c2/k_B/T_nu0= 353.71 (mnu= 0.060 eV)
Reion opt depth = 0.0900
Age of universe/GYr = 17.719
zstar = 1088.72
rs(zstar)/Mpc = 146.38
100*theta = 0.930036
DA(zstar)/Gpc = 15.73971
zdrag = 1059.70
rs(zdrag)/Mpc = 149.01
kD(zstar) Mpc = 0.1392
100*theta_D = 0.143348
zEQ (if vnu=1) = 3216.47
kEQ Mpc (if vnu=1) = 0.009817
100*theta_EQ = 0.758217
100*theta_rs_EQ = 0.417776
tau_recomb/Mpc = 284.76 tau_now/Mpc = 16024.4

Then it takes a very large amount of time to end.
Is there a mistake with the way I modified the code? or is this normal

Thanks
Nishanth

Post Reply