Физики рассчитали толщину нейтронной оболочки в ядре изотопа свинца-208 — она находится в диапазоне от 0,14 до 0,20 фемтометра. Для этого учёные объединили расчёт из первых принципов, статистические методы и машинное обучение. Такая схема работает быстрее прямых симуляций и в будущем её можно будет использовать для описания как лёгких, так и тяжёлых ядер, пишут учёные в своём исследовании, опубликованном в журнале Nature Physics.
В тяжёлых ядрах, в которых число нейтронов превышает число протонов, нуклоны распределены неравномерно: избыток нейтронов концентрируется снаружи, из-за чего у ядра появляется «нейтронная кожа» — внешняя оболочка, в которой нейтронов больше, чем протонов. Этот эффект связан с обменной энергией между протонами и нейтронами, от него зависит размер и форма тяжёлых ядер. Похожее взаимодействие нейтронов происходит и в нейтронных звёздах, поэтому экспериментальные исследования и моделирование нейтронных оболочек в ядре помогут учёным разобраться и со структурой нейтронных звёзд.
При моделировании тяжёлых ядер физикам обычно приходится выбирать между точностью и скоростью расчётов. Приближённые методы, например, теория самосогласованного поля, экономят вычислительные ресурсы, но дают систематические ошибки, а более точные методы ab initio — расчёта из первых принципов — редко применяются для систем, содержащих более 100 частиц, в первую очередь из-за вычислительной сложности алгоритмов даже для суперкомпьютеров. Под термином «ab initio» в физике твёрдого тела подразумевается группа методов, в которых систему моделируют, описывая первичные квантовые принципы взаимодействия между частицами — для получения её параметров необходимо решить уравнение Шрёдингера. И, если структуру нейтронной оболочки магического изотопа кальция-48, состоящего из 28 нейтронов и 20 протонов, таким способом обсчитали в 2015 году, то значительно более тяжёлое ядро олово-132 (p = 50, n = 82) — только в 2020 году. А для свинца-208 (p = 82, n = 126) ab initio расчётов до сих пор не было.
Чтобы смоделировать столь тяжёлое ядро, учёные под руководством Байшаня Ху (Baishan Hu) из канадского ускорительного центра TRIUMF использовали эмуляторы. Методы байесовской оптимизации имитируют работу других алгоритмов и выдают соответствующие результаты, при этом не выполняя все рутинные вычисления. В данном случае эмуляторы помогли заменить решатели задачи многих тел для ab initio расчётов.
Для описания ядерных сил учёные использовали теорию хиральных возмущений, согласно которой энергия нуклонов определяется несколькими порядками взаимодействия. Дальнодействующая компонента соответствует пионному обмену между нейтроном и протоном. Ближнее взаимодействие второго порядка описывается 17 константами низких энергий, 4 из которых фиксированы, а значения 13 оставшихся учёные определили в работе.
Варьируя эти 13 параметров в ab anitio симуляциях и сравнивая результаты с известными значениями фазовых сдвигов нуклон-нуклонного рассеяния на низких энергиях (до 200 МэВ), а также энергией и радиусом дейтерия, трития, гелия-4, кислорода-16, учёные нашли 34 набора констант, которые согласуются с экспериментальными данными. При сравнении экспериментальных данных с предсказаниями моделирования, авторы учитывали экспериментальную ошибку, точность эмулятора, ошибку аппроксимации метода и модели поля. Затем модель откалибровали, определив значимость каждого набора констант с помощью функции правдоподобия.
В результате учёным удалось рассчитать энергию основного состояния, радиус точечных протонов и поляризуемость диполей для кальция-48 и свинца-208. Толщина нейтронного слоя для свинца-208 оказалась в диапазоне от 0,14 до 0,20 фемтометра — это хорошо согласуется с большинством экспериментальных работ (правда, противоречит результатам коллаборации PREX). Физики говорят, что использованную ими схему моделирования в будущем можно использовать и для расчёта параметров других ядер, поскольку полученный набор констант не зависит от числа нуклонов. Однако неясно, потребуется ли перерасчёт констант и насколько точно они будут описывать другие ядра.
Учёные отмечают, что расчёт нуклонного взаимодействия из первых принципов на основании теории хиральных возмущений можно применять как для небольших, так и для тяжёлых атомов. Принципиально ограничивать точность вычислений будут только данные рассеяния нуклонов.