Double Chooz публикует первое количественное измерение остаточного потока антинейтрино

Остановленный ядерный реактор не умолкает. Продукты деления, накопившиеся в его активной зоне и бассейнах, продолжают распадаться в течение месяцев или лет, а их цепочки бета-распадов испускают антинейтрино. Коллаборация Double Chooz публикует первое количественное измерение этого остаточного потока, полученное благодаря одновременным остановкам обоих реакторов АЭС Шооз.

Источник: journals.aps.org

Простыми словами

Когда гасят дровяной костёр, угли ещё долго тлеют. С ядерным реактором происходит то же самое, но его свечение нельзя увидеть глазами: это поток частиц-призраков, антинейтрино, возникающий в цепочках бета-распадов продуктов деления в облучённом топливе, которое остаётся в активных зонах, и в отработавших сборках, хранящихся в бассейнах.

За 17,2 дня живого времени, когда оба реактора в Шооз, в Арденнах, были остановлены, ближний детектор — расположенный примерно в 400 метрах от активных зон и защищённый 115 метрами водного эквивалента — выделил 106 ± 18 остаточных кандидатов после вычитания фона. Это около шести в день. Подробная модель предсказывает, что 44 % сигнала поступает из бассейнов, где охлаждается отработавшее топливо; детектор напрямую не разделяет эти доли.

Предполагаемая практическая ценность связана с контролем ядерных материалов: антинейтрино проходят сквозь вещество с очень малой вероятностью взаимодействия, позволяя наблюдать реактор без доступа к его активной зоне. Однако этот результат получен с использованием нескольких десятков кубических метров сцинтиллятора и подземной защиты в 115 метров водного эквивалента. Он подтверждает расчёт остаточного потока, но не доказывает, что переносной прибор сможет работать на поверхности атомной электростанции.

Техническое объяснение

1. Как поймать антинейтрино и почему их требуется много. Используемая реакция — обратный бета-распад (IBD): $\bar{\nu}_e + p \rightarrow e^+ + n$νˉe​+p→e++n. Электронное антинейтрино взаимодействует со свободным протоном — ядром водорода в сцинтилляторе — и порождает нейтрон и позитрон. Сигнатура состоит из двух частей, и именно это делает измерение возможным при частоте всего в несколько событий в день. Позитрон сразу отдаёт свою кинетическую энергию, а затем аннигилирует: это быстрый сигнал. Нейтрон рассеивается до захвата; метод «полного захвата нейтронов» Double Chooz использует захваты на гадолинии, водороде и углероде во внутренних объёмах, формируя задержанный сигнал. Требование двух сигналов в окне от 0,5 до 800 µs и на расстоянии менее 1,2 m друг от друга устраняет основную часть окружающей радиоактивности. Энергия антинейтрино восстанавливается по быстрому сигналу как $E_{\bar{\nu}_e} \simeq E_{\text{prompt}} + 0{,}78\ \text{MeV}$Eνˉe​​≃Eprompt​+0,78 MeV, а кинематический порог реакции составляет около 1,8 MeV: ниже него этот канал антинейтрино не регистрирует.

2. Откуда берётся остаточное свечение: три цепочки, трое часов. Работающий реактор получает основную часть потока антинейтрино от короткоживущих продуктов деления, вклад которых обрушивается в течение нескольких часов после остановки. Регистрируемый остаток поступает главным образом от трёх пар «родитель–дочерний нуклид»: во время облучения, а затем охлаждения долгоживущий родитель непрерывно питает гораздо более короткоживущий дочерний нуклид. Поскольку последний быстро исчезает после каждого образования, его активность в основном следует активности родителя; остановка прежде всего убирает доминирующий поток, который скрывал эти устойчивые цепочки.

   Пара	Период полураспада родителя	Дочерний нуклид: период полураспада и максимальная энергия $Q_\beta$Qβ​	Роль
   $^{144}\mathrm{Ce} \rightarrow {}^{144}\mathrm{Pr}$144Ce→144Pr	285 дней	17,3 min · 3,00 MeV	≈ 54 % ожидаемого моделью сигнала 1–3 MeV
   $^{106}\mathrm{Ru} \rightarrow {}^{106}\mathrm{Rh}$106Ru→106Rh	372 дня	30,1 s · 3,54 MeV	≈ 38 % ожидаемого моделью сигнала
   $^{90}\mathrm{Sr} \rightarrow {}^{90}\mathrm{Y}$90Sr→90Y	28,9 года	3,19 h · 2,28 MeV	становится доминирующим в долгосрочном предсказании модели (> 90 % через десятилетие)

   Эти энергии определяют всё остальное: при $Q_\beta$Qβ​ от 2 до 3,5 MeV и пороге IBD 1,8 MeV регистрируется лишь высокоэнергетический хвост каждого бета-спектра. Поэтому авторы ожидают 98,7 % регистрируемого остаточного потока ниже 3 MeV — этим обосновано ограничение анализа окном 1–3 MeV по энергии быстрого сигнала, и это же делает задачу трудной: именно здесь фон наиболее плотен.

3. Что именно измеряет окно 1–3 MeV и какой контроль подтверждает результат. Если ограничить анализ окном, где ожидается сигнал, возникает классическое возражение: искомое находят там, где его ищут. Авторы отвечают отрицательным контролем. Выше 3 MeV выходы после вычитания фона совместимы с нулём. Модель действительно предсказывает небольшой остаточный хвост, но менее 0,1 события в день, то есть ниже нынешней чувствительности. Отсутствие измеримого избытка в этой области поддерживает модель фона, но не доказывает, что физический поток там строго равен нулю.

4. Фон и значение подземной защиты. При шести событиях в день важно всё. Три семейства фона рассматриваются раздельно. Случайные совпадения — два не связанных друг с другом энерговыделения, одно из которых принимается за быстрый сигнал, а другое за задержанный, — подавляются нейронной сетью с коэффициентом отбраковки более 400; остаток измеряется in situ по частоте одиночных энерговыделений. Коррелированный фон порождают космические мюоны: это либо быстрые нейтроны, возникающие при расщеплении ядер в породе, которые имитируют быстрый сигнал и затем захватываются, либо β-n-изотопы, такие как $^{9}\mathrm{Li}$9Li, образованные при мюонном расщеплении $^{12}\mathrm{C}$12C сцинтиллятора, чей распад сам даёт пару «электрон–нейтрон». Последние отмечаются по временной и пространственной корреляции с родительским мюоном. Защита ближнего детектора в 115 m водного эквивалента уменьшает мюонный поток и, следовательно, образование этих изотопов, которые невозможно полностью устранить никаким вето.

5. Предсказание: расчёт инвентаря для каждой сборки. Ожидаемое число — не оценка порядка величины, а результат подробного моделирования. Предсказанный спектр IBD записывается как
   $\frac{dN_{\bar{\nu}_e}}{dE_{\text{vis}}} = \sum_i \frac{N_p}{4\pi L_i^2}\iint \varphi_i(E,t)\,\sigma_{\text{IBD}}(E)\,P_{ee}(E,L_i)\,\mathcal{D}(E_{\text{vis}}|E)\,dE\,dt,$dEvis​dNνˉe​​​=i∑​4πLi2​Np​​∬φi​(E,t)σIBD​(E)Pee​(E,Li​)D(Evis​∣E)dEdt,
   где $N_p$Np​ — число протонов-мишеней, $L_i$Li​ — расстояние до источника $i$i, $\varphi_i$φi​ — зависящий от времени спектр излучения, $\sigma_{\text{IBD}}$σIBD​ — сечение, $P_{ee}$Pee​ — вероятность выживания с учётом осцилляций, а $\mathcal{D}$D — отклик детектора. Активности продуктов деления рассчитаны кодами APOLLO-2.8.4 и DARWIN-3 с библиотекой ядерных данных JEFF-3.1.1; бета-спектры взяты из библиотеки BESTIOLE, которая применяет теорию Ферми с поправками на электромагнитные эффекты, конечный размер ядра, атомное экранирование и форм-факторы для разрешённых и запрещённых переходов. История облучения каждой топливной сборки моделировалась отдельно, с отслеживанием её тепловой мощности в течение всего срока службы вплоть до начала каждого периода остановки.

6. Почему модель приписывает 44 % бассейнам. Два бассейна хранения расположены близко к активным зонам, поэтому для детектора примерно в 400 m их геометрический фактор $1/L^2$1/L2 остаётся сопоставимым. Каждая сборка поступает туда после нескольких циклов облучения и охлаждается годами. Остановленная активная зона содержит свежий инвентарь, который быстро распадается; бассейн накапливает сборки разного возраста, в которых некоторые изотопы сохраняются долго. При длительной остановке суммарный предсказанный вклад бассейнов становится сопоставимым с вкладом активных зон. Интегральный расчёт даёт 56 % для активных зон и 44 % для бассейнов; существующая геометрия не позволяет измерить эти доли раздельно.

7. Что ограничивает предсказание и каких выводов это не позволяет делать. Неопределённость нормировки предсказанного числа достигает 7,4 % и складывается из 6,0 % спектрального моделирования, 2,9 % геометрии расстояний до источников, 2,1 % инвентаря продуктов деления, 2,0 % инвентаря сборок в бассейнах и менее 0,8 % для детектора. Доминирующий член связан с формой бета-перехода $^{144}\mathrm{Pr}$144Pr; его неопределённость принимается равной различию между подробным расчётом ядерной структуры и приближением, называемым $\xi$ξ. Эти 7,4 % не измеряют ни точность согласия данных с моделью, ни качество вывода об инвентаре: исследование сопоставляет подсчёт с предсказанием, не оценивая в обратном направлении массу или состав топлива.

Обозначения на схеме: $N_{\rm obs}$Nobs​ — измеренное остаточное число событий, $N_{\rm pred}$Npred​ — предсказание модели, а $H_0(B)$H0​(B) — гипотеза одного лишь фона. Величина $5{,}9\,\sigma$5,9σ сопоставляет $N_{\rm obs}$Nobs​ с $H_0(B)$H0​(B), а не с $N_{\rm pred}$Npred​.

Почему это сработало

Измерение стало возможным благодаря редкому обстоятельству и геометрии. Обстоятельство: в 2017 году оба реактора Шооз одновременно останавливались четыре раза, включая период продолжительностью 20,8 дня — это необходимое условие, поскольку ближний детектор при полной мощности регистрирует около 900 событий IBD в день против 6,2 ± 1,1 остаточного события в день здесь. Геометрия: уже построенный, откалиброванный и хорошо изученный детектор, установленный в 400 m, систематические неопределённости которого были определены за годы измерений осцилляций.

Численное сравнение прямое. Предыдущий анализ данных 2011–2012 годов собрал лишь около двадцати кандидатов — слишком мало для выделения остаточного спектра. Здесь после вычитания фона 106 ± 18 остаточных кандидатов сравниваются с 88 ± 7 ожидаемыми: разность равна $\Delta = 18 \pm 19$Δ=18±19 событиям, то есть наблюдается согласие. На дальнем детекторе получено 27 ± 13 при 14 ± 1 ожидаемых, и из-за большой погрешности вывод гораздо менее строгий.

Ограничения сформулированы явно. Сигнальный анализ относится к окну 1–3 MeV. Выше 3 MeV выходы после вычитания фона совместимы с нулём, тогда как модель сохраняет небольшой остаточный хвост примерно до 4,5 MeV, менее 0,1 события в день, который эта выборка не разрешает. Неопределённость нормировки 7,4 % характеризует предсказание; эксперимент не доказывает точность вывода об инвентаре.

Разницу между доказанным и предполагаемым важно обозначить. Доказано следующее: остаточный поток отличается от фона на уровне $5{,}9\,\sigma$5,9σ, а число $106\pm18$106±18 совместимо с предсказанием $88\pm7$88±7, поскольку $\Delta=18\pm19$Δ=18±19. Не доказана эксплуатационная возможность контроля топлива. Для демонстрации потребовалось более 30 m³ сцинтиллятора под 115 m водного эквивалента, чтобы получать около шести событий в день на расстоянии 400 m. Компактному наземному прибору ещё предстоит справиться с космическим фоном и продемонстрировать чувствительность к изменениям инвентаря.

Причинная цепочка

Деление в активной зоне → накопление при облучении долгоживущих родителей ($^{144}\mathrm{Ce}$144Ce, $^{106}\mathrm{Ru}$106Ru, $^{90}\mathrm{Sr}$90Sr), питающих короткоживущие бета-распадающиеся дочерние нуклиды → остановка реакторов и быстрое исчезновение доминирующего потока → сохранение уже установившихся цепочек в активных зонах и сборках, перенесённых в бассейны → четыре одновременные остановки общей длительностью 24,4 дня в 2017 году → 17,2 дня живого времени ближнего детектора → $106\pm18$106±18 остаточных кандидатов между 1 и 3 MeV, несовместимых с одним лишь фоном на уровне $5{,}9\,\sigma$5,9σ → небольшой предсказанный хвост выше 3 MeV, менее 0,1 события в день, который не разрешён → совместимость с моделью для каждой сборки ($88\pm7$88±7, $\Delta=18\pm19$Δ=18±19) → первый количественный ориентир для проверки расчётов остаточного потока → будущие исследования разделения активных зон и бассейнов и специализированных детекторов.

Историческая деталь

Double Chooz уже пыталась провести это измерение по данным 2011–2012 годов. Тогда коллаборация собрала лишь около двадцати остаточных кандидатов: ближний детектор ещё не был доступен, а одновременная остановка длилась слишком мало для выделения спектра. В 2017 году ближний детектор уже работал, а четыре одновременные остановки, включая одну длительностью 20,8 дня, дали 17,2 дня живого времени; благодаря этим двум изменениям количественное измерение стало возможным.

Предыстория и актуальные данные

Измерение даёт первую опубликованную точку количественного сравнения расчёта методом суммирования остаточного потока с реальным подсчётом. Предсказание имеет неопределённость нормировки 7,4 %; наблюдаемое сравнение равно $\Delta=18\pm19$Δ=18±19 событиям. Эти две величины нельзя смешивать.

С прикладной стороны наиболее непосредственно полезна новая информация о предсказанном распределении: модель приписывает 44 % сигнала остановленной площадки бассейнам, а не активным зонам. Поэтому метод контроля, моделирующий только топливо, остающееся в активных зонах, сильно недооценил бы ожидаемый сигнал. Однако публикация не позволяет количественно оценить чувствительность компактного наземного детектора: в статье нет таких данных, а авторы не выполняют экстраполяцию от активных объёмов подземной установки.

Взгляд исследователя — открытые вопросы

Это возможные направления продолжения работы, а не результаты исследования.

• Решающий спектральный тест. При большей статистике наблюдение в течение нескольких месяцев за относительной эволюцией $^{144}\mathrm{Ce}$144Ce и $^{106}\mathrm{Ru}$106Ru проверило бы временную форму модели. Для наблюдения доминирования пары $^{90}\mathrm{Sr}$90Sr–$^{90}\mathrm{Y}$90Y потребуются годы: по предсказанию её вклад превышает 90 % лишь более чем через десять лет.
• Геометрически разделить бассейны и активные зоны. Значения 56 %/44 % здесь являются предсказанием модели, а не результатом измеренного разделения: на расстоянии 400 m разность в 38 m не разрешается. Детектор, размещённый значительно ближе, или два детектора под разными углами могли бы экспериментально разделить вклады вместо их расчётного распределения.
• Ограничение со стороны ядерной структуры. Доминирующая неопределённость связана не с детектором, а с формой бета-перехода $^{144}\mathrm{Pr}$144Pr. Специализированное измерение этого бета-спектра могло бы уменьшить вклад 6,0 %, но исследование не определяет конечную точность модели в таком случае.

Источники

Ссылки проверены в ходе фактчекингового аудита 16 августа 2026 года.

• Double Chooz Collaboration, «First Measurement of Neutrino Emissions from Spent Nuclear Fuel by the Double Chooz Experiment», Physical Review Letters 137, 061803, 4 августа 2026 года — рецензированная статья. DOI: 10.1103/dr26-j19g
• Препринт той же работы: arXiv:2510.04869 — изучен полный текст (аннотация, методы, таблицы систематических неопределённостей).

Справочные источники

1. P. Vogel и J. F. Beacom — Angular distribution of neutron inverse beta decay — эталонная кинематика и сечение обратного бета-распада.

Декларация уверенности

Высокая уверенность в приведённых числах событий, значимостях, энергетических окнах, геометрии и составляющих неопределённости, воспроизведённых из рецензированной статьи. Средняя уверенность в прикладной значимости: публикация упоминает ядерные гарантии, но не испытывает ни компактный наземный детектор, ни вывод об инвентаре. Доля 44 %, приписанная бассейнам, является результатом модели, а не измеренным разделением. У этого первого количественного анализа пока нет независимого воспроизведения. Проверка уникальности: ни один DOI, идентификатор arXiv или эквивалентный центральный результат не встречается в реестре A001–A018; вердикт — уникально. DOI, даты, формулы, ссылки и отсутствие дубликатов проверены.