Как рассчитать возраст Вселенной с помощью астропии?

Я хочу преобразовать значение постоянной Хаббла H0 в Gyr, а также рассчитать возраст Вселенной t0 с помощью python. Код, который я использовал, показан ниже. Преобразование H0 выглядит нормально, но значение t0 должно быть 13,78, но оно показывает 0,013. Но если я произведу расчет t0 аналитически, то легко найду 13,7. Почему это происходит? Ожидая вашей помощи.

import numpy as np

Om=0.3
Ol=0.7
H0 = 72  #km/s/Mpc
from astropy.cosmology import WMAP9 as cosmo
from astropy import units
H0 = cosmo.H(0)
H0.value, H0.unit  
print H0.to('Gyr**-1')
H0 = H0.value

# inverse sinh (Ol/Om)^0.5 = 1.2099
t0 = (2./(3.*H0*np.sqrt(Ol))) * 1.2099
print t0

astropy
person Photon    schedule 05.12.2018    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Ниже приводится наивный возраст Вселенной (без учета ускорения):

from astropy.cosmology import WMAP9 as cosmo
import astropy.units as u
H0 = cosmo.H(0)
t0 = 1.0/H0
t0_Gyr = t0.decompose().to(u.Gyr)
print (t0_Gyr)

<Quantity 14.10548502 Gyr>
person astrosnapper    schedule 05.12.2018
comment
Большое спасибо за ответ. Хотя я сделал это успешно, задав этот вопрос здесь. - person Photon; 12.12.2018