COSMOMC: OMEGAV~0.68 ?!

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Minghua Li
Posts: 15
Joined: August 25 2009
Affiliation: Institue of High Energy Physics, Chinese Academy of Sciences

COSMOMC: OMEGAV~0.68 ?!

Post by Minghua Li » October 25 2009

Settings:

Use_SN = T and all other Use_CMB, Use_XX, ... = F
action = 0
Samples_to_get = 20,000


My params.ini file:
#parameter start center, min, max, start width, st. dev. estimate
#ombh2
param1 = 0.0223 0.005 0.1 0.001 0.001
#omdmh2
param2 = 0.105 0.01 0.99 0.01 0.01
#theta
param3 = 1.04 0.5 10 0.002 0.002
#tau
#param4 = 0.09 0.01 0.8 0.03 0.03
param4 = 0.00 0.00 0.0 0.00 0.00


#omk
param5 = 0 0 0 0 0
#nufrac
param6 = 0 0 0 0 0
#w
param7 = -1 -1 -1 0 0

#n_s
param8 = 0.95 0.5 1.5 0.02 0.01
#n_t
param9 = 0 0 0 0 0
#n_run
param10 = 0 0 0 0 0

#log[10^10 A_s]
param11 = 3 2.7 4 0.01 0.01
#amp_ratio
param12 = 0 0 0 0 0
#SZ amplitude, as in WMAP analysis
param13 = 1 0 2 0.4 0.4


Reslults:
param mean sddev lower1 upper1 lower2 upper2
1 0.363199E-01 0.102716E-01 0.255479E-01 0.470955E-01 0.150799E-01 0.528502E-01 \Omega_b h^2
2 0.916378E-01 0.413039E-01 0.545598E-01 0.122560E+00 0.415733E-01 0.209161E+00 \Omega_c h^2
3 0.948984E+00 0.638113E-01 0.877551E+00 0.994789E+00 0.839910E+00 0.110133E+01 \theta
8 0.112733E+01 0.219944E+00 0.910778E+00 0.138105E+01 0.669724E+00 0.148456E+01 n_s
11 0.323718E+01 0.372499E+00 0.270021E+01 0.365387E+01 0.270021E+01 0.393914E+01 log[10^{10} A_s]
13 0.977879E+00 0.578001E+00 0.000000E+00 0.200000E+01 0.000000E+00 0.200000E+01 A_{SZ}
14 0.679480E+00 0.458794E-01 0.633649E+00 0.725960E+00 0.587615E+00 0.762129E+00 Omega_lambda
15 0.152279E+02 0.240025E+01 0.127401E+02 0.177274E+02 0.103971E+02 0.192567E+02 Age/GYr
16 0.320520E+00 0.458794E-01 0.274063E+00 0.366354E+00 0.237894E+00 0.412569E+00 \Omega_m
20 0.626215E+02 0.111770E+02 0.519260E+02 0.731843E+02 0.467391E+02 0.893054E+02 H_0

Limits are: 0.68;0.95
1 two tail; \Omega_b h^2
2 two tail; \Omega_c h^2
3 two tail; \theta
8 two tail; n_s
11 two tail; log[10^{10} A_s]
13 one tail; A_{SZ}
14 two tail; Omega_lambda
15 two tail; Age/GYr
16 two tail; \Omega_m
20 two tail; H_0


Why didn't I get the OMEGAV~0.73 ?!
Is there anything wrong with my params.ini settings? Or the did not have enough samples ?
PS: The dashed line in test_1D.ps in OMEGA_V and OMEGA_M are pretty close to the solid line, Dose this mean that it was converged?

Thanks for all help, limh

Post Reply