[CosmoMC] Rombint failed to converge (Out of range finding H0)

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Mason Ng
Posts: 25
Joined: May 17 2017
Affiliation: The University of Auckland

[CosmoMC] Rombint failed to converge (Out of range finding H

I'm actually running it with ModeCode, but the error is still from CosmoMC. I set param[w] = 0.33333 (radiation dominated; i.e. just after inflation).

However, I get the error as written in the title. This is probably a basic question - I looked at http://cosmocoffee.info/viewtopic.php?t=2593 which brought up the idea of finding the starting value of theta that commensurates with w -
how can I derive an 'intelligent' start value for theta?

In the radiation-dominated universe, I can write

$D_A = \frac{1}{H_0(1+z)} \int_0^z \frac{dz}{\sqrt{\Omega_r(1+z)^4}}$, where
$a = \frac{1}{1+z} = \sqrt{2H_0\sqrt{\Omega_r}}$; I'm not sure where to integrate to though. Additionally, I'm wondering where the equation of state parameter $w$ comes in. Knowing D_A, I was wondering if there was something I need to consider for getting the sound horizon distance $r_s$. From there, I want to get $\theta$. Thanks in advance.

For now, I have (conservatively)

Code: Select all

H0_min = 1
H0_max = 150
param&#91;theta&#93; = 1.0411 1.00 1.08 0.0004 0.0004
param&#91;w&#93; = 0.33333 &#40;but plan to put a prior on it&#41;. 
[/b]

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

Re: [CosmoMC] Rombint failed to converge (Out of range findi

Theta for any model that fits the CMB data will always be close to the LCDM best-fit value (by construction). Better to use a more sensible value for w (or other degenerate parameters).