Кейс: наблюдаем световую кривую транзита экзопланеты с глубиной 1% и периодом 10 суток у звезды типа G2V — какие параметры планеты можно оценить из этих данных и какие дополнительные наблюдения нужны для определения плотности
Можно оценить из транзитной кривой и заданных данных: - Радиус планеты RpR_pRp (относительно звезды): - глубина транзита δ\deltaδ даёт δ=(RpR∗)2\displaystyle \delta=\left(\frac{R_p}{R_*}\right)^2δ=(R∗Rp)2, отсюда RpR∗=δ\displaystyle \frac{R_p}{R_*}=\sqrt{\delta}R∗Rp=δ. - при δ=0.01\delta=0.01δ=0.01 и для звезды G2V (R∗≈R⊙R_*\approx R_\odotR∗≈R⊙) получаем Rp≈0.1 R∗≈0.1 R⊙∼1 RJR_p\approx0.1\,R_*\approx0.1\,R_\odot\sim1\,R_{\rm J}Rp≈0.1R∗≈0.1R⊙∼1RJ. - Полуось орбиты aaa (при известной массе звезды M∗M_*M∗) через 3-й закон Кеплера: - a3=GM∗P24π2\displaystyle a^3=\frac{G M_* P^2}{4\pi^2}a3=4π2GM∗P2. - при M∗≈M⊙M_*\approx M_\odotM∗≈M⊙ и P=10P=10P=10 сут. a≈0.09a\approx0.09a≈0.09 AU. - Температура равновесия (оценочно, при заданной альбедо ABA_BAB): - Teq=TeffR∗2a (1−AB)1/4\displaystyle T_{\rm eq}=T_{\rm eff}\sqrt{\frac{R_*}{2a}}\,(1-A_B)^{1/4}Teq=Teff2aR∗(1−AB)1/4 — для Teff∼5770T_{\rm eff}\sim5770Teff∼5770 K и AB∼0A_B\sim0AB∼0 даёт Teq∼900T_{\rm eq}\sim900Teq∼900 K (порядок). - Геометрические параметры орбиты и звезды из формы транзита: - по длительности транзита и длительностям вх/выхода можно оценить масштабированную полуось aR∗\displaystyle \frac{a}{R_*}R∗a, параметр сближения bbb и наклонение iii. Для круговой орбиты длительность TTT связана с этими величинами через T=Pπarcsin (R∗a(1+k)2−b2sini)\displaystyle T=\frac{P}{\pi}\arcsin\!\Big(\frac{R_*}{a}\frac{\sqrt{(1+k)^2-b^2}}{\sin i}\Big)T=πParcsin(aR∗sini(1+k)2−b2), где k=Rp/R∗k=R_p/R_*k=Rp/R∗. - из aR∗\displaystyle \frac{a}{R_*}R∗a и PPP можно получить среднюю плотность звезды: ρ∗=3πGP2(aR∗)3\displaystyle \rho_*=\frac{3\pi}{G P^2}\Big(\frac{a}{R_*}\Big)^3ρ∗=GP23π(R∗a)3. Что дополнительно нужно для определения плотности планеты ρp\rho_pρp: - Для плотности нужен массовый параметр: ρp=Mp43πRp3\displaystyle \rho_p=\frac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp=34πRp3Mp. Т.о. нужно измерить массу MpM_pMp. - Основной способ — радиальная скорость (RV): измерить полуширину KKK спектроскопически и получить MpM_pMp (поскольку транзит даёт i≈90∘i\approx90^\circi≈90∘, sini≈1\sin i\approx1sini≈1): K=(2πGP)1/3Mpsini(M∗+Mp)2/311−e2\displaystyle K=\Big(\frac{2\pi G}{P}\Big)^{1/3}\frac{M_p\sin i}{(M_*+M_p)^{2/3}}\frac{1}{\sqrt{1-e^2}}K=(P2πG)1/3(M∗+Mp)2/3Mpsini1−e21. Для P=10P=10P=10 сут. и Mp∼1 MJM_p\sim1\,M_{\rm J}Mp∼1MJ ожидаемый K∼90K\sim90K∼90 м/с (порядок), значит RV с точностью десятки м/с легко даст массу; для меньших масс нужна лучшая точность (м/с). - Альтернативы: транзитные вариации времени (TTV) при наличии соседних возмущающих планет — требуют многолетних, очень точных измерений времени транзитов. - Точные параметры звезды M∗M_*M∗ и R∗R_*R∗: - радиус влияет на RpR_pRp, масса — на перевод aaa и RV-оценку массы; нужны спектроскопия, SED+параллакс (Gaia) или астросейсмология для уменьшения систематической погрешности в плотности планеты. - Доп. проверки для корректности радиуса: - высокоразрешающая визуальная/адаптивная съёмка для исключения источников разбавления (blend/dilution), корректировка ламбда-темп. и лимб-дакаринга при фиттинге световой кривой. Коротко: из глубины и периода вы получите Rp/R∗R_p/R_*Rp/R∗, aaa (при известном M∗M_*M∗), a/R∗a/R_*a/R∗, bbb, iii и приблизительную TeqT_{\rm eq}Teq. Чтобы получить плотность планеты, необходимо измерить массу MpM_pMp (обычно RV или TTV) и улучшить знания о M∗M_*M∗ и R∗R_*R∗ (спектроскопия, параллакс/модели или астросейсмология).
- Радиус планеты RpR_pRp (относительно звезды):
- глубина транзита δ\deltaδ даёт δ=(RpR∗)2\displaystyle \delta=\left(\frac{R_p}{R_*}\right)^2δ=(R∗ Rp )2, отсюда RpR∗=δ\displaystyle \frac{R_p}{R_*}=\sqrt{\delta}R∗ Rp =δ .
- при δ=0.01\delta=0.01δ=0.01 и для звезды G2V (R∗≈R⊙R_*\approx R_\odotR∗ ≈R⊙ ) получаем Rp≈0.1 R∗≈0.1 R⊙∼1 RJR_p\approx0.1\,R_*\approx0.1\,R_\odot\sim1\,R_{\rm J}Rp ≈0.1R∗ ≈0.1R⊙ ∼1RJ .
- Полуось орбиты aaa (при известной массе звезды M∗M_*M∗ ) через 3-й закон Кеплера:
- a3=GM∗P24π2\displaystyle a^3=\frac{G M_* P^2}{4\pi^2}a3=4π2GM∗ P2 .
- при M∗≈M⊙M_*\approx M_\odotM∗ ≈M⊙ и P=10P=10P=10 сут. a≈0.09a\approx0.09a≈0.09 AU.
- Температура равновесия (оценочно, при заданной альбедо ABA_BAB ):
- Teq=TeffR∗2a (1−AB)1/4\displaystyle T_{\rm eq}=T_{\rm eff}\sqrt{\frac{R_*}{2a}}\,(1-A_B)^{1/4}Teq =Teff 2aR∗ (1−AB )1/4 — для Teff∼5770T_{\rm eff}\sim5770Teff ∼5770 K и AB∼0A_B\sim0AB ∼0 даёт Teq∼900T_{\rm eq}\sim900Teq ∼900 K (порядок).
- Геометрические параметры орбиты и звезды из формы транзита:
- по длительности транзита и длительностям вх/выхода можно оценить масштабированную полуось aR∗\displaystyle \frac{a}{R_*}R∗ a , параметр сближения bbb и наклонение iii. Для круговой орбиты длительность TTT связана с этими величинами через
T=Pπarcsin (R∗a(1+k)2−b2sini)\displaystyle T=\frac{P}{\pi}\arcsin\!\Big(\frac{R_*}{a}\frac{\sqrt{(1+k)^2-b^2}}{\sin i}\Big)T=πP arcsin(aR∗ sini(1+k)2−b2 ), где k=Rp/R∗k=R_p/R_*k=Rp /R∗ .
- из aR∗\displaystyle \frac{a}{R_*}R∗ a и PPP можно получить среднюю плотность звезды:
ρ∗=3πGP2(aR∗)3\displaystyle \rho_*=\frac{3\pi}{G P^2}\Big(\frac{a}{R_*}\Big)^3ρ∗ =GP23π (R∗ a )3.
Что дополнительно нужно для определения плотности планеты ρp\rho_pρp :
- Для плотности нужен массовый параметр: ρp=Mp43πRp3\displaystyle \rho_p=\frac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp =34 πRp3 Mp . Т.о. нужно измерить массу MpM_pMp .
- Основной способ — радиальная скорость (RV): измерить полуширину KKK спектроскопически и получить MpM_pMp (поскольку транзит даёт i≈90∘i\approx90^\circi≈90∘, sini≈1\sin i\approx1sini≈1):
K=(2πGP)1/3Mpsini(M∗+Mp)2/311−e2\displaystyle K=\Big(\frac{2\pi G}{P}\Big)^{1/3}\frac{M_p\sin i}{(M_*+M_p)^{2/3}}\frac{1}{\sqrt{1-e^2}}K=(P2πG )1/3(M∗ +Mp )2/3Mp sini 1−e2 1 .
Для P=10P=10P=10 сут. и Mp∼1 MJM_p\sim1\,M_{\rm J}Mp ∼1MJ ожидаемый K∼90K\sim90K∼90 м/с (порядок), значит RV с точностью десятки м/с легко даст массу; для меньших масс нужна лучшая точность (м/с).
- Альтернативы: транзитные вариации времени (TTV) при наличии соседних возмущающих планет — требуют многолетних, очень точных измерений времени транзитов.
- Точные параметры звезды M∗M_*M∗ и R∗R_*R∗ :
- радиус влияет на RpR_pRp , масса — на перевод aaa и RV-оценку массы; нужны спектроскопия, SED+параллакс (Gaia) или астросейсмология для уменьшения систематической погрешности в плотности планеты.
- Доп. проверки для корректности радиуса:
- высокоразрешающая визуальная/адаптивная съёмка для исключения источников разбавления (blend/dilution), корректировка ламбда-темп. и лимб-дакаринга при фиттинге световой кривой.
Коротко: из глубины и периода вы получите Rp/R∗R_p/R_*Rp /R∗ , aaa (при известном M∗M_*M∗ ), a/R∗a/R_*a/R∗ , bbb, iii и приблизительную TeqT_{\rm eq}Teq . Чтобы получить плотность планеты, необходимо измерить массу MpM_pMp (обычно RV или TTV) и улучшить знания о M∗M_*M∗ и R∗R_*R∗ (спектроскопия, параллакс/модели или астросейсмология).