Поверхность. Рентгеновские, синхротронные и нейтронные исследования, 2020, № 12, стр. 70-77
Двухпотоковая модель обратного рассеяния моноэнергетического пучка электронов: приложение к задачам диагностики тонкопленочных структур
Н. Н. Михеев *
ФНИЦ “Кристаллография и фотоника” Российской академии наук
248640 Калуга, Россия
* E-mail: kmikran@spark-mail.ru
Поступила в редакцию 20.01.2020
После доработки 25.03.2020
Принята к публикации 28.03.2020
Аннотация
Применительно к задачам метода спектроскопии обратно рассеянных электронов, выполнены расчеты энергетических спектров электронов, отраженных от отдельных слоев, расположенных внутри многослойной тонкопленочной мишени. Уточнена структура аналитических выражений, описывающих распределение коэффициента обратного рассеяния электронов по глубине образца. Получена формула для двумерного распределения обратно рассеянных электронов по полярному углу после их выхода из пленочной мишени заданной толщины. Разработан алгоритм расчетов, необходимых для построения энергетических распределений обратно рассеянных электронов при заданном угловом положении входной диафрагмы спектрометра анализатора. Это позволяет приблизить выполняемые расчеты к реальным условиям регистрации обратно рассеянных электронов. Приводятся результаты всесторонней проверки полученных соотношений на основе сопоставления проведенных модельных расчетов энергетических распределений обратно рассеянных электронов с экспериментальными спектрами от свободных пленок, от массивных образцов и от многослойного тонкопленочного образца на массивной подложке.
ВВЕДЕНИЕ
Спектроскопия обратно рассеянных электронов (ОРЭ) моноэнергетического пучка РЭМ является перспективным электронно-зондовым методом для послойной, бесконтактной и неразрушающей диагностики многослойных тонкопленочных структур микро- и нано-электроники с высоким пространственным разрешением по глубине [1–3]. Практическая реализация такой электронно-зондовой томографии обусловлена возможностями выбора соответствующей глубины зондирования исследуемого образца (за счет варьирования энергии первичного пучка электронов) и использования режимов фильтрации и дискриминации регистрируемого энергетического спектра ОРЭ. Однако, как правило, общий вид энергетического спектра весьма слабо отражает внутреннюю структуру исследуемого образца. И определить из него, в каком диапазоне энергий следует расположить энергетическое “окно” анализатора для наилучшего выявления того или иного слоя, затруднительно. Поэтому для успешного и эффективного применения этого метода требуется детальная информация о вкладе в общий энергетический спектр структуры каждого отдельного слоя образца. А это предполагает умение решать прямую задачу по расчету сигнала ОРЭ от многослойной структуры с известной геометрией слоев и знанием их элементного состава. Полуэмпирическая модель расчета энергетических распределений ОРЭ, которая основывалась на учете вклада в общий энергетический спектр каждого отдельного тонкого слоя известного состава, расположенного внутри образца, была предложена ранее в работах [4, 5]. Однако вопросы соотношения двух групп первичных электронов: а) испытавших однократное рассеяние на угол θ > π/2; и б) изменивших направление своего движения в результате только малоуглового многократного рассеяния, – тогда окончательно решены не были. Эта проблема была решена только в работе [6], в результате использования статистической теории многократного рассеяния заряженных частиц в веществе [7–9].
Целью данной работы являлось решение прямой задачи аналитического описания пространственного и энергетического распределения обратно рассеянных электронов от многослойной тонкопленочной мишени на основе использования статистической теории многократного рассеяния заряженных частиц в веществе. В ходе ее реализации решалась задача по учету влияния углового распределения ОРЭ на их энергетические спектры. Получена формула для двумерного распределения ОРЭ по углу после их выхода из пленочной мишени заданной толщины. Разработан алгоритм расчетов, необходимых для построения энергетических распределений ОРЭ при заданном значении угла их регистрации. Это позволяет приблизить выполняемые расчеты к реальным условиям регистрации ОРЭ, то есть при фиксированном угловом положении входной диафрагмы спектрометра анализатора относительно поверхности образца. Приводятся результаты всесторонней проверки полученных соотношений на основе сопоставления проведенных модельных расчетов энергетических распределений ОРЭ с экспериментальными спектрами для свободных пленок, для массивных образцов и для образцов типа “пленочная мишень–подложка”.
БАЗОВЫЕ РЕЗУЛЬТАТЫ ПО ОРЭ
Известно, что ОРЭ представляют собой совокупность двух групп первичных электронов [11]: а) испытавших однократное рассеяние на угол θ > π/2 (первая группа с коэффициентом отражения η1) и б) изменивших направление своего движения в результате только малоуглового многократного рассеяния (вторая группа с коэффициентом отражения η2). Вопросы количественного соотношения вкладов этих групп в коэффициент η обратного рассеяния, как для пленок, так и для массивных мишеней представляет теоретический и практический интерес и в целом был решен в работе [6]. Коэффициент обратного рассеяния η первичных электронов пучка в двух потоковой модели определяется как η = η1 + η2, где коэффициент обратного рассеяния второй группы ОРЭ η2 = 0.24α ×η [6].
Здесь:
Величина параметра α определяется отношением максимального пробега к транспортной длине первичных электронов, и как показано в работе [14], зависит от величины среднего атомного номера Z мишени следующим образом. Для образцов с Z ≤ 16 величина α ≈ 0, так как расстояние от поверхности, на котором для электронов второй группы при многократном и малоугловом рассеянии все направления становятся равновероятными, располагается слишком глубоко. Так, что на возвратном пути в сторону поверхности они почти полностью теряют энергию, то есть термализуются. Например, для электронов пучка с энергией 20 кэВ в бериллии Ltr = 5.0 мкм, а Re = = 7.0 мкм, соответственно, α = 0.008. В алюминии Ltr = 1.7 мкм, а Re = 3.03 мкм, и α = 0.048. На интервале значений Z: 16 < Z ≤ 29, параметр α возрастает до α ≈ 1.0, и при Z > 29 параметр α = 1.0. Например, для электронов пучка с энергией 20 кэВ в серебре Ltr = 0.22 мкм, а Re = 0.86 мкм, α = 1.0. В золоте Ltr = 0.10 мкм, а Re = 0.51 мкм, соответственно, α = 1.0.
Распределение ОРЭ первой группы по глубине образца [4]:
где:(2)
$\begin{gathered} \frac{{d{{\eta }_{1}}}}{{dt}} = \frac{{1.15{{\eta }_{1}}}}{{\sqrt \pi {{z}_{{ss}}}}}\exp \left\{ { - \frac{{{{{\left( {t - 0.77{{z}_{{ss}}}} \right)}}^{2}}}}{{z_{{ss}}^{2}}}} \right\}, \\ {{z}_{{ss}}} = {\text{ }}{{L}_{{{\text{tr}}}}}{{Z}^{{--1/3}}}. \\ \end{gathered} $Распределение ОРЭ второй группы по глубине образца [5]:
где:(4)
$\frac{{d{{\eta }_{2}}}}{{dt}} = \frac{{1.40{{\eta }_{2}}}}{{\sqrt \pi {{z}_{{ss}}}}}\exp \left\{ { - \frac{{{{{\left( {t - 0.77{{L}_{{tr}}}} \right)}}^{2}}}}{{1.25{{z}_{{ss}}}\left( {{{z}_{{ss}}} - 0.2t} \right)}}} \right\}.$Примеры достигаемого соответствия расчетов распределений коэффициентов η1 и η2 по глубине в алюминии и меди по формулам (1)–(4) с экспериментом приведены на рис. 1а и 1б.
Из представленных результатов следует, во-первых, что в тонких слоях толщиной h ≤ zss присутствие ОРЭ второй группы пренебрежимо мало. Во-вторых, что для любой однородной пленочной мишени, представленной в виде суммы n тонких одинаковой толщины слоев, число ОРЭ первой группы ∆η1i, которые отразились в i-ом слое толщиной ∆zi, расположенном на глубине zi, всегда может быть определено с приемлемой точностью. Поэтому, если применить к этой группе ОРЭ результаты разработанной в [7–10] статистической модели неупругого многократного рассеяния, то энергетический спектр Ni(E) группы ∆η1i ОРЭ, которые отразились в i-ом слое ∆zi, расположенном на глубине zi, и которые характеризуются наиболее вероятной энергией Epi, представляется в аналитическом виде как:
(5)
${{N}_{i}}\left( E \right) = \Delta {{\eta }_{{1i}}}\frac{{{{A}_{i}}}}{{{{E}_{0}}}}\exp \left\{ { - \frac{{{{{\left( {E - {{E}_{{pi}}}} \right)}}^{2}}}}{{0.125{{E}_{0}} \times \left( {{{E}_{0}} - E} \right)}}} \right\},$И если далее провести суммирование по всем n слоям, то получим энергетический спектр NB(E) ОРЭ всей первой группы электронов, вышедших в полупространство над поверхностью пленочного образца:
Этот алгоритм расчета энергетического спектра электронов соответствует случаю регистрации ОРЭ с помощью сеточного квазисферического энергетического анализатора [4, 18]. В современных РЭМ для регистрации энергетических спектров ОРЭ, как правило, используются спектрометры анализаторы дисперсионного типа, обеспечивающие сбор и анализ ОРЭ, распространяющихся в сторону входной диафрагмы спектрометра в пределах небольшого телесного угла ∆Ω. Вследствие того, что отраженные в каждом тонком слое электроны характеризуются и определенным распределением по полярному углу θ (с углом наиболее вероятного направления θpi), а входная диафрагма спектрометра-анализатора в камере образцов растрового электронного микроскопа располагается под некоторым углом θ0 к месту падения пучка электронов на поверхность образца, то необходимо учесть влияние этого пространственного распределения Ni = Ni(E, θ0, θpi) на регистрируемый сигнал NB(E) ОРЭ.
РАСПРЕДЕЛЕНИЕ ОРЭ ПО ПОЛЯРНОМУ УГЛУ θ ПРИ НОРМАЛЬНОМ ПАДЕНИИ ЭЛЕКТРОНОВ НА ПЛЕНКУ, ТОЛЩИНОЙ h
При нормальном падении на поверхность образца пучка заряженных частиц распределение F(θ, φ) частиц в пространстве зависит явно только от полярного угла θ и не зависит явно от азимутального угла φ. Как следствие, экспериментально регистрируемые спектры распределения частиц, как прошедших пленочную мишень, так и отраженных от нее, строго симметричны относительно первоначального направления их движения. Это дает основание использовать полученное в работе [9] решение транспортного уравнения для описания двумерного процесса многократного рассеяния частиц в веществе. Пусть θ – угловое отклонение частицы от первоначального направления пучка после прохождения пленки известной толщины h, а величина θmax – максимальное вероятное угловое отклонение частицы в результате однократного рассеяния. Тогда для прошедших пленку электронов, в случае малоуглового рассеяния частиц, их функция распределения по углу FT(θ, θmax) хорошо описывается выражением работы [7]:
(6)
$\begin{gathered} {{F}_{T}}\left( {{\theta ,}{{{\theta }}_{p}}} \right) = \frac{{d{{N}_{T}}}}{{d\Omega }} = \frac{{d{{N}_{T}}}}{{2{\pi \theta }d{\theta }}} = \\ = A\exp \left[ { - \frac{{{{{\theta }}^{2}}}}{{0.75{{{\theta }}_{{\max }}}\left( {{{{\theta }}_{{\max }}} + {\theta }} \right)}}} \right]. \\ \end{gathered} $Применительно к ОРЭ первой группы при нормальном падении пучка электронов на образец θmax = π/4 [7] и для них функция распределение по углу θ – FB(θ, θp, θmax) – принимает вид:
(7)
$\begin{gathered} {{F}_{B}} = \frac{{d{{\eta }_{1}}}}{{2\pi \sin \theta d\theta }} = {{\eta }_{1}}\left( h \right)\frac{{1.0826}}{\pi } \\ \times \,\,\exp \left\{ { - \frac{{{{{\left( {{{\theta }_{p}} - \theta } \right)}}^{2}}}}{{0.75\pi /4\left( {\pi /4 + \theta } \right)}}} \right\}\Phi \left( {\frac{{\theta - \pi /2}}{{\sqrt 2 \pi /8}}} \right). \\ \end{gathered} $Среди научных публикаций число экспериментальных работ, посвященных исследованию угловых распределений ОРЭ от свободных пленок, – единицы. Для проверки формулы (7) были использованы экспериментальные данные по угловому распределению ОРЭ от поликристаллических пленок золота, полученные в работе [19]. На рис. 2б приведены экспериментальные спектры из этой работы. На рис. 2а представлены результаты модельных расчетов, выполненные по формуле (7) для пленок золота такой же толщины и для такой же энергии первичных электронов E0= = 30 кэВ. В целом, можно говорить о приемлемом соответствии модельных расчетов экспериментальным данным. Особенностью формулы (7) является возможность ее применения к объектам в виде массивных мишеней. Действительно, при значительном возрастании толщины пленки h → Re изменения в (7) коснутся лишь введения в нее двух вполне определенных величин: π как θp и η1 = = η1(h). При этом возможности более детальной проверки формулы существенно расширяются из-за наличия большого массива выполненных экспериментальных измерений. Итоговые результаты представлены на рис. 3 и рис. 4. Следует отметить широкий диапазон материалов мишеней, привлеченных к апробации полученной формулы: от бериллия до урана. А также то, что формула хорошо описывает не только экспериментальные угловые распределения η1 ОРЭ первой группы (бериллий и алюминий), но и спектр суммарного коэффициента η обратного рассеяния электронов (медь, серебро, золото). Это свидетельствует об изоморфности угловых распределений ОРЭ первой и второй групп при нормальном падении первичных электронов на поверхность образца.
ФУНКЦИЯ ЭНЕРГЕТИЧЕСКОГО РАСПРЕДЕЛЕНИЯ ОРЭ, ОТ ПЛЕНОЧНОЙ МИШЕНИ ИЗВЕСТНОГО СОСТАВА И ЗАДАННОЙ ГЕОМЕТРИИ
С учетом того, что входная диафрагма энергетического анализатора располагается под известным нам углом θ0, а угол наиболее вероятного направления θpi этих ОРЭ также известен, то вклад ∆Ni(E, θ0, θpi) этого слоя в спектр регистрируемых ОРЭ будет:
Проводя суммирование по всем слоям – получаем искомый энергетический спектр ОРЭ первой группы NB:
Результат такого подхода и достигаемое при этом соответствие экспериментальному спектру, на примере свободной алюминиевой пленочной мишени, представлен на рис. 5. Штрихпунктирные линии соответствуют вкладам в общий спектр первой группы ОРЭ отдельных распределений от четырех тонких слоев, на которые была разбита пленка. На рис. 6 приведены результаты расчета энергетического спектра ОРЭ для массивной золотой мишени. Энергия первичных электронов 15 кэВ и угол регистрации ОРЭ 155°. Пунктирная кривая 1 представляет рассчитанное распределение ОРЭ первой группы, как результат суммирования 24 отдельных распределений ∆Ni(E, θ0, θpi). Это соответствует представлению всей области ОРЭ первой группы в образце в виде 24 одинаковых слоев. Пунктирная кривая 2 представляет рассчитанный вклад распределения ОРЭ второй группы. Итоговое энергетическое распределение в виде сплошной линии хорошо согласуется с экспериментальным спектром работы [18]. И наконец, итоговые результаты решения прямой задачи по расчету спектра NB(E, θ0) ОРЭ как совокупности энергетических спектров электронов, отраженных от отдельных слоев структуры с заданной геометрией и известным элементным составом, представлены на рис. 7. Расчеты выполнялись для трехслойной тонкопленочной структуры на массивной подложке: Cu (толщина слоя 0.09 мкм, кривая 1)–Al (толщина слоя 0.15 мкм, кривая 2)–Cu (толщина слоя 0.09 мкм, кривая 3)–Si (подложка, кривая 4). Спектр от каждого слоя структуры рассчитывался путем разбиения каждого слоя на определенное число одинаковых подслоев: медные слои – на 6, алюминиевый слой – на 2, подложка – на 24 подслоя. При этих расчетах предполагалось, что число ОРЭ первой группы, отраженных в каждом тонком подслое, такое же, как в однородной мишени, а меняются лишь их энергетические потери при прохождении слоев иного элементного состава. Использованная в эксперименте и при расчете величина энергии первичных электронов E0 = 20 кэВ обеспечивает формирование в анализируемых слоях структуры лишь ОРЭ первой группы. Приведенное на рис. 7 сопоставление результатов расчета с экспериментальным спектром из работы [3] показывает удовлетворительное соответствие расчета и эксперимента практически во всем диапазоне энергий ОРЭ.
ЗАКЛЮЧЕНИЕ
В результате проведенных исследований получена формула, описывающая двумерное распределение обратно рассеянных электронов по полярному углу после их выхода из пленочной мишени заданной толщины. Разработан алгоритм расчетов, необходимых для построения энергетических распределений ОРЭ первой группы при заданном значении угла θ0 входной диафрагмы анализатора при их регистрации. Проведены модельные расчеты функции энергетического распределения частиц для пучка электронов, вышедших из пленочной мишени известного состава и заданной геометрии. Показано, что при использовании двухпотоковой модели обратного рассеяния электронов, достигается хорошее соответствие модельных расчетов экспериментальным спектрам.
Список литературы
Спивак Г.В., Рау Э.И., Филиппов М.Н., Сасов А.Ю. // Современная электронная микроскопия в исследовании вещества. М.: Наука, 1982. С. 5.
Аристов В.В, Дремова Н.Н., Зайцев С.И., Рау Э.И. и др. // ДАН СССР. 1988. Т. 301. С. 611.
Дремова Н.Н., Зайцев С.И., Рау Э.И., Якимов Е.Б. // Известия РАН. Сер. Физ.1993. Т. 57. № 8. С. 9.
Михеев Н.Н, Петров В.И., Степович М.А. // Известия РАН. Сер. Физ. 1995. Т. 59. № 2. С. 144.
Михеев Н.Н., Степович М.А. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 1998. № 1. С. 6.
Михеев Н.Н., Никифорова Н.А., Ганчев А.С. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2015. № 9.С. 64. https://doi.org/10.7868/S0207352815090115
Михеев Н.Н. // Поверхность. Рентгеновские, синхротронные и нейтронные исследования. 1998. № 10. С. 37.
Михеев Н.Н. // Известия РАН. Сер. Физ. 2000. Т. 64. № 11. С. 2137.
Михеев Н.Н., Степович М.А., Юдина С.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2009. № 3. С. 53.
Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2010. № 4. С. 25.
Niedrig H. // J. Appl. Phys. 1982. V. 53. № 4. R 15.
Михеев Н.Н. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2019. № 8. С. 56. https://doi.org/10.1134/S0207352819080109
Михеев Н.Н., Колесник А.С. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2017. № 12. С. 88. https://doi.org/10.7868/S0207352817120083
Михеев Н.Н., Степович М.А., Широкова Е.В. // Известия РАН, Сер. Физ. 2010. Т. 74. № 7. С. 1049.
Рау Э.И., Сенов Р.А. // Известия РАН, Сер. Физ. 2004. Т. 68. № 9. С. 1342.
Niedrig H., Sieber P. // Z. Angew. Phys. 1971. B. 31. S. 327.
Cosslett V.E., Thomas R.N. // Brit. J. Appl. Phys. 1965. V. 16. P. 779.
Kulenkampff H., Spira W. // Z. Phys. 1954. B 137. S. 414.
Hohn F.J., Niedrig H. // Optik. 1972. B. 35. S. 290.
Balladore J.L., Soum G., Martinez J.P., Arnal F. // J. Microsc. Spectrosc. Electron.1980. V. 5. P. 371.
Бишоп Г. // Основы рентгеноспектрального локального анализа, перевод с англ. / Под ред. Боровского И.Б. М.: Наука, 1973. С. 109.
Kanter H. // Ann. Phys. 1957. B. 20. S. 144.
Рязанов М.И., Тилинин И.С. Исследование поверхности по обратному рассеянию частиц. М.: Энергоатомиздат, 1985. 152 с.
Kulenkampff H., Ruttiger K. // Z. Phys. 1958. B. 152. S. 249.
Купреенко С.Ю., Орликовский Н.А., Рау Э.И., Тагаченков А.М., Татаринцев А.А. // ЖТФ. 2015. Т. 85. С. 101.
Дополнительные материалы отсутствуют.
Инструменты
Поверхность. Рентгеновские, синхротронные и нейтронные исследования