Дан транзитный световой кривой звезды и её спектроскопические наблюдения: как по глубине и форме транзита, временным вариациям (TTV) и потоку поглощения в линиях возбуждённых элементов можно реконструировать радиус, плотность, возможную атмосферу и наличие дополнительных планет в системе; приведите пример расчёта радиуса планеты для звезды солнечного типа и глубины транзита 0,1 %
Кратко — по каким наблюдаемым параметрам что можно получить и как это вычислять; затем пример расчёта радиуса. 1) Радиус планеты (Rp) - Основная формула из глубины транзита: глубина δ\deltaδ примерно равна отношению площадей: δ≈(RpR∗)2\delta \approx \left(\frac{R_p}{R_*}\right)^2δ≈(R∗Rp)2. Отсюда RpR∗=δ,Rp=R∗δ.
\frac{R_p}{R_*} = \sqrt{\delta},\qquad R_p = R_* \sqrt{\delta}. R∗Rp=δ,Rp=R∗δ.
- Форма транзита (ингресс/эгресс, длительность) даёт a/R∗a/R_*a/R∗ и параметр смещения bbb. Для круговой орбиты aR∗=(1+k)2−b2sin (πTdur/P),k=RpR∗,
\frac{a}{R_*} = \frac{\sqrt{(1+k)^2-b^2}}{\sin\!\left(\pi T_\mathrm{dur}/P\right)},\quad k=\frac{R_p}{R_*}, R∗a=sin(πTdur/P)(1+k)2−b2,k=R∗Rp,
где TdurT_\mathrm{dur}Tdur — длительность транзита, PPP — период. Из a/R∗a/R_*a/R∗ можно получить плотность звезды: ρ∗≈3πGP2(aR∗)3.
\rho_* \approx \frac{3\pi}{G P^2}\left(\frac{a}{R_*}\right)^3. ρ∗≈GP23π(R∗a)3.
Поэтому из точной формы транзита можно проверить и улучшить оценку R∗R_*R∗ и, соответственно, RpR_pRp. 2) Масса планеты и плотность (ρp\rho_pρp) - Массу дают: а) радиальные скорости (RV) — полуширина кривой: K=(2πGP)1/3Mpsini(M∗+Mp)2/3,
K=\left(\frac{2\pi G}{P}\right)^{1/3}\frac{M_p\sin i}{(M_*+M_p)^{2/3}}, K=(P2πG)1/3(M∗+Mp)2/3Mpsini,
для транзитной системы sini≈1\sin i\approx1sini≈1. Решив относительно MpM_pMp получаем массу. б) альтернативно — TTVs: временные вариации моментов транзита позволяют оценить массу возмущающей планеты через моделирование N--тельной динамики. - Плотность: ρp=Mp43πRp3.
\rho_p = \frac{M_p}{\tfrac{4}{3}\pi R_p^3}. ρp=34πRp3Mp. 3) Атмосфера (трансмиссионная спектроскопия) - В спектре во время транзита глубина в линиях или в полосах увеличивается из–за поглощения в атмосфере. Приближённая амплитуда спектрального сигнала: Δδλ≈(Rp+NHR∗)2−(RpR∗)2≈2RpNHR∗2,
\Delta \delta_\lambda \approx \left(\frac{R_p+N H}{R_*}\right)^2 - \left(\frac{R_p}{R_*}\right)^2 \approx 2\frac{R_p N H}{R_*^2}, Δδλ≈(R∗Rp+NH)2−(R∗Rp)2≈2R∗2RpNH,
где HHH — шкала высоты атмосферы, N∼3 − 7N\sim 3\!-\!7N∼3−7 — число шкал, видимых в линии. - Шкала высоты: H=kBTμmpg,g=GMpRp2,
H=\frac{k_B T}{\mu m_p g},\qquad g=\frac{G M_p}{R_p^2}, H=μmpgkBT,g=Rp2GMp,
где TTT — температура на терминирующей слое, μ\muμ — средняя молярная масса (в а. е.), mpm_pmp — масса протона. - По глубинам и форме линий (Na, K, Hα\alphaα, He 10830 Å и др.) можно оценить состав (легкая/тяжёлая атмосфера), наличие расширенной/уходящей атмосферы (широкие, смещённые вектора, асимметрия), скорость ветров и истечение (через профили и доплеровские сдвиги). 4) TTV и поиск дополнительных планет - Сдвиги моментов транзита (TTV) указывают на гравитационные возмущения соседних планет. Амплитуда и период TTV зависят от масс и близости к резонансу; по серии измерений проводится N‑тельная подгонка для извлечения масс и орбит возмущающих тел. - Дополнительно: долгосрочные изменения длительности/склонения (TDV) могут указывать на наклонённые компаньоны или прецессию. 5) Ограничения и систематики - Лимб‑дамминг, активные участки звезды (споты), систематические ошибки на RV и на времени наблюдений влияют на оценки. Нужна совместная модель световой кривой + RV + спектроскопии. Пример: звезда солнечного типа и глубина транзита δ=0.1%\delta=0.1\%δ=0.1% - Запишем δ= 0.1% = 0.001\delta=\;0.1\%\;=\;0.001δ=0.1%=0.001. - Отсюда отношение радиусов: k=RpR∗=δ=0.001≈0.03162.
k=\frac{R_p}{R_*}=\sqrt{\delta}=\sqrt{0.001}\approx 0.03162. k=R∗Rp=δ=0.001≈0.03162.
- Для звезды типа Солнца положим R∗=R⊙R_*=R_\odotR∗=R⊙. Тогда Rp=R⊙0.001≈0.03162 R⊙.
R_p = R_\odot\sqrt{0.001}\approx 0.03162\,R_\odot. Rp=R⊙0.001≈0.03162R⊙.
Подставим численно (R⊙=696 340 kmR_\odot=696\,340\ \mathrm{km}R⊙=696340km, R⊕=6 371 kmR_\oplus=6\,371\ \mathrm{km}R⊕=6371km, RJ=71 492 kmR_J=71\,492\ \mathrm{km}RJ=71492km): Rp≈0.03162×696 340 km≈22 020 km.
R_p \approx 0.03162\times 696\,340\ \mathrm{km}\approx 22\,020\ \mathrm{km}. Rp≈0.03162×696340km≈22020km.
В земных и юпитерианских радиусах: Rp≈22 0206 371 R⊕≈3.46 R⊕,Rp≈22 02071 492 RJ≈0.308 RJ.
R_p\approx \frac{22\,020}{6\,371}\,R_\oplus\approx 3.46\,R_\oplus,\qquad R_p\approx \frac{22\,020}{71\,492}\,R_J\approx 0.308\,R_J. Rp≈637122020R⊕≈3.46R⊕,Rp≈7149222020RJ≈0.308RJ.
Итого: при глубине δ=0.1%\delta=0.1\%δ=0.1% у солнечного аналога радиус планеты примерно ∼3.5 R⊕\sim 3.5\ R_\oplus∼3.5R⊕ (супер‑Земля/мини‑Нептун). Если нужна оценка плотности — требуется масса MpM_pMp (из RV или TTV). После получения MpM_pMp вычисляете ρp=Mp43πRp3\rho_p=\dfrac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp=34πRp3Mp и по значению ρp\rho_pρp и по признакам в трансмиссионной спектре делаете выводы об атмосфере (газовый/скалистый, наличие лёгких компонентов, расширенная атмосфера). Если хотите, могу привести пример полного расчёта массы из измеренной RV‑полуширины KKK и затем плотности для конкретных чисел.
1) Радиус планеты (Rp)
- Основная формула из глубины транзита: глубина δ\deltaδ примерно равна отношению площадей: δ≈(RpR∗)2\delta \approx \left(\frac{R_p}{R_*}\right)^2δ≈(R∗ Rp )2. Отсюда
RpR∗=δ,Rp=R∗δ. \frac{R_p}{R_*} = \sqrt{\delta},\qquad R_p = R_* \sqrt{\delta}.
R∗ Rp =δ ,Rp =R∗ δ . - Форма транзита (ингресс/эгресс, длительность) даёт a/R∗a/R_*a/R∗ и параметр смещения bbb. Для круговой орбиты
aR∗=(1+k)2−b2sin (πTdur/P),k=RpR∗, \frac{a}{R_*} = \frac{\sqrt{(1+k)^2-b^2}}{\sin\!\left(\pi T_\mathrm{dur}/P\right)},\quad k=\frac{R_p}{R_*},
R∗ a =sin(πTdur /P)(1+k)2−b2 ,k=R∗ Rp , где TdurT_\mathrm{dur}Tdur — длительность транзита, PPP — период. Из a/R∗a/R_*a/R∗ можно получить плотность звезды:
ρ∗≈3πGP2(aR∗)3. \rho_* \approx \frac{3\pi}{G P^2}\left(\frac{a}{R_*}\right)^3.
ρ∗ ≈GP23π (R∗ a )3. Поэтому из точной формы транзита можно проверить и улучшить оценку R∗R_*R∗ и, соответственно, RpR_pRp .
2) Масса планеты и плотность (ρp\rho_pρp )
- Массу дают: а) радиальные скорости (RV) — полуширина кривой:
K=(2πGP)1/3Mpsini(M∗+Mp)2/3, K=\left(\frac{2\pi G}{P}\right)^{1/3}\frac{M_p\sin i}{(M_*+M_p)^{2/3}},
K=(P2πG )1/3(M∗ +Mp )2/3Mp sini , для транзитной системы sini≈1\sin i\approx1sini≈1. Решив относительно MpM_pMp получаем массу. б) альтернативно — TTVs: временные вариации моментов транзита позволяют оценить массу возмущающей планеты через моделирование N--тельной динамики.
- Плотность:
ρp=Mp43πRp3. \rho_p = \frac{M_p}{\tfrac{4}{3}\pi R_p^3}.
ρp =34 πRp3 Mp .
3) Атмосфера (трансмиссионная спектроскопия)
- В спектре во время транзита глубина в линиях или в полосах увеличивается из–за поглощения в атмосфере. Приближённая амплитуда спектрального сигнала:
Δδλ≈(Rp+NHR∗)2−(RpR∗)2≈2RpNHR∗2, \Delta \delta_\lambda \approx \left(\frac{R_p+N H}{R_*}\right)^2 - \left(\frac{R_p}{R_*}\right)^2 \approx 2\frac{R_p N H}{R_*^2},
Δδλ ≈(R∗ Rp +NH )2−(R∗ Rp )2≈2R∗2 Rp NH , где HHH — шкала высоты атмосферы, N∼3 − 7N\sim 3\!-\!7N∼3−7 — число шкал, видимых в линии.
- Шкала высоты:
H=kBTμmpg,g=GMpRp2, H=\frac{k_B T}{\mu m_p g},\qquad g=\frac{G M_p}{R_p^2},
H=μmp gkB T ,g=Rp2 GMp , где TTT — температура на терминирующей слое, μ\muμ — средняя молярная масса (в а. е.), mpm_pmp — масса протона.
- По глубинам и форме линий (Na, K, Hα\alphaα, He 10830 Å и др.) можно оценить состав (легкая/тяжёлая атмосфера), наличие расширенной/уходящей атмосферы (широкие, смещённые вектора, асимметрия), скорость ветров и истечение (через профили и доплеровские сдвиги).
4) TTV и поиск дополнительных планет
- Сдвиги моментов транзита (TTV) указывают на гравитационные возмущения соседних планет. Амплитуда и период TTV зависят от масс и близости к резонансу; по серии измерений проводится N‑тельная подгонка для извлечения масс и орбит возмущающих тел.
- Дополнительно: долгосрочные изменения длительности/склонения (TDV) могут указывать на наклонённые компаньоны или прецессию.
5) Ограничения и систематики
- Лимб‑дамминг, активные участки звезды (споты), систематические ошибки на RV и на времени наблюдений влияют на оценки. Нужна совместная модель световой кривой + RV + спектроскопии.
Пример: звезда солнечного типа и глубина транзита δ=0.1%\delta=0.1\%δ=0.1%
- Запишем δ= 0.1% = 0.001\delta=\;0.1\%\;=\;0.001δ=0.1%=0.001.
- Отсюда отношение радиусов:
k=RpR∗=δ=0.001≈0.03162. k=\frac{R_p}{R_*}=\sqrt{\delta}=\sqrt{0.001}\approx 0.03162.
k=R∗ Rp =δ =0.001 ≈0.03162. - Для звезды типа Солнца положим R∗=R⊙R_*=R_\odotR∗ =R⊙ . Тогда
Rp=R⊙0.001≈0.03162 R⊙. R_p = R_\odot\sqrt{0.001}\approx 0.03162\,R_\odot.
Rp =R⊙ 0.001 ≈0.03162R⊙ . Подставим численно (R⊙=696 340 kmR_\odot=696\,340\ \mathrm{km}R⊙ =696340 km, R⊕=6 371 kmR_\oplus=6\,371\ \mathrm{km}R⊕ =6371 km, RJ=71 492 kmR_J=71\,492\ \mathrm{km}RJ =71492 km):
Rp≈0.03162×696 340 km≈22 020 km. R_p \approx 0.03162\times 696\,340\ \mathrm{km}\approx 22\,020\ \mathrm{km}.
Rp ≈0.03162×696340 km≈22020 km. В земных и юпитерианских радиусах:
Rp≈22 0206 371 R⊕≈3.46 R⊕,Rp≈22 02071 492 RJ≈0.308 RJ. R_p\approx \frac{22\,020}{6\,371}\,R_\oplus\approx 3.46\,R_\oplus,\qquad
R_p\approx \frac{22\,020}{71\,492}\,R_J\approx 0.308\,R_J.
Rp ≈637122020 R⊕ ≈3.46R⊕ ,Rp ≈7149222020 RJ ≈0.308RJ . Итого: при глубине δ=0.1%\delta=0.1\%δ=0.1% у солнечного аналога радиус планеты примерно ∼3.5 R⊕\sim 3.5\ R_\oplus∼3.5 R⊕ (супер‑Земля/мини‑Нептун).
Если нужна оценка плотности — требуется масса MpM_pMp (из RV или TTV). После получения MpM_pMp вычисляете ρp=Mp43πRp3\rho_p=\dfrac{M_p}{\tfrac{4}{3}\pi R_p^3}ρp =34 πRp3 Mp и по значению ρp\rho_pρp и по признакам в трансмиссионной спектре делаете выводы об атмосфере (газовый/скалистый, наличие лёгких компонентов, расширенная атмосфера).
Если хотите, могу привести пример полного расчёта массы из измеренной RV‑полуширины KKK и затем плотности для конкретных чисел.