Моделюва́ння нуклеосинтезу — розділ ядерної астрофізики, що займається чисельним описом утворення хімічних елементів у зорях, наднових, кілоновах та інших астрофізичних об'єктах. Воно ґрунтується на інтегруванні системи диференціальних рівнянь для реакційних мереж, які описують еволюцію ізотопного складу залежно від температури, густини, електронної частки та потоків частинок.
Моделювання нуклеосинтезу поєднує ядерні дані (січення реакцій, періоди напіврозпаду, енергії зв'язку), астрофізичні моделі (структуру зір, гідродинаміку вибухів, акреційні диски) та чисельні алгоритми розв'язання жорстких систем рівнянь. Для цього використовують спеціалізовані програмні пакети, зокрема MESA, NuGrid, SkyNet, PRISM, XNet, Torch, WinNet, Webnucleo, pynucastro та інші.
Предмет і завдання
Основне завдання моделювання нуклеосинтезу — обчислити зміни хімічного та ізотопного складу речовини у заданому астрофізичному середовищі. Це дає змогу:
пояснити походження хімічних елементів і їх ізотопів у Всесвіті; інтерпретувати спостережні дані: спектри зір, наднових, кілонових, метеоритів; знайти чутливість до невідомих ядерних параметрів (січень, мас, періодів напіврозпаду) і обґрунтувати необхідні експерименти; тестувати моделі будови та еволюції зір, нейтронних зір, акреційних дисків тощо.
На практиці моделювання нуклеосинтезу майже завжди зводиться до розв'язання системи рівнянь для молярних часток (або масових часток) ізотопів: , де — молярна частка ядра з номером , а — швидкості ядерних реакцій чи розпадів. Кількість рівнянь може сягати тисяч і більше.
Типи нуклеосинтезу в контексті моделювання
Різні фізичні процеси нуклеосинтезу вимагають різних реакційних мереж і підходів до моделювання:
Первинний нуклеосинтез — Big Bang nucleosynthesis (високі температури, але прості мережі для легких елементів). Гідростатичний нуклеосинтез у зорях — протон-протонний ланцюжок, CNO-цикл, потрійна альфа-реакція, вуглецеве, кисневе, неонове, кремнієве горіння в ядрах та оболонках зір головної послідовності, гігантів та надгігантів. Тут важлива еволюція структури зорі. s-процес (повільний захват нейтронів) — характерний для AGB-зір та певних стадій еволюції масивних зір; моделюється великими мережами, але при відносно помірних потоках нейтронів. r-процес (швидкий захват нейтронів) — відбувається, ймовірно, в нейтронно-зоряних злиттях, кілоновах, можливо в гіпернових та інших екстремальних середовищах; потребує тисяч нестабільних ізотопів і складних моделей. rp-процес (rapid proton capture) та пов'язані з ним процеси p-нуклеосинтезу — важливі для моделювання рентгенівських спалахів на поверхні нейтронних зір, термонуклеарних вибухів на білих карликах тощо. Вибуховий нуклеосинтез — в оболонках наднових (Ia та core-collapse), у навколодискових вітрах кілонов та акреційних дисків.
У кожному випадку моделі нуклеосинтезу мають бути узгоджені з гідродинамічною моделлю джерела (структурою зорі, профілями температури та густини в часі, нейтринним полем, геометрією викидів тощо).
Чисельні методи та реакційні мережі
Реакційні мережі
Реакційна мережа — це набір ізотопів і реакцій між ними (захвати частинок, фотодисоціація, бета-розпади, поділ тощо), що визначає праву частину рівнянь для . Обираючи мережу, моделювальник має балансувати між:
повнотою (достатня кількість ізотопів і реакцій для опису потрібного процесу); чисельною складністю (розмірність системи та час інтегрування).
Для простого водневого горіння достатньо десятків реакцій, для r-процесу в нейтронно-зоряних злиттях застосовують мережі з кількома тисячами ядер та сотнями тисяч реакцій.
Жорсткі системи та інтегратори
Системи рівнянь реакційних мереж є жорсткими: існують дуже різні часові шкали (від мікросекунд для деяких реакцій до років і більше для розпадів). Тому використовують неявні або напівнеявні методи інтегрування:
методи родини backward Euler / BDF; неявні методи Рунге — Кутти; спеціалізовані алгоритми з адаптивним часом, побудовані під нуклеосинтетичні мережі.
Більшість сучасних пакетів (MESA, XNet, SkyNet, Torch, WinNet) реалізують власні варіанти неявних інтеграторів, оптимізованих під конкретні розміри мереж і архітектури обчислювальних кластерів.
Вбудований нуклеосинтез та постпроцесінг
Є два основні підходи до моделювання:
Вбудовані реакційні мережі (inline burning): нуклеосинтез інтегрують одночасно з рівняннями гідродинаміки чи зоряної структури (як у MESA). Це потрібно, коли зворотний вплив вивільнення енергії на динаміку істотний. Постпроцесінг (post-processing): спочатку розв'язують гідродинамічну/зоряну задачу зі спрощеною фізикою або малою мережею, зберігають термодинамічні траєкторії, а потім проганяють їх крізь велику реакційну мережу (NuGrid, XNet, PRISM, SkyNet). Це дозволяє використовувати дуже великі мережі, не обмежуючи гідродинамічний розрахунок.
Джерела ядерних даних
Для моделювання нуклеосинтезу необхідні:
таблиці січень ядерних реакцій і збуджень; термоядерні реакційні ставки як функції температури; маси ядер, періоди напіврозпаду, гілки розпаду, параметри поділу; рівняння стану густої матерії для внутрішніх областей зір і нейтронних зір.
Найпоширеніші бібліотеки реакційних ставок:
JINA REACLIB — велика бібліотека аналітично апроксимованих реакційних ставок; NACRE, CF88, BRUSLIB, NON-SMOKER та інші набори, що часто комбінуються в «гібридні» бібліотеки; локальні бібліотеки, побудовані в рамках конкретних проєктів.
Для автоматизованої побудови та трансформації таких бібліотек існують інфраструктурні платформи, наприклад Computational Infrastructure for Nuclear Astrophysics (CINA), яка дозволяє через вебінтерфейс формувати власні бібліотеки, запускати постпроцесінг з XNet тощо.
Програмне забезпечення для моделювання нуклеосинтезу
У ядерній астрофізиці сформувався багатий ландшафт відкритих та напіввідкритих кодів. Нижче подано основні з них.
MESA
MESA (Modules for Experiments in Stellar Astrophysics) — один з найпоширеніших пакетів для одномірного моделювання еволюції зір з урахуванням нуклеосинтезу.
Основні риси:
модульна бібліотека на Fortran 95: рівняння стану, непрозорості, ядерні реакції, атмосфера, транспорт енергії; модуль MESA star розв'язує повністю зв'язану систему рівнянь структури та складу зорі з адаптивною сіткою і кроком по часу; вбудована реакційна мережа, яку можна розширювати: від кількох десятків до сотень/тисяч ізотопів; різні джерела реакційних ставок: власна бібліотека, комбінована з CF88, NACRE, REACLIB та ін.; підтримка паралелізму (OpenMP) і масштабованість до великих кластерів.
MESA часто використовують у зв'язці з постпроцесінговими кодами (напр., NuGrid) для детального аналізу нуклеосинтезу в AGB-зір, класичних нових, масивних зорях, зорях Вольфа — Рає тощо.
NuGrid (PPN, MPPNP, інші коди)
NuGrid — колаборація, що розробляє набір інструментів для постпроцесінгового нуклеосинтезу, орієнтованого на зоряні моделі.
Основні коди:
PPN (Post-Processing Network) — однозонний постпроцесінговий код для заданих термодинамічних траєкторій. MPPNP — багатозонний (multi-zone) код, здатний обробляти повний вивід одномірних еволюційних кодів (MESA, інші), з динамічним збільшенням/зменшенням розміру мережі залежно від складу оболонок. набір Python-скриптів і Jupyter-ноутбуків (NuGridPy, NuPyCEE) для аналізу результатів та моделювання галактичної хімічної еволюції.
Особливості:
мережі з тисячами ізотопів (s-, r-, p-процеси); реакційна бібліотека, що комбінує NACRE, CF88, BRUSLIB, ранні версії REACLIB та інші джерела; використання MPI для паралельного опрацювання великої кількості зон.
SkyNet
SkyNet — модульна бібліотека ядерної реакційної мережі загального призначення, орієнтована на екстремальні сценарії, зокрема r-процес у злиттях нейтронних зір та кілонових.
Характерні риси:
реалізована на C++, з можливістю виклику з інших кодів та через обгортки; підтримує великі мережі (до тисяч ізотопів), у тому числі зосереджені на r-процесі; включає розв'язувач для ядерної статистичної рівноваги (NSE) та самонагрівну мережу (self-heating), де температура еволюціонує разом зі складом; модульна структура для підключення різних рівнянь стану, реакційних бібліотек, моделей електронного екранування; велику увагу приділено нейтрино-індукованому нуклеосинтезу (ν-процес).
SkyNet широко використовується у моделях кілонових, злиттів нейтронних зір та акреційних дисків навколо компактних об'єктів.
PRISM
PRISM (Portable Routines for Integrated nucleoSynthesis Modeling) — сучасний код реакційних мереж, створений у Лос-Аламоській національній лабораторії для детального моделювання r-процесу та пов'язаних явищ.
Особливості:
орієнтація на важкі та надважкі елементи, включно з поділом; використання теоретичних реакційних ставок, обчислених, зокрема, кодом CoH (Hauser–Feshbach); реалізація nucleosynthesis tracing — методики відстеження траєкторій ядер у мережі, що дозволяє кількісно оцінити внесок окремих реакцій і каналів розпаду у формуванні кінцевого розподілу r-процесу; інтеграція з інфраструктурою FIRE (Fission In R-process Elements) для дослідження впливу ядерних даних на r-процес.
XNet, Torch, WinNet, Webnucleo та інші реакційні мережі
У ядерній астрофізиці історично розроблено низку реакційних мереж, які й досі активно використовуються:
XNet — повністю неявний постпроцесінговий код, здатний відстежувати до ≈5000 ізотопів; широко використовується в симуляціях core-collapse наднових (у зв'язці з гідрокодами FLASH, CHIMERA). Torch — більш «старий» код загального призначення для термоядерних мереж із типовою мережею ≈500 ізотопів і можливістю розширення. WinNet — мережа, орієнтована на r-процес, здатна відстежувати понад 5800 ізотопів; застосовується для моделювання r-процесу в наднових і злиттях нейтронних зірReichert M. et al. The Nuclear Reaction Network WinNet. // Astrophysical Journal Supplement Series. — 2023. — Vol. 266, 21.. Webnucleo / NucNet Tools — набір відкритих бібліотек і вебінтерфейсів для побудови реакційних мереж, обробки ядерних даних та виконання нуклеосинтетичних розрахунків у середовищі Jupyter або на локальній машині.
До цієї ж групи належать спеціалізовані коди на кшталт r-Java 2.0 для моделювання r-процесу з графічним інтерфейсом.
Pynucastro та бібліотеки для конструкторів мереж
pynucastro — відкрита Python-бібліотека, що полегшує створення та аналіз нуклеосинтетичних мереж.
Можливості:
завантаження реакцій з REACLIB та інших бібліотек; автоматична побудова графу мережі та перевірка її придатності для конкретного завдання; генерація коду реакційної мережі на Fortran або C++ для включення у великі гідродинамічні коди; візуалізація мереж (графи реакцій, потоки, чутливості).
Такі бібліотеки відіграють роль «конструкторів мереж», дозволяючи дослідникам швидко збирати мережі для конкретних сценаріїв (наприклад, rp-процесу в рентгенівських спалахах або s-процесу в AGB-зірках).
Узагальнювальна таблиця пакетів
Нижче наведено спрощений огляд найуживаніших пакетів: