CosmoCoffee

 FAQ   Search  SmartFeed   Memberlist    Register Profile   Log in Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

 Constraining a cosmology with two new coupled parameters
Author Message
Albin Nilsson

Joined: 05 Oct 2016
Posts: 13
Affiliation: NCBJ, Warsaw

 Posted: October 17 2016 I'm sorry if this is a trivial question, but I am just getting started with cosmomc and I am a little confused. I have a cosmological model which basically introduces two new parameters into the Friedmann equation, call them param1 and param2. These parameters are not independent, but depend on each other in some way. I would like to constrain these parameters using Planck data. I assume that I need to calculate and predict some parameter values which I can then constrain using cosmomc. Could anyone tell me how to go about doing this? I know that there are a few threads about adding new parameters, but I have not been able to gather how to go about coding my specific case. Thanks in advance
Antony Lewis

Joined: 23 Sep 2004
Posts: 1303
Affiliation: University of Sussex

 Posted: October 17 2016 If param1 is a function of param2, then just add the one parameter param1, and calculate param2 in your theory/likelihood calculation.
Albin Nilsson

Joined: 05 Oct 2016
Posts: 13
Affiliation: NCBJ, Warsaw

 Posted: October 17 2016 By adding you mean following the procedure in: http://cosmocoffee.info/viewtopic.php?t=1149&postdays=0&postorder=asc&start=15? My problem is really that I don't understand how cosmomc then knows what param1 is, and how the equations for cosmological evolution depend on it. In this case, it might be easier to redefine ΩM, ΩΛ, Ωcurv and Ωrad, since my new parameters attach to all of those. Where are they defined?
Antony Lewis

Joined: 23 Sep 2004
Posts: 1303
Affiliation: University of Sussex

 Posted: October 17 2016 It's the "Modify Calculator_CAMB.f90" step of the readme FAQ where you interface the cosmomc variable with CAMB which determines what it actually means.
Albin Nilsson

Joined: 05 Oct 2016
Posts: 13
Affiliation: NCBJ, Warsaw

Posted: October 18 2016

Thanks for the help so far. This is what I have done:

1. Added line 115 in CosmologyTypes.f90, added 2 new parameters, param1 and param2
2. Changed the number of hard parameters on line 71 in CosmologyParameterization
s.f90 from 16 to 18.
3. Changed line 64 to call myparams.paramnames instead of params_CMB.paramnames
4. Created myparams.paramnames and added param1 and param2 at the end
5. Added lines 85 and 86 in Calculator_CAMB.f90 (real(dl) param1, and real(dl) param2)

The last step is then to modify the code to account for the structure changed by these new parameters; this is where I am stumped (partly because I don't speak Fortran).
My parameters attach to Ωm, ΩΛ, Ωrad and Ωcurv in different ways.
If my new Ωm is (for example) Ω = Ωm * param1, can I express that like this in Calculator_CAMB:
 Code: P%omegab = CMB%omb * param1

and
 Code: P%omegac = CMB%omc * param1

also, where do the other ΛCDM-parameters (curv and rad) come in?. I assume Ωm = Ωc + Ωb

Thanks again for answering my questions
 Display posts from previous: All Posts1 Day7 Days2 Weeks1 Month3 Months6 Months1 Year Oldest FirstNewest First
 All times are GMT + 5 Hours Page 1 of 1

 Jump to: Select a forum Arxiv paper discussion----------------arXiv papers Topic discussion----------------Early UniverseCosmological ModelCosmological Observations  Projects and Resources----------------Computers and softwareTeaching, Papers and PresentationsResearch projectsiCosmo Coming up----------------Job vacanciesConferences and meetings Management----------------CosmoCoffee
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum