Известия РАН. Серия физическая, 2019, T. 83, № 4, стр. 473-478

Роль механизма срыва тяжелого кластера в реакции 27Al(р, α)24Mg

Л. И. Галанина 1*, Н. С. Зеленская 1, В. Ю. Огнев 1

1 Федеральное государственное бюджетное образовательное учреждение высшего образования “Московский государственный университет имени М.В. Ломоносова”, Научно-исследовательский институт ядерной физики имени Д.В. Скобельцына
Москва, Россия

* E-mail: galan_lidiya@mail.ru

Поступила в редакцию 01.10.2018
После доработки 15.10.2018
Принята к публикации 19.11.2018

Полный текст (PDF)

Аннотация

Рассчитан вклад в сечение реакции 27Al(р, α)24Mg(0+, 2+) при Ep = 7.4 МэВ в интервале углов вылета α-частиц 25°–160° (лаб.) механизма срыва тяжелого кластера 23Na с начального ядра 27Al(5/2+). Расчет сечения этого механизма проведен по программе FRESCO с учетом вклада основного 3/2+ состояния 23Na. Спектроскопические амплитуды в вершинах распада 27Al(5/2+) → 23Na(3/2+) + α и 24Mg(0+, 2+) → 23Na(3/2+) + р определялись с использованием волновых функций модели Нильссона с учетом их зависимости от деформации. Рассчитанные сечения для уровней 0+ и 2+ ядра 24Mg меньше экспериментальных на два порядка. Сечение для уровня 0+ имеет подъем при θα ∼ 180° в задней полусфере, а для уровня 2+ – максимум при θα > 140°. Незначительный вклад механизма срыва тяжелого кластера связан с малыми спектроскопическими амплитудами в вершинах распада начального и конечного ядра вследствие слабого перекрывание волновых функций ядер 27Al, 23Na и 24Mg(0+, 2+) в модели Нильссона.

ВВЕДЕНИЕ

В работе [1] получены экспериментальные угловые распределения (УР) дифференциального сечения реакции 27Al(p, α)24Mg(0+, 2+) при Ep = 7.4 МэВ в интервале углов вылета α-частиц 25°–160° (лаб.) и сопоставлены с теоретическими. При расчетах в [1] волновые функции (ВФ) участвующих в реакции ядер выбирались в модели Нильссона [2], в качестве механизмов реакции рассматривались подхват тритона в рамках метода связи каналов (МСК) [3] (программный код CHUCK [4]) и в модель составного ядра (СЯ) (программный код CNDENSI [5]). Было показано, что в рамках выбранных моделей УР более или менее корректно описываются только до углов вылета α-частиц θα < 120° с образованием 24Mg и в основном состоянии, и в состоянии 2+. При бóльших значениях углов θα расчеты не описывают экспериментальные сечения, особенно максимум при θα ≈ 150°.

В настоящей работе с целью описать экспериментальные УР α-частиц при больших углах их вылета рассмотрен вклад в сечение реакции 27Al(p, α)24Mg(0+, 2+) механизма срыва [6] тяжелого кластера 23Na с начального ядра 27Al(5/2+). Расчет сечения этого механизма проведен по программе FRESCO [7] с учетом вклада основного 3/2+ состояния 23Na [8]. Спектроскопические амплитуды (СА) в вершинах распада 27Al(5/2+) → → 23Na(3/2+) + α и 24Mg(0+, 2+) → 23Na(3/2+) + р, как и в [1], рассчитывались с использованием ВФ модели Нильссона [2] с учетом их зависимости от деформации.

ВОЛНОВЫЕ ФУНКЦИИ ЯДЕР 27Al(5/2+), 24Mg(0+, 2+) и 23Na(3/2+, 5/2+) В МОДЕЛИ НИЛЬССОНА

ВФ ядра $\left| {J\,\Omega } \right\rangle $ с полным моментом J и проекцией Ω на ось его симметрии в модели Нильссона [2] может быть представлена в виде суперпозиции конфигураций ВФ нильссоновских орбиталей каждая из которых может быть разложена по ортонормированному базису $\left| {2{{l}_{i}}\,{{\mu }_{i}}\,{{\sigma }_{i}}} \right\rangle $ одночастичных собственных функций (заполнение орбиталей нуклонами начинается с их низших собственных энергий)

(1)

В (1) ${{l}_{i}}({{\mu }_{i}})$ – значение орбитального момента нуклона (его проекции), ${{\sigma }_{i}}$ – проекции его спина, Ωi = μi + σi; $\sum\nolimits_i {{{\Omega }_{i}}} = \Omega .$ Четыре нуклона с различными проекциями спинов и изоспинов образуют замкнутую конфигурацию орбиталей с Ω = 0.

Ядра 27Al(JA = 5/2+), 24Mg(JB = 0+, 2+) и 23Na(JС = = 3/2+) принадлежат 1d–2s-оболочке. ВФ в модели Нильссона 27Al(JA = 5/2+), 24Mg(JB = 0+, 2+) с положительной квадрупольной деформацией (β2 = = +0.25 для 27Al и +0.4 для 24Mg) [9]) имеют вид [2]:

(2)
(3)
(4)

Основное состояние 23Na(JС = 3/2+) в модели Нильссона имеет протонную дырку в конфигурации (№7)4, так что его ВФ определяется выражением

(5)

КАЧЕСТВЕННЫЕ ОЦЕНКИ ВЛИЯНИЯ ДЕФОРМАЦИИ ЯДЕР НА СПЕКТРОСКОПИЧЕСКИЕ АМПЛИТУДЫ МЕХАНИЗМА СРЫВА ТЯЖЕЛОГО КЛАСТЕРА В РЕАКЦИИ 27Al(р, α)24Mg(0+, 2+)

Механизм срыва тяжелого кластера 23Na в реакции 27Al(р, α)24Mg иллюстрируется полюсной диаграммой рис. 1. Вероятность распада в каждой из вершин этой диаграммы характеризуется СА.

Рис. 1.

Диаграмма, иллюстрирующая механизм срыва тяжелого кластера 23Na(3/2+) в реакции 27Al(р, α)24Mg.

Основная сложность заключается в расчете СА в вершине 27Al → 23Na + α. Согласно (3), в основном состоянии 27Al содержит три $1{{d}_{{{5 \mathord{\left/ {\vphantom {5 2}} \right. \kern-0em} 2}}}}$ нуклона (один протон и два нейтрона в орбитали (№ 5)), ВФ которых не зависит от деформации. Чтобы сформировать α-частицу, из 27Al необходимо отделить еще один $1{{d}_{{3/2}}}$ протон из конфигурации (№ 7)4 с компонентой его ВФ $\left| {222 - } \right\rangle $, при этом весовой коэффициент зависит от деформации 27Al и при β2 = +0.25 меньше единицы [2].

Вероятность распада 24Mg(0+, 2+) → 23Na(3/2+) + + р также определяется компонентой $\left| {222 - } \right\rangle $ в 23Na(3/2+) и 24Mg(0+, 2+). Из (3–5) следует, что ВФ 24Mg и 23Na во всех рассматриваемых состояниях содержат компоненту $\left| {222 - } \right\rangle $. Ее вклад в ВФ (3)–(5) при большой деформации 24Mg и 23Na относительно мал, что уменьшает величину СА в вершине распада 24Mg → 23Na + р.

Таким образом, учет деформации начального, конечного ядер и передаваемого тяжелого кластера приводит к уменьшению СА в вершинах рассматриваемой диаграммы. В следующем разделе проведены расчеты конкретных СА для механизма срыва тяжелого кластера в реакции 27Al(р, α)24Mg(0+, 2+).

РАСЧЕТ СА ДЛЯ МЕХАНИЗМА СРЫВА ТЯЖЕЛОГО КЛАСТЕРА 23Na(3/2+) В РЕАКЦИИ 27Al(р, α)24Mg(${{J}_{B}}$ = 0)

Согласно [6] СА распада в вершине 1 – $B \to C + x$ – определяется выражением

(6)
$\Theta _{{{{\Lambda }_{1}}{{J}_{C}}}}^{{B \to C + x}} = \sum {\left\langle {{{{L}_{B}}{{S}_{B}}}} \mathrel{\left | {\vphantom {{{{L}_{B}}{{S}_{B}}} {{{J}_{B}}}}} \right. \kern-0em} {{{{J}_{B}}}} \right\rangle \left\langle {{{{L}_{C}}{{S}_{C}}}} \mathrel{\left | {\vphantom {{{{L}_{C}}{{S}_{C}}} {{{J}_{C}}}}} \right. \kern-0em} {{{{J}_{C}}}} \right\rangle } \left\langle {{{{L}_{x}}{{S}_{x}}}} \mathrel{\left | {\vphantom {{{{L}_{x}}{{S}_{x}}} {{{J}_{x}}}}} \right. \kern-0em} {{{{J}_{x}}}} \right\rangle \left\langle {{{{L}_{C}}{{L}_{1}}}} \mathrel{\left | {\vphantom {{{{L}_{C}}{{L}_{1}}} {{{L}_{B}}}}} \right. \kern-0em} {{{{L}_{B}}}} \right\rangle \left\langle {{{{S}_{C}}{{S}_{x}}}} \mathrel{\left | {\vphantom {{{{S}_{C}}{{S}_{x}}} {{{S}_{B}}}}} \right. \kern-0em} {{{{S}_{B}}}} \right\rangle \left\langle {{{{\Lambda }_{1}}{{L}_{x}}}} \mathrel{\left | {\vphantom {{{{\Lambda }_{1}}{{L}_{x}}} {{{L}_{1}}}}} \right. \kern-0em} {{{{L}_{1}}}} \right\rangle \tilde {\Theta }_{{{{\Lambda }_{1}}{{L}_{B}}{{L}_{C}}{{L}_{1}}{{L}_{x}}{{S}_{B}}{{S}_{C}}{{S}_{x}}}}^{{B \to C + x}},$
где суммирование ведется по всем моментам (их проекциям), разрешенным правилами отбора. При использовании кода FRESCO в случае передачи тяжелого кластера С в СА (6) необходимо выделить связь моментов в виде ❬Λ1JC | I1❭❬JxI1 | JB❭, где ${{\Lambda }_{1}}$ – орбитальный момент функции относительного движения кластеров С и х.

Выполняя соответствующие алгебраические преобразования в (6), получим [6]:

(7)
$\Theta _{{\Lambda {{J}_{C}}}}^{{B \to C + x}} = \left\langle {{{{\Lambda }_{1}}{{J}_{C}}}} \mathrel{\left | {\vphantom {{{{\Lambda }_{1}}{{J}_{C}}} {{{I}_{1}}}}} \right. \kern-0em} {{{{I}_{1}}}} \right\rangle \left\langle {{{{J}_{x}}{{I}_{1}}}} \mathrel{\left | {\vphantom {{{{J}_{x}}{{I}_{1}}} {{{J}_{B}}}}} \right. \kern-0em} {{{{J}_{B}}}} \right\rangle W({{\Lambda }_{1}}{{I}_{1}}),$
где

(8)
$\begin{gathered} W({{\Lambda }_{1}}{{I}_{1}}) = \tilde {\Theta }_{{{{\Lambda }_{1}}{{L}_{B}}{{L}_{C}}{{L}_{1}}{{L}_{x}}{{S}_{B}}{{S}_{C}}{{S}_{x}}}}^{{B \to C + x}}{{( - 1)}^{{{{S}_{C}} + {{S}_{x}} - {{S}_{B}}}}}\sqrt {\frac{{LL{{S}_{C}}{{J}_{B}}{{J}_{x}}{{J}_{C}}}}{{{{L}_{C}}{{I}_{1}}}}} \times \\ \times \,\,\sum\limits_{{{I}_{1}}{{M}_{{{{I}_{1}}}}}} {\left\{ {\begin{array}{*{20}{c}} {{{L}_{x}}}&{{{S}_{x}}}&{{{J}_{x}}} \\ {{{L}_{1}}}&{{{S}_{C}}}&{{{I}_{1}}} \\ {{{L}_{B}}}&{{{S}_{B}}}&{{{J}_{B}}} \end{array}} \right\}} \cdot U({{L}_{x}}{{\Lambda }_{1}}{{L}_{B}}{{L}_{C}}:{{L}_{1}}L) \cdot U({{\Lambda }_{1}}{{L}_{1}}{{J}_{C}}{{S}_{C}}:{{L}_{C}}{{I}_{1}}). \\ \end{gathered} $

В (6) приведенная СА $\tilde {\Theta }_{{{{\Lambda }_{1}}{{L}_{B}}{{L}_{C}}{{L}_{1}}{{L}_{x}}{{S}_{B}}{{S}_{C}}{{S}_{x}}}}^{{B \to C + x}}$ определяется выражением

$\tilde {\Theta }_{{{{\Lambda }_{1}}{{L}_{B}}{{L}_{C}}{{L}_{1}}{{L}_{x}}{{S}_{B}}{{S}_{C}}{{S}_{x}}}}^{{B \to C + x}} = {{\left( {\frac{B}{C}} \right)}^{{{{{{N}_{x}}} \mathord{\left/ {\vphantom {{{{N}_{x}}} 2}} \right. \kern-0em} 2}}}}{{\Im }_{{1{{J}_{B}}}}}\left\langle {{{T}_{x}}{{T}_{C}}} \right|\left. {{{T}_{B}}} \right\rangle ,$
где ${{N}_{x}}$ – число квантов, уносимых отделяемыми нуклонами, ${{\left( {\frac{B}{C}} \right)}^{{{{{{N}_{x}}} \mathord{\left/ {\vphantom {{{{N}_{x}}} 2}} \right. \kern-0em} 2}}}}$ – множитель отдачи, ${{\Im }_{{1{{J}_{B}}}}}$ – интеграл перекрытия ВФ ядер В (в состоянии JB) и С:

$\begin{gathered} {{\Im }_{{1{{J}_{B}}}}} = \left\langle {{{\Psi }_{B}}({{J}_{B}})} \right.\left| {{{\Psi }_{C}}} \right\rangle = \\ = {{K}_{x}}(l,\Lambda )Г К \left\langle {{{\varphi }_{x}}} \right|{{\Psi }_{{{{N}_{x}}\,\Lambda }}}({{{\vec {r}}}_{x}} - {{{\vec {R}}}_{C}}), \\ \end{gathered} $

для расчета которого в LS-связи мы использовали нильссоновскую модель, $\left\langle {{{T}_{x}}\,{{T}_{C}}\,} \right|\left. {{{T}_{B}}} \right\rangle $ – изоспиновый коэффициент Клебша–Гордана, ГК – генеалогический коэффициент для ядер 1d-оболочки [10].

Аналогичным образом определяется СА распада в вершине 2 – $A \to C + y.$

10 для вершины 24Mg(0+) = 23Na(3/2+) + р

Определим интеграл перекрытия ВФ ядер 24Mg(0+) (В) и 23Na(3/2+) (С) в модели Нильссона. Для этого преобразуем ВФ (3) основного состояния ядра 24Mg с параметром квадрупольной деформации β2 = +0.4.

(9)
где $\left\langle {{{d}^{4}}\,[4]\,{}^{{11}}S} \right|\left. {{{d}^{3}}\,[4]\,{}^{{22}}D,\,\,d} \right\rangle $ – одночастичный ГК [10], равный 1. Компонента $\left| {222 - } \right\rangle $ ВФ ядра 24Mg имеет вес = 0.174 при β2 = +0.4. Подставляя в (9) эти величины для ℑ10, получаем

(10)
${{\Im }_{{10}}} = 0.348.$

12 для вершины 24Mg(2+) = 23Na(3/2+)+

Для расчета интеграла перекрытия ℑ12 в модели Нильсона ВФ (4) ядра 24Mg(2+) (В) с параметром квадрупольной деформации β2 = +0.4 представим в виде:

(11)

В (11) = 0.985, = 0.174, =  –0.996, а $\left| {221{\kern 1pt} - } \right\rangle $ = –$\left| {221 - } \right\rangle {\text{*}}.$ ГК $\left. {\left\langle {{{d}^{3}}\,[3]{}^{{22}}D} \right|{{d}^{3}}\,[2]{{\,}^{{13}}}S,d} \right\rangle \,$ = $\frac{{\sqrt 7 }}{{\sqrt {15} }}.$ В результате (11) принимает вид

(12)

Аналогичным образом преобразуем ВФ 23Na(3/2+):

(13)

Объединяя (12, 13), для интеграла перекрытия получаем

(14)
${{\Im }_{{12}}} = \frac{{3\sqrt 2 }}{{\sqrt 5 }}0.174 \cdot 0.981.$

2 для вершины 27Al(5/2+) = 23Na(3/2+) + α

Интеграл перекрывания ВФ ядер 27Al(5/2+) (β2 = = +0.25) и 23Na(3/2+) (β2 = +0.4) равен

(15)
где = 0.252 [2] при β2(27Al) = +0.25; ${{K}_{t}}({{d}^{3}}\,[3]\,{{l}_{2}} = 2)$ – обобщенные коэффициенты Тальми [11]. Их конкретные значения, выделяющие в ВФ трех d-нуклонов внутреннюю ВФ тритона и ВФ относительного движения трития и ядра 24Mg, приведены в [1]. $\left\langle {6222:{{\Lambda }_{2}}} \right|\left. {3,1} \right|\left. {8{{\Lambda }_{2}}00:{{\Lambda }_{2}}} \right\rangle $ – коэффициенты Тальми с разными массами [11], выделяющие в системе четырех d-нуклонов внутреннюю ВФ α-частицы и ВФ относительного движения α-частицы и ядра 23Na. Произведения ${{K}_{t}}({{d}^{3}}[3]\,\,{{l}_{2}} = 2)$·$\left\langle {6222:{{\Lambda }_{2}}} \right|\left. {3,\,\,1} \right|\left. {8{{\Lambda }_{2}}00:{{\Lambda }_{2}}} \right\rangle $ равны $\frac{{\sqrt {33} }}{{32\sqrt {14} }}$ = 0.058 для ${{\Lambda }_{2}}$ = 2 и $\frac{{\sqrt {11 \cdot 13} }}{{64\sqrt 7 }}$ = 0.071 для ${{\Lambda }_{2}}$ = 4. В результате

(16)
${{\Im }_{2}} = 0.504\sqrt {165} \left\{ \begin{gathered} 0.06\,\,\,\,{{\Lambda }_{2}} = 2 \hfill \\ 0.07\,\,\,\,{{\Lambda }_{2}} = 4. \hfill \\ \end{gathered} \right.$

Все рассчитанные нами СА для левой и правой вершин диаграммы рис. 1 приведены в табл. 1.

Таблица 1.  

Спектроскопические амплитуды для механизма срыва тяжелого кластера 23Na(3/2+) в реакции 27Al(р, α)24Mg

  Λi Ii $W({{\Lambda }_{i}}{{I}_{i}})$
24Mg(0+) → 23Na(3/2+) + p 2 1/2 –0.363
24Mg(2+) → 23Na(3/2+) + p 0 3/2 0.113
2 3/2 0.201
2 5/2 0.113
27Al(5/2+) → 23Na(3/2+) + α 2 5/2 0.33
4 5/2 0.943

ВКЛАД МЕХАНИЗМА СРЫВА ТЯЖЕЛОГО КЛАСТЕРА 23Na(3/2+) В УГЛОВОЕ РАСПРЕДЕЛЕНИЕ α-ЧАСТИЦ – ПРОДУКТОВ РЕАКЦИИ 27Al(р, α)24Mg

Расчет сечения механизма передачи тяжелого кластера 23Na(3/2+), иллюстрируемого диаграммой на рис. 1, проведен по программе FRESCO [7]. СА в вершинах распада ядер брались из табл. 1. Значения параметры оптических потенциалов во входном и выходном каналах реакции взяты такими же, как и в [1], и приведены в табл. 2.

Таблица 2.  

Параметры оптических потенциалов во входном и выходном каналах реакции 27Al(р, α)24Mg

  V, МэВ rV, фм aV, фм W, МэВ rWD, фм аWD, фм Vso, МэВ rso, фм aso, фм rC, фм Ссылка
p+27Al(5/2+) 55.317 1.169 0.674 7.62 1.295 0.533 5.644 0.97 0.59 1.33 [12]
α+24Mg(0+,2+) 194. 1.4 0.65 7.25 1.4 0.65       1.4 [13]

Результаты расчетов вместе с экспериментальным сечением, полученным в [1], показаны на рис. 2.

Рис. 2.

Рассчитанные дифференциальные сечения механизма срыва тяжелого кластера 23Na(3/2+) в реакции 27Al(р, α)24Mg (сплошные кривые) с образованием 24Mg в основном (а) и возбужденном $2_{1}^{ + }$ (б) состояниях. Экспериментальные данные (черные точки) взяты из [1].

Как видно из рисунка, расчетное угловое распределение α-частиц на два порядка меньше экспериментального для обоих состояний 24Mg. Такое различие экспериментальных и теоретических сечений можно объяснить малыми значениями СА в обеих вершинах распада начального и конечного ядер, обусловленных тем, что коэффициенты разложения ВФ нильссоновской орбитали существенно меньше единицы.

Более того, форма рассчитанного углового распределения не характерна для механизма срыва тяжелого кластера: за счет большого переданного импульса сечение этого механизма должно иметь осцилляции в задней полусфере и подъем при углах θα, близких к 180°. При образовании 24Mg в основном состоянии рассчитанное сечение имеет практически симметричную форму угловой зависимости и узкий максимум при θα ~ 180°, в состоянии $2_{1}^{ + }$ – максимум при θα > 140°, но их величина на два порядка меньше экспериментального. Другими словами, механизм срыва тяжелого кластера не дает вклад в дифференциальное сечение реакции 27Al(р, α)24Mg(0+, 2+).

ЗАКЛЮЧЕНИЕ

В настоящей работе рассмотрен вклад механизма срыва тяжелого кластера 23Na в сечение реакции 27Al(р, α)24Mg(0+, 2+).

Расчет сечения этого механизма проведен по программе FRESCO [7], СА в вершинах распада 24Mg → 23Na + р и 27Al → 23Na + α вычислялись в модели Нильссона с учетом деформации всех ядер.

Рассчитанные сечения для уровней 0+ и 2+ ядра 24Mg меньше экспериментальных на два порядка, сечение с образованием уровня 0+ имеет узкий подъем к углу θα ∼ 140°. Для уровня 2+ в расчетном сечении имеется максимум при θα > 140°, но его величина также меньше экспериментального на два порядка. Такое поведение вклада механизма срыва тяжелого кластера можно, по крайней мере, качественно объяснить аномально малыми СА в вершинах распада начального и конечного ядер вследствие слабого перекрывания ВФ ядер 27Al, 23Na и 24Mg в модели Нильссона. Малые значения СА приводят к тому, что вклад механизма срыва тяжелого кластера 23Na в угловое распределение протонов в реакции 27Al(p, α)24Mg занижен по сравнению с экспериментом на два порядка. Это означает, что рассматриваемый механизм не играет существенной роли в реакциях передачи на ядрах 1d-оболочки.

Список литературы

  1. Галанина Л.И., Зеленская Н.С., Лебедев В.М. и др. // Изв. РАН. Сер. физ. 2016. Т. 80. С. 338; Galanina L.I., Zelenskaya N.S., Lebedev V.M. et al. // Bull. Russ. Acad. Sci. Phys. 2016. V. 80. P. 304.

  2. Нильссон С.Г. Связанные состояния индивидуальных нуклонов в сильно деформированных ядрах. В книге “Деформация атомных ядер”. М.: Иностр. лит-ра, 1958. С. 232.

  3. Tamura T. // Rev. Mod. Phys. 1965. V. 37. P. 679.

  4. Kunz P.D., Rost E. // Comp. Nucl. Phys. 1993. V. 2. P. 88.

  5. Belyaeva T.L., Zelenskaya N.S., Odintzov N.V. // Comp. Phys. Comm. 1992. V. 73. P. 161.

  6. Зеленская Н.С., Теплов И.Б. Обменные процессы в ядерных реакциях. М.: Изд-во МГУ, 1985. 167 с.

  7. Tompson I.J. // Comput. Phys. Rep. 1988. V. 7. P. 167; http:// www.fresco.org.uk/.

  8. Endt P.M., Van Der Leun C. // Nucl. Phys. A 1978. V. 310. P. 94.

  9. http://cdfe.sinp.msu.ru/services/radchart/radmain.html/.

  10. Jahn H.A., Van Wieringer H. // Pros. Roy. Sos. 1951. V. A209. P. 502.

  11. Неудачин В.Г., Смирнов Ю.Ф. Нуклонные ассоциации в легких ядрах. М.: Изд-во “Наука”, 1969. 414 с.

  12. Koning A.J., Delaroche J.P. // Nucl. Phys. A. 2003. V. 713. P. 231.

  13. Thompson W.J., Grawford G.E., Davis R.H. // Nucl. Phys. A. 1967. V. 98. P. 228.

Дополнительные материалы отсутствуют.