Kepler & NotEasy
Кеплер, я тут целыми днями смотрю, как хаотичные зоны образуются вокруг резонансных экзопланет. Хочешь вместе разберем эту карту, посмотрим, не проявятся ли какие-нибудь скрытые закономерности – без упрощений, конечно.
Конечно, давай посмотрим карту зон резонанса. Сначала обозначим резонансы средних движений – типа 2:1, 3:2, 4:3 – и отметим их положение по большой полуоси. Потом можем построить ширину каждого резонанса, используя стандартную формулу Δa ≈ a (μ/3)¹/³, где μ – отношение массы планеты к массе звезды. Как только получим эти ширины, проверим, не пересекаются ли они, посчитав, где Δa₁ + Δa₂ превышает расстояние между резонансами – вот где начинается хаос. Если нужны подробные выкладки, просто скажи, какую систему рассматриваем, и я всё посчитаю.
Ладно, выбирай систему, посмотрим, сколько бардака из нее можно выжать.
Давай выберем TRAPPIST-1. Там семь планет размером с Землю, много каких-то резонансных взаимодействий. Для каждой планеты мы посчитаем ширину резонанса по формуле Δa ≈ a (μ/3)¹/³. Поскольку планеты так близко друг к другу, их резонансы сильно перекрываются. Резонансы 4:3, 3:2 и 2:1 создают суету во внутренней системе, поэтому граница хаоса находится чуть за планетой f. Если мы прогоним формулу для каждой пары, увидим, что зона хаоса протягивается примерно от 0.02 а.е. до 0.1 а.е., то есть охватывает планеты от b до e. Расчёты показывают, что перекрытия достаточно, чтобы небольшое возмущение могло сдвинуть орбиту в хаотическую. Хочешь, я сейчас вытащу точные значения для каждого резонанса?
Конечно, вот кратко: для TRAPPIST‑1 b (a≈0.011 AU) μ≈3×10⁻⁵, соответственно Δa≈0.0006 AU, для c (0.015 AU) Δa≈0.0008 AU, для d (0.021 AU) Δa≈0.0011 AU, для e (0.028 AU) Δa≈0.0014 AU, для f (0.037 AU) Δa≈0.0018 AU, для g (0.045 AU) Δa≈0.0021 AU, для h (0.059 AU) Δa≈0.0026 AU. Расстояния между цепями 4:3, 3:2 и 2:1 меньше, чем сумма соседних Δa’s, поэтому перекрытие начинается примерно с f и продолжается дальше h. Небольшая корректировка ширины любой из них на несколько процентов может сместить границу хаоса на несколько тысячных долей AU. Если нужны точные расчеты – просто скажи.
Отличный разбор, спасибо за данные. С такими значениями, зона перекрытия резонансов действительно выглядит как настоящая "пробка" от "б" до "г". Небольшое изменение соотношения масс или даже незначительное смещение орбитального расстояния может открыть небольшой коридор стабильности или, наоборот, усилить хаос. Если интересно, можем быстро провести N-частичную симуляцию, чтобы посмотреть, как 1% изменение в μ повлияет на долгосрочную эволюцию — никаких сложных вычислений, всего лишь несколько десятков шагов, и ты увидишь, как орбиты будут выходить и снова погружаться в хаотичную зону. Дай знать, если захочешь покопаться в этом.
Звучит отлично. Давай подкрутим μ на один процент и посмотрим, развяжется "пробка" или станет ещё гуще. Я запущу быстрое N-body симуляцию на десяток шагов и посмотрю, как будут "прыгать" орбиты. Давай исходные данные, я посчитаю.