Поверхность. Рентгеновские, синхротронные и нейтронные исследования, 2023, № 1, стр. 20-24

Тормозная способность вещества для пучка моноэнергетических альфа-частиц

Н. Н. Михеев a*, И. Ж. Безбах a

a Институт кристаллографии им. А.В. Шубникова ФНИЦ “Кристаллография и фотоника” РАН
119333 Москва, Россия

* E-mail: kmikran@spark-mail.ru

Поступила в редакцию 25.06.2022
После доработки 14.07.2022
Принята к публикации 14.07.2022

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

Аннотация

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

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

ВВЕДЕНИЕ

Тормозная способность вещества для направленных пучков моноэнергетических заряженных частиц представляет собой исключительно важную характеристику сложных процессов взаимодействия частиц с веществом. Она необходима во многих областях фундаментальной и прикладной физики [13]. Например, для количественных методов локальной диагностики современных наноматериалов и структур (сканирующей гелиевой ионной микроскопии [4, 5]) необходимы надежные знания о тормозной способности исследуемых образцов по отношению к альфа-частицам с энергией E0 от 100 эВ до 30 кэВ. Для пучка моноэнергетических протонов задача аналитического описания зависимости тормозной способности вещества от энергии в широком диапазоне первичной энергии частиц была решена в [6]. Получены формулы, учитывающие влияние зависимости энергетических потерь частиц от их скорости в четырех различных энергетических диапазонах. Особенность транспорта протонов в твердотельном материале состоит в том, что заряд иона водорода Z1 остается неизменным и равным единице электрического заряда вплоть до энергии частиц E0 несколько кэВ. Ситуация существенно иная для ионов с большим электрическим зарядом начиная с альфа-частиц. Общепринято [7], что при скорости V0 пучка ионов меньше, чем $5Z_{1}^{{{2 \mathord{\left/ {\vphantom {2 3}} \right. \kern-0em} 3}}} \times 2.178 \times {{10}^{8}}\,\,\frac{{{\text{см}}}}{{\text{с}}},$ необходимо уже учитывать зависимость зарядового состояния Z1ef ионов от их скорости. Первоначальная оценка среднего заряда иона, двигающегося со скоростью V0, была сделана Бором [8], который определил средний заряд иона Z1ef как:

$\begin{gathered} \,\,\,\,\,\,{{Z}_{{1{\text{ef}}}}} \approx Z_{1}^{{{1 \mathord{\left/ {\vphantom {1 3}} \right. \kern-0em} 3}}}({{{{V}_{0}}} \mathord{\left/ {\vphantom {{{{V}_{0}}} {{{V}_{{\text{B}}}}}}} \right. \kern-0em} {{{V}_{{\text{B}}}}}}), \\ {\text{где}}\,\,{{V}_{{\text{B}}}} = 2.18 \times {{10}^{8}}\,\,{{{\text{см}}} \mathord{\left/ {\vphantom {{{\text{см}}} {\text{с}}}} \right. \kern-0em} {\text{с}}},\,\,{\text{а}}\,\,{{V}_{{\text{B}}}} < {{V}_{0}} < {{Z}^{{{2 \mathord{\left/ {\vphantom {2 3}} \right. \kern-0em} 3}}}}{{V}_{{\text{B}}}}. \\ \end{gathered} $
А в настоящее время наиболее часто используется выражение для Z1ef [9]:
${{Z}_{{{\text{1ef}}}}} = {{Z}_{1}}\left[ {1 - C\exp \left( { - Z_{1}^{\gamma }\frac{{{{V}_{0}}}}{{{{V}_{{\text{B}}}}}}} \right)} \right],$
где C ≈ 1, γ ≈ 2/3 для большинства материалов, изученных в [9]. Детальный обзор существующих аналитических выражений, используемых для описания зарядового состояния низкоэнергетических ионов, проведен в [7].

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

Цель настоящей работы заключалась в разработке новой методики расчета зарядового состояния Z1ef пучка альфа-частиц в тонких слоях вещества на основе изучения зависимости Z1ef от соотношения скоростей V0 частиц пучка и средней скорости V1i самых медленных электронов вещества образца.

МЕТОДИКА РАСЧЕТА Z1ef ДЛЯ АЛЬФА-ЧАСТИЦ

Базовая формула для Z1ef

Как показано в [6, 10], для пучка заряженных частиц массой M с энергией меньше E01 = = (M/me)J 2/8ε1i в тонких пленках наиболее вероятная потеря энергии ε1i соответствует энергии связи самых медленных электронов образца, средняя скорость которых V1i (J – усредненное значение потенциальной энергии электронов мишени). Поэтому в статистической модели многократного рассеяния этот параметр является ключевым при описании энергетических потерь частиц с энергией меньше E01.

Ранее в [11] учет уменьшения вероятности неупругих потерь пучка заряженных частиц на электронах мишени, когда скорость пучка становится сопоставимой со средней скоростью VZ/2 электронов области экранирования, был осуществлен путем введения функции $F_{M}^{N}{\text{:}}$

$F_{M}^{N} = 2.143\left[ {1 - \exp \left\{ { - 0.1714{{{\left( {\frac{{{{m}_{e}}V_{0}^{2}}}{{{{\varepsilon }_{{{Z \mathord{\left/ {\vphantom {Z 2}} \right. \kern-0em} 2}}}}}}} \right)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}} \right\}} \right],$
где εZ/2 – энергия электронов атома, соответствующих граничной области экранирования заряда ядра в модели Томаса–Ферми [12]. Она успешно используется в прикладных задачах [1316], связанных с описанием энергетических потерь пучков электронов и протонов средних и низких энергий в веществе. Функцию $F_{M}^{N}$ можно представить в виде:
$F_{M}^{N} = 2.314\left[ {1 - \exp \left\{ { - 0.2424\frac{{{{V}_{0}}}}{{{{V}_{{{Z \mathord{\left/ {\vphantom {Z 2}} \right. \kern-0em} 2}}}}}}} \right\}} \right].$
По существу, эта формула отражает зависимость вероятности неупругого рассеяния пучка частиц на электронах, расположенных вне области экранирования, когда V0/VZ/2 < 4.7. И поэтому вполне разумно, используя метод аналогий, предположить, что уменьшение заряда альфа-частиц с Z1 = 2 до заряда Z1ef может быть обусловлено двумя возможными факторами. Во-первых, уменьшением вероятности однократной ионизации атомов вещества, вызванным снижением отношения V0/V1i – оно становится меньше некоторой величины 4.7–5.0. И, во-вторых, изменением заряда некоторой части альфа-частиц пучка после захвата ими слабо связанных атомных электронов мишени. Поэтому для Z1ef возможно записать базовую формулу в виде:
(1)
${{Z}_{{1{\text{ef}}}}} = {{Z}_{1}}\left[ {1 - \exp \left( { - {\kern 1pt} \beta \frac{{{{V}_{0}}}}{{{{V}_{{1i}}}}}} \right)} \right],$
где β – некоторая функция переменной V0/V1i, опосредовано связанная с вероятностью присоединения электрона вещества к иону гелия в зависимости от соотношения скоростей иона V0 и слабо связанного электрона V1i.

Определение функции $\beta $ для альфа-частиц

Для изучения зависимости функции β от отношения V0/V1i в широком диапазоне от 1.5 до 30 были проведены расчеты тормозной способности алюминия как модельного материала для альфа-частиц с Z1 = 2 и энергией E0 в диапазоне от 10 кэВ до 10 МэВ. Расчеты проводили в рамках концепции дискретного процесса торможения заряженных частиц [6], используя значения параметров этого материала ε1i = 6 эВ и εZ/2 = 80 эВ, надежно определенные по результатам торможения пучка протонов в алюминии. Далее полученные результаты сравнивали с экспериментальными значениями тормозной способности S, представленными в [17, 18]. По результатам проведенного анализа была предложена формула, устанавливающая функциональную зависимость величины β от переменной V0/V1i для моноэнергетических пучков альфа-частиц в алюминии:

(2)
$\begin{gathered} \beta = 0.805 + 0.25\exp \left[ { - \frac{{{{{\left( {4.25 - \frac{{{{V}_{0}}}}{{{{V}_{{1i}}}}}} \right)}}^{2}}}}{{1.2211}}} \right] \\ {\text{при}}\,\,1.5 \leqslant {{{{V}_{0}}} \mathord{\left/ {\vphantom {{{{V}_{0}}} {{{V}_{{1i}}}}}} \right. \kern-0em} {{{V}_{{1i}}}}} \leqslant 4.25, \\ \end{gathered} $
(3)
$\beta = 1.055\,\,{\text{при}}\,\,{{{{V}_{0}}} \mathord{\left/ {\vphantom {{{{V}_{0}}} {{{V}_{{1i}}}}}} \right. \kern-0em} {{{V}_{{1i}}}}} > 4.25.$
Полученные формулы (1)(3) определяют поведение зарядового состояния Z1ef пучка альфа-частиц с энергией E0, соответствующей интервалу скоростей частиц V1i 1.55.0. При V0 < 1.5V1i заряд Z1ef пучка альфа-частиц достигает своего равновесного значения 2√2/2 = √2.

РЕЗУЛЬТАТЫ ПРИМЕНЕНИЯ МЕТОДИКИ РАСЧЕТА Z1ef

Результаты применения полученной функциональной зависимости для β при расчете тормозной способности S алюминия для альфа-частиц с использованием формулы (1) представлены на рис. 1. Как видно из рис. 1а, учет Z1ef позволяет сделать следующие выводы: для альфа-частиц с энергией E0 более 1.0 МэВ заряд е Z1ef = 2 в алюминии остается неизменным; в интервале значений энергии E0 1.0–0.1 МэВ Z1ef в алюминии плавно уменьшается до равновесного значения √2; на графике зависимости тормозной способности S от энергии E0 стабилизация Z1ef (достижение равновесного значения) сопровождается появлением характерного перегиба.

Рис. 1.

Зависимость тормозной способности S алюминия для пучка альфа-частиц с энергией от 10 кэВ до 100 МэВ (He++13Al): а – расчет с Z1 = 2 по формулам [6] (пунктир) и Z1ef по формулам (1)(3) (сплошная кривая); б – расчет с Z1ef по формулам (1)(3) (сплошная кривая) и экспериментальные данные [17] (символы).

Результаты сравнения рассчитанных по методике значений S с результатами экспериментальных измерений [17, 18] представлены на рис. 1б и свидетельствуют о хорошем соответствии расчетов и экспериментальных данных для алюминия во всем приведенном энергетическом диапазоне E0.

Результаты использования новой методики учета зарядового состояния ионов гелия при расчете тормозной способности S углерода, серебра и золота для пучка альфа-частиц с использованием формулы (1) представлены на рис. 2–4. Сказанное выше о результатах, представленных на рис. 1, в полной мере относится и к результатам, приведенным на рис. 2–4, которые подтверждают значимость и возможности новой методики учета зарядового состояния альфа-частиц для образцов, относящихся к широкому диапазону значений среднего атомного номера Z материала мишени.

Рис. 2.

Зависимость тормозной способности S углерода для пучка альфа-частиц с энергией от 6 кэВ до 50 МэВ (He++6C): пунктирная линия – расчет с Z1 = 2 по формулам [6]; сплошная кривая – расчет с Z1ef по формулам (1)(3).

Рис. 3.

Зависимость тормозной способности S серебра для пучка альфа-частиц с энергией от 2 кэВ до 100 МэВ (He++47Ag): пунктирная линия – расчет с Z1 = 2 по формулам [6]; сплошная кривая – расчет с Z1ef по формулам (1)(3).

Рис. 4.

Зависимость тормозной способности S золота для пучка альфа-частиц с энергией от 10 кэВ до 100 МэВ (He++79Au): пунктирная линия – расчет с Z1 = 2 по формулам [6]; сплошная кривая – расчет с Z1ef по формулам (1)(3).

Отметим, что присутствие характерного перегиба на рассчитанных кривых S = S(E0) не является случайным результатом проведенных вычислений, а только подтверждает и объясняет особенность экспериментальных измерений тормозной способности. Эта особенность, однако, слабо выявляется из-за значительного разброса экспериментально измеренных величин S в этом диапазоне значений энергии E0.

ЗАКЛЮЧЕНИЕ

Разработана новая методика учета зарядового состояния Z1ef альфа-частиц. Предложены формулы для расчета, устанавливающие функциональную зависимость Z1ef от отношения V0/V1i.

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

Установлено хорошее соответствие выполненных расчетов тормозной способности S экспериментальным данным во всем исследованном энергетическом диапазоне E0. Показано, что до пятикратного превышения V0 над V1i величина Z1ef = Z1 = 2. При уменьшении отношения V0/V1i до 1.5 происходит плавное изменение заряда ионов гелия Z1ef до равновесного значения √2.

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

  1. ICRU Report 49. Stopping Powers and Ranges for Protons and Alpha Particles. International Commission on Radiation Units and Measurements. 1993.

  2. Andersen H.H., Ziegler J.F. Hydrogen Stopping Powers and Ranges in All Elements. N.Y.: Pergamon Press, 1977.

  3. ICRU Report 37. Stopping Powers for Electrons and Positrons. International Commission on Radiation Units and Measurements. 1984.

  4. Ward E.W., Notte J., Economou N.P. // J. Vac. Sci. Technol. B. 2006. V. 24. P. 2871.

  5. Van Gastel R., Hlawacek G., Zandvliet H.J.W., Poelsema B. // Microelectron. Reliability. 2012. V. 52. № 9–10. P. 2104. https://doi.org/10.1016/j. microrel. 2012.06.130

  6. Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2022. № 3. С. 94. https://doi.org/10.31857/S1028096022030141

  7. Cruz S.A. // Rad. Effects. 1986. V. 88. P. 159.

  8. Bohr N. // K. Dan Vidensk. Mat. Fys. Medd. 1948. V. 18. P. 8.

  9. Betz H.D., Hortig G, Leischner E., Schmelzer Ch., Stadler B., Weihrauch J. // Phys. Lett. A. 1966. V. 22. P. 643.

  10. Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2022. № 6. С. 54. https://doi.org/10.31857/S102809602

  11. Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2010. № 4. С. 25.

  12. Ландау Л.Д., Лившиц Е.М. Квантовая механика. Нерелятивистская теория. М.: Наука, 1974. 607 с.

  13. Михеев Н.Н., Колесник А.С. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2017. № 12. С. 88. https://doi.org/10.7868/S0207352817120083

  14. Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2020. № 12. С. 70. https://doi.org/10.31857/ S1028096020120120201

  15. Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2020. № 3. С. 77. https://doi.org/10.31857/S1028096020030127

  16. Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2019. № 8. С. 56. https://doi.org/10.1134/S0207352819080109

  17. Paul H. IAEA. NDS. https://www-nds.iaca.org/stopping

  18. Ziegler J.F. Handbook Stopping Powers and Ranges in All Elemental Matters. N.Y.: Pergamon Press, 1978. 368 p.

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