Конец погони за стерильным нейтрино: как одно математическое упрощение 30 лет заставляло физиков искать несуществующую частицу

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

На протяжении тридцати лет в экспериментах по физике нейтрино наблюдалось явление, которое полностью удовлетворяло этому строгому критерию, но при этом вступало в прямое противоречие с общепринятой теорией строения материи. Это явление вошло в научную литературу под названием «галлиевая аномалия».

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

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

Хроника дефицита: от проекта SAGE до эксперимента BEST

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

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

Калибровочные измерения проводились тремя независимыми научными коллаборациями:

1. SAGE (Soviet-American Gallium Experiment) — советско-американский эксперимент в Баксанской нейтринной обсерватории, расположенной в подземных выработках под горой Андырчи на Северном Кавказе на глубине двух километров.
2. GALLEX (Gallium Experiment) — европейский проект в подземной лаборатории Гран-Сассо в Италии, защищенный толщей скальных пород от космических лучей.
3. BEST (Baksan Experiment on Sterile Transitions) — современная российская установка повышенной точности, развернутая в Баксанской обсерватории в 2019 году с использованием двух независимых концентрических зон детектирования.

Результаты всех трех исследовательских программ совпали: количество реально извлеченных атомов германия-71 неизменно составляло около восьмидесяти процентов от теоретически рассчитанного объема. Разница была системной, повторялась при использовании разных источников и на разных установках, а совокупная статистическая значимость расхождения со временем превысила пять сигма.

Тупик для гипотезы о стерильных частицах

Естественной физической попыткой объяснить дефицит стали нейтринные осцилляции — квантовый процесс, при котором нейтрино одного типа в полете превращается в нейтрино другого типа.

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

Чтобы объяснить столь быстрое исчезновение частиц на ультракоротких расстояниях, ученые предложили гипотезу о существовании четвертого типа частиц — стерильного нейтрино с массой около одного электронвольта. В отличие от трех известных типов, стерильное нейтрино не должно обладать зарядами слабого взаимодействия и способно проявлять себя только через квантовое смешивание с обычными нейтрино. Предполагалось, что часть электронных нейтрино на пути в один метр переходит в это невидимое состояние и перестает реагировать с галлием.

Однако в последние годы серия прецизионных экспериментов опровергла эту гипотезу:

• Эксперимент MicroBooNE в Национальной лаборатории имени Ферми в США исследовал потоки ускорительных нейтрино на коротких дистанциях и не обнаружил признаков перехода частиц в стерильные состояния в заявленном диапазоне масс.
• Детекторы STEREO, PROSPECT и DANSS, расположенные в непосредственной близости от активных зон ядерных реакторов, с высокой точностью измерили спектры электронных антинейтрино и исключили пространственные колебания потока, обязательные для гипотезы стерильного нейтрино.
• Эксперимент KATRIN в Германии, анализирующий энергетический спектр электронов при бета-распаде трития, наложил жесткие прямые ограничения на массы нейтрино и закрыл область параметров, необходимую для объяснения дефицита на галлии.

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

Принцип детального равновесия: математический обходной путь

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

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

В конце 1990-х годов американский астрофизик Джон Бакал предложил использовать для вычислений фундаментальный физический принцип — принцип детального равновесия. Он опирается на инвариантность микроскопических процессов относительно обращения времени: вероятность прямого процесса математически строго связана с вероятностью обратного процесса.

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

Главное достоинство этого подхода заключалось в том, что вероятность электронного захвата не требовалось вычислять теоретически. Период полураспада изотопа германия-71 многократно и надежно измерен в лабораторных условиях: он составляет ровно 11,465 суток с погрешностью всего в три тысячных доли дня. Опираясь на эту экспериментальную константу, физики вывели значение ядерного матричного элемента — параметра, определяющего вероятность внутренней перестройки ядра, — и перенесли его на реакцию захвата нейтрино.

Где возникла погрешность: допущение о факторизации

Чтобы связать период полураспада германия с сечением поглощения нейтрино галлием, теоретики применили стандартный математический прием — факторизацию волновых функций.

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

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

Это допущение позволяло вынести значения функций лептонов за знак интеграла, зафиксировав их в условном центре ядра, и разделить расчет на два независимых множителя:

1. Лептонный множитель, описывающий поведение электрона и нейтрино.
2. Ядерный матричный элемент, описывающий перестройку протонов и нейтронов внутри ядра и считающийся одинаковым для прямой и обратной реакций.

Именно в этой точке была допущена ошибка, исказившая итоговый результат на двадцать процентов.

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

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

Квантовое строение ядра и узловая структура плотности

Итальянские физики провели расчет сечения реакции без применения факторизации, выполнив строгое прямое интегрирование всех компонентов процесса.

Сначала исследователи рассчитали точные релятивистские волновые функции электронов в электрическом поле ядра конечного размера, решив систему квантовых уравнений Дирака — Хартри — Фока — Слэтера. Полученные пространственные распределения электрона и нейтрино были напрямую объединены под знаком интеграла с функцией плотности слабого перехода внутри ядра.

Плотность слабого перехода описывает, как именно распределена вероятность превращения нейтрона в протон по объему ядра от его центра к периферии под действием оператора ядерного перехода Гамова — Теллера.

В ранних работах предполагалось, что плотность перехода распределена в объеме ядра монотонно: она достигает максимума в центре и плавно спадает к краям, в целом повторяя распределение электрического заряда.

Однако атомные ядра со сложной оболочечной структурой устроены иначе. В ядрах галлия и германия квантовый переход Гамова — Теллера затрагивает нуклоны, находящиеся на различных квантовых оболочках и подоболочках. Основной вклад вносит переход протонов и нейтронов между базовыми состояниями, но одновременно происходят переходы между более глубокими орбиталями с идентичным моментом импульса.

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

Результаты расчетов: устранение расхождения

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

Авторы работы построили математические модели переходной плотности ядра с учетом узловых точек и проверили их на строгое соответствие экспериментальным данным. К моделям предъявлялось абсолютное граничное условие: они обязаны точно воспроизводить экспериментальный период полураспада германия-71 (11,465 суток).

Результаты прямого расчета показали следующее:

• Модели с гладкой однородной плотностью ядра без смены знака не дают необходимого снижения вероятности и сохраняют статистическое противоречие с экспериментом.
• Модели со знакопеременной плотностью перехода, учитывающие квантовую интерференцию внутриядерных состояний, автоматически снижают теоретическое сечение поглощения нейтрино ровно на двадцать процентов.
• При этом модель полностью сохраняет соответствие измеренному периоду полураспада германия-71: расчетное значение времени распада составило от 11,463 до 11,465 суток, что укладывается в пределы экспериментальной погрешности.

Для источника нейтрино на основе хрома-51 старая теория с использованием факторизации давала завышенное значение сечения около 5,53 на десять в минус сорок пятой степени квадратных сантиметров. Новый строгий расчет снизил эту величину до диапазона от 4,39 до 4,50 на десять в минус сорок пятой степени квадратных сантиметров. Это значение находится в полном согласии с усредненным экспериментальным результатом детекторов SAGE, GALLEX и BEST, составляющим 4,44 с погрешностью 0,33 на десять в минус сорок пятой степени квадратных сантиметров.

Аналогичное согласие получено и для источника на основе аргона-37: старый расчет давал значение 6,62, в то время как новый расчет показал величину от 5,10 до 5,24 на десять в минус сорок пятой степени квадратных сантиметров, что идеально соответствует экспериментально измеренному значению 5,21 с погрешностью 0,72.

Статистическое расхождение между теорией и экспериментом снизилось с критических пяти сигма до нуля, полностью устранив аномалию.

Значение открытия для фундаментальной физики

Разрешение проблемы галлиевого дефицита имеет прямое значение для всей структуры современной физики частиц.

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

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

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

Источник: arXiv