The magicity, the radii of neutron orbits 1f7/2, 2p3/2 and halo-like structure of 52,54Ca nuclei
- Authors: Bespalova O.V.1, Klimochkina A.A.2
-
Affiliations:
- Lomonosov Moscow State University, Skobeltsyn Institute of Nuclear Physics
- Lomonosov Moscow State University
- Issue: Vol 87, No 2 (2024)
- Pages: 102-107
- Section: МАТЕРИАЛЫ 73-Й МЕЖДУНАРОДНОЙ КОНФЕРЕНЦИИ ПО ЯДЕРНОЙ ФИЗИКЕ. Ядра. Теория
- URL: https://journals.rcsi.science/0044-0027/article/view/265627
- DOI: https://doi.org/10.31857/S0044002724020065
- EDN: https://elibrary.ru/KRSPAP
- ID: 265627
Cite item
Full Text
Abstract
The evolution of neutron single-particle spectra of isotones with N = 32 and 34 new magic neutron numbers in the region 16 ≤ Z ≤ 32 was calculated in the dispersive optical model. It was shown that the minimum of the difference between the Fermi energy and the half-sum of the energy levels of the last predominantly occupied state and the first predominantly unoccupied state is achieved in the magic isotones with N = 32 and 34. The calculated root-mean-square radius of the neutron halo-like state 2p3/2 in double magic 52Ca nucleus exceeded the radius of the underlying 1f7/2 state by 0.8 fm. It is consistent with the recent experimental data and theoretical predictions that explain ‟unexpectedly” large root-mean-square charge radius of this nucleus.
Full Text
1. ВВЕДЕНИЕ
Экспериментальные данные об энергии возбуждения [1] свидетельствуют о новой магичности изотопов Са с числами нейтронов N = 32, 34 (рис. 1, а), которой посвящено большое количество исследований (см., например, [2–4] и ссылки в них). Для традиционных магических чисел N = 20 и 28 одновременно с увеличенными значениями наблюдается уменьшение зарядовых радиусов rch соответствующих изотопов по сравнению с соседними. Однако для N = 32 в изотопах Са такого уменьшения не происходит [5], и, начиная с 48Са, радиусы rch демонстрируют неуклонный рост (рис. 1, б). “Неожиданно” большие радиусы rch нейтронно-избыточных изотопов Са, полученные в [5], поставили вопрос о степени проявления магических свойств изотопов 52Са с числом N = 32.
Рис. 1. Экспериментальные энергии (а) и среднеквадратичные зарядовые радиусы rch (б) четных изотопов Са.
Рост радиуса rch в ядрах с N > 28 был приписан [6] галоподобным особенностям структуры нейтронных орбит 2p3/2 и 2p1/2, которые заполняются в изотопах Са с N > 28, и сопровождающим это заполнение расширением распределения протонов, ‟следующих” за нейтронами. В [6] была предсказана разность в 0.7 Фм между среднеквадратичными радиусами rrms нейтронных орбит 1f7/2 и 2p3/2, которая соответствует наблюдаемому росту rch для 52Са при сохранении этим ядром дважды магической природы. Недавние экспериментальные данные, полученные в реакции выбивания нейтрона из 52Са [7], показали превышение rrms нейтронов на орбите 2p3/2 на 0.61(23) Фм по сравнению с орбитой 1f7/2. Так, данные составили 4.74(18) Фм и 4.13(14) Фм для нейтронных состояний 2p3/2 и 1f7/2 соответственно. Расчеты по методу Хартри—Фока—Боголюбова с взаимодействием Скирма (SKM), выполненные в [7], привели к rrms = 4.12 Фм для состояния 1f7/2 и к 4.49 Фм для состояния 2p3/2. Если первое значение хорошо согласуется с экспериментальным, то второе заметно меньше его. В настоящей работе выполнены расчеты эволюции нейтронных одночастичных спектров изотонов с магическим числом нейтронов N = 32 и 34 в дисперсионной оптической модели (ДОМ), определены среднеквадратичные радиусы rrms нейтронных состояний 1f7/2, 2p3/2 и распределение нейтронной плотности в изотопах 52,54Са, согласующиеся с экспериментальными данными [7] о галоподобной структуре состояния 2p3/2 в ядре 52Са.
2. ОСНОВЫ МОДЕЛИ И ОПРЕДЕЛЕНИЕ ЕЕ ПАРАМЕТРОВ
Дисперсионная оптическая модель [8], развитая в конце прошлого века, в отличие от традиционной оптической модели опирается на дисперсионное соотношение, связывающее действительную и мнимую части потенциала ДОМ:
(1)
где Р обозначает главное значение. После перехода от нелокального потенциала к локально эквивалентному действительная центральная часть нейтронного потенциала V(r, E) представляется в виде суммы плавно зависящей от энергии составляющей хартри-фоковского (ХФ) типа и дисперсионной составляющей, которая резко зависит от энергии вблизи энергии Ферми EF:
(2)
где f(r, ri, ai) — функция Вудса–Саксона. Индексы HF, s, d относятся соответственно к ХФ, объемной и поверхностной составляющим потенциала ДОМ. Дисперсионная составляющая эффективно учитывает корреляции, которые испытывает нуклон в ядре, как распределенные по объему, так и сконцентрированные на поверхности, и приводит к уменьшению частично-дырочной щели G между последним преимущественно занятым состоянием и первым преимущественно свободным состояниями нуклона в ядре. В настоящей работе мнимая часть дисперсионного оптического потенциала (ДОП) предполагалась симметричной относительно энергии EF и была параметризована следующими выражениями для объемной Ws и поверхностной Wd частей:
(3)
где w1, w2, d1, d2, d3 — параметры. Энергия Ферми EF определялась по данным об энергиях отделения Sn нейтрона от ядра с (N, Z) и (N + 1, Z) [9]:
(4)
Энергетическая зависимость ХФ-составляющей определялась выражением
(5)
Параметр находился из условия согласия суммарного числа нейтронов в состояниях с квантовыми числами nlj с числом N изотопа. Для этого вероятность заполнения Nnlj одночастичных орбит определялась по формуле теории Бардина–Купера–Шриффера (БКШ):
(6)
с эмпирическим значением параметра щели спаривания:
(7)
Таким образом, условие Nn = N соответствует согласию с эмпирическим значением энергии Ферми.
Энергии Enlj, использованные в (6), вычислялись при решении уравнения Шредингера с центральной действительной частью нейтронного ДОП
(8)
где
(9)
Радиальная часть unlj(r) полной волновой функции (r) была скорректирована в целях учета эффекта нелокальности:
(10)
где отношение ХФ эффективной массы нуклона к его полной массе m определяется выражением
(11)
Коэффициент Cnlj находился при нормировке на единицу. Среднеквадратичные радиусы rms нейтронных состояний определялись согласно соотношению
(12)
а распределение нейтронной плотности — в одночастичном подходе, аналогично [10]:
(13)
3. ЭВОЛЮЦИЯ МАГИЧЕСКИХ ЧИСЕЛ N = 32, 34
Для расчета характеристик нейтронных одночастичных состояний в качестве параметров мнимого, спин-орбитального потенциалов и геометрических параметров rHF, aHF были взяты значения, экстраполированные в соответствии с параметрами глобальных потенциалов Конинга–Делароша (KD) [11]. В результате глубина потенциала VHF(EF) являлась единственным свободным параметром, который находился из условия согласия Nn с числом N изотопа с точностью не хуже 0.5%. Погрешности определения параметра γ практически не влияли на энергии нейтронных состояний вблизи EF. Этот параметр для 52Са полагался равным 0.42. Расчетные одночастичные нейтронные энергии Enlj изотонов с магичным N = 32, вычисленные с такими глобальными параметрами ДОП, показаны на рис. 2 в сравнении с экспериментальными данными [12, 13]. Обращает на себя внимание взаимное расположение энергии Ферми EF и энергий состояний 2p, которые являются последним занятым и первым свободным состояниями в рамках модели независимых частиц (МНЧ). Значение EF соответствует почти половине суммарной энергии уровней 2p3/2 и 2p1/2 только в изотонах Ca и Ti. Это иллюстрирует вставка к рис. 2, на которой представлено различие между энергией EF и средним положением состояний 2р в изотонах с N = 32. Оно достигает наименьших значений в 0.1 МэВ в 52Са и 54Ti, что составляет около 10% от щели G между энергиями состояний 2p3/2 и 2p1/2.
Рис. 2. Нейтронные одночастичные энергии изотонов с N = 32. Светлые значки, соединенные сплошными линиями, — расчет с ДОП, штрихпунктирная линия — энергия EF, штриховые — энергии –Sn(A), –Sn(A + 1). Темные значки — результат совместной оценки данных реакций срыва и подхвата нейтрона на одном и том же ядре [12, 13].
В ядрах с традиционными магическими числами последнее состояние ниже EF с энергией практически полностью заполнено, а первое состояние выше EF с энергией — свободно. Поэтому для таких ядер энергии Sn (с противоположным знаком) отделения нуклона от него и ядра с A + 1 близки соответственно к энергиям и , энергия EF (4) — к их среднему, а основные характеристики ядра хорошо описываются МНЧ. В ядрах 52Са, 54Ti с новым магическим числом N = 32 эмпирические энергии отделения нейтрона –Sn от ядер с (N = 32, Z) и (N = 33, Z) отличаются от расчетных по ДОМ энергий Enlj состояний 2p3/2 и 2p1/2, однако среднее между энергиями Sn близко к среднему между рассчитанными в ДОМ энергиями Enlj. Такое расположение энергии EF характерно и для других ядер с традиционными и новыми магическими числами [14] и может рассматриваться в качестве маркера магичности. Для того чтобы при этом ядро проявляло себя как магическое, необходимо наличие значительной энергетической щели G между и .
В ДОМ не учитывается вклад тензорного взаимодействия, который важен для количественного объяснения магических свойств 52Ca и 54Ti [4]. За счет тензорных сил протоны заполненной оболочки 1d3/2 в изотопах Са притягивают нейтроны подоболочки 2p3/2 и отталкивают — 2p1/2, расширяя тем самым нейтронную щель N = 32 по сравнению с расчетом по ДОМ. В 54Ti такое действие тензорного взаимодействия ослабевает, так как начинает заполняться протонная оболочка 1f7/2. Несмотря на отсутствие тензорного взаимодействия, ДОМ хорошо описывает динамику необходимого условия для возникновения магического числа N = 32. При увеличении Z среднее между энергиями состояний 2p3/2 и 2p1/2 удаляется от EF (см. рис. 2). При этом минимум разности ∆F достигается для двух ядер: 52Са и 54Ti. Такая картина соответствует исчезновению магических свойств числа N = 32 (рис. 3) с увеличением Z, которое можно проследить по данным о . Высокие значения этой энергии характерны для ядер 52Са и 54Ti. И хотя при переходе от дважды магического ядра 52Са к 54Ti наблюдается существенное уменьшение энергии , ее значение в 54Ti остается выше, чем в изотоне 60Ni с магическим числом Z = 28.
Рис. 3. Экспериментальные энергии изотонов с N = 32 (квадраты) и N = 34 (кружки).
Аналогично можно проследить динамику нового магического числа N = 34. В МНЧ числу N = 34 соответствует заполнение всех подоболочек ниже 2p1/2 включительно, вышележащая подоболочка 1f5/2 свободна. В стабильных изотонах 62Ni и 64Zn ситуация другая. Согласно результатам совместной оценки данных реакций срыва и подхвата нейтрона [13] состояние 2p1/2 расположено выше энергии EF и заполнено лишь на четверть, что не предполагает замыкания подоболочки N = 34 в этих ядрах. Согласно расчету по ДОМ (рис. 4) при уменьшении Z уровень 2p1/2 опускается ниже EF, так что последняя оказывается между уровнями 2p1/2 и 1f5/2. В этом проявляется магичность числа N = 34 в изотонах 54Ca и 52Ar, которая находит свое отражение в энергиях (см. рис. 3). Отличие энергии EF от среднего значения энергий последнего преимущественно занятого и первого преимущественно свободного состояний 2p1/2 и 1f5/2 показано на врезке к рис. 4 для N = 34. Оно достигает наименьшего значения ≈0.2 МэВ в 54Са и 52Ar, что составляет около 17% от щели между энергиями состояний 2p1/2 и 1f5/2. Таким образом, в расчетных одночастичных спектрах магичность числа N = 34 менее выражена по сравнению с N = 32. Это качественно согласуется с уменьшением экспериментальной энергии (см. рис. 1, а) в 54Са по сравнению с 52Са. Значение ∆F для соседнего с 54Са изотона 56Ti с N = 34 по абсолютной величине близко к значению для 56Cr с N = 32, для которого магические свойства числа N = 32 уже не проявляются. Такая картина согласуется с энергиями изотонов с N = 32, 34 в ядрах 56Ti и 56Cr, которые отличаются друг от друга лишь на 10% (см. рис. 3).
Рис. 4. То же, что и на рис. 2 для изотонов с N = 34.
4. РАДИУСЫ НЕЙТРОННЫХ СОСТОЯНИЙ 1f7/2 И 2p3/2
Ранее в [15] нами был выполнен расчет нейтронных и протонных радиусов изотопов Са по ДОМ. Расчетный протонный радиус 48Са оказался меньше, чем у соседних изотопов, а при дальнейшем росте числа N он продемонстрировал рост, качественно согласующийся с экспериментальными данными [5].
Среднеквадратичные радиусы rrms 52Са для нейтронных орбит 1f7/2 и 2p3/2, вычисленные нами с теми же параметрами, которые определены при вычислении энергий нейтронных орбит, равны соответственно 4.22 и 5.04 Фм. Разность этих rrms составляет 0.82 Фм, что согласуется с экспериментальными данными 0.61(23) Фм с учетом их погрешности. Радиус rrms состояния 1f7/2 находится в согласии с экспериментальным радиусом = 4.13(14) Фм, в то время как rrms состояния 2p3/2 завышен по сравнению с = 4.74(18) Фм.
Индивидуальные особенности ядер можно учесть, изменяя, в частности, геометрические параметры ХФ-составляющей ДОП. Следует отметить, что зависимость глобального параметра диффузности систематики KD от массового числа А имеет значительную неопределенность [16], что может быть связано с отсутствием зависимости этого параметра от нейтрон-протонной асимметрии (N – Z) / A [17]. Для нового магического ядра 52Са параметр диффузности aHF был уменьшен до 0.58 Фм (KD aV = 0.67 Фм). Такой параметр диффузности aHF соответствует более резкой границе потенциала и может быть сопоставлен с параметром aHF для 48Са [12], являющегося кором для ядра 52Са. С параметрами aHF = 0.58 Фм и rHF = 1.2 Фм (KD rV = 0.195 Фм) были получены значения rrms = 4.13 Фм и 4.94 Фм для состояний 1f7/2 и 2p3/2 в 52Са. Последнее значение лучше согласуется с экспериментальными данными [5]. В ядре 54Са, согласно МНЧ, состояние 2p1/2 — последнее заполненное. С измененными по сравнению с KD параметрами (aHF = 0.58 Фм, rHF = 1.2 Фм, VHF(EF) = 46.9 МэВ) для состояния 2p1/2 в 54Са получено значение rrms = 5.48 Фм, бóльшее, чем для состояния 2p3/2.
Существенное улучшение согласия rrms состояния 2p3/2 в 52Са с экспериментом возможно при учете тензорного взаимодействия, притягивающего нейтроны в этом состоянии и, следовательно, уменьшающего rrms. Учет тензорного вклада во взаимодействие Скирма (SGII + T) в приближении ХФ + БКШ для 52Са привел к уменьшению радиуса rrms орбит 1f7/2 и 2p3/2 соответственно на 0.09 и 0.1 Фм. Отметим, что и в этих расчетах полученный радиус rrms = 4.44 Фм нейтронного состояния 2p3/2, так же как и в [5], занижен по сравнению с экспериментальным.
На рис. 5, а представлены нейтронные плотности ρn(r) ядер 52,54Са и 48Са, вычисленные по ДОМ с параметрами aHF = 0.58 Фм и rHF = 1.2 Фм. Как видно из рисунка, для ядер 52Са и 54Са характерен более протяженный хвост распределения плотности, напоминающий структуру гало. В [18] предложен модельно-независимый метод выделения области, в которой нуклоны гало пространственно не коррелируют с кором. Метод основан на изменении наклона плотности ρn(r) на поверхности ядра, где она сформирована как нуклонами кора, так и гало. Изменение наклона приводит к появлению характерного максимума при r = rмакс в радиальном распределении второй логарифмической производной плотности по радиусу (рис. 5, б). Согласно [18] при rмакс плотность кора и гало дают одинаковый вклад в суммарную плотность ядра. Функция для 52Са имеет слабо выраженный максимум при rrms = 6.8 Фм, в отличие от 48Са. Аналогичный максимум характерен и для ядра 54Ca. Этот результат находится в согласии с представлением о галоподобной структуре изотопов 52, 54Са.
Рис. 5. Нейтронная плотность ρn(r) (а) и вторая производная от ее логарифма (б) для изотопов 48Са (штриховая кривая),52Са (сплошная) и 54Са (точечная).
Мы оценили также количество нейтронов N(rмакс) и N(r0) в ядрах 52,54Ca, которые расположены в области r ≥ rмакс и в области гало r ≥ r0. Последняя определена в [18] как область за пределами радиуса r0, при котором плотность гало в 10 раз превышает плотность кора. Значение r0 для 52,54Са близко к 9.5 Фм. В ядре 52Са в области r ≥ rмакс расположено около N(rмакс) ≈ 0.5 нейтрона, а в области r ≥ r0 всего лишь N(r0) = 0.03 ± 0.01 нейтронов. Теоретическая погрешность N(r0) найдена в соответствии со способом оценки погрешности r0, предложенным в [18]. Аналогичные значения для 54Са составили N(rмакс) ≈ 0.8 нейтрона и N(r0) = 0.07 ± ± 0.03 нейтрона. Это соответствует представлению о 52,54Ca как о ядрах с галоподобной структурой, но не собственно с гало в силу малости числа N(r0).
5. ЗАКЛЮЧЕНИЕ
Выполнены расчеты, предсказывающие одночастичные нейтронные характеристики изотонов с новыми магическими числами N = 32, 34 в рамках дисперсионной оптической модели. Показано, что для изотонов 52Ca, 54Ti с N = 32 и 52Ar, 54Ca с N = 34 характерна близость середины между последним преимущественно заполненным и первым преимущественно свободным нейтронными состояниями в МНЧ к определяемой эмпирически энергии Ферми EF, что отражает магические свойства этих ядер. Для соседних изотонов эта близость быстро исчезает. Такая особенность расчетных спектров может быть полезна для предсказаний возможных кандидатов в новые магические ядра.
Рассчитанный среднеквадратичный радиус rrms последнего преимущественно заполненного нейтронного состояния 2p3/2 в 52Са превысил радиус нижележащего состояния 1f7/2 на 0.8 Фм, что согласуется с представлением о галоподобной структуре состояния 2p3/2, которая рассматривается в качестве причины “неожиданно” большого зарядового радиуса этого ядра. Однако rrms нейтронного состояния 2p3/2 в ДОМ завышен по сравнению с экспериментальным, что можно объяснить отсутствием в этой модели учета притягивающего тензорного взаимодействия нейтронов в 2p3/2 состоянии с протонами в заполненном состоянии 1d3/2.
Рассчитанные нейтронные плотности ρn(r) ядер 52,54Са и 48Са показали, что для ядер 52Са и 54Са присутствует, в отличие от 48Са, максимум второй производной от логарифма плотности , характерный для ядер со структурой гало. Однако этот максимум слабо выражен, а число нейтронов в области периферии ядра, где они пространственно не коррелированы с нуклонами кора, несущественно. Такая картина соответствует представлению лишь о галоподобных особенностях распределения нейтронов в состояниях 2p в ядрах 52Са и 54Са, которые заключаются в бóльшем среднеквадратичном радиусе этих состояний по сравнению с нижележащим состоянием 1f7/2.
Авторы благодарят аспиранта физического факультета МГУ С.В. Сидорова за предоставленные результаты выполненного им расчета радиусов нейтронных состояний 52Са в приближении ХФ + БКШ с учетом и без учета тензорного вклада во взаимодействие Скирма.
About the authors
O. V. Bespalova
Lomonosov Moscow State University, Skobeltsyn Institute of Nuclear Physics
Author for correspondence.
Email: besp@sinp.msu.ru
Russian Federation, Moscow
A. A. Klimochkina
Lomonosov Moscow State University
Email: besp@sinp.msu.ru
Faculty of Physics
Russian Federation, MoscowReferences
- База данных ENSDF, http://www.nndc.bnl.gov/ensdf
- A. Gade, R. V. F. Janssens, D. Bazin, R. Broda, B. A. Brown, C. M. Campbell, M. P. Carpenter, J. M. Cook, A. N. Deacon, D.-C. Dinca, B. Fornal, S. J. Freeman, T. Glasmacher, P. G. Hansen, B. P. Kay, P. F. Mantica, et al., Phys. Rev. C 74, 021302(R) (2006).
- D. Steppenbeck, S. Takeuchi, N. Aoi, P. Doornenbal, M. Matsushita, H. Wang, H. Baba, N. Fukuda, S. Go, M. Honma, J. Lee, K. Matsui, S. Michimasa, T. Motobayashi, D. Nishimura, T. Otsuka, et al., Nature 502, 207 (2013).
- M. Honma, T. Otsuka, B.A. Brown, and T. Mizusaki, Eur. Phys. J. A 25, 499 (2005).
- R. F. Garcia Ruiz, M. L. Bissell, K. Blaum, A. Ekström, N. Frömmgen, G. Hagen, M. Hammen, K. Hebeler, J. D. Holt, G. R. Jansen, M. Kowalska, K. Kreim, W. Nazarewicz, R. Neugart, G. Neyens, W. Nörtershäuser, et al., Nature Phys. 12, 594 (2016).
- J. Bonnard, S. M. Lenzi, and A. P. Zuker, Phys. Rev. Lett. 116, 212501 (2016).
- M. Enciu, H. N. Liu, A. Obertelli, P. Doornenbal, F. Nowacki, K. Ogata, A. Poves, K. Yoshida, N. L. Achouri, H. Baba, F. Browne, D. Calvet, F. Château, S. Chen, N. Chiga, A. Cors, et al., Phys. Rev. Lett. 129, 262501 (2022).
- C. Mahaux and R. Sartor, Adv. Nucl. Phys. 20, 1 (1991).
- M. Wang, W. J. Huang, F. G. Kondev, G. Audi, and S. Naimi, Chin. Phys. C 45, 030003 (2021).
- M. Jaminon and C. Mahaux, Nucl. Phys. A 440, 228 (1985).
- A. J. Koning and J. P. Delaroche, Nucl. Phys. A 713, 231 (2003).
- О. В. Беспалова, И. Н. Бобошин, В. В. Варламов, Т. А. Ермакова, Б. С. Ишханов, Е. А. Романовский, Т. И. Спасская, Т. П. Тимохина, ЯФ 71, 37 (2008) [O. V. Bespalova, I. N. Boboshin, V. V. Varlamov, T. A. Ermakova, B. S. Ishkhanov, E. A. Romanovsky, T. I. Spasskaya, and T. P. Timokhina, Phys. At. Nucl. 71, 36 (2008)].
- О. В. Беспалова, И. Н. Бобошин, В. В. Варламов, Т. А. Ермакова, Б. С. Ишханов, А. А. Климочкина, С. Ю. Комаров, Ч. Коура, Е. А. Романовский, Т. И. Спасская, Изв. РАН. Сер. физ. 74, 575 (2010) [O. V. Bespalova, I. N. Boboshin, V. V. Varlamov, T. A. Ermakova, B. S. Ishkhanov, A. A. Klimochkina, S. Yu. Komarov, H. Koura, E. A. Romanovsky, and T. I. Spasskaya, Bull. Russ. Acad. Sci.: Phys. 74, 542 (2010)].
- О. В. Беспалова, Е. А. Романовский, Т. И. Спасская, ЯФ 78, 123 (2015) [O. V. Bespalova, E. A. Romanovsky, and T. I. Spasskaya, Phys. At. Nucl. 78, 118 (2015)].
- О. В. Беспалова, А. А. Климочкина, ЯФ 80, 516 (2017) [O. V. Bespalova and A. A. Klimochkina, Phys. At. Nucl. 80, 919 (2017)].
- C. D. Pruitt, J. E. Escher, and R. Rahman, Phys. Rev. C 107, 014602 (2023).
- О. В. Беспалова, А. А. Климочкина, ЭЧАЯ 53, 428 (2022) [O. V. Bespalova and A. A. Klimochkina, Phys. Part. Nucl. 53, 476 (2022)].
- V. Rotival and T. Duguet, Phys. Rev. C 79, 054308 (2009).
Supplementary files
