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

Зарождение концепции: капельная модель Бора-Уилера и классический порог деления

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

Именно этот, казалось бы, грубый классический подход позволил Нильсу Бору и Джону Уилеру в 1939 году создать первую строгую теорию деления ядер, которая до сих пор служит фундаментом для преподавания ядерной физики.

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

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

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

Поверхностная энергия ( $E_s$ ). Нуклоны на границе ядра имеют меньше соседей, чем нуклоны внутри. Возникает аналог поверхностного натяжения, стремящийся стянуть ядро в идеальную сферу — форму с минимальной площадью поверхности.
Кулоновская энергия ( $E_c$ ). Протоны отталкиваются друг от друга электростатически. Это дальнодействующая сила, которая стремится разорвать каплю на части.

В сферическом ядре эти энергии имеют базовые значения $E_s^{(0)}$ и $E_c^{(0)}$ . Для тяжелых ядер, таких как уран-238, кулоновское отталкивание огромно (около 1000 МэВ), но поверхностное натяжение все еще удерживает систему в равновесии. Вопрос в том, насколько устойчиво это равновесие при деформациях.

Деформация и энергетический баланс

Чтобы понять, как ядро делится, необходимо математически описать изменение его формы. Бор и Уилер разложили радиус-вектор поверхности ядра по полиномам Лежандра. Для простоты мы рассмотрим только главную моду — квадрупольную деформацию, описываемую параметром $\alpha_2$ . При $\alpha_2 > 0$ ядро вытягивается, превращаясь в эллипсоид вращения (форму сигары).

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

E_s(\alpha_2) \approx E_s^{(0)} \left(1 + \frac{2}{5}\alpha_2^2\right)

В то же время, при вытягивании протоны в среднем отдаляются друг от друга. Кулоновское отталкивание ослабевает:

E_c(\alpha_2) \approx E_c^{(0)} \left(1 - \frac{1}{5}\alpha_2^2\right)

Полное изменение потенциальной энергии ядра $\Delta E$ при малой деформации — это сумма этих двух вкладов:

\Delta E = \left(\frac{2}{5}E_s^{(0)} - \frac{1}{5}E_c^{(0)}\right) \alpha_2^2

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

Параметр делимости

Из уравнения выше видно, что абсолютный предел стабильности наступает, когда кулоновское падение энергии превышает поверхностный рост:

\frac{1}{5}E_c^{(0)} > \frac{2}{5}E_s^{(0)}

Бор и Уилер ввели безразмерный параметр делимости $x$ , который является фундаментальной мерой близости ядра к нестабильности:

x = \frac{E_c^{(0)}}{2E_s^{(0)}}

Поскольку кулоновская энергия пропорциональна $Z^2/A^{1/3}$ , а поверхностная — $A^{2/3}$ (где $Z$ — заряд, а $A$ — массовое число), параметр $x$ прямо пропорционален отношению $Z^2/A$ .

Критическое значение, при котором $x=1$ , достигается при $(Z^2/A)_{\text{crit}} \approx 50$ . В природе таких ядер нет. Для урана-238 ( $Z=92$ , $A=238$ ) отношение $Z^2/A \approx 35.5$ , что дает параметр делимости $x \approx 0.71$ .

Поскольку для урана $x < 1$ , сферическая форма стабильна по отношению к малым деформациям. Но что происходит при больших?

Классический барьер деления

При сильном вытягивании (больших значениях $\alpha_2$ ) ядро начинает приобретать форму гантели. В этот момент кулоновское отталкивание двух формирующихся «половинок» начинает доминировать над поверхностным натяжением перешейка между ними.

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

Максимум на этом графике называется седловой точкой, а высота этого максимума над энергией основного сферического состояния — барьером деления ( $E_f$ ). Для урана-238 высота этого барьера составляет около 6 МэВ.

В рамках классической механики, чтобы ядро разделилось, ему нужно сообщить извне энергию, превышающую $E_f$ . Именно это происходит при вынужденном делении: ядро урана-235 захватывает нейтрон, получает энергию связи (около 6.5 МэВ при образовании составного ядра U-236), переваливает через барьер и делится.

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

Однако в 1940 году советские физики Георгий Флёров и Константин Петржак экспериментально зафиксировали, что уран-238 делится сам по себе, находясь в полной темноте глубоко под землей, без всякого облучения нейтронами. Классическая капельная модель столкнулась с непреодолимым противоречием. Чтобы объяснить спонтанное деление, физике потребовалось отказаться от классических траекторий и обратиться к квантовому туннелированию, о котором мы поговорим на следующем шаге.

Проблема спонтанности: квантовое туннелирование и первые оценки периодов полураспада

Для классической физики ядро урана-238 абсолютно стабильно. Его барьер деления составляет около 6 МэВ, а энергия основного состояния лежит глубоко на дне потенциальной ямы. Чтобы ядро разделилось, ему нужно откуда-то получить эту энергию активации. Однако в 1940 году Георгий Флёров и Константин Петржак экспериментально зафиксировали спонтанное деление урана глубоко под землей на станции московского метро «Динамо», где влияние космических лучей было исключено. Ядро делилось само по себе, находясь в основном энергетическом состоянии.

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

От альфа-распада к коллективной динамике

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

Но перенос этой логики на спонтанное деление столкнулся с фундаментальной концептуальной трудностью. При альфа-распаде мы рассматриваем движение одной, уже сформированной частицы с фиксированной массой внутри статического потенциала ядра-остатки. При спонтанном делении движется не отдельная частица, а форма всего ядра в целом.

Туннелирует сама геометрическая конфигурация капли, эволюционируя от сферы к гантели. Это означает, что нам нужно применять аппарат квантовой механики не к координатам отдельных нуклонов, а к макроскопическому параметру деформации, такому как $\alpha_2$ .

Квазиклассическое приближение (ВКБ) и интеграл действия

Для расчета вероятности подбарьерного перехода тяжелых ядер традиционно используется квазиклассическое приближение Вентцеля — Крамерса — Бриллюэна (ВКБ). Оно отлично работает для широких и высоких барьеров, характерных для спонтанного деления.

Вероятность прохождения через барьер $P$ определяется экспоненциальным фактором:

P \approx \exp\left(-\frac{2}{\hbar} S\right)

где $S$ — это интеграл действия вдоль траектории под барьером:

S = \int_{\alpha_{in}}^{\alpha_{out}} \sqrt{2 B(\alpha) \left[ V(\alpha) - E_{gs} \right]} \, d\alpha

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

Пределы интегрирования ( $\alpha_{in}$ и $\alpha_{out}$ ): Это классические точки поворота. Точка входа находится в основном состоянии, а точка выхода — по ту сторону барьера, где потенциальная энергия $V(\alpha)$ снова сравнивается с энергией основного состояния $E_{gs}$ .
Форма барьера ( $V(\alpha) - E_{gs}$ ): Разность между локальной высотой барьера и энергией покоящегося ядра. В капельной модели эта кривая представляет собой гладкий перевернутый параболоид вблизи седловой точки.
Массовый параметр ( $B(\alpha)$ ): Это не масса ядра в килограммах. Это эффективная инерция коллективного движения.

Практический пример: представьте, что ядро — это воздушный шарик с водой, который мы пытаемся растянуть в гантель. $V(\alpha) - E_{gs}$ — это сила упругости резины, мешающая растяжению, а $B(\alpha)$ — инерция самой перетекающей внутри воды.

Проблема гидродинамической массы

Если форма потенциала $V(\alpha)$ была относительно понятна из баланса поверхностной и кулоновской энергий, то определение массового параметра $B$ стало главным вызовом.

Что сопротивляется изменению формы ядра? В первых оценках физики (например, Яков Френкель) использовали гидродинамическую аналогию: ядро рассматривалось как капля идеальной несжимаемой жидкости с безвихревым течением. В такой модели массовый параметр $B_{hydro}$ зависит от числа нуклонов $A$ и геометрии течения.

Ошибка в массовом параметре всего на 20% приводит к изменению расчетного периода полураспада в $10^3$ — $10^5$ раз. Это связано с тем, что $B$ находится под знаком квадратного корня в показателе экспоненты.

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

Чтобы перейти от вероятности туннелирования $P$ к наблюдаемой величине — периоду полураспада $T_{1/2}$ , необходимо ввести понятие «частоты попыток» (assault frequency) $n$ .

Ядро в основном состоянии совершает нулевые квантовые колебания (вибрации формы). Каждый раз, когда ядро максимально деформируется в ходе такого колебания, оно «ударяется» о барьер деления.

T_{1/2} = \frac{\ln 2}{n \cdot P}

Частота попыток $n$ оценивается из энергии нулевых колебаний $\hbar \omega_0 \approx 1$ МэВ, что дает $n \approx 10^{20}$ ударов в секунду.

Практический пример: ядро урана-238 стучится в барьер сто квинтиллионов раз каждую секунду ( $10^{20}$ Гц), и лишь спустя $10^{16}$ лет одна из попыток заканчивается успешным туннелированием.

От расчетов на ENIAC к эмпирическим систематикам: триумф и провал

В 1947 году Стэнли Франкель и Николас Метрополис, используя один из первых компьютеров ENIAC, провели масштабные расчеты барьеров деления на основе капельной модели. Опираясь на их результаты и приближение ВКБ, в 1952 году Гленн Сиборг (а также независимо от него Уильям Уайтхаус и Уильям Гэлбрейт) сформулировали красивый и логичный вывод.

Поскольку эффективная высота барьера деления в капельной модели убывает с ростом параметра $Z^2/A$ , логарифм периода полураспада должен демонстрировать строгую линейную зависимость от этого параметра:

\log_{10} T_{1/2} = c_1 + c_2 \frac{Z^2}{A}

В этой формуле $T_{1/2}$ — период полураспада, $Z$ — заряд ядра (число протонов), $A$ — массовое число (общее число нуклонов), а $c_1$ и $c_2$ — эмпирические константы.

Практический пример: при переходе от урана-238 ( $Z^2/A \approx 35.6$ ) к плутонию-240 ( $Z^2/A \approx 36.8$ ) параметр делимости увеличивается, эффективный барьер снижается, и период полураспада резко падает с $10^{16}$ лет до $1.1 \times 10^{11}$ лет.

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

Но по мере синтеза всё более тяжелых элементов (калифорний, фермий) стройная картина рухнула. Экспериментальные значения периодов полураспада отклонялись от предсказаний капельной модели на десятки порядков. Более того, для одного и того же элемента изотопы с нечетным числом нейтронов жили в тысячи раз дольше своих четно-четных соседей (эффект нечетного нуклона), что гладкая капля жидкости объяснить в принципе не могла.

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

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

В 1962 году в лаборатории Объединенного института ядерных исследований в Дубне группа физиков под руководством С. М. Поликанова облучала изотопы плутония и урана ионами неона и кислорода. Среди продуктов реакции исследователи зафиксировали спонтанное деление с периодом полураспада всего 14 миллисекунд. Анализ показал, что делится ядро америция-242. Но согласно устоявшимся расчетам вероятности туннелирования, америций-242 в основном состоянии должен был делиться с периодом около $10^{14}$ лет. Разница почти в 24 порядка означала, что макроскопическая теория туннелирования формы, блестяще работавшая для урана и плутония, столкнулась с явлением, которое в принципе не могла описать.

Парадокс 14 миллисекунд

Обнаруженное ядро находилось не в основном, а в возбужденном (изомерном) состоянии, которое получило обозначение $^{242f}\text{Am}$ (индекс f от англ. fission — деление). Само по себе существование возбужденных состояний ядер не было новостью. Кризис заключался в сочетании двух фактов: способа распада этого изомера и его времени жизни.

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

Поместим ядро в возбужденное состояние внутри этой потенциальной ямы. У него есть два пути:

Электромагнитный переход (гамма-распад). Ядро сбрасывает излишек энергии, излучая гамма-квант, и возвращается в основное состояние. Характерное время такого процесса для типичных изомеров колеблется от наносекунд до микросекунд.
Деление. Если энергия возбуждения выше барьера, деление происходит мгновенно (за ядерное время $\sim 10^{-20}$ с). Если энергия ниже барьера, ядро туннелирует. Чем выше ядро "поднято" в потенциальной яме, тем барьер у́же и ниже, а значит, проницаемость экспоненциально возрастает.

Если $^{242f}\text{Am}$ живет 14 миллисекунд, значит, его энергия возбуждения лежит ниже вершины барьера (иначе он разделился бы мгновенно). Но почему тогда он не сбрасывает энергию через гамма-квант за наносекунды? Гамма-распад должен был опередить подбарьерное туннелирование в миллиарды раз, сделав спонтанное деление из этого состояния ненаблюдаемым.

Чтобы изомер мог жить 14 миллисекунд и при этом распадаться именно путем деления, он должен быть заперт в какой-то «ловушке», из которой путь назад (в основное состояние) сильно затруднен, а путь вперед (к делению) — облегчен. Классическая капельная парабола с одним барьером такой ловушки не имеет.

Асимметрия осколков: след структуры в точке разрыва

Проблема изомеров стала самым ярким, но не единственным провалом капельной модели. Второй фундаментальный кризис назревал в анализе масс осколков, на которые распадается ядро.

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

Радиохимический анализ показал прямо противоположную картину.

При спонтанном (и тепловом) делении тяжелых ядер, таких как уран или плутоний, ядро почти никогда не делится пополам. График выхода осколков имеет два ярко выраженных пика: легкий осколок группируется вокруг массы $A \approx 95$ , а тяжелый — строго вокруг $A \approx 140$ .

Более того, при переходе от урана-233 к калифорнию-252 масса легкого осколка растет, а пик тяжелого осколка остается неизменным — он «заякорен» на массе 140. Капельная макроскопическая физика не имеет механизмов для фиксации абсолютной массы осколка.

Разгадка кроется в магических числах. Осколок с массой около 140 близок к конфигурации изотопа олова-132, который имеет $Z=50$ (число протонов) и $N=82$ (число нейтронов). Это дважды магическое ядро с полностью заполненными протонными и нейтронными оболочками. Оно обладает повышенной энергией связи и исключительной стабильностью. Формирующееся ядро «чувствует» эту квантовую выгоду еще до момента физического разрыва шейки.

Необходимость синтеза моделей

К середине 1960-х годов стало очевидно, что макроскопическое описание деформации исчерпало себя.

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

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

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

Метод оболочечных поправок Струтинского: формирование двугорбого барьера

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

Решение, предложенное в 1966 году, оказалось элегантным математическим трюком, который навсегда изменил ландшафт ядерной физики.

Проблема двойного счета

К середине 1960-х годов физики понимали, что полная энергия ядра зависит как от его макроскопических свойств (поверхностное натяжение, кулоновское отталкивание), так и от микроскопических квантовых эффектов (заполнение нуклонных оболочек).

Самый очевидный путь синтеза — просто сложить энергию жидкой капли $E_{mac}$ и сумму энергий одночастичных квантовых уровней $\sum \epsilon_i$ , рассчитанных в каком-либо потенциале. Однако этот подход приводил к катастрофическим ошибкам в расчетах масс.

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

Теорема Струтинского и процедура сглаживания

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

E_{total} = E_{mac} + \delta E_{shell} + \delta E_{pairing}

Где $E_{mac}$ — энергия по капельной модели, $\delta E_{pairing}$ — поправка на спаривание нуклонов, а $\delta E_{shell}$ — та самая оболочечная поправка. Например, для ядра урана-238 капельная энергия составляет сотни МэВ, а оболочечная поправка всего около 2–3 МэВ, но именно она определяет сложную форму барьера деления.

Суть метода Струтинского заключается в определении $\delta E_{shell}$ как разности между истинной (дискретной) суммой энергий уровней и «сглаженной» суммой:

\delta E_{shell} = \sum_{i=1}^N \epsilon_i - \int_{-\infty}^{\lambda} \epsilon \tilde{g}(\epsilon) d\epsilon

Здесь $\epsilon_i$ — энергии дискретных уровней, $\lambda$ — энергия Ферми, а $\tilde{g}(\epsilon)$ — сглаженная плотность уровней. Представьте, что вы измеряете высоту ступенек на лестнице (дискретные уровни). Разность между реальной высотой каждой ступеньки и гладким пандусом, проведенным через них (сглаженная плотность), и даст вам локальную поправку.

Чтобы получить $\tilde{g}(\epsilon)$ , Струтинский применил математическую процедуру размазывания каждого дискретного уровня в распределение Гаусса с определенной шириной $\gamma$ (параметром сглаживания). К гауссиану добавляются полиномы Эрмита, чтобы скорректировать кривизну и гарантировать, что при сглаживании не потеряется ни один нуклон.

Оболочечная поправка $\delta E_{shell}$ — это не абсолютная энергия квантовых уровней, а мера их локальной флуктуации относительно равномерного распределения. Если уровни сгруппированы плотно (оболочка заполнена), поправка отрицательна (ядро стабильнее). Если уровни разрежены — поправка положительна.

Деформированные оболочки Нильссона

В сферическом ядре оболочечная поправка максимальна по модулю (и отрицательна) для магических чисел (50, 82, 126). Но при делении ядро меняет форму. Что происходит с квантовыми уровнями при вытягивании ядра?

Свен Гёста Нильссон показал, что при деформации вырожденные сферические уровни расщепляются. Одни уровни идут вверх по энергии, другие — вниз. В результате при определенных, сильно вытянутых формах ядра (например, при отношении осей $2:1$ ) уровни снова группируются, образуя деформированные магические числа.

Это означает, что оболочечная поправка $\delta E_{shell}$ не является константой. Она осциллирует по мере того, как ядро деформируется на пути к разрыву. В одних точках деформации она уменьшает потенциальную энергию, в других — увеличивает.

Рождение двугорбого барьера

Теперь соединим два компонента вместе.

Макроскопическая капельная энергия $E_{mac}$ дает гладкую, широкую параболу барьера деления, достигающую максимума в классической седловой точке.
Микроскопическая поправка $\delta E_{shell}$ накладывает на эту параболу осциллирующую волну с амплитудой около 2–3 МэВ.

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

Так возникает двугорбый барьер деления. Он состоит из:

Внутреннего барьера (Барьер А): преодоление начального сопротивления сферической формы.
Второго минимума (изомерной ямы): состояния с сильной деформацией (отношение осей $\approx 2:1$ ), где ядро временно стабилизируется благодаря деформированным оболочкам.
Внешнего барьера (Барьер Б): последнего препятствия перед необратимым спуском к точке разрыва.

Разрешение парадокса изомеров

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

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

Почему изомер Am-242 делится за 14 миллисекунд, а для основного состояния этот процесс занял бы сотни триллионов лет? Из основного состояния ядру нужно туннелировать сквозь весь массивный двугорбый барьер. А из второго минимума ядру остается преодолеть только внешний Барьер Б, который значительно тоньше и ниже. Вероятность туннелирования возрастает на десятки порядков.

Почему изомер не сбрасывает энергию через гамма-квант, возвращаясь в основное состояние? Для этого ядру пришлось бы изменить форму обратно на сферическую, то есть туннелировать «назад» сквозь внутренний Барьер А. Для многих актиноидов Барьер А обычно выше и толще Барьера Б, поэтому ядру энергетически выгоднее разделиться, чем вернуться в исходную форму.

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

Микроскопический подход: теория Хартри-Фока с эффективными силами и адиабатическое приближение

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

Переход к полностью микроскопическому описанию деления требует решения квантовой задачи многих тел для системы из тяжелого ядра (например, 240 нуклонов). Точное решение уравнения Шрёдингера для такого количества частиц невозможно даже на современных суперкомпьютерах. Поэтому фундаментальный подход строится на серии строгих физических приближений, центральным из которых является теория среднего поля.

От парных взаимодействий к среднему полю

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

\hat{H} = \sum_{i} \hat{T}_i + \frac{1}{2} \sum_{i \neq j} \hat{V}_{ij}

Здесь $\hat{T}_i$ — кинетическая энергия $i$ -го нуклона, а $\hat{V}_{ij}$ — потенциал парного взаимодействия между нуклонами $i$ и $j$ . Множитель $1/2$ исключает двойной учет пар. В качестве практического примера масштаба проблемы: для ядра урана-238 количество таких парных связей превышает 28 тысяч. Метод Хартри-Фока (Hartree-Fock, HF) радикально упрощает задачу: он заменяет хаотичную сеть индивидуальных взаимодействий на концепцию самосогласованного среднего поля.

Предполагается, что каждый нуклон движется не в поле тысяч отдельных соседей, а в едином усредненном потенциале $U_{HF}$ , который создается всеми остальными нуклонами. Гамильтониан разбивается на сумму одночастичных гамильтонианов:

\hat{H}_{HF} = \sum_{i} (\hat{T}_i + U_{HF}(\mathbf{r}_i))

Здесь $U_{HF}(\mathbf{r}_i)$ — усредненный потенциал (среднее поле), действующий на $i$ -й нуклон со стороны всех остальных. Практический пример: представьте движение человека в плотной толпе. Вместо того чтобы рассчитывать столкновения с каждым отдельным прохожим (парные взаимодействия), мы описываем его движение сквозь общий усредненный «поток» людей (среднее поле).

Проблема в том, что сам потенциал $U_{HF}$ зависит от того, как распределены нуклоны, а распределение нуклонов (их волновые функции) определяется потенциалом $U_{HF}$ . Возникает замкнутый круг, который разрывается методом итераций:

Задается начальное (пробное) распределение плотности нуклонов.
На его основе вычисляется форма среднего поля $U_{HF}$ .
Решается уравнение Шрёдингера для нуклонов в этом поле, получаются новые волновые функции.
Из новых волновых функций строится новая плотность.
Цикл повторяется до тех пор, пока плотность на текущем шаге не совпадет с плотностью на предыдущем. Это состояние называется самосогласованным.

Проблема жесткого кора и эффективные силы

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

Причина кроется в природе сильного взаимодействия. На расстояниях менее 0.4 фемтометра между нуклонами возникает мощное отталкивание — так называемый «жесткий кор» (hard core). В свободном пространстве нуклоны просто отскакивают друг от друга, но в плотной ядерной среде метод Хартри-Фока, оперирующий независимыми частицами, не может адекватно описать такие резкие короткодействующие корреляции.

Чтобы сохранить элегантность метода среднего поля, физики перешли к использованию эффективных взаимодействий.

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

Наиболее популярными в физике деления стали два типа эффективных сил:

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

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

Метод Хартри-Фока с ограничением (Constrained Hartree-Fock)

Стандартный метод Хартри-Фока ищет абсолютный минимум энергии системы. Если запустить самосогласованный цикл для ядра плутония-240, алгоритм сойдется к основному состоянию — дну первого минимума потенциального барьера. Но для описания спонтанного деления нам нужно знать энергию ядра при сильных деформациях, вплоть до точки разрыва.

Чтобы заставить ядро деформироваться вопреки его желанию оставаться в минимуме энергии, вводится математическое принуждение — метод Хартри-Фока с ограничением (CHF).

В гамильтониан добавляется оператор ограничения (например, оператор квадрупольного момента $\hat{Q}_{20}$ , отвечающий за вытягивание ядра) с множителем Лагранжа $\lambda$ . Мы минимизируем новый функционал:

\delta \langle \Psi | \hat{H} - \lambda \hat{Q}_{20} | \Psi \rangle = 0

Здесь $\delta$ означает поиск минимума, $\Psi$ — многочастичная волновая функция ядра, $\hat{H}$ — исходный гамильтониан, $\hat{Q}_{20}$ — оператор деформации, а $\lambda$ — множитель Лагранжа. Множитель $\lambda$ имеет ясный физический смысл — это внешняя «сила», которую нужно приложить к ядру, чтобы удержать его в заданном деформированном состоянии. Практический пример: представьте, что вы растягиваете тугую пружину эспандера. Множитель $\lambda$ — это усилие ваших рук, необходимое для удержания пружины на заданной длине (заданной деформации $\langle \hat{Q}_{20} \rangle$ ) вопреки ее стремлению сжаться обратно в минимум энергии.

Алгоритм построения барьера деления выглядит так:

Выбирается сетка значений деформации (от сферы до гантели).
Для каждой точки сетки подбирается такое значение $\lambda$ , чтобы ожидаемое значение деформации $\langle \hat{Q}_{20} \rangle$ равнялось заданному.
Проводится полный самосогласованный цикл Хартри-Фока.
Вычисляется полная микроскопическая энергия ядра $E$ в этой точке.

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

Адиабатическое приближение

Построение барьера деления методом CHF неявно опирается на фундаментальное допущение о соотношении скоростей процессов внутри ядра. Это допущение называется адиабатическим приближением.

В процессе спонтанного деления происходят два типа движения:

Быстрое одночастичное движение нуклонов по своим орбиталям внутри потенциальной ямы (характерное время $\sim 10^{-22}$ с).
Медленное коллективное изменение формы самого ядра, то есть макроскопическая деформация (характерное время подбарьерного туннелирования $\sim 10^{-20}$ с и более).

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

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

Именно адиабатическое приближение дает нам право использовать статическую потенциальную поверхность $E(q)$ , рассчитанную методом Хартри-Фока, как барьер, сквозь который туннелирует система. Если бы ядро деформировалось слишком быстро, нуклоны не успевали бы перестраиваться, их орбитали пересекались бы (эффект Ландау-Зинера), система переходила бы в возбужденные состояния, и эффективный барьер деления резко возрос бы.

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

Динамика процесса: тензор инерции, траектории наименьшего действия и современные вычислительные методы

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

От эффективной массы к тензору инерции

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

Деление — это процесс, описываемый не одной переменной (удлинением), а целым набором коллективных координат. Обозначим их как $q_i$ . Это может быть удлинение ядра, толщина шейки, параметр асимметрии масс (формирование тяжелого и легкого осколков). В многомерном пространстве кинетическая энергия коллективного движения $T$ уже не описывается одним числом. Она задается матрицей — тензором инерции:

T = \frac{1}{2} \sum_{i,j} B_{ij}(q) \dot{q}_i \dot{q}_j

Здесь $B_{ij}(q)$ — компоненты массового тензора, зависящие от текущей формы ядра $q$ .

Тензор инерции $B_{ij}$ показывает, какую массу «ощущает» ядро при попытке изменить форму вдоль координаты $i$ , если оно одновременно движется вдоль координаты $j$ .

Почему инерция зависит от формы? Вспомним адиабатическое приближение: при медленном растяжении ядра нуклоны остаются на своих орбиталях. Но по мере деформации энергетические уровни нуклонов пересекаются. В точках пересечений (или сильного сближения) уровней ядру приходится кардинально перестраивать свою внутреннюю волновую функцию. В этих зонах сопротивление деформации резко возрастает — компоненты тензора $B_{ij}$ демонстрируют гигантские пики, превышающие гидродинамические оценки в десятки раз.

Для расчета этих пиков сегодня применяется модель крэнкинга (cranking model) или более строгий метод адиабатического зависящего от времени Хартри-Фока (ATDHF). Суть модели крэнкинга изящна: мы математически «тянем» ядро с заданной скоростью и смотрим, как квантовая система сопротивляется этому внешнему возмущению, вычисляя отклик нуклонных орбиталей.

Многомерное туннелирование и интеграл действия

Имея многомерную потенциальную энергию $V(q)$ и тензор инерции $B_{ij}(q)$ , мы можем рассчитать вероятность спонтанного деления. Для этого используется многомерное обобщение квазиклассического приближения (ВКБ).

Вероятность туннелирования определяется интегралом действия $S(L)$ , вычисленным вдоль некоторой траектории $L$ в пространстве деформаций:

S(L) = \int_{L} \sqrt{2 B_L(q) (V(q) - E)} \, ds

В этой формуле $E$ — энергия основного состояния ядра, $ds$ — элемент длины дуги вдоль траектории, а $B_L(q)$ — эффективная инерция вдоль выбранного пути $L$ , которая собирается из проекций тензора $B_{ij}$ на касательную к траектории.

Главный вопрос динамики деления: какую именно траекторию $L$ выберет ядро?

Согласно принципам квантовой механики, доминирующей будет та траектория, для которой интеграл действия $S(L)$ минимален. Именно она даст наибольший вклад в вероятность туннелирования.

Статика против динамики: почему ядро «срезает углы»

Исторически физики искали так называемую статическую траекторию. Это путь, который проходит через минимумы потенциальной энергии и самые низкие седловые точки (перевалы). Статическая траектория минимизирует подкоренное выражение $V(q) - E$ . Кажется логичным: чтобы легче пройти сквозь гору, нужно найти самый низкий перевал.

Но квантовое ядро «видит» не только высоту барьера $V(q)$ , но и массу $B_L(q)$ .

Представьте, что вы идете через горный хребет. Самый низкий перевал проходит через глубокое вязкое болото (огромная эффективная инерция $B_L$ ). Каждый шаг дается с колоссальным трудом. Рядом есть другой перевал — он физически выше (больше $V(q)$ ), но путь пролегает по твердой скале (малая инерция $B_L$ ). Квантовая система, минимизируя общее действие $S(L)$ , выберет путь по скале.

Этот путь называется динамической траекторией.

Динамическая траектория наименьшего действия часто отклоняется от дна потенциальной долины. Ядро способно «срезать углы» в пространстве деформаций, проходя через области с более высокой потенциальной энергией, если это позволяет обойти зоны с аномально высокой эффективной инерцией (те самые пики от пересечения нуклонных уровней).

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

Современные вычислительные методы (GCM)

Расчет многомерного действия — сложнейшая вычислительная задача. Сегодня стандартом де-факто в микроскопической теории деления является метод генераторных координат (Generator Coordinate Method, GCM) в сочетании с гауссовым приближением перекрытия (GOA).

Метод Хартри-Фока, который мы обсуждали ранее, дает нам одно лучшее состояние для каждой формы ядра. Но квантовая механика говорит, что истинное состояние ядра — это суперпозиция (смешивание) различных форм. Ядро не имеет строго фиксированной деформации, оно постоянно «дышит», испытывая квантовые флуктуации.

Метод GCM математически собирает волновую функцию ядра как интеграл по всем возможным деформациям (генераторным координатам). Применяя приближение GOA к уравнениям GCM, физики получают дифференциальное уравнение типа уравнения Шрёдингера для коллективного движения. Из этого уравнения естественным образом, без дополнительных искусственных допущений, выводятся и потенциальная энергия $V(q)$ , и тензор масс $B_{ij}(q)$ , и квантовые поправки к энергии нулевых колебаний.

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