Турбулентность у черной дыры: новые модели проливают свет на ускорение частиц

В последние годы интерес к исследованию высокоэнергетических астрофизических явлений, таких как процессы, происходящие в окрестностях сверхмассивных черных дыр, стремительно растет. Наблюдения, проведенные с помощью телескопов Event Horizon Telescope и Gravity, позволили получить изображения (квази)стационарных плазменных структур, например, дисков, окружающих черные дыры, что подчеркивает важность моделирования потоков плазмы в сильных гравитационных полях. Традиционно для моделирования подобных явлений использовались методы магнитогидродинамики (МГД), которые описывают плазму как сплошную среду. Однако общепризнано, что плазма вблизи малоактивных сверхмассивных черных дыр является бесстолкновительной, то есть ее поведение определяется не столкновениями частиц, а их взаимодействием с электромагнитными полями. Это означает, что для корректного моделирования таких систем необходимо учитывать кинетические эффекты, которые не могут быть полностью учтены в рамках МГД-подхода.

Одним из ключевых процессов, определяющих динамику аккреционных дисков, является магнитовращательная неустойчивость (MRI). MRI — это механизм, который дестабилизирует вращающийся поток плазмы в присутствии магнитного поля, приводя к развитию турбулентности и эффективному переносу углового момента, необходимого для аккреции вещества на черную дыру. Для изучения MRI в бесстолкновительной плазме используются кинетические симуляции методом частиц в ячейках (PIC), которые позволяют моделировать движение отдельных частиц в самосогласованных электромагнитных полях. Однако, из-за высоких вычислительных затрат, до сих пор такие симуляции проводились в ограниченном диапазоне масштабов, что не позволяло достоверно исследовать переход от МГД- к кинетическому режиму турбулентности.

В новой работе, опубликованной на сервисе препринтов arXiv, команда ученых из Бельгии и США представила результаты самых масштабных на сегодняшний день трехмерных PIC-симуляций MRI-турбулентности в плазме, состоящей из электронов и позитронов. Благодаря использованию мощного суперкомпьютера и усовершенствованного численного кода ZELTRON, исследователи смогли достичь беспрецедентного разделения макроскопических и микроскопических масштабов в моделируемой системе. Это позволило им детально изучить динамику MRI-турбулентности, начиная от линейной стадии развития неустойчивости и заканчивая нелинейным насыщением и развитием каскада энергии вплоть до кинетических масштабов.

В симуляциях использовалась локальная модель «сдвигового ящика» (shearing-box), которая представляет собой небольшой участок аккреционного диска, вращающегося вокруг черной дыры с кеплеровской скоростью. Для подавления нефизического ускорения частиц на начальной стадии развития MRI, связанного с пересоединением магнитных силовых линий в макроскопических токовых слоях, ученые впервые ввели в PIC-симуляции MRI-турбулентности эффект радиационного охлаждения. Этот механизм эффективно «гасил» высокоэнергетические частицы в начальной фазе, позволяя изучать ускорение частиц, обусловленное исключительно турбулентностью MRI.

Результаты симуляций показали, что на больших масштабах плазма ведет себя в соответствии с предсказаниями МГД, демонстрируя экспоненциальный рост магнитного поля и формирование характерных каналов течения. Однако на кинетических масштабах спектры турбулентности становятся более крутыми, что свидетельствует о переходе к кинетическому режиму. Кроме того, было обнаружено, что значительный вклад в перенос углового момента вносит анизотропия давления плазмы, связанная с различием в движении частиц вдоль и поперек магнитного поля. Важно отметить, что это согласуется с результатами более ранних работ, основанных на гибридных моделях, которые сочетали МГД-описание для ионов с кинетическим описанием для электронов.

Одним из главных результатов исследования стало однозначное подтверждение того, что турбулентность MRI способна ускорять частицы до нетепловых энергий. Это открытие имеет важное значение для понимания механизмов генерации высокоэнергетического излучения в аккреционных дисках черных дыр. Включение радиационного охлаждения в симуляции позволило впервые четко отделить этот эффект от нефизического ускорения частиц на начальной стадии развития MRI.

Полученные данные могут быть использованы для построения более точных моделей аккреционных дисков, учитывающих кинетические эффекты. В частности, результаты симуляций могут быть использованы для калибровки МГД-моделей с анизотропным давлением, которые применяются для интерпретации наблюдений с помощью телескопов нового поколения. В будущем авторы планируют провести аналогичные симуляции для более реалистичной плазмы, состоящей из электронов и ионов, что позволит получить еще более полную картину процессов, происходящих вблизи черных дыр.