О некоторых оценках возможности необратимого термического воздействия электронного зонда на анализируемый микрообъем. Результаты математического моделирования

Язык труда и переводы:
УДК:
53.043:537.533.9:54.062:004.942
Дата публикации:
04 декабря 2022, 00:35
Категория:
Математическое моделирование физических процессов и технических систем
Авторы
Аннотация:
Рассмотрены процессы потерь энергии в конденсированном веществе, обусловленные облучением мишени остро сфокусированным электронным пучком, и, как следствие этого, нагрев облучаемого образца. Рассмотрение проведено для математической модели, основанной на возможности раздельного количественного описания вклада в потери энергии поглощенных в мишени и обратно рассеянных электронов. Решение стационарного уравнения теплопроводности получено с использованием функции Грина. Показано, что для тяжёлых элементов зависимость перегрева мишени от энергии электронов пучка имеет немонотонный характер, что для таких материалов может привести к необратимым процессам в исследуемых материалах.
Ключевые слова:
электронный зонд, нагрев мишени, зависимость от энергии электронов зонда, уравнение теплопроводности
Основной текст труда

Введение

При использовании электронно-зондовых приборов для диагностики материалов плотность энергии, выделяемой электронным пучком в мишени, может колебаться примерно от 10 до 1000 Вт/см3, а характерный размер микрообъема потерь энергии электронами пучка может быть около 1 или 1,5 микрометра [1–4]. Таким образом, когда мишень облучается остро сфокусированным электронным пучком, нагрев образца может быть весьма значительным, что может привести к необратимым процессам в мишени, например, плавлению или сублимации компонентов в вакуум и изменению состава мишени. Поэтому осо­бую ценность на прак­тике приоб­ре­тают рас­чет­ные оцен­ки ве­ли­чи­ны нагрева элек­трон­ным зондом, основанные на ре­ше­нии урав­не­ния теп­ло­про­вод­нос­ти, позволяющие оценивать максимальную температуру нагрева образца.

Математическая модель, описывающая потери энергии электронным пучком в мишени

При взаимодействии киловольтных (до 50 кэВ) электронов с конденсированным веществом большая часть энергия электронов идет на нагрев, а распределение тепловых источников будет совпадать с распределением потерь энергии первичными электронами [5, 6]. Большинство из использовавшихся ранее математических моделей описывают энергетические потери в виде функции Гаусса и потому дают лишь качественное соответствие с экспериментальными результатами, причем применение каждой из них ограничивается сравнительное узким кругом материалов и энергий электронов пучка. Анализ существующих моделей довольно подробно выполнен ранее [5, 6] и литературу там же. На основе проведённого анализа предложена новая модель, основанная на возможности раздельного количественного описания вклада в потери энергии поглощенных в мишени и обратно рассеянных электронов. В этой модели в качестве функции генерации источника тепла использовано выражение

P(x,\ y,\ z)={\frac {1,085(1-\eta )P_{0}}{\pi ^{3/2}a_{1}^{2}z_{ms}(1-\eta +\eta {\frac {z_{ss}}{z_{ms}}})}}\left\{\exp \left\{-\left[{\frac {x^{2}+y^{2}}{a_{1}^{2}}}+\left({\frac {z-z_{ms}}{z_{ms}}}\right)^{2}\right]\right\}+{\frac {\eta a_{1}^{2}}{\left(1-\eta \right)a_{2}^{2}}}\exp \left\{-\left[{\frac {x^{2}+y^{2}}{a_{2}^{2}}}+\left({\frac {z-z_{ss}}{z_{ss}}}\right)^{2}\right]\right\}\right\}.

Здесь начало координат совпадает с точкой падения электронного зонда на образец, координаты x и y расположены на плоской поверхности мишени, а ось z направлена в глубину мишени; P_{0} — мощность пучка первичных электронов; z_{ms} — глуби­на максимальных потерь энергии первичными электронами, испытавшими малоугловое рассеяние и поглощенными мишенью;  z_{ss} — глубина максимальных потерь энергии обратно рассеянными электронами, при этом z_{ss}=Z^{-1/3}z_{ms} , где Z — эффективный порядковый но­мер материала мишени. В формуле также используются следующие параметры: \eta — коэффициент обратного рас­сеяния электронов пучка, \eta =\left(1,085ae/\pi ^{1/2}z_{ms}\right)\int _{0}^{z}\exp \left[-\left(z-z_{ss}\right)^{2}/z_{ss}^{2}\right]dz , постоянная a=0,024\ \left(Z^{2}/A\right) , e — основание натуральных логарифмов, при этом для массивной мишени \eta \ =\ \int _{0}^{\infty }{\eta \ \left(z\right)dz\ \approx \ 0,024eZ^{1,67}/A}, а A — эффективная относительная атомная масса материала мишени. Параметры a_{1}  и a_{2} , входящие в выражения для z_{ms}  и z_{ss} , определяются согласно [5, 6]. Сравнение с экспериментальными результатами [5, 6] показало, что эта модель может быть успешно использована при проведении практических расчетов для широкого класса материалов (практически от Al по Au) в широком диапазоне энергий первичных электронов (практически от 2 до 50 кэВ).

Математическая модель, описывающая нагрев мишени электронным пучком

Время, в течение которого электронный зонд позиционируется на поверхности образца и установление стационарного температурного режима анализируемого микрообъема, как правило, значительно меньше времени процессов, формирующих регистрируемый информативный сигнал [5]. Это позволяет искать распределение температуры в микрообъёме мишени на основе решения стационарного уравнения теплопроводности

\left({\frac {\partial ^{2}\Delta \ T\left(x,\ y,\ z\right)}{\partial x^{2}}}+{\frac {\partial ^{2}\Delta \ T\left(x,\ y,\ z\right)}{\partial y^{2}}}+{\frac {\partial ^{2}\Delta \ T\left(x,\ y,\ z\right)}{\partial z^{2}}}\right)=\ -{\frac {P(x,\ y,\ z)}{k}}.

Здесь \Delta T\left(x,\ y,\ z\right)=T\left(x,\ y,\ z\right)-T_{0} , где T  – температура образца в точке (x,\ y,\ z)  после установления стационарного режима под воздействием электронного пучка, T_{0} — температура образца до воздействия электронного пучка; k — коэффициент теплопроводности, а \Delta T(x,\ y,\ z) удовлетворяет следующим граничным условиям: \lim \limits _{x\rightarrow -\infty }{\Delta T=0} , \lim \limits _{x\rightarrow +\infty }{\Delta T=0} , \lim \limits _{y\rightarrow -\infty }{\Delta T=0} , \lim \limits _{y\rightarrow +\infty }{\Delta T=0} , \lim \limits _{z\rightarrow +\infty }{\Delta T=0} .

Поскольку в электронно-зондовых технологиях облучение мишеней электронным пучком в основном проводят в условиях, близких к вакууму, то можно задать следующее граничное условие, означающее отсутствие теплообмена с внешней средой: \left.{\frac {\partial \Delta T}{\partial z}}\right|_{z=0}=0.

Уравнение теплопроводности можно решить с использованием функции Грина G.  Со­от­вет­ст­ву­ю­щее решение будет иметь вид: \Delta T\ \left(x,\ y,\ z\right)=\ \int _{D}{G\left(x_{0},\ y_{0},\ z_{0},\ x,\ y,\ z\right){\frac {P\left(x_{0},\ y_{0},\ z_{0}\right)}{k}}dV_{\varrho }},

где D — область распространения тепла, (-\infty <x_{0}<+\infty , -\infty <y_{0}<+\infty , 0\leq z_{0}<+\infty ), dV_{\varrho }=dx_{0}dy_{0}dz_{0} , а G={\frac {1}{4\pi }}\left({\frac {1}{\sqrt {\left(x-x_{0}\right)^{2}+\left(y-y_{0}\right)^{2}+\left(z-z_{0}\right)^{2}}}}+{\frac {1}{\sqrt {\left(x-x_{0}\right)^{2}+\left(y-y_{0}\right)^{2}+\left(z+z_{0}\right)^{2}}}}\right)

Здесь \left(x_{0},\ y_{0},\ z_{0}\right) — координата точечного источника.

Результаты расчетов

На рисунке представлены рассчитанные зависимости перегрева мишеней \Delta T\ \left(x,\ y,\ z\right)  от энергии электронов пучка E_{0} для алюминия (легкий металл), меди (металл со средним атомным номером) и золота (тяжелый металл). Показан также вклад поглощенных в мишени и отраженных электронов. Расчеты проведены при x  и y , равными нулю, при токе зонда 10–7 A. Параметры материалов приведены в [7].

Зависимость перегрева мишеней от энергии электронов пучка (кривые 1) для алюминия (а), меди (б), золота (в),  вклад поглощенных электронов (кривые 2) и отраженных электронов (кривые 3)

Анализ результатов расчетов показал, что для тяжелых металлов характерным оказалась немонотонная зависимость \Delta T\ \left(E_{0}\right) . Для того чтобы избежать возможных необратимых процессов в мишени, для ряда материалов это может приводить к необходимости учета этого явления при выборе режима облучения мишени.

Заключение

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

Грант
Работа выполнена при финансовой поддержке гранта РФФИ, проект № 19-03-00271
Литература
  1. Филиппов М.Н. Оценка теплового воздействия электронного зонда в растровой электронной микроскопии и рентгеноспектральном анализе. Известия. РАН. Серия физическая, 1993, № 8, с. 163–171.
  2. Масловская А.Г. Моделирование теплового воздействия электронного зонда в растровой электронной микроскопии. Моделирование систем и процессов, 2007, № 2 (14), с. 40–51.
  3. Кузин А.Ю., Степович М.А., Митюхляев В.Б., Тодуа П.А., Филиппов М.Н. Тепловые эффекты при низковольтном электронно-зондовом рентгеноспектральном микроанализе с нанометровой локальностью. Измерительная техника, 2016, № 10, с. 27–29.
  4. Амрастанов А.Н., Кузин А.Ю., Митюхляев В.Б., Серегина Е.В., Степович М.А., Тодуа П.А., Филиппов М.Н. Тепловое воздействие электронного зонда при рентгеноспектральном наноанализе. Измерительная техника, 2017, № 6, с. 13–15.
  5. Михеев Н.Н., Петров В.И., Степович М.А. Количественный анализ материалов полупроводниковой оптоэлектроники методами растровой электронной микроскопии. Известия АН СССР. Серия физическая, 1991, т. 55, № 8, с. 1474–1482.
  6. Михеев Н.Н., Степович М.А. Распределение энергетических потерь при взаимодействии электронного зонда с веществом. Заводская лаборатория. Диагностика материалов, 1996, т. 62, № 4, с. 20–25.
  7. Амрастанов А.Н., Серегина Е.В., Степович М.А. Оценка нагрева поверхности однородной металлической мишени электронным зондом. Известия РАН. Серия физическая, 2019, т. 83, № 11, с. 1455–1460.
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.