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

Исправление — 2026-09-09 — В таблице теперь указано около 64,05 h для основного состояния 90Y вместо 3,19 h — значения для метастабильного состояния 90mY. Исправлена идентификация ядерного состояния; это не доказывает, что в моделировании исследования использовалось ошибочно напечатанное значение. Уточнены модельное распределение вкладов активных зон и бассейнов, экспериментальные ограничения и предположения о неопределённостях.

Исправление и перевод исправления выполнены с помощью ИИ с целевыми проверками источников и согласованности; авторство исходной статьи не изменено.

ENSDF: 90Y / 90mY

Редакционная методика News — это исправление

Остановленный ядерный реактор не умолкает. Продукты деления, накопившиеся в его активной зоне и бассейнах, продолжают распадаться в течение месяцев или лет, а их цепочки бета-распадов испускают антинейтрино. Коллаборация 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β​ от 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 записывается как

где $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​.

Схема без соблюдения масштаба. N_obs — наблюдаемое число событий, N_pred — предсказание. B1/B2 — активные зоны, P1/P2 — бассейны; распределение 56 %/44 % предсказано моделью.

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

Измерение стало возможным благодаря редкому обстоятельству и геометрии. Обстоятельство: в 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 событиям. Эти две величины нельзя смешивать.

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

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

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

• Решающий спектральный тест. При большей статистике наблюдение в течение нескольких месяцев за относительной эволюцией $^{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 %, но исследование не определяет конечную точность модели в таком случае.

Источники

Источники изучены 9 сентября 2026 года: окончательная публикация и препринт v1 от 6 октября 2025 года.

• 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 — эталонная кинематика и сечение обратного бета-распада.

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

Приведённые числа событий и предсказания модели сопоставлены с первичной публикацией. Исходные данные не опубликованы в открытом доступе; их анализ и моделирование здесь не воспроизводились. Наличие независимого воспроизведения этим аудитом не установлено. Научную и языковую проверку выполняет ИИ, без независимой проверки человеком.

Примечание к расчёту: Δ = 106 − 88 = 18 событий. При предположении о некоррелированных приведённых ошибках σΔ = √(18² + 7²) ≈ 19,3 события. Эта приближённая проверка согласия отличается от указанной в статье значимости сигнала 5,9σ. Активность дочернего нуклида приблизительно следует за активностью долгоживущего родительского нуклида лишь после переходного периода, при допущении равновесия в цепочке распада.

Physical Review Letters — PDF

Редакционная методика News — это исправление

Для этого исправления числа событий и предсказания сопоставлены с первичной статьёй, а период полураспада — с оценёнными данными ENSDF. Простой расчёт неопределённости повторён с явным указанием допущения. Английская версия и семь переводов проверены на численную, научную и языковую согласованность; координаты схемы и отображение формул проверены отдельно. Это редакционные проверки с помощью ИИ, а не независимая экспертиза человеком или повторение эксперимента. Исходное авторство сохранено. Датированное примечание описывает существенные изменения; предыдущая редакция сохраняется для прослеживаемости. Неразрешённые утверждения остаются с оговорками.