
книги из ГПНТБ / Физические основы рентгеноспектрального локального анализа
..pdfПРОГРЕСС В ВЫЧИСЛЕНИЙ ТОРМОЗНОЙ СПОСОБНОСТИ И ФАКТОРА ОБРАТНОГО РАССЕЯНИЯ
П. Дапкамб, С. Рид
Введение. Кастен [1] показал, что весовая концентра ция СА анализируемого элемента А в образце и отношение интенсивностей анализируемого излучения (обычно К а линин) от образца и эталона — чистого элемента А — в первом приближении связаны соотношением
кА = сл- |
(1) |
При этом предполагается, что перед определением |
кА |
в интенсивности внесены аппаратурные поправки, в их чис ле — на мертвое время и фон, и что интенсивности излуче ния от образца и стандарта измерены при одинаковых условиях опыта. Отклонения от простой пропорционально сти (1) могут быть большими, и при необходимости точно го количественного анализа приходится вводить поправоч ный множитель:
кА — сА • (поправочный множитель). |
(2) |
Этот множитель введен в правую часть соотношения (2) по той причине, что он сам является функцией сА , а также концентраций других элементов, присутствующих в об разце.
Следуя Кастену, большинство исследователей разби вало поправочный фактор на три отдельных множителя, соответствующих поправкам на поглощение, флуоресцен цию и атомный номер. Методы расчета поправок на погло щение характеристического рентгеновского излучения в образце и на дополнительное возбуждение анализируемо го излучения флуоресценцией других элементов, присут ствующих в образце, были предложены Кастеном и впо следствии совершенствовались многими авторами. Третья из поправок, так называемая «поправка на атомный но мер», учитывает нелинейность соотношения между А'А и СА, остающуюся после введения других двух поправок.
Более правильным было бы этот множитель называть «поправкой на возбуждение», поскольку его величина оп ределяется соотношением между интенсивностью харак теристических рентгеновских лучей, возбужденных в об разце, и составом образца. Название же «поправка на атомный номер» связано с тем, что излучение элементом характеристических лучей в первую очередь определяется его атомным номером.
В рентгеновском микроанализе электроны, ускоренные до энергии 5—40 кэв, фокусируются в пучок диаметром примерно 1 мкм или несколько меньше; ток пучка обычно заключен в пределах от 10~6 до 10~9 а. Падающие электро ны взаимодействуют с атомами, входящими в состав об разца. Эти взаимодействия можно разделить на два широ ких класса — упругое рассеяние, при котором электроны существенно изменяют направления своего движения, но теряют лишь пренебрежимо малую энергию, и неупругое рассеяние, при котором происходят существенные потери энергии. Основной вклад в упругое рассеяние на большие углы вносят взаимодействия электронов с ядрами атомов, а неупругое рассеяние в основном связано с взаимодейст виями падающих и атомных электронов. Пространствен ное распределение падающих электронов в образце опре деляется главным образом их упругим рассеянием, тогда как неупругое рассеяние определяет их торможение и путь, который они пройдут в веществе до остановки. Одним из процессов неупругого рассеяния является ионизация внут ренней электронной оболочки атома, например, К- или L-оболочки, сопровождающаяся в некоторой определен ной доле возбуждением квантов характеристического рент геновского излучения К- или L-серии. Однако ионизация внутренней электронной оболочки — сравнительно редкое событие, и основной вклад в торможение падающих элект ронов вносят гораздо более частые неупругие столкнове ния с электронами внешних оболочек атома. До остановки или вылета из образца падающие на него электроны испы тывают многократное рассеяние (рис. 1). Вследствие этого физическая картина явления отличается большой слож ностью, и это вызывает особые трудности при попытках представить «поправку на атомный номер» в общем и при том достаточно простом виде.
Эффективное сечение ионизации Q определяется как ве роятность на единицу пути электрона с данной энергией ионизации им внутренней электронной оболочки (К-, L -
или М-) атома образца. Значительный накопленный экс периментальный материал позволяет с достаточной точ ностью установить зависимость Q от энергии электронов Е. Поэтому, зная число и энергетическое распределение электронов в образце, можно вычислить интенсивность
характеристического рентгенов |
|
|
|
|
|
|
||||||||
ского излучения. Для этой цели |
|
|
|
|
|
|
||||||||
необходимо |
|
располагать |
мо |
|
|
|
|
|
|
|||||
делью физической картины |
рас |
|
|
|
|
|
|
|||||||
сеяния, |
которая |
позволяла |
бы |
|
|
|
|
|
|
|||||
провести |
математический |
рас |
|
|
|
|
|
|
||||||
чет. Все |
модели подобного рода |
|
|
|
|
|
|
|||||||
неизбежно содержат упрощения |
|
|
|
|
|
|
||||||||
и тем |
самым |
требуют |
экспери |
|
|
|
|
|
|
|||||
ментальной |
проверки |
получен |
|
|
|
|
|
|
||||||
ных с их помощью результатов. |
|
|
|
|
|
|
||||||||
Сравнительно |
простая |
диффу |
|
|
|
|
|
|
||||||
зионная |
модель Арчарда |
[2] не |
|
|
|
|
|
|
||||||
настолько |
близка |
к |
реальной |
Рис. 1. |
Схематическое |
изобра |
||||||||
жение объема, в котором рассе |
||||||||||||||
картине, чтобы служить надеж |
иваются падающие |
на |
образец |
|||||||||||
ной основой |
для |
расчета |
по |
электроны. |
Характеристическое |
|||||||||
излучение |
может возбуждаться |
|||||||||||||
правки на атомный номер. Более |
вдоль |
всей |
части |
траектории |
||||||||||
сложные |
расчеты |
Брауна |
[31, |
электрона вплоть до точки, в |
||||||||||
которой его энергия становится |
||||||||||||||
базирующиеся |
на транспортном |
равной |
энергии |
|
возбуждения |
|||||||||
ЯдНекоторые |
электроны вы |
|||||||||||||
уравнении |
Больцмана в приме |
ходят |
обратно |
через |
поверх |
|||||||||
нении |
к |
электронам, |
приводят |
ность образца и тем самым |
||||||||||
уменьшают среднее |
число актов |
|||||||||||||
к поправкам на поглощение |
ха |
ионизации, |
приходящееся на |
|||||||||||
рактеристического |
излучения, |
|
один электрон. |
|
||||||||||
|
|
|
|
|
|
|||||||||
вполне |
удовлетворительно |
со |
|
|
|
|
|
|
||||||
гласующимся |
с опытом, однако не в состоянии точно вос |
произвести наблюдаемое энергетическое распределение об ратно рассеянных электронов. Наиболее удовлетворитель ный с физической точки зрения подход заключается в рас чете траекторий электронов методом Монте-Карло с помо щью вычислительных машин при использовании известных закономерностей индивидуальных актов рассеяния элект ронов. Грин [4] и Бишоп [5] применили этот метод для расчета характеристик процесса рассеяния, в том числе коэффициентов обратного рассеяния электронов и рас пределения возбуждения рентгеновского излучения по глубине образца; некоторые из этих характеристик могут быть измерены на опыте. Полученное при этом согласие с опытом оказывается все еще недостаточным для удовлет ворительного расчета поправки на атомный номер. Отча-
сти это связано с неточностью исходных данных расчета для сечений рассеяния, а также с тем, что даже самые большие современные вычислительные машины не позво ляют смоделировать весь процесс рассеяния достаточно точно. В любом случае использование метода Монте-Карло при практических анализах, разумеется, совершенно не практично. Наконец, даже в этом, наиболее точном сейчас методе неизбежно приходится использовать непрерыв ную математическую функцию для описания потерь энер гии электронами при их движении в образце, тогда как
вдействительности потери энергии являются дискретными
ихарактеризуются значительным разбросом для электро нов, прошедших одинаковый путь в образце.
Располагая подходящим соотношением для потери энер гии на единицу пробега электронов и внося его в выраже ние для Q, можно без труда и независимо от каких бы то ни было моделей или допущений вычислить интенсивность рентгеновского излучения, возбужденного теми из падаю щих электронов, которые остановились в образце. Такой подход игнорирует статистический характер потерь энер гии при неупругом рассеянии, но он ничуть не хуже ука занного выше модельного подхода, главное значение ко торого заключается в даваемой им информации об обрат ном рассеянии электронов. Вместе с тем на сегодняшний день собран большой экспериментальный материал по обратному рассеянию, и на практике предпочтительнее использовать именно этот материал, а не результаты, полученные с помощью теоретических моделей, хотя по следние ценны, например, в предсказании характеристик обратного рассеяния от образцов сложного состава.
Способность вещества замедлять падающие на него электроны обычно характеризуется удельной тормозной способностью S, имеющей следующее определение:
где р — плотность вещества, х — путь, пройденный элект роном, olE — среднее изменение энергии при прохожде нии им пути dx. Определенная таким образом тормозная
способность |
оказывается |
приблизительно постоянной |
для каждого |
химического |
элемента безотносительно к |
его физическому или химическому состоянию. Из аддитив ного характера потерь энергии вытекает, что тормозная
способность сложного образца |
АВ... |
равна |
|
|
||||||
|
|
|
SAB |
= 2 С А . |
|
|
|
(4) |
||
|
|
|
|
г |
|
|
|
|
|
|
где |
с; — весовая |
концентрация |
г-го |
элемента, |
S T — его |
|||||
тормозная способность. |
|
|
|
|
|
|
|
|||
|
Наилучшее выражение для тормозной способности бы |
|||||||||
ло |
получено Бете |
[6]: |
|
1 . |
1,1667? |
|
/ с . |
|||
|
с |
= |
. |
Z |
, |
|||||
|
S |
const • — |
— |
In |
j |
|
(5) |
|||
где |
Z ж А — соответственно |
атомный |
номер и |
атомный |
вес элемента. В применении к электронам более правиль ные значения S дает выбор постоянной 1,166 вместо часто встречающегося в литературе значения 2, которое при менимо скорее в случае протонов и других положительных ионов (Бете и Ашкин [7]).
Формула Бете была первоначально выведена для ато мов водорода, и для ее распространения на более тяжелые атомы потребовалось введение среднего потенциала иони зации J , который должен быть выражен в тех же едини цах, что и энергия электронов Е в формуле (5). Блох [8], исходя из статистической модели атома Томаса — Ферми, заключил, что / должно быть прямо пропорционально Z, так что JIZ = const, и что значение постоянной должно быть около 13,5 т. Некоторые экспериментальные исследо вания подтвердили по крайней мере приближенную спра ведливость этого заключения, хотя и область энергий, охваченная указанными опытами, в общем лежит значи тельно выше используемой при микроанализе; кроме того, согласие между разными экспериментами не было хоро шим. Уилсон [9] получил значение J/Z = 11,5 эв для алю миния; это значение широко используется в литературе, и его можно внести в выражение (5), предполагая, разумеет ся, что J/Z постоянно по всей периодической системе. Для некоторых задач такое предположение не приводит к су щественной ошибке, но мы позже увидим, что лучше ис пользовать несколько иное соотношение между / и Z.
Из определения эффективного сечения ионизации Q сле дует, что среднее число п ионизации і^-оболочки в расчете на один падающий электрон с энергией Е0 равно
N0pcA I' Q
n - — A - \ i m d E ^ |
( 6 ) |
где N0 — число |
Авогадро, |
СА — весовая концентрация |
|
элемента А, |
ЕК |
— энергия ионизации его ./t-уровня. Та |
|
же формула |
при подстановке |
в нее соответствующих зна |
|
чений Q и |
энергий возбуждения применима и к другим |
спектральным сериям. С учетом (3) выражение (6) прини мает вид
п = |
N u C * |
Г |
Q лі? |
( 7 ) |
А |
|
-%-dE. |
|
|
|
|
|
|
к
Величина п прямо пропорциональна интенсивности /Ка линин, измеряемой при анализе. Однако^ при расчете последней должно быть принято во внимание обратное рассеяние. Оно учитывается введением множителя R , рав ного отношению действитель но возбуждаемой интенсивно сти к той, которая возбуж далась бы, если бы все падаю щие электроны остались в об разце; из этого определения ясно, что R всегда меньше 1.
Таким образом, мы приходим к следующей формуле для ин тенсивности возбужденного в образце характеристического рентгеновского излучения:
SO дО 10 SO 60 70 80 fO ISO |
|
|
||
|
1 |
/ А = const • CAR ^ |
dE. |
|
Рис. 2. Доля |
обратно рассеянных |
|||
|
|
|||
электронов из |
пучка, падавшего |
|
|
|
нормально на поверхность образца, |
|
( о ) |
||
в зависимости |
от его атомного но- |
|
||
мера (Бишоп [10]). |
|
|
Вычисление фактора обратного рассеяния R . Величи на R связана с коэффициентом обратного рассеяния элект ронов г], равным отношению обратно рассеянного и падаю щего электронных потоков. На рис. 2 показана зависи мость т) (Z), измеренная Бишопом [10] на электронном пучке с энергией Е0 = 20 кэв, падавшем нормально на образцы. Аналогичные результаты получены другими ис следователями [11]. Можно считать, что п является плавной функцией Z и почти не зависит от Е0 в рассматри ваемой области энергий. Если бы все обратно рассеянные электроны покидали образец с начальной энергией Е0, то
1 — і? было бы равно г]. В действительности эти электроны имеют распределение по энергиям, меньшим Е0, так что 1 — R меньше т| на величину, которая зависит от этого распределения.
Вебетер с сотрудниками [12] показал, что R можно вы числить из энергетического распределения обратно рас сеянных электронов, располагая выражениями для Q и S,
Рис. 3. а) Энергетическое распределение обратно рассеянных электронов, нор мально падавших на поверхность углерода, меди, серебра и золота (Бишоп [10]); б) интегральная форма энергетического распределения; дополнительно приведены найденные интерполированием кривые для алюминия, титана и
плутония.
Альтернативные методы расчета R, предложенные Томас [13] и Томлином [14], недостаточно строги. Согласно Вебстеру R вычисляется по формуле
где W = Е/Е0, d\]/dW — энергетическое распределение обратно рассеянных электронов. Экспериментальные кри вые dr\ldW для углерода, меди, серебра и золота, получен ные Бишопом [10], приводятся на рис. 3, а. Из опыта най-
1?:5
дено, что от Е0 практически не зависит не только г), но и форма dx\/dW в функции W.
Энергетическое распределение обратно рассеянных электронов можно также записать в интегральной форме
г| (И7 ), дающей |
долю этих |
|
электронов с |
энергией, |
боль |
||||||||
шей WE0, |
что |
приводит |
к |
существенно |
более |
простой |
|||||||
формуле для R: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
Е„ |
|
|
|
|
|
|
|
1_д |
= |
J n(W)-§-d\VI |
jj |
-%-dE. |
|
|
(10) |
|||||
На рис. З, б приведены значения R, |
вычисленные |
по |
|||||||||||
соотношению (10) с |
использованием данных |
Бишопа для |
|||||||||||
|
|
|
|
|
|
|
т) (W), |
|
нормированных |
||||
|
|
|
|
|
|
|
на и (Z) |
из |
рис. 2, и с |
||||
|
|
|
|
|
|
|
использованием плавной |
||||||
|
|
|
|
|
|
|
интерполяции |
для эле |
|||||
|
|
|
|
|
|
|
ментов, |
|
располагаю |
||||
|
|
|
|
|
|
|
щихся |
между теми, для |
|||||
|
|
|
|
|
|
|
которых |
эксперимен |
|||||
|
|
|
|
|
|
|
тально |
найдены энерге |
|||||
|
|
|
|
|
|
|
тические |
распределе |
|||||
|
|
|
|
|
|
|
ния. Для S использова |
||||||
|
|
|
|
|
|
|
лось |
выражение |
(5) |
с |
|||
|
|
|
|
|
|
|
J/Z = |
11,5 зб, для Q — |
|||||
|
|
|
|
|
|
|
выражение, |
полученное |
|||||
|
|
|
|
|
|
|
Вебстером |
и |
др. [12]. |
||||
|
|
|
|
|
|
|
Изменение J/Z в преде |
||||||
|
|
|
|
|
|
|
лах погрешности, с ко |
||||||
|
|
|
|
|
|
|
торой |
известна эта |
ве |
||||
Рис. 4. Доля ионизации, вызванных |
в |
об-, |
личина, оказывает |
пре |
|||||||||
разце электронами, за вычетом потерь на |
небрежимо малое |
влия |
|||||||||||
обратное рассеяние. Вычислена по данным |
ние на значение R. Ана |
||||||||||||
Бишопа, .приведенным на |
рис. 2 |
и |
3, б. |
||||||||||
|
|
|
|
|
|
|
логично, зависимость R |
||||||
от Е0, возникающая |
из-за изменения S с энергией, также |
||||||||||||
весьма слаба, и ею можно |
|
пренебречь. |
|
В результате |
R |
||||||||
можно выразить в функции только Z и избытка |
напряже |
||||||||||||
ния U = |
Е0/Ец, |
чт |
сделано на рис. 4. |
|
|
|
|
|
Приведенные на этом графике кривые R (Z, U) очень близки к полученным Грином [15], который использовал аналогичный метод обработки более ранних данных по обратному рассеянию; они также хорошо согласуются с данными прямых экспериментальных измерений Дерья-
на [16]. Однако от них существенно отличаются кривые, опубликованные Томас [17]. Это различие можно припи сать, во-первых, неточному методу интегрирования, при нятому Томас, во-вторых, сомнительной точности исполь зованных ею данных по обратному рассеянию и, в-третьих, что наименьший атомный номер, для которого Томас вы числила R, был равен 29, а кривые экстраполировались к R — 1 при Z = 0 для всех U. Рис. 4, однако, свидетель ствует о перегибе кривых в этой области, что связано с об щей тенденцией изменения формы энергетических распре
делений обратно |
рассеянных электронов |
в |
зависимости |
|||
от Z. Следовательно, |
использованная Томас |
линейная |
эк |
|||
страполяция является неточной и привносит |
значительные |
|||||
ошибки в том диапазоне атомных номеров, |
в |
котором |
рас |
|||
полагаются наиболее |
распространенные |
элементы. |
|
|||
При введении |
поправки на атомный номер возникает |
|||||
проблема вычисления R для образцов сложного состава. |
||||||
Экспериментальные |
исследования г\ |
показывают, |
что |
весьма точные значения для сплавов и соединений дает среднее значение по весовым содержаниям компонентов:
ц =21 |
с ^ 1 - |
( И ) |
|
|
Поэтому разумно считать, что та же процедура будет при менима и к В, так что для состава А В
(12)
Теоретические расчеты Бишопа [18] и Брауна [3] нахо дятся в весьма хорошем согласии с этим допущением, хотя Бишоп и предположил, что лучше усреднять форму энер гетического распределения обратно рассеянных электро нов независимо от п. Этот вопрос требует дальнейшего изучения, а пока для практических целей можно приме нять выражение (12).
Более удобно выражать R в зависимости не от U, а от i/U: как отметил Бишоп [18], при данном Z величина/? изменяется почти линейно с 1/U, что существенно облегча ет интерполирование. Для практических целей удобно иметь табличное представление R. Соответствующие дан ные рис. 4 приведены в табл. 1, откуда для произвольных Z и MU величину R можно получить линейной интерпо ляцией с интерполяционной ошибкой менее 0,2%. Эту
Т а б л и ц а 1
Значения фактора обратного рассеяния R в зависимости
от |
1/17 и |
Z |
|
|
|
|
|
|
\ l |
/ U |
0,01 |
0,10 |
1 |
0,20 |
0,30 |
0,40 |
0,50 |
z |
|
|||||||
|
|
|
|
|
|
|
|
|
|
0 |
1,000 |
1,000 |
|
1,000 |
1,000 |
1,000 |
1,000 |
|
10 |
0,934 |
0,944 |
|
0,953 |
0,961 |
0,968 |
0,975 |
|
20 |
0,856 |
0,873 |
|
0,888 |
0,903 |
0,917 |
0,933 |
|
30 |
0,786 |
0,808 |
|
0,828 |
0,847 |
0,867 |
0,888 |
|
40 |
0,735 |
0,760 |
|
0,782 |
0,804 |
0,827 |
0,851 |
|
50 |
0,693 |
0,718 |
|
0,741 |
0,764 |
0,789 |
0,817 |
|
70 |
0,662 |
0,688 |
|
0,713 |
0,737 |
0,764 |
0,793 |
|
90 |
0,635 |
0,663 |
|
0,687 |
0,713 |
0,740 |
0,770 |
|
80 |
0,611 |
0,639 |
|
0,665 |
0,691 |
0,718 |
0,750 |
|
60 |
0,592 |
0,613 |
|
0,639 |
0,665 |
0,695 |
0,730 |
|
99 |
0,578 |
0,606 |
|
0,634 |
0,661 |
0,691 |
0,725 |
|
0,60 |
0,70 |
0,80 |
0,50 |
1,00 |
0 |
1,000 |
1,000 |
1,000 |
1,000 |
1,000 |
10 |
0,981 |
0,988 |
0,993 |
0,997 |
1,000 |
20 |
0,948 |
0,963 |
0,977 |
0,990 |
1,000 |
30 |
0,911 |
0,935 |
0,959 |
0,981 |
1,000 |
40 |
0,878 |
0,907 |
0,938 |
0,970 |
1,000 |
50 |
0,847 |
0,881 |
0,919 |
0,959 |
1,000 |
60 |
0,829 |
0,862 |
0,904 |
0,950 |
1,000 |
70 |
0,805 |
0,844 |
0,889 |
0,941 |
1,000 |
80 |
0,785 |
0,826 |
0,874 |
0,932 |
1,000 |
90 |
0,767 |
0,811 |
0,862 |
0,924 |
1,000 |
99 |
0,763 |
0,806 |
0,858 |
0,921 |
1,000 |
таблицу можно легко ввести в машинную программу рас чета поправок к данным микроанализа.
Вычисление тормозной способности S. Из выражения
8)следует, что отношение кл интенсивностей харак
теристического |
рентгеновского |
излучения от образца |
|
и эталона, в |
пренебрежении |
его |
поглощением и |