2024 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 11 (69). Вып. 1
МАТЕМАТИКА. МЕХАНИКА. АСТРОНОМИЯ
АСТРОНОМИЯ
УДК 521.135 МБС 70F07
Орбитальная прецессия в ограниченной задаче трех тел: точные представления*
А. А. Березина
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9 Институт прикладной астрономии РАН,
Российская Федерация, 197110, Санкт-Петербург, ул. Ждановская, 8
Для цитирования: Березина А. А. Орбитальная прецессия в ограниченной задаче трех тел: точные представления // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 1. С. 185-197. https://doi.org/10.21638/spbu01.2024.113
Рассмотрены аналитические представления для скорости апсидальной прецессии в эллиптической плоской ограниченной задаче трех тел в случае, когда орбита возмущающего тела — внешняя по отношению к орбите тестовой частицы. Аналитические выражения сопоставлены с численными данными, полученными в ходе массовых расчетов скорости прецессии. Получено аналитическое выражение для скорости апсидальной прецессии в виде степенного ряда по параметру, равному отношению больших полуосей орбит тестовой частицы и возмущающей планеты. Показано, что аналитические выражения для скорости апсидальной прецессии частицы надежны лишь на удалении от зоны нестабильности в окрестности орбиты возмущающей планеты. Вблизи же люка Виздома линейная вековая теория перестает работать. Предложена поправочная эмпирическая формула для вычисления скорости апсидальной прецессии при относительно высоких (но менее 0.5) эксцентриситетах частицы и возмущающей планеты. Предложенные формулы применены к описанию прецессии орбит в реальных экзо-планетных системах.
Ключевые слова: ограниченная задача трех тел, прецессия линии апсид, люк Виздома, экзопланеты, планетные системы.
* Работа выполнена при финансовой поддержке Российского научного фонда (грант №22-2200046).
© Санкт-Петербургский государственный университет, 2024
1. Введение. В данной работе рассмотрены аналитические представления скорости изменения долготы перицентра w тестовой пассивно гравитирующей частицы на внутренней орбите в эллиптической плоской ограниченной задаче трех тел.
Хеппенхаймер [1], Уитмайр и др. [2], Тибо и др. [3] исследовали движение тестовой частицы в двойной звездной системе в плоском случае, когда частица обращается вокруг одного из компонентов двойной системы. Предполагалось, что эксцентриситет частицы ei ^ 1, отношение больших полуосей орбит частицы и двойной а = ai/a,2 ^ 1. Также учитывалось, что эксцентриситеты двойных, как правило, не малы, а лежат в широком диапазоне: от 0 до 1 со средним значением 0.55 [1, 4, 5]. Пусть массы измеряются в массах Солнца, расстояния — в астрономических единицах, время — в годах. Таким образом, гравитационная постоянная G = 4п2. Далее будем использовать эти единицы измерения.
Если взять в начальный момент времени ei = 0, тогда из уравнений движения можно получить [1-3]
sin ut 3п Ш2 OJ"1 , 2\ -3/2
tanUJi = — --, U=—--777-—о— (1 — во ) , (1)
1 — cos ut 2 тоу2 а\ v 2> к '
где mo — масса главного компонента двойной системы; m2 — масса меньшего компонента двойной системы; wi — долгота перицентра тестовой частицы. Вводя переменную y = ut/2, можно записать (см. формулу (13.6) в [6], а также [7]):
w1 = y + const, (2)
где const — константа, зависящая от того, к какой четверти принадлежит угол y. Таким образом, в некотором смысле величину u можно рассматривать как скорость прецессии линии апсид частицы.
Мюррей и Дермотт [8] рассмотрели случай плоской круговой ограниченной задачи трех тел, в которой внешнее тело возмущает внутреннее тело пренебрежимой массы. В отсутствие резонансов формула для скорости изменения долготы перицентра имеет вид
1 db(1) db(0) 1 m2 2w w o l(0) , obW , 3/2 2 3/2 , ,
w\ =--nia Wi, W\ = —Aab\L + ¿b\ L + a----a -—, (3)
4 m0 3/2 3/2 ¿a da
где ni — среднее движение тестовой частицы; а = ai/a2, b(+)i/2 — коэффициент Лапласа, которые определяются как
1ЬЩа) = ±Г_COS^__(4)
2 8 к > 2тг У0 (l-2acosф + a2y, [ '
где s — полуцелое число [8].
о /2
Заметим, что при подстановке п\ = 27г-у/тоо/а/ коэффициент перед скобками в формуле (1) отличается от соответствующего в формуле (3) на множитель 3а. В дальнейшем можно видеть (см., например, раздел 4), что формула (1) без учета эксцентриситета e2 является первым членом разложения по а формулы (3).
Применяя метод усреднения Гаусса, Мюррей и Дермотт [8] также получили представление апсидальной скорости прецесии в эллиптической плоской ограниченной задаче трех тел. Эксцентриситет частицы предполагался малым. Усредненное
влияние внешнего тела на внутреннее определялось путем распределения массы внешнего тела вдоль его орбиты. После вычисления текущей скорости прецессии она усреднялась по истинной аномалии внутреннего тела. Результат данного подхода тождественен расчету усредненной скорости прецессии через уравнения Лагран-жа с возмущающей функцией, в которой отбрасываются все члены, не зависящие от средних долгот. Усредненная скорость прецессии в этом случае (с точностью до членов второго порядка в разложении возмущающей функции по эксцентриситетам)
{т i)
1 m-2 2
—-п\о.
4 m°
W2 = -3-
db
da
^2 a
H--— W2 cos(m2 - mi)
ei 2
(0) 3/2
2b
(i) 3/2
db
2a-
(i) 3/2
db!
da
+
(2) 3/2
(5)
da
где m2 — долгота перицентра возмущающего тела.
Расчеты скоростей прецессии важны для выявления вековых резонансов в планетных системах. Наличие и взаимодействие вековых резонансов влияет на долговременную устойчивость системы и даже способно приводить к дезинтеграции системы (см., например, [9]).
Нашей целью является вывод компактных и удобных для приложений представлений для скорости прецессии в зависимости от параметров задачи.
2. Численное моделирование. Мы провели массовые расчеты с целью численного определения скоростей апсидальной прецессии. Вычисления проводились в Межведомственном суперкомпьютерном центре РАН. Использовался несимплекти-ческий интегратор с адаптивным временным шагом 15-го порядка [10] в программном пакете REBOUND [11]. Центральное тело (звезда) имело единичную массу, возмущающее тело (планета) находилось на единичном расстоянии от звезды, долгота перицентра планеты задавалась нулевой. Время моделирования составляло 1000 оборотов планеты вокруг звезды. Варьировались масса возмущающей планеты m-2 от Ю-7 до 10~3, большая полуось орбиты тестовой частицы ai от 0.01 • ^ff Д° 0.95 • Y^ff > эксцентриситет орбиты планеты е2 от 0.01 до 0.1, а также начальный эксцентриситет частицы ei от 0.01 до 0.1 и начальная долгота перицентра орбиты частицы mi от 0 до 2^.
Заметим, что входящие в формулы (3) и (5) эксцентриситет частицы ei и долгота перицентра mi меняются с течением времени. Рассмотрим, как это влияет на скорость прецессии mi. На рис. 1 показаны численные и теоретические данные изменения ei, mi и mi со временем t, а также зависимость эксцентриситета от долготы перицентра на типичном примере при длительном моделировании (50 тыс. оборотов возмущающей планеты вокруг звезды). Время измеряется в периодах обращения внешней планеты T2. Теоретические кривые соответствуют уравнениям движения в линейной вековой теории [8]:
т i
1 2
—m2nia
db
(°) d2 b(°) 1/2 a__V£ __
da da2 aei
e2
2b/ - 2a
db
(i) i/2
e\ = — m2niae2
db(i)
26V2 " da
da
da2
d2b(i) x2Q °l/2
da2
2^ • - SI
sin mi.
cos mi
(6) (7)
2
1
Рис. 1. На верхних графиках показаны численные и теоретические данные зависимости эксцентриситета ei и долготы перицентра wi от времени t для типичной системы выборки. T2 — период обращения внешней планеты. На нижнем левом графике — зависимость эксцентриситет — долгота перицентра; на нижнем правом — зависимость скорости прецессии орбиты тестовой частицы w 1 от времени. Кривыми показаны теоретические данные согласно линейной вековой модели, точками — модельные данные, пунктиром — скорость прецессии согласно формуле (5), штрихпунк-тиром — скорость прецессии согласно формуле (3).
£ 0.000500.00025 -
в nnnnnnJ . • ---,_
0.0 0.1 0.2 0.3 0.4 0.5 ai/a2
Рис. 2. Численные данные и аналитические значения скорости апсидальной прецессии тестовой частицы Wi в зависимости от отношения больших полуосей орбит частицы и возмущающей планеты ai/a2 при эксцентриситете орбиты тестовой частицы ei = 0.02 и начальной долготе перицентра орбиты частицы wi = п. На левом графике эксцентриситет возмущающей планеты e2 = 0.04, масса планеты m2 = 10-4 (в единицах массы центральной звезды); на правом: e2 = 0.06, m2 = 10-5. Точки обозначают данные численных экспериментов, синяя кривая соответствует формуле (5), зеленая — формуле (3), оранжевая — формуле (1).
ai/a2
Видим хорошее соответствие теории и результатов численного моделирования, но с течением времени их расхождения возрастают. Что касается скорости апсидальной прецессии, формула (5) соответствует данным в начальный момент времени, а формула (3) отражает среднее значение на всем промежутке времени. Действительно, усреднение сов(ет2 — по времени в формуле (5) дает ноль, и формула (5) приводится к виду (3).
На рис. 2 представлены типичные графики скорости прецессии, согласно результатам численного моделирования и аналитическим представлениям (2), (4), (6), в зависимости от отношения а = ах/а2 при различных параметрах системы. Видно, что наилучшее соответствие имеет аппроксимация (5).
3. Люки Виздома. Согласно Виздому [12] и Дункану и др. [13], в ограниченной круговой задаче трех тел при в\ ~ 0 вблизи орбиты планеты существует хаотическая область вследствие перекрытия резонансов средних движений вида где р> 1. Радиальный размер люка Виздома (см. формулу (11.53) в книге [6]) определяется по формуле
есг ~ 1.62М2/7, (8)
где (1 = "|2 .
" п Ш0+Ш2
Рис. 3. Зависимость скорости апсидальной прецессии тестовой частицы Шi от отношения ai/a2. Если иное не указано на графике, параметры системы составляют: ei = 0.02, roi = п, в2 = 0.04, m.2 = 10_Точки — данные численного моделирования, синяя кривая соответствует формуле (5), зеленая вертикальная линия отмечает границу люка Виздома согласно формуле (9), красная линия — согласно формуле (8).
Полуширину области нестабильности можно также получить используя двумерное столкновительное отображение [6]; в итоге она определяется формулой
есг ~ 1.67(Ме1)1/5, (9)
где в1 — начальный эксцентриситет тестовой частицы (см. формулу (11.50) в [6]).
На рис. 3 представлены типичные графики зависимости скорости прецессии от а при различных параметрах системы, указано положение люков Виздома согласно формулам (8) и (9). Мы нашли, что во всех рассматриваемых нами случаях численного моделирования формула (9) задает зону нестабильности шире, чем формула (8). Рис.3 иллюстрирует важный факт: отклонения теории от результатов численного моделирования резко возрастают вблизи и внутри хаотической области (люка Виздома).
Таким образом, рассмотренные нами аналитические выражения для скорости апсидальной прецессии частицы надежны лишь на удалении от зоны нестабильности в окрестности орбиты возмущающей планеты.
4. Новые представления. Формула (5) хорошо аппроксимирует результаты численного моделирования при параметрах системы, указанных в разделе 2, за исключением больших значений а, но не всегда удобна для приложений. Представим ее в более удобном и компактном виде как ряд по степеням а. Заметим, что формулы (3) и (5) получены из разложений по эксцентриситету, который предполагается малым; коэффициенты этих разложений представляют собой точные (через коэффициенты Лапласа) зависимости от а.
Коэффициенты Лапласа представим в виде ряда:
160, = s(s+l)...(s+j-l)^
о 3 у ' I •::..../
1++
S(S + 1)(S + J)(S+J + 1) 4 1-2(J + 1)(J+2)
(10)
(см., например, [8]). При этом их производные по а:
ОЪ^ = а (ъЛ — 2о$Ъ + Ъ^) . (11)
Рассмотрим сумму, обозначенную нами ранее через и подставим в нее значения производных коэффициентов Лапласа:
аъ(Ц аъ(0)
= -3 аЪУ3/>2 + 2Ь\2 + а--а'
= -3аЬ<% + 2ЪЦ + За ^ + а2) - 6а2б/ + ^аб/. (12)
Из разложений коэффициентов Лапласа в ряды можно видеть, что в ряде для у нас будут только нечетные степени о; при этом при первой степени будет коэффициент 3, при третьей — Обозначим через № (а) отношение вида
л{1)/ , _ В{8 + 1) ■ ■ ■ {8 + | ~ 1) • {3+3){8+3 + 1) • • • {8 + 3 + | ~ 1) 1-2...|.(, + !)(,+2)...(, + §) •
Тогда в разложении (12) для при г-й степени о коэффициент будет равен
-б4°/2(г -1)+б4у2(г -1)+з^^
8 (14)
Нетрудно видеть, что
'2А<%Ла) = (з+ ¡У =
42)(а) = 2 + АР(а), А^(г - 1) = + ^ ~ ^^(г - 3). 8 К ' (в + 1) (2 + |) > ((г — 1)/2) 8 ^ '
С учетом этих соотношений получаем коэффициент при г-й степени о, равный
-М а<°> и - 1) - (г + V (г!!)2 (16)
Теперь рассмотрим ^2, и аналогично исключим производные:
цг2 = -^111 + 260/2 + 2а^ + =
ао 3/2 ао ао
1,(1) ^ 1_ о Л> , „Л . (1) , 3,(3)
26/ + 12аб/ - 3 (Ё + 2а2) б/ + | (17)
Здесь также получаются только нечетные степени о. При первой степени коэффи-
15 »105
при третьей —
15 105
циент будет равен --у-, ПРИ третьей — при г-и степени:
б4°/2(г - 1) + 244°/2(г - 1) - - 1) - 30АЦ(г - 3) + - 3). (18)
При этом
8 1 ; (2+§)(3+§)(5 + 1)(5 + 2) 8 1 ; 1 ;
Тогда коэффициент при г-й степени а приводится к виду
( 2г + 5 \ „(0) , ( 2г + 5 \ (¿!!)2 , ,
"4 I1 + ) ^ " :1) = Л1 + ^ПШ) 2*-1((4)!)3- (20)
Тогда разложение будет давать в общей сумме нечетные степени а, а разложение W^2 — четные, на единицу большие (с учетом а перед Wг2 в формуле (5)).
Выпишем в явном виде разложение для скорости прецессии до 21-го порядка по а:
, > 1 т2 о
4 то
45 2 525 4 11025 6 218295 8 2081079 10
3 H--а Л--а4 H--а6 H--а8 H--а10+
8 64 1024 16384 131072
19324305 12 703956825 14 25264228275 16 224009490705 18
H--а H--а14 H--а1Ь H--а18+
1048576 33554432 1073741824 8589934592
1967210618373 20 е2 /15 105 3 4725 5 24255 7
Н--а"'--—аН--а Ч--а Ч--а +
68719476736 2ei \ 2 8 256 1024
945945 9 4459455 11 328513185 13 1486131075 15
H--а H--а11 H--а13 H--а15 +
32768 131072 8388608 33554432
106109758755 17 468383480565 19 16393421819775 21\
H--а H--а H--а cos( ujn - m-\ )
2147483648 8589934592 274877906944 ) v " ;
+ 7—m3 [3 + 5.6250а2 + 8.2031а4 + 10.7666а6 + 13.3237а8 + 15.8774а10+ 4 т0
+ 18.4291а12 + 20.8795а14 + 23.5291а16 + 26.0781а18 + 28.6267а20 - (7.5000а+
2ei
+13.1250а3 + 18.4570а5 + 23.6865а7 + 28.8680а9 + 34.0229а11 + 39.1618а13+
+44.2902а15 + 49.4112а17 + 54.5270а19 + 59.6389а21) cos(i2 - ш1)] . (21)
Выясним, насколько точно полученное разложение представляет численные данные. Построим графики относительной разности скорости прецессии согласно формуле (5) (обозначим ее как inexact) и согласно разложению (21) до 21-й степени а с рациональными коэффициентами (обозначим как iexpansion). На рис. 4 показаны
ro,xact-ro,xpanBion ра3дИЧНЫХ Па_
типичные графики относительной разности 6 = ш
раметрах системы. Мы видим, что при больших а точность разложения существенно
ai/эз aila2
Рис. 4. Зависимость относительной разности скоростей апсидальной прецессии 6 согласно формуле (5) и разложению (21) от отношения больших полуосей орбит тестовой частицы и возмущающей планеты «1/0,2 при эксцентриситете орбиты тестовой частицы вх = 0.02 и начальной долготе перицентра орбиты частицы шх = п. На левом графике эксцентриситет возмущающей планеты в2 = 0.06, масса планеты т2 = 10_5; на правом: в2 = 0.04, т2 = 10_4.
падает. При а < 0.75 относительная разность находится в пределах примерно 0.01, при а < 0.5 — в пределах 10~5.
5. Случай сильно вытянутых орбит. Формула (5) для скорости апсидальной прецессии выведена в предположении малых эксцентриситетов орбит возмущающего тела в2 и тестовой частицы в\. Рассмотрим случай умеренных эксцентриситетов: в! < 0.5, в2 < 0.5. При в1,в2 ~ 0.5 формула (5) гораздо хуже описывает значение скоростей прецессии даже при небольших а. При разложении возмущающей функции до четвертого порядка по эксцентриситетам (см. Приложение Б в книге [8]) можно уточнить разложение (5), но при этом оно станет громоздким. На левых графиках рис. 5 показаны скорости апсидальной прецессии в зависимости от начального эксцентриситета орбиты частицы в! при разных эксцентриситетах орбиты возмущающего тела в2. Видим, что формула (5) хорошо описывает данные лишь при в2 < 0.1, а разложение до четвертого порядка хорошо работает вплоть до в2 = 0.3. Но при больших в2 оба представления существенно отклоняются от численных данных.
Дальнейшее разложение возмущающей функции будет чрезвычайно громоздким. Рассмотрим зависимость скорости прецессии от эксцентриситета орбиты планеты в2 при разных начальных эксцентриситетах частицы в1, приведенную на правых графиках рис. 5. В данном случае с изменением в! и в2 характер функции меняется. Рассмотрим разность модельных данных формулы (5) и обозначим ее как Лги. На рис. 6 показаны типичные графики зависимости Лги от параметров планетной
Рис. 5. Зависимость скорости апсидальной прецессии тестовой частицы Ш 1 от начального эксцентриситета частицы вх при различных значениях эксцентриситета возмущающей планеты в2 (слева) и от эксцентриситета планеты в2 при различных начальных эксцентриситетах частицы вх (справа). Масса возмущающей планеты т2 = 10_5, начальная долгота перицентра орбиты частицы П71 = 7Г и а = 0.532Точки — данные численного моделирования, синяя кривая соответствует формуле (5), оранжевая — разложению возмущающей функции до четвертого порядка по эксцентриситетам частицы и возмущающей планеты (см. Приложение Б в книге [8]), зеленая — формуле (5) с поправкой (22).
0.1 0.2 0.3 0.4 0.5 0.6 0.7
ег
Рис. 6. Зависимость погрешности формулы (5) от параметров системы. Точками обозначены результаты расчетов, кривыми представлена аппроксимация (22). При варьировании одного из параметров, остальные были равны а = 0.5321 Тв2 ■ т1 = Ю-Б, П71 = 7Г, е\ = 0.2, еп = 0.3.
системы. Разности обозначены точками. Видим, что зависимость от о в логарифмическом масштабе имеет степенной вид, от в2 — линейный, от ех — гиперболический, от т2 — линейный, от шх — вид косинуса. Кривыми на рисунке обозначена эмпирическая формула следующего вида:
= (в— + С—со8(то-т1)\ (22)
V 1 — е2) \ е1 ех )
где А = 3.67 ± 0.03; В = 1.34 ± 0.02; С = —2.38 ± 0.04 — коэффициенты, полученные с использованием нелинейного метода наименьших квадратов по данным в области: 0.05 < ех < 0.5, 0.05 < е2 < 0.5, 10-7 < т2 < 10-4, 0 < < 2п, о < 0.3. Под ех и шх подразумеваются начальные эксцентриситет и долгота перицентра орбиты частицы. Значение о ограничено до 0.3 с целью исключить выбросы, не ложащиеся на общую закономерность, которые появляются вблизи люка Виздома. В качестве ошибок определения коэффициентов А, В и С указаны значения одного стандартного отклонения.
На рис.5 зеленая кривая соответствует формуле (5) с поправкой (22). Такой подход с гораздо лучшей точностью описывает значение скоростей прецессии при умеренных ех, е2.
6. Применение к реальным системам. Отберем реальные экзопланетные системы, к которым можно применить рассмотренные формулы. Воспользуемся базой данных ехор1апе1.еи. Выберем системы с двумя планетами, у которых известны массы и большие полуоси, с отношением масс < 0.1, где менее массивная планета расположена на внутренней орбите. Среди этих систем отбирались такие, у которых небольшая или неизвестная разность наклонений орбит. Таких систем было найдено девять (см. таблицу). В таблице указаны массы звезд, известные параметры орбит планет и вычисленная скорость прецессии для внутренней планеты. В системах звезд 55 Спс и Кер1ег-167 есть и другие планеты, но с неизвестными
СО
Список планетных систем, к которым применимы формулы для скорости прецессии в рамках плоской ограниченной задачи трех тел
Система Масса звезды, Мзип Планета Масса планеты, М]ир а, а. е. е г, град. ■го, град. Т ■ -'-реп? юлианские дни ■го, рад/Тг
Кер1ег-167 0.77 Кер1ег-167 с1 0.006 0.1405 0.12 89.352 — — 4.77 X 10~4
Кер1ег-167 е 4 1.89 0.062 89.976 201 - —
Кер1ег-93 1.09 Кер1ег-93 Ь 0.0143 0.053 0 89.18 — — 4.47 X 10~Б
Кер1ег-93 с 8.5 4.5 - - - - -
55 Спс 1.015 55 Спс е 0.02703 0.015439 0.028 90.36 170 2450000 3.81 X Ю-4
55 Спс Ь 0.84 0.11339 0.0023 89.73 110 2453035 -
НБ 38529 1.48 НБ 38529 Ь 0.93 0.131 0.248 57 95.9 50020.19 9.00 X Ю-4
НБ 38529 с 23.7 3.695 0.36 48.3 22.52 48002 —
К2-290 1.19 К2-290 Ь 0.0343 0.0923 0 88.14 90 — 5.82 X Ю-4
К2-290 с 0.774 0.305 0 89.37 90 — —
КО 1-12 1.5 К01-12 Ь 1.1 0.154 0.34 88.95 — — 4.63 X 10~4
К01-12 с 22 4.2 0.65 — -89 2455911 —
Кер1ег-117 1.13 Кер1ег-117 Ь 0.094 0.1445 0.0493 88.74 254.3 — 5.74 X 10~3
Кер1ег-117 с 1.84 0.2804 0.0323 89.64 305 - —
Т01-1670 1.21 Т01-1670 Ь 0.0434 0.103 0.59 86.87 163.6 — 1.63 х 10~3
Т01-1670 с 0.63 0.249 0.09 88.84 105.5 - —
К01-984 0.928 К01-984 Ь 0.0544 0.0504 0.12 — 285 — —5.52 X КГ4
К01-984 с 0.6601 0.1467 0.38 45 84 - -
о> Примечание. Массы звезд указаны в массах Солнца, массы планет — в массах Юпитера. Т2 обозначает период обращения внешней
V планеты; г — наклон орбиты планеты (угол между плоскостью орбиты планеты и картинной плоскостью).
Со
с*
массами. Если известны эксцентриситеты и долготы перицентров и если эксцентриситет внутренней планеты не нулевой, скорость прецессии рассчитывается по формуле (5) с поправкой (22); в противном случае — по формуле (3). Для внутренних планет HD 38529 b, K2-290 b, Kepler-117 b, T0I-1670 b и KOI-984 b известна долгота перицентра, при этом эпоха прохождения перицентра известна для 55 Cnc e и HD 38529 b. Эти данные были получены посредством использования доплеровской спектрометрии (см., например, [14]). Зная скорости прецессии, мы можем прогнозировать изменение долготы перицентра и эпохи его прохождения. При несоответствии прогноза и наблюдения можно предполагать существование других планет в системе, возмущающих движение маломассивной, и оценивать массы и параметры орбит этих планет.
7. Заключение. Знание скоростей апсидальной прецессии и умение их предсказывать важны для исследований динамики планетных систем. Нами рассмотрены и проанализированы известные аналитические представления для скорости апси-дальной прецессии в эллиптической плоской ограниченной задаче трех тел. При этом аналитические представления сопоставлены с данными, полученными в ходе массовых расчетов скоростей апсидальной прецессии. Даже лучшая из представленных формул не позволяет предсказать скорость прецессии вблизи и внутри области нестабильности (люка Виздома).
Представлено аналитическое выражение для скорости апсидальной прецессии в виде ряда по a. Полученное выражение при разложении до 21-го порядка показывает хорошую точность, за исключением областей a ~ 1. Количественно определены пределы по a для точного представления скорости прецессии с помощью найденного разложения.
Предложена эмпирическая поправка к формуле (5) для определения скорости апсидальной прецессии при умеренных эксцентриситетах частицы и возмущающей планеты (менее 0.5).
Рассчитаны скорости апсидальной прецессии для ряда реальных экзопланетных систем.
Автор благодарит И. И. Шевченко за ценные советы и рекомендации при написании данной статьи, а также искренне признателен рецензентам за подробные полезные замечания.
Литература
1. Heppenheimer T.A. On the formation of planets in binary star systems. Astronomy and Astrophysics 65, 421-426 (1978).
2. Whitmire D., Matese J., Criswell L., Mikkola S. Habitable planet formation in binary star systems. Icarus 132, 196-203 (1998).
3. Thebault P., Marzari F., Scholl H. Relative velocities among accreting planetesimals in binary systems: The circumprimary case. Icarus 183, 193-206 (2006). https://doi.org/10.1016/j.icarus.2006.01.022
4. Mayor M., Duquennoy A., Halbwachs J.-L., Mermilliod J.-C. CORAVEL Surveys to Study Binaries of Different Masses and Ages. Complementary Approaches to Double and Multiple Star Research, ASP Conference Series 32 IAU Colloquium 135, 73-81 (1992).
5. Larson R. B. Implications of Binary Properties for Theories of Star Formation. The Formation of Binary Stars, Proceedings of IAU Symp. 200, 93-106 (2001). https://doi.org/10.48550/arXiv.astro-ph/0006288
6. Shevchenko I.I. Dynamical Chaos in Planetary Systems. Springer Nature (2020). https://doi.org/10.1007/978-3-030-52144-8
7. Demidova T. V., Shevchenko I.I. Spiral patterns in planetesimal circumbinary disks. Astrophys-ical Journal 805 (38), 8 (2015). https://doi.org/10.1088/0004-637X/805/V38
8. Мюррей К., Дермотт С. Динамика солнечной системы,, пер. с англ. под ред. И. И. Шевченко. Москва, Физматлит (2009).
9. Lithwick Y., Wu Y. Theory of Secular Chaos and Mercury's Orbit. Astrophysical Journal 739 (1), 17 (2011). https://doi.Org/10.1088/0004-637X/739/1/31
10. Rein H., Spiegel D. S. IAS15: a fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits. Monthly Notices of the Royal Astronomical Society 446 (2), 1424-1437 (2015). https://doi.org/10.1093/mnras/stu2164
11. Rein H., Liu S.-F. REBOUND: an open-source multi-purpose N-body code for collisional dynamics. Astronomy & Astrophysics 537 (A128), 10 (2012). https://doi.org/10.1051/0004-6361/201118085
12. Wisdom J. The resonance overlap criterion and the onset of stochastic behavior in the restricted three-body problem. Astronomical Journal 85, 1122-1133 (1980). https://doi.org/10.1086/112778
13. Duncan M., Quinn T., Tremaine S. The long-term evolution of orbits in the solar system. A mapping approach. Icarus 82 (2), 402-418 (1989). https://doi.org/10.1016/0019-1035(89)90047-X
14. Маров М.Я., Шевченко И. И. Экзопланеты: природа и модели. Успехи физических наук 190, 897-932 (2020). https://doi.org/10.3367/UFNr.2019.10.038673
Статья поступила в редакцию 20 мая 2023 г.;
доработана 30 августа 2023 г.; рекомендована к печати 31 августа 2023 г.
Контактная информация:
Березина Алена Александровна — лаборант-исследователь; st055707@student.spbu.ru
Orbital precession in the restricted three-body problem: Exact representations*
A. A. Berezina
St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation Institute of Applied Astronomy of the Russian Academy of Sciences, 8, ul. Zhdanovskaya, St. Petersburg, 197110, Russian Federation
For citation: Berezina A. A. Orbital precession in the restricted three-body problem: Exact representations. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11(69), issue 1, pp. 185-197. https://doi.org/10.21638/spbu01.2024.113 (In Russian)
Analytical representations for the rate of apsidal precession in the planar elliptical restricted three-body problem are considered in the case when the orbit of the disturbing body is external to the orbit of the test particle. Analytical expressions are compared with numerical data obtained during numerous calculations of the precession rate. An analytical expression is obtained for the apsidal precession velocity in the form of a power series with a parameter equal to the ratio of the semi-major axes of the orbits of the test particle and the disturbing planet. It is shown that the analytical expressions for the velocity of the apsidal precession of the particle are reliable only at distances not close to the instability zone near the orbit of the disturbing planet. Near the Wisdom gap, the linear secular theory stops working. A correction empirical formula is proposed for calculating the apsidal precession rate at relatively high (less than 0.5) eccentricities of the particle and the disturbing planet. The proposed formulas are applied to the description of the precession of orbits in real exoplan-etary systems.
Keywords: restricted three-body problem, apsidal precession, Wizdom gap, exoplanets, planetary systems.
*The work is supported by Russian Science Foundation (grant no. 22-22-00046).
References
1. Heppenheimer T. A. On the formation of planets in binary star systems. Astronomy and Astrophysics 65, 421-426 (1978).
2. Whitmire D., Matese J., Criswell L., Mikkola S. Habitable planet formation in binary star systems. Icarus 132, 196-203 (1998).
3. Thebault P., Marzari F., Scholl H. Relative velocities among accreting planetesimals in binary systems: The circumprimary case. Icarus 183, 193-206 (2006). https://doi.org/10.1016Zj.icarus.2006.01.022
4. Mayor M., Duquennoy A., Halbwachs J.-L., Mermilliod J.-C. CORAVEL Surveys to Study Binaries of Different Masses and Ages. Complementary Approaches to Double and Multiple Star Research, ASP Conference Series 32 IAU Colloquium 135, 73-81 (1992).
5. Larson R. B. Implications of Binary Properties for Theories of Star Formation. The Formation of Binary Stars, Proceedings of IAU Symp. 200, 93-106 (2001). https://doi.org/10.48550/arXiv.astro-ph/0006288
6. Shevchenko I. I. Dynamical Chaos in Planetary Systems. Springer Nature (2020). https://doi.org/10.1007/978-3-030-52144-8
7. Demidova T. V., Shevchenko 1.1. Spiral patterns in planetesimal circumbinary disks. Astrophys-ical Journal 805 (38), 8 (2015). https://doi.org/10.1088/0004-637X/805/1/38
8. Murray C.D., Dermott S. F. Solar ¡System Dynamics. Cambridge University Press (1999). [Rus. ed.: Murray C.D., Dermott S. F. Dinamika solnechnoi sistemy. Moscow, Fizmatlit Publ. (2009)].
9. Lithwick Y., Wu Y. Theory of Secular Chaos and Mercury's Orbit. Astrophysical Journal 739 (1), 17 (2011). https://doi.org/10.1088/0004-637X/739/V31
10. Rein H., Spiegel D. S. IAS15: a fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits. Monthly Notices of the Royal Astronomical Society 446 (2), 1424-1437 (2015). https://doi.org/10.1093/mnras/stu2164
11. Rein H., Liu S.-F. REBOUND: an open-source multi-purpose N-body code for collisional dynamics. Astronomy & Astrophysics 537 (A128), 10 (2012). https://doi.org/10.1051/0004-6361/201118085
12. Wisdom J. The resonance overlap criterion and the onset of stochastic behavior in the restricted three-body problem. Astronomical Journal 85, 1122-1133 (1980). https://doi.org/10.1086/112778
13. Duncan M., Quinn T., Tremaine S. The long-term evolution of orbits in the solar system. A mapping approach. Icarus 82 (2), 402-418 (1989). https://doi.org/10.1016/0019-1035(89)90047-X
14. Marov M. Ya., Shevchenko I.I. Exoplanets: Nature and models. Uspekhi fizicheskikh nauk 190, 897-932 (2020). https://doi.org/10.3367/UFNr.2019.10.038673 (In Russian) [Eng. transl.: Phys. Uspekhi 63, 837-871 (2020). https://doi.org/10.3367/UFNe.2019.10.038673].
Received: May 20, 2023 Revised: August 30, 2023 Accepted: August 31, 2023
Author's information:
Alena A. Berezina — st055707@student.spbu.ru