ОБЗОРЫ
АСТРОНОМИЯ, АСТРОФИЗИКА И КОСМОЛОГИЯ
Общая модель образования небесных тел от начальной конденсации газопылевых частиц до формирования «зародышей» планет
Т. Р. Абдульмянова
Казанский государственный энергетический университет, кафедра инженерной кибернетики. Россия, 420066, Казань, ул. Красносельская, д. 51.
Поступила в редакцию 11.03.2019, после доработки 23.04.2019, принята к публикации 29.04.2019.
В настоящей работе анализируются основные идеи формирования небесных тел в аккреционных и протопланетных дисках. Рассматриваются механизмы радиальной и орбитальной миграции небесных тел. Согласно результатам моделирования движения малых тел Солнечной системы формирование малых тел и рост их размеров и масс происходит в устойчивых точках либрации орбитальных резонансов. Появление единственного «зародыша» в протопланетных кольцах объясняется результатом действия механизма орбитального резонанса 1/1. При помощи модели волновых возмущений пылевых оболочек протозвезд и модели орбитального резонанса построена общая модель формирования и роста «зародышей» планет для дисков, имеющих массы менее 0.1М^.
Ключевые слова: протопланетный диск, протопланетные кольца, орбитальный резонанс, центры либрации.
УДК: 52-5. РЛСБ: 96.12.Bc.
ВВЕДЕНИЕ
Газопылевые диски вокруг молодых звезд впервые были обнаружены в 1984 г. при помощи спутниковых наблюдений в инфракрасном диапазоне (IRAS). Результаты этих и более поздних наблюдений показывают, что крупномасштабные структуры газопылевых дисков подразделяются на отдельные группы: некоторые диски имеют зоны низкой плотности круговой формы (круговые люки) [19, 20]. Некоторые другие имеют спиралевидную структуру или содержат пустоты определенных форм [16, 23, 35, 38, 50, 58]. Появление круговых люков объясняется тем, что такие пустоты появляются как зона формирования и роста планетезималей. Возникновение локальных пустот (холлов) объясняется [40] гравитационным влиянием формирующихся больших планет. Газопылевые диски вокруг звезд TW Hya, HD 163296, HD 16 9142, HD 97048 [23, 24, 52, 57] имеют кольцеобразную фрагментацию. Формирование такой структуры дисков объясняется устойчивостью равновесия газопылевых оболочек протозвезд [9]. Газопылевые диски вокруг звезд MWC 758, HD 100453, HD 135344B [23, 24, 41] имеют спиралевидную структуру. Появление таких структур у дисков может быть следствием изменения условия устойчивости и смены режима равновесия оболочки на режим сжатия протозвезды.
В 2008 г. астрофизики института Макса Планка обнаружили в окрестности звезды TW Гидры b пылевой диск с формирующейся планетой. В 2013 г. изображение этого пылевого диска было получено при помощи космического телескопа «Хаббл». Анализ результатов наблюдений газопылевых дисков показывает [21, 32-34, 53], что для интерпретации этих наблюдений необходимы математические модели процессов, происходящих в газопылевых дисках.
Динамика газопылевых сгущений (клампы) в массивных дисках, имеющих массу более 0.2МЯ, исследована в работе [55]. Показано, что на интервалах времени порядка 103 лет газовые клампы на звезду не падают и из них могут формироваться планеты как типа планет земной группы, так и подобные планетам-гигантам. В статьях [42, 47] исследована возможность радиальной миграции зародышей планет в «дебрис-дисках». Механизм формирования самих «дебрис-дисков» и тел метровых, десятиметровых размеров к настоящему времени до конца не выяснен [11].
В статьях [44-46] рассматривается возможность миграции зародышей планет земной группы из зоны планет-гигантов на ранней стадии эволюции Солнечной системы. В статье [22] показано, что в зоне планет земной группы условие локальной гравитационной неустойчивости не выполняется.
Результаты моделирования [9] показывают, что на расстоянии от 5 а. е. до 50 а. е. от центра протозвезды находится зона плохой обусловленности численного интегрирования уравнений газовой динамики. В таком случае моделирование эволюции газопылевых дисков с массой менее 0.2М^ может быть выполнено только при помощи численно-аналитических методов.
К настоящему времени найдены решения многих задач классической планетной космогонии. Однако до сегодняшнего дня не определены соответствия между стадиями эволюции протопланетных дисков и известными стадиями эволюции Солнца. Неизвестно также, росли ли планетезимали от пылинок с размерами 10-5 см до десятикилометровых тел путем слипания при взаимных столкновениях или на определенной стадии они прошли через стадию джинсовской неустойчивости в пылевом субдиске [11, 18].
а E-mail: abdulmyanov.tagir@yandex.ru
В настоящей работе рассматривается механизм формирования тел и роста их масс в дисках имеющих массы менее 0.2М^. В таких дисках на процесс формирования тел существенное влияние оказывает механизм волновых возмущений плотности газопылевой оболочки протозвезды и механизм орбитального резонанса [1, 2, 4-9]. Появление «зародышей» планет в протопланетных дисках может быть результатом действия механизма орбитального резонанса 1/1 [4], что подтверждается данными наблюдений астероидов вблизи орбиты Юпитера и вблизи орбит других планет Солнечной системы. Согласно модели, рассматриваемой в настоящей работе, центрами начального формирования тел и роста масс этих тел являются устойчивые точки либрации орбитальных резонансов.
1. МОДЕЛЬ ФОРМИРОВАНИЯ ЯДРА
ПРОТОЗВЕЗДЫ: УСЛОВИЕ ГАЗОДИНАМИЧЕСКОГО РАВНОВЕСИЯ ОБОЛОЧКИ ПРОТОЗВЕЗДЫ
Рассмотрим газопылевое облако, масса которого равна массе Солнца Мя, с соотношением газа — 99% и пыли — 1%. Такое облако будет гравитационно неустойчивым [56], если будет иметь радиус 2 • 106 Дя, среднюю плотность р0 = 10-19 г/см3, температуру Т = 15 К и будет сжиматься. За время сжатия (свободное падение) ^ = 0.2 млн лет радиус облака уменьшится в 100 раз. Выберем прямоугольную систему координат (ж, у, г), вращающуюся с невозмущенной (кеплеровой) угловой скоростью И так, чтобы ось г совпадала с осью вращения облака. Начало координат в центре протозвездного облака. Вращающееся облако представим как объединение тонких осесимметричных дисков, в которых движение межзвездного газа и пыли происходит без внутреннего трения (начальный этап формирования протозвезд). Пылевая оболочка протозвезды будет находиться в состоянии газодинамического равновесия, если скорость V = «хех + ууеу движения ее частиц, плотность р, давление р и потенциальная функция Ф будут представлять собой решение известных уравнений гидродинамики [11, 51]:
д
——+ V • grad р = —р V,
дV
Ж + ^ V =
= — 1УР — УФ — 2(П х V) + П2(жеж + уеу), У2Ф = 4пСр
(1) (2) (3)
в уравнения (1)-(3) и учитывая, что р0, Р0, Ф0, Vo — решение системы (1)-(3), с точностью до величин второго порядка получим следующие уравнения для малых возмущений [51]:
+ VoVрl + р^0 = —У1Ур0 — (4)
£ + (^У) • V! + • Vo =
—1— У(р0 + Р1) + — Ур0— (5) Р0 + Р1 Р0
У2Ф1
— УФ1 — 2(П х V1), 4пСр1.
(6)
Из уравнения (2) следует, что равновесие возможно только в том случае, когда р0, Р0, Ф0 являются решением следующего уравнения:
—1 Ур0 — УФ0 + П2(жвж + уву) = 0.
Р0
Вычисляя дивергенцию правой и левой части этого уравнения и учитывая уравнение Пуассона (3), получим окончательный вид уравнения газодинамического равновесия пылевой оболочки протозвезды:
У2Ф0 = —V у Ур^ + 2П2 = 4пСр0. (7)
Для межзвездного газа выполняется уравнение состояния идеального газа: р = (к/р)рТ. Откуда получим (к/р)Т = р0/р0. Учитывая, что во вращающейся системе координат невозмущенная скорость Vo равна нулю, вместо уравнений (4)-(6) получим следующие уравнения:
дР1 = —■У1Ур0 — (8)
дV1 1 _ 1 ^
=--,—у(р0 + Р0 + — уР0—
д£ Р0 + Р1 Р0 (9)
— УФ1 — 2(П х V1), У2Ф1 = 4пСр1. (10)
Вблизи устойчивого положения равновесия малым изменениям координат будут соответствовать малые изменения плотности. Тогда выражение VlVрo в уравнении (8) будет величиной второго порядка малости. Кроме того, с точностью до величин второго порядка получим
-1-У(р0 + р1) + — Ур0 =
Р0 + Р1 Р0
и если будут выполнены следующие условия: = — — (1 — — + .. Л У(р0+р1)+ — Ур0 « — — Ур1.
Р0 \ Р0 / Р0 Р0
дУ _ др _ дР _ дФ _ 0
Ж = Ж = Ж = Ж = .
В уравнении (2) через ех, еу обозначены векторы единичной длины, направленные вдоль оси ж и у. Обозначим через р0, р0, Ф0, V0 — равновесное решение системы (1)-(3), а через р1, р1, Ф1, V1 — малые возмущения равновесного решения. Подставляя р = Р0 + р1, р = р0 + р1, Ф = Ф0 + Ф1, V = Vo + V
В результате вместо уравнения (9) получим следующее:
^ = — — Ур1 — УФ1 — 2(П х Vl). (11) д£ р0
Вычисляя дивергенцию правой и левой части равенства (11), получим
- (УУ,) = -= дг( 1) дгр0вг
= - — У2Р1 - 4пОрх - 2 аып х VI). (12) Р0
В равновесии пылевая оболочка протозвезды будет находиться до тех пор, пока не произойдет уменьшение ее температуры за счет инфракрасного излучения. При уменьшении температуры равновесие оболочки нарушается и происходит сжатие протозвезды, в течение которого температура пылевого слоя существенно не изменяется. Поэтому предположим, что в равновесии пылевая оболочка будет находиться при постоянной температуре. В этом случае давление будет пропорциональным плотности: р\ = с2рь где с — скорость звука для газа (водород). Подставляя это выражение для плотности р1 в уравнение (12), получим уравнение для возмущений плотности пылевой оболочки:
д2р1
= с2У2р1 + 4пСр0Р1 + 2р0&у(П х Vl). (13)
дг2
Рассмотрим сначала уравнение (13) без учета вращения протозвездного облака:
д2р1
дг2
С2У2Р1 + 4пСр0Р1,
(14)
затем найдем решение уравнения (13), учитывающее вращение протозвездного облака.
1.1. Волновые возмущения плотности оболочки протозвезды
Уравнение (14) является линейным дифференциальным уравнением. Тогда его аналитическое решение можно найти при помощи метода разделения переменных. В цилиндрической системе координат х = г еов(у), у = г вт(у), г = г уравнение (14) имеет следующий вид:
д2Р1 = с2 (д2р1 +1 др1 +1 д2Р1 + д2р1 ^ +
дг2 у дг2 г дг г2 ду2 дг2 у
+ 4пСр0Р1. (15)
Согласно методу разделения переменных решение Р1(г, у, г) будем искать в виде произведения двух новых искомых функций Т и V: р1 (г, у, г, г) = Т(г) х х V(г, у, г). Подставляя это выражение для функции Р1 в уравнение (15), получим следующее уравнение для функций Т и V:
удт = /Тд?у + Т1 ду + Т &Г + Тд!Х V
с2 дг2 у дг2 г дг г2 ду2 дг2 у
+ 4пСр0ТУ/с2. (16)
Уравнение (16) можно заменить равносильной ему системой из двух более простых уравнений. Для этого необходимо правую и левую части уравнения (16) разделить на произведение функций V и Т:
1 д2Т (1 1 ^ 1 1 д^А
+ Т7
+
+
+
с2Т дг2 VV дг2 ' Vг дг ' Vг2 ду2 ' V дг2 )
+ 4пОр0/с2 = -А2. (17)
Левая часть уравнения (17) является функцией только переменной г, а правая — зависит только от г, у и г. В этом случае как правая, так и левая
части равенства (17) будут равны одной и той же константе. Обозначим эту константу через — А2. В результате вместо уравнения (17) получим следующие два уравнения:
ж + А2А2Т
0,
1 д^ 1 дV V дг2 Vг дг Vг2 дф2
1 д 2V 1 д ^А
+
V дг2
+4пОр0/с2 = -А2
+
(18)
(19)
Уравнение (18) имеет простое решение. Поэтому рассмотрим только уравнение (19). Решение уравнения (19) будем искать также методом разделения переменных, представив искомую функцию V в виде произведения функций Ш и Z: V(г, у, г) = Ш(г, у)^(г). Подставляя это выражение для функции V в уравнение (19), получим следующее уравнение:
д2Ш Z дШ Z д2Ш ттrд2Z\ Z--1----1----+ Ш- =
дг2 г дг г2 ду2 дг2 у
= -[4пСр0/с2 + А2]Ш • Z. (20)
Разделяя переменные Ш и X, вместо уравнения (20) получим следующие два уравнения:
d2Z 27
0,
(21)
д2Ш 1 дШ 1 д2Ш г. „ /2 ,2 2т, п - + —— + "2 1^ + [4пСр0/с2 + А2+72]Ш = 0.
дг2
дг
ду'
2
(22)
Уравнение (21) имеет простое решение. Для решения уравнения (22) также воспользуемся методом разделения переменных. Представим функцию Ш в виде произведения двух новых функций: Ш = К(г) • Б (у). Подставляя это выражение для функции Ш в уравнение (22), получим:
„д2К Б дК Кд2 Б 2 2 2
Б—^ I---т,—I—о +[4пСр0/с2 + А2 + 72]К• Б = 0.
ог2 г ог г2 оу2
(23)
После разделения переменных К, Б вместо уравнения (23) получим следующие два уравнения:
¿у2
^ + v2S = 0,
¿2 К 1 ¿К
¿г2 г ¿г
{4пОр0/с2 + А2 + 72} - ^
(24)
К = 0. (25)
Уравнение (24) имеет простое решение, а уравнение (25) можно представить в виде уравнения Бесселя. Обозначим через х выражение [4пСр0/с2 + + А2 + 72]1/2г. Тогда уравнение (25) примет классическую форму уравнения Бесселя:
¿2 К 1 ¿К
¿х2
+ --Г + I 1 - - )К = 0.
¿х
(26)
Частными решениями уравнения (26) являются функции Бесселя порядка V:
К(х) = (х) = (г^4пОр0/с2 + А2 + 72).
2
Величину константы Л определим так, чтобы выполнялось граничное условие р1(Й0, у, г, 4) = 0:
Й0^4пСр0/с2 + Л2 + 72 = Лк„.
Откуда получим константу Л:
Л2 = ЛЪ /Д — 4пСр0/е2 — 72.
Согласно решению уравнения (18) равновесие пылевой оболочки будет устойчивым, если будет выполнено следующее условие:
At/Д2 - 4nGpo/c2 - Y2 > 0.
(27)
Произвольная константа y в уравнении (21) равна: Y = [1/(zi - zo)] • ln{p(zo)/p(zi)}, где zo, zi — начальный и конечный радиусы протозвездного облака, а константа v в уравнении (25) является порядком функции Бесселя. Учитывая это, для случая устойчивого равновесия получим следующие частные решения уравнения (15):
pi(r, у, z, t) = (t) • Rkv(r) • Sv(y) • Z7(z),
где Tkv (t) = akv cos(cAt) + bkv sin(cAt), Z7 (z) = = p(zi) x exp(-Yz), Sv(у) = cos(vy + yo), Rkv(r) = = Jv(r), Jv — функция Бесселя порядка v. Коэффициенты akv, bkv определяются при помощи начальных условий краевой задачи. Уравнение (i5) является линейным, поэтому сумма полученных частных решений также будет решением уравнения (i5). Эта сумма будет общим решением уравнения (i5):
pi(r,у, z,t) =
СО СО
= p(zi) exp(-Yz) Y Sv (у) Y Tkv (t)Rkv (r). (28)
v=0 k=0
Если условие (27) не будет выполняться, то равновесие будет неустойчивым:
A2 = 4nGpo/c2 - A\v/R0 + Y2 > 0. (29)
В этом случае в общем решении (28) функция Tkv(t) будет другой:
Tkv(t) = akv • exp[cAt] + • exp[-cAt].
Начальные и граничные условия краевой задачи для волнового уравнения (i5) выбираются исходя из величины начальной плотности p(zo) = = p0 = i0~i9 г/см3 протозвездного облака. Практически нулевая начальная плотность p(zo) дает основание для выбора модели краевой задачи с нулевыми начальными и граничными условиями:
pi (r, zo, 0) = ^ Е Jo(Akr/Ro), (r, zo, 0) = 0,
n + 1
k=0
|р1 (0, г, ¿)| < то, р1(Й0, г,¿) = 0. (30)
В случае почти нулевой начальной плотности и однородности протозвездного облака индекс V в функции Б1>(у) будет равен нулю. Учитывая начальные и граничные условия (30), получим, что коэффициенты ак
будут равны ак = А = р(^) • ехр^0)/[р(г1) • (п + 1)], Ьк = 0. В этом случае решение краевой задачи будет следующее:
р1 (г, г, 4) = р(^1) ехр(—■7^) • Ах
• оо® (Лк/йс)2 — 4пСр0/е2 — 7Л х к=0 ^ '
х 70(Лкг/Й0). (31)
Рассмотрим теперь решение волнового уравнения (13) с учетом вращения протозвезды (П = 0). Дивергенцию векторного произведения П х V можно вычислить по следующей формуле [14]:
div(n х V) = V • rot П - П • rot V, (32)
где
rot V
dz
1 dVR
r дф
1 d(rVv
r dr
er +
AÖVr ÖVz
rot П = —
3 /GM
r5 eV>
\ dz
1 dVv
r дф
0,0,
dr
GM
ev+
h h V =(VR, К,0), Vr = -e • sin v, Vv = -, p R
p = a(1 - e2), - = \JGMp, R = p/(1 + e • cos v),
a, e — большая полуось и эксцентриситет кеплеров-ской орбиты, er, ev, ez — векторы единичной длины по направлениям r, z и по направлению касательной к орбите. Радиус-вектор R эллиптической орбиты и полярный радиус r будут совпадать только при e = 0. Также будут отличаться полярный угол у и истинная аномалия эллиптической орбиты v. Вычисляя частные производные и подставляя в формулу дивергенции (32), получим:
3 GM • (1 + e • cos v)
divíO x V) = --, -1 к,
2 yja( 1 - e2)r5/2
3 GM • (1 + e • cos v)
2
= f (r, v).
Решение с учетом вращения облака будем искать в виде суммы общего решения р1(г, г, 4) уравнения (14) и частного решения р^(г, 4) уравнения (13). Для этого функции / (г, V) и р^(г, 4) разложим в ряд Фурье по функциям Бесселя 70(Лкг/До):
1 Г п
/(г) = 2" /(г, = £ дк70(Лкг/Д0),
^ _Л
k=0
Pd(r, t) = Y dk(t)J0(Akr/R>),
1
k=0 R 2n
g0
gk =
n[A0Ji (A0)]2 1
n[AkJi (Ak)]2
J Jf (r, v) • r • J0(A0r/ñ0)dvdr,
Z1 0
R 2n
/// (r, v)t^(A, r/R0)dvdr,
z 1 0
k > 0.
z
3
r
Подставляя функцию р1 (г, г, г) + р(1 (г, г) в уравнение (13) и учитывая, что р1 является решением уравнения (14), для нахождения функций ¿к(г) получим следующие дифференциальные уравнения второго порядка:
d2 dk (t) dt2
4nGpo - c2
Xk Ro
dk (t) + 2pogk. (33)
Частными решениями уравнений (33) являются следующие функции:
dk (t) = Ak cos (в1/21) + 2pogk/в, в = (cXk/Ro)2 - 4nGpo,
где ck = -p(zo)/[(n+1)(1 +2pogk/в)]• Следовательно, решение уравнения (13), учитывающего вращение протозвездного облака, будет иметь следующий вид:
pi (r, z, t) + pd(r, t) = p(zi) exp(-yz) • Ax
X ]T cos (tc^(Xk/Ro)2 - 4nGpo/c2 - x
k=0 ^ ' n
X Jo (Xk r/Ro) + ^ dk (t) Jo (Xk r/Ro). (34) k=o
Коэффициенты ck определяются из условия p1 (r, Ro,o) + pd(r,o) = o краевой задачи для уравнения (13).
Рассмотрим графики функций p1 (r, z, t), pd(r, t) и выясним влияние вращения протозвездного облака на плотность газопылевой оболочки прото-звезды. На рис. 1, а изображены графики функций p1 (r, Ro,o) и pd(r, t) для параметра Ro = 5 a.e. Согласно рис. 1, а максимумы (минимумы) плотности p1 расположены на расстояниях 2 • Ro, кратных значению параметра Ro = 5 a. e. Такая особенность расположения максимумов характерна для любых значений параметра Ro и связана со свойствами бесселевых функций. В зонах максимальной плотности p1 пылевые частицы будут аккумулироваться быстрее, чем в окружающих их зонах. Аккумуляция пылевых частиц будет ускорять их оседание на экваториальную плоскость. При оседании этих частиц на экваториальную плоскость особенности в распределении частиц пылевой оболочки протозвезды будут проявляться также и на экваториальной плоскости протозвездного облака. На рис. 1, б представлены графики функций p1 (r, Ro,o) и pd(r, o) для значения параметра Ro = o. 18 а. е.
Сравнивая эти графики (рис. 1), заметим, что вариации плотности pd(r,o) составляют около 2% от вариаций p1 (r, Ro, o) для Ro = 5 а. е. (рис. 1, а), а для Ro = o. 18 а. е. — 12o% (рис. 1, б). Для значения параметра Ro = 1.45 а. е. — 6%.
Таким образом, при больших радиусах Ro ядра протозвезды вращение протозвездного облака не оказывает существенного влияния на плотность p1 (r, Ro,o). При Ro = o. 18 а. е. величина pd(r, o) намного превосходит величину плотности p1 (r, Ro,o). Однако это влияние не будет распространяться на вариации плотности p1 (r, Ro,o), на расположение максимумов плотности и на перемещение волн плотности (при t = o).
0.10 0.05-0.05 -0.10-1 -0.15
10
15
20 25
30
35 40 r, а.е.
0.30
0.25
о 0.20
а 0.15
а
0.10
0.05
2
3
4
r, а.е.
Рис. 1. Профили плотности р1 (г, Я0, Ь) + ра(г, V, Ь) для значений параметра Я0: а — для Я0 = 5 а. е. на интервале 1 ^ г ^ 40 а.е.; б — для Я0 = 0.18 а.е. на интервале 1 ^ г ^ 5 а.е. (е = 0.3)
а 0.10
+
а 0.05 0
-0.05 -0.10
10
2
4
6
8
r, а.е.
10
Рис. 2. Профили плотности р1 (г, Я0, Ь) + ра(г, V, Ь) для значения параметра Яо = 1.45 а.е.: а — для Ь = 0; 1.5; 3; б — для Ь = 4.5; 6; 7.5 (рад)
На рис. 2 изображены профили плотности р1 (г, у, К0, г) + р,1 (г, V, г) пылевой оболочки для значений параметров К0 = 1.45 а.е., для значений г = 0; 1.5; 3 (рис. 2, а) и для значений г = 4.5; 6; 7.5 (рис. 2, б), V = п/2, у = 0. Профили плотности р1 (г, у, К0, г) + р(1 (г, V, г), изображенные на рис. 2, показывают квазипериодическое движение плотност-ных волн для значения параметра К0 = 1.45 а. е.
2
5
б
5
0.15 0.1 0.05 0.0
-0.05 -0.1 -0.15 -I
0
30
40
40
0.008 0.006
0.004
0.002 0.0 -0.002 -0.004 -0.006 -0.008
-0.008
I " '' 111
-0.004
0.0
11 ■111 ■ ■ ■ ■ I ■ ■ ■ ■ I ■
0.004 0.008
Рис. 3. Плотность дисков и орбиты миграции частиц и малых тел: а — 3Э-график плотности р1 (г, Я0,0) + ра(г,0) для параметра Я0 = 5 а. е. и полярных расстояний 0.18 ^ г ^ 40 (а. е.); б — либрационные орбиты малых небесных тел вблизи орбитального резонанса 1/1 для резонансных параметров а2 = 0.6; 0.8; 1.7; 2.2 и наклонов орбит г = 25°,
г* = г — 0.817 (проекции на плоскость орбиты зародыша планеты)
б
Скорость бегущей волны, согласно решению (34), приближенно равна сАк/Я0 = Ак/(Я0 )11/8, где с = 1/(^0)3/8. Тогда период возвращения бегущей волны для параметра Я0 будет приближенно равен 2п(Я0)11/8/А1 = 4.36 рад. Принимая скорость звука равным с = 1286 м/с (для водорода), получим оценку периода времени возвращения бегущей волны в исходное положение: Р = 115.36/Ак лет. Наименьший период имеет волна для наибольшего значения Ак (2 < Ак < 20). Следовательно, для периода Р бегущей волны получим следующую оценку: 5 < Р < 60 лет.
На рис. 3, а представлены 3Э-профили плотности р1 (г, Я0, £) + ра(г, V, £), вычисленные для £ = 0, Я0 = 5 а. е. Полярная плоскость (г, преобразована в декартову плоскость (х, у), г = (х2 + у2)1/2, 0.5 ^ х < 40, 0.5 ^ у ^ 40 (а. е.), 0 ^ ^ ^ п/2.
1.2. Условие устойчивости равновесия оболочек протозвезд и границы зоны гравитационной неустойчивости
В общем решении (31) предполагалось, что движение частиц газопылевой оболочки протозвезд происходит без внутреннего трения. Однако такое предположение справедливо только для газовой составляющей оболочки. В газопылевой части оболочки, особенно в экваториальном диске протозвезды, уже будет влиять кинематическая вязкость. Для учета влияния кинематической вязкости необходимо уравнение (2) заменить уравнением Навье-Стокса [62]:
^ = -1VP - УФ - 2(П х V) + О2(хеж + уеу)+
+ V
V2 V +1 У(а1у V) 3
где V — кинематическая вязкость. Динамика газопылевых частиц экваториальных дисков (аккреционных дисков, «дебрис-дисков», протопланетных дисков) в настоящее время исследуется при помощи численного интегрирования уравнений газовой динамики [33, 37, 39, 49].
В работе Снытникова [55] численными методами исследована динамика газопылевых сгущений (клам-пы) для дисков массы более 0.2М8. Показано, что формирующиеся клампы мигрируют в радиальном направлении к ядру протозвезды по сложным орбитам. Было обнаружено также, что на интервалах интегрирования порядка 103 лет клампы не аккре-цируют на ядро протозвезды. Однако не выяснена причина, по которой клампы не будут аккрецировать на более длительных интервалах интегрирования.
В статьях Морбидели [42] и Несворного [47] исследована возможность радиальной миграции зародышей планет в «дебрис-дисках». Механизм образования самих «дебрис-дисков», как и тел метровых и десятиметровых размеров, к настоящему времени не выяснен [11]. Следовательно, полученная в этих статьях миграция зародышей планет имеет гипотетический характер,
В статьях Наякшина [44-46] рассматривалась возможность миграции зародышей планет земной группы из зоны планет-гигантов. В зоне планет земной группы условие локальной гравитационной неустойчивости не выполняется [22]. На основании этого предположено, что зародыши планет земной группы могли мигрировать из зоны планет-гигантов. В монографии [54] было показано выполнение условия гравитационной неустойчивости в зоне планет-гигантов.
В работе [9] показано, что газопылевые оболочки и их экваториальные диски неустойчивы для всех масс протозвездного облака только на расстояниях более 50 а. е. от центра протозвезды. На расстояниях от 5 а. е. до 50 а. е. почти все экваториальные диски, за исключением сверхмассивных, могут быть как гравитационно неустойчивыми, так и устойчивыми. Эта зона расположена близко к границе устойчивости [9]. В этом случае задача численного интегрирования уравнений газовой динамики становится плохо обусловленной задачей.
Локальная гравитационная неустойчивость в дисках малой массы (< 0.2М8) может наступить в ре-
зультате действия волновых механизмов. Плотност-ные волны, согласно решению (34) краевой задачи для волнового уравнения (13), могут обеспечить локальное выполнение условия гравитационной неустойчивости в результате фрагментации экваториальных дисков на протопланетные кольца [9]. Механизм плотностных волн является вторичным по отношению к механизму гравитационного сжатия протозвездного облака, так как является следствием аккреции межзвездного газа и газовых уплотнений на ядро протозвезды. Эти два механизма действуют независимо друг от друга, несмотря на то, что они взаимосвязаны. Признаком, по которому можно определить завершение или временную остановку процесса аккреции, является выполнение условия гравитационной устойчивости (27).
Интервал полярных расстояний, меньших 5 а. е., находится (за исключением сверхмассивных межзвездных газовых облаков) в зоне устойчивости [9]. Согласно статьям [44-46] на таких расстояниях невозможно формирование зародышей планет и зародыши должны мигрировать из внешней зоны. Однако радиальная миграция по какой-то причине должна быть остановлена. Причиной прекращения радиальной миграции может быть установление устойчивого равновесия и завершение процесса аккреции или временная его остановка (разд. 1.1). Согласно решению (34) в зоне устойчивости экваториального диска формирование зародышей планет возможно за счет действия механизма плотностных волн и кольцевой фрагментации этой зоны. Следовательно, зародыши планет земной группы могли не мигрировать из зоны планет-гигантов, а могли формироваться в зоне планет земной группы. В этом состоит отличие выводов Наякшина от выводов, полученных в настоящей работе.
Таким образом, в экваториальных дисках малой массы (< 0.2Ы,Э), согласно модели (34), зародыши планет могут формироваться в условиях устойчивости газопылевых оболочек протозвезд. Учитывая, что экваториальный диск является частью газопылевой оболочки протозвезды, нетрудно заметить противоположность локальной неустойчивости экваториальных дисков и устойчивости оболочки про-тозвезды. Однако волновой механизм и механизм гравитационного сжатия протозвездного облака действуют независимо друг от друга. Следовательно, такая противоположность не является противоречивой. В дальнейшем будем исходить из предположения о кольцевой фрагментации оболочек протозвезд и их экваториальных дисков. Для выбора модели формирования малых тел в экваториальных дисках далее предварительно рассматриваются результаты наблюдений и исследований динамики малых тел Солнечной системы [1, 2, 15, 28-31, 43].
2. РЕЗОНАНСНАЯ МОДЕЛЬ ОРБИТАЛЬНОЙ МИГРАЦИИ МАЛЫХ ТЕЛ
В статье [4] рассматривалась динамика малых тел Солнечной системы в модели плоской круговой
ограниченной задачи трех тел (случай плоских экваториальных дисков). Для плоских экваториальных дисков была определена плотность распределения
газопылевых частиц в кольце, при достижении которого начинает действовать механизм орбитального резонанса и формирование зародыша планеты. Однако экваториальные диски протозвезд не являются плоскими. Орбиты частиц и орбиты малых тел в таких дисках будут иметь ненулевые наклоны. Это обстоятельство необходимо учитывать в соответствующей модели. Поэтому далее рассмотрим проблему формирования и динамики малых тел в пространственной модели ограниченной задачи трех тел.
Рассмотрим движение астероида под действием гравитационных сил со стороны Солнца и одной из больших планет. Массу астероида будем считать бесконечно малой по сравнению с массами больших планет. В качестве основной системы возьмем прямоугольную гелиоцентрическую систему координат, в которой орбита большой планеты является круговой. Основной плоскостью выбранной системы координат будем считать плоскость орбиты планеты. Единицей измерения масс тел — сумму масс Солнца и планеты, а единицей измерения длины — величину большой полуоси планеты. Предположим также, что единица времени выбрана так, что постоянная Гаусса равна единице. Движение астероида вблизи соизмеримости (р + q)/p средних движений астероида и планеты означает выполнение условия: |п - (р + (1)п1/р)1 ^ 0(п1л/ш), где т — масса планеты в выбранной системе единиц, п и П1 — средние движения астероида и планеты. Канонические уравнения круговой ограниченной задачи трех тел в системе координат, вращающейся с угловой скоростью п1, имеют следующий вид [1, 2]:
¿С _ЗГ ¿Г _ дГ ¿(Н - С) _ дГ = Ж', = ~31, ¿г = Ж,
¿К _ дГ ¿А' _ дГ ¿1 _ дГ
Иг = Ж, Ж = - дС, ¿г = - дГ,
dh ~dt
F =
dF
d(H - G)'
dk lit
dF dK,
1
2(G + Г)2
+ n1G - П1К + mR,
R = (1 + r2 - 2r • cos N)-1/2 - r • cos N, cos N = c • cos(t - T1) + s • cos(t + T1 - 2П), F = Fo + F1 + F2 + F3 + O(me4), Fo = 1/[2(G + Г)2]+ mG + mfo(X'), fo(X') = (1 + r2 - 2roccos A')-1/2 - roAcos X',
(35)
L = va, G = L\J 1 - e2,
Г = L - G,
H = G • cos(i), K =
c = cos2(«/2), i < 4oo,
dF , dkdt,
где a, e, i, I, g, t, A, П — большая полуось, эксцентриситет, угол наклона орбиты, средняя аномалия, аргумент перигелия, истинная и средняя долгота астероида, восходящий узел, A1, т1 — средняя и истинная долгота планеты, ro — величина полярного радиуса r в центре либрационных движений астероида, A' = A - A1, k = n1t + ko. Предположим, что наибольшая величина эксцентриситета не будет превосходить o.3. Тогда составляющие Fo, F1 + F2, F3
0.008 0.006 0.004 0.002 0.0 -0.002 -0.004 -0.006 -0.008
100
0.008 -0.004 0.0..... 0.004.....0.008
150
200
250
300
350
Рис. 4. Либрационные орбиты малых тел вблизи резонанса 1/1 для значений параметра а2 = 0.6; 0.8; 1.7; 2.2: а — для наклона орбиты г = 25° в полярных координатах (г*, А'), г* = г — 0.817 (проекция орбит на плоскость орбиты планеты); б — распределение наиболее крупных нумерованных «греков» и «троянцев» (резонанс 1/1) в полярной системе координат
гамильтониана Р будут величинами порядка т, тоз/2, т2 соответственно. Для того, чтобы найти решение системы канонических уравнений (35) при помощи методов теории возмущений, необходимо сначала найти решение задачи для промежуточного гамильтониана р). Гамильтониан Р системы уравнений (35) является функцией канонических переменных {С, Г, Н — С, К; А', 1, Л., к}. Для определения значения С0 канонической переменной С в центре либрации необходимо ввести систему координат, вращающуюся с угловой скоростью (р + я)"-1/р, и вычислить частную производную р [1, 2]:
Ж = -(С + Г) +— Приравнивая к нулю эту частную производную и решая полученное уравнение относительно С, найдем С0:
Go
Р
1/3
- Го, ro = GO + O(m). (36)
(р + я)«1.
Угловые координаты центров либрации определяются как точки экстремумов усредненной функции /0(А') [2]. Решение системы канонических уравнений (19) для промежуточного гамильтониана р имеет следующий вид [1, 28-31]:
С = С0 + С2р1 + 2С0р2/3 + ..., Р1 = — 1/3V6m[a2 - /0(А')]1/2 sgn(dА'/dt),
t
V6mG0
4 G3(A, А' ) [a2 - /о(А')]1/2 - 9G3(A - Лl),
^ = V2Tcos(1 + 10), n = V2Tsin(1 + 10),
H — G = const, K = const,
h = —nit + h-0, k = nit + k0, (37)
где А1 — значение переменной А' = А — Ai = = А* + qA1/p на границе либрации, А* = А — — (p + qMi/p, ti — момент прохождения границы либрации, а2 — резонансный параметр, I0 — начальное значение средней аномалии астероида.
На рис. 4, а представлены либрационные орбиты резонанса 1/1 для значений резонансного параметра а2 = 0.6; 0.8; 1.7; 2.2 и наклона орбиты i = 25°
в полярных координатах (г*, А'), е = 0.2, г* = = г - 0.817. Полярный радиус г определяется по формуле г = р2 + 2С0р + С0, р = С - С0, каноническая переменная С — при помощи решения (37). На рис. 4, б показано распределение наиболее крупных нумерованных астероидов резонанса 1/1 («троянцы» и «греки» Юпитера), которые сконцентрированы в окрестности двух устойчивых точек Лагранжа: Ь4 (С0, п/3) и Ь5 (С0,5п/3). При помощи решения (37) канонических уравнений получим формулу для оценки ширины резонансной зоны:
Д :
22 Gmax — G0
—-Gq 1ас + O(m),
2 2 2 а = а — а
По этой формуле можно определить ширину резонансной зоны, за пределами которой расположена внешняя часть зоны резонанса и внешняя часть протопланетного кольца.
Согласно результатам современных наблюдений (рис. 4, б), вблизи устойчивых точек Лагранжа ¿4 и ¿5 располагаются десятки тысяч астероидов больших и малых размеров. При появлении пылевых частиц в таких точках произойдет их аккумуляция. По этой причине еще на ранних этапах эволюции протопланетного диска в таких точках могли формироваться малые тела до метровых и больших размеров. Если это так, то по какой причине не происходит дальнейшая интеграция «троянцев» и «греков» и вместо этих групп астероидов не образуются два больших астероида? По какой причине останавливается процесс роста масс тел? Для того, чтобы найти ответы на эти вопросы, рассмотрим более детально либрационные орбиты малых тел и пылевых частиц внутри резонансной зоны. Сравнивая либрационные орбиты, изображенные на рис. 4, а, заметим, что с увеличением резонансного параметра а2 и наклона орбиты г малые тела и пылевые частицы будут приближаться к той точке гравитационной системы Солнце—астероид—планета, в которой располагается планета. Однако в решении (37) резонансный параметр а2 является константой. Увеличение этого параметра может происходить только под влиянием
б
внешних механизмов. Таким механизмом может быть механизм роста масс малых тел. Оседание пылевых частиц на экваториальную плоскость протозвезды приведет к росту масс малых тел. В результате этого кеплеровы орбиты малых тел будут преобразованы в орбиты с вращающимися перигелиями [3]. Такое изменение орбит малых тел приводит к орбитальной миграции малых тел и их аккреции на зародыш планеты. Увеличение массы зародыша планеты приводит, согласно интегралу «энергии» (38), к росту величины постоянной «энергии» или резонансного параметра а2. Следовательно, практически весь приток пылевых частиц из оболочки протозвезды в экваториальную плоскость будет обеспечивать формирование и рост массы малых тел, их миграцию по либрационным орбитам (рис. 4, а). Рост массы «троянцев» и «греков» останавливается с окончанием притока пыли в зону резонанса 1/1. В это же время завершается и рост массы самого зародыша планеты и зародыш становится планетой.
2.1. Механизм формирования малых тел в экваториальных дисках протозвезд
Аккреционные диски от протопланетных дисков отличаются прежде всего сравнительно большой скоростью аккреции газа на ядро протозвезды. Если скорость аккреции будет оставаться постоянной и сравнительно большой до завершения аккреции, то диск в этом случае будет аккреционным. Однако результаты моделирования [9] показывают, что аккреционными будут диски только при сжатии сверхмассивных газопылевых облаков, имеющих массу более 103Ms. К таким случаям относятся формирование спиральных галактик и эволюция тесных двойных звезд [59]. В остальных случаях скорость аккреции будет уменьшаться. Экваториальный диск будет аккрециннным только до тех пор, пока радиус ядра протозвезды будет больше 50 а. е. С уменьшением радиуса ядра протозвезды (< 50 а. е.) скорость аккреции будет уменьшаться. При этом оболочка протозвезды будет находиться вблизи граничной зоны равновесия, где режим устойчивости может смениться на режим неустойчивости. Экваториальный диск в этой зоне становится протопланетным диском. На этом этапе в экваториальных дисках будут формироваться малые тела и «дебрис-диски». Так называемые «дебрис-диски» представляют собой ранний этап эволюции протопланетных дисков.
Формирование тел малых размеров и рост их масс в «дебрис-дисках» может происходить в устойчивых точках либрации орбитальных резонансов (p + q)/p [4]. Кроме того, в разд. 2 было замечено, что либрационные орбиты малых тел с ростом их масс трансформируются в орбиты аккреции на зародыш планеты. Выясним возможные изменения этих орбит при помощи первых интегралов ограниченной задачи трех тел. Для этого рассмотрим интегралы канонических уравнений (35) для промежуточного гамильтониана Fo [1]: Г = const, H — G = const, K = const. Полный дифференциал усредненного гамильтониана Fo будет равен нулю. Следовательно, вдоль решения системы канонических уравнений усредненный гамильтониан будет оставаться константой:
Fo
+ — mG + mf (X*) = h,
2(G + Г)2 ' p
2n
f(X*) = ^ f fo(X — Xi)dl, X* = X — X1, 2W p
(38)
где к — константа интегрирования. Интеграл (38) называется интегралом энергии [30]. Значение резонансного параметра а2 в решении (37) равно значению усредненной функции /о(А') в точках поворота прямого к обратному либрационному движению [30]. Усреднение в (38) выполнено по схеме Моисеева [27]. Выбор такой схемы связан с выбором соответствующей порождающей орбиты.
В работе Пуанкаре [48] была отмечена важность выбора порождающих орбит при определении топологии либрационных орбит. Полная классификация кеплеровых орбит в системе координат, вращающейся с угловой скоростью П1, проведена в монографии [15]. Показано, что порождающие орбиты в этом случае замыкаются через д оборотов вокруг начала координат. В точной соизмеримости эти орбиты представляют собой окружности, а в остальных случаях — лепестковые розетки. Синодический период таких орбит равен 2пр/щ. При возмущенном движении малых тел вблизи соизмеримостей средних движений траектории их движения будут представлять собой раскачивающиеся розетки.
Однако в системе координат, вращающейся с угловой скоростью [(р + д)/р]п1, эти же кеплеровы орбиты будут иными. В точной соизмеримости круговая орбита вырождается в точку. Все некруговые орбиты в этом случае замыкаются после одного оборота вокруг точки, в которую вырождается круговая орбита. Начало координат при этом остается во внешней стороне орбиты. Период таких орбит равен 2пр/[(р + д)п1]. При возмущенном движении малых тел вблизи соизмеримостей средних движений их траектории будут представлять собой орбиты либра-ционных движений гамильтоновой системы (35).
Полярный радиус г0 центра либрации определяется в формуле (36). Для того, чтобы определить полярный угол точки равновесия канонической системы, необходимо вычислить производную усредненной функции /о(Х') и приравнять эту производную к нулю. Решая полученное уравнение, определим угловые координаты центров либрации. На рис. 5 представлены графики усредненной функции /о(Х') для соизмеримостей 1/1, 4/3, 3/2 (рис. 5, а) и соиз-меримостей 2/1, 7/3, 3/1 (рис. 5, б). Для резонансов 1/1, 4/3, 3/2 (рис. 5, а) усредненные функции имеют два минимума и два максимума. В случае люковых резонансов 2/1 и 7/3 усредненные функции имеют три минимума и три максимума, а соизмеримость 3/1 — два минимума и два максимума (рис. 5, б). Либрационный тип движения будет вблизи всех устойчивых точек равновесия (минимумы) [2].
Сравнение результатов усреднения функции /0(Х') по схеме Моисеева, по схеме Фату и по времени показывает [10], что схема Моисеева является усреднением по наименьшему периоду, при котором либрационный тип движения будет для всех соиз-меримостей. Увеличение периода усреднения сгла-
Рис. 5. Графики функции /0(А'), усредненной по схеме Моисеева: а — для соизмеримостей 1/1, 4/3, 3/2; б — для
соизмеримостей 2/1, 7/3, 3/1
живает существующие особенности гамильтоновой системы [10]. Адекватность гамильтоновой системы, полученной при помощи усреднения по схеме Моисеева, можно проверить при помощи теоремы Боголюбова [13]. Несмотря на то, что теорема Боголюбова и ее е-оценка справедливы только для стандартных систем, при учете особенностей резонансных систем [10] эту теорему можно применить и к системам, усредненным по порождающей траектории. В других обоснованиях схем усреднения, кроме теоремы Боголюбова, нет необходимости. Следовательно, классификацию либрационных орбит внутри каждого в отдельности взятого резонанса можно проводить при помощи классических методов теории возмущения [36, 43]. Общая классификация либрационных орбит, принадлежащих разным резонансам, возможна только при помощи методов теории дискретных структур.
Подставляя в интеграл (38) выражение (36) для первой координаты С центра либрации, получим:
F0
1
2 ([p/((p + q)ni) - Го]1/3 +г„)2
г +
(p + q)
-ni
(p + q)ni_
1/3
- Го +
+ 2та2 = Ь. (39)
Учитывая связь между гамильтонианами во вращающейся и неподвижной системах координат [31], получим следующее выражение, связывающее постоянную Якоби С с резонансным параметром а2 для резонансов общего вида (р + я)/р:
C = 2Р0 + m(1 - m)
1
([p/((p + q)n1 ) - Го]1/3 + Г„)2 1/3
> +
, 2(p + q)
+--П1
p
p
- Г0 +
(р + я)п1_ + 2та2 + т(1 - т). (40)
Выражение (38) интеграла энергии было получено для случая орбитального резонанса, то есть в предположении, что линия апсид и линия узлов возмущаемого тела являются неподвижными. В случае их
вращения для аргумента перицентра и долготы узла орбиты возмущаемого тела в линейном приближении получим: ш = W0t + С1 и П = ^t + c2. Интеграл энергии (38) в этом случае будет иметь следующий вид:
Fo =
1
2(G + Г)2 ' p
p + q
+ -— (n1 + wo + Oo)G+
+ mf (A*) = h. (41)
Существование интеграла энергии (41) показывает, что орбитальные и вековые резонансы представляют собой единый механизм эволюции орбит малых небесных тел. Выполнение условия векового или орбитального резонанса может играть стабилизирующую роль в движении малых тел, если движение этих тел происходит вблизи точек минимума усредненного гамильтониана (рис. 5). В таких точках может происходить формирование тел начиная с малых до метровых и километровых размеров. Вблизи неустойчивых точек либраций движение становится стохастичным. Однако влияние такой стохастично-сти ограничено величиной вариации констант в интегралах канонических уравнений (35). Величина этих констант будет медленно меняться по мере роста масс малых тел и массы m зародыша планеты. При малых возмущениях движения малых тел, эти тела будут совершать либрационные движения вблизи устойчивых точек либраций [2]. При больших возмущениях малые тела аккрецируют в зародыш планеты (резонанс 1/1) либо выходят из резонанса.
Таким образом, радиальная миграция тел в модели идеального резонанса может быть представлена в виде последовательности переходов малых тел из одних резонансов в другие. При этом аккрециия малых тел на зародыш планеты будет происходить при значительно меньшей величине постоянной энергии h, чем при выходе из резонанса 1/1. Следовательно, рост массы зародышей планет начинается только после завершения радиальной миграции и завершения аккреции газа и малых тел на ядро протозвезды. В этом отношении кометы с большими эксцентриситетами орбит являются остаточными признаками завершения аккреции на ядро протозвезды, а астероиды резонанса 1/1 — признаками завершения аккреции малых тел на зародыш планеты.
p
p
2.2. Механизм формирования и роста зародышей планет
Если рост массы зародышей планет начинается только после завершения радиальной миграции малых тел и завершения аккреции газа и газовых уплотнений на ядро протозвезды, то возможны варианты эволюции протозвезды. Газодинамическое равновесие ядра и оболочки протозвезды может наступить только после ее вступления на главную последовательность. Тогда звезда может не иметь планетную систему. В этом случае экваториальный диск будет аккреционным и диск не эволюционирует в протопланетный. Если же равновесие наступит до вступления протозвезды на главную последовательность, то в условиях устойчивого равновесия будут формироваться зародыши планет и в целом планетная система. Состояние равновесия в газопылевой оболочке может наступить несколько раз. В этом случае будут формироваться различающиеся друг от друга планетные группы.
При выполнении условия устойчивого равновесия (27) в пылевой оболочке протозвезды будут формироваться плотностные волны. Вид этих волн представлен на рис. 3, а. Согласно рис. 3, а максимумы волн плотности имеют кольцеобразную форму. Вариация плотности (и давления) может быть причиной начала переходных процессов: конденсации, кристаллизации, коагуляции и др. Конденсирующие газопылевые частицы будут оседать на экваториальную плоскость вращения протозвезды, формируя протопланетные кольца.
Дальнейший процесс аккумуляции частиц пыли и газа будет происходить в экваториальных дисках. Процесс аккумуляции можно представить при помощи модели ограниченной задачи трех тел с двумя частицами малой массы (модель неупругих столкновений частиц) [3] или в более сложной модели метода крупных частиц [12]. При определенной плотности распределения частиц [4] в процессе формирования и роста тел будет действовать механизм орбитального резонанса (разд. 2.1). В случае большой скорости аккреции (аккреционные диски), орбиты частиц экваториального диска будут иметь большие эксцентриситеты. По таким орбитам будет обеспечиваться рост массы ядра протозвезды. При малой скорости аккреции (протопланетные диски) эксцентриситеты орбит частиц будут малыми. Тогда частицы будут захвачены в орбитальные резонансы и будут участвовать в формировании и росте малых тел и зародышей планет в экваториальной плоскости.
Периодические перемещения максимумов плотности (рис. 3, а) могут формировать протопланетные кольца, расположенные на расстоянии й (й = 2Д0), Д0 = 5 а. е. (рис. 3, а) друг от друга. При гравитационном сжатии протозвезды радиус ядра Д0 будет уменьшаться. Тогда будет сокращаться и расстояние й между протопланетными кольцами. Согласно модели (34) волновой фрагментации пылевой оболочки протозвезды сначала будут формироваться протопланетные кольца планет более удаленных от центра протозвезды, а затем те, которые расположены ближе к ее центру. Приток частиц для быстрого роста зародышей всех планет может обеспечиваться
волновым механизмом (34) и быстрым оседанием частиц пылевой оболочки на экваториальную плоскость.
Согласно модели Сафронова [54] формирование и рост масс малых тел происходит в результате неупругих взаимных столкновений малых тел. Наибольшие из них, вследствие выпадения на них значительно более мелких тел, образовали зародыши планет. Предполагается, что зародыши двигались по круговым орбитам, расположенным в наиболее плотной центральной плоскости роя. Орбиты падающих на зародыши тел — это орбиты пересечения с круговыми орбитами зародышей. Вероятность столкновения с зародышами определяется по формуле пт2/<^00 [54]. В отличие от модели Сафронова предположим, что падение тел на «зародыши» планет происходит не по случайным траекториям, а по регулярным орбитам резонанса 1/1 (рис. 3, б). В такой модели зародыши будут формироваться [4] на одной из вершин равностороннего треугольника БЬЕ, где Б — центр протозвезды, Ь — одна из двух треугольных точек Лагранжа Ь4 и Ь5, Е — зародыш. В соответствии с барометрической формулой [54] при отсутствии турбулентности твердые частицы за короткое время будут оседать на экваториальную плоскость, образуя в ней слой повышенной плотности. При критических значениях плотности диска будут возникать тела километровых размеров [54]. В статье [49] показано, что гравитационная неустойчивость протопланетного диска и появление «зародышей» возможна на расстоянии 30-100 а. е. от центра протопланетного диска. В газовых дисках с массой менее 0.1М^, согласно модели (34) фрагментации пылевых оболочек протозвезд, гравитационная неустойчивость наступает в результате кольцевой фрагментации диска и диссипации оболочки протозвезды.
ЗАКЛЮЧЕНИЕ
Стандартная модель происхождения и эволюции Солнечной системы сформировался к настоящему времени в результате усилий многих исследователей: Сафронова [54], Витязева и Печерниковой [17, 18], Горькавого и Фридман [25] и многих других. Однако до настоящего времени не выяснен механизм формирования тел метровых, километровых размеров в дисках с массой менее 0.1 М^ [11, 15].
Согласно модели (34) волновых возмущений плотности пылевых оболочек протозвезд, планетные системы начинают формироваться еще до вступления звезд на главную последовательность. Еще на ранней стадии формирования звезд происходит разделение пылевых оболочек протозвезд и их экваториальных дисков на кольцевые фрагменты. Выполнение условия устойчивости газодинамического равновесия в пылевой оболочке протозвезд является необходимым условием для аккумуляции газопылевых частиц и формирования тел в экваториальных дисках с массой менее 0.1М^. Достаточное условие для формирования тел метровых и километровых размеров в таких дисках — существование в протопла-нетных кольцах устойчивых точек гравитационного равновесия. Такими точками равновесия являются устойчивые точки либрации [4].
Механизм орбитального резонанса 1/1 дает объяснение причин появления зародышей планет и быстрого роста их масс. Эффективность действия такого механизма подтверждается результатами наблюдений астероидов троянской группы вблизи орбиты Юпитера.
Результаты наблюдений молодых звезд TW Hya, HD 163296, HD 16 9142, HD 97048 [23, 24,
52, 57] показывают, что газопылевые диски вокруг этих звезд имеют кольцеобразную структуру. Такая структура может быть связана с устойчивостью пылевых оболочек этих звезд. Газопылевые диски звезд MWC 758, HD 100453, HD 135344B [23, 24] имеют спиралевидную структуру. Такая структура могла быть сформирована в результате нарушения устойчивости пылевых оболочек этих звезд.
Автор настоящей работы выражает благодарность академику АН РТ Н. А. Сахибуллину за замечания, касающиеся миграции «зародышей» планет, профессору СПбГУ К. В. Холшевникову за замечания, связанные с обоснованием методов усреднения, а также всем сотрудникам отдела небесной механики ГАИШ МГУ и кафедры астрономии и космической геодезии КФУ за участие в обсуждениях работы.
СПИСОК ЛИТЕРАТУРЫ
1. Абдульмянов Т. Р., Загретдинов Р. В. // Кинем. и физ. неб. тел. 1994. 11, №4. С. 34.
2. Абдульмянов Т. Р. // Астрон. вестн. 2001. 35, №5. С. 449.
3. Абдульмянов Т. Р. // Вестник КГЭУ. 2011. №2 (9). С. 64.
4. Абдульмянов Т. Р. // Вестн. Моск. Ун-та. Физ. Астрон. 2014. №6. С. 115. (Abdul'myanov T.R. // Moscow Univ. Phys. Bull. 2014. 69, N 6. P. 558.)
5. АЬды-lmyanov T.R. // In: Rep. of the Intern. Symp. «Lunar Explor. and Space Techn. Heritage»/ Kazan: Kazan federal university, 2016. P. 41.
6. Abdulmyanov T.R. // Abstracts EPSC2017 - 6, 17-27 Sep. Riga. Latvia. 2017.
7. Abdulmyanov T.R. // Abstracts PFE - 17, 25-27 Sep. Jena, Germany. 2017. P. 30.
8. Abdulmyanov T.R. // Abstracts EPSC2018 - 12, 16-21 Sep. Berlin. Germany. 2018.
9. Абдульмянов Т. Р. // Современная звездная астрономия 2018. 22-26 октября. Москва, ГАИШ МГУ. 2018. Т. 1. С. 68.
10. Абдульмянов Т. Р. // В сб. матер. Междун. науч. семинара «GRACOS-18», Казань: Казан. федер. университет, 2018. С. 103.
11. Armitage P.J. Lecture notes on the formation and early evolution of planetary systems. In: Astrophysics of Planet Formation. Cambridge University Press. 2010.
12. Белоцерковский О.М., Давыдов Ю.М. Метод крупных частиц в газовой динамике. М., 1982.
13. Боголюбов Н.Н., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. М., 1958.
14. Бронштейн И.Н., Семендяев К. А. Справочник по математике. М., 1986.
15. Брюно А. Д. Ограниченная задача трех тел. М., 1990.
16. Wilner D. J., Holman M. J., Kuchner M. J., Ho P. T. P. // Astrophys. J. 2002. 569. L115.
17. Витязев А. В., Печерникова Г. В., Сафронов В. С. Планеты земной группы: происхождение и ранняя эволюция. М., 1990.
18. Витязев А. В., Печерникова Г. В. // Активность звезд и Солнца. 17-18 декабря 2010. Москва, ГАИШ МГУ. 2010. 1. С. 161.
19. Wolf S., Gueth F., Henning T. et al. // Astrophys. J. 2002. 566, L97.
20. Wolf S, Hillenbrand L.A. // Astrophys. J. 2003. 596. P. 603.
21. Wolf S. // Astrophys. Space Sci. 2008. 313. P. 109.
22. Vorobyov E., Basu S. // Astrophys. J. 2006. 714. P. 133.
23. Garufi, Meeus, Benisty et al. // Astron. and Astroph. 2017. 603.
24. Garufi A. // Abstracts PFE-2017, Jena 25-27 Sep., Germany, 2017, 12
25. Горькавый Н.Н., Фридман А.М. Физика планетных колец. М., 1994.
26. Goto M., Usuda T., Dullemond C.P. et al. // Astrophys. J. 2006. 652. P. 758.
27. Гребеников Е.А. Методы усреднения в прикладных задачах. М.: Наука, 1986.
28. Garfinkel B., Jupp A., Williams I. // Astron. J. 1971. 76, № 2.
29. Garfinkel B. // Celest. mech. 1976. 13, N 2.
30. Garfinkel B. // Astron. J. 1977. 82, N 5. P. 368.
31. Garfinkel B. // Celest. mech. 1980. 22. P. 267.
32. Goto M., Usuda T., Dullemond C.P. et al. // Astrophys. J. 2006. 652. P. 758
33. Gräfe C., Wolf S., Guilloteau S. et al. // Astron. and Astroph. 2013. 553. A69.
34. Greaves J. S., Mannings V., Holland W. S. // Icarus. 2000. 143. P. 155.
35. Dent W.R.F., Walker H.J., Holland W.S., Greaves J. S. // Mon. Not. R. Astron. Soc. 2000. 314. P. 702.
36. Джакалья Г. Е. О. Методы теории возмущений для нелинейных систем. М., 1979.
37. Johansen A., Oishi J., Mac Low M.-M. et al. // Nature. 2007. 448. P. 1022.
38. Calvet N., D'Alessio P., Hartmann L. et al. // Astrophys. J. 2002. 568. 1008.
39. Klahr H., Rozyczka M., Dziourkevitch N. et al. // Planet Formation: Theory, Observations and Experiments. Cambridge University Press. 2006. P. 42
40. Kley W. // Mon. Not. R. Astron. Soc. 1999. 303. P. 696.
41. Meru F. // Abstracts PFE-2017, 25-27 Sep., Jena, Germany. P. 19.
42. Morbidelli A., Raymond S.N. Challenges in Planet Formation. ArXiv: 1610.07202v2 [astro-ph.EP] 25 Okt 2016.
43. Мюррей К., Дермотт С. Динамика Солнечной системы. М.: ФИЗМАТЛИТ, 2009.
44. Nayakshin S. // Mon. Not. R. Astron. Soc. 2010. 408, L36.
45. Nayakshin S. // Mon. Not. R. Astron. Soc. 2011. 413, 1462.
46. Nayakshin S. // Mon. Not. R. Astron. Soc. 2014. 441, 1380.
47. Nesvorny D. Dynamical Evolution of the Early Solar System. ArXiv: 1807.06647v1 [astro-ph.EP] 17 Jul 2018.
48. Пуанкаре А. Избранные труды. Т. 1. Новые методы небесной механики. М., 1971.
49. Rafikov R.R. // Astrophys. J., 621, 2005, L69.
50. Rice W.K.M., Wood K, Armitage P.J. et al. // Mon. Not. R. Astron. Soc. 2003. 342. P. 79.
51. РольфсК. Лекции по теории волн плотности. М., 1980.
52. Ruge J.P., Wolf S., Uribe A.L., Klahr H.H. // Astron. and Astroph. 2013. 549. A97
53. Sauter J., Wolf S. // Astron. and Astroph. 2011. 527. A27.
54. Сафронов В. С. Эволюция допланетного облака и образование Земли и планет. М., 1969.
55. Снытников В.Н., Стояновская О. П. // Земля на ранних этапах развития солнечной планетной системы. 28-30 ноября 2016. М., ГАИШ МГУ. 2017. С. 43.
56. Сурдин В. Г., Ламзин С. А. Протозвезды: где, как и из чего формируются звезды. М., 1992.
57. Facchini S. 2017, Abstracts PFE-2017, 25-27 Sep., Jena, Germany, 11.
58. Holland W.S., Greaves J.S., Dent W.R.F. et al. // Astrophys. J. 2003. 582. P. 1143.
59. Shakura N.I., Sunyaev R.A. // Ap. J. 1973. 24. 337.
60. Schegerer A.A., Wolf S., Ratzka Th., Leinert Ch. // Astron. and Astroph. 2008. P. 779.
61. Schegerer A.A., Wolf S. // Astron. and Astroph. 2010. 517. A87.
62. Яворский Б.М., Детлаф А. А. Справочник по физике. М.: Наука, 1965.
A General Model of Celestial Body Formation from Initial Condensation of Gas and Dust Particles to the «Embryos» of Planets
T. R. Abdulmyanov
Department of Engineering Cybernetics, Kazan State Power Engineering University, Krasnoselskaya st. 51.
Kazan 420066, Russia.
E-mail: abdulmyanov.tagir@yandex.ru.
This paper analyzes the main ideas of the formation of celestial bodies in accretion and protoplanetary disks. Mechanisms of the radial and orbital migration of celestial bodies are considered. According to the results of modeling, the formation of small bodies of the Solar System and the growth in their size and mass occurs at stable libration points of orbital resonances. The appearance of a single embryo in protoplanetary rings is explained by of the action of the 1/1 orbital resonance mechanism. Using a model of wave perturbations of the dust shells of protostars and a model of orbital resonance, a general model of the formation and growth of "embryos" of planets for disks with masses less than 0.1MS was constructed. Keywords: protoplanetary disk, protoplanetary rings, orbital resonance, libration centers. PACS: 96.12.Bc Received 11 March 2019.
English version: Moscow University Physics Bulletin. 2019. 74, No. 4. Pp. 309-322.
Сведения об авторе
Абдульмянов Тагир Раисович — канд. физ.-мат. наук, профессор; тел.: (843) 562-43-36, e-mail: abdulmyanov.tagir@yandex.