Гравитация. От хрустальных сфер до кротовых нор

Все это однозначно определяет эволюцию масштабного фактора a(t), плотность и давление p. Становится известной геометрия космологического пространства-времени, а значит, тип и динамика 3-мерного пространства, в котором звезды, галактики, скопления галактик ведут себя в соответствии с современными наблюдениями. То есть гравитационного парадокса не возникает.

Теперь зададимся, возможно, провокационным вопросом. А можно ли описать расширяющуюся Вселенную с помощью гравитации Ньютона? Оказывается, можно! Обратимся к опыту Зельдовича. В своих лекциях он всегда старался представить материал простейшим способом, а рассказывая о космологических решениях, по возможности ограничивался теорией Ньютона.

Рис. 9.4. Схема расчета ускорений

Рассмотрим шар радиуса R₀ и элемент массы m внутри шара на расстоянии от центра R R₀ (рис. 9.4). Такая задача рассматривалась еще Ньютоном. Он же и установил, что сила, действующая на m, определяется массой материи внутри сферы радиуса R, а гравитационное действие внешних областей взаимно компенсируется:

F = GMm/R ²,

здесь: M = (4/3)R ³ – масса материи внутри сферы, – плотность материи, распределенная однородно. Следовательно, ускорение

a = (4/3)GR = CR.

Таким образом, ускорение элемента в точке R пропорционально его расстоянию от центра и не зависит от радиуса шара R₀ при любом, сколь угодно большом R₀. Это фактически означает, что соотношение справедливо для бесконечной однородной вселенной.

Но остается вопрос: решение было найдено для некоторого центра, в котором ускорение равно нулю, а в других точках имеет вполне определенную величину и направлено к центру. А где такой центр в бесконечной однородной вселенной? На самом деле, никакого выделенного центра нет или, если угодно, таким центром может быть любая точка. Возьмем произвольную точку О', находящуюся на некотором расстоянии R_О' от «нашего» центра, ускорение которой a_О' = – СR_О'. Радиус-вектор и ускорение направлены, разумеется, в разные стороны. Перейдем в систему координат с центром в точке О' (рис. 9.5). Величины в этой новой системе координат будем обозначать штрихом.

Рис. 9.5. Переход к другой системе координат

Ускорения в старой и новой системах координат связаны правилом Галилея, которое, если кто забыл, справедливо не только для скоростей, но и для ускорений:

a'_A = a_A – a_О'.

Подставляя в это соотношение выражение для ускорений a_A = – CR_A и a_О' = – СR_О' и используя правило сложения векторов R '_A = R_A – R_O' (рис. 9.5), получим:

a'_A = – CR'_A.

Следовательно, наблюдатель в точке О' будет видеть ту же картину – все частицы материи имеют ускорение, направленное к нему. Ситуация несколько непривычная – ускорение направлено к центру, но центр «виртуальный», им всегда является точка, в которой находится наблюдатель. Такая ситуация концептуально отличается от ньтоновой, в которой предполагается наличие выделенного пространства, общего для всех наблюдателей.

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

Рассмотрим ситуацию, когда есть некоторые начальные скорости, направленные от наблюдателя (от «центра»). Для сохранения однородности в постановке задачи необходимо, чтобы начальная скорость была пропорциональна расстоянию от наблюдателя:

V = HR,

здесь H – коэффициент пропорциональности.

Вселенная будет расширяться, но скорость расширения будет падать. Из-за расширения будет уменьшаться плотность, а, следовательно, и ускорение. Что «пересилит»? Если начальная плотность достаточно велика, или, если угодно, мала начальная скорость, расширение через некоторое время сменится сжатием. При достаточно большой начальной скорости расширение будет продолжаться вечно. Качественно ситуация аналогична, например, рассмотрению стартовавшей с Земли ракеты. При скорости, большей второй космической, ракета может преодолеть притяжение и улететь на бесконечность.

В нашем случае также можно определить критическое распределение скоростей, в данном случае это параметр H_к, при превышении которого сжатие никогда не наступит. Его значение определяется соотношением:

3H_к² = 8G.

Но точно так же, можно оперировать с критической величиной плотности, рассчитывая ее по отношению к параметру H. Именно так делается при анализе решений Фридмана. Мало того, это соотношение для определения критической плотности полностью совпадает с фридмановским, см. Дополнение 8.