Вычисляем возраст Вселенной в R
Как узнали возраст Вселенной?
В 1929 году американский астроном Эдвин Хаббл обнаружил, что галактики удаляются друг от друга. Поделив расстояние между ними на скорость удаления, можно вычислить, как давно они были в одной точке. Это грубая оценка возраста Вселенной.
График зависимости скорости от расстояния в оригинальной работе Хаббла
“Как видите, он был хорошим ученым, потому что провел прямую линию через набор точек” (см. “Возникновение Вселенной”)
Вычисляем возраст Вселенной в R
Исходные данные
Попробуем почувствовать себя настоящими учеными, проведя прямую линию через набор точек. Для начала загрузим данные Freedman et al. (2001) по расстоянию (Distance, Mpc) до 36 сверхновых звезд типа Ia, а также по скорости удаления (Velocity, km/s) этих сверхновых, полученные при помощи космического телескопа “Хаббл”. Данные доступны в пакете lava:
library(lava)
## lava version 1.5.1
data(hubble) str(hubble)
## 'data.frame': 36 obs. of 3 variables: ## $ v : num 9065 12012 15055 16687 9801 ... ## $ D : num 135 159 199 239 117 ... ## $ sigma: num 2.3 3.1 2.8 2.8 3.4 2.9 3.1 2.2 3.4 2.7 ... ## - attr(*, "labels")= chr "Velocity (km/s)" "Distance (Mpc)" "sigma"
Строим диаграмму Хаббла
plot(hubble$D, hubble$v, axes = FALSE, pch = 20, col = "blue",
xlim= c(0,500), ylim = c(0, 30000),
main = "", xlab = "", ylab = "")
axis(1, at = seq(0, 500, by = 100), labels = seq(0, 500, by = 100))
axis(2, at = x <- seq(0, 3e4, by = 1e4),
labels = expression(0, 1 %*% 10^4, 2 %*% 10^4, 3 %*% 10^4))
title(main = "Диаграмма Хаббла \n для сверхновых типа Ia",
xlab = "Расстояние, Мпк",
ylab = "Скорость, км/с")
abline(lm(hubble$v~hubble$D - 1), col = "red", lwd = 2)
Вычисляем постоянную Хаббла при помощи простой линейной регрессии
fit <- lm(hubble$v~hubble$D - 1)
summary(fit)
## ## Call: ## lm(formula = hubble$v ~ hubble$D - 1) ## ## Residuals: ## Min 1Q Median 3Q Max ## -2748.6 -329.6 127.4 703.3 3020.8 ## ## Coefficients: ## Estimate Std. Error t value Pr(>|t|) ## hubble$D 70.6672 0.8385 84.28 <2e-16 *** ## --- ## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 ## ## Residual standard error: 1012 on 35 degrees of freedom ## Multiple R-squared: 0.9951, Adjusted R-squared: 0.995 ## F-statistic: 7103 on 1 and 35 DF, p-value: < 2.2e-16
Как видим, оцененное значение постоянной Хаббла составило 70.667 км/с на мегапарсек. Это означает, что галактики, находящиеся на расстоянии 100 мегапарсек, удаляются от нас со скоростью 7000—8000 км/с. Постоянная Хаббла имеет размерность, обратную времени и выражается обычно в км/с на мегапарсек. Один мегапарсек - это 3.09∗10^19 км.
mpc <- 3.09 * 10^19 mpc
## [1] 3.09e+19
ysec <- 60^2*24*365.25 ysec
## [1] 31557600
age <- 1/(fit$coefficients[[1]]/mpc) age
## [1] 4.372608e+17
age/ysec
## [1] 13855957106
age/ysec/10e8
## [1] 13.85596
Итак, мы определили, что возраст Вселенной составляет ~ 13.86 миллиардов лет.
Из биографии Эдвина Хаббла: Будущий астроном мало занимался математикой. Трудности с правописанием не оставляли его до конца жизни. В возрасте 25 лет Эдвин не мог финансово помогать семье, не знал, чем хочет заниматься, у него даже не было девушки. Только в возрасте 40 лет Хаббл написал свою первую статью, в которой сформулировал свой закон. Учился на юриста, но стал выдающимся астрономом, работающим у телескопа в военной форме. Высокомерный Хаббл враждовал с коллегами, но дружил с одним астрономом – представительницей Гарварда Сесилией Пейн, которую называл “лучшим мужчиной Гарварда”.
Вселенная Кеплера соответствовала Солнечной системе, Гершель расширил ее до размеров галактики, а во вселенной Хаббла наша галактика стала всего лишь маленькой точкой среди миллиардов других галактик. Он первым в мире осознал, насколько огромна Вселенная.
“Наблюдения всегда подразумевают под собой теорию” (с) Эдвин Хаббл