Гравитационные поля, возникающие при стохастическом самоподобном распределении материи во Вселенной Текст научной статьи по специальности «Физика»

2015

ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА

Сер. 4. Том 2 (60). Вып. 1

ФИЗИКА

УДК 53.01+51-71

А. И. Шерстюк, О. М. Бритаева

ГРАВИТАЦИОННЫЕ ПОЛЯ, ВОЗНИКАЮЩИЕ ПРИ СТОХАСТИЧЕСКОМ САМОПОДОБНОМ РАСПРЕДЕЛЕНИИ МАТЕРИИ ВО ВСЕЛЕННОЙ

Санкт-Петербургский государственный экономический университет, Российская Федерация, 191023, Санкт-Петербург, ул. Садовая, 21

В статье анализируются возможности использования мультифрактальной модели для описания гравитационных полей, возникающих при стохастическом самоподобном распределении материи во Вселенной. Модель фрактального распределения принципиально отличается от модели однородного распределения вещества в стандартной космологии и учитывает тонкие особенности фрактальных распределений, связанных с кластеризацией масс, что имеет принципиальное значение при вычислении распределения гравитационного поля при наличии флуктуаций плотности материи в пространстве. В случае обычно предполагаемого изотропного и однородного характера крупномасштабного распределения масс в пространстве напряжённость создаваемого ими поля в заданной точке в среднем равна нулю. Отличие от нуля значения напряжённости целиком обусловлено флуктуациями распределения масс. Предложенная бифрактальная модель учитывает фрактальную кластеризацию распределения материи в пространстве, а также то, что флуктуации на больших расстояниях вносят заметный вклад в значение напряжённости поля в заданной точке. Использование математического аппарата теории фракталов открывает новые возможности исследования самогравитирующих систем. Библиогр. 15 назв. Ил. 1.

Ключевые слова: теория фракталов, самоподобное распределение материи, гравитационные поля.

A. I. Sherstyuk, O. M. Britaeva

THE GRAVITATIONAL FIELDS ARISING FROM A STOCHASTIC SELF-SIMILAR DISTRIBUTION OF MATTER IN THE UNIVERSE

Saint-Petersburg State University of Economics, 21, ul. Sadovaya, St. Petersburg, 191023, Russian Federation

The article analyzes the possibility of using multi-fractal models for the description of the gravitational field arising from a stochastic self-similar distribution of matter in the Universe. The fractal distribution model is fundamentally different from the model of homogeneous distribution of matter in the standard cosmology. This model takes into account the subtle features of fractal distributions associated with mass clustering, which is of fundamental importance when calculating the distribution of the gravitational field in the presence of fluctuations of the density of matter in space. In the event of the typically assumed isotropic and homogeneous nature of

large-scale mass distribution in space, the intensity of the field they create is, at a given point, on average equal to zero. Non-zero intensity is entirely due to fluctuations in the distribution of masses. Proposed bifractal model takes into account the fractal clustering of the distribution of matter in space, and the fact that the fluctuations at large distances make a significant contribution to the value of the field strength at a given point. The use of mathematical apparatus of the theory of fractals opens new possibilities in the study of self-gravitating systems. Refs 15. Figs 1.

Keywords: fractal theory, self-similar distribution of matter, gravitational fields.

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

Данные глубоких обзоров распределения галактик с определением красных смещений [2, 3] и с учётом наличия тёмной материи, существованием которой ранее пренебрегали, указывают на наличие крупномасштабных неоднородностей вплоть до расстояний в сотни мегапарсек. Установлено, в частности, что парная корреляционная функция распределения галактик и их скоплений обладает свойством масштабной инвариантности до расстояний порядка 800 Мпк независимо от направления в пространстве. Вытекающее из космологического принципа [4] свойство локальной изотропии при условии неоднородности распределения масс предполагает нарушение требования аналитичности (регулярности) распределения материи в пространстве. Хорошим приближением для описания подобного распределения является модель изотропного статистического фрактала [5, 6] в области расстояний от Ri « 1 Мпк до верхней границы R2 ~ 800 Мпк. Нижняя граница определяется размером типичной галактики. До расстояний с параметрами R > R2 достоверные данные в настоящее время отсутствуют.

Поскольку модель фрактального распределения [7] принципиально отличается от однородного распределения вещества в стандартной космологии, с точки зрения изучения звёздной динамики представляет интерес, в частности, оценка флуктуаций напряжённости гравитационного поля F от всех источников в пространстве. Для однородного распределения вещества эта задача была решена в классических работах Чанд-расекара и фон Неймана [8, 9], где показано, что плотность распределения поля W(F) удовлетворяет соотношению Хольцмарка [10]. Пример применения этих соотношений к исследованию полного уширения спектральных линий, обусловленного совокупным воздействием различных механизмов, содержится в [11].

Случай степенного характера убывания плотности вещества

р - r-p ,p> 0 (1)

рассмотрен в работах [12, 13]. Зависимость (1) воспроизводит, в частности, усреднённое распределение вещества, характерное для монофрактальной модели с размерностью D = 3 — p («миксированный фрактал»). Однако такая модель не учитывает более тонкие особенности фрактальных распределений, связанных с кластеризацией масс, что имеет принципиальное значение при вычислении распределения гравитационного поля при наличии флуктуаций плотности материи в пространстве. В общем случае трёхмерное распределение масс может быть описано в рамках мультифрактальной модели

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

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

Согласно принципу самоподобия, число источников поля в сферическом объёме радиуса Я удовлетворяет соотношению

N (Я) = БЕВя, (2)

где П^ — фрактальная размерность по числу источников поля; Б — постоянный множитель, определяющий морфологию распределения и зависящий от внутреннего радиуса обрезания. Если Я0 — радиус образующего элемента (нижняя граница зоны самоподобия), то Б = N(Ео)/Е^м. Согласно [2], « 1,23.

Учёт кластеризации объектов предполагает введение аналогичного соотношения для радиального распределения суммарной массы в объёме радиуса Я:

М(Н) = АН*>", А=(3) Я0

где Бм — размерность по отношению к распределению масс.

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

В случае обычно постулируемого изотропного и однородного характера крупномасштабного распределения масс в пространстве напряжённость создаваемого ими поля в заданной точке в среднем равна нулю. Отличие от нуля значения напряжённости целиком обусловлено флуктуациями распределения масс. Величина таких флуктуаций определяется наличием локальных неоднородностей распределения. В предположении однородного в среднем крупномасштабного распределения источников вклад удалённых масс быстро зануляется с увеличением расстояния, что связано прежде всего со взаимной компенсацией влияния многочисленных флуктуаций. Зачастую достаточно бывает учесть доминирующий вклад только одного ближайшего источника (звезды или галактики) [10, 11].

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

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

п(г, т) = 4пп(г, т) = по(г)/Т (т),

где по (г) — плотность числа объектов на поверхности радиуса г; /т(т) — нормированная на единицу функция распределения масс на расстоянии г. Поскольку характер распределения поля в точке наблюдения несущественно зависит от конкретного вида функции /т (т), можно положить

/т(т) = Ь(т — тср.(г)),

(4)

где тср.(г) — средняя масса объектов на расстоянии г от точки наблюдения (г = 0), т. е. будем считать, что массы детерминированным образом зависят от расстояния, причём средняя масса источника тяготения на расстоянии г

тср.(г) = J /(г,т)тйт.

Число объектов в сфере радиуса Я

Г *

N (Я) =4п п0(г)г2 ¿г, (5)

■! 0

а радиальное распределение масс: р(г) = п0(г)тср (г). При этом

гя

0

Из выражений (2) и (5) вытекает:

, я

М (Я) = 4п р(г)г2 ¿г.

0

п0(г)

Р(г)

1 ¿Ж

4пг2 ¿Я 1 ¿Ж

БВ

Я=т

4п

АРМ

4ж

N г-(3-£>лг)

-(3-Бм)

4пг2 ¿Я я=т

Таким образом, эффективная средняя масса источника тяготения

Р(г)

тср (г)

по(г)

тога,

(6)

где то = АВм/БВN, а = Вм — ВN.

Далее будем считать, что фрактальное распределение имеет место во всей области изменения г(Я\ ^ 0,Я2 ^ те). В «ближайшей» области (в пределах галактики) имеет место почти однородное распределение объектов. Вкладом области г > Я2 в распределение поля в точке наблюдения (г = 0) можно пренебречь из-за быстрого убывания напряжённости с увеличением расстояния.

С учётом выражения (6) в рамках ньютоновской теории тяготения напряжённость гравитационного поля в точке с координатами г

т = а)

где О — гравитационная постоянная.

Плотность распределения вероятностей значений гравитационного поля каждой из N частиц в объёме радиуса Е будет определяться формулой (7). Согласно общей теории [14], плотность распределения вероятностей совокупности N частиц даётся обобщённой свёрткой плотностей распределения вероятностей каждой из частиц. Следовательно, характеристическая функция Ф(к) суммы независимых случайных величин равна произведению характеристических функций отдельных величин. Учитывая, что плотность распределения вероятности каждой из величин равна по(г),

Ф(к)

1 ]N

— exp(ikF(r))no(r)dr

(8)

Таким образом, для плотности распределения поля W(F), создаваемого N источника-

W(F) = Т—[ Ф(к)ехр(-1кР)0к. (9)

(2п)3 J

Подставляя (7) в (8) и переходя к пределу R ^ ro(N(R) ^ ж), с учётом сферической симметрии системы после некоторых преобразований получаем

Ф(к) = exp(-ако), (10)

где к = |к|,

Dn

о =-,

2 -Dm + DN'

f ОО

a = (Gm0)aoB (z - sin z)z-2-adz. (11)

J 0

Выражение (11) имеет смысл только при о < 2,т. е. DM < D* = DN/2+2, где D* — критическое значение фрактальной размерности по массам. Подставляя (10) в (9), получаем окончательное выражение для распределения напряжённости гравитационного поля о

W(F) = 9 о з J^ xsinxexp |— ж0 jete, (12)

где в силу сферической симметрии W(F) = W(F), F = |F|, F0 = а1/о — «стандартная» величина напряжённости поля.

Далее, следуя Хольцмарку [10], удобно ввести стандартную нормированную на единицу функцию распределения Но(в) = 4nF2F0W(F), зависящую от безразмерного параметра в = F/F0,

На($) = ъ [ ехР | — (тт) Ixsinxífa;. (13)

ПРJ0 { \Р/ J

При однородном распределении (DN = DM = 3) получаем: о = 3/2, а функция Но(в) совпадает с функцией Хольцмарка: Н3/2(в) = Н(в).

Следует отметить, что в ряде случаев интеграл в формуле (13) может быть выражен в элементарных функциях. При Вм = 2 будем иметь о = 1 независимо от значения DN. В этом случае формула (13) даёт

Нг(в)

4в2

п(в2 + 1)2

(14)

При о = 1 Г0 -При о = 2

пОт0Б/4 = пАО/2Вм, т. е. в = 2ВМЯ/пАС.

Н2(в)

2у/л

е — в2/4.

(15)

При о = 2п , где п — целое число, асимптотика Но(в) при больших значениях в имеет вид

Но(в)

2Г(о

2) пв

^ , .— эш-.

2

При целых чётных о функция Но убывает экспоненциально с ростом в.

При а = 3/2 распределение (13) совпадает с распределением Хольцмарка, где При о > 3/2 максимум функции Но(в) сдвигается в сторону больших в по сравнению с классическим распределением Хольцмарка, т. е. имеет место относительное увеличение вероятности сильных полей («гиперхольцмарковское» распределение) за счёт вклада далёких галактик и их скоплений. При о < 3/2 происходит относительное увеличение вероятности слабых полей («субхольцмаровское» распределение) вследствие более быстрого убывания плотности распределения масс с расстоянием.

В большинстве работ, посвящённых анализу наблюдательных данных по крупномасштабному распределению галактик, приводятся значения ВN = 1,2 + 1,4. Значение Вм может быть получено из анализа светимости галактик в видимом и инфракрасном диапазонах с использованием известного соотношения между массой и светимостью. Согласно последним оценкам, Вм « 2,0 [5, 6]. Однако это значение требует дальнейшего уточнения с учётом значительного вклада тёмной материи, масса которой, по последним данным, составляет значительную часть массы галактик.

На рисунке представлена плотность распределения модуля гравитационного поля Но(в) (13) при ВN = 1,3 и различных значениях параметра о.

Рис. 1. Стандартная функция распределения безразмерной напряжённости поля

0,40,35 0,3 0,25 0,2 0,15 0,1 0,05

0

0,40,81,21,6 2 2,42,83,23,6 4 4,44,!

в

5,25,6 6 6,4 6,8

а

2

С учётом значений DM: если DM = 2,43 — распределение Хольцмарка (жирная линия); если DM = 2 — формула (14); если DM = D* — формула (15); если DM < 2,43 — субхольцмарковское распределение; если 2,433 < DM < 2,65 — гиперхольц-марковское распределение. Суб- и гиперхольцмарковское распределение представлены со стандартными характеристиками.

Следует отметить, что полученные распределения устанавливаются при значительно больших «радиусах обзора» R, чем в случае однородного в среднем распределения [15] источников поля.

Литература

1. Peebles P. J.E. The large scale structure of the Universe. Princeton: Princeton Univ. Press., 1990.

2. Perdang J. Self-gravitational fractal configuranions // Vistas in Astron. 1990. Vol. 33. P. 171—191.

3. Барышев Ю. В., Теерикорпи П. Фрактальная структура Вселенной. Очерк развития космологии. Нижний Архыз: Изд-во САО РАН, 2005. 410 с.

4. Weinberg S. Gravitation and Cosmology. N.-Y., 1972.

5. Baryshev Yu. V. On the fractal nature of the large scale of the universe // Astron. Astrophys. Transactions. 1994. Vol. 5. P. 15-23.

6. Baryshev Yu. V., Sylos LabiniF., MontuoriM., Pietronero L. Facts and ideas in modern cosmology // Vistas in Astronomy. 1994. Vol. 38. P. 419-500.

7. Mandelbrot B. B. Fractals and multifractals. New York: Springer, 1988.

8. Chandrasekhar S., von Neumann S. // Astrophys. J. 1943. Vol. 97. P. 1.

9. Chandrasekhar S. // Rev. Mod. Phys. 1943. Vol. 15. P. 1.

10. Holtsmark J. // Ann. Phys. Ser. 4. 1919. Vol. 58. P. 365.

11. Каули Ч. Теория звёздных спектров. М.: Мир, 1974.

12. Шерстюк А. И. и др. // Тез. докл. 9-й Рос. гравитационной конф. Новгород, 1996. C. 131.

13. Sylos LabiniF. // Astrophys. J. 1994. Vol. 443. P. 464.

14. Золотарёв В. М. Одномерные устойчивые распределения. М.: Наука, 1983.

15. Del Popolo A. // Astron. and Astrophys. 1996. Vol. 305. P. 999.

Статья поступила в редакцию 20 августа 2014 г.

Контактная информация

Шерстюк Алексей Иванович — доктор физико-математических наук, профессор. Бритаева Ольга Муратовна — старший преподаватель; e-mail: olga-britaeva@yandex.ru

Sherstyuk Alexey Ivanovich — Doctor of Physics and Mathematics, Professor. Britaeva Olga Muratovna — Senior Lecturer; e-mail: olga-britaeva@yandex.ru