На примере транзитной кривой короткопериодной экзопланеты оцените методы определения радиуса и плотности планеты и обсудите неопределённости при отсутствии допплеровских измерений скорости
Кратко — что даёт транзит и как на этой основе получить радиус и плотность, с примером и оценкой неопределённостей при отсутствии допплеровских измерений. Что извлекается из транзитной кривой - Глубина транзита: δ≃(RpR⋆)2\delta \simeq \left(\dfrac{R_p}{R_\star}\right)^2δ≃(R⋆Rp)2 → отношение радиусов RpR⋆=δ\dfrac{R_p}{R_\star}=\sqrt{\delta}R⋆Rp=δ. - Форма (длительность T14T_{14}T14, длительность погружения, скос bbb) даёт масштаб орбиты a/R⋆a/R_\stara/R⋆ и косвенно плотность звезды ρ⋆\rho_\starρ⋆ через закон Кеплера: ρ⋆≈3πGP2(aR⋆)3\displaystyle \rho_\star \approx \frac{3\pi}{G P^2}\left(\frac{a}{R_\star}\right)^3ρ⋆≈GP23π(R⋆a)3. - Для точного RpR_pRp требуется абсолютный радиус звезды R⋆R_\starR⋆: Rp=R⋆δR_p = R_\star\sqrt{\delta}Rp=R⋆δ. Простой числовой пример (короткопериодная планета) - Пусть P=1.5P=1.5P=1.5 суток и наблюдаем δ=0.01\delta=0.01δ=0.01. Тогда RpR⋆=0.1\dfrac{R_p}{R_\star}=0.1R⋆Rp=0.1. - Если по спектроскопии/парсингу Гайи R⋆=1.0 R⊙R_\star=1.0\,R_\odotR⋆=1.0R⊙, то Rp=0.1 R⊙≃1.0 RJR_p=0.1\,R_\odot\simeq 1.0\,R_JRp=0.1R⊙≃1.0RJ (приблизительно, RJ≈0.10045 R⊙R_J\approx0.10045\,R_\odotRJ≈0.10045R⊙). - Ошибка в R⋆R_\starR⋆ прямо переводится в ошибку в RpR_pRp: если σ(R⋆)/R⋆=10%\sigma(R_\star)/R_\star=10\%σ(R⋆)/R⋆=10%, то и σ(Rp)/Rp≈10%\sigma(R_p)/R_p\approx10\%σ(Rp)/Rp≈10%. Как получить массу и плотность без RV — методы и их ограничения 1. Массо-радиусные эмпирические соотношения (statistical M–R): - дают оценку MpM_pMp по RpR_pRp, но имеют большую разбросанность (особенно для переходных и малых планет). - Формула для плотности: ρp=Mp43πRp3\rho_p=\dfrac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp=34πRp3Mp. Погрешность σ(ρp)\sigma(\rho_p)σ(ρp) агрегирует σ(Mp)\sigma(M_p)σ(Mp) и σ(Rp)\sigma(R_p)σ(Rp) как ∼\sim∼ пропорционально σ(Mp)\sigma(M_p)σ(Mp) и 3σ(Rp)3\sigma(R_p)3σ(Rp). - Для газовых гигантов эмпирическая разбросанность в массе при заданном радиусе может быть фактором нескольких → плотность может быть неопределена на порядок. 2. TTV (транзитные временные вариации) в многопланетных системах: - При подвижном возмущении можно извлечь массы компаньонов с хорошей точностью, но нужен заметный TTV-сигнал и дополнительная планета в резонансе/близко. 3. Фотофазовая фотометрия: эллипсоидальные вариации и доплеровское ближение - Амплитуды дают зависимость от Mp/M⋆M_p/M_\starMp/M⋆ и (R⋆/a)3(R_\star/a)^3(R⋆/a)3. При благоприятном S/N (яркая цель, большая масса, короткий aaa) возможно ограничить массу без RV. - Примерно: амплитуда эллипсоидальных вариаций Aellip∼αMpM⋆(R⋆a)3A_{\rm ellip}\sim \alpha \dfrac{M_p}{M_\star}\left(\dfrac{R_\star}{a}\right)^3Aellip∼αM⋆Mp(aR⋆)3 — обычно мала (ppm–тысячи ppm), требует очень точной фотометрии (Kepler/TESS для ярких звёзд). 4. Вторичные минимуми/фазовые кривые: - Для очень горячих короткопериодных гигантов можно получить ограничения на тепловой вклад и иногда на массу через отражение/фазовую кривую, но это косвенно и модельно-зависимо. 5. Статистические/галактические приоритеты: - Использовать массовые приоритеты из популяций (например, вероятность того, что это суперземля/нептун/юпитер) — даёт широкий диапазон масс. Основные источники неопределённостей без RV - Неопределённость звёздных параметров (R⋆,M⋆R_\star, M_\starR⋆,M⋆): напрямую влияет на RpR_pRp и косвенно на массу при использовании a/R⋆a/R_\stara/R⋆. Типично σ(R⋆)/R⋆\sigma(R_\star)/R_\starσ(R⋆)/R⋆ = few % (с хорошей спектроскопией + Gaia) или ≳10%\gtrsim10\%≳10% иначе. - Разброс массово-радиусных соотношений: для малых планет — масса может варьировать на фактор ∼3−10\sim 3-10∼3−10; для гигантов — тоже большой разброс из‑за раздувания/надувания. - Смешивание (dilution, блендинг) и ложные сигналы (э клипсинг бинарные звёзды) могут занижать глубину → занижать RpR_pRp. - Неизвестная эксцентриситет eee: влияет на связь между a/R⋆a/R_\stara/R⋆ и длительностью транзита; неправильное предположение e=0e=0e=0 даёт систематическую ошибку в ρ⋆\rho_\starρ⋆ и далее в интерпретации. - Грациозный транзит (grazing) и неопределённый bbb: плохая оценка радиуса и больших ошибок в RpR_pRp. Количество ошибок — числовая иллюстрация - Допустим Rp=1 RJR_p=1\,R_JRp=1RJ с σ(Rp)=10%\sigma(R_p)=10\%σ(Rp)=10%. Если у нас нет прямой меры MpM_pMp и мы предполагаем MpM_pMp в диапазоне 0.3–3 MJ0.3\text{–}3\,M_J0.3–3MJ, то ρp/ρJ=Mp/MJ(Rp/RJ)3\rho_p/\rho_J = \dfrac{M_p/M_J}{(R_p/R_J)^3}ρp/ρJ=(Rp/RJ)3Mp/MJ. При Rp=1 RJR_p=1\,R_JRp=1RJ плотность варьирует в пределах ∼0.3–3\sim 0.3\text{–}3∼0.3–3 по отношению к ρJ\rho_JρJ — фактор ∼10\sim10∼10 разброса массы даёт почти порядок неопределённости в плотности. Даже при более узком предположении масс 0.5–2 MJM_JMJ вы получите фактор 444 в плотности. Рекомендации для минимизации неопределённости - Получить качественную спектроскопию и паралаксу (Gaia) для снижения ошибок R⋆,M⋆R_\star, M_\starR⋆,M⋆. - Искать TTV в длительных сериях наблюдений (если система мультипланетная). - Проверить и исправить блендинг/третье светило (high-resolution imaging). - Если нет RV, попытаться измерить эллипсоидальные/бимпинг/фазовые сигналы при высокой точности фотометрии и моделировании. - Всегда приводить плотность с учётом систематической погрешности от метода оценки массы; при отсутствии RV лучше указывать диапазон, а не единственное значение. Короткая итоговая формула - Радиус: Rp=R⋆δ\displaystyle R_p = R_\star\sqrt{\delta}Rp=R⋆δ. - Плотность (при известной массе): ρp=Mp43πRp3\displaystyle \rho_p=\frac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp=34πRp3Mp. Без независимых измерений MpM_pMp погрешности в ρp\rho_pρp обычно доминируют и могут достигать порядка единиц — факторов нескольких до порядка величины. Если нужно, могу привести расчёт конкретного числового примера с оценкой ошибок по заданным входным данным (P, δ, σ(R★), предположения по M_p).
Что извлекается из транзитной кривой
- Глубина транзита: δ≃(RpR⋆)2\delta \simeq \left(\dfrac{R_p}{R_\star}\right)^2δ≃(R⋆ Rp )2 → отношение радиусов RpR⋆=δ\dfrac{R_p}{R_\star}=\sqrt{\delta}R⋆ Rp =δ .
- Форма (длительность T14T_{14}T14 , длительность погружения, скос bbb) даёт масштаб орбиты a/R⋆a/R_\stara/R⋆ и косвенно плотность звезды ρ⋆\rho_\starρ⋆ через закон Кеплера:
ρ⋆≈3πGP2(aR⋆)3\displaystyle \rho_\star \approx \frac{3\pi}{G P^2}\left(\frac{a}{R_\star}\right)^3ρ⋆ ≈GP23π (R⋆ a )3.
- Для точного RpR_pRp требуется абсолютный радиус звезды R⋆R_\starR⋆ : Rp=R⋆δR_p = R_\star\sqrt{\delta}Rp =R⋆ δ .
Простой числовой пример (короткопериодная планета)
- Пусть P=1.5P=1.5P=1.5 суток и наблюдаем δ=0.01\delta=0.01δ=0.01. Тогда
RpR⋆=0.1\dfrac{R_p}{R_\star}=0.1R⋆ Rp =0.1.
- Если по спектроскопии/парсингу Гайи R⋆=1.0 R⊙R_\star=1.0\,R_\odotR⋆ =1.0R⊙ , то
Rp=0.1 R⊙≃1.0 RJR_p=0.1\,R_\odot\simeq 1.0\,R_JRp =0.1R⊙ ≃1.0RJ (приблизительно, RJ≈0.10045 R⊙R_J\approx0.10045\,R_\odotRJ ≈0.10045R⊙ ).
- Ошибка в R⋆R_\starR⋆ прямо переводится в ошибку в RpR_pRp : если σ(R⋆)/R⋆=10%\sigma(R_\star)/R_\star=10\%σ(R⋆ )/R⋆ =10%, то и σ(Rp)/Rp≈10%\sigma(R_p)/R_p\approx10\%σ(Rp )/Rp ≈10%.
Как получить массу и плотность без RV — методы и их ограничения
1. Массо-радиусные эмпирические соотношения (statistical M–R):
- дают оценку MpM_pMp по RpR_pRp , но имеют большую разбросанность (особенно для переходных и малых планет).
- Формула для плотности: ρp=Mp43πRp3\rho_p=\dfrac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp =34 πRp3 Mp . Погрешность σ(ρp)\sigma(\rho_p)σ(ρp ) агрегирует σ(Mp)\sigma(M_p)σ(Mp ) и σ(Rp)\sigma(R_p)σ(Rp ) как ∼\sim∼ пропорционально σ(Mp)\sigma(M_p)σ(Mp ) и 3σ(Rp)3\sigma(R_p)3σ(Rp ).
- Для газовых гигантов эмпирическая разбросанность в массе при заданном радиусе может быть фактором нескольких → плотность может быть неопределена на порядок.
2. TTV (транзитные временные вариации) в многопланетных системах:
- При подвижном возмущении можно извлечь массы компаньонов с хорошей точностью, но нужен заметный TTV-сигнал и дополнительная планета в резонансе/близко.
3. Фотофазовая фотометрия: эллипсоидальные вариации и доплеровское ближение
- Амплитуды дают зависимость от Mp/M⋆M_p/M_\starMp /M⋆ и (R⋆/a)3(R_\star/a)^3(R⋆ /a)3. При благоприятном S/N (яркая цель, большая масса, короткий aaa) возможно ограничить массу без RV.
- Примерно: амплитуда эллипсоидальных вариаций Aellip∼αMpM⋆(R⋆a)3A_{\rm ellip}\sim \alpha \dfrac{M_p}{M_\star}\left(\dfrac{R_\star}{a}\right)^3Aellip ∼αM⋆ Mp (aR⋆ )3 — обычно мала (ppm–тысячи ppm), требует очень точной фотометрии (Kepler/TESS для ярких звёзд).
4. Вторичные минимуми/фазовые кривые:
- Для очень горячих короткопериодных гигантов можно получить ограничения на тепловой вклад и иногда на массу через отражение/фазовую кривую, но это косвенно и модельно-зависимо.
5. Статистические/галактические приоритеты:
- Использовать массовые приоритеты из популяций (например, вероятность того, что это суперземля/нептун/юпитер) — даёт широкий диапазон масс.
Основные источники неопределённостей без RV
- Неопределённость звёздных параметров (R⋆,M⋆R_\star, M_\starR⋆ ,M⋆ ): напрямую влияет на RpR_pRp и косвенно на массу при использовании a/R⋆a/R_\stara/R⋆ . Типично σ(R⋆)/R⋆\sigma(R_\star)/R_\starσ(R⋆ )/R⋆ = few % (с хорошей спектроскопией + Gaia) или ≳10%\gtrsim10\%≳10% иначе.
- Разброс массово-радиусных соотношений: для малых планет — масса может варьировать на фактор ∼3−10\sim 3-10∼3−10; для гигантов — тоже большой разброс из‑за раздувания/надувания.
- Смешивание (dilution, блендинг) и ложные сигналы (э клипсинг бинарные звёзды) могут занижать глубину → занижать RpR_pRp .
- Неизвестная эксцентриситет eee: влияет на связь между a/R⋆a/R_\stara/R⋆ и длительностью транзита; неправильное предположение e=0e=0e=0 даёт систематическую ошибку в ρ⋆\rho_\starρ⋆ и далее в интерпретации.
- Грациозный транзит (grazing) и неопределённый bbb: плохая оценка радиуса и больших ошибок в RpR_pRp .
Количество ошибок — числовая иллюстрация
- Допустим Rp=1 RJR_p=1\,R_JRp =1RJ с σ(Rp)=10%\sigma(R_p)=10\%σ(Rp )=10%. Если у нас нет прямой меры MpM_pMp и мы предполагаем MpM_pMp в диапазоне 0.3–3 MJ0.3\text{–}3\,M_J0.3–3MJ , то
ρp/ρJ=Mp/MJ(Rp/RJ)3\rho_p/\rho_J = \dfrac{M_p/M_J}{(R_p/R_J)^3}ρp /ρJ =(Rp /RJ )3Mp /MJ .
При Rp=1 RJR_p=1\,R_JRp =1RJ плотность варьирует в пределах ∼0.3–3\sim 0.3\text{–}3∼0.3–3 по отношению к ρJ\rho_JρJ — фактор ∼10\sim10∼10 разброса массы даёт почти порядок неопределённости в плотности. Даже при более узком предположении масс 0.5–2 MJM_JMJ вы получите фактор 444 в плотности.
Рекомендации для минимизации неопределённости
- Получить качественную спектроскопию и паралаксу (Gaia) для снижения ошибок R⋆,M⋆R_\star, M_\starR⋆ ,M⋆ .
- Искать TTV в длительных сериях наблюдений (если система мультипланетная).
- Проверить и исправить блендинг/третье светило (high-resolution imaging).
- Если нет RV, попытаться измерить эллипсоидальные/бимпинг/фазовые сигналы при высокой точности фотометрии и моделировании.
- Всегда приводить плотность с учётом систематической погрешности от метода оценки массы; при отсутствии RV лучше указывать диапазон, а не единственное значение.
Короткая итоговая формула
- Радиус: Rp=R⋆δ\displaystyle R_p = R_\star\sqrt{\delta}Rp =R⋆ δ .
- Плотность (при известной массе): ρp=Mp43πRp3\displaystyle \rho_p=\frac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp =34 πRp3 Mp .
Без независимых измерений MpM_pMp погрешности в ρp\rho_pρp обычно доминируют и могут достигать порядка единиц — факторов нескольких до порядка величины.
Если нужно, могу привести расчёт конкретного числового примера с оценкой ошибок по заданным входным данным (P, δ, σ(R★), предположения по M_p).