Астрономический журнал, 2023, T. 100, № 2, стр. 133-143
Устранение хаббловского несоответствия при наличии взаимосвязи темной энергии и материи в современной Вселенной
Г. С. Бисноватый-Коган 1, *, А. М. Никишин 2, **
1 Институт космических исследований РАН
Москва, Россия
2 Московский инженерно-физический институт (МИФИ)
Москва, Россия
* E-mail: gkogan@iki.rssi.ru
** E-mail: nikishin-5@yandex.ru
Поступила в редакцию 31.08.2022
После доработки 14.11.2022
Принята к публикации 08.12.2022
- EDN: CJENUQ
- DOI: 10.31857/S0004629923020032
Аннотация
В современной космологии принято, что скалярное поле, ответственное за инфляционную стадию ранней Вселенной, полностью превращается в вещество. Предполагается, что ускоренное расширение в настоящее время управляется темной энергией (DE), которая, по всей вероятности, определяется космологической постоянной Эйнштейна, не связанной со скалярным полем, ответственным за инфляцию. Мы рассматриваем здесь космологическую модель, в которой DE может иметь в настоящее время два компонента, один из которых – постоянная Эйнштейна ($\Lambda $), а другой, меньший переменный компонент DEV (${{\Lambda }_{V}}$), связан с остатком скалярного поля, вызвавшего инфляцию, после того, как основная часть скалярного поля превратилась в вещество. Мы рассматриваем здесь только стадии эволюции Вселенной после рекомбинации ($z \lesssim 1100$), где DM – преобладающий компонент вещества. Предполагается, что превращение скалярного поля в вещество продолжается в настоящее время и сопровождается обратным процессом превращения DM в скалярное поле. Рассматривается связь между DM и DEV, которая приводит к линейному соотношению между плотностями энергии этих компонентов после рекомбинации ${{\rho }_{{{\text{DM}}}}} = \alpha {\kern 1pt} {{\rho }_{{{\text{DEV}}}}}$. Рассматриваются также варианты с зависимостью от красного смещения $z$ коэффициента $\alpha (z)$. Одна из возникших в современной космологии проблем, получившая название Hubble Tension (HT) – Несоответствие Хаббла, состоит в расхождении значений постоянной Хаббла в настоящее время (${{H}_{0}}$), измеренных по наблюдениям Вселенной на малых красных смещениях ($z \lesssim 1$) и по наблюдениям флуктуаций реликтового излучения во Вселенной при больших красных смещениях ($z \approx 1100$). В рассматриваемой модели это несоответствие может быть объяснено отклонением существующей космологической модели от использованной общепринятой $\Lambda $CDM модели плоской Вселенной действием добавочного компонента темной энергии DEV на стадиях после рекомбинации. В рамках этой расширенной модели мы рассматриваем различные функции $\alpha {\kern 1pt} (z)$, которые могут устранить HT. Чтобы поддерживать близким к постоянному соотношение плотностей энергии DEV и DM на протяжении интервала $0 \leqslant z \lesssim 1100$, необходимо допустить существование широкого спектра масс частиц темной материи.
1. ВВЕДЕНИЕ
Закон Хаббла ${v} = Hr$ является одним из важнейших законов космологии, задающим скорость разбегания космических объектов ${v}$ в зависимости от расстояния между ними $r$. Измерение параметра ${{H}_{0}}$ в настоящую эпоху является весьма трудной и нетривиальной задачей, которой многие годы занимаются различные научные группы. Скорость разбегания определяется по красному смещению линий, наблюдаемых в спектрах улетающих галактик, но наибольшая трудность состоит в измерениях расстояний до них. При этом используется метод “лестницы расстояний” (cosmic distance ladder), в котором используются все более яркие стандартные свечи – объекты с известной светимостью. Этот метод позволяет измерить расстояния до относительно близких объектов (в пределах тысячи парсек), и не свободен от систематических ошибок, связанных с неизбежным разбросом светимостей используемых стандартных свечей. Использование разных ступеней в этой лестнице разными группами привело в 1972–1974 гг. к существенному расхождению значений ${{H}_{0}}$: от $ \sim {\kern 1pt} 50$ км/с/Mпс группой Сэндеджа-Таммана до $ \sim {\kern 1pt} 100$ км/с/Mпс группой Вокулера [1]. Со временем измерения ${{H}_{0}}$ на красных смещениях $ \lesssim {\kern 1pt} 1$ были существенно улучшены благодаря созданию больших телескопов, включая измерения на телескопах миссии Хаббла, что позволило ограничить интервал значений до ${{H}_{0}} \approx 72{\kern 1pt} - {\kern 1pt} 75$ км/с/Mпс.
Измерения спектра флуктуаций реликтового излучения на спутниках WMAP и Planck привели к возможности независимого измерения параметра Хаббла ${{H}_{r}}$ в эпоху рекомбинации. Используя теоретическую зависимость $H(z)$ в рамках принятой космологической модели, предположительно $\Lambda $CDM, было получено современное значение параметра Хаббла $H_{0}^{{{\text{Dist}}}}$, которое отличалось от значения, полученного в локальных измерениях $H_{0}^{{{\text{Loc}}}}$ на статистически значимую величину в пределах $(4.5{\kern 1pt} - {\kern 1pt} 6.3)\sigma $ [2]. В этом расхождении и заключается проблема HT (см., однако, [3]).
В космологических моделях, принятых в настоящее время, предполагается, что скалярное поле, ответственное за существование ранних стадий быстрого экспоненциального расширения, полностью превращается в вещество в процессе инфляции. Однако современное ускоренное разбегание галактик определяется темной энергией (DE), предположительно постоянной Эйнштейна $\Lambda $, которая никак не связана со скалярным полем, которое привело к инфляции.
Согласно измерениям флуктуаций реликта (WMAP, Planck), наша Вселенная с точностью $ < {\kern 1pt} 1\% $ является плоской, и ее средняя плотность равна критической ${{\rho }_{c}}$. Современная $\Lambda $CDM модель Вселенной характеризуется следующими параметрами:
(1)
$\begin{gathered} {\text{вклад}}\;{\text{нерелятивистского}}\;{\text{вещества}} \\ {{\Omega }_{{\text{M}}}} = {{\Omega }_{{{\text{DM}}}}} + {{\Omega }_{{\text{B}}}} \approx 0.26 + 0.04 \approx 0.3, \\ {\text{вклад}}\;{\text{темной}}\;{\text{энергии}}\quad {{\Omega }_{\Lambda }} \approx 0.7. \\ \end{gathered} $В данной работе, следуя [4, 5], для объяснения несоответствия HT мы рассматриваем космологическую модель, которая является расширением принятой $\Lambda $CDM модели, и включает два компонента темной энергии, один из которых – постоянная Эйнштейна ($\Lambda $), а другой, малый переменный компонент DEV (${{\Lambda }_{V}}$), является остатком скалярного поля, создавшего инфляцию. Появление проблемы HT здесь связывается с пересчетом параметра Хаббла ${{H}^{{{\text{Dist}}}}}({{z}_{r}})$, измеренного на момент рекомбинации, к современному значению ${{H}^{{{\text{Dist}}}}}(0)$, используя стандартную $\Lambda $CDM модель, которая предполагается неполной, что приводит после пересчета к заниженному значению современной постоянной Хаббла. В данной работе получено более корректное решение космологического уравнения, где учитывается действие переменного добавочного компонента темной энергии DEV, что существенно уточняет результаты, полученные в [4, 5], для значений плотности DEV, при которых проблема HT не возникает.
2. НЕСООТВЕТСТВИЕ ХАББЛА
Проблема, возникшая в космологии в последние годы, состоит в различии значений постоянной Хаббла в настоящую эпоху, получаемых в разных экспериментах. Анализ наблюдений миссии Planck, измеряющей флуктуации реликтового излучения в период рекомбинации, приводит к современному значению постоянной Хаббла [6–8]:
В то же время измерения с использованием сверхновых типа Ia (SNIa) с калибровкой расстояния по цефеидам [9–13], дают значение Измерения с использованием временны́х задержек линзированных квазаров [14] дают значение ${{H}_{0}} = 73.3_{{ - 1.8}}^{{ + 1.7}}$ км/с/Мпк. В работе [15] было найдено значение ${{H}_{0}} = 72.4 \pm 1.9$ км/с/Мпк, используя ветвь красных гигантов, приложенную к SNIa, которая не зависит от шкалы расстояний цефеид. Анализ набора этих и других недавних измерений на малых и больших красных смещениях показывает [16], что несоответствие между результатами Planck [8], и любыми тремя независимыми измерениями в поздней Вселенной лежит в интервале между $4\sigma $ и $6\sigma $. Несколько новых космических экспериментов были предложены для проверки достоверности этого расхождения значений постоянной Хаббла [17, 18].Было предложено много различных вариантов объяснения возникновения HT, некоторые из которых были опровергнуты наблюдательными данными. Общепринятого и доказанного экспериментально решения проблемы HT в настоящее время пока нет. Предложенные решения можно разбить на дорекомбинационные и послерекомбинационные, которые неким образом меняют процесс эволюции Вселенной либо в период до рекомбинации, либо после нее соответственно. Рассматриваются также варианты с использованием теорий гравитации, основанных на модификациях ОТО (см. работы [19–41]). Подробный обзор большинства предложенных методов решения данной проблемы приведен в работе [18]. В ней изложены предлагаемые способы объяснения данного феномена, а также возможное влияние их присутствия на другие космологические параметры.
Наблюдается значимое расхождение между экспериментальными значениями ${{H}_{0}}$, полученными по флуктуациям реликтового излучения миссией Planck, и значениями, полученными в локальных измерениях. В связи с этим мы рассмотрим возможность решения проблемы HT, как расхождения между результатами этих экспериментов, рассматривая усредненную величину, полученную в локальных измерениях, как истинную. В дальнейших численных расчетах мы используем значение $H_{0}^{{{\text{Dist}}}}$, полученное в результате измерений на больших красных смещениях, и значение $H_{0}^{{{\text{Loc}}}}$, полученное в результате измерений в локальной Вселенной.
3. ПРЕДЛАГАЕМОЕ РЕШЕНИЕ ПРОБЛЕМЫ HT
Темная материя и темная энергия составляют около 96% всей энергии во Вселенной [6, 9, 10], но их природа до сих пор неизвестна. Современное значение плотности DE может быть представлено космологической постоянной Эйнштейна $\Lambda $ [42], но также может быть результатом действия Хиггского скалярного поля, которое предлагается в качестве причины инфляции в ранней Вселенной [43–46]. Величина индуцированной ${{\Lambda }_{V}}$, относящейся к инфляции, на много порядков больше, чем ее современное значение и не было попыток найти связь между ними. Физическая природа DM остается неясной. Существует много предположений о ее происхождении [34, 47, 48], но ни одно из них не подтверждено экспериментально или наблюдательно, в то время как многие из них были опровергнуты.
Чтобы объяснить происхождение Hubble Tension в [4, 5], была введена переменная часть темной энергии (DEV) как переменная составляющая космологической постоянной ${{\Lambda }_{V}}$. Ее вклад в процесс расширения Вселенной растет с увеличением $z$, и превышает вклад космологической постоянной Эйнштейна $\Lambda $ при красных смещениях $z \gg 1$. Таким образом, DE может иметь два компонента, один из которых является постоянной Эйнштейна $\Lambda $, и другой, малый в настоящее время переменный компонент DEV – ${{\Lambda }_{V}}$, который может иметь происхождение из остатков инфляционного скалярного поля, являющегося источником происхождения существующего вещества. Предполагалось, что процесс перехода энергии поля в вещество в настоящее время сопровождается обратным процессом превращения массы в энергию поля, что приводит к динамической связи между плотностями материи и ${{\Lambda }_{V}}$.
Мы рассматриваем только послерекомбинационный период ($z \lesssim 1100$) расширения Вселенной, где DM представляет собой самый значительный компонент вещества во Вселенной. Поэтому для простоты считаем, что между плотностями энергии темной материи и добавочной части темной энергии существует связь вида
Здесь $\alpha (z)$ – функция красного смещения, форма которой ограничивается требованием того, чтобы она позволила устранить противоречие HT, не создавая при этом дополнительных трудностей при интерпретации космологических наблюдений. Как показано в [4, 5], этому требованию удовлетворяет, например, постоянное значение $\alpha $, устраняющее HT. Для поддержания динамического равновесия (4) темная материя должна состоять из спектра частиц различной массы, включая очень легкие частицы. Рождение таких частиц в процессе взаимных превращений DM и скалярного поля на всем интервале $z \lesssim 1100$ позволит поддерживать динамическое равновесие между плотностями DM и DEV типа (4) при уменьшении энергии всех переменных составляющих Вселенной при ее расширении.Как следует из численного моделирования процессов, ведущих к формированию современной крупномасштабной структуры Вселенной, наилучшее согласие с наблюдениями показывает модель с холодной DM, т.е. частицы темной материи являются нерелятивистскими. Для этого часто рассматриваются массивные частицы из суперсимметричной теории поля (нейтралино, фотино и др. [49]). C другой стороны, в качестве DM рассматриваются очень легкие частицы аксионы [47], существование которых следует из некоторых теоретических моделей, объясняющих наблюдаемое нарушение CP инвариантности в ядерных процессах. В нашей модели Вселенной для устранения HT требуется более сложная структура DM, где могут присутствовать частицы промежуточных масс, причем масса легких частиц (аксионов) также может быть представлена целым спектром.
В присутствии DEV постоянная Хаббла уменьшается со временем медленнее, чем без нее. Это создает большее современное значение ${{H}_{0}}$ для одного и того же значения ${{H}_{{{\text{rec}}}}}$ в эпоху рекомбинации. Поэтому мы предполагаем, что значение $H_{0}^{{{\text{Dist}}}}$, измеренное телескопом Planck, было получено с помощью экстраполирования значения ${{H}_{{{\text{rec}}}}}$ с момента рекомбинации ${{z}_{r}} \approx 1100$ до настоящего времени $z = 0$ в модели Фридмана плоской пылевой Вселенной с учетом космологической постоянной $\Lambda $. В случае равнораспределенной Вселенной экстраполяция эта должна быть произведена в модели, учитывающей добавочные компоненты темной энергии DEV. В нашей интерпретации HT связано с неточной экстраполяцией данных Planck в рамках модели, не учитывающей DEV.
Примем, что истинное современное значение постоянной Хаббла определяется локальными измерениями, т.е. $H_{0}^{{{\text{Loc}}}}$, и оба измерения корректны, но значение $H_{0}^{{{\text{Dist}}}}$ возникло при неточной экстраполяции. Зная это, можно найти такую функцию $\alpha (z)$, чтобы при пересчете данных Planck получалось значение постоянной Хаббла, совпадающее с локальными измерениями. Таким образом, задача сводится к нахождению такой функции $\alpha (z)$, при которой в процессе пересчета постоянной Хаббла в настоящее время HT не возникает. Чтобы вычислить влияние малой добавки DEV на расчет современного значения ${{H}_{0}}$ из измерений ${{H}_{{{\text{rec}}}}}$ на момент рекомбинации, мы должны построить космологическую модель с учетом компонента DEV на интервале $z \in [0,{{z}_{r}}]$.
4. РАСШИРЕНИЕ ВСЕЛЕННОЙ
Эволюция Вселенной описывается уравнениями ОТО Эйнштейна [1]:
(5)
${{R}_{{ik}}} - \frac{1}{2}{{g}_{{ik}}}{\kern 1pt} R - \Lambda {{g}_{{ik}}} = \frac{{8\pi G}}{{{{c}^{4}}}}{{T}_{{ik}}}.$(6)
${{\left( {\frac{{\dot {a}}}{a}} \right)}^{2}} = \frac{{8\pi G}}{3}\rho + \frac{{\Lambda {{c}^{2}}}}{3} - \frac{{\kappa {{c}^{2}}}}{{{{a}^{2}}}}{\kern 1pt} ,$5. СКАЛЯРНОЕ ПОЛЕ – ИСТОЧНИК DM И DEV
Скалярное поле принимается в качестве основной причины инфляционной стадии и рождения вещества во Вселенной [43–45]. Рассмотрим скалярное поле интенсивности $\phi $, находящееся в потенциале $V(\phi )$. В однородной изотропной расширяющейся Вселенной зависимость $\phi $ от времени определяется уравнением [51]:
Плотность энергии ${{\rho }_{V}}$ и давление ${{P}_{V}}$ скалярного поля определяются как в работе [51] (здесь и в большинстве дальнейших уравнений принимается $c = 1$)(8)
${{\rho }_{V}} = \frac{{{{{\dot {\phi }}}^{2}}}}{2} + V(\phi ),\quad {{P}_{V}} = \frac{{{{{\dot {\phi }}}^{2}}}}{2} - V(\phi ).$Рассмотрим Вселенную с начальной интенсивностью скалярного поля ${{\phi }_{{{\text{in}}}}}$, начальным потенциалом ${{V}_{{{\text{in}}}}}$ и нулевой начальной производной ${{\dot {\phi }}_{{{\text{in}}}}} = 0$. Производная интенсивности скалярного поля растет на начальных стадиях инфляции при уменьшении потенциала $V$. Предположим, что после достижения справедливости соотношения
оно продолжает выполняться на дальнейших стадиях расширения. Кинетическая часть энергии скалярного поля превращается в вещество, предположительно, главным образом в темную материю, а функция $\alpha (z)$ определяет связь плотности темной энергии (DE), определяемой величиной $V$, и плотностью вещества, определяемой кинетическим членом. Как следует из наблюдений, основная часть DE в настоящее время может быть связана с постоянной Эйнштейна $\Lambda $. На предшествующих стадиях расширения постоянная $\Lambda $ была меньше переменной части ${{\Lambda }_{V}}$ в широком интервале допустимых функций $\alpha (z)$. При выполнении условия ((9)) введем следующие обозначения(10)
${{\rho }_{\phi }} = V,\quad {{P}_{\phi }} = - V,\quad {{\rho }_{m}} = \frac{{{{{\dot {\phi }}}^{2}}}}{2}{\kern 1pt} ,\quad {{P}_{m}} = \beta {\kern 1pt} \frac{{{{{\dot {\phi }}}^{2}}}}{2}{\kern 1pt} ,$(11)
$\begin{gathered} \rho = {{\rho }_{\phi }} + {{\rho }_{m}} = (1 + \alpha ){\kern 1pt} V, \\ P = {{P}_{\phi }} + {{P}_{m}} = - (1 - \alpha \beta ){\kern 1pt} V. \\ \end{gathered} $(12)
$\dot {\rho } = - 3\alpha {\kern 1pt} \frac{{1 + \beta }}{{1 + \alpha }}\frac{{\dot {a}}}{a}\rho ,\quad \frac{\rho }{{{{\rho }_{ \star }}}} = {{\left( {\frac{{{{a}_{ \star }}}}{a}} \right)}^{{3\alpha (1 + \beta )/(1 + \alpha )}}}.$6. КОСМОЛОГИЧЕСКИЕ МОДЕЛИ И РЕШЕНИЕ ПРОБЛЕМЫ HT
Решение космологических уравнений (6) и (12) для плоской пылевой Вселенной после рекомбинации при $\kappa = \beta = 0$, и постоянном $\alpha $ рассматривалось в [4, 5]. Приведем решения для двух случаев.
6.1. $\Lambda = 0$, $\alpha (z) = {\text{const}}$
В этом случае получаем следующие зависимости функций от времени:
(13)
$\begin{gathered} \frac{a}{{{{a}_{*}}}} = {{(6\pi G\rho _{*}^{{}}{{t}^{2}})}^{{(1 + \alpha )/3\alpha }}}{{\left[ {\frac{\alpha }{{1 + \alpha }}} \right]}^{{2(1 + \alpha )/3\alpha }}} = \\ \, = {{\left( {\frac{{\rho _{*}^{{}}}}{\rho }} \right)}^{{(1 + \alpha )/3\alpha }}} = {{\left( {\frac{t}{{t_{*}^{{}}}}} \right)}^{{2(1 + \alpha )/3\alpha }}}, \\ \end{gathered} $(14)
$\rho = {{\left( {\frac{{1 + \alpha }}{\alpha }} \right)}^{2}}\frac{1}{{6\pi G{{t}^{2}}}}{\kern 1pt} .$(15)
${{H}_{\alpha }} \equiv \frac{{\dot {a}}}{a} = \frac{{2(1 + \alpha )}}{{3\alpha {\kern 1pt} t}} = \frac{2}{3}\sqrt {6\pi G{{\rho }_{0}}} {{(z + 1)}^{{3\alpha /2(1 + \alpha )}}}.$(16)
${{H}_{\alpha }}(z) = \frac{2}{3}\sqrt {6\pi G{{\rho }_{0}}} {\kern 1pt} {{({{z}_{r}} + 1)}^{{3/2(1 + \alpha )}}}{{(z + 1)}^{{3\alpha /2(1 + \alpha )}}},$(18)
$\frac{{H_{0}^{{{\text{Loc}}}}}}{{H_{0}^{{{\text{Dist}}}}}} = ({{z}_{r}} + {{1)}^{{3/2(1 + \alpha )}}},\quad {{\alpha }_{{{\text{HT}}}}} \approx 133.$Таблица 1.
Параметр | Значение |
---|---|
Локальное значение постоянной Хаббла | $H_{0}^{{{\text{Loc}}}} \approx 73$ км/с/Мпк |
Постоянная Хаббла, измеренная по реликтовому излучению | $H_{0}^{{{\text{Dist}}}} \approx 67.5$ км/с/Мпк |
Полная плотность вещества в плоской Вселенной | ${{\rho }_{{{\text{tot}}}}} \approx 1.066 \times {{10}^{{ - 29}}}$ г/см3 |
Измеренное локально значение плотности космологической постоянной [10] (статистика 2σ) | ${{\rho }_{\Lambda }} = (0.44{\kern 1pt} - {\kern 1pt} 0.96){\kern 1pt} {{\rho }_{{{\text{tot}}}}}$ |
Измеренное дистанционно значение плотности космологической постоянной | ${{\rho }_{\Lambda }} \approx 0.7{{\rho }_{{{\text{tot}}}}}$ |
Космологическая постоянная из дистанционных измерений | $\Lambda = \frac{{8\pi {\kern 1pt} G{{\rho }_{\Lambda }}}}{{{{c}^{2}}}} \approx 1.40 \times {{10}^{{ - 56}}}$ см‒2 |
Асимптотика постоянной Хаббла | ${{H}_{{{\text{ac}}}}} = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3}} \approx 63.2$ км/с/Мпк |
Усредненный возраст Вселенной [52] | ${{t}_{0}} \approx 4.35 \times {{10}^{{17}}}$ c |
Красное смещение, соответствующее эпохе рекомбинации | ${{z}_{r}} \approx 1100$ |
6.2. $\Lambda \ne 0$, $\alpha (z) = {\text{const}}$
Современные наблюдения (1) указывают на преобладание темной энергии, определяемой космологической постоянной $\Lambda $. С учетом дополнительного переменного члена DEV, параметры пылевой Вселенной после рекомбинации ($z < 1100$) описываются формулами, получающимися из решения системы (6), (12) в виде [4, 5]:
(19)
$\begin{gathered} {{\left( {\frac{a}{{a_{*}^{{}}}}} \right)}^{{3\alpha /2(1 + \alpha )}}} = \sqrt {\frac{{8\pi G\rho _{*}^{{}}}}{{\Lambda {{c}^{2}}}}} \sinh \left( {\sqrt {\frac{\Lambda }{3}} \frac{{3\alpha }}{{2(1 + \alpha )}}ct} \right) = \\ = \sqrt {\frac{{\rho _{*}^{{}}}}{\rho }} , \\ \end{gathered} $(20)
$\sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G\rho }}} = \sinh \left( {\sqrt {\frac{\Lambda }{3}} \frac{{3\alpha }}{{2(1 + \alpha )}}ct} \right).$(21)
$\begin{gathered} z + 1 \equiv \frac{{{{a}_{0}}}}{a} = \\ = {{\left[ {\sqrt {\frac{{8\pi G{{\rho }_{0}}}}{{\Lambda {{c}^{2}}}}} \sinh \left( {\sqrt {\frac{\Lambda }{3}} \frac{{3\alpha }}{{2(1 + \alpha )}}ct} \right)} \right]}^{{ - 2(1 + \alpha )/3\alpha }}} = \\ \, = {{\left( {\frac{\rho }{{{{\rho }_{0}}}}} \right)}^{{(1 + \alpha )/3\alpha }}}. \\ \end{gathered} $(22)
${{H}_{\alpha }}(t) \equiv \frac{{\dot {a}}}{a} = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3}} {\kern 1pt} \coth \left( {\sqrt {\frac{\Lambda }{3}} \frac{{3\alpha }}{{2(1 + \alpha )}}ct} \right).$(23)
$\begin{gathered} {{\left( {\frac{a}{{{{a}_{0}}}}} \right)}^{{3/2}}} = \sqrt {\frac{{8\pi G{{\rho }_{0}}}}{{\Lambda {{c}^{2}}}}} \sinh \left( {\sqrt {\frac{\Lambda }{3}} \frac{3}{2}ct} \right) = \sqrt {\frac{{{{\rho }_{0}}}}{\rho }} , \\ \sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G\rho }}} = \sinh \left( {\sqrt {\frac{\Lambda }{3}} \frac{3}{2}ct} \right), \\ \end{gathered} $(24)
$\begin{gathered} z + 1 \equiv \frac{{{{a}_{0}}}}{a} = {{\left[ {\sqrt {\frac{{8\pi G{{\rho }_{0}}}}{{\Lambda {{c}^{2}}}}} \sinh \left( {\sqrt {\frac{\Lambda }{3}} \frac{3}{2}ct} \right)} \right]}^{{ - 2/3}}} = {{\left( {\frac{\rho }{{{{\rho }_{0}}}}} \right)}^{{1/3}}}, \\ H(t) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3}} \coth \left( {\sqrt {\frac{\Lambda }{3}} \frac{3}{2}ct} \right). \\ \end{gathered} $(25)
$\begin{gathered} {{H}_{\alpha }}(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3}} \times \\ \times \;\coth \left( {{{{\sinh }}^{{ - 1}}}\left[ {\sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G{{\rho }_{0}}}}} {\kern 1pt} {{{(z + 1)}}^{{ - 3\alpha /2(1 + \alpha )}}}} \right]} \right), \\ \end{gathered} $(26)
$H(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3}} {\kern 1pt} \coth \left( {{{{\sinh }}^{{ - 1}}}\left[ {\sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G{{\rho }_{0}}}}} {\kern 1pt} {{{(z + 1)}}^{{ - 3/2}}}} \right]} \right).$(27)
$\begin{gathered} {{H}_{{r\alpha }}} \approx \frac{2}{3}\sqrt {6\pi G{{\rho }_{0}}} {\kern 1pt} {{({{z}_{r}} + 1)}^{{3\alpha /2(1 + \alpha )}}}, \\ {{H}_{r}} \approx \frac{2}{3}\sqrt {6\pi G{{\rho }_{0}}} {\kern 1pt} {{({{z}_{r}} + 1)}^{{3/2}}}. \\ \end{gathered} $(28)
$\begin{gathered} {{H}_{\alpha }}(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3}} \coth \left( {{{{\sinh }}^{{ - 1}}}\left[ {\sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G{{\rho }_{0}}}}} \times } \right.} \right. \\ \left. {\left. {\mathop \times \limits_{\mathop {}\limits_{_{{_{{_{{}}}}}}} }^{} \;{{{(z + 1)}}^{{ - 3\alpha /2(1 + \alpha )}}}{{{({{z}_{r}} + 1)}}^{{ - 3/2(1 + \alpha )}}}} \right]} \right), \\ \end{gathered} $(29)
$H(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3}} \coth \left( {{{{\sinh }}^{{ - 1}}}\left[ {\sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G{{\rho }_{0}}}}} {\kern 1pt} {{{(z + 1)}}^{{ - 3/2}}}} \right]} \right).$(30)
$\begin{gathered} {{(z + 1)}^{{ - 3\alpha /2(1 + \alpha )}}}{{({{z}_{r}} + 1)}^{{ - 3/2(1 + \alpha )}}}\sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G{{\rho }_{0}}}}} = \\ \, = \sinh \left[ {{{{\coth }}^{{ - 1}}}\left( {\frac{{{{H}_{\alpha }}}}{{\sqrt {\Lambda {{c}^{2}}/3} }}} \right)} \right], \\ \end{gathered} $(31)
${{(z + 1)}^{{ - 3/2}}}\sqrt {\frac{{\Lambda {{c}^{2}}}}{{8\pi G{{\rho }_{0}}}}} = \sinh \left[ {{{{\coth }}^{{ - 1}}}\left( {\frac{H}{{\sqrt {\Lambda {{c}^{2}}{\text{/}}3} }}} \right)} \right],$(32)
$\begin{gathered} \frac{{\sinh \left[ {{{{\coth }}^{{ - 1}}}\left( {\frac{{H_{0}^{{{\text{Loc}}}}}}{{\sqrt {\Lambda {{c}^{2}}{\text{/}}3} }}} \right)} \right]}}{{\sinh \left[ {{{{\coth }}^{{ - 1}}}\left( {\frac{{H_{0}^{{{\text{Dist}}}}}}{{\sqrt {\Lambda {{c}^{2}}{\text{/}}3} }}} \right)} \right]}} = ({{z}_{r}} + {{1)}^{{ - 3/2(1 + \alpha )}}}, \\ ^{{^{{}}}}{{\alpha }_{{{\text{HT}}}}} \approx 24. \\ \end{gathered} $7. РЕШЕНИЕ ЗАДАЧИ С ИСПОЛЬЗОВАНИЕМ ФУНКЦИИ $H(z)$
Если ограничить поставленную задачу, исключив из нее вопрос построения космологической модели в виде зависимости $a(t)$ и оставив только нахождение поведения параметра Хаббла в виде зависимости $H(z)$, то задача о решении проблемы HT существенно упрощается. В частности удается аналитически решить ее для целого набора функций $\alpha (z)$. Космологическое уравнение (6) и уравнение адиабаты (12), при $\kappa = \beta = 0$, запишутся в виде:
(33)
${{H}^{2}} = \frac{{8\pi G}}{3}\rho + \frac{{\Lambda {{c}^{2}}}}{3},\quad \frac{{\dot {\rho }}}{\rho } = - \frac{{3\alpha }}{{1 + \alpha }}\frac{{\dot {a}}}{a}.$7.1. $\alpha = {\text{const}}$, $\Lambda \ne 0$
Аналогично предыдущему методу, сначала рассмотрим случай постоянной функции $\alpha (z)$. Второе уравнение в (33) для этого случая проинтегрировано в (12). Запишем это решение в виде
(34)
$\frac{{{{\rho }_{\alpha }}}}{{{{\rho }_{{0\alpha }}}}} = {{\left( {\frac{{{{a}_{0}}}}{{{{a}_{{0\alpha }}}}}} \right)}^{{3\alpha /(1 + \alpha )}}} = (1 + z{{)}^{{3\alpha /(1 + \alpha )}}}.$(35)
${{H}_{\alpha }}(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3} + \frac{{8\pi G}}{3}{{\rho }_{{0\alpha }}}{{{(1 + z)}}^{{3\alpha /(1 + \alpha )}}}} {\kern 1pt} .$(36)
$H(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3} + \frac{{8\pi G}}{3}{{\rho }_{0}}{{{(1 + z)}}^{3}}} .$(37)
$\begin{gathered} {{\rho }_{{0\alpha }}}{{(1 + {{z}_{r}})}^{{3\alpha /(1 + \alpha )}}} = {{\rho }_{0}}{{(1 + {{z}_{r}})}^{3}}, \\ {{\rho }_{{0\alpha }}} = {{\rho }_{0}}{{(1 + {{z}_{r}})}^{{3/(1 + \alpha )}}}{\kern 1pt} . \\ \end{gathered} $Отождествляя вычисленные современные значения $H$ с наблюдаемыми
получаем, используя (35)–(37), уравнение для определения значения $\alpha $, устраняющего HT, в виде(39)
$\frac{{{{{(H_{0}^{{{\text{Loc}}}})}}^{2}} - \frac{{\Lambda {{c}^{2}}}}{3}}}{{{{{(H_{0}^{{{\text{Dist}}}})}}^{2}} - \frac{{\Lambda {{c}^{2}}}}{3}}} = \frac{{{{\rho }_{{0\alpha }}}}}{{{{\rho }_{0}}}} = (1 + {{z}_{r}}{{)}^{{3/(1 + \alpha )}}}{\kern 1pt} .$(40)
$\mathop {\sinh }\nolimits^2 \left[ {\mathop {\coth }\nolimits^{ - 1} \left( {\frac{{H_{0}^{{{\text{Loc}}}}}}{{\sqrt {\Lambda {{c}^{2}}{\text{/}}3} }}} \right)} \right] = {{\left[ {{{{\left( {\frac{{H_{0}^{{{\text{Loc}}}}}}{{\sqrt {\Lambda {{c}^{2}}{\text{/}}3} }}} \right)}}^{2}} - 1} \right]}^{{ - 1}}},$(41)
$\mathop {\sinh }\nolimits^2 \left[ {\mathop {\coth }\nolimits^{ - 1} \left( {\frac{{H_{0}^{{{\text{Dist}}}}}}{{\sqrt {\Lambda {{c}^{2}}{\text{/}}3} }}} \right)} \right] = {{\left[ {{{{\left( {\frac{{H_{0}^{{{\text{Dist}}}}}}{{\sqrt {\Lambda {{c}^{2}}{\text{/}}3} }}} \right)}}^{2}} - 1} \right]}^{{ - 1}}}.$(42)
$\begin{gathered} \sinh ({{\coth }^{{ - 1}}}x) = \frac{1}{{\sqrt {{{x}^{2}} - 1} }}, \\ \mathop {\coth }\nolimits^{ - 1} x = {{\sinh }^{{ - 1}}}\frac{1}{{\sqrt {{{x}^{2}} - 1} }}. \\ \end{gathered} $(43)
$\begin{gathered} {{\coth }^{{ - 1}}}x = \frac{1}{2}\left[ {\ln \left( {1 + \frac{1}{x}} \right) - \ln \left( {1 - \frac{1}{x}} \right)} \right], \\ {{\sinh }^{{ - 1}}}y = \ln (y + \sqrt {{{y}^{2}} + 1} ). \\ \end{gathered} $7.2. $\alpha (z) = \frac{\alpha }{{{{{(1 + z)}}^{\gamma }}}},\;\Lambda \ne 0$
Рассмотрим модель с зависящей от красного смещения связью плотностей ${{\rho }_{m}}$ и ${{\rho }_{{{\text{DEV}}}}}$. Считаем параметр $\gamma \ll 1$, чтобы значения функции $\alpha (z)$ не слишком различались в эпоху рекомбинации и в настоящее время. Тогда вместо второго уравнения (33) получаем с учетом (21)
(44)
$\frac{{\dot {\rho }}}{\rho } = - \frac{{3\alpha (z)}}{{1 + \alpha (z)}}\frac{{\dot {a}}}{a} = \frac{{3\alpha }}{{\alpha + {{{(1 + z)}}^{\gamma }}}}\frac{{\dot {z}}}{{1 + z}}.$(45)
$\frac{{{{\rho }_{\alpha }}}}{{{{\rho }_{{0\alpha }}}}} = {{\left[ {\frac{{1 + \alpha }}{{1 + \alpha {{{(1 + z)}}^{{ - \gamma }}}}}} \right]}^{{3/\gamma }}}.$(47)
$3\alpha \int \frac{{dz}}{{(1 + z)[\alpha + {{{(1 + z)}}^{\gamma }}]}} = - \frac{3}{\gamma }\ln (1 + \alpha {{(1 + z)}^{{ - \gamma }}}).$(48)
$\begin{gathered} {{H}_{\alpha }}(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3} + \frac{{8\pi G}}{3}{{\rho }_{{0\alpha }}}{{{\left[ {\frac{{1 + \alpha }}{{1 + \alpha {{{(1 + z)}}^{{ - \gamma }}}}}} \right]}}^{{3/\gamma }}}} , \\ H(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3} + \frac{{8\pi G}}{3}{{\rho }_{0}}{{{(1 + z)}}^{3}}} . \\ \end{gathered} $(49)
${{\rho }_{{0\alpha }}}{{\left[ {\frac{{1 + \alpha }}{{1 + \alpha {{{(1 + {{z}_{r}})}}^{{ - \gamma }}}}}} \right]}^{{3/\gamma }}} = {{\rho }_{0}}{{(1 + {{z}_{r}})}^{3}}.$(50)
$\begin{gathered} {{H}_{\alpha }}(z) = \\ = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3} + \frac{{8\pi G}}{3}{{\rho }_{0}}{{{\left[ {\frac{{1 + \alpha {{{(1 + {{z}_{r}})}}^{{ - \gamma }}}}}{{1 + \alpha {{{(1 + z)}}^{{ - \gamma }}}}}} \right]}}^{{3/\gamma }}}{{{(1 + {{z}_{r}})}}^{3}}} , \\ H(z) = \sqrt {\frac{{\Lambda {{c}^{2}}}}{3} + \frac{{8\pi G}}{3}{{\rho }_{0}}{{{(1 + z)}}^{3}}} {\kern 1pt} . \\ \end{gathered} $(51)
$\frac{{{{{(H_{0}^{{{\text{Loc}}}})}}^{2}} - \frac{{\Lambda {{c}^{2}}}}{3}}}{{{{{(H_{0}^{{{\text{Dist}}}})}}^{2}} - \frac{{\Lambda {{c}^{2}}}}{3}}}{\kern 1pt} {{\left[ {\frac{{1 + \alpha }}{{1 + \alpha {{{(1 + {{z}_{r}})}}^{{ - \gamma }}}}}} \right]}^{{3/\gamma }}} = (1 + {{z}_{r}}{{)}^{3}}.$Возраст Вселенной ${{t}_{0}}$, где отсутствует HT, зависит от выбора пары ($\gamma ,{\kern 1pt} {\kern 1pt} {{\alpha }_{{{\text{HT}}}}}$), где только параметр $\gamma $ является независимым. Только при одном значении $\gamma $ этот возраст совпадает с возрастом из табл. 1, полученным на основе наблюдательных данных [52]. Для расчета возраста в модели используются следующие соотношения, рассмотренные выше
(52)
${{t}_{0}} = \int\limits_0^{{{a}_{0}}} \frac{{da}}{{a{\kern 1pt} H}} = \int\limits_0^1 \frac{{dx}}{{x{\kern 1pt} {{H}_{\alpha }}(x)}}{\kern 1pt} ,$8. ЗАКЛЮЧЕНИЕ
Мы рассмотрели космологическую модель, в которой вклад темной энергии DE состоит из двух компонентов: космологической постоянной Эйнштейна $\Lambda $ и малой переменной составляющей ${{\Lambda }_{V}}$ (DEV), связанной со скалярным полем как остаток от инфляционной стадии в ранней Вселенной. Было предположено, что энергия переменного компонента однозначно связана с плотностью вещества, образованного из скалярного поля на ранних стадиях эволюции Вселенной. Для простоты было рассмотрено расширение стандартной $\Lambda $CDM космологической модели с линейной связью между плотностями энергии в виде ${{\rho }_{m}} = \alpha (z){\kern 1pt} {{\rho }_{{{\text{DEV}}}}}$, где коэффициент пропорциональности, в общем случае, зависит от красного смещения. В данной модели нет проблемы HT, которая, предположительно, является следствием неточной экстраполяции при расчете современного значения постоянной Хаббла по наблюдениям флуктуаций реликтового излучения. Мы решили уравнение Фридмана при условии существования добавочного компонента темной энергии и его связи с веществом, и рассмотрели несколько функций $\alpha (z)$, при которых проблема HT устраняется. Современная плотность энергии DEV, необходимая для объяснения феномена HT, мала по отношению к космологической постоянной $\Lambda $, поэтому влияет на расширение Вселенной только при больших $z$, когда вклад постоянной Эйнштейна уменьшается.
В настоящее время ситуация противоположная, $\Lambda \gg {{\Lambda }_{V}}$, потому что уменьшение плотности вещества в процессе космологического расширения определяет переход от квази-Фридмановской стадии расширения к квази-Ситтеровской стадии. В случае постоянного коэффициента ${{\alpha }_{{{\text{HT}}}}}$ противоречие устраняется при ${{\alpha }_{{{\text{HT}}}}} \approx 24$. При этом плотность ${{\rho }_{{{\text{DEV}}}}}$ в настоящую эпоху соответствует относительному вкладу ${{\Omega }_{{{\text{DEV}}}}} \approx 0.014$. Наилучшее совпадение современного возраста Вселенной с модельным значением имеет место при слабом росте вклада ${{\rho }_{{{\text{DEV}}}}}$ с красным смещением и большим современным значением ${{\alpha }_{{{\text{HT}}}}} \approx 61$, что соответствует меньшему современному вкладу в плотность ${{\Omega }_{{{\text{DEV}}}}} \approx 0.0055$.
Современные значения параметров $\Lambda $CDM модели были определены из измерений флуктуаций реликтового излучения экспериментами WMAP и Planck. Данная процедура весьма сложна и основывается на нахождении экстремумов в многопараметрическом пространстве. В нашей модели, устраняющей проблему HT, если она действительно существует, возникают некоторые изменения процедуры нахождения космологических параметров, поэтому космологические параметры, полученные при этом, могут слегка измениться.
В нашей модели DM должна быть представлена широким спектром масс частиц, в отличие от модели холодной темной материи (CDM) с частицами одинаковой массы, которая обычно рассматривается. Принимая, что на ранних стадиях расширения существовало термодинамическое равновесие, число легких частиц DM должно быть близко к числу реликтовых фотонов. Зная относительный вклад в энергию CMB, можно грубо оценить, что наименьшая масса частиц DM не должна в настоящее время превышать значения ${{({{\Omega }_{{{\text{DEV}}}}}/{{\Omega }_{{{\text{CMB}}}}})}^{{1/4}}}k{{T}_{{{\text{CMB}}}}} \sim 5 \times {{10}^{{ - 4}}}$ эВ.
Список литературы
Y. B. Zeldovich and I. D. Novikov, Structure and Evolution of the Universe (Moscow: Nauka, 1975).
G. Riess, Nature Rev. Phys. 2 (1), 10 (2020).
W. L. Freedman, Astrophys. J. 919 (1), id. 16 (2021).
G. Bisnovatyi-Kogan, arXiv:2002.05602 [astro-ph.CO] (2020).
G. S. Bisnovatyi-Kogan, Universe 7 (11), 412 (2021).
D. N. Spergel, L. Verde, H. V. Peiris, E. Komatsu, et al., Astrophys. J. Suppl. 148 (1), 175 (2003).
P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, et al., Astron. and Astrophys. 594, id. A13 (2016).
N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, et al., Astron. and Astrophys. 641, id. A6 (2020).
A. G. Riess, A. V. Filippenko, P. Challis, A. Clocchiatti, et al., Astron. J. 116 (3), 1009 (1998).
S. Perlmutter, G. Aldering, G. Goldhaber, R. A. Knop, et al., Astrophys. J. 517 (2), 565 (1999).
A. G. Riess, L. M. Macri, S. L. Hoffmann, D. Scolnic, et al., Astrophys. J. 826 (1), 56 (2016).
A. G. Riess, S. Casertano, W. Yuan, L. Macri, et al., A-strophys. J. 861 (2), 126 (2018).
A. G. Riess, S. Casertano, W. Yuan, L. M. Macri, and D. Scolnic, Astrophys. J. 876 (1), 85 (2019).
K. C. Wong, S. H. Suyu, G. C.-F. Chen, C. E. Rusu, et al., Monthly Not. Roy. Astron. Soc. 498 (1), 1420 (2020).
W. Yuan, A. G. Riess, L. M. Macri, S. Casertano, and D. M. Scolnic, Astrophys. J. 886 (1), 61 (2019).
L. Verde, T. Treu, and A. G. Riess, Nature Astron. 3, 891 (2019).
C. A. Bengaly, C. Clarkson, and R. Maartens, J. Cosmology and Astroparticle Phys. № 05, id. 053 (2020).
E. Di Valentino, O. Mena, S. Pan, L. Visinelli, et al., Classical and Quantum Gravity 38, id. 153001 (2021).
T. Karwal and M. Kamionkowski, Phys. Rev. D 94 (10), id.103523 (2016).
E. Mörtsell and S. Dhawan, J. Cosmology and Astroparticle Phys. № 09, id. 025 (2018).
V. Poulin, T. L. Smith, T. Karwal, and M. Kamionkowski, Phys. Rev. Letters 122 (22), id. 221301 (2019).
W. Yang, S. Pan, E. Di Valentino, R. C. Nunes, S. Vagnozzi, and D. F. Mota, J. Cosmology and Astroparticle Phys. № 09, id. 019 (2018).
S. Vagnozzi, Phys. Rev. D 102 (2), id.023518 (2020).
E. Di Valentino, A. Melchiorri, O. Mena, and S. Vagnozzi, Phys. Dark Universe 30, id. 100666 (2020).
C. Umiltá, M. Ballardini, F. Finelli, and D. Paoletti, J. Cosmology and Astroparticle Phys. 2015 (08), id. 017 (2015).
M. Ballardini, F. Finelli, C. Umiltá, and D. Paoletti, J. Cosmology and Astroparticle Phys. № 05, id. 067 (2016).
M. Rossi, M. Ballardini, M. Braglia, F. Finelli, D. Paoletti, A. A. Starobinsky, and C. Umiltá, Phys. Rev. D. 100 (10), id. 103524 (2019).
L. Knox and M. Millea, Phys. Rev. D. 101 (4), id. 043533 (2020).
V. V. Luković, B. S. Haridasu, and N. Vittorio, Monthly Not. Roy. Astron. Soc. 491, 2075 (2020).
W. Kenworthy, D. Scolnic, and A. Riess, Astrophys. J. 875, id. 145 (2019).
E. Mörtsell and S. Dhawan, J. Cosmology and Astroparticle Phys. № 09, id. 025 (2018).
J. Sakstein and M. Trodden, Phys. Rev. Letters 124 (16), id. 161301 (2020).
A. Gogoi, R. Kumar Sharma, P. Chanda, and S. Das, Astrophys. J. 915, id. 132 (2021).
G.-B. Zhao, M. Raveri, L. Pogosian, Y. Wang, et al., N-ature Astron. 1, 627 (2017).
M. Mortonson, W. Hu, and D. Huterer, Phys. Rev. D. 80 (6), id. 067301 (2009).
X. Li and A. Shafieloo, Astrophys. J. Letters 883 (1), id. L3 (2019).
L. Parker and D. A. Vanzella, Phys. Rev. D. 69 (10), id. 104009 (2004).
G. Steigman, D. N. Schramm, and J. E. Gunn, Phys. Letters B 66 (2), 202 (1977).
L. Amendola, Phys. Rev. D. 62 (4), id. 043511 (2000).
M.-X. Lin, M. Raveri, and W. Hu, Phys. Rev. D. 99 (4), id. 043514 (2019).
W. Hu and I. Sawicki, Phys. Rev. D. 76 (6), id. 064004 (2007).
Einstein, Kosmologische Betrachtungen zur allgemeinen Relativitätstheorie, Sitzungsberichte der Königlich Preussichen Akademie der Wissenschaften zu Berlin. Reprinted in The Collected Papers of Albert Einstein, 1914–1917 6 (Princeton University Press, 1996).
A. Guth, The Inflationary Universe (Reading, Massachusetts: Perseus Books, 1998).
A. D. Linde, Phys. Letters B. 129 (3/4), 177 (1983).
A. Starobinsky, Phys. Letters B. 117 (3/4), 175 (1982).
V. F. Mukhanov, G. V. Chibisov, JETP 56 (2),258 (1982).
K. Arun, S. Gudennavar, and C. Sivaram, Adv. Space Research 60, 166 (2017).
D. Samart and P. Channuie, European Phys. J. C 79 (4), id. 347 (2019).
Д. С. Горбунов, В. А. Рубаков, Введение в теорию ранней Вселенной. Теория горячего Большого взрыва (M.: Ин-т ядерных исследований РАH, 2007).
А. А. Фридман, Успехи физ. наук 80 (7), 439 (1963).
P. J. E. Peebles, Principles of physical cosmology (Princeton University Press, 1993).
Age of the Universe, WikipediA (2021), in press https://en.wikipedia.org/wiki/Ageoftheuniverse .
И. С. Градштейн, И. М. Рыжик, Таблицы интегралов, сумм, рядов и произведений (М.: Гос. изд.-во Физ.-Мат. литературы, 1962).
Дополнительные материалы отсутствуют.
Инструменты
Астрономический журнал