Аккреционный диск — где гравитация превращается в свет | Интерстеллар (Interstellar)

Астрофизика компактных объектов. Аккреционный диск

Вращающийся газовый поток, закручивающийся внутрь, теряет угловой момент из-за вязкости и преобразует гравитационную энергию в излучение — это двигатель каждого квазара, рентгеновской двойной системы и протозвезды. Аккреционный диск представляет собой уплощенную, вращающуюся структуру, в которой газ медленно закручивается на центральное тело. Вязкие напряжения переносят угловой момент наружу, высвобождая гравитационную потенциальную энергию в виде тепла и света — с эффективностью до 42% от массы покоя для вращающихся черных дыр, что в десятки раз лучше, чем ядерный синтез.

Стандартная модель: Шакура и Сюняев, 1973

– α-вязкость: 0,01–0,3

– Внутренний стабильный круговой орбит (ISCO) Шварцшильда: 6 GM/c²

– Максимальная эффективность Керра: 42% от mc²

– Масштабирование температуры: T ∝ r⁻³ᐟ⁴

Почему образуются диски, а не сферы

Сферическое свободное падение на компактный объект встречается редко в реальной Вселенной. Почти любой участок газа, падающий к звезде, нейтронной звезде или черной дыре, прибывает с как минимум небольшим угловым моментом — полученным от двойного компаньона, галактического потенциала или просто турбулентности в молекулярном облаке. Даже небольшое вращение имеет последствия. Когда газ падает внутрь, сохранение углового момента требует, чтобы его скорость вращения увеличивалась как обратная величина расстояния. На некотором радиусе центробежный барьер совпадает с гравитацией, и газ не может падать дальше радиально; вместо этого он оседает на вращательно поддерживаемый диск в экваториальной плоскости. Движение вне плоскости быстро гасится столкновениями, в то время как вращение сохраняется. Получающийся диск представляет собой своего рода “клиринговый дом”. Материя не может просто закручиваться внутрь самостоятельно — она должна передать свой угловой момент чему-то. В реальных дисках этим “чему-то” является турбулентный стресс, который эффективно действует как вязкость: сдвиг между соседними кольцами переносит угловой момент наружу. Когда участок теряет L, он опускается на меньшую орбиту, отдает гравитационную потенциальную энергию и излучает избыток. Это фундаментальный цикл: вязкость перемещает угловой момент наружу, масса дрейфует внутрь, связанная энергия выходит в виде света.

Модель α-диска Шакуры-Сюняева

Базовая теоретическая основа — статья 1973 года Николая Шакуры и Рашида Сюняева. Они записали уравнения для стационарного, геометрически тонкого, оптически толстого диска, в котором неизвестная вязкость ν параметризуется как ν = α c_s H, где c_s — местная скорость звука, H — вертикальный масштаб высоты, а α — безразмерная величина, предполагаемая порядка единицы или меньше. Эмпирически, значения α, соответствующие наблюдениям, группируются между 0,01 и 0,3. Модель дает чистые прогнозы для радиальных профилей поверхностной плотности Σ, температуры в средней плоскости T и эффективной температуры T_eff. Ключевой результат для наблюдателя заключается в том, что T_eff(r) = [3GM Ṁ / (8π σ r³)]^(1/4) × [1 − (r_in/r)^(1/2)]^(1/4) ∝ r^(−3/4) (далеко от внутреннего края). Каждое кольцо излучает примерно как черное тело при своей местной T_eff. Общий спектр представляет собой сумму над кольцами — “многоцветный диск черного тела”, который заканчивается возле температуры самого горячего, внутреннего кольца. Эта одна кривая, с ее характерным плоским вершинным F_ν ∝ ν^(1/3) средним участком и высокочастотным срезом Вина, является диагностическим отпечатком тонкого диска.

Что на самом деле движет вязкостью

α-предписание было намеренно сделано “черным ящиком”. Молекулярная вязкость горячей ионизированной плазмы на несколько порядков меньше, чем необходимо для объяснения наблюдаемых скоростей аккреции: она предсказывала бы α ~ 10^(-12), тогда как выводимые значения составляют 10^(-2) до 10^(-1). В течение двух десятилетий источник недостающего стресса был известной нерешенной проблемой. Прорыв произошел благодаря Стивену Балбусу и Джону Хоули в 1991 году, которые поняли, что магниторотационная нестабильность (MRI) — известная физикам плазмы с 1959 года (Велихов) и 1960 года (Чандрасекар), но никогда не связанная с дисками — является общей в магнитизированной, дифференциально вращающейся жидкости. Слабое магнитное поле дестабилизируется сдвигом, вызывает турбулентность, которая извлекает угловой момент из внутренних колец и откладывает его на внешние, и самоподдерживается через действие динамо. Турбулентность MRI воспроизводит α ~ 0,01 до 0,1 в симуляциях коробки сдвига и сейчас является канонической картиной для ионизированных дисков. Холодные, частично нейтральные протопланетарные диски имеют “мертвые зоны”, где MRI подавляется; в этих регионах другие механизмы переноса (гравитационные нестабильности, магнитные дисковые ветры) берут верх.

Сколько энергии высвобождается

Для частицы, движущейся от бесконечности к круговой орбите на радиусе r вокруг невращающейся массы M, высвобождаемая энергия связи составляет E_bind = (1/2) GM/r (на единицу массы, ньютоновская). Половина гравитационной потенциальной энергии превращается в орбитальную кинетическую энергию; другая половина излучается, когда частица движется внутрь. К моменту достижения ISCO на 6 GM/c² для черной дыры Шварцшильда полностью релятивистский расчет дает общую эффективность η = 1 − √(8/9) ≈ 0,057 от массы покоя — примерно 6%. Если считать энергию покоя материи, в конечном итоге падающей через горизонт, входной, то доля, выходящая в виде света, составляет примерно 5,7%, обычно округляемая до “около 10%”. Для максимально вращающейся черной дыры Керра (a = M, проградная орбита) ISCO отступает до 1,235 GM/c², а эффективность повышается до 0,423 — 42% от mc². Это более чем в 50 раз превышает эффективность водородного синтеза (0,7%) и объясняет, почему активные галактические ядра могут затмить свои родительские галактики из области размером с Солнечную систему.

Температуры диска и масштабирование массы

Для тонкого диска, излучающего на пределе Эддингтона, пиковая T_eff масштабируется как M^(-1/4). Подставив числа в формулу Шакуры-Сюняева при r = 10 GM/c²:

|—|—|—|—|—|

Это одно масштабирование объясняет огромное количество наблюдаемой феноменологии. УФ-спектр активного галактического ядра, мягкий избыток рентгеновских двоек, ИК-избыток звезд T Tauri — все это разные кольца одного и того же многоцветного диска, наблюдаемого в огромном диапазоне масс.

Пример: сколько материи питает квазар?

Рассмотрим сверхмассивную черную дыру массой 10⁸ M☉, излучающую на пределе Эддингтона. Лимит Эддингтона составляет L_Edd = 1,26 × 10³⁸ (M/M☉) эрг/с = 1,26 × 10⁴⁶ эрг/с для M = 10⁸ M☉. Принимаем радиационную эффективность η = 0,1 (невращающаяся дыра, консервативно). Связь между скоростью аккреции и светимостью: L = η Ṁ c² → Ṁ = L / (η c²). Подставив L = 1,26 × 10⁴⁶ эрг/с, η = 0,1, c² = 9 × 10²⁰ см²/с²: Ṁ = 1,26 × 10⁴⁶ / (0,1 × 9 × 10²⁰) = 1,4 × 10²⁵ г/с ≈ 1,4 × 10²⁵ × (3,156 × 10⁷ с/год) г/год ≈ 4,4 × 10³² г/год ≈ 2,2 M☉/год (M☉ = 2 × 10³³ г). Таким образом, яркий квазар потребляет около двух солнечных масс газа каждый год. За типичный цикл активности 10⁷ лет это накапливается до примерно 2 × 10⁷ M☉ — сопоставимо с общей массой сверхмассивной черной дыры. Это размерное рассуждение стоит за соотношением Солтана, связывающим плотность светимости квазаров с плотностью массы сверхмассивных черных дыр: интегрированное излучение активных галактических ядер показывает, сколько массы современных черных дыр образовалось за счет яркой аккреции. Если вместо этого предположить η = 0,42 (максимально вращающаяся Керр), то необходимая Ṁ составляет только 0,5 M☉/год — это иллюстрирует, почему вращение сверхмассивных черных дыр за счет выровненной аккреции является важным наблюдаемым фактором, влияющим на оценки скоростей питания.

Режимы: тонкие, толстые, тонкие, ADAF

Решение Шакуры-Сюняева применимо в диапазоне скоростей аккреции примерно 0,01 ≲ Ṁ/Ṁ_Edd ≲ 0,3. За пределами этого диапазона берут верх другие режимы, каждый со своими наблюдаемыми сигнатурами.

|—|—|—|—|—|

Два крайних случая — тонкие диски и ADAF — возникают, когда радиационное охлаждение не может идти в ногу с вязким нагревом. В тонком диске фотоны захвачены в потоке, потому что диффузия медленнее, чем адвекция; диск вздуется, и светимость насыщается возле нескольких L_Edd. В ADAF газ настолько разрежен, что он безударен на соответствующих временных масштабах; ионы не могут передать свою тепловую энергию электронам, газ нагревается до ~10⁹ K, и большая часть энергии адвектирована через горизонт, а не излучена. Sgr A* и центры почти всех массивных эллиптических галактик являются ADAF.

Варианты и расширения

– Диск Новикова-Торна. Общерелятивистская версия Шакуры-Сюняева, введенная в 1973 году Игорем Новиковым и Кипом Торном. Добавляет поправки общей теории относительности к гравитации, красному смещению и внутреннему граничному условию; воспроизводит эффективности ISCO Шварцшильда и Керра.

– Тонкий диск (Абрамович и др., 1988). Отменяет приближение "тонкого". При скоростях аккреции Ṁ близких или выше L_Edd диск становится геометрически толстым, захват фотонов делает охлаждение неэффективным, и светимость насыщается. Модель, питающая современную картину сверхлюминесцентных рентгеновских источников.

– ADAF / RIAF (Нараян-И, 1994; Юань-Нараян). Адвекционно-доминирующий поток аккреции / радиационно-неэффективный поток аккреции. При очень низких скоростях аккреции ионы и электроны декуплируются, и большая часть тепла адвектирована через горизонт. Основание для моделирования Sgr A и M87.

– Усеченный диск + горячая корона. Составная геометрия, используемая для рентгеновских двоек в их твердом состоянии: тонкий диск простирается от внешнего края до ~10–100 r_g, где он переходит в горячую, оптически тонкую Комптоновскую корону, производящую хвост закона мощности.

– Магнитно-арестованный диск (MAD). Когда магнитный поток, накапливающийся возле дыры, становится достаточно сильным, чтобы остановить аккрецию устойчивым образом, диск переходит в режим MAD. Наблюдается в поляриметрии EHT M87* и может коррелировать с запуском релятивистских джетов.

Где наблюдаются аккреционные диски

– Активные галактические ядра. Сверхмассивные черные дыры массой 10⁶–10¹⁰ M☉, аккрецирующие при 10⁻²–1 L_Edd. Температуры диска 10⁴–10⁵ K производят УФ-“Большой синий бугор”; Комптоновское излучение короны дает рентгеновский закон мощности 2–100 кэВ. Квазары (самые яркие АГЯ) достигают L ~ 10⁴⁷ эрг/с.

– Рентгеновские двойные системы. Звездные черные дыры или нейтронные звезды массой 1,4–30 M☉, аккрецирующие от звездного компаньона. Диски пикуются при ~1 кэВ, дают “мягкое” состояние спектра; в “твердом” состоянии внутренний диск заменяется горячей Комптоновской короной. Cygnus X-1 и GX 339-4 являются каноническими примерами.

– Катаклизмические переменные. Белые карлики (~1 M☉), аккрецирующие от низкомассивных компаньонов через переполнение лобовой полости Роша. Температуры диска 10⁴–10⁵ K, пик в УФ. Карликовые новые, такие как SS Cygni, проходят через термально-вязкие вспышки каждые ~50 дней, отличную местную лабораторию для физики дисков.

– Протопланетарные диски. 0,001–0,1 M☉ газа и пыли вокруг формирующейся звезды T Tauri или Herbig Ae/Be. Холодные (~10–1500 K), частично нейтральные. Пик в ИК и субмиллиметровом диапазоне; ALMA разрешает зазоры и кольца (HL Tau, TW Hya), интерпретируемые как формирующиеся планеты, очищающие свои зоны питания.

– События разрыва приливов. Когда звезда заходит внутрь приливного радиуса сверхмассивной черной дыры, она разрывается, и примерно половина ее массы возвращается на высокоэксцентрические орбиты, чтобы закруглиться в транзиторный аккреционный диск. Результатом является вспышка рентгеновского/УФ излучения в течение месяцев-лет, например ASASSN-14li, PS1-10jh.

Временная изменчивость — диск не стационарен

Реальные диски наблюдаются изменчивыми на всех временных масштабах от миллисекунд (высокочастотные QPO рентгеновских двоек) до десятилетий (долгосрочные события изменения вида АГЯ). Релевантные временные масштабы диска на радиусе r:

t_dyn = 1 / Ω_K ≈ 10⁻⁴ с × (r/r_g)^(3/2) для 10 M☉

t_th = t_dyn / α тепловой временной масштаб

t_visc = t_dyn / (α (H/r)²) вязкий (радиальный дрейф)

Для 10 M☉ черной дыры звездной массы при r = 10 r_g и α = 0,1: t_dyn ~ 3 мс, t_th ~ 30 мс, t_visc ~ несколько секунд (при условии H/r ~ 0,1). Для 10⁸ M☉ сверхмассивной черной дыры те же орбиты масштабируются до часов и веков соответственно. Наблюдаемая изменчивость УФ-излучения АГЯ на временных масштабах лет-декад значительно быстрее канонического временного масштаба тонкого диска t_visc, открытая загадка, которая мотивировала пересмотренные модели дисков (например, магнитно-индуцированный поток).

Распространенные ошибки

– Спутать ISCO с горизонтом событий. ISCO — это внутренний край диска; горизонт событий (радиус Шварцшильда 2 GM/c²) намного меньше. Фотоны, испущенные внутри ISCO, все еще могут вырваться; материя внутри него падает почти свободно.

– Рассматривать α как константу природы. α — это параметризация, а не фундаментальная величина. Она варьируется с состоянием диска, магнитной топологией, степенью ионизации. Симуляции МГД показывают, что “α” повышается с 0,01 в стратифицированных коробках MRI до 0,1+ в конфигурациях MAD.

– Забыть поправку граничного условия внутреннего края. Радиальный профиль T_eff включает фактор [1 − (r_in/r)^(1/2)]^(1/4), который исчезает на внутреннем крае — без него T кажется расходящимся при r → r_in. Фактор отражает тот факт, что на внутренней границе не действует крутящий момент.

– Неправильно применять L = η Ṁ c² выше Эддингтона. Эффективность тонкого диска η — это свойство геометрии. Когда Ṁ превышает Ṁ_Edd, захват фотонов уменьшает эффективную эффективность, и берет верх режим тонкого диска; использование η тонкого диска значительно завышает светимость.

– Игнорировать корону. Оптически тонкая горячая корона отвечает за 10–80% рентгеновской светимости в реальных системах и не входит в стандартное решение диска. Моделирование широкополосных рентгеновских спектров без нее приводит к выводам о параметрах диска, противоречащим наблюдениям.

Часто задаваемые вопросы

Почему газу нужно терять угловой момент, чтобы упасть внутрь?

Частица в круговой кеплеровской орбите на радиусе r имеет специфический угловой момент L = √(GMr). Чтобы опуститься на меньший радиус, она должна потерять этот L — и угловой момент сохраняется, поэтому L должен куда-то деться. В диске вязкие напряжения переносят его наружу к внешнему диску, который поэтому постепенно расширяется, в то время как внутренний диск истощается. Без этого переноса газ просто вечно бы вращался и не аккрецировался.

Что такое α-вязкость, и почему это просто параметр?

Молекулярная вязкость в астрофизических дисках на несколько порядков меньше, чем необходимо для объяснения наблюдаемых скоростей аккреции. Шакура и Сюняев параметризовали неизвестный источник стресса как ν = α c_s H, где c_s — местная скорость звука, H — масштаб высоты, а α обычно составляет от 0,01 до 0,3. Микрофизика теперь, как считается, является магниторотационной нестабильностью (MRI), открытой для аккреционных дисков Балбусом и Хоули в 1991 году, которая генерирует магнитогидродинамическую турбулентность, действующую как эффективная вязкость.

Почему АГЯ пикуются в УФ, а рентгеновские двойные системы в рентгене?

Потому что пиковая температура диска масштабируется как T ∝ M^(-1/4) для фиксированной дроби Эддингтона. Диск черной дыры звездной массы 10 M☉ пикуется возле 10⁷ K — мягкий рентген. Диск сверхмассивной черной дыры 10⁸ M☉ пикуется возле 10⁵ K — УФ. Это объясняет, почему аккреция на галактических ядрах производит УФ-“Большой синий бугор”, а звездные рентгеновские двойные системы светятся в 0,1–10 кэВ рентгене.

Насколько эффективна аккреция по сравнению с ядерным синтезом?

Водородный синтез преобразует около 0,7% энергии покоя. Аккреция на невращающуюся черную дыру Шварцшильда — ограниченная ISCO на 6 GM/c² — высвобождает около 5,7% (или около 10%, если считать энергию покоя аккрецируемой материи входной). Максимально вращающаяся черная дыра Керра, с ISCO на 1,235 GM/c², извлекает до 42%. Аккреция поэтому является самым эффективным источником энергии во Вселенной, примерно в десять раз эффективнее синтеза.

Что такое предел Эддингтона?

Предел Эддингтона — это верхний предел, при котором давление излучения на свободные электроны (через рассеяние Томсона) точно уравновешивает гравитацию на полностью ионизированной водородной плазме. L_Edd = 4πGM m_p c / σ_T ≈ 1,26 × 10³⁸ (M/M☉) эрг/с. Выше его излучение выталкивает ветер, который нарушает устойчивую аккрецию. Реальные диски могут кратковременно превышать L_Edd, если геометрия не сферическая (тонкие диски, сверх-Эддингтоновские воронки), но долгосрочное устойчивое состояние — подсветимое.

Являются ли протопланетарные диски той же физикой, что и диски черных дыр?

Бухгалтерия переноса углового момента идентична: вязкость перемещает L наружу, масса дрейфует внутрь, гравитационная связанная энергия высвобождается в виде излучения. Различия заключаются в масштабе и микрофизике. Протопланетарные диски холодные (10–1000 K), частично нейтральные (так что MRI подавляется в “мертвых зонах”), пыльные, самогравитирующие в внешних регионах. Диски АГЯ горячие, полностью ионизированные, и доминируются давлением излучения в своих внутренних частях. Основные уравнения тонкого диска работают для обоих; границы нагрева, охлаждения и стабильности различаются.

Что такое ISCO, и почему оно задает внутренний край?

Внутренний стабильный круговой орбит (ISCO) — это наименьший радиус, на котором тестовая частица может поддерживать круговую орбиту вокруг черной дыры. Внутри ISCO возмущенные орбиты нестабильны, и частица падает. Для Шварцшильда ISCO = 6 GM/c² ≈ 9 км × M/M☉. Для максимально вращающейся Керр (a = M, проградная орбита) ISCO отступает до 1,235 GM/c². ISCO обычно рассматривается как внутренний край диска; материя там практически находится в свободном падении и испускает мало дополнительного излучения.