Элементы современной космологии и теории бран. Учебное пособие
.pdf
Постановка задачи
имела нулей между инфляционной фазой, когда вторая производная от масштабного фактора положительна, и собственно фридмановским режимом. В противном случае фаза инфляции и фаза эпохи доминирования пылевой материи окажутся разделенными точкой сингулярности, т.е. окажутся двумя причинно не связанными вселенными. Именно это свойство позволяет обеспечить выход из инфляции.
Во-вторых, ПД можно использовать для построения инфляционных моделей. Чтобы это показать, перепишем уравнение (1.2) в виде
φ&2 = V(1) −V ,
где V(1) — потенциал, построенный по обычным крамовским формулам из V.
Предположим, что эта функция имеет N нулей, тогда V(1) будет иметь по крайней мере N особенностей. Будем считать, что все эти особенности отсутствуют в исходном потенциале
V. Если в окрестности одной из этих особенностей V(1) не является положительно определенной функцией (или функция V имеет свою особенность, окрест которой она положительно определена), то существуют такие t1 и t2 , лежащие слева и
справа относительно особенности, что
φ&(t1 )= φ&(t2 )= 0, φ&2 < 0, t1 < t < t2 .
Доказательство этого утверждения мы оставляем читателю в качестве упражнения. Теперь ясно, что при выполнении указанных условий вселенная непременно войдет в инфляционную фазу при t, стремящемся к t1 . Действительно, предполо-
жим, что вселенная рождается в некоторый момент времени, предшествующий t1 , и допустим, что существует такой промежуток времени, в течение которого выполняются все энергетические условия. Можно показать, что при стремлении к t1
вначале нарушится сильное энергетическое условие при действующем слабом. В результате вселенная начнет раздуваться
31
Глава 1
вплоть до момента t1 . В момент времени t1 мы получаем в
точности фазу Деситтера, но лишь на один момент. Далее следует «запретная зона», которую следует удалить, сшив решения в точках с t1 и t = t2 . Заметим, что такая сшивка требует
более жестких условий, чем в аналогичных задачах нерелятивистской квантовой механики, основанных на уравнении Шредингера. Здесь необходимо потребовать совпадения не только волновых функций и их первых производных в точке сшивки, но и вторых производных. Последнее необходимо во избежание скачков в давлении и плотности. Общее решение уравнения (1.1) содержит две константы интегрирования, поэтому не очевидно, что этим трем условиям сшивки всегда можно удовлетворить. Тем не менее легко убедиться, что третье условие (совпадение вторых производных) автоматически выполняется при выполнении первых двух! Это счастливое обстоятельство является следствием того, что сшивка происходит в точке перегиба logψ . Цена, которую приходится
заплатить, — наличие скачка у третьей производной масштабного фактора — не выглядит чрезмерной, по крайней мере пока справедливо классическое (т. е. не квантовое) рассмотрение. Таким образом, контролируя нули у опорной функции ПД, можно строить точные космологии с инфляцией. Если при этом использовать в качестве «затравочных» решения с правильными асимптотиками (например, с динамикой по закону двух третей), то вселенная будет автоматически выходить из инфляции без каких-либо подгонок. Ниже мы проиллюстрируем, как это работает, на простых примерах.
Простейшая модель |
|
|
|
|
Закон a ≈ t2/ 3 приводит к потенциалу |
|
|||
V = |
2 |
, |
(1.3) |
|
3t2 |
||||
|
|
|
||
32
Простейшая модель
Общее решение уравнения (1.1) с потенциалом (1.3) и нулевой космологической постоянной имеет вид
ψ(−) = − |
c1 |
−c2 t2 , |
(1.4) |
|
t |
||||
|
|
|
где c1,2 — постоянные интегрирования. Решение описывает две не связанные друг с другом вселенные, (рис. 1).
Рис. 1. Зависимость масштабного фактора от времени для общего решения уравнений Эйнштейна с потенциалом (1.3)
Первая неограниченно продолжается в прошлое и имеет конечную сингулярность, а вторая выглядит достаточно странно: рождается из сингулярности, затем расширяется и переживает суперинфляцию, т. е. к моменту t = 0 масштабный фактор становится бесконечно большим. Затем происходит сжатие с переходом к конечному расширению, причем масштабный фактор меняет знак. Разумеется, это решение нельзя воспринимать буквально. Зона суперинфляции (слева и справа
33
Глава 1
от момента t = 0 ) является запретной, так как в ней нарушается слабое энергетическое условие. Легко убедиться, что давление и плотность имеют вид
p = − |
c1 |
(c1 + 4c2t3 ) |
, |
ρ = |
(2c2t3 −c1 )2 |
||
|
|
|
|
, |
|||
t2 |
(c1 +c2t3 )2 |
3t2 (c1 +c2t3 )2 |
|||||
а потому запрещенная зона, о которой шла речь выше, будет находиться в интервале от t1 , до t2 , где
t1 = |
4 −3 |
2 |
|
c1 |
≈ −0,12 |
c1 |
, |
t2 |
≈ 4,12 |
c1 |
. |
2 |
|
c2 |
c2 |
|
|||||||
|
|
|
|
|
|
|
c2 |
||||
Ясно, что решение справа должно определяться другими константами интегрирования:
ψ(+) = bt1 + b2 t2 .
Можно показать, что сшивка возможна при специальном выборе констант в этой формуле. Мы настоятельно советуем читателю в качестве упражнения вычислить эти константы. Вместе с тем такая сшивка оставляет неудовлетворенность, поскольку мы сшиваем два решения в разные моменты времени. Однако, несложно произвести сдвиг по времени так, чтобы выполнялись два условия:
1)начальная сингулярность соответствовала бы t = 0 ;
2)t1 = t2 , то есть отсутствует скачок по времени.
Опуская вычисления, приведем ответ:
1) при 0 < t < t1 динамика описывается формулами
|
|
|
c3 |
|
|
|
2 |
|
2 / 3 |
|
|
|
|
|
|||
ψ(−)(t)= − |
|
1 |
|
−c3 |
(t −c / c |
|
) |
, V(−) = |
|
|
|
|
|
|
; |
||
|
|
|
2 |
(t −c / c |
|
)2 |
|
||||||||||
|
|
t −c1 / c2 |
2 |
1 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
2 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
||
2) при t1 < t |
динамика описывается формулами |
|
|
|
|
|
|
|
|||||||||
ψ(+)(t)= |
c13 (3 − 2 2 ) |
+ c32 |
(3 + 2 2 )(t −g)2 , V(+) = |
2 / 3 |
|
|
, |
||||||||||
|
|
|
(t −g)2 |
|
|||||||||||||
|
|
|
t −g |
|
|
|
|
|
|
|
|
|
|||||
34
Простейшая модель
где
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
3 2 −4 |
3 |
|
c |
|
|
|
|||||
t1 = |
|
|
1 |
1− |
|
|
|
|
≈ 0,51 |
|
|
1 |
, |
|||||||
|
|
|
|
|
2 |
c |
|
|||||||||||||
|
|
|
c |
|
|
|
|
|
|
|
|
2 |
|
|
||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
3 2 − |
4 3 |
|
|
c |
|
|
|||||||
g = |
|
|
|
1 |
|
|
|
|
|
|
≈ 0,01 |
|
|
1 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
. |
||||||||||
|
c |
2 |
|
|
|
|
2 |
|
|
|
c |
2 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Здесь введены новые константы c1,2 |
(для удобства) и t2 = t1 . |
|||||||||||||||||||
Сшивка уже осуществлена: |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
ψ(−)(t1 )= ψ(+)(t2 ), |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
& |
(−) |
|
& |
(+) |
(t2 ). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ψ |
|
(t1 )= ψ |
|
|
|
|
|
|
|
|||
График a(t) для второй вселенной приведен на рисунке 2.
Рис. 2. Зависимость масштабного фактора от времени после сшивки
35
Глава 1
Отметим, что при t < t1 кривая описывается формулой
a(−) = ψ(−) (t) , а при t > t |
1 |
: a(+) = ψ(+) (t) |
( c = c |
2 |
=1 ). На ри- |
|
|
1 |
|
сунке 3 изображен график зависимости ускорения &a& как функции времени. Излом расположен в точке t = t1 . Как видно,
имеет место инфляция на конечном интервале (там, где выступает зубец над осью t). Разумеется, инфляция возникает не при всяких t. Можно показать, что для ее возникновения должно выполняться неравенство
c2 |
> |
2 |
|
(12 2 −16)1/ 3 (3 −2 2 ) |
≈ |
0,38 . |
||
9 |
c12 (10 −7 2 ) |
|
||||||
|
|
|
|
c12 |
||||
Рис. 3. Зависимость ускорения масштабного фактора от времени после сшивки
Выбор c1 = c2 =1 как раз этому условию удовлетворяет. Отметим, что при достаточно больших t потенциал ведет себя как
V ≈ |
2 |
e− 3φ, φ ≈ |
2 |
log t . |
|
3 |
3 |
||||
|
|
|
36
Простейшая модель
Таким образом, мы построили простейшую модель инфляции с выходом и асимптотическим переходом к режиму фридмановского расширения вселенной находящейся в фазе преобладания вещества. Разумеется, этот пример не описывает нашу вселенную поскольку тут нет фазы осцилляции. Нет и вторичного ускорения. Тем не менее ясно, что, применяя предложенную методику, можно без принципиальных затруднений строить модельные космологии с выходом без использования подгонок параметров. Также нет необходимости использовать гибридные модели. Напомним, что решение описывает еще одну, сжимающуюся, вселенную, которая не имеет начала, но имеет конец. Такие решения выглядят весьма необычно, но их существование неизбежно, поскольку уравнения общей теории относительности обратимы по времени.
Преобразование Дарбу
Ясно, что любое точное решение уравнения (1.1), согласованное с энергетическими условиями, порождает точное решение уравнений Эйнштейна в метрике ФЛРУ. Поскольку не все решения (1.1) будут порождать положительно определенное выражение после подстановки в правую часть (1.2), то точные решения уравнений Эйнштейна являются подмножеством точных решений уравнения (1.1). С другой стороны, фактически все точные решения уравнения Шредингера (1.1) с «физическими» потенциалами могут быть построены с помощью ПД (мы не рассматриваем потенциалы типа прямоугольной ямы, ограничиваясь непрерывными функциями.). В качестве примера можно привести гармонический осциллятор, многосолитонные и даже конечнозонные потенциалы. Мы будем использовать это как «экспериментальный факт», не пытаясь придать утверждениям типа «все потенциалы, для которых уравнение (1.1) может быть явно решено для всех значений спектрального параметра, строятся с помощью ПД и,
37
Глава 1
возможно, еще некоторого конечного набора симметрий из других потенциалов с тем же свойством» точный математический смысл. Тогда, неформально говоря, можно заключить, что все точные решения уравнений Эйнштейна в плоской вселенной ФЛРУ, заполненной вещественным скалярным полем, могут быть построены с помощью ПД. В качестве простейшего примера можно привести потенциал (1.3), который стро-
ится однократным ПД на фоне нулевого потенциала V(0) = 0 . Результат N-кратного ПД можно компактно записать с помощью формул Крама:
ψ |
(N) |
= |
N+1 |
, |
V |
(N) |
= V − |
2 d2 ln N |
, |
(1.5) |
|||||
|
N |
|
|
3 |
|
|
dt2 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dkψ |
|
|
|
|
|
|
|
|
|
|
|
||
N = det |
|
i |
, |
k = 0,K, N −1; |
i =1,KN |
|
|||||||||
dtk |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dkψ |
|
|
|
|
|
|
|
|
|
|
||
N+1 = det |
|
|
i |
, |
k = 0,K, N; i =1,KN +1 . |
|
|||||||||
dtk |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В этих выражениях, ψi — частные решения (1.1) при не-
котором общем исходном потенциале V и разных собственных значениях 3Λ . Теперь ясно, как можно использовать ПД для построения моделей с двумя ускорениями. Мы уже обсу-
ждали необходимость наличия в потенциале V(N) особенностей. Из формул Крама следует, что все особенности совпадают с нулями функций, входящих в крамовский определитель N , поэтому задача сводится к контролю за этими ну-
лями. Будем считать, что исходный потенциал определяется выражениям (1.3). Решение (1.4) отвечает нулевому космологическому члену. Решения уравнения (1.1) с положительным и отрицательным значениями Λ-члена имеют вид
38
Преобразование Дарбу
ψ |
|
= μ |
|
− |
1 |
|
kt |
+μ |
|
|
+ |
|
1 |
|
−kt |
|
|
|
|
1 |
|
e |
|
|
1 |
|
|
e |
|
, |
|
||||||
|
kt |
|
|
|
|
|
|
|||||||||||
|
+ |
1 |
|
|
|
|
|
2 |
|
|
kt |
|
|
|
(1.6) |
|||
|
|
|
sin kt |
|
|
|
|
|
|
cos kt |
||||||||
ψ− |
= ν1 |
|
|
|
|
|
|
|||||||||||
|
kt |
−cos kt +ν2 |
|
kt |
|
+sin kt , |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где μ1,2 и ν1,2 — произвольные константы.
Используя эти функции для разных значений k в формулах (1.5) и считая, что N + 1-кратно одетое решение задано выражением (1.4) для всех номеров N , можно построить множество решений, асимптотически выходящих на фридмановский режим двух третей. Ясно, что N-кратно одетые функции будут содержать множество особенностей, окрест которых нарушается слабое энергетическое условие и, значит, будет возникать режим ускорения. Как и в примере из предыдущего раздела, решения нужно оборвать в точках перегиба логарифма N-кратно одетой волновой функции, где реализуется точный режим Деситтера, после чего произвести сшивку в этих точках. При этом необязательно сшивать решения в двух «соседних» точках. Под соседними мы понимаем точки, между которыми лежит лишь одна особенность N-кратно одетой волновой функции. Как уже отмечалось, всегда можно сшить любые две такие точки, что дает известную свободу для построения нужного нам решения. Заметим, что с ростом N явный вид решений быстро становится чрезвычайно сложным. Приведем несколько примеров.
Всюду считается, что «нормировочные» константы равны единице. Однократное ПД с опорной функцией из (1.6) приводит к выражениям
ψ |
(1) |
= |
|
μ[(x −1)3 +k3 +1]ex −[(x +1)3 +k3 −1]e−x |
||||||||
|
|
|
|
|
|
|
|
|
|
, |
||
|
|
|
|
|
|
μ(x −1)ex +(x +1)e−x |
|
|||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
V |
(1) |
= |
2k |
2 μ2e2x +e−2x −4μx2 −2μ |
, |
|
|||
|
|
|
|
3 |
|
|
(μ(x −1)ex +(x +1)e−x )2 |
|
||||
39
Глава 1
где μ1 =1, μ2 = μ, x = kt . При μ < 0 график для ψ(1) (t) изображен на рисунке 4.
Рис. 4. Зависимость ψ(1) от времени при μ < 0
Здесь уже вырезаны области, в которых нарушается слабое энергетическое условие. Граничные точки (там, где кинетический член скалярного поля равен нулю) следует сшить. Тогда становится ясно, что такое решение описывает две вселенные. Первая сжимается (причем в определенный момент, предшествующий коллапсу) в финальную сингулярность, она переживает фазу обращенной по времени инфляции. Динамика второй вселенной фактически является обращенной по времени динамикой первой: вселенная рождается в сингулярности, некоторое время расширяется, затем переживает фазу раздувания, выходит из инфляции и асимптотически переходит к фридмановскому расширению, характерному для вселенной с преобладанием вещества. Новая ситуация возникает при μ > 0 . На рисунке 5 приведен
график ψ(1) (t) для μ = 2 .
40
