Предложите модель аккреции вокруг молодой протозвезды, объясните механизм формирования планетезималей и оцените влияние радиации и магнитных полей на миграцию формирующихся планет; какие наблюдения ALMA помогут подтвердить модель?
Краткая модель аккреции (концептуальная схема) - Диск типа α (внутренняя аккреция + внешний поток пепла): поверхность плотности и температура Σ(r)=Σ0(rr0)−p,T(r)=T0(rr0)−q\displaystyle \Sigma(r)=\Sigma_0\left(\frac{r}{r_0}\right)^{-p},\qquad T(r)=T_0\left(\frac{r}{r_0}\right)^{-q}Σ(r)=Σ0(r0r)−p,T(r)=T0(r0r)−q
(обычно p∼1 − 1.5, q∼0.5 − 0.75p\sim1\!-\!1.5,\; q\sim0.5\!-\!0.75p∼1−1.5,q∼0.5−0.75). Вязкость задаётся ν=αcsH\nu=\alpha c_s Hν=αcsH, H=cs/ΩH=c_s/\OmegaH=cs/Ω. - Зоны: активная MRI-зона в верхних слоях, «мертвая» зона (низкая ионизация, малое α\alphaα) в средней плоскости → образование градиентов давления и петель (zonal flows). - Внешние механизмы: магнитно-дувные ветры (MHD winds) и фотоэвапорация формируют внутр. пустоты и изменяют массовый расход. Механизм образования планетезималей (последовательность) 1. Рост и отлёты пыли: фракционный рост до пепла/гранул, затем радиальный дрейф со скоростью vr≈−2ηvKStSt2+1,\displaystyle v_r\approx -2\eta v_K\frac{\mathrm{St}}{\mathrm{St}^2+1},vr≈−2ηvKSt2+1St,
где St=tsΩ\mathrm{St}=t_s\OmegaSt=tsΩ — число Стокса, η∼(cs/vK)2\eta\sim (c_s/v_K)^2η∼(cs/vK)2 — параметр подпружинивания. 2. Концентрация в ловушках давления: локальные максимумы давления (пояса/щели/вихри) останавливают дрейф и повышают локальный solid-to-gas. 3. Стриминговая неустойчивость (SI) при достаточной металличности/концентрации: требуется локальная доля твёрдых частиц ZZZ выше критической (Zcrit∼0.02 − 0.03Z_{\rm crit}\sim 0.02\!-\!0.03Zcrit∼0.02−0.03) и St∼10−2−1 \mathrm{St}\sim10^{-2}-1St∼10−2−1. SI приводит к образованию плотных сгустков, которые гравитационно коллапсируют в планетезимали (размеры ~10–100 km). 4. Дальнейший рост: пепел-аккреция (pebble accretion) на ядра; скорость аккреции быстро растёт при оптимальном St\mathrm{St}St. Ключевые количественные условия - Устойчивость к гравитационному коллапсу: Toomre Q=csΩπGΣ\displaystyle Q=\frac{c_s\Omega}{\pi G\Sigma}Q=πGΣcsΩ, планетообразование через гравитац. нестабильность требует Q≲1.5−2Q\lesssim1.5-2Q≲1.5−2. - Время радиального дрейфа для частиц St≪1\mathrm{St}\ll1St≪1: τdrift∼r∣vr∣∼1ηΩ1St\displaystyle \tau_{\rm drift}\sim \frac{r}{|v_r|}\sim \frac{1}{\eta\Omega}\frac{1}{\mathrm{St}}τdrift∼∣vr∣r∼ηΩ1St1. Влияние радиации и магнитных полей на миграцию формирующихся планет - Общая миграция: - Type I (меньшие массы): типичное время τI∼(M∗Mp)(M∗Σr2)(hr)2Ω−1.\displaystyle \tau_{\rm I}\sim\left(\frac{M_*}{M_p}\right)\left(\frac{M_*}{\Sigma r^2}\right)\left(\frac{h}{r}\right)^2\Omega^{-1}.τI∼(MpM∗)(Σr2M∗)(rh)2Ω−1.
- Type II (гигантские планеты, gap opening): миграция привязана к вязкому времени τII∼r2/ν\tau_{\rm II}\sim r^2/\nuτII∼r2/ν. - Эффекты излучения (звездное излучение, фотоэвапорация): - Нагрев и фотоэвапорация уменьшают поверхностную плотность газа Σ\SigmaΣ, повышают h/rh/rh/r в ближней части и сокращают миграционные торки → миграция замедляется или останавливается при образовании внутренней щели. Типичные массовые потери: M˙wind∼10−9−10−8 M⊙ yr−1\dot M_{\rm wind}\sim10^{-9}-10^{-8}\,M_\odot\,\mathrm{yr}^{-1}M˙wind∼10−9−10−8M⊙yr−1 (зависит от UV/X‑ray). - Радиационное нагревание меняет градиенты плотности/температуры, что может менять знак короткодиапазонного (corotation) торка и приводить к «пасткам» миграции (полевая точка останова). - Магнитные поля: - MRI-турбулентность задаёт эффективную α\alphaα (активная зона α∼10−3−10−2\alpha\sim10^{-3}-10^{-2}α∼10−3−10−2, мёртвая зона α≲10−4\alpha\lesssim10^{-4}α≲10−4) → это определяет ν\nuν и τII\tau_{\rm II}τII. - Наличие «мертвых зон» и магнитных зон (zonal flows) породит постоянные давления-ловушки → планеты/ядра могут остановиться или мигрировать медленнее. - MHD-ветры отводят угловой момент, по сути уменьшают плотность газа и могут приводить к внутреннему «перекрытию» миграции (меньше газа → слабее торки). - Нон-идеальная MHD (амбиполярное рассеяние, Холл-эффект) меняют структуру турбулентности и поэтому эффективность коротационных торков; в некоторых условиях магнитные силы могут вызвать нерегулярные ускоренные дрейфы (случайные «киксы»). - Итог: сильное излучение и сильные магнитные эффекты обычно замедляют или локально останавливают миграцию (через уменьшение Σ\SigmaΣ, создание ловушек давления или уменьшение торков); в активных турбулентных зонах миграция может быть stochastic. Наблюдения ALMA, подтверждающие модель (что именно смотреть) 1. Пыльные кольца/щели и асимметрии (continuum high‑res, λ∼0.8–3 mm): - Наличие узких колец/клиней и ярких «лис» (crescents) → давление-ловушки/вихри. - Многодорожечные карты (разные длины волны) дают карту спектрального индекса αν\alpha_\nuαν — снижение αν\alpha_\nuαν указывает на рост зерен (большие частицы, St↑). 2. Радиационно/газовое распределение (CO, C18O, HCO+): - Сопоставление газа и пыли: зоны с дефицитом газа при наличии плотной пыли указывают на фильтрацию и gap-opening. - Кинематические особенности: «kinks» и локальные отклонения от Кеплеровского потока → наличие планет и невязкости/локальной массы. 3. Турбулентность и α: - Линии ширины (CO isotopologues, CS) для оценки turbulent linewidth → ограничение α\alphaα. 4. Ионизация и магнитные поля: - Линии-ионы (HCO+, N2H+, CN) как трейсеры ионизации; дефицит/номенклатура указывает на мёртвую зону. - Поляризация континуума и линий (dust polarization, aligned grains, Goldreich-Kylafis) для карты магнитного поля (морфология B‑поля). 5. Доказательства ветров/выбросов: - Высоко-скоростные крылья в CO → диск‑ветер (MHD wind / photoevaporative flow). 6. Пространственно‑разрешённые карты размера зерен: - Радіальные профили размера зерен (multi-λ) показывают концентрацию пепла в трапах; сопоставление с моделями St и τ_drift. 7. Индикаторы стриминговой неустойчивости / перегрузки: - Очень высокая локальная dust-to-gas (через continuum / C18O) → условие для SI (ZlocalZ_{\rm local}Zlocal превышает ZcritZ_{\rm crit}Zcrit). Коротко: ожидаемые ALMA-признаки модели — совпадение пылевых колец/асимметрий с локальными изменениями газа и кинематикой, снижение спектрального индекса в ловушках, малые скорости турбулентности в «мертвых зонах», следы ветров и поляризация, соответствующая конфигурация магнитного поля. Эти данные позволяют связать: (i) места концентрации пепла → планетезимали, (ii) уменьшение газа/градиенты → влияние радиации, (iii) морфологию B‑поля и уровни турбулентности → влияние магнитных эффектов на миграцию.
- Диск типа α (внутренняя аккреция + внешний поток пепла): поверхность плотности и температура
Σ(r)=Σ0(rr0)−p,T(r)=T0(rr0)−q\displaystyle \Sigma(r)=\Sigma_0\left(\frac{r}{r_0}\right)^{-p},\qquad
T(r)=T_0\left(\frac{r}{r_0}\right)^{-q}Σ(r)=Σ0 (r0 r )−p,T(r)=T0 (r0 r )−q (обычно p∼1 − 1.5, q∼0.5 − 0.75p\sim1\!-\!1.5,\; q\sim0.5\!-\!0.75p∼1−1.5,q∼0.5−0.75). Вязкость задаётся
ν=αcsH\nu=\alpha c_s Hν=αcs H, H=cs/ΩH=c_s/\OmegaH=cs /Ω.
- Зоны: активная MRI-зона в верхних слоях, «мертвая» зона (низкая ионизация, малое α\alphaα) в средней плоскости → образование градиентов давления и петель (zonal flows).
- Внешние механизмы: магнитно-дувные ветры (MHD winds) и фотоэвапорация формируют внутр. пустоты и изменяют массовый расход.
Механизм образования планетезималей (последовательность)
1. Рост и отлёты пыли: фракционный рост до пепла/гранул, затем радиальный дрейф со скоростью
vr≈−2ηvKStSt2+1,\displaystyle v_r\approx -2\eta v_K\frac{\mathrm{St}}{\mathrm{St}^2+1},vr ≈−2ηvK St2+1St , где St=tsΩ\mathrm{St}=t_s\OmegaSt=ts Ω — число Стокса, η∼(cs/vK)2\eta\sim (c_s/v_K)^2η∼(cs /vK )2 — параметр подпружинивания.
2. Концентрация в ловушках давления: локальные максимумы давления (пояса/щели/вихри) останавливают дрейф и повышают локальный solid-to-gas.
3. Стриминговая неустойчивость (SI) при достаточной металличности/концентрации: требуется локальная доля твёрдых частиц ZZZ выше критической (Zcrit∼0.02 − 0.03Z_{\rm crit}\sim 0.02\!-\!0.03Zcrit ∼0.02−0.03) и St∼10−2−1 \mathrm{St}\sim10^{-2}-1St∼10−2−1. SI приводит к образованию плотных сгустков, которые гравитационно коллапсируют в планетезимали (размеры ~10–100 km).
4. Дальнейший рост: пепел-аккреция (pebble accretion) на ядра; скорость аккреции быстро растёт при оптимальном St\mathrm{St}St.
Ключевые количественные условия
- Устойчивость к гравитационному коллапсу: Toomre
Q=csΩπGΣ\displaystyle Q=\frac{c_s\Omega}{\pi G\Sigma}Q=πGΣcs Ω , планетообразование через гравитац. нестабильность требует Q≲1.5−2Q\lesssim1.5-2Q≲1.5−2.
- Время радиального дрейфа для частиц St≪1\mathrm{St}\ll1St≪1:
τdrift∼r∣vr∣∼1ηΩ1St\displaystyle \tau_{\rm drift}\sim \frac{r}{|v_r|}\sim \frac{1}{\eta\Omega}\frac{1}{\mathrm{St}}τdrift ∼∣vr ∣r ∼ηΩ1 St1 .
Влияние радиации и магнитных полей на миграцию формирующихся планет
- Общая миграция:
- Type I (меньшие массы): типичное время
τI∼(M∗Mp)(M∗Σr2)(hr)2Ω−1.\displaystyle \tau_{\rm I}\sim\left(\frac{M_*}{M_p}\right)\left(\frac{M_*}{\Sigma r^2}\right)\left(\frac{h}{r}\right)^2\Omega^{-1}.τI ∼(Mp M∗ )(Σr2M∗ )(rh )2Ω−1. - Type II (гигантские планеты, gap opening): миграция привязана к вязкому времени τII∼r2/ν\tau_{\rm II}\sim r^2/\nuτII ∼r2/ν.
- Эффекты излучения (звездное излучение, фотоэвапорация):
- Нагрев и фотоэвапорация уменьшают поверхностную плотность газа Σ\SigmaΣ, повышают h/rh/rh/r в ближней части и сокращают миграционные торки → миграция замедляется или останавливается при образовании внутренней щели. Типичные массовые потери: M˙wind∼10−9−10−8 M⊙ yr−1\dot M_{\rm wind}\sim10^{-9}-10^{-8}\,M_\odot\,\mathrm{yr}^{-1}M˙wind ∼10−9−10−8M⊙ yr−1 (зависит от UV/X‑ray).
- Радиационное нагревание меняет градиенты плотности/температуры, что может менять знак короткодиапазонного (corotation) торка и приводить к «пасткам» миграции (полевая точка останова).
- Магнитные поля:
- MRI-турбулентность задаёт эффективную α\alphaα (активная зона α∼10−3−10−2\alpha\sim10^{-3}-10^{-2}α∼10−3−10−2, мёртвая зона α≲10−4\alpha\lesssim10^{-4}α≲10−4) → это определяет ν\nuν и τII\tau_{\rm II}τII .
- Наличие «мертвых зон» и магнитных зон (zonal flows) породит постоянные давления-ловушки → планеты/ядра могут остановиться или мигрировать медленнее.
- MHD-ветры отводят угловой момент, по сути уменьшают плотность газа и могут приводить к внутреннему «перекрытию» миграции (меньше газа → слабее торки).
- Нон-идеальная MHD (амбиполярное рассеяние, Холл-эффект) меняют структуру турбулентности и поэтому эффективность коротационных торков; в некоторых условиях магнитные силы могут вызвать нерегулярные ускоренные дрейфы (случайные «киксы»).
- Итог: сильное излучение и сильные магнитные эффекты обычно замедляют или локально останавливают миграцию (через уменьшение Σ\SigmaΣ, создание ловушек давления или уменьшение торков); в активных турбулентных зонах миграция может быть stochastic.
Наблюдения ALMA, подтверждающие модель (что именно смотреть)
1. Пыльные кольца/щели и асимметрии (continuum high‑res, λ∼0.8–3 mm):
- Наличие узких колец/клиней и ярких «лис» (crescents) → давление-ловушки/вихри.
- Многодорожечные карты (разные длины волны) дают карту спектрального индекса αν\alpha_\nuαν — снижение αν\alpha_\nuαν указывает на рост зерен (большие частицы, St↑).
2. Радиационно/газовое распределение (CO, C18O, HCO+):
- Сопоставление газа и пыли: зоны с дефицитом газа при наличии плотной пыли указывают на фильтрацию и gap-opening.
- Кинематические особенности: «kinks» и локальные отклонения от Кеплеровского потока → наличие планет и невязкости/локальной массы.
3. Турбулентность и α:
- Линии ширины (CO isotopologues, CS) для оценки turbulent linewidth → ограничение α\alphaα.
4. Ионизация и магнитные поля:
- Линии-ионы (HCO+, N2H+, CN) как трейсеры ионизации; дефицит/номенклатура указывает на мёртвую зону.
- Поляризация континуума и линий (dust polarization, aligned grains, Goldreich-Kylafis) для карты магнитного поля (морфология B‑поля).
5. Доказательства ветров/выбросов:
- Высоко-скоростные крылья в CO → диск‑ветер (MHD wind / photoevaporative flow).
6. Пространственно‑разрешённые карты размера зерен:
- Радіальные профили размера зерен (multi-λ) показывают концентрацию пепла в трапах; сопоставление с моделями St и τ_drift.
7. Индикаторы стриминговой неустойчивости / перегрузки:
- Очень высокая локальная dust-to-gas (через continuum / C18O) → условие для SI (ZlocalZ_{\rm local}Zlocal превышает ZcritZ_{\rm crit}Zcrit ).
Коротко: ожидаемые ALMA-признаки модели — совпадение пылевых колец/асимметрий с локальными изменениями газа и кинематикой, снижение спектрального индекса в ловушках, малые скорости турбулентности в «мертвых зонах», следы ветров и поляризация, соответствующая конфигурация магнитного поля. Эти данные позволяют связать: (i) места концентрации пепла → планетезимали, (ii) уменьшение газа/градиенты → влияние радиации, (iii) морфологию B‑поля и уровни турбулентности → влияние магнитных эффектов на миграцию.