В ядерной физике и физике материалов тормозная способность — это замедляющая сила, действующая на заряженные частицы , как правило, альфа- и бета-частицы , из-за взаимодействия с веществом , что приводит к потере кинетической энергии частицы . [1] [2] Тормозная способность также интерпретируется как скорость, с которой материал поглощает кинетическую энергию заряженной частицы . Ее применение важно в широком спектре термодинамических областей, таких как защита от радиации , ионная имплантация и ядерная медицина . [3]
Как заряженные, так и незаряженные частицы теряют энергию при прохождении через вещество. В большинстве случаев ниже рассматриваются положительные ионы . Тормозная способность зависит от типа и энергии излучения и от свойств материала, через который оно проходит. Поскольку для образования пары ионов (обычно положительного иона и (отрицательного) электрона) требуется фиксированное количество энергии (например, 33,97 эВ в сухом воздухе [4] : 305 ), число ионизаций на длину пути пропорционально тормозной способности. Тормозная способность материала численно равна потере энергии E на единицу длины пути, x :
Знак минус делает S положительным.
Сила обычно увеличивается к концу диапазона и достигает максимума, пика Брэгга , незадолго до того, как энергия падает до нуля. Кривая, описывающая силу как функцию глубины материала, называется кривой Брэгга . Это имеет большое практическое значение для лучевой терапии .
Уравнение выше определяет линейную тормозную способность , которая в международной системе выражается в Н , но обычно указывается в других единицах, таких как МэВ / мм или подобных. Если вещество сравнивается в газообразной и твердой форме, то линейные тормозные способности двух состояний сильно различаются просто из-за разной плотности. Поэтому часто делят силу на плотность материала, чтобы получить массовую тормозную способность , которая в международной системе выражается в м4 / с2 , но обычно находится в единицах, таких как МэВ/(мг/см2 ) или подобных. Тогда массовая тормозная способность очень мало зависит от плотности материала.
На рисунке показано, как тормозная способность альфа-частиц с энергией 5,49 МэВ увеличивается по мере того, как частица проходит через воздух, пока не достигает максимума. Эта энергия соответствует энергии альфа-частиц, излучаемых естественным радиоактивным газом радоном ( 222 Rn), который присутствует в воздухе в ничтожных количествах.
Средний диапазон можно рассчитать путем интегрирования обратной тормозной способности по энергии: [5]
где:
Выделившуюся энергию можно получить путем интегрирования тормозной способности по всей длине пути иона при его движении в материале.
Электронное торможение относится к замедлению иона-снаряда из-за неупругих столкновений между связанными электронами в среде и ионом, движущимся через нее. Термин неупругое используется для обозначения того, что энергия теряется в ходе процесса (столкновения могут приводить как к возбуждению связанных электронов среды, так и к возбуждению электронного облака иона). Линейная электронная тормозная способность идентична неограниченной линейной передаче энергии .
Вместо передачи энергии некоторые модели рассматривают электронную тормозную способность как передачу импульса между электронным газом и энергичным ионом. Это согласуется с результатом Бете в диапазоне высоких энергий. [6]
Поскольку число столкновений иона с электронами велико, а зарядовое состояние иона при прохождении через среду может часто меняться, очень сложно описать все возможные взаимодействия для всех возможных зарядовых состояний ионов. Вместо этого электронная тормозная способность часто задается как простая функция энергии, которая является средним значением, взятым по всем процессам потери энергии для различных зарядовых состояний. Ее можно теоретически определить с точностью до нескольких % в диапазоне энергий выше нескольких сотен кэВ на нуклон из теоретических расчетов, наиболее известной из которых является формула Бете . При энергиях ниже примерно 100 кэВ на нуклон становится сложнее определить электронное торможение с помощью аналитических моделей. [7] Недавно теория функционала плотности, зависящая от времени, в реальном времени была успешно использована для точного определения электронного торможения для различных систем ионов-мишеней в широком диапазоне энергий, включая режим низкой энергии. [8] [9]
Графические представления экспериментальных значений электронной тормозной способности для многих ионов во многих веществах были даны Полом. [10] Точность различных таблиц торможения была определена с помощью статистических сравнений. [11]
Ядерная тормозная способность относится к упругим столкновениям между ионом-снарядом и атомами в образце (устоявшееся обозначение «ядерный» может сбивать с толку, поскольку ядерное торможение не обусловлено ядерными силами, [12] но оно призвано отметить, что этот тип торможения включает взаимодействие иона с ядрами в мишени). Если знать форму потенциальной энергии отталкивания между двумя атомами (см. ниже), можно рассчитать ядерную тормозную способность . На рисунке тормозной способности, показанном выше для ионов алюминия в алюминии, ядерное торможение пренебрежимо мало, за исключением самой низкой энергии. Ядерное торможение увеличивается с увеличением массы иона. На рисунке, показанном справа, ядерное торможение больше электронного торможения при низкой энергии. Для очень легких ионов, замедляющихся в тяжелых материалах, ядерное торможение слабее электронного при всех энергиях.
Особенно в области радиационного повреждения в детекторах термин « неионизирующие потери энергии » (NIEL) используется как термин, противоположный линейной передаче энергии (LET), см., например, ссылки [13] [14] [15] . Поскольку по определению ядерная тормозная способность не включает в себя электронные возбуждения, NIEL и ядерное торможение можно считать одной и той же величиной при отсутствии ядерных реакций.
Таким образом, общая нерелятивистская тормозная способность представляет собой сумму двух членов: . Было разработано несколько полуэмпирических формул тормозной способности. Модель, предложенная Циглером, Бирсаком и Литтмарком (так называемая «ZBL»-остановка, см. следующую главу), [16] [17], реализованная в различных версиях кодов TRIM/SRIM , [18] используется сегодня чаще всего.
Радиационную тормозную способность , которая обусловлена испусканием тормозного излучения в электрических полях частиц в проходимом материале, необходимо учитывать при чрезвычайно высоких энергиях ионов. [12] Для электронных снарядов радиационная тормозная способность всегда важна. При высоких энергиях ионов могут также иметь место потери энергии из-за ядерных реакций, но такие процессы обычно не описываются тормозной способностью. [12]
Вблизи поверхности твердого материала мишени как ядерное, так и электронное торможение может привести к распылению .
В начале процесса замедления при высоких энергиях ион замедляется в основном электронным торможением и движется почти по прямой. Когда ион достаточно замедлится, столкновения с ядрами (ядерное торможение) становятся все более и более вероятными, в конечном итоге доминируя над замедлением. Когда атомы твердого тела получают значительные энергии отдачи при ударе иона, они будут удалены из своих положений в решетке и произведут каскад дальнейших столкновений в материале. Эти каскады столкновений являются основной причиной возникновения повреждений при ионной имплантации в металлы и полупроводники.
Когда энергии всех атомов в системе падают ниже пороговой энергии смещения , производство новых повреждений прекращается, и концепция ядерного торможения больше не имеет смысла. Общее количество энергии, переданной ядерными столкновениями атомам в материалах, называется ядерной депонированной энергией.
Вставка на рисунке показывает типичное распределение пробега ионов, осажденных в твердом теле. Показанный здесь случай может быть, например, замедлением иона кремния с энергией 1 МэВ в кремнии. Средний пробег иона с энергией 1 МэВ обычно находится в микрометровом диапазоне.
На очень малых расстояниях между ядрами отталкивающее взаимодействие можно считать по сути кулоновским. На больших расстояниях электронные облака экранируют ядра друг от друга. Таким образом, отталкивающий потенциал можно описать, умножив кулоновское отталкивание между ядрами на экранирующую функцию φ(r/a),
где φ(r/a) → 1 при r → 0. Здесь и — заряды взаимодействующих ядер, а r — расстояние между ними; a — так называемый параметр экранирования.
За эти годы было предложено большое количество различных отталкивающих потенциалов и экранирующих функций, некоторые из которых были определены полуэмпирически, другие — теоретическими расчетами. Часто используемый отталкивающий потенциал — это потенциал, предложенный Циглером, Бирсаком и Литтмарком, так называемый отталкивающий потенциал ZBL. Он был построен путем подгонки универсальной экранирующей функции к теоретически полученным потенциалам, рассчитанным для большого разнообразия пар атомов. [16] Параметр и функция экранирования ZBL имеют вид
и
где x = r/a u , а a 0 — радиус атома Бора = 0,529 Å.
Стандартное отклонение соответствия универсального отталкивающего потенциала ZBL теоретически рассчитанным парно-специфическим потенциалам, которым он соответствует, составляет 18% выше 2 эВ. [16] Еще более точные отталкивающие потенциалы можно получить из самосогласованных расчетов полной энергии с использованием теории функционала плотности и приближения локальной плотности (LDA) для электронного обмена и корреляции. [19]
В кристаллических материалах ион в некоторых случаях может «каналироваться», т. е. фокусироваться в канале между кристаллическими плоскостями, где он почти не испытывает столкновений с ядрами. [20] Кроме того, электронная тормозная способность может быть слабее в канале. [21] Таким образом, ядерная и электронная тормозная способность зависят не только от типа и плотности материала, но также от его микроскопической структуры и поперечного сечения.
Методы компьютерного моделирования для расчета движения ионов в среде разрабатывались с 1960-х годов и в настоящее время являются доминирующим способом теоретической обработки тормозной способности. Основная идея в них заключается в том, чтобы отслеживать движение иона в среде, моделируя столкновения с ядрами в среде. Электронная тормозная способность обычно учитывается как сила трения, замедляющая ион.
Традиционные методы, используемые для расчета пробегов ионов, основаны на приближении бинарных столкновений (BCA). [22] В этих методах движение ионов в имплантированном образце рассматривается как последовательность отдельных столкновений между ионом отдачи и атомами в образце. Для каждого отдельного столкновения классический интеграл рассеяния решается путем численного интегрирования.
Параметр удара p в интеграле рассеяния определяется либо из стохастического распределения, либо способом, учитывающим кристаллическую структуру образца. Первый метод подходит только для моделирования имплантации в аморфные материалы, поскольку не учитывает каналирование.
Самая известная программа моделирования BCA — TRIM/SRIM ( аббревиатура от TRansport of Ions in Matter, в более поздних версиях называется Stopping and Range of Ions in Matter), которая основана на электронном торможении и межатомном потенциале ZBL . [16] [18] [23] Она имеет очень простой в использовании пользовательский интерфейс и имеет параметры по умолчанию для всех ионов во всех материалах вплоть до энергии ионов 1 ГэВ, что сделало ее чрезвычайно популярной. Однако она не учитывает кристаллическую структуру, что серьезно ограничивает ее полезность во многих случаях. Несколько программ BCA преодолевают эту трудность; некоторые из них довольно известны — MARLOWE, [24] BCCRYS и crystal-TRIM.
Хотя методы BCA успешно применялись для описания многих физических процессов, у них есть некоторые препятствия для реалистичного описания процесса замедления энергичных ионов. Основное предположение о том, что столкновения являются бинарными, приводит к серьезным проблемам при попытке учесть множественные взаимодействия. Кроме того, при моделировании кристаллических материалов процесс выбора следующего сталкивающегося атома решетки и параметр удара p всегда включают несколько параметров, которые могут не иметь идеально определенных значений, что может повлиять на результаты на 10–20% даже для вполне разумных, на первый взгляд, выборов значений параметров. Наилучшая надежность в BCA достигается путем включения в расчеты множественных столкновений, что нелегко сделать правильно. Однако, по крайней мере, MARLOWE делает это.
Принципиально более простой способ моделирования множественных атомных столкновений обеспечивается моделированием молекулярной динамики (МД), в котором временная эволюция системы атомов рассчитывается путем численного решения уравнений движения. Были разработаны специальные методы МД, в которых число взаимодействий и атомов, участвующих в моделировании МД, было уменьшено, чтобы сделать их достаточно эффективными для расчета пробегов ионов. [25] [26] Моделирование МД автоматически описывает ядерную тормозную способность. Электронную тормозную способность можно легко включить в моделирование молекулярной динамики либо как силу трения [25] [27] [28] [29] [26] [30] [31] [32] или более продвинутым способом, также отслеживая нагрев электронных систем и связывая электронные и атомные степени свободы. [33] [34] [35]
За пределами максимума тормозная способность уменьшается примерно как 1/v 2 с ростом скорости частицы v , но после минимума она снова увеличивается. [36] Минимальная ионизирующая частица (MIP) — это частица, средняя скорость потери энергии которой через вещество близка к минимуму. Во многих практических случаях релятивистские частицы (например, мюоны космических лучей ) являются минимально ионизирующими частицами. Важным свойством всех минимально ионизирующих частиц является то, что приблизительно верно, где и — обычные релятивистские кинематические величины. Более того, все MIP имеют почти одинаковую потерю энергии в веществе, значение которой равно: . [36]
{{cite book}}
: CS1 maint: location missing publisher (link)