EduKG
← Назад

Динамика Солнечной системы и экзопланеты

Цели занятия

  • Итог

    Применять третий закон Кеплера T2=4π2GMa3T^2 = \dfrac{4\pi^2}{GM}a^3 для вычисления периода обращения по большой полуоси и обратно, а также для оценки массы центрального тела.

  • Итог

    Вычислять минимальную массу экзопланеты msin⁡im\sin i по амплитуде радиальных скоростей KK и периоду PP, используя функцию масс.

  • Итог

    Определять радиус планеты по глубине транзита δ=(Rp/R∗)2\delta = (R_p/R_*)^2, а по транзиту и радиальным скоростям вместе — среднюю плотность и сделать вывод о типе планеты.

Понять идею

Текст

Движение планет, их спутников и экзопланет определяется всемирным тяготением: F=GMmr2F = G\dfrac{Mm}{r^2}. Из этого закона выводятся законы Кеплера, и для круговой орбиты вокруг тела массы MM период связан с большой полуосью: T2=4π2GM a3.T^2 = \frac{4\pi^2}{GM}\,a^3. Обратите внимание: период зависит от aa и MM, но не от массы самой планеты, поэтому по наблюдаемым TT и aa можно найти массу центрального тела. Реальные системы не идеальны: взаимные возмущения планет, резонансы, приливы и миграция дают малые отклонения от кеплеровского движения — именно эти отклонения позволяют измерять массы планет Солнечной системы. Экзопланеты обнаруживают в основном косвенными методами: по доплеровскому смещению линий звезды (радиальные скорости), по ослаблению блеска при транзите, а также микролинзированием и прямыми изображениями. Метод радиальных скоростей даёт проекцию массы msin⁡im\sin i, транзит даёт отношение радиусов δ=(Rp/R∗)2\delta = (R_p/R_*)^2, а совместные наблюдения — плотность планеты и вывод о её составе. При этом методы чувствительны к разным планетам: горячие гиганты находятся легче, чем небольшие холодные планеты, поэтому известная выборка экзопланет смещена.

Разобранные примеры

Напрямую

Задача

Планета обращается по круговой орбите вокруг звезды массой M=1 M⊙M = 1\,M_\odot с большой полуосью a=1.5a = 1.5 а.е. Найдите период обращения TT в годах и в сутках.

Решение

Запишем третий закон Кеплера в удобной форме, где aa — в а.е., а масса — в массах Солнца: T=a3M/M⊙ лет.T = \sqrt{\frac{a^3}{M/M_\odot}}\ \text{лет}. Подставляем a=1.5a = 1.5 и M/M⊙=1M/M_\odot = 1: T=(1.5)3=3.375≈1.84 года.T = \sqrt{(1.5)^3} = \sqrt{3.375} \approx 1.84\ \text{года}. В сутках: 1.84×365≈6721.84 \times 365 \approx 672 сут. Проверка: для Земли a=1a = 1 а.е. и M=1 M⊙M = 1\,M_\odot дают ровно T=1T = 1 год, значит форма записи выбрана верно.

Применение

Задача

У звезды массой M=1 M⊙M = 1\,M_\odot спектральные линии периодически смещаются с амплитудой K=12K = 12 м/с и периодом P=4P = 4 сут. Оцените минимальную массу планеты msin⁡im\sin i.

Решение

Исходим из функции масс: m3sin⁡3i(M+m)2=PK32πG.\frac{m^3\sin^3 i}{(M+m)^2} = \frac{P K^3}{2\pi G}. Так как m≪Mm \ll M, принимаем M+m≈MM + m \approx M и получаем msin⁡i≈(PK3M22πG)1/3.m\sin i \approx \left(\frac{P K^3 M^2}{2\pi G}\right)^{1/3}. Переводим период: P=4×86400=3.456×105P = 4 \times 86400 = 3.456\times10^5 с. Считаем произведение в числителе: K3=123=1728K^3 = 12^3 = 1728, тогда PK3=3.456×105×1728=5.97×108P K^3 = 3.456\times10^5 \times 1728 = 5.97\times10^8; масса в квадрате M2=(1.989×1030)2=3.96×1060M^2 = (1.989\times10^{30})^2 = 3.96\times10^{60}; итого 2.36×10692.36\times10^{69}. Знаменатель: 2πG=4.19×10−102\pi G = 4.19\times10^{-10}. Дробь равна 5.63×10785.63\times10^{78}, кубический корень даёт msin⁡i≈1.8×1026 кг≈30 M⊕≈0.09 MЮп.m\sin i \approx 1.8\times10^{26}\ \text{кг} \approx 30\,M_\oplus \approx 0.09\,M_\text{Юп}. Это планета типа Нептуна. Полученное значение — нижняя оценка: при наклонении i<90∘i < 90^\circ истинная масса больше.

Перенос

Задача

Транзит планеты по диску звезды радиусом R∗=1 R⊙R_* = 1\,R_\odot уменьшает блеск на δ=1%\delta = 1\%. Метод радиальных скоростей при почти рёберном наклонении (i≈90∘i \approx 90^\circ) даёт массу m≈0.9 MЮпm \approx 0.9\,M_\text{Юп}. Найдите радиус и среднюю плотность планеты и решите, каменистая она или газообразная.

Решение

Радиус находим из глубины транзита: δ=(RpR∗)2⇒RpR∗=0.01=0.1.\delta = \left(\frac{R_p}{R_*}\right)^2 \Rightarrow \frac{R_p}{R_*} = \sqrt{0.01} = 0.1. Значит Rp=0.1 R⊙=0.1×6.96×108 м=6.96×107R_p = 0.1\,R_\odot = 0.1 \times 6.96\times10^8\ \text{м} = 6.96\times10^7 м, то есть примерно 1111 радиусов Земли — размер газового гиганта. Объём: V=43πRp3=43π(6.96×107)3≈1.41×1024 м3.V = \frac{4}{3}\pi R_p^3 = \frac{4}{3}\pi (6.96\times10^7)^3 \approx 1.41\times10^{24}\ \text{м}^3. Переводим массу: m=0.9×1.898×1027=1.71×1027m = 0.9 \times 1.898\times10^{27} = 1.71\times10^{27} кг. Плотность: ρ=mV=1.71×10271.41×1024≈1.2×103 кг/м3=1.2 г/см3.\rho = \frac{m}{V} = \frac{1.71\times10^{27}}{1.41\times10^{24}} \approx 1.2\times10^3\ \text{кг/м}^3 = 1.2\ \text{г/см}^3. Это близко к плотности Юпитера (1.331.33 г/см3^3) и намного меньше плотности камня (33–55 г/см3^3), поэтому планета газообразная, а не каменистая.

Практика

Текст

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

  1. Откройте онлайн-симулятор транзитов и кривых радиальных скоростей либо возьмите архивные данные какой-нибудь экзопланеты из открытой базы и постройте графики в таблице.
  2. Задайте звезду 1 M⊙1\,M_\odot и планету на круговой орбите; запишите полученные PP и aa.
  3. По формуле M=4π2a3GT2M = \dfrac{4\pi^2 a^3}{G T^2} восстановите массу звезды из «наблюдений» и сравните её с заданной.
  4. По глубине транзита найдите Rp/R∗R_p/R_*, по амплитуде KK найдите msin⁡im\sin i, вычислите плотность и решите, каменистая планета или газообразная.
  5. Меняя радиус и период планеты, отметьте, при каких параметрах сигнал падает ниже уровня шума — это наглядно показывает наблюдательную селекцию.

Безопасность. Никогда не смотрите на Солнце без специального фильтра; любые наблюдения ярких объектов ведите только через проверенный солнечный фильтр, установленный перед объективом.

Визуальная подсказка

Назначение

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

Представление

Схема из двух графиков с общей осью времени. Верхний — кривая радиальных скоростей vr(t)=Ksin⁡(2πt/P)v_r(t) = K\sin(2\pi t/P) с подписанными амплитудой KK и периодом PP. Нижний — кривая блеска с провалами глубиной δ=(Rp/R∗)2\delta = (R_p/R_*)^2 в моменты транзита; под одним из провалов схематично показаны звезда и проходящая по её диску планета. Моменты провалов совпадают с нулями синусоиды.

Примеры из жизни

Навигация GPS в автомобиле

Сценарий

Приёмник GPS в автомобиле использует сигналы со спутников, обращающихся на высоте около 20 200 km20\,200\ \text{km}. Каждый спутник совершает один оборот примерно за 11 h 58 min11\ \text{h}\ 58\ \text{min}, двигаясь со скоростью около 3.9 km/s3.9\ \text{km/s}. Приёмник измеряет время прохождения сигнала с точностью до нескольких наносекунд, чтобы определить положение автомобиля с точностью около 5 m5\ \text{m}.

Связь

Радиус орбиты и период спутника связаны третьим законом Кеплера, T2∝a3T^2 \propto a^3, который следует из закона всемирного тяготения Ньютона. Те же законы позволяют вычислять орбиты экзопланет по наблюдаемым периодам их транзитов. Время прохождения сигнала также использует скорость света, что напрямую связано с измерением расстояний в астрономии.

Приливы на побережье

Сценарий

В прибрежном городе разница между приливом и отливом обычно составляет около 2.0 m2.0\ \text{m}. Во время полнолуния, когда Солнце и Луна действуют вдоль одной линии, амплитуда прилива может возрасти до 3.0 m3.0\ \text{m} (сизигийный прилив). Луна находится на расстоянии 384 400 km384\,400\ \text{km}; её приливное воздействие на океаны Земли примерно вдвое больше солнечного, хотя Солнце гораздо массивнее, потому что приливная сила убывает как 1/r31/r^3.

Связь

Это применяет закон всемирного тяготения Ньютона, F=Gm1m2r2F = G \frac{m_1 m_2}{r^2}, и понятие приливных сил, которые возникают из-за разницы гравитационного ускорения в разных точках тела. Та же дифференциальная гравитация заставляет звезду «качаться» из-за обращающейся вокруг неё экзопланеты — сигнал, обнаруживаемый методом радиальных скоростей.

Наблюдение транзита экзопланеты

Сценарий

Тринадцатилетний школьник использует телескоп с диаметром 100 mm100\ \text{mm} и камеру смартфона, чтобы записывать яркость звезды HD 209458 в течение 44 часов. Яркость звезды падает примерно на 1.5%1.5\% в течение 33 часов, когда планета HD 209458 b проходит перед ней. Планета обращается с периодом 3.53.5 дня.

Связь

Глубина падения яркости даёт отношение радиуса планеты к радиусу звезды. Длительность транзита и орбитальный период вместе с третьим законом Кеплера позволяют найти большую полуось и наклонение орбиты. Это транзитный метод — один из основных способов обнаружения и characterization экзопланет.