О природе «черных дыр», «темной материи» и динамике Вселенной
Григорян Самвел Самвелович (1930-2015), доктор физико-математических наук, профессор, академик РАН, заведующий Отделом механики природных процессов Института механики МГУ имени М.В. Ломоносова, заслуженный профессор МГУ, действительный член Российской академии естественных наук, Российской инженерной академии, Национальной академии наук Республики Армения, Инженерной академии Украины, Международной академии астронавтики, Международной академии информатизации
В статье сформулировано положение, согласно которому теоретические количественные модели, предназначенные для описания природных процессов, не должны содержать сосредоточенных особенностей с бесконечными значениями физических параметров. На этой основе автор утверждает о непригодности теории гравитации А. Эйнштейна для описания астрофизических и космологических явлений. Автором построена количественная модель нейтронных звёзд и «чёрных дыр», основанная на классической механике и теории гравитации И. Ньютона, на современной теории элементарных частиц и на представлении о силовом воздействии гравитационного поля на фотоны.
Ключевые слова: астрофизика, космология, «темная материя», галактики, нейтронные звёзды, нейтрон-но-кварковые звёзды, «чёрные дыры», теория гравитации И. Ньютона, решение Шварцшильда.
Введение
В 1916 году Альберт Эйнштейн опубликовал свою теорию гравитации - так называемую общую теорию относительности (ОТО), имевшую целью устранить небольшие количественные отклонения предсказаний теории гравитации Исаака Ньютона от фактов, регистрируемые в астрономических наблюдениях. Основной «физический» эффект, содержащийся в ОТО, состоит в том, что присутствие в пространстве материи (массы) приводит к искривлению континуума пространство-время, к изменению его метрики, причём это происходит не только в местах, где располагается материя, но и вне этих мест - в пустоте. Этот эффект значительно усиливается при приближении к областям, где имеется материя.
Одно из первых точных решений задач для математической модели ОТО - простейшей стационарной центрально-симметричной задачи о метрике пространства-времени вокруг сосредоточенной конечной массы, построенное Карлом Шварцшильдом в том же 1916 году, показало, что модель ОТО вместо того, чтобы устранить «малые» несоответствия теории Ньютона данным наблюдений, приводит к весьма существенным и неожиданным качественным и количественным эффектам вокруг «точечной» массы - порождает сложную математическую особенность («сингулярность») в точке и определяет «гравитационный радиус» Rg - расстояние от центра симметрии в задаче, на котором имеется другая особенность, в частности, невозможность выхода фотонов за пределы сферы с этим радиусом. Внутри этой сферы, кроме того, кардинально меняется метрика пространства- времени и возникают странные проблемы причинно-следственного характера. Построенные позднее различные другие точные решения уравнений ОТО («космологические» решения А.А. Фридмана1, обобщения решения Шварцшильда на случаи вращающейся массы (точечной!), наличия магнитного поля и др.) - все без исключения также содержат математические сингулярности. Впоследствии С. Хокингом2 была доказана теорема, устанавливающая, что эволюция решения математической задачи для модели ОТО при весьма общих и мало ограничительных предположениях об исходных данных приводит к возникновению в решении сингулярностей!
С моей точки зрения, этого всего было достаточно, чтобы сделать простой вывод о непригодности модели ОТО для описания гравитационных процессов в природе, ибо в математических построениях, претендующих на описание любых естественных явлений, не должно быть особенностей (сингулярностей) с бесконечными значениями физических величин. Однако такой вывод почему-то в науке сделан не был. Более того, решения А.А. Фрид-
1 Friedmann A.A. "Über die Krümmung des Raumes." Zeitschrift für Physik 10.1 (1922): 377-386; Idem. "Über die Möglichkeit einer Welt mit konstanter negativer Krümmung des Raumes." Zeitschrift für Physik 21.1 (1924): 326-332.
2 Hawking S.W., Penrose R. "The Singularities of Gravitational Collapse and Cosmology." Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 314.1519 (1970): 529-548; Idem. "Black Holes in General Relativity." Communications in Mathematical Physics 25.2 (1972): 152-166; Idem. A Brief History of Time from the Big Bang to Black Holes. New York: Bantam Dell Publishing Group, 1988. 256 p.
мана, содержащие «вселенские» сингулярности, с привлечением некоторых наблюдательных фактов (закона «раз-бегания» галактик Э. Хаббла, факта существования «реликтового» излучения) были «канонизированы» и стали основой современной космологии...
Важно отметить ещё одну теоретическую работу - американских физиков Р. Оппенгеймера и Г. Волкова1, в которой изучался аналог задачи Шварцшильда, но для массы, распределённой в конечном объёме. В этой работе при использовании уравнений ОТО и конкретного уравнения состояния вида p = f (р) (p и р - давление и плотность) для распределённой массы было установлено, что равновесное решение задачи без особенностей существует только для значений суммарной массы m, не превышающих некоторую конечную величину т.. Удивительный результат, всецело определяемый сильнейшей нелинейностью уравнений ОТО! Можно показать, что этот результат не зависит от деталей уравнения состояния (то есть от вида функции f, при определённых естественных ограничениях на поведение f при больших p). Вывод Оппенгеймер и Волков сделали такой: поскольку при т > т. равновесной конфигурации без сингулярностей не существует, то, значит, в реальности происходит «коллапс» с собиранием всей массы в точку, то есть осуществляется решение Шварцшильда. И даже была выполнена работа2, в которой количественно описывалась динамика такого коллапса. Правда, при этом было принято, что давлением материи можно пренебречь, что, конечно, недопустимо, поэтому эта работа в принципе неудовлетворительна, даже в рамках ОТО. Курьёзность ситуации ещё и в том, что при т < т. одна и та же задача имеет два решения, и существенно разных: решение Оппенгеймера и Волкова без сингулярностей и решение Шварцшильда, имеющее сингулярность. Какое из них предпочесть? Такой вопрос ни у кого из лиц, занимавшихся этой проблематикой, даже не возникал.
Начиная с 1971 г., когда в космос была выведена аппаратура для регистрации источников рентгеновского излучения «на небе» (американский спутник «Ухуру»), было обнаружено множество (уже много сотен в наше время) таких источников - в составе двойных звёздных систем с одной видимой (оптической) и одной невидимой компонентами, а также в активных ядрах галактик. Анализ рентгеновских сигналов от таких объектов, возникающих, по существующим представлениям, при ускоряющемся натекании (аккреции) массы на гравитирующий невидимый объект из его окрестностей, и теоретические рассуждения о возможной их физической природе привели к заключению, что невидимые компоненты двойных систем с массой, меньшей трёх солнечных, - это нейтронные звёзды (компактные образования из «спрессованных» гравитацией нейтронов, формирующиеся на последних стадиях эволюции достаточно массивных светящих звёзд), а объекты с большими массами, включая и тёмные объекты с массами порядка 106-1010 солнечных масс в активных ядрах галактик, - это сингулярные объекты ОТО, так называемые «чёрные дыры» (наименование связано с тем, что, как выше отмечалось, свет из центральной области решения Шварцшильда не может «вырваться» наружу). В качестве доказательства приводятся два факта: 1) по особым характеристикам рентгеновского излучения (его периодичности и т.п.) делается вывод, что маломассивные объекты имеют твёрдую внешнюю поверхность, на которую натекает газ, «высасываемый» гравитацией этого объекта из оптической компоненты, поэтому делается заключение, что это - нейтронные звёзды; 2) отсутствие таких характеристик у рентгеновского излучения - это свидетельство отсутствия у тёмного объекта твёрдой поверхности, и, значит, этот объект есть сколлапсировавшая материя, собранная в точку, это - сингулярность ОТО. Разделяющая эти две разновидности объектов величина массы по порядку величины согласуется с результатом работы Оппенгеймера и Волкова (уточнённым позднее), целиком обязанного своим существованием модели ОТО, и из этого делается вывод о торжестве ОТО в интерпретации феномена «чёрных дыр». Таково современное положение дел в этой области астрофизики.
Однако приведённые выше выводы некорректны, - просто потому, что они основаны на ОТО с её физически недопустимыми сингулярностями, и они - не единственно «возможная» интерпретация наблюдательных данных. Мной предлагается другая их интерпретация, минимально «нагруженная» недоказуемыми допущениями типа ОТО.
Количественная модель нейтронных звёзд и «чёрных дыр», основанная на классической механике и теории гравитации И. Ньютона
Предлагаемая мною модель состоит в принятии простейшей теории гравитации - закона всемирного тяготения Ньютона и предположения о том, что фотоны также испытывают воздействие этого тяготения (как и «массовая» материя). Это допущение основано на зарегистрированном прямыми наблюдениями во время полного солнечного затмения факте искривления световых лучей в гравитационном поле Солнца, интерпретированном блестящим астрономом сэром Артуром Эддингтоном (правда, на базе ОТО).
Предположения этой модели были сделаны впервые очень давно - Джоном Митчеллом в 1783 г. и, несколько позднее, Пьером Симоном Лапласом, в 1798 г. Из этих предположений следует существование гравитационного радиуса, определяемого простой формулой
2 Gm
Rg =— ' (1)
где G - гравитационная постоянная Ньютона, т - масса гравитирующего тела, c - скорость света. Удивительным является совпадение выражения (1) для Rg у Митчелла - Лапласа и в решении Шварцшильда, несмотря на существенно разные, физически и математически, исходные положения в двух теоретических подходах.
Приняв отмеченные два допущения, можно решить задачу Оппенгеймера - Волкова для любого уравнения состояния типа p = f (р) и убедиться в том, что в этой модели нет эффекта несуществования решения задачи при «больших» величинах суммарной массы объекта т. В частности, при простейшем виде уравнения состояния р = const = р., что является хорошим приближением в статических задачах рассматриваемого типа (но не в динамических - см. ниже), решение получается в явном виде в элементарных функциях. Если в этом решении в качестве р.
1 Oppenheimer J. R., Volkoff G.M. "On Massive Neutron Cores." Physical Review 55 (1939): 374-381.
2 Oppenheimer J. R., Snyder H. "On Continued Gravitational Contraction." Physical Review 56 (1939): 455-459.
принять плотность нейтронного вещества (р. = рп ~ 1014 г/см3), то получается, что радиус Rb такого тела при его массе т < т. оказывается больше гравитационного (по Митчеллу - Лапласу)
3 т 2 Gm
Rb = 3 -- >—— = Rg , т < т.. (2)
\4лр. с2 g
В области т > т. верно обратное неравенство, т.е. такое тело есть «тёмная звезда» Митчелла - Лапласа. Когда Rb = Rg , р. = р.. , и для р.. получается формула
_ _3_
р«.«._
4 л
( ^ V
2 G
m2 V m
2
1 _ 2 • 1016 [ I (г/см3)' (3)
где т., - масса Солнца. Отсюда видно, что при р.. = рп ~ 1014 г/см3 получается т = т. ~ 14т,. Таким образом, в области т < 14т, рентгеновские источники должны бы быть нейтронными звёздами, в то время как признаками таковых обладают лишь объекты с т < 2т, ^3т, .
Моё объяснение этого факта состоит в следующем. В современной физике элементарных частиц известно, что нейтроны (и протоны) состоят из «более элементарных» частиц - кварков, которые, как и нейтроны, являются фер-мионами1 и «газ» из них описывается уравнением состояния «холодного» вида р = f (р) (конечно, с другой функцией /, чем у нейтронного «газа»). Понятно, что при сжатии гравитационными силами достаточно большой массы т произойдёт «коллапс» нейтронного «газа» - при определённом пороговом значении давления р. произойдёт скачкообразное сжатие почти всего нейтронного вещества с возрастанием плотности от р. = рп до плотности кваркового вещества р. = рд. А дальнейший рост давления будет сжимать уже кварковое вещество2. Для объекта заданной массы его радиус Rbq в «кварковом» состоянии будет значительно меньше радиуса в «нейтронном» состоянии, Rbq < Rbn, поскольку рq > рп. А это означает, что при р = р. = рп с ростом массы т будет достигнуто такое его значение, при котором резко возрастёт плотность и уменьшится радиус тела - произойдёт кварковый коллапс. При этом гравитационный радиус этого тела Rg не изменится Я линейно зависит от т, но не зависит от р). Таким образом, произойдёт скачкообразное падение Rb от значения Rbn > Rg до значения Rbq < Rg , т.е. тело из состояния «нормальной» нейтронной звезды скачком превратится в невидимую звезду Митчелла-Лапласа. Детальное исследование этой проблемы, приведённое ниже, обнаруживает удивительные особенности решения, имеющие очень важное значение в качественном и количественном отношениях.
Итак, принимаем уравнения состояния для нейтронного и кваркового веществ в простейшем виде
¡Рп, р < р.; Р =\ (4)
{РЧ , р > р.;
где, как отмечалось, р. - «пороговое» значение давления, при котором нейтронное вещество скачком переходит в сильно уплотнённое кварковое. Конечно, представление (4) является идеализацией, ибо и при р < р. будет какая-то зависимость рп = рп(р), но с малым изменением рп при росте р , а при р > р. будет р = р^) с еще меньшей сжимаемостью кваркового вещества. Но эти малые эффекты не изменят основных свойств решаемой задачи, для которой представление (4) вполне приемлемо. Уравнение равновесия в радиально-симметричной задаче для вещества при ньютоновом законе тяготения имеет вид
У
dr r
(r), M (r)_ 4 л J p(4 2d4 , (5)
Г - J
где г - радиальная координата, отсчитываемая от центра системы, р (г) - значение плотности на этом расстоянии, М (г) - суммарная масса вещества внутри сферы радиуса г. Для представления (4) имеем из (5) (если всюду р = рп)
р = р )-^ G рП г 2, (6)
и, принимая р(Я) = 0, R - внешний радиус тела (в действительности, р(Я) = р0 > 0, где р0 - давление, при котором любое вещество сжимается до нейтронного состояния; однако, как будет ясно из дальнейшего, р0 намного меньше интересующих нас давлений, близких к значению р., при котором происходит «смятие» нейтронной массы в квар-ковую), получаем
р (0 )= ^ G рП Я2. (7)
Это - максимальное давление в недрах тела и именно оно будет достигать значения р. при Я = Я. и массе М.,
1 Сисакян А.Н. Избранные лекции по физике частиц. Дубна: ОИЯИ, 2004. 377 с.
2 Glendenning N.K., Weber F. "Nuclear Solid Crust on Rotating Strange Quark Stars." Astrophysical Journal 100.2 (1992): 647658; Negreiros R.P., Mishustin I.N., Schramm S., Weber F. "Properties of Bare Strange Stars Associated with Surface Electric Fields." Physical Review D82.10 (2010):103010-103016.
поэтому
Я„ —
Г \ 1/2
Г 3 р. ^
2л Gp
п У
4л 3 М. — — РпЯ. .
(8)
Дальнейшие изменения состояния тела будут характеризоваться образованием вокруг его центра однородного ядра из кваркового вещества с радиусом г. и плотностью рд и «мантии» из нейтронного вещества, заполняющей сферический слой в интервале г.< г < Я. При этом для М(г) получается формула
М (г ) —
4л 3
— рдг , 0 < г < г. ;
4л
3 4л / 3 3\ 3 рдг» + —Рп(г -г. А г.<г < Я;
(9)
а для распределения давлений - формула
Р —
р.+ ^ ^^р¡2 (г.2 - г2), 0 < г < г. ;
2л „ 2
р ~3Г GРq
(
г.2 - г2 - 2 г.3
Р
Р
Л
- 1
V п у
^ - 1 гг
г < г < Я
(10)
Давление р на границе ядра и мантии непрерывно и равно р.. При г = Я из (9) получаем
М , X — (1 + V,3 )у3,
м. у г
(11)
а из (10) при г = Я и Я (р) = 0 - соотношение
я
Я
■ = У — [1 + (2v - 1)и2 - 2vu3]-
- 1) и2 - 2vu3 1-1/2
(12)
Здесь приняты обозначения
Рд , V — —- - 1
Рп
и
Я
(13)
Соотношения (11), (12) задают в параметрическом виде зависимость У = У (X, V), которая для однородной (без скачка плотности) среды переходит в такое выражение:
У — X1/3 (у (и, 0) — (1 - и2 )1/2, X (и,0) — У3 — (1 - и2) 3/2).
[и, 0) — - и
Асимптотика зависимостей (11), (12) при и ^ 0, X ^ У имеет вид
У — 1 - 1 [^ - 1 )и2 - 2vu3 ]+ о (и4)— 1 -IV - 1 | и2 + vu3 + ...
1
X — Л - 3 [IV - 1 )и2 - 2vu3 ]+ о (и4 М + vu3)— 1 - V - 1
- — |и2 + 4vu3 + ... .
(14)
(15)
Из этих представлений следует, что для 0 < V <1/2 с ростом и от нуля X и У растут от их значений 1 при и = 0, причём кривые У = У^, V) в точке X = У=1 имеют общую касательную с предельной при V = 0 кривой У = X . В области 1/2 < V < да весь пучок этих кривых имеет такое же касание. При V = 1/2 асимптотика имеет вид
X — 1 + 2и3 + ..., У — 1 + 1 и3 ,
2
то есть эта единственная кривая пучка отходит от точки X =У =1 под конечным углом к кривой У = X1'3 . При и^1 получается асимптотика
У ^ [2 (1 + V )(1 - и )]-1/2 , X ^(1 + V )у 3 ,
2
г
то есть
Y =
что при v = 0 совпадает с Y = X , а при v ^ да имеем:
X
1 + v
1/3
Y ^
Г Л1 p. x
q У
Картина распределения семейства кривых Y = Y(X, v) в показана на рис. 1.
Удивительным фактом является то, что при v > 1/2 росту радиуса r. кваркового ядра (росту u) соответствует не только уменьшение Y, т.е. общего радиуса тела R, но и уменьшение X, т.е. массы тела! Таким образом, зависимость Y = Х1/3 чисто нейтронного тела с переходом к нейтронно-кварковому образованию с ростом u сначала совершает «попятное» движение в область X <1 и лишь затем выходит на монотонную ветвь с приведённой выше асимптотикой. Всё это означает некую неустойчивость построенного статического решения в области 0 < X < 1: сечение прямой X = const < 1 с кривыми семейства определяет три точки - одну на Y = Х1/3 < 1 (чисто нейтронная конструкция) и две на ветви Y = Y(X, v) , которая при некотором X. = X.(v) <1 имеет вертикальную касательную. А это означает, что достаточно интенсивное динамическое воздействие на чисто нейтронную звезду в диапазоне X.(v) < X < 1 может привести к образованию кваркового ядра конечных размеров, сопровождающемуся резким (скачкообразным) уменьшением общего размера тела R , причём такой скачок может произойти в две возможные точки - одну на ветви, идущей от точки X = Y = 1 до точки X = X.(v) , другую - на нижней ветви, идущей от точки X = X.(v) вниз и затем, пройдя точку X=X..(v), где касательная к кривой Y = Y(X, v) горизонтальна, вверх к асимптотике:
N1/3
Рис. 1. Распределение семейства кривых Y = Y(X, v).
Y =
X 1 + v
при X ^ да.
Скачкообразный переход в первую из этих точек следует считать нереализуемым, ибо из него дальнейший рост ядра (параметра и) соответствует уменьшению общей массы (уменьшается X). Для второй точки такой переход вполне реален, ибо квазистатический рост ядра из этого состояния определяется ростом общей массы системы (например, за счёт её плавного нарастания в процессе спокойной аккреции вещества извне). Принципиально важно, что при таких скачках радиус тела может стать меньше гравитационного радиуса Митчелла-Лапласа, и световые сигналы из окрестностей возникшей скачком нейтронно-кварковой конструкции будут не в состоянии выйти за пределы этого радиуса, то есть конструкция будет представлять из себя «чёрную дыру»!
Скачкообразное уменьшение У из точки X = У = 1 при постоянном X = 1 будет происходить обязательно, ибо точка X = У = 1 выделена тем, что в центре нейтронного тела давление уже достигло значения р. - давления перехода нейтронной материи в кварковую, поэтому для возникновения скачка не требуется сильных (конечной интенсивности) воздействий - для этого будет достаточно малых динамических воздействий. Таким образом, после достижения системой состояния X = У = 1 (М = М.(р.), Я = Я.(р.)) существование чисто нейтронных звёзд становится невозможным, и дальнейший, даже очень малый, прирост массы скачком переводит систему в нейтронно-кварковую звезду с кварковым ядром конечного (не малого) радиуса (именно этот случай отмечен выше, в общей части статьи). При этом радиус такого образования может быть и больше гравитационного радиуса для общей массы конструкции, тогда по всем признакам наблюдательных данных она будет классифицироваться как «нейтронная» звезда. Понятно, что это будет в случае небольших значений V, т.е. небольшого превышения плотностью кваркового вещества плотности нейтронного. Но поскольку реальные оценки, основанные на современной теории элементарных частиц1, дают для V величину в несколько (много) десятичных порядков, то, как будет показано ниже, скачок из точки X = У =1 дает величину У = Ус, существенно меньшую величины Уg = Яg/Я..
А из этого следует, что именно величина массы М = М. ограничивает сверху диапазон значений масс нейтронных звёзд. Для всех значений М > М. нейтронно-кварковая конструкция будет «чёрной дырой» («нормальным» небесным телом без сингулярностей), ибо для них будет соблюдаться условие Я = Яс = Я.Ус < Я»У^1) = Я^М»).
Определим теперь, при каких значениях параметров задачи (р.у и т.п.) будет соблюдаться условие Ус^) < Уг. Из формулы (1) и формул (8) для М., Я. следует, что
= R^ = 2GM, = g R c2 R pc2
1 Сисакян А.Н. Указ. соч. 96
Легко показать, что асимптотика величины Ус = У(Х = 1^) при V > 1 есть Ус =у 1/3 +..., и, таким образом, условие У: > Ус представляется при V > 1 в виде
V-"3 < 4р
РС*' (17)
Из данных наблюдений для нейтронных звёзд и «чёрных дыр» известно, что при критической массе М = М* = 3М* будет < R*, т.е. У: <1. Действительно, при R = Rg плотность тела выражается простой формулой (3), и чтобы было р = рп ~ 10 г/см , должно быть М ~10М , что втрое больше М*, т.е. при М = М* = 3М* будет R = R* > Rg. Таким образом, окончательно при V > 1 получаем ограничения на р*
. 2 2 .1/3 , 4 р.. Рпс _ ^ рпс
<
пп
< 1 или —Ц3 < р. <—-— . (18)
Рпс2 4у 4
Полагая рп = 1014 г/см3 = 1016 кРсек2/м4, при с2 = 9^1016 м2/сек2, получаем оценку:
РпС 2
При V = 103 имеем:
,2
6
а при V = 10
Рпс!
4у1/3
2,5 • 10 28 кГ/см2 . Р4т - 2,5 -1027 кГ/см2, 2,5 -1026 кГ/см2 .
Таким образом, для р* получились оценки:
2,5 • 1026 ■ 2,5 • 1027 кГ/см2 < р. < 2,5 • 1028 кГ/см2,
ибо оценки 103 < V < 106 для V = рд /рп - 1 представляются разумными1. Полученные числа позволяют оценить характерную величину давления перехода нейтронного вещества в кварковое. Здесь мы имеем замечательную ситуацию: «измерение» р* осуществляется построенной рациональной моделью макрокосмических объектов очень больших масс и сведениями из современной теории элементарных частиц, а имея оценку р* из других данных о макрокосмических событиях рассматриваемого масштаба (например, о динамических процессах), можно по полученным здесь соотношениям оценить характеристику микромира рд (где рд = (1 + V) рп), что другими средствами, по-видимому, и невозможно сделать. В этом отношении мир нейтронных звёзд и «чёрных дыр» - это лаборатория для исследований по физике элементарных частиц в условиях очень больших, конечно, недостижимых в земных условиях, давлений и плотностей материи.
Подводя итог изложенному выше, можно заключить, что регистрируемое изменение характеристик рентгеновского излучения, исходящего от тёмных объектов (нейтронных звёзд и «чёрных дыр»), в области значений масс М > 3М* , связано не с эффектами ОТО, а со скачкообразной перестройкой нейтронной звезды в намного более компактное образование - не имеющую никаких сингулярностей нейтронно-кварковую звезду, приобретающую основную характеристику «чёрной дыры» - способность «запирать» фотоны в её окрестности, ограниченной классическим гравитационным радиусом Митчелла - Лапласа. Таким образом, построенная модель является подтверждающей гипотезу Митчелла - Лапласа теорией, опирающейся вовсе не на ОТО с её физически неприемлемыми математическими сингулярностями, а на классическую механику, на закон тяготения Ньютона и представления современной физики элементарных частиц.
Для завершения исследования вопроса о финальных стадиях эволюции звёзд нужно рассмотреть динамические задачи для нейтронных и нейтронно-кварковых звёзд, что и сделано ниже.
Уравнение состояния вида р = /(р) характеризует сплошную среду с разделяющейся внутренней энергией. Для такой среды термодинамическая модель имеет вид
е(у,*)=еи+*2И, р = /И—ОУ>0,
аУ
(Т) а (19) * = s(т)=\СyTT)dT , е2 = ^2 (т), су = ОТ = Су (т),
где е - внутренняя энергия, * - энтропия, Т - абсолютная температура, У =1/р- удельный объём2. Рассматриваемая нейтронно-кварковая среда есть частный случай среды (19), у которойр = /(У) имеет вид, изображенный на рис. 2А.
Для такой среды рассмотрим задачи об «ударном» нагружении и «ударном» разгружении, когда к поверхности среды скачкообразно прикладывается высокое давление или, наоборот, - начальное высокое давление в среде скачкообразно снимается на её поверхности. Эти задачи автомодельны и допускают точные решения. Не приводя здесь дета-
Там же.
1
2 Григорян С.С. О некоторых специальных вопросах термодинамики сплошных сред // Приладная математика и механика. 1960. Т. 24. Вып. 4. С. 651-662.
лей построения (очень простых) этих решений и результирующих их явных формульных описаний, представим интересующие нас здесь результаты графически на рис. 2В и рис. 2С, соответственно.
1 А
• ! 1 ! 1 1 Ъгт • \ • \
1 1 1 1 1 1
К*
Vqi Vq*
ln*VnO
Рис. 2. Диаграммы состояния нейтронно-кварковой среды.
К* К
ft* пз
Приложение к поверхности полупространства из нейтронного вещества, находящегося под давлением p0, давления pi, превышающего «критическое» давление p», порождает автомодельное решение с бегущими вглубь полупространства двумя ударными волнами - в одной происходит скачок из точки V0, p0 в точку Vn(p»), р» (быстрая волна), в другой - скачок из точки Vn(p»), р» в точку Vq(pi), pi. Это - единственное построение, удовлетворяющее второму закону термодинамики. По этому решению, в среде происходит нагрев, мера которого показана в виде затемнённых областей на рис. 2B (нагрев происходит в два этапа на фронтах двух ударных волн). Величина нагрева (рост температуры в скачках) определяется соотношениями:
e2(T») - e2(T0) = величине затемнённой площади для первого скачка, e2(Ti) - e2(T») = величине затемнённой площади для второго скачка.
Снятие начального давления pi, имеющегося в кварковом полупространстве, за счёт скачкообразного понижения давления на его поверхности до давления ps < p» приводит к картине, изображенной на рис. 2С: по среде распространяется непрерывная волна разрежения из точки Vq(pi), pi до точки Vq(p»), p» - это тоже быстрая волна, за ней следует скачок разрежения из точки Vq(p»), p» до точки Vn(ps), ps. Этот скачок удовлетворяет требованию второго закона термодинамики о росте энтропии s и температуры T на скачке, и нагрев на скачке определяется соотношением:
e2(Ts) - e2(Ti) = площади затемнённой области на рис. 2C.
Величины нагревов как в случае ударного нагружения, так и в случае ударного разгружения нейтронно-кварковой среды весьма велики, и их легко оценить по построенным решениям. Действительно, из-за того, что кривые Vq = Vq(p) и Vn = Vn(p) пологи к оси давлений, их просто можно заменить прямыми линиями Vq = const, Vn = const, и тогда затемнённая площадь Fc для второй (главной) ударной волны на рис. 2B будет равна
Fc = 0,5(p1 - p.)[Vn(p.) - Vq(p.)]= 0,5(p1 - p.Xv - Vq), (20)
а для рис. 2C затемнённая площадь Fr будет равна
F = 0,5 (p. - ps )Vn (p.)-Vq (p.)]= 0,5(p.- ps )(Vn - Vq ), (2i)
Особый интерес представляет величина нагрева в скачке разрежения, определяемая по (21) приps = 0 и Vq < Vn:
Pn
и если для оценок принять cVn(T) = const =cVn, cVq(T) = const =cVq, из (22) получится
p.
e2T)-e2(T)« 0,5p.Vn = 0,5^,
(22)
c T - c T « c T = 0,5
vn s vq i vn s '
Pncv,
(23)
Однако, поскольку мало что известно о теплоёмкостях нейтронного и кваркового веществ Су„ и cVq , формула (23) не позволяет сделать оценку нагрева в величинах абсолютных температур.
Правая часть (22), определяющая рост тепловой составляющей внутренней энергии среды, при оценках р* и ря, приведенных выше, составляет:
0,5-2,5-1027 кГ/см2:1014 г/см3 ~ 1027 кГ/см210-16 м4/кГсек2 = 1015 м2/сек2
- при средней оценке для р*. Для крайних оценокр* результат будет отличаться на один порядок в каждую сторону, т.е. 10 м /сек < 0,5 р*/ря < 10 м /сек . Эти величины составляют от 10- с до 0,1с , где с ~ 10 м /сек - квадрат скорости света. Таким образом, в скачке разрежения получается «чудовищный» нагрев нейтронной среды, который эту холодную при её образовании материю превращает в очень сильно нагретый, уже разреженный нейтронный газ, разлетающийся в пустоту; у реальных нейтронно-кварковых звёзд из-за трёхмерности процесса - со значительным дальнейшим адиабатическим расширением и охлаждением. В результате этого нейтроны «развалятся» на протоны,
электроны и множество других лептонов и фотонов, и при достижении в разлетающемся облаке температуры рекомбинации водорода последний будет массово синтезироваться!
Таким образом, описанный динамический процесс ставит точку в конце эволюции звёзд (из водорода), когда их холодные нейтронно-кварковые «останки» при столкновениях полностью разрушаются и очень энергично разлетаются в сильно нагретом состоянии (за счёт накопленной в них до динамического разрушения «холодной» внутренней энергии, «взятой» из собственного гравитационного поля «умершей» звезды) и на определённой стадии адиабатического охлаждения «приготовленной» таким образом плотной газовой среды вновь превращаются в исходное вещество - газообразный водород и смесь лептонов и фотонов.
Интенсивность нагрева холодного нейтронно-кваркового вещества в случае столкновения двух объектов из подобной материи и, следовательно, энергетика последующего его разлёта в ударных волнах разрежения и адиабатического охлаждения с образованием водорода будет весьма высокой, если скорость «лобового» столкновения нейтрон-но-кварковых останков звёзд будет достаточно велика. В этом случае нагрев будет двухступенчатым (если пренебречь «слабым» нагревом в предвестниковой быстрой ударной волне) и определяться величинами затемнённых областей на рис. 2В и рис. 2С (конечно, детали «газодинамических» расчётов будут уже не столь простыми, как в рассмотренных выше двух простейших задачах, но для оценок результирующей интенсивности нагрева можно нагревы на рис. 2В и рис. 2С сложить, связав предварительно расчётом величину результирующего максимального давления соударения р1 (рис. 2С) со скоростью v0 столкновения рассматриваемых тел, что с помощью описанных выше автомодельных решений делается очень просто - с явными формулами. Эти формулы таковы (при р = р„(р) и р = рд(р) для нейтронной и кварковой материй, соответственно)
Рг =
Povo
1 _ р^ Рт
< Р.
v < v„ =
(
Ро
Ро = Рт (0) , Рт. = Рт (р. ) ;
Ро
Л
1
V Рт. J
(24)
Рг = Р. +
Рт.^о _ V0.)2
1 _
Рт
А
(Рг )
< Рг. , vo. < V < V
о*.
(25)
o«
1 _
Ро
р
v =_ь_
q.. J
1 _
Ро
Р q**= Ра (Рг.) ,
(26)
Р
т.. J
где V)««, р,-» - значения Vo и р¡, соответствующие условиям, когда скорости предвестниковой и основной ударных волн сравниваются (исчезает двухвол-новая конфигурация);
Рг =
Ро*о
1 _
Ро_
Раг
Рог = Ра (Рг ) .
(27)
Рис. 3. Пятно контакта двух соударяющихся выпуклых тел.
По этим формулам можно определить давление p¡ на пятне контакта двух соударяющихся выпуклых тел на самых ранних стадиях соударения, когда образующиеся в них ударные волны представляют собой очень пологие поверхности и присоединены к границе пятна контакта (рис. ЗА; АА- проекция пятна контакта, АВ1А', АВ2А' - проекции ударных волн на плоскость фигуры). С течением времени наступает момент, когда эти поверхности отрываются от границы пятна контакта, и образуется кольцевая быстро расширяющаяся свободная поверхность (с нулевым давлением на ней), через которую сжатое ударными волнами вещество разлетается в виде тонкой дискообразной струи с огромными скоростями (рис. ЗВ; А1СА2, А/СА2' - проекция струи).
Такая задача о начальной стадии динамики лобового соударения выпуклых тел из обычных материалов с образованием струи была рассмотрена в моей работе1. Скорости в струе получаются очень большими в сравнении со скоростью соударения v0 - возникает своеобразный кумулятивный эффект. Для рассматриваемого случая соударения нейтронных и нейтронно-кварковых тел этот эффект, очевидно, имеет важное значение - он, в частности, определяет скорости и конфигурацию облака разлетающейся сильно нагретой материи.
По формулам (24) - (26) можно оценить порядок скоростей соударения, нужных для возникновения удар волн сжатия и последующего разрежения. Из (24) для v0* получаем при полученных выше оценках для р.:
1 Григорян С.С. О динамике начальной стадии соударения тел с большими скоростями // Некоторые вопросы механики сплошной среды: Сб. научных статей / Под ред. С.С. Григоряна. М.: Изд-во МГУ, 1978. С.157-172.
1 - А = 1(ГкУ0* = 1 0-к/^2,5-1027 кГ/сМ :1с14 г\смМ = 10-к/2 • 1,6-10 км/сек.
Рп*
И при к = 2, 4, 6 получаем, соответственно, у0. ~ 1,6 103, 1,6402, 16 км/сек.
Поскольку зависимость рп = рп(р) неизвестна, приведённые значения к можно принять для оценки порядков величин. Принимая во внимание, что
Ро „ Рс
< 1.
Р а * * р
(Рг*)
по (26) имеем для у0 = у0..
V,
0** 1 п.к
Vо*
= 10к •
( \ 0
Ра
10к
V г а** J
т.е. получаются значения у0.. ~ 1,6405, 1,6^ 106, 1,6-107 км/сек.
Уже при к = 2 получается у0.. порядка скорости света с. Поскольку полученные значения у0.. очень велики, это означает, что при реальных скоростях соударений состояния с одноволновой конфигурацией недостижимы. Если же у0 = ау0., то для р1 по (24) и (25) получается
( „Л
Рг - Р* ~а2Р>
1 -
Р
0
V Рп* J
= 10гкр * а2
и чтобы при этом было р1 - р. < р., нужно а < 10к/2. Для принятых выше оценок к имеем а ~10, 100, 1000, а для у0 получаем у0 ~1,6-104 км/сек (независимо от к!). Полученные числа позволяют оценить нагрев в ударных волнах сжатия, предшествующий нагреву в ударной волне разрежения, - чтобы эти нагревы были одного порядка, нужно иметь скорость соударения у0 ~ 1,6-10 км/сек.
Рассмотренные выше динамические эффекты, возможные при тесных взаимодействиях нейтронных и нейтрон-но-кварковых звёзд (их непосредственных столкновениях друг с другом или при близком прохождении, вызывающем интенсивные приливные эффекты, выводящие их из устойчивого состояния), казалось бы, маловероятны из-за разреженности звёздного «населения» Вселенной. Однако это население собрано в сгустки - звёздные скопления, галактики, скопления галактик, в которых распределение звёзд и газопылевых облаков материи существенно неоднородно - плотность распределения возрастает при приближении к центру скопления и достигает там значительных величин, образуя ядро скопления. Ядра многих галактик активны - наблюдаются мощные радиационные потоки из них - проявления протекающих там интенсивных физических процессов, которые, бесспорно, связаны с тесными взаимодействиями плотных сгустков массы (светящих и тёмных нейтронных и нейтронно-кварковых звёзд и газопылевых облаков). Факт значительного возрастания со временем плотности распределения звёзд в ядрах скоплений в процессе их эволюции был теоретически установлен ещё в ранней работе В.А. Амбарцумяна1. Таким образом, существует эволюционный механизм образования у звёздных скоплений густонаселённых ядер, где рассмотренные выше динамические эффекты становятся реальными. Об этом свидетельствует, в частности, открытый В.А. Амбар-цумяном в конце 40-х годов прошлого столетия факт группового рождения звёзд из единых материнских образований - открытие звёздных ассоциаций, состоящих из очень молодых звёзд. Это открытие привело В.А. Амбарцумяна к формулировке гипотезы о существовании протозвёзд2 - неких невидимых объектов, которые, разрушаясь взрывным процессом, рождают обычные (светящие) звезды. Таким образом, представленная здесь моя теория подтверждает эту гипотезу В.А. Амбарцумяна, а существование ассоциаций из молодых одновозрастных звёзд подтверждает сделанный из этой теории вывод о реальности динамического распада нейтронно-кварковых звёзд, приводящего к восстановлению из их материи водорода, из которого вновь образуются «обычные» звёзды. Эта теория одновременно является подтверждением и исходной гипотезы Митчелла-Лапласа о существовании «тёмных» звёзд. Из этой теории также следует объяснение природы так называемой тёмной материи во Вселенной, то есть факта необходимости предположения о существовании в звёздных системах (скоплениях, галактиках) невидимой массы, гравитационное воздействие которой на видимые (светящие) звёзды позволяет количественно «свести концы с концами» в динамике этих систем (при этом оказывается, что тёмная материя составляет около 90% всей массы системы).
В самом деле, поскольку, с одной стороны, в основном объёме звёздной системы (вне её ядра) звёзды расположены очень далеко друг от друга, так что их столкновения маловероятны, а с другой стороны, продолжительность «жизни» светящих звёзд относительно мала в сравнении с таковой для самих звёздных систем, то из этого следует, что основной объём этих систем «заселён» уже «умершими» звёздами - нейтронными и нейтронно-кварковыми звёздами, которые «живут» там неопределённо долго, лишенные возможности тесных взаимодействий с последующей регенерацией водорода. Новые же светящие звёзды в этот объём заносятся коллективной динамикой в плотном по населению ядре, где процессы регенерации водорода из звёздных останков возможны и постоянно протекают. Здесь уместно ещё отметить, что выявляемое наблюдениями существование в активных ядрах галактик «чёрных дыр» с массами порядка 10-10 масс Солнца не обязательно означает, что эти образования представляют собой единые объекты, - это могут быть скопления нейтронных и нейтронно-кварковых звёзд, и коллективное их воздействие на соседние объекты (обычные звёзды, газовые и пылевые облака), порождающее множественную аккрецию
1 Амбарцумян В.А. К вопросу о динамике открытых скоплений // Ученые записки Ленинградского гос. университета. 1938. № 22. Серия математических наук (астрономия). Вып. 4. С. 19-22.
2 Амбарцумян В.А. О протозвёздах //Доклады АН Арм.ССР. 1953. Т. 16. № 4. С. 97-102.
материи на них, формирует суммарный эффект рентгеновского и иных излучений, регистрируемых наблюдениями. К числу этих «иных» излучений следует отнести космические лучи - потоки частиц с очень высокими энергиями. Таким образом, построенная здесь теория объясняет также и природу генерирования космических лучей. Здесь же отметим, что она может дать объяснение и вероятной природы квазаров.
Наконец, нужно отметить, что ещё много лет назад в одной из моих работ1 я высказал гипотезу о возможности существования «чёрных дыр» малой плотности - разреженных «островных» систем во Вселенной, у которых гравитационный радиус больше радиуса (размера) самого материального скопления. В такой системе излучение (фотоны) будет «заперто», оно будет в нём неограниченно «блуждать», поглощаться сгустками массы и переизлучаться, в результате образуется некий средний фон равновесного излучения. В той же работе высказана гипотеза о том, что наблюдаемая часть Вселенной, вероятно, и есть такая островная система, а так называемое реликтовое излучение -такой фон запертого в «чёрной дыре» равновесного излучения. Это заключение хорошо согласуется с установленным наблюдениями фактом очень малой величины равновесной чернотельной температуры (единицы градусов Кельвина) реликтового излучения, поскольку основное население этой системы - холодные тёмные останки светящих звёзд. Построенная мною теория подтверждает высокую вероятность реальности положений также и этих двух гипотез.
В той же работе мною высказано также предположение о том, что так называемый Большой взрыв - это событие в конечной наблюдаемой части Вселенной, в такого рода островной системе. Здесь к этому можно добавить, что это событие могло быть результатом сжатия островной системы на определённой стадии её эволюции, когда основная её масса собралась к центру в виде плотного скопления «умерших» звёзд, которое по описанному здесь механизму «взорвалось» и превратило всю эту массу в разлетающуюся сильно нагретую «элементарную» материю, а далее всё происходило по существующей ныне теории эволюции Большого взрыва.
Изложенное выше представляет собой простейшее построение, основанное на теории тяготения Ньютона и современных представлениях о физике микромира, показывающее, что «чёрные дыры» это вполне «нормальные» физические образования без каких-либо сингулярностей, согласующееся с существующими представлениями о строении материи и управляющих её «жизнью» силах и не требующее интеллектуальных сверхусилий для размещения этих образований в прокрустовом ложе ОТО, которая приводит к физически неприемлемым предсказаниям и поэтому непригодна для описания природных явлений. Конечно, построенная здесь конструкция может претерпеть какие-то количественные изменения в деталях, если использовать теорию гравитации, приемлемым образом уточняющую закон тяготения Ньютона и не порождающую сингулярностей. Но всё это уже «не принципиально».
Заключение
В астрономии - астрофизике - астромеханике попытки определить массы звёздных объединений - галактик по числу (надёжно оцениваемому) светящих звёзд и оценке индивидуальных звёздных масс (по разнообразию оптических свойств звёзд) привели к жёсткой необходимости сделать предположение, что общая масса всех таких звёзд в данном звёздном кластере - галактике существенно (на один порядок величины) недостаточна для обеспечения (по закону всемирного тяготения Ньютона) величин, определяемых в наблюдениях характеристик её динамики - скоростей звёзд, их распределения в объёме галактики и т. п. И это «заставило» предположить, что в галактиках существует несветящая -тёмная материя, также создающая гравитационное (тоже по Ньютону) силовое воздействие на светящие звёзды и определяющая их индивидуальную и «коллективную» динамику. Это представление о «тёмной» материи получило всеобщее признание в науке о мироздании, однако никаких вразумительных гипотез и наблюдательных результатов, позволяющих судить о её природе, до сих пор в этой науке нет, если не считать таковыми чисто умозрительные высказывания о том, что должны существовать некие «элементарные» носители «тёмного» гравитационного феномена.
В последние годы я построил2 рациональную количественную теорию механики и физики регистрируемых современными средствами наблюдений сверхплотных небесных тел Митчелла - Лапласа - нейтронных и нейтронно-кварковых звёзд - останков «умерших» звёзд и сформулировал гипотезу о том, что «создателями» тёмного гравитационного феномена в галактиках являются эти сверхплотные холодные несветящие макротела, составляющие подавляющую часть населения галактик из-за того, что время жизни галактик значительно больше времени жизни звёзд (особенно звёзд массивных), так что в «великовозрастных» галактиках накапливается очень много этих сверхплотных останков умерших звёзд, которые и определяют «гравитационную погоду» в таких галактиках. Так это или нет, должны, как всегда, решить результаты астрономических наблюдений. Ибо именно наблюдения играют в астрономии - астрофизике - астромеханике две главные роли - «добытчика» фактов и верховного судьи для всех теоретических построений (в их числе - и бесплодных математических измышлений, ставших в последние десятилетия доминирующим жанром в астрофизике - астромеханике и в теоретической физике вообще). Именно поэтому великий астрофизик академик В.А. Амбарцумян часто повторял: «Наблюдайте, наблюдайте, наблюдайте!», - ибо в наблюдениях выявляются реальные факты, и только они определяют истинность теоретических построений.
Так вот, относительно недавно были обнародованы3 результаты последних наблюдений особенностей строения и радиационных свойств огромной популяции чрезвычайно далёких (exceedingly distant) окутанных пылевыми «ко-
\ Григорян С.С. О динамике и строении Вселенной // Проблемы современной механики. Ч. \. Сб. научных статей / Под ред. акад. Л.И. Седова. М.: Изд. МГУ, \983. С. 38-43.
2 Григорян С.С. О природе наблюдаемых во Вселенной «чёрных дыр» // Ломоносовские чтения. Тезисы докладов научной конференции. Секция механики. Москва, апр. 2005 г., МГУ имени М.В. Ломоносова. М.: Изд-во Моск. ун-та, 2005. С. 82; Он же. О природе наблюдаемых во Вселенной «чёрных дыр» и динамике Вселенной // Известия Национальной Академии наук Армении. Серия Механика. 2007. Т. б0. № 4. С. 3-Î9; Он же. Гипотеза С.С. Григоряна о природе «тёмной материи» в галактиках и новые данные астрономических наблюдений, подтверждающие её // Ломоносовские чтения. Тезисы докладов научной конференции. Секция механики. Москва, \б-25 апр. 20\2, МГУ имени М.В. Ломоносова. М.: Изд-во Моск. ун-та, 20\2. С. 57-58; Он же. О природе «тёмной материи» в галактиках // Доклады Национальной Академии Наук Армении. 20\2. Т. \ \2. № 3. С.273-27б.
3 Amblard A., Cooray A., Blain A., Bock J., Dowell C.D., Levenson L., Lu N., Nguyen H.T., Schulz B., Shupe D.L., Vieira J.D., Xu C.K., Zemcov M. "Submillimetre Galaxies Reside in Dark Matter Haloes with Masses Greater than 3* Ю^ Solar Masses." Nature 470.7335 (20\\): 5\0-5\2; ESA. "Hershel Finds Less Dark Matter But More Stars (From ESA News Release, \б Febr. 20\\)." Space Research Today: COSPAR's Information Bulletin. \80 20\\: 5\-52.
конами» галактик, в которых идёт (точнее, шёл - это было очень давно на световом конусе) бурный процесс звездообразования. Наблюдения были проведены в дальней инфракрасной субмиллиметровой части спектра космическим телескопом «Гершель» (Европейское Космическое Агентство). Каждая из этих галактик имеет массу порядка 3-1011 масс Солнца, чего явно недостаточно, по существующим теоретическим представлениям, для «зажигания» интенсивного звездообразования (для этого требуется масса галактики в 5-1012 масс Солнца). Из этого факта был сделан вывод, что в этих галактиках нет нужного преобладания (на порядок величины) массы тёмной материи над светящей массой, что наблюдается в «обычных» галактиках, и что он делает серьёзный вызов существующим представлениям теории насчёт природы «тёмной» материи. Более того, авторами обнаруженного эффекта было сделано и заключение, что именно столь маломассивные прегалактические скопления газа и пыли наиболее благоприятны для интенсивного охлаждения в них газа, его результирующего уплотнения и запуска в нём очень активного (в 3-5 раз интенсивнее, чем «обычно», в видимых длинах волн) звездообразования, но никаких предположений о природе обнаруженного феномена и попыток его объяснения авторами (их - 74 персоны!) предложено не было.
Обнаруженный в инфракрасных лучах и обследованный класс галактик, без огромного избытка тёмной материи очень энергично «делающих» звёзды и находящихся на экстремально далёких расстояниях от нас, - это класс очень молодых сгустков обычной материи (чем дальше объект, тем он «моложе» на световом конусе). Из этого с несомненностью следует вывод, что отсутствие избытка тёмной материи в них и их молодость подтверждают мою гипотезу о том, что тёмная материя - это холодные плотные останки обычных звёзд - нейтронные и нейтронно-кварковые звёзды1, ибо в молодых сгустках газа (и пыли), в которых очень интенсивно и впервые образуются обычные (светящие) звёзды, таких плотных тёмных образований ещё не накопилось в столь большом количестве, чтобы влиять на гравитационную динамику обычной материи в сгустке, в частности, и на сам процесс и интенсивность звездообразования, как это происходит в «пожилых» галактиках, «захламлённых» холодными останками давно умерших и продолжающих умирать в них обычных звёзд.
ЛИТЕРАТУРА
1. Амбарцумян В.А. К вопросу о динамике открытых скоплений // Ученые записки Ленинградского гос. университета.
1938. № 22. Серия математических наук (астрономия). Вып. 4. С. 19-22.
2. Амбарцумян В.А. О протозвёздах // Доклады АН Арм.ССР. 1953. Т.16. №4. С. 97-102.
3. Григорян С.С. Гипотеза С.С. Григоряна о природе «тёмной материи» в галактиках и новые данные астрономических
наблюдений, подтверждающие её // Ломоносовские чтения. Тезисы докладов научной конференции. Секция механики. Москва, 16-25 апр. 2012, МГУ имени М.В. Ломоносова. М.: Изд-во Моск. ун-та, 2012. С. 57-58.
4. Григорян С.С. О динамике и строении Вселенной // Проблемы современной механики. Ч. 1. Сб. научных статей / Под
ред. акад. Л.И. Седова. М.: Изд. МГУ, 1983. С. 38-43.
5. Григорян С.С. О динамике начальной стадии соударения тел с большими скоростями // Некоторые вопросы механики
сплошной среды. Сб. научных статей / Под ред. С.С. Григоряна. М.: Изд-во МГУ, 1978. С. 157-172.
6. Григорян С.С. О некоторых специальных вопросах термодинамики сплошных сред // Прикладная математика и меха-
ника. 1960. Т. 24. Вып. 4. С. 651-662.
7. Григорян С.С. О природе наблюдаемых во Вселенной «чёрных дыр» // Ломоносовские чтения. Тезисы докладов научной
конференции. Секция механики. Москва, апр. 2005 г., МГУ имени М.В. Ломоносова. М.: Изд-во Моск. ун-та, 2005. С. 82.
8. Григорян С.С. О природе наблюдаемых во Вселенной «чёрных дыр» и динамике Вселенной // Известия Национальной
Академии наук Армении. Серия Механика. 2007. Т. 60. № 4. С. 3-19.
9. Григорян С.С. О природе «тёмной материи» в галактиках // Доклады Национальной Академии Наук Армении. 2012.
Т. 112. № 3. С. 273-276.
10. Сисакян А.Н. Избранные лекции по физике частиц. Дубна: ОИЯИ, 2004. 377 с.
11. Amblard A., Cooray A., Blain A., Bock J., Dowell C.D., Levenson L., Lu N., Nguyen H.T., Schulz B., Shupe D.L., Vieira J.D.,
Xu C.K., Zemcov M. "Submillimetre Galaxies Reside in Dark Matter Haloes with Masses Greater than 3*10u Solar Masses." Nature 470.7335 (2011): 510-512.
12. ESA. "Hershel Finds Less Dark Matter But More Stars (From ESA News Release, 16 Febr. 2011)." Space Research Today:
COSPAR's Information Bulletin. 180 2011: 51-52.
13. Friedmann A.A. "Über die Krümmung des Raumes." Zeitschrift für Physik 10.1 (1922): 377-386.
14. Friedmann A.A. "Über die Möglichkeit einer Welt mit konstanter negativer Krümmung des Raumes." Zeitschrift für Physik 21.1
(1924): 326-332.
15. Hawking S.W., Penrose R. "The Singularities of Gravitational Collapse and Cosmology." Proceedings of the Royal Society A:
Mathematical, Physical and Engineering Sciences 314.1519 (1970): 529-548.
16. Hawking S.W., Penrose R. "Black Holes in General Relativity." Communications in Mathematical Physics 25.2 (1972): 152-166.
17. Hawking S.W., Penrose R. A Brief History of Time from the Big Bang to Black Holes. New York: Bantam Dell Publishing
Group, 1988. 256 p.
18. Oppenheimer J. R., Volkoff G.M. "On Massive Neutron Cores." Physical Review 55 (1939): 374-381.
19. Oppenheimer J. R., Snyder H. "On Continued Gravitational Contraction." Physical Review 56 (1939): 455-459.
20. Glendenning N.K., Weber F. "Nuclear Solid Crust on Rotating Strange Quark Stars." Astrophysical Journal 100.2 (1992): 647-658.
21. Negreiros R.P., Mishustin I.N., Schramm S., Weber F. "Properties of Bare Strange Stars Associated with Surface Electric
Fields." Physical Review D82.10 (2010):103010-103016.
Цитирование по ГОСТ Р 7.0.11—2011:
Григорян, С. С. О природе «чёрных дыр», «тёмной материи» и динамике Вселенной / С.С. Григорян // Пространство и Время. — 2015. — № 3(21). — С. 92—102. Стационарный сетевой адрес: 2226-7271provr_st3-21.2015.26.
1 Григорян С.С. О природе наблюдаемых во Вселенной «чёрных дыр»; Он же. О природе наблюдаемых во Вселенной «чёрных дыр» и динамике Вселенной; Он же. Гипотеза С.С. Григоряна...