Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Физические основы рентгеноспектрального локального анализа

..pdf
Скачиваний:
17
Добавлен:
25.10.2023
Размер:
15.53 Mб
Скачать

ПРОГРЕСС В ВЫЧИСЛЕНИЙ ТОРМОЗНОЙ СПОСОБНОСТИ И ФАКТОРА ОБРАТНОГО РАССЕЯНИЯ

П. Дапкамб, С. Рид

Введение. Кастен [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)следует, что отношение кл интенсивностей харак­

теристического

рентгеновского

излучения от образца

и эталона, в

пренебрежении

его

поглощением и

Соседние файлы в папке книги из ГПНТБ