Новые классы точных решений и согласованные модели космологической эволюции

УДК 524.834

А. В. Асташёнок, А. В. Юров, А. В. Япарова

НОВЫЕ КЛАССЫ ТОЧНЫХ РЕШЕНИЙ

И СОГЛАСОВАННЫЕ МОДЕЛИ КОСМОЛОГИЧЕСКОЙ ЭВОЛЮЦИИ

Развит метод построения точных космологических решений уравнений Эйнштейна, основанный на их представлении в форме линейного дифференциального уравнения второго порядка. Метод позволяет, в частности, использовать произвольное известное решение для построения более общего, параметризованного двумя константами. Показано, что в определенном случае возникающие новые классы точных решений содержат особенности, обладающие следующим свойством: геодезическая, начинающаяся или заканчивающаяся в сингулярности, имеет бесконечную длину. Такая сингулярность приводит к отсутствию горизонтов событий (если это сингулярность будущего). В этом случае, вероятно, можно построить космологическую модель, удовлетворяющую требованиям согласованности в течение всего времени существования вселенной.

15

A method for construction of exact cosmological solutions of Einstein equations is developed. The method is based on representation of these equations as a second-order linear differential equation. Particularly, the method allows to construct a more general solution, parametrized by two constants, using any known solution. Also it is shown, that in a specific case new classes of exact solutions have singularities with the following property: a geodesic, which begins or ends in the singularity, has an infinite length. Such a singularity leads to an absence of event horizons (if this is a future singularity). In this case it might prove possible to construct cosmological model, satisfying every consistency requirement during a universe lifetime.

Ключевые слова: уравнения Эйнштейна — Фридмана, скалярное поле, преобразование Дарбу.

Keywords: Einstein — Friedman equations, scalar field, Darboux transformation.

Введение

В работе [1] предложен подход для нахождения точных решений уравнений Фридмана. Основная идея выглядит следующим образом: в плоской вселенной функция, равная масштабному фактору в кубе у/= a3, удовлетворяет линейному дифференциальному уравнению второго порядка:

Ц- = 12*0 [р- X

М2 Г с

= 12xG\p-^\W.

© Асташёнок А. В., Юров А. В., Япарова А. В., 2018

Вестник Балтийского федерального университета им. И. Канта. Сер.: Физико-математические и технические науки. 2018. № 4. С. 15 — 48.

Обозначения стандартны: р — плотность, а р — давление материи. Если вселенная заполнена скалярным полем с минимальной связью ф с потенциалом самодействия У(ф) и вакуумной энергией с плотностью Л, то (1) оказывается в точности стационарным одномерным уравнением Шредингера

^ = {и-ХУг,

аг

где и(г) = 24лСУ, а «энергия» Я = -24яСЛ. Разумеется, потенциал V следует считать функцией времени: V(t) = V(ффt)) [1]. Конечно, для явного 16 представления потенциала в таком виде необходимо знать точные ре——' шения полевых уравнений, однако можно поступить иначе — задать явный вид и(г) как функции времени и начальные условия, после чего решить (2) и определить масштабный фактор. Если не фиксировать начальное условие, то масштабный фактор автоматически окажется зависящим от двух произвольных констант, одна из которых может быть поглощена трансляцией по времени, а вторая может дать весьма нетривиальную динамику. Эта идея частично реализована в работе [1] на примере трех моделей: гармонического осциллятора (и ~ г2), степенного (и ~ г-2) и односолитонного (безотражательного) (и ~ 1/сЬ2г) потенциалов. Главным выводом обсуждаемой работы является то, что для всех трех типов моделей возникает инфляция. Вывод: инфляция возникает при самых общих видах потенциалов.

Однако с выходом из инфляции ситуация обстоит куда проблематичнее1. Уравнение (2) в космологии исследовалось в работах [2] и [3] с помощью формализма преобразования Дарбу (ПД)2. Показано, в частности, что выбирая соответствующим образом функции, входящие в определитель Крама, можно получать решения с интересным асимптотическим поведением. Например, известно, что решения уравнения Шрёдингера, принадлежащие пространству функций, интегрируемых с квадратом I2, ПД преобразуют в функции из этого же пространства, если опорные функции не имеют особенностей и демонстрируют наличие перемежающихся нулей. Правильная асимптотика позволяет проконтролировать выход из инфляции, и потому это свойство ПД вызывает большой интерес в контексте конструирования интегрируемых космологических моделей в плоском пространстве.

Одним из самых интересных и неожиданных результатов, полученных при изучении космологических моделей с помощью (2), является то, что решения этого уравнения демонстрируют нарушение слабого

1 Так, авторы цитируемой работы предложили такую модификацию потенциалов, чтобы параметры оказались зависящими от температуры. Тогда выход на фридмановский режим можно описывать как фазовый переход в материи ранней вселенной.

В работе [4] аналогичная техника использовалась для построения точных решений на бране и в объемлющем пространстве, снабженном структурой ор-биобразия.

энергетического условия3, причем оно имеет место как следствие динамической эволюции в процессе расширения вселенной. В работе [5] этот удивительный феномен был назван эффектом гладкой фантоми-зации, а позднее его стали называть эффектом пересечения фантомной зоны. В упомянутой работе было показано, что пересечение фантомной зоны (ПФЗ) вполне естетственно может быть описано в популярных космологических моделях со скалярным полем. Другими словами, уравнения Фридмана демонстрируют совершенно новый физический эффект или (и такая точка зрения тоже имеет право на существование) патологическое поведение. В любом случае очевидно, что предложенный в [1] и развитый в последующих работах метод является весьма эффективным и плодотворным.

В настоящем исследовании мы предлагаем модификацию данного метода. Из уравнений Эйнштейна следует, что функция более общего вида ц/п = ап тоже удовлетворяет уравнению Шрёдингера с зависящим от п потенциалом (мы будем его обозначать ип), который оказывается некоторой линейной комбинацией давления (деленного на квадрат скорости света) и плотности. Задавая историю «потенциала», можно получать решения для функции ц% и, таким образом, для масштабного фактора. Ясно, что даже для одной временной зависимости ип(Ь) спектр решений существенно расширяется.

Мы упоминали выше, что наш метод применялся в плоском случае. При наличии ненулевой кривизны данный подход оказывается намного менее эффективным просто потому, что для произвольной степени п функции ц% = ап удовлетворяют уже нелинейному уравнению, а не уравнению типа (2). Для нелинейноых уравнений, получаемых путем такой процедуры, общие или эффективные методы точного интегрирования неизвестны. Это приводит к тому, что в общем случае не удается построить решение, параметризованное двумя константами интегрирования, а именно это необходимо для реализации процедуры, аналогичной проведенной в работе [1]. Однако для п = 1 ситуация упрощается.

В этом случае роль «пси-функции» играет непосредственно масштабный фактор. Действительно, одно из уравнений Эйнштейна в метрике Фридмана оказывается линейным уравнением второго порядка, причем роль потенциала (с точностью до множителя, пропорционального ньютоновской константе) играет величина р + 3р/с2, взятая с отрицательным знаком. Тогда физический смысл «потенциала» тоже становится понятным: при нулевой вакуумной энергии знак потенциала определяет, выполняется или нарушается сильное энергетическое условие. Это особенно полезно ввиду того, что нарушение сильного энергетического условия является необходимой (но не достаточной) причиной генерации инфляции. Из-за важности вопроса мы посвятим случаю п = 1 отдельный раздел. Изучение этих моделей приведет нас

17

3 В работе принята следующая терминология: мы говорим, что выполнено сильное энергетическое условие, если справедливы неравенства р + р/с2 > 0, р + 3р/с2 > 0. Слабое же энергетическое условие или условие энергодоминантности подразумевает выполнение р > 0, р + р/с2 > 0.

18

к еще одному замечательному наблюдению: построенные нашим методом двухпараметрические точные решения в чрезвычайно общем случае демонстрируют отсутствие горизонтов событий на фазе сжатия. Другими словами, геодезические в такой вселенной продолжаются неограниченно. Мы считаем, что это свойство является чрезвычайно важным для решения следующей амбициозной задачи — построение космологической модели, физически согласованной в течение всего времени своего существования.

Работа организована следующим образом. В разделе 1 описывается сам метод. Далее, действуя аналогично работе [1], мы рассматриваем примеры точных решений, используя для этого известные (и немногочисленные) потенциалы, для которых уравнение Шредингера может быть решено точно в аналитическом виде. В разделе 3 отдельно рассматриваются упомянутые выше решения, не содержащие горизонтов событий, и поясняется, почему с нашей точки зрения такие космологические модели представляют огромный интерес. В разделе 4 мы кратко описываем метод преобразований Дарбу, позволяющий строить новые решения. В последнем разделе подводятся итоги работы.

1. Процедура линеаризации

Система уравнений Эйнштейна — Фридмана имеет вид

а2 8я"С кс2 . .

-2 = —^Р--Т,

а 2 а

а 4жв ( 3р Л

Пусть функция а = а(г) для заданных р = р(г), р = ррг) является решением (3), (4), причем к = 0. Тогда ц/п = ап — решение уравнения

У = ип

то есть уравнения Шрёдингера со сложным «потенциалом» вида ип(г) = 8лв(п2р - (3п/2)(р + р/с2))/3.

В случае, когда вселенная заполнена самодействующим скалярным

ф2

полем ф и лагранжиан имеет стандартную форму I = —- - V (ф),

2с

и (п (п - 3)ф 2т/(ф)Л

ип =—ф + пТ (ф.

Случай, когда п = 3 и = 72лС/3, подробно исследован в [1].

Замечание 1. Если одновременно с полями материи во вселенной присутствует ненулевая космологическая постоянная (с рЛс2), то уравнение (5) выглядит как классическая задача на собственные значения:

У =

со спектральным параметром, определенным соотношением Хп = - 8лСп2рЛ/3. Постановка задачи на собственные значения и соб-

ственные функции для (6) требует дополнительно наличия однородных начальных (независимая переменная здесь время, поэтому мы и говорим о начальном, а не граничном условии) условий. Тот факт, что вакуумная энергия в такой форме выступает как собственное значение, оказался достаточно неожиданным, и не исключено, что на нашем пути удастся прояснить загадку вычисления этой величины (эта идея была впервые высказана в [1]). Очевидно, подобная гипотеза правомочна и в случае уравнения (6), и в разделе 2 мы обсудим эту интригующую возможность. Пока что заметим, что если известно частное решение (6), то используя (3) и (6), можно найти скалярное поле и потенциал в виде простой квадратуры:

Ф

V (* ) =

8яТ

2

V " " ЧУ п у у

п п2 V —

Соотношения (7), (8) позволяют вычислить V = У(ф). Разумеется, в явном виде это можно сделать лишь в редких случаях, поскольку интегралы в общем виде выражаются с помощью специальных функций, а процедура обращения функции вообще редко осуществима в явном виде.

Общее решение (5) имеет вид

% = — + фп,

где —/п — другое линейно независимое решение, то есть вронскиан, построенный на этих двух функциях, равен единице:

7 .1-11

Фп

; ¥2

Отсюда вытекает следующая теорема:

Линеаризационная теорема. Пусть а = а

1^

1/п л

С1 + С21 а2п V а У

удовлетворяет системе (3), (4) с новой плотностью энергии рп и новым давлением рп, причем

п2Рп - | Р + *") = п2р- | [р+

Замечание 2. Линеризационная теорема справедлива для случая к = 0. Если к = + 1, то следует положить п = 0,1.

Приведем несколько примеров точных решений.

19

ап = а

20

2. Генерация решений при фиксированном U„

В качестве простых «потенциалов» Un

m

Un

u"

ch2

и исследуем решения для потенциалов (a), (b), (c) на предмет существования инфляционных режимов и выхода из них.

(a) Решением уравнения (6) c нулевыми граничными условиями при t ^ + ю являются функции

ц/n = AnHs(/Ut) exp(-/t2/2n),

где A — постоянная; Hs(¿ut) — полиномы Эрмита порядка s. Соответственно динамика масштабного фактора имеет вид

a(() = AH1/" (/t )exp(- / 2¡2n).

В простейшем случае y/n = exp(-/t2/2). Рассмотрим этот случай (n = 1) подробнее, принимая во внимание, что при n = 1 k необязательно равно нулю.

В общем виде решение для масштабного фактора:

a

где t0 — момент времени, отвечающий наблюдению, сделанному в настоящее время; a0 = 1028 см — текущее значение масштабного фактора; параметр u выражается через to; наблюдаемое значение постоянной Хаббла Ho = 24,3 х 10-19 секунда-1: / = - Ho/to.

Это решение представляет самостоятельный интерес, поскольку содержит свободный параметр t0. Как следует из (3), верно соотношение

^ = H°

а (го) ^

С другой стороны, наблюдения свидетельствуют об ускорении вселенной. Выражая вторую производную,

;1=¥*<0<7рс>=^

Сравнивая (14) и (15), легко выразить г0 через Н0:

10

'0 -- н

Возраст вселенной T будем отмерять от момента времени tpl, когда a(tpi) = api = 10-33 см. За этот интервал вселенная совершила

N = ln = 61 ln 10 = 140,458

api

e-расширений. Полагая tpl = t0 - T и учитывая (16), получаем квадратное уравнение

3H02T2 + 20H0T - 20N = 0, откуда находим возраст вселенной в годах:

T = 2

3H0 H0

Это на два порядка больше реальной величины 15 — 20 х 109 лет. Тем не менее надо иметь в вицу, что формально возраст такой вселенной равен бесконечности, а 358 миллиардов лет — это промежуток времени, прошедший от момента, когда масштабный фактор имел планков-ский размер, до настоящего момента. Определение именно этого момента времени как «начала» выглядит естественным, но это не означает, что такой выбор действительно правилен.

Промежуток времени между современной эпохой

At1 = t1 -10 = 10 = 19,67 х 109,

110 3H0

а время от текущего момента до момента времени t = t2, когда имеет место большая остановка (a(t2 = 0)) составляет

At2 = 10 / (3H0) = 43,5 X 109 лет4.

Отдельно необходимо исследовать выполнение слабого и сильного энергетических условий. Очевидно,

p 3H02 A0 + 10kc 2e3H01V10

P+^T = ——-5-,

c2 40жСА^

p+ 3p =" 9 H2 ($H0¡t2 -10) c2 400жС ,

где Ао = аов5/3.

Из (18) следует, что слабое энергетическое условие нарушается при достаточно поздних временах 171 во вселенных с отрицательной кривизной (к = -1) и всегда выполняется в плоских и закрытых моделях.

4 Заметим, что дискуссии по поводу величины современного значения пара-

метра Хаббла Н0 = к х 0,324 х 10-19 обратных секунд пока не прекратились.

Наиболее распространена точка зрения, согласно которой к = 75. Именно этот

случай рассматривается здесь и в дальнейшем. Если использовать нижнее воз-

можное значение к = 45, то Д71 = 32,79 х 109 лет, а Д72 = 72,5 х 109 лет.

21

22

Ускорение имеет место в течение бесконечных интервалов в прошлом и будущем, на временах г > л/10/(Н(,л/3), а отсутствует лишь на конечном интервале времени. В современную эпоху г = 'о плотность составляет

( ) з н2 зкс2

Ро =Р('О Ь^ТГ +

8жС 8лСа2

Используя численные значения, находим, что р» = - 0,123 х 10-28 г/ см3 при к = - 1, ро = 0,38 х 10-29 г/см3 при к = 0 и ро = 0,2 х 10-28 г/см3 при к = + 1. Отсюда становится очевидным, что случай отрицательной кривизны неизбежно приводит к нарушению условия энергодоминантности в настоящее время.

Несложно вычислить и величину вакуумной энергии. Предположим, что вселенная заполнена самодействующим скалярным полем ф с минимальной связью и потенциалом V = У(ф) и — дополнительно — космологической постоянной с плотностью рл (часто вакуумную энергию отождествляют с минимумом потенциала, но мы ввели ее как самостоятельную величину). Тогда

1 = - (рф+ %1 + = < 2 - 2,.

а 3 ^ с ) 3

Отсюда следует, что, во-первых, рл = -9Н02/ (80^ ) = -0,3рс, а во-вторых, форма потенциала V(ф) имеет вид

V (ф) = 9Н0Ф1.

20с2

Значит, скалярное поле в нашей модели описывает невзаимодействующие скалярные частицы с очень маленькой массой — т ~ 0,3 х 10-60 г, что на 33 порядка меньше массы электрона. Кроме того, плотность вакуумной энергии оказывается отрицательной.

Но согласуется ли описанная простая модель с фактом существования галактик? На первый взгляд, наличие очень большой величины ускорения в прошлом должно препятствовать образованию наблюдаемой крупномасштабной структуры, что указывает на неадекватность нашей модели. Удивительно, но это не так.

Пусть рм — плотность барионного вещества, которая в настоящее время равна рм0. Эти два значения связаны друг с другом известным соотношением

рм = рмо

где 2 — величина красного смещения.

Ограничения на модуль ускорения получаются из требования, чтобы величина рп = 3а/(8жСа) была меньше правой части (19) при 2 = 2т, что соответствует моменту, когда формировались самые ранние галактики:

ру <Рмо

При хт = 4,5 получаем ру < 170рм0, а при хт = 10 (известная величина красного смещения, часто рассматривающаяся как верхний наблюдаемый предел) рр < 4000рм0. В дальнейшем мы будем использовать именно такую величину красного смещения.

Формула (13) позволяет вычислить промежуток времени, который потребовался бы свету, испущенному галактикой с величиной красного смещения 2т = 10, чтобы достигнуть нас в настоящий момент времени:

Т = 2 - 5 + У 25 +151п 11

= 3 Н0 ,

что дает Т = 24,43 х 109 лет. В момент испускания величина ру составляла

~ = 3Н02

р = 80яС ,

то есть 0,226 х 10-28 г/см3. В то же время плотность барионов — рм0 = 10-31 г/см3, следовательно рм0 (1 + 2т)3 = 0,133 х 10-27 г/см3. Значит, условие (20) выполняется для галактик с красным смещением 2 = 2т = 10.

Выше мы исходили из того, что «галактическая материя» — это ба-рионы. Если повторить вычисления для случая темной материи, составляющей примерно 0,3 от критической плотности, то получится рм0 (1 + 2т)3 = 0,42 х 10-26 г/ см3. Другими словами, и в этом случае условие (20) выполняется.

Можно убедиться, что (20) не приводит к верхнему пределу величины красного смещения, однако для темной материи существует нижний предел, который определяется условием выполнения неравенства

Дг) = 4 + 61п

Единственный вещественный корень трансцендентного уравнения Дг) = 0 равен 2 = 2» = 0,4556, а условие (21) выполняется при 2 > 2». Переходя от величины красного смещения к переменной «время» с помощью соотношения

= ^100 + 601п

7 = — ,

3 Н0

находим, что в модели (13) галактики перестали образовываться совсем недавно, примерно 0,356/Н = 4,6 х 109 лет назад. Интересно, что в рамках данной модели образование галактик прекращается не из-за слишком большой величины космологического ускорения, а из-за чрезмерной разреженности темной материи в настоящее время.

Разумеется, выбор решения (6) в виде волновой функции основного состояния произволен и иллюстративен. Можно, например, попробовать рассмотреть интегрируемые с квадратом (лежащие в пространстве I2) функции, отвечающие возбужденным состояниям. К сожалению, такие решения буквально неприменимы на всем интервале своей области определения. Действительно, в соответствии с осцилляционной теоремой волновая функция в-го возбужденного уровня имеет в нулей. Каждый из этих нулей означает наличие космологической сингулярности, причем для положительных значений п > 0 масштабный фактор

23

будет обращаться в нуль (сингулярность типа Big Bang или Big Chrunch), а если n < 0, то масштабный фактор будет взрываться, за конечное время приводя к сингулярности типа Big Rip. Разумеется, можно использовать решения для возбужденных уровней с номерами s > 1, ограничивая динамику определенным интервалом времения (между двумя соседними нулями). Такая вселенная будет, очевидно, начинаться и/или заканчиваться в точках сингулярности. Например, если в качестве решения (6) с «потенциалом» Un = i2t2 выбрать функцию ц/n ~ Hi(^i)exp(- mut2/2) и положить n = -1, то эволюция масштабного фактора в такой вселенной определится формулой

a

Это простое решение описывает нетривиальную динамику вселенной: в «начальный» момент времени t = 0 масштабный фактор равен бесконечности (имеем сингулярность типа Big Rip), затем происходит сжатие вплоть до момента t = ¡-г1/2, а после начинается неограниченно долгое инфляционное расширение. Используя (7) и (8), несложно вычислить асимптотику скалярного поля и потенциала при t ~ 0 и t ^ œ. Можно убедиться в том, что кинетический член плотности энергии скалярного поля отрицателен, то есть эволюция (22) соответствует фантомному полю, что неудивительно, поскольку сингулярности типа «большой разрыв» означают наличие фантомной компоненты. Асимптотика фантома и его потенциала самодействия выглядит так: а) t ^ 0:

W(t )'i sh ta Ь 1

У (ф) "о exP ("F70) , 7 = ;

2hGt0

б) t ^ œ:

W) "i Jr',

y (Фф t

~2еГ

Решению (22) соответствует отрицательная космологическая константа Л = -9i/ 8nG.

Можно рассмотреть решение вида

a(t) = а.Ьф^-^, 1 - ¡t

которое описывает вселенную, возникающую в момент времени t = -i1/2 из сингулярности Big Rip. Инфляция начинается в момент времени t = 0 и продолжается вплоть до коллапса вселенной, который имеет место при t = ¡г1/г. Эта модель приводит к очень похожим асимптотикам фантома его и потенциала при приближении к сингулярности Big Rip.

(b) В случае потенциалов типа (b) решение уравнения (6) для произвольного значения т > 0 и нулевого собственного числа Xn = 0 имеет вид

Wn = dfs + % + afs + V2,

где s = V1 + 4т/2 . В зависимости от выбора параметра n возникают различные варианты динамики масштабного фактора:

а) c1 > 0, c2 > 0, n > 0. Функция, описывающая масштабный фактор, положительна при t > 0 и обладает минимумом в некоторой точке t0 > 0, после прохождения которого динамика описывается степенной зависимостью с выходом на асимптотический режим a ~ t (2s+1)/2n. Инфляция возникает, если значения т превышают значения n(n - 1);

б) c1 > 0, c2 < 0, n > 0. Эволюция начинается из точки сингулярности при t = (- c2/c1)1/2s и сразу выходит на стадию степенной инфляции при условии m > n(n - 1);

в) ci > 0, С2 > 0, n < 0. Вселенная возникает при t = 0. Если т > n(n - 1), то имеет место фаза инфляции вплоть до определенного момента времени, после чего динамика переходит на фридмановский режим. В течение этого времени масштабный фактор достигает максимальной величины и начинает уменьшаться, асимптотически стремясь к нулю.

Изучим случай, соответствующий выбору c2 = 0 и n >0. Без потери общности ограничимся выбором n = 1. Решение имеет вид

a = Cts+1/2,

где C — положительная константа.

Решения для скалярного поля и потенциала в случае с положительной кривизны имеют вид

ф{\. )=^0 ±

cV 1 + 2s

(2s - 1)У2жС

it

2s-1 + а2 1

+—ln

,s-1/2

•Tt4t

2s-1 + а2 - ts

ЛЛ

Jt-ft

2s-1 +a2 + ts

V =

8nG

12s2 + 8s +1

4t2

25

где фо — константа интегрирования; также обозначено а2 = ^ 2 '

Видно, что в явном виде получить зависимость V = У(ф) невозможно. Однако можно рассмотреть асимптотику для № ~ 1 << а2 и № ~ 1 >> а2. Имеем:

а) при № " 1 << а2:

Ы1 + 2э

Фа ±

(2s -1)2kG t

а

.s-1/2 ,

V

4kCG

(2s - 1)^lnG c (1 + 2 s )а

(±(Ф-Фф )2 s-1;

2

2

б) при t2s 1 >> а2.

26

ф ф х cj 2s +1 t

ф *Фо ln Т.'

V(ф)-( + 8s + ^ e^ Ф" У ' 32 kG

где введены обозначения t. = I — I и ß =

а \ 2s-1 „ 4yllxG

. 2 ) cV 1 + 2s '

В плоском случае легко найти явную зависимость потенциала от скалярного поля V = V( ф).

, (t) = + С>/1 + 2s j^t

ф () + 2V2KG t0'

, ч (12s2 + 8s +1)) V(Ф )=( 8kG ' e'

Здесь Ь0 — константа интегрирования, а 3 определена в предыдущем примере.

Рассмотрим случай другого знака спектрального параметра. Пусть Яп = ¿р- > 0. Используя подстановку у/п (^, можно свести соответствующее уравнение Шредингера для функции у/п к уравнению Бесселя для новой функции Д^. Как и раньше, рассмотрим простейший случай п = 1. Нас будет интересовать решение, содержащее сингулярность в начальный момент времени, который выберем за нуль, то есть я(0) = 0. Получаем

я = СЛ],

где ] — функция Бесселя первого рода порядка в, а параметр в связан с / так же, как и в случае Я = 0.

Решение (26) описывает вселенную, которая начинается с сингулярности и ускоренно расширяется до момента времени г = V452 —1/2^. На следующем этапе масштабный фактор достигает максимума, после чего вселенная коллапсирует в финальную сингулярность.

Рассмотрим (26) для в = 3/2 (это означает, что / = 42). Как известно, функция Бесселя порядка 3/2 может быть выражена через тригонометрические функции. Обозначая ^ = х, запишем решение следующим образом:

я = С

Масштабный фактор достигает максимума при х « 2,75, финальная сингулярность наступает при х = х/ « 4,5.

Интересно рассмотреть асимптотическое поведение скалярного поля (как функции времени) и потенциала (как функции поля) вблизи начальной и финальной сингулярностей:

а) X ~ 0:

(х) ~ —7~т== 1пх/

У

"2^/ЛG

5Л

7г =

77с '

б) X ^ X/:

ф

-1п(п х/х - 008 X),

2\ лО

0,241 7 7 4Л

72 =-

лО с

(с) Решение для потенциала (с) при Лп = -Л2 ^ 0 (это соответствует неотрицательному значению космологической константы) имеет вид

ц/п = сг

Если Л = Ы, то (28) существенно упрощается:

= С

Цп = Мы).

Если п > 0, то, как и в случае «потенциала» (я), получаем инфляцию на полупрямой (- <», ^), причем ^ — точка перегиба функции сК^Ы).

Интересный набор решений можно получить в простейшем случае нулевого собственного значения Лп = 0:

Цп = СШ

Если п > 1, то решение описывает «рождение» вселенной в момент t = 0 с последующим асимптотическим переходом в стационарное состояние, причем я < 0 в течение всей эволюции. Если же п лежит в интервале 0 < п < 1, то реализуется инфляция с естественным выходом, после чего вселенная асимптотически стремится к стационарному состоянию. Используя формулы (7) и (8), определим асимптотическое поведение скалярного поля и потенциала:

а) t ~ 0:

ф

24плй

3 - п

1п (Л 2),

32лGn2

2 сТ70ф

ЛЛ е

70 =-

4л/ лпв

б) Ы >> 1:

Ф( )* +

2с

ЛплС

-Ы

У (Ф)*-Лф

27

е

28

В случае Л > Ло, сг > 0, с2 = 0 имеется решение Цп = сг

Такое всюду регулярное решение достаточно необычно, потому что может описывать две последовательные инфляции. Можно показать, что при определенных значениях параметров Л0, Л и п (30) демонстрирует вселенную, испытывающую стадию инфляции на некотором конечном интервале (0, £1). После этого реализуется выход из стадии инфляционного расширения, но начиная с некоторого момента времени Ь2 > Ь вселенная вновь оказывается в режиме инфляции. Чтобы конкретизировать условия, при которых реализуется такая ситуация, рассмотрим я в начальный момент времени, соответствующий Ь = 0. Используя (30), находим

я(0) ~ (1 -п)у4 -2у2 +1,

где введено обозначение у = Л/ Л< 1. Значит, вторая производная в нулевой момент времени будет неотрицательной, если выполняется условие

_ 2 2 1 -л/и 0 < у2 < Уо =-7-^ • 1-п

Дальнейшее рассмотрение этого вопроса приводит к заключению, что в случае п > 1 вселенная сразу выходит на стадию инфляционного расширения. Если у2 > у^, то решение (30) демонстрирует появление инфляции начиная с определенного момента времени Ь0 > 0. Если п < 1, то в достаточно узком интервале значений у0 - А < у2 < у^, А << у0 имеет место быть область расширения с отрицательным ускорением, но после этого вселенная вновь становится инфляционной. Скажем, при п = 0,25, Ло = 0,8Л вселенная расширяется с ускорением в течение конечного интервала времени (0, 0,08/Л), затем следует конечная область степенного расширения с длительностью порядка 0,4/Л. На временах Ь > 5/Л масштабный фактор уже меняется фактически по деситтеров-скому закону я ~ е4^.

Заключая этот раздел, вычислим асимптотическое поведение потенциала скалярного поля на самой ранней стадии. Используя (7), (8) и ограничиваясь точностью до первого порядка (по времени), находим при Ь << 1/ Л:

Ф(£ )~Ф0 —т^л ль ,

2ы пл

и с)" 4

Как мы видим, если «первичная инфляция» реализуется, то она соответствует режиму медленного скатывания (что естественно), а сам потенциал медленно падает по линейному закону.

3. Космологические модели без горизонтов событий 3.1. Метод

Сейчас мы продемонстрируем, что изученный выше метод позволяет строить при п = 1 космологические модели, не содержащие горизонтов событий.

Очевидно, при п = 1 уравнение (5) совпадает с (4). Пусть а(4) — некоторое известное решение (4). Будем метить «крышечкой» второе, линейно независимое решение с тем же потенциалом — Я(4) = v(t~)а(), причем

4 й4'

(И а2

я2(Г)'

Значит, общее решение (4), которое мы обозначим А(4), имеет вид

А

где С1,2 — две произвольные константы.

Такие функции являются двухпараметрическими обобщениями исходных решений (случай С2 = 0 возвращает нас к исходному решению). На втором шаге подставим (31) в (3) и вычислим плотность энергии в такой двухпараметрической модели. Поскольку при всех этих преобразованиях величина р + -31т не меняется (помним, что это функция вре-

с

мени)5, легко найти давление р = р(4) и окончательно получить обобщенное уравнение состояния.

Пусть момент времени 4/ является нулем функции (31), причем будем предполагать, что 4/ не является нулем исходной функции я(4): а(4/) = а/ ф 0. Тогда

4V/=—

С2

В эпсилон-окрестности точки 4/ (4 = 4/ — е при е~ 0(1)) имеем простое соотношение (используем обозначения А/ = А(4/), А/ = А (4) при 4 = 4/):

С

А

Поскольку линейный элемент имеет вид

йз2 = с2й42 - А2

для направленных в будущее нулевых радиальных геодезических получаем

сй4 = са/ / й4 (а-Л

Г сй4 сиг г й4

ЛТ = —тг =-— -= - да X

4^А(4) С2 4^ 4-4( ё

С

29

4{-е '} V-2

5 Важно отметить, что если в исходной модели выполнялось сильное энергетическое условие, то оно выполняется и в новой.

Отсюда следует, что если sgn(я// С2) = - 1, то Дг = + да: нулевые геодезические совершают бесконечно много «оборотов-осцилляций» (в замкнутой вселенной) и лишь потом попадают в финальную сингулярность при Ь = Ь/. Поскольку вселенная с линейным элементом (33) однородна и изотропна, можно доказать, что единственной сингулярностью будущего будет с-граница в вице точки, достигаемой при Ь = Ь/. Значит, такое пространство-время не содержит горизонтов в будущем. Кроме того, параметр уравнения состояния т = р/(рс2) ^ -1/3 при Ь ^ Ь/. То, что параметр оказывается отрицательным, — хорошо известный факт, следующий из анализа динамики коллапсирующих вселенных с использованием глобальных методов.

30 Если вселенная заполнена скалярным полем ф с потенциалом У(ф) и

минимальной связью (причем вклад ф в полный тензор энергии-импульса материи является доминирующим при Ь ^ Ь/), то можно показать, что в окрестности точки финальной сингулярности Ь = Ь/ выполняются следующие соотношения:

1

ф = ф(0) + — 1п

V

1 - ±л

V ] ,

V (ф) = -г1те-2т (ф-ф(0)),

с V Ь ]

где

2«2Л/жС 2 С2 v = +—. , и =---.

с^ / + кс2 а/

Очевидно, при Ь ^ Ь/ выполняется vф^ - да, а значит, V ^ + да.

Будем считать, что масштабный фактор а(Ь) обращается в нуль в некоторый момент времени, предшествующий Ь/. Не теряя общности, выберем этот момент за Ь = 0. Нас будет интересовать случай вселенной, возникающей из начальной сингулярности при Ь = 0 и заканчивающей существование в финальной сингулярности при Ь = Ь/. Это имеет место быть при условии регулярности выражения а(Ь)и(Ь) при Ь ^ 0. Сделаем два дополнительных предположения: (г) при достаточно малых положительных значениях времени Ь работает баротропное уравнение состояния р = тс2р; (и) выполняется слабое энергетическое условие или условие энергодоминантности т > - 1. Тогда выражение а(Ь)и(Ь) остается регулярным при Ь ^ 0, если т >- 1/3 (например, модель вселенной, заполненной барионным веществом (т = 0) или электромагнитным излучением (т = 1/3)). В случае т >1/3 при Ь ^ 0 в (31) преобладает первое слагаемое, а если - 1/3 < т < 1/3, то второе. Например, если начинать нашу процедуру со вселенной, заполненной холодным веществом (ба-рионами или темной материей), то общее решение при Ь ^ 0 будет вести себя как А ~ Ь1/3, то есть описывать динамику вселенной, в которой на начальной стадии доминирует материя с так называемым предельно жестким уравнением состояния т = 1, а не с уравнением состояния пылевой материи с т = 0.

3.2. Условие самосогласованности и пространство-время без горизонтов

На этом этапе может возникнуть вопрос, почему мы акцентируем внимание на моделях вселенных без горизонтов событий и демонстрирующих смену расширения на режим сжатия с коллапсом в финальную сингулярность? Ответ: эти модели удовлетворяют «условию самосогласованности». Поясним этот момент.

Традиционно для космологии считать, что обнаружение каких-либо несогласованностей, возникающих при сравнении известных физических закономерностей с данными космологических наблюдений, свидетельствует о существовании пока не открытых фундаментальных принципов. Скажем, необходимость чрезвычайно тонкой настройки начальных условий в ранней фридмановской космологии рассматривалась как свидетельство того, что традиционная физика просто непреме-нима вблизи космологической сингулярности. Аналогично выглядела проблема однородности и изотропности реликтового фона. Позднее все эти проблемы были успешно разрешены в рамках инфляционной парадигмы (включая объяснение происхождения крупномасштабной структуры).

Тем не менее можно предложить другой, значительно менее принятый подход, а именно считать основные физические закономерности уже заданными и использовать упомянутые несогласованности для отбора космологических моделей, свободных от этих несогласованностей. Выразимся яснее: принцип самосогласованности означает, что глобальная геометрия пространства-времени определяется не только распределением вещества, но и отсутствием несогласованностей, отсутствием нестыковок, возникающих при применении различных физических методов для описания космологических процессов. При этом изучение глобальной эволюции вселенной становится принципиально важным просто потому, что если динамика вселенной происходит согласованно в течение какого-то времени, нет гарантий, что несогласованности не проявятся позднее. Например, с точки зрения квантовой теории поля фантомы оказываются полями с отрицательными кинетическими членами, что приводит (на квантовом уровне) к потере унитарности. Вместе с тем, как мы показали выше, даже если фантомов не существует в начальный момент времени, они вполне могут появиться динамически. Значит, если исследования покажут, что фантомы действительно несовместимы с теорией поля, следует наложить на динамику вселенной условие, препятствующее их образованию в будущем. Важно отметить, что, по-видимому, почти все пространство-время вселенной лежит в будущем относительно текущего момента, поэтому изучение глобальной эволюции означает главным образом изучение будущего.

О каких еще несогласованностях (кроме фантомов) идет речь? Данные, полученные при наблюдениях над удаленными сверхновыми типа 1а, которые являются «стандартными свечами», показали, что вселенная испытывает ускоренное расширение [8; 9]. В рамках примени-

31

32

мости уравнений Эйнштейна и выполнения слабого энергетического условия следует предположить наличие материи особого вида (темная энергия), которая может представлять собой или скалярное поле неизвестной природы (квинтэссенция), или космологическую постоянную, или комбинацию этих двух объектов. Измерение параметра адиабатич-ности довольно уверенно свидетельствует о втором варианте [10].

В этом случае мы наблюдаем процесс смены динамического закона расширения от фридмановского (степенного) к деситтеровскому (экспоненциальному). В случае точной вселенной де Ситтера (дБ) вся будущая эволюция наблюдаемой вселенной должна быть ограничена горизонтом событий Я = с/Н (напомним, что Н — параметр Хаббла, который в дБ-вселенной является постоянным). Этот факт влечет за собой ряд крайне сложных проблем в совершенно неожиданной области. К примеру, активно обсуждается вопрос о том, как формулировать фундаментальную теорию (теорию суперструн / М-теорию) в конечном объеме [11 — 13]. Дело в том, что традиционный аппарат квантовой теории поля основан на формализме Б-матрицы, а Б-матрица, как представляется, не может быть математически корректно определена в этой ситуации. Очень вероятно, что непротиворечивое описание появляется лишь при Я = да. Таким образом, появляется пример несогласованности между космологией и теорией поля.

Остроумный выход (в духе определения согласованной модели) был предложен в [14]. Автор пришел к заключению, что темная энергия должна не быть космологической постоянной, а реализовываться в рамках моделей с «голографической» фантомной энергией. Отталкиваясь от модели, описанной в [15], ученый показал, что вполне возможны ситуации, когда Я = да. Другими словами, голографическая фантомная энергия (но весьма специального вида) оказывается более предпочтительным кандидатом на роль темной энергии с позиций возможности непротиворечивой формулировки фундаментальной М-тео-рии. Это прекрасный пример, иллюстрирующий смысл введенного нами выше условия согласованности или самосогласованности.

Однако фантомы выглядят (по крайней мере, пока) слишком экзотически, чтобы так легко сделать вывод об их существовании, поэтому разумно оставить голографический фантомный сценарий на «черный день» и поискать другие модели. Очевидно, кандидатом (разумеется, если мы принимаем принцип согласованности) остается только квинтэссенция. Приведем еще один аргумент против космологической постоянной: неограниченное расширение вселенной приводит к парадоксальному выводу, впервые отмеченному в работе [16]. Оказывается, если возраст деситтеровской вселенной с современным темпом расширения превысит Т = 1060 лет, то в такой вселенной на протяжении всей ее бесконечной эволюции будут доминировать так называемые больцманов-ские мозги, а не обычные наблюдатели. Можно было бы сделать вывод, что вселенная должна существовать конечное время. Собственно, теория струн предсказывает конечное время жизни, связанное с тем, что состояние «нашего вакуума» метастабильно. Однако приведенная выше [16] оценка максимального времени жизни вселенной чрезвычайно

отличается от оценки времени жизни тметастабильной деситтеровской фазы. Согласно исследованиям [17; 18], эта величина имеет характерный

порядок т < е°'5х10 лет. В некоторых моделях это число удается существенно понизить до е и даже до е лет (например, в модели, содержащей КРУ-инстантоны [19]), но и в этом случае время жизни вселенной чудовищно больше, чем 1060! Аналогичные трудности будут, видимо, возникать и в фантомных моделях, для которых = да6. Мы заключаем, что эти результаты могут служить веским аргументом против моделей вселенных с вечным расширением. Отсюда следует, что более предпочтительными являются сценарии, в которых современная фаза расширения должна смениться на фазу сжатия с последующим коллапсом, причем время жизни вселенной не должно превышать предельных значений, полученных в [16]. Таким образом, наиболее привлекательными представляются модели, для которых выполняются два условия: (г) вселенная оканчивает свое существование в финальной сингулярности на стадии коллапса; (и) вселенные не содержат горизонта событий К < да. В дальнейшем будем называть согласованными в будущем (пишется без кавычек) или просто «согласованными» (пишется в кавычках) космологические модели, удовлетворяющие этим двум условиям. Это означает, в частности, что темная энергия должна быть квинтэссенцией, а не космологической постоянной.

Необходимо отметить чрезвычайно важное обстоятельство: во вселенной, заполненной веществом, для которого выполняются все энергетические условия, приближение к финальной сингулярности будет иметь осцилляторно-хаотический характер. В наиболее интересной и изученной модели — Бианки-1Х — сжатие демонстрирует постоянные смены «казнеровских эпох», причем характер этих смен быстро приобретает стохастическое поведение7. В таком режиме главный вклад вносят зависящие от пространственной кривизны члены, отвечающие плотности энергии, убывающей как шестая степень а6, поэтому при а ^ 0 имеем бесконечное число смен «казнеровских эпох» вплоть до достижения вселенной точки финальной сингулярности. Это приводит к существенному осложнению: в этих условиях вероятность случайной реализации вселенной с единственной с-границей в виде точки (неизбежное следствие условия отсутствия горизонтов) оказывается пренебрежимо малой величиной. Вместе с тем эти выводы кардинально меняются при наличии безмассового скалярного поля (например, дила-тона), в присутствии которого колебательный режим приближения к сингулярности сменяется монотонным. Значит, принцип согласованности в будущем требует наличия переменного безмассового скаляр-

6 Интересно отметить, что парадокс из [16] не появляется в фантомных космологических моделях общего вида, в которых Кь ^ 0 при неограниченном росте масштабного фактора [20]. Однако эти модели нас не интересуют по вышеизложенной причине.

7 В течение каждой эпохи расстояния вдоль двух осей осциллируют, а расстояние вдоль третьей монотонно уменьшается.

33

ного поля. При этом условие отсутствия горизонтов в замкнутой вселенной выполняется, если при приближении к сингулярности параметр т > -1/3.

Это замечательный вывод! Мы можем заключить, что описанный в разделе 3.1 математический подход построения точных решений космологических уравнений Эйнштейна — Фридмана приобретает весьма глубокий смысл, поскольку в общем случае он дает решения, обладающие именно таким поведением, которое требуется для выполнения «условий согласованности». Ниже, в разделах 3.3 — 3.6 мы продемонстрируем конкретные примеры таких решений.

__В качестве начальных (затравочных) решений будем выбирать

34 масштабные факторы, соответствующие различным значениям параметра уравнения состояния т и кривизны к. Наиболее интересны и физически осмысленны ситуации, когда исходное решение соответствует барионному веществу и излучению. Для таких моделей сильное энергетическое условие, очевидно, всегда выполнено, а значит, как мы уже отмечали, оно будет автоматически выполняться и для их двухпара-метрических обобщений. Для пространств с положительной и нулевой кривизной это приводит к выполнению и условия энергодоминантности, поэтому в дальнейшем мы не будем подробно исследовать слабое энергетическое условие. Удобно ввести следующие обозначения: если в качестве затравочного решения берется решение уравнений (3), (4) для пыли или/и излучения, то построенную по описанной схеме модель будем называть В-, И- или ВИ-моделью соответственно.

3.3. Б-модели

Для вещества т = 0, и уравнения (3), (4) принимают следующий вид:

а2 8жСС кс2 ,„_.

„2 а„3 2

а 3а а а 4nGC

а=~ юг,

а за

где C — положительная константа.

Замкнутое пространство

В случае положительной кривизны (k = 1) решения уравнений (35) и (36) имеют виц

а = х2 sin2 в,

t = Х-

где х = — константа. Опуская численные коэффициенты,

приведем линейно независимое решение:

а ~ sin2 0Í—d0--sin0cos0 .

J sin2 0

Отсюда находим общее двухпараметрическое решение:

А = 8тв8т

где константа 8 удовлетворяет условию 0 < 8 < л. Вторая константа входит мультипликативным образом, перенормируя величину %, и может быть поглощена за счет переопределения С. Параметр в в решении (39) связан со временем по формуле (38). Заметим, что положив в уравнении (39) 8= л, мы придем к (37). Масштабный фактор достигает макси-

,2

мума Amax = %2 sin2[ — | при tmax = -%(8 - sin 8). В момент времени

V 2 ) 4c

tf =-2-

масштабный фактор обращается в нуль. Поучительно вычислить долю времени от полного времени своего существования, в которой вселенная расширяется у= tmax /tf. Получаем

8- sin8 ....

у =-.

28 - sin 28

При малых 8^ 0 вселенная расширяется приблизительно У 8 часть своей жизни. При 8^ 7г получаем 50 % (у= Уг). Величины плотности ри давления p в параметрическом вице:

C 4sin4 в + sin2

Р=-

%6sin6 в 4sin2

Р =

c 2С

f 4sin4 в + sin2 (8-2в)) 4sin2 (8-в

3% sin6 в

1 -

Параметр уравнения состояния меняется следующим образом: в окрестности начальной сингулярности w ~ 1, а при приближении к финальной сингулярности (как и должно быть) w —1/3.

На этом примере покажем, что пространство-время действительно не содержит горизонтов. Рассмотрим распространение светового луча в B-модели для случая замкнутого пространства. Простоты ради выберем радиальную изотропную геодезическую

ds2 = c2dt2 - A(t)2dr2 = 0.

Расстояние, которое луч света пройдет за время t:

\ dt r = c — .

Jo A

Соответствующий путь, который луч успевает пройти за все время (конечное) жизни вселенной, равен

— ' в

r = 2 f de —sin-- = - 2(в cos 8 + sin 8 ln sin(8 - в))— .

0 sin

35

36

Таким образом, при 8= я длина пути конечна, но это специальный случай, поскольку для всех остальных 8ф я интеграл расходится. В этом, по сути, общем случае светоподобные геодезические огибают вселенную бесконечное число раз, вплоть до коллапса вселенной. Из-за высокой симметрии (изотропности и однородности) модели всегда можно выбрать координатную систему таким образом, чтобы две произвольные пространственные точки были соединены посредством пространственного сегмента изотропной геодезической. Это означает, что все бесконечные времениподобные кривые определяют одну и ту же с-границу, которая представляет собой точку, и горизонты событий отсутствуют.

Открытое пространство

Необычно, что в открытых моделях с к = -1 существуют решения с финальной сингулярностью. Действуя по схеме, описанной в 3.1, получаем в параметрическом виде решение уравнений (35) и (36):

А = ех2$ЬвМ8-в),

%2,л \

Ь = А- (Ьв- 20),

причем параметр е = 1 при 8> 0 и е = - 1 при 8= 0. В последнем случае получаем известную открытую фридмановскую модель. Если параметр 8 достаточно велик, то большая часть времени существования вселенная переживает стадию коллапса. Уравнение состояния:

р = 1 (_ 48Ь2 (8-0

т = 2 рс2 3

8Ь2 (8-20)-48Ь4 в

Из (42) следует картина, аналогичная случаю пространства с положительной кривизной, а именно: т ^ 1 при Ь ^ 0 и т ^ -1/3 при Ь ^ В некоторый момент времени, предшествующий достижению максимального значения масштабного фактора, параметр т терпит разрыв, поскольку в нуль обращаются плотности энергии р. Горизонты событий в такой вселенной тоже отсутствуют, в чем легко убедиться, вычислив путь луча света вдоль радиальной геодезической.

Плоское пространство

Это самый простой случай, для которого решение выписываются явно: А =

1 4

т =---1—

3 3

( ^1/3 - £ 1/3 А Ь1/3 - 2Ь1/3

V 1

Из (43) следует, что в момент времени Ьтах = Ь// 8 значение масштабного фактора достигает максимального значения и в этот же самый момент параметр т обращается в бесконечность (см. (44)).

2

Поведение параметра уравнения состояния тоже неоригинально: в начальный момент времени т = 1, а в финальной сингулярности т = -1/3.

Несложно понять, что обращение функции т в бесконечность будет общим свойством всех моделей, для которых кривизна пространства равна нулю или отрицательна и для которых существует финальная сингулярность, поэтому ниже мы специально это не оговариваем.

3.4. Я-модели

Для излучения т = 1/3, и уравнения (3), (4) принимают следующий вид:

а2 8 л ОС кс2

37

„2 п „4 „2

а 3а а

а 4лО

^ = ~~Т" +,

а 3 ^ ас )

где С — положительная константа. Для удобства мы использовали те же самые обозначения, что применялись в разделе 3.3, несмотря на то, что размерности этих констант различны. Как и раньше, последовательно рассмотрим три случая: к = 1, -1, 0.

Замкнутое пространство

В случае положительной кривизны (к = 1) общее решение системы (45) и (46) в параметрическом виде имеет следующий вид:

А ¡1 , 1 - сс8в= —,

где а (здесь и далее) — положительная константа. Конечная, на стадии сжатия, сингулярность возникает в момент времени

= / = ехр(2/а) , / с 1 + ехр(2/ а)'

а параметр уравнения состояния

1 2 (1

т = —+-

3 3 ( ( в

1 + 1 со8в1 1 -а1ntg— -а

ведет себя так: на ранних временах (малые в) т ^ 1/3, то есть во вселенной доминирует излучение, затем т в определенный момент времени обращается в нуль и наконец достигает значения -1/3 в финальной сингулярности. Путь луча света вдоль радиальной геодезической за время существования вселенной равен

38

2arctg(exp(l/a)) ^д

r = J т~гд=+"

о 1 -aln — 2

Вновь убеждаемся, что горизонты отсутствуют.

Открытое пространство

Случай k = -1 наиболее богат и интересен. Дело в том, что из него можно получить по крайней мере три класса решений, в том числе и стационарное. Все эти решения не содержат горизонтов событий и потому относятся к классу «согласованных». Стартовое решение уравнений (45) и (46) для k = - 1 имеет вид

a = zshв, chg-1 = —.

Линейно независимое решение a находится по (48) следующим образом:

a ~ shgf-dg= shginthд J shg 2

Общее решение:

]

A = 2shgin-^, chg-1 = —.

a 2

Без потери общности можно считать, что s = 0 для исходного решения и s = 1 во всех остальных случаях, к рассмотрению которых мы и переходим.

Сингулярные R-модели в открытом пространстве

Таким решениям соответствуют значения а, превышающие единицу. В момент времени

t = 22 1

f c а2 -1

масштабный фактор обращается в нуль, а при приближении к финальной сингулярности параметр w ^ -1/3.

R-модели с вечным расширением

При а< 1 вселенная вечно расширяется. При больших t масштабный фактор растет линейно со временем:

А = -1п( а)с*.

Линейность расширения указывает, что ш асимптотически стремится к -1/3. Очевидно, горизонт событий отсутствует, так как путь света вдоль радиальной геодезической

=1-

йв

Г вЛ

0 1п

2

а

стремится к бесконечности, поскольку при больших величинах в подынтегральное выражение стремится к постоянной величине.

Стационарная И-модель

Положив в (49) а = 1, получим, что при t ^ ж А ^ %, то есть имеем стационарное решение. Любопытно, что, как и в случае открытых решений, на больших временах ш ^ -1/3. В самом деле, при t ^ ж с с2С

р ^---, а р ^ —4-. Горизонт событий также отсутствует (расчет ана-

2 3%

логичен случаю открытого решения).

Плоское пространство

Случай нулевой кривизны как всегда наиболее прост для анализа. Решение для масштабного фактора записывается в виде

Г I л

А = У 2%а 1 - а1п — ,

V ^о)

где Ь0 — константа, введенная из соображений размерности. Без потери обшцости можно принять а = 1, так как качественно это не меняет поведение решения. Момент времени Ьтах, когда масштабный фактор имеет наибольшее значение, и время достижения финальной сингулярности tf связаны простым соотношением

тх = ехр(-2^

из которого следует, что в фазе расширения вселенная проводит порядка V 7 части своего существования. Горизонт событий отсутствует.

3.5. ЕЯ-модели в плоском пространстве

Теперь рассмотрим случай, когда в качестве исходного взято решение, описывающее плоскую вселенную, заполненную и излучением, и веществом:

С

а = С 8Ь2 в,

Р ) С

сЬв(сЬ2 в/3 -1) =—%с*.

Р1

39

40

При получении решения (51) мы полагали, что плотность энергии вещества меняется как р = р]/я3, где р — константа, а излучения ри = С/а4. Применив нашу схему, находим общее решение для данного р + 3р/с2 в виде

С сШ£-А = С эЬ2 — 1п-2,

р1 а

причем время связано с параметром в, как и в решении (51). При е= 0 получаем исходное, затравочное решение, а в остальных случаях без потери общности можно положить е= 1. Вновь видим, что имеется три класса таких решений.

БИ-модели с асимптотикой

В этом случае а < 1. На достаточно больших временах асимптотическое поведение масштабного фактора имеет вид

А ^ 1паС8Ь2 в~ t2/3. Р1

Для простоты мы опустили все несущественные численные коэффициент^!. С точностью до констант асимптотика (52) при а < 1 совпадает с асимптотикой (51). Отличие между решениями проявляется лишь на начальном этапе, когда существенен логарифмический множитель.

БИ-модель с асимптотикой

В этом случае а = 1, А ^ 2Ср эЬ в/р ~ № при t ^ да.

Интересное поведение в этих двух моделях имеет функция т. В самом начале т ^ 1/3, затем увеличивается, а потом уменьшается до 0 для первого случая и до 1 во втором случае. Заметим, что горизонты событий в этих двух моделях отсутствуют, так как путь луча света вдоль геодезической за бесконечное время стремится к бесконечности:

г ~ I ав

1пйЬ—- 1па 2

Однако в данном случае это не так важно с точки зрения «согласованных» моделей, поскольку в затравочной модели горизонты событий уже отсутствовали.

Сингулярные БИ-модели

Для а > 1 решение будет заканчиваться сингулярностью. Нетрудно убедиться, что горизонт событий в таких моделях будет отсутствовать.

3.6. Л-модели

Особый интерес представляют модели, построенные на стартовых решениях, содержащих ненулевой Л-член (ю = -1). Ниже мы представим общую схему задачи и рассмотрим один частный случай.

Пусть вселенная с ненулевой вакуумной энергией Л заполнена еще и «материей», плотность которой рЕ убывает по степенному закону:

С

Ре =—,

ап

где п — произвольное положительное число. Уравнения (3) и (4) принимают вид

2 Ч„гг( г \

а2 8яв

а2

С д — + Л

А у

кс2

2

а = -4*С

Исследование системы уравнений (53) и (54) в общем виде довольно сложно, однако в ряде случаев удается найти достаточно простые решения. Например, для случая п = 4 (излучение) плоского пространства к = 0 и положительной вакуумной энергии затравочное решение имеет вид

а =

«Потенциал» в такой модели равен

Гр+ц. 1=м лГ 1 -

3 Г с2 У 3 I 8Ь2 в

то есть при Ь > ^ = I—3— агсэт 1 нарушается сильное энергетическое 32жСЛ

условие. Общее решение для такого потенциала:

Г С )1/4 с&е —

А = [С| 8Ь1/2 в 1п-2,

^лу а

где, как и ранее, мы ввели параметр е, нулевое значение которого дает нам исходное решение и который во всех остальных случаях можно полагать равным единице. Нетрудно заметить, что в последнем случае имеется три класса решений. При а < 1 вселенная открыта и на больших временах расширяется по экспоненциальному закону. Разумеется, в такой модели есть и горизонт событий.

41

42

Случай а> 1 более интересен. Если строго а> 1, то в момент времени tf = J arcc^h а имеем конечную сингулярность. Любопытно отметить, что при a > 1 + 42 сильное энергетическое условие не нарушается (вселенная успевает сколлапсировать раньше момента времени tv).

Когда а = 1, получающееся решение можно назвать квазисингулярным, так как

lim А ~ limsh1/2 lncth— ~ exp(- (9/2).

t>>0 9>>0 2

Значит, сингулярность достигается лишь асимптотически.

И в сингулярной, и квазисингулярной моделях горизонты событий отсутствуют. Действительно, для квазисингулярной вселенной путь светового луча вдоль радиальной геодезической

г d9

J 9 '

о sh1/2 9lncth— 2

Этот интеграл расходится, поскольку на верхнем пределе подынтегральное выражение имеет экспоненциальное асимптотическое поведение. Аналогично доказывается отсутствие горизонтов и для сингулярной модели.

Наконец, вычислим уравнение состояния:

с&

w =---+-

1 2 (1 - sh2 ()ln2-2

3 3 ' cth( ' ch9ln-2 - 2

Анализируя это соотношение, получаем, что и для открытой, и для квазисингулярной моделей вселенная стремится к режиму дБ (т ^ - 1) с ростом t. Финальная сингулярность соответствует т = - 1/3, а начальная (для всех моделей) — радиационно-доминированной вселенной с т = 1/3.

Замечание 3. Разумеется, все приведенные решения являются только модельными хотя бы потому, что они не содержат инфляционной стадии (ни быстрой, ни медленной). Однако можно использовать метод, описанный в [1], вводя (разумеется, вручную) слабую температурную зависимость параметров. Например, если допустить, что а = а(Т) в (56) (причем модуль а возрастает при убывании температуры (Т ^ 0)), то в будущем вселенная, находящаяся сейчас на стадии ускоренного расширения (а < 1), может сколлапсировать (если параметр вырастет так, что будет выполняться неравенство а > 1, то коллапс становится неизбежным) в финальную сингулярность, причем горизонт событий будет отсутствовать.

r

а

Замечание 4. Использование условия согласованности в будущем выглядит очень необычно. Традиционный физический подход подразумевает, что какие-либо предсказания делаются с использованием уравнений движения и начальных условий. Условие же «согласованности», анонсированное в данной работе, кажется парадоксальным, поскольку основано главным образом на различных предположениях о будущей эволюции вселенной. С нашей точки зрения мы имеем дело не с логически противоречивой, а просто с непривычной ситуацией. Более того, такого рода обстоятельства не являются чем-то совершенно новым в физике.

Скажем, знаменитая проблема «убегающих решений» в классической электродинамике, связанная с учетом сил радиационного трения, решается наложением дополнительных граничных условий в будущем. Еще более ярким примером является физика черных дыр: поскольку горизонт событий черной дыры определяется всем веществом, падающим на черную дыру на протяжении ее существования, динамика горизонта определяется путем задания граничных условий не в прошлом, а в будущем [24]. Поэтому подход, основанный на применении «принципа согласованности в будущем», выглядит даже менее странно, чем использование «телеологических» граничных условий в физике черных дыр. В конце концов, вся остальная физика разворачивается на арене пространства-времени не какой-то абстрактной, а совершенно реальной вселенной, и не кажется таким уж странным, что требование согласованности накладывает условие не только на локальные физические законы, но и на геометрию пространства-времени этой вселенной.

Если серьезно отнестись к «принципу согласованности в будущем», то, как мы видели, его использование накладывает достаточно жесткие условия на динамику вселенной (см., например, условия (г), (и) из раздела 3.2). Это означает, что такой принцип вовсе не носит какого-то «общефилософского» значения, а является чисто физическим и проверяемым. Резонно предположить, что последовательное и систематическое изучение данного принципа позволит получить еще более строгие ограничения на то, какой должна быть вселенная. Хорошо известно, что очень многие параметры наблюдаемой вселенной оказываются тщательно «подогнанными», и единственным существующим объяснением этого удивительного обстоятельства является антропныи принцип. Антропныи принцип имеет смысл только при наличии космологического мультиверса [25], но даже при этом он не объясняет всех совпадений. Не исключено, что использование принципа «согласованности в будущем» позволит усилить антропные предсказания.

4. Преобразование Дарбу

Выше мы использовали уравнение (4) в качестве генератора точных «согласованных» космологических решений. Разумеется, уравнение Шредингера можно решить лишь для очень специальных видов потенциала и(Ь). Эффективным способом генерации новых интегрируемых

43

44

моделей является метод преобразования Дарбу (ПД) [26; 27], который уже использовался одним из авторов этой статьи в [3] для построения новых классов точных решений уравнений Эйнштейна — Фридмана в рамках метода линеаризации, использованного в [1]. Ключевым свойством, делающим ПД полезным в космологическом контексте, является то, что эти преобразования переводят функции Йоста исходного уравнения Шрёдингера (6) с заданным потенциалом и в функции Йоста уравнения с преобразованным потенциалом и(1) [22]. Другими словами, ПД сохраняют асимптотическое поведение исходных решений при условии, что опорные функции, входящие в крамовские определители (см. [23]), имеют перемежающиеся нули. С точки зрения построения точных решений, описывающих «согласованные» космологические модели, это свойство можно переформулировать следующим образом: корректно выбирая опорные функции, с помощью которых реализуется ПД, можно получать бесконечные последовательности таких потенциалов №), что решения уравнения (4) будут описывать согласованные в будущем космологические модели, то есть модели, удовлетворяющие условиям (г), (гг) из раздела 3.2. Мы не будем приводить точное математическое доказательство этого утверждения, а ограничимся демонстрацией одного конкретного примера.

Пусть т = - 1/3. Легко убедиться, что уравнение (6) выполняется при Я = и = 0, а решение имеет вид

/а , 3(а02 + с21 с2К + с21 ,сгг.

а

8жСа^ 8жGa0í

Для реализации ПД необходимо второе, отличное от (57) решение щ (6) с тем же потенциалом и = 0, но другим значением спектрального параметра Я. Будем считать, что значение спектрального параметра отрицательно: Я = -к2, где к — вещественное число. Регулярное решение уравнения (6) очевидно: щ = сЬ К. Однократное ПД определяется формулами

а ^ а« = ащ1- ащ1 ,

Щ ¿2 (58) и = 0 ^ и(1) = и - 2 —- 1пщ1, Я^Я.

дЛ

Следовательно,

аЩ) = а0

с х = К. Несложно убедиться, что

(Х) 3 а2к(сЬ х яЬ х) + с2 сЬ4 х > 0 8жСа° сЬ2 х(сЬ х - х яЬ х)2

и

р<^+ = К >0.

х

Решение (59) описывает вселенную, которая рождается при хг — 1,1997 и коллапсирует в финальную сингулярность при х/ = - хг. Следовательно, в окрестности точки х = х/

а«1) ~ аок2

при Ь ^ /. То есть построенная нами простейшая модель подобно исходной не содержит горизонтов событий в будущем, как мы и говорили выше. Кроме того,

а*1) ~ ао^ - Ь)

при Ь ^ Ь, поэтому для нашего решения не выполняются условия применимости знаменитой теоремы Борде — Гута — Виленкина (см. [6]), поскольку параметр Хаббла

к

Н = -

2сЬ2 х

и Н ^ - ж при Ь ^ íl.

Таким образом, можно заключить, что использование ПД может оказаться чрезвычайно эффективным методом генерации новых классов точных космологических решений, демонстрирующих свойство «согласованности».

45

Заключение

В данной работе мы предложили простой метод генерации точных решений уравнений Фридмана, при котором задача сводится к решению уравнения Шрёдингера для функции, степенным образом зависящей от масштабного фактора ап, и с «потенциалом», пропорциональным п2р - 3п(р + р/с2)/2.

Показано, что данный подход позволяет получать широкий набор интересных решений, включая решения, содержие инфляционную фазу и естественный (без подгона параметров) выход из нее.

Численно проанализирован пример интегрируемого потенциала в виде гармонического осциллятора. Продемонстрировано, что такая модель может быть определенным образом согласована с данными астрономических наблюдений.

Проведен достаточно подробный анализ двухпараметрических семейств решений в пространствах с нулевой (см [1; 3]) кривизной к = 0, но для к = + 1 при п = 1.

Выявлено, что такие двухпараметрические решения в общем случае не содержат будущих горизонтов событий, что позволяет говорить о «согласованности» данных моделей. По мере приближения к финальной сингулярности параметр уравнения состояния ю ^ - 1/3, а потенциал самодействующего скалярного У(ф) ^ е-^, что характерно для струнных моделей.

Доказано, что использование преобразований Дарбу позволяет строить богатые семейства точных решений, описывающих «согласованные» космологические модели.

Нет сомнений, что описанный метод может оказаться весьма плодотворным не только в изучении режимов инфляции или при построении моделей пространства-времени без горизонтов, но и в решении других космологических задач. В качестве интересного примера приведем следующий результат, недавно вызвавший активные дебаты и широкий резонанс в научном сообществе. В работе [28] (см. также [29—34]) изучался вопрос о существовании метастабильного dS-вакуума в рамках теории струн. Авторы предположили, что из многообразия вакуу-мов (логарифм числа которых по оценке «сверху» достигает трех (!) порядков) лишь немногие могут быть физически реализуемы. Физические потенциалы скалярных полей, по мнению авторов этого исследования, должны удовлетворять следующему критерию: модуль логарифмического градиента от потенциала должен быть ограничен «снизу» положительной константой. Этот критерий (если он верен) делает практически невозможным реализацию деситтеровского вакуума в рамках фундаментальной теории и однозначно указывает на то, что физическая темная энергия должна быть квинтэссенцией. Этот вывод вполне согласуется с принципом согласованности в будущем. Замечательно и то, что простые лиувиллевские потенциалы вроде тех, что были построены выше, как раз удовлетворяют этому критерию. Разумеется, все это не может служить непосредственным аргументом в пользу принципа «согласованности», но важно отметить, что и результаты работы [28], и использование этого принципа приводят к одинаковым и неожиданным предсказаниям касаемо природы темной энергии — вероятно, главной загадки современной фундаментальной физики.

Список литературы

1. Журавлев В. М., Червон С. В., Щиголев В. К. Новые классы точных решений в инфляционной космологии // Журнал экспериментальной и теоретической физики. 1998. Т. 114, вып. 2. С. 406—417.

2. Yurov A. V. Phantom Scalar Fields Result in Inflation Rather Than Big Rip // Eur. Phys. J. Plus. 2011. Vol. 126, iss. 12. 132.

3. Верещагин С. Д., Юров А. В. Преобразование Дарбу и точно решаемые космологические модели // Теоретическая и математическая физика. 2004. Т. 139, № 3. С. 405—422.

4. Yurov A. V., Yurov V.A. Nonsingular Brane Solutions via the Darboux Transformation // Phys. Rev. D. 2005. Vol. 72, iss. 2. 026003.

5. Andrianov A.A., Cannata F., Kamenshchik A.Y. Smooth Dynamical Crossing of the Phantom Divide Line of a Scalar Field in Simple Cosmological Models // Phys. Rev. D. 2005. Vol. 72, iss. 4. 043531.

6. Borde A., Guth A., Vilenkin A. Inflationary Space-Times are Incompletein Past Directions // Phys. Rev. Lett. 2003. Vol. 90, iss. 15. 151301.

7. Долгов А. Д., Зельдович Я. Б., Сажин М. В. Космология ранней вселенной. М., 1988.

8. Riess A. G., Filippenko A. V., Challis P. et al. Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant // Astron. J. 1998. Vol. 116, iss. 3. P. 1009—1038.

9. Perlmutter S., Aldering G., Goldhaber G. et al. Measurements of Q and Л from 42 High-Redshift Supernovae // Astron. J. 1999. Vol. 517, iss. 2. P. 565 — 586.

10. Чернин А. Д. Космический вакуум // Успехи физических наук. 2001. Т. 171, № 11. С. 1154-1175.

11. Banks T. Cosmological Breaking of Supersymmetry? arXiv:hep-th/0007146.

12. Witten E. Quantum Gravity in De Sitter Space. arXiv:hep-th/0106109.

13. He X.-G. Accelerating Universe and Event Horizon. arXiv:astro-ph/0105005.

14. Gonzalez-Diaz P. F. Fundamental Theories in a Phantom Universe. arXiv:hep-th/0411070.

15. Li M. A Model of Holographic Dark Energy // Phys. Lett. B. 2004. Vol. 603, iss. 1. P. 1-5.

16. Page D.N. The Lifetime of the Universe // J. Korean Phys. Soc. 2006. Vol. 49, iss. 9 (2). P. 711-714.

17. Goheer N., Kleban M., Sussind L. The Trouble with de Sitter Space / / JHEP. 2003. Vol. 7. 056.

18. Kachru S., Kallosh R., Linde A., Trivedi S. R. de Sitter Vacua in String Theory // Phys. Rev. D. 2003. Vol. 68, iss. 4. 046005.

19. Kachru S., Pearson J., Verlinde H. Brane / Flux Annihilation and the String Dual of a Nonsupersymmetric Field Theory // JHEP. 2002. Vol. 6. 021.

20. Yurov A. V., Yurov V. A. One More Observational Consequence of Many-Worlds Quantum Theory. arXiv:hep-th/0511238.

21. Tipler F.J. Genesis: How the Universe Began According to Standard Model Particle Physics. arXiv:astro-ph/0111520.

22. Matveev V.B., Salle M.A. Darboux Transformation and Solitons. Berlin ; Heidelberg. 1991.

23. Crum M.M. Associated Sturm-Liouville Systems // Quart. J. Math. 1955. Vol. 6, iss. 2. P. 121-127.

24. Новиков И. Д., Фролов В. П. Черные дыры во Вселенной // Успехи физических наук. 2001. Т. 171, № 3. С. 307-324.

25. Линде А. Д. Физика элементарных частиц и инфляционная космология. М., 1990.

26. Darboux G. Sur une proposition relative aux équations linéaires // C. R. Acad. Sci. 1882. Vol. 94. P. 1456-1459.

27. Darboux G. Theorié generale des surfaces. N.Y., 1972.

28. Obied G., Ooguri H., Spodyneiko L., Vafa C. De Sitter Space and the Swampland. arXiv:1806.08362.

29. Ooguri H., Palti E., Shiu G., Vafa C. Distance and de Sitter Conjectures on the Swampland // Phys. Lett. B. 2019. Vol. 788. P. 180-184.

30. Agrawal P., Obied G., Steinhardt P. J., Vafa C. On the Cosmological Implications of the String Swampland // Phys. Lett. B. 2018. Vol. 784. P. 271-276.

31. Denef F., Hebecker A., Wrase T. de Sitter Swampland Conjecture and the Higgs Potential // Phys. Rev. D. 2018. Vol. 98. 086004.

32. Akrami Y., Kallosh R., Linde A., Vardanyan V. The Landscape, the Swampland and the Era of Precision Cosmology // Fortsch. Phys. 2018. Vol. 66. 1800075.

33. Cicoli M., De Alwis S., Maharana A. et al. De Sitter vs Quintessence in String Theory // Fortsch. Phys. 2018. Vol. 66. 1800079.

34. David Marsh M. C. The Swampland, Quintessence and the Vacuum Energy // Phys. Lett. B. 2019. Vol. 789. P. 639-642.

Об авторах

47

Артём Валерьевич Асташёнок — д-р физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Россия. E-mail: artyom.art@gmail.com

48

Артём Валерианович Юров — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Россия. E-mail: artyom_yurov@mail.ru

Анна Валентиновна Япарова — канд. физ.-мат. наук, ассист., Балтийский федеральный университет им. И. Канта, Россия. E-mail: AYAparova@kantiana.ru

The authors

Prof. Artyom V. Astashenok, I. Kant Baltic Federal University, Russia. E-mail: artyom.art@gmail.com

Prof. Artyom V. Yurov, I. Kant Baltic Federal University, Russia. E-mail: artyom_yurov@mail.ru

Dr Anna V. Yaparova, Assistant, I. Kant Baltic Federal University, Russia. E-mail: AYAparova@kantiana.ru