книги из ГПНТБ / Электронно-зондовый микроанализ [сборник]
..pdfЭлектронно-зондовый микроанализ |
129 |
поглощения для углерода в работе [109] отличаются на 20% отдан ных Генриха; в работах [103, 108] выявлен ряд расхождений с дан ными Генриха для поглощения Al, Zr и Nb, но данные Генриха для А1 ближе всего к значениям ц/р [109]. Келли сопоставил полу ченные им значения с данными [103, 108] и получил хорошее совпа дение результатов; в частности, хорошее согласие (в противополож ность данным Генриха) получено для А1. Детальное полное сравнение результатов Генриха и Келли не проводилось. Можно предпо ложить, что различия будут сравнительно незначительными и любая совокупность значений могла бы быть использована как обширный и в общем случае надежный источник. Обзорная работа Тайзена и других авторов не сравнивалась с работами Генриха и Келли, но, по-видимому, эта работа превосходит их по точности измерений.
Следует упомянуть также работу [112] по расчету коэффициента поглощения, так как этот метод был использован несколькими авто рами. Полностью он описан в работе [48]; метод полезен при опре делении приближенного значения данного коэффициента и уступа ет по точности данным Генриха и Келли.
Имеется очень мало данных по коэффициентам поглощения для 12 Â . Наиболее обширными и, вероятно, наиболее надежными являются обзоры Хенке и сотр. [113, 114]. В работе [113] приведе ны коэффициенты для длин волн вплоть до 44 Â , что соответствует СКя-излучению; рассмотрены поглотители с атомным номером не более 33, постепенно убывающим до 16 при %> 8 Â . В работе [114]
диапазон длин волн расширен до 200 Â , но рассмотрены поглотите ли с Z <1 18. Поскольку легкие элементы часто встречаются в ком
бинации с более тяжелыми, необходимы измерения для более тяже лых поглотителей. Успешной в этом направлении была работа [115], где даны значения для поглотителей с 6 Z 82 для К а-линий
С , N и О.
Массовый коэффициент поглощения для многокомпонентной мишени может быть получен из коэффициентов для составляющих
элементов по уравнению |
|
|
(Ы?)а - sІ |
И ([х/р)j . |
(50) |
|
В . Обзор теоретически обоснованных уравнений для поправок на атомный номер и поглощение
В этом разделе рассмотрены уравнения из предшествующих разделов, практически применяемые для внесения поправок, вклю чая табличные данные, программы для электронно-вычислительных машин и различные упрощения. Точность этих уравнений и калиб ровочных формул оценена в разд. III.
5—584
130 |
П . Мартин, Д . П ул |
1. Метод П у л а и Томас [11]
Зависимость между интенсивностью рентгеновского излучения и концентрацией, найденная Пулом и Томас, уступает по точности методу Томас [4, 5], где используется более надежный и менее под верженный ограничениям фактор обратного рассеяния.
2. Метод Томас [4, 5]
В этом методе поправка на атомный номер по Томас комбиниру ется с модифицированной поправкой на поглощение по Филиберу [69—71]. Концентрация выражается через отношение интенсивнос тей излучения
(51)
где a t = S f / R f . Значение S f , соответствующее Е = (£0+ £ с)/2, находится по графику зависимости 5 от £ для атомного номера Z,
определенного из рис. 3. Последний построен на основе таблиц 123, 24] с использованием закона Бете [уравнение (7)]. По мнению Генриха [91], использование функции S A IZ , изменяющейся более плавно по сравнению с изменением S от Z , может повысить точность
интерполяции значений для элементов, не включенных в таблицы Нельмса [23, 24] (в области легких элементов может быть полезен
обзор 1116]). Значение R t взято из графика зависимости R от Z
в работе Томас, основанной на данных обратного рассеяния [58] при соответствующем значении E 0IE f . Функция і л І'/а) рассчитана
по уравнению (37); |
f£{y*) также рассчитана по |
уравнению (37), |
но с подстановкой у' |
[см. уравнение (43)] вместо х и h — h а • |
|
Так как для вычисления поправок требуется |
знание состава |
образца, сначала необходимо использовать в правой части уравне ния (51) приближенное значение са ■ Значение са , получаемое в
левой части уравнения, вновь подставляется в правую часть, и процедура повторяется. После нескольких итераций значения са
сходятся и принимают постоянное значение; Томас [5] привела пример такого вычисления. Обычно метод итераций необходим во всех случаях, когда поправки зависят от концентрации.
Несмотря на то что расчет по уравнению (51) может быть прове ден вручную, он становится трудоемким при большом числе ите раций или в случае анализа многокомпонентных сплавов. Разумеет ся, значительно быстрее расчет можно выполнить на ЭВМ . Шоу [137] разработал программу ФОРТРАН IV для расчетов на ЭВМ поправки по методу Томас и поправки на флуоресценцию по Риду [118] (см. разд. I, Д , 1). С помощью этой программы производится расчет поправок для системы, содержащей 8 элементов, при условии,
Электронно-зондовый микроанализ |
131 |
что измерены отношения интенсивностей для всех элементов, кро ме одного. Все параметры, необходимые для поправок, рассчиты ваются машиной на основе использования эмпирических уравне ний для получения согласия с табличными значениями а по Фили
беру, значениями 5 по Нельме, кривыми Томас для а и с данными обзора Генриха [102] по коэффициентам массового поглощения. Основные данные, например длины волн краев полос поглощения рентгеновских лучей, длины волн излучения главных линий К, L и М и их потенциалы возбуждения, содержатся в программе или
выдаются машиной и их не требуется включать во входные данные. В дальнейшем программа была расширена для предварительного отбора входных данных [119]; в настоящее время используется ме
тод Данкамба и да Каза [12, |
13] вместо метода Томас. Другая гиб |
||
кая программа Ф ОРТРАН , |
которая может включать |
поправку |
|
на атомный номер Томас, составлена Броуном [120]. |
|
||
При отсутствии ЭВМ можно обратиться к составленным Колби |
|||
[121, 122] таблицам значений / (х), рассчитанным по |
уравнению |
||
Филибера, принимая х = х |
для образца и используя данные, соот |
||
ветствующие h = ііа |
для эталона и образца. |
|
|
3. |
П ервы й |
метод П ейсскер а [123] |
|
Метод представляет, по существу, модификацию метода Томас [4, 5], в котором используется поправка на поглощение Филибера [69— 71] со значением ас по Данкамбу и Шилдсу [68] [уравнение (39)] вместо а. Значение R получено из кривых, подобных кривым
в работе [60], основанной на данных [59]. Следует подчеркнуть, что в этом методе не предусматривается применяемая Томас замена
X на х ' , а Ха вычисляется из суммы 2 ctxf\ значение h вычисляет-
І
С Я ИЗ с у м м ы E C j / i j .
4. Второй метод П ейсскера [123]
Разница между этим методом и методом, описанным выше, со
стоит в том, что отношение аА /af |
получают из выведенного Пейс- |
|
скером эмпирического уравнения для бинарной системы AB в виде |
||
Л- |
0,26 |
0,075 ln ( Z A / ZB ) |
Е А |
(52) |
|
А |
|
|
|
с 0 |
|
5. |
Метод Д анка м ба и да К азы [12, 13] |
1) Зависимость между концентрацией и отношением интенсив ностей может быть представлена по Данкамбу и да Казу следующим
5 *
132 П . Мартин, Д . П ул
уравнением:
|
/ л |
r a |
S c ; s f |
f i ( - d ) |
|
|
Сл~ |
' S* * S c ; Rf |
t i W |
} |
|
Значение R |
определяется по экспериментальным данным Бишопа |
||||
[61], S — из |
закона Бете [19, |
20, 22] |
[уравнение (7)] при Е |
= |
|
= {Е0-\-Ес)/2 и расчете J |
из уравнения |
(8); f (%) вычисляется |
по |
формуле Филибера [69*—71] [уравнение (37)], модифицированной Данкамбом и Шилдсом [68] с заменой а на ас [уравнение (39)].
Расчет са производится методом итераций на ЭВМ . Программа
Шоу [119] |
для ЭВМ |
составлена на основе |
метода |
Данкамба и |
да Каза и, |
несмотря на то что она специально разработана для нового |
|||
варианта (см. ниже), |
необходимо, очевидно, |
лишь |
незначительно |
изменить ее для применения с первым методом Пейсскера [123]. Программа Броуна [120] также может быть приспособлена для при менения с этим методом. Хобби и Вуд [124] составили программу для модифицированной поправки на поглощение Филибера на язы ках АЛ ГОЛ , автокод АТЛ АС и в варианте АТЛ АС для автокода М ЕРК УРИ Й ; программа включает также поправку на характерис тическую флуоресценцию Кастена [2] (см. разд. I, Д).
Колби [125] дополнил таблицы значений / (у) (рассчитанные по
уравнению Филибера), принимая во внимание модификацию Дан камба и Шилдса (<ус). Аналогичные таблицы составлены Адлером
иГольдстейном [126].
2)Позднее было использовано значение ас по Генриху [97]
[уравнение (40)], которое дает, по-видимому, более высокую точ ность [13].
6. П ервы й метод Ш прингера [15]
Простая формула Шпрингера весьма сходна с уравнением Дан камба и да Каза [12], но вместо S используется S' = SE:
Ra |
(54) |
|
X ici Rf f a ( l a ) |
Значение R определено Шпрингером [60] из данных по обратному
рассеянию [59]; S — из закона Бете [19, 20, 22] в форме, незначи тельно отличающейся от уравнения (7), для Е = Е 0\ f (у) найдено
по формуле Филибера [69—71] [уравнение (37)]. В дополнение к таб лицам Колби [121, 122] Хечгессон [127] составил программу для ЭВМ в случае формулы Филибера, а Броун [128] — для формулы Филибера с поправкой на флуоресценцию по Кастену (см. разд. I, Д). Эти программы могут быть полезны при использовании формулы: Шпрингера.
Электронно-зондовый микроанализ |
133 |
7.Метод Филибера и Тиксье [9, 10]; второй метод Шпрингера [8]
Вдополнение к расчету / {у) Филибер совместно с Тиксье рас
смотрел отдельно поправку на |
атомный номер. |
Выражение для |
||
са , выведенное в этой работе и объединенное с / (у), |
имеет вид |
|||
сА |
|
f U l i ) |
(55) |
|
Rn |
/aAU ) ’ |
|||
|
|
где S получено из уравнения (18). Филибер и Тиксье [9] полагают, что формула Филибера [69—71] (с модификацией Данкамба и Шилд са или без нее) может быть использована для определения / (у). Значение R определяется из кривых зависимости R от Z, построен
ных Данкамбом и сотр. [12, 13] по данным Бишопа и приведенных Филибером и Тиксье.
Шпрингер в более поздней работе приходит, по существу, к такой же формуле, что и уравнение (55), с определением S по урав нению (18). Он рекомендует применение формулы Филибера, моди фицированной Данкамбом и Шилдсом для f (у), и значений R по
Бишопу [62].
8. Первый метод Тайзена [98]; первый метод Тонга
При первом рассмотрении влияния атомного номера и погло щения Тайзен получил уравнение для F (у), подобное уравнению
Филибера [уравнение (36)], но |
со значением h = Ъ ,Ь А И %, а не |
|||
1,2/1/Z2. Таким образом, |
|
|
|
|
сА |
/А |
Р І Ы ) |
(56) |
|
? аА { г аА ) |
||||
|
' А |
|
||
Идентичное уравнение для F |
(у) |
приписывается Тонгу. |
|
9.Второй метод Тайзена [99—101]
Позднее Тайзен вывел уравнение для сА в виде
сА |
r a |
. Га Ы ) |
(57) |
|
‘1А < |
П { г Аа ) ' |
|||
|
|
Значение R получают из кривых зависимости R от Z по данным Томас [4, 5]. Значение R A находится путем определения среднего
атомного номера по уравнению
Z = 2 Ci Z;. |
(58) |
134 |
П . Мартин, Д . Пул |
Значение / (х) вычисляется по модифицированному уравнению Тайзена — Тонга (разд. I, В, 7), в котором а заменено на ас [уравне
ние (41)] и h вычисляется из уравнения (42). В этом методе 5 — до минирующий фактор в поправке на атомный номер — отсутствует. Тайзен составил таблицы значений различных величин, входящих в уравнение (57), однако часто требуется весьма сложная интерпо ляция данных. Следует отметить, что приведенные Тайзеном таб лицы основных данных имеют ряд ошибок.
10. Второй метод Тонга [129]
Бимен [130] сообщает о разработанном Тонгом втором методе. Однако во время работы над данным обзором авторы не располага ли какими-либо данными относительно этого метода.
11. Метод Ильина и Лосевой [14]
Основное различие между этим методом и приведенными выше заключается в поправке на поглощение fx. Зависимость между сд и /л //л имеет вид
(59)
І
Значение R вычисляется из уравнения (16), а 5 — по закону Бете
[19—21] [уравнение (7)]; второй подстрочный индекс «Л» является лишним для R и 5 при расчете по методу Ильина и Лосевой, так
как зависимость этих величин от критического потенциала возбуж дения излучения не учитывается. Для значений у 1000 поправка fx определяется из уравнения (44).
12.Метод Биркса [82, 55]
Первое приближение по Бирксу имеет вид
СА = !L |
Ң Ы ) |
(60) |
|
Ѵ Ы )
F (х) найдено по универсальным кривым Биркса при соответствую щем значении Е 0. Для большей точности Бирке [82] принял, что
сА |
F a ( I x P z Ü ) |
(61) |
|
Ft dxPzt f ) |
|||
|
|
Электронно-зондовый микроанализ |
135 |
Значение P z получено из таблиц, основанных на данных Филибера
[69—71]. Полагают, что [j P z\a = 2 И [цРг^> хотя у Биркса нет
І
такой формулировки. В любом случае существуют сомнения отно сительно применения параметра P z. Бирке [131] впоследствии со гласился, что использование величины P z неоправданно.
13. Метод Арчарда и Мулви [18, 53]
Расчеты по методу Арчарда и Мулви позволили получить ряд кривых, учитывающих поправку на атомный номер для ограничен ного числа бинарных систем, а также поправку на поглощение для ряда элементов. Пул и Томас [132] путем экстраполирования и интерполирования данных, полученных при помощи кривых Ар чарда и Мулви, распространили предложенные последними поправ ки на некоторые другие системы, однако метод не внушает полного доверия. Метод Арчарда и Мулви должен рассматриваться как имеющий только ограниченное применение. Для оценки точности этого метода Пул и Томас сравнили его с другими методами.
14. Метод Белка [133]
Белк использовал поправку на атомный номер, предложенную Пулом и Томас [11], а также поправку на поглощение по Филиберу [69—71] и, упростив их, получил следующее уравнение для Са -
Знаменатель в выражениях для фУр )а и Z отражает то обстоя тельство, что сумма значений Ра11\ обычно не равна 1; использова
ние отношения lallt вместо ct исключает необходимость в итера
ции. Несмотря на то что это уравнение — результат серьезного упрощения, полученное выражение можно легко применять для мно гокомпонентных систем даже при отсутствии ЭВМ . Как показано в разд. III, сравнение результатов, рассчитанных по данному ме
136 |
П .- Мартин, Д . Пул |
|
толу и с помощью более сложных выражений, в ряде случаев |
дает |
|
одинаковую точность. |
Однако применение его ограничено; |
метод |
обладает достаточной точностью тогда, когда эффект поглощения мал. В работе [133] установлены границы, за пределами которых точность понижается.
15. Первый метод Смита [134]
Метод включает большой объем уже рассмотренных теоретиче ских положений. Вместе с тем для S принимается эмпирическое уравнение, основанное на предположениях Тонга и Рида:
S = [4,7 + 0,9 In ( Е 0- Е.?) - ln Z] 4 “ • |
(63) |
Зависимость между интенсивностью и концентрацией имеет вид
с |
Kg |
, f i ( l i ) |
(64) |
|
* A |
fa ( і а ) |
|||
|
|
Значение R получено из работы Грина [28]; значения / (у) для Mg,
Al и Si определяются многочленом, коэффициенты для которого приводит Смит. В случае элементов с Z > 14 Смит предлагает применять для f (у) формулы Филибера [69—71] или Биркса [82, 84, 85].
16.Метод Боровского и Рыдника [31]
Врезультате совместного применения поправок по уравнениям
(47а) и (476) зависимость между отношением интенсивностей и концентрацией но Боровскому и Рыднику принимает следующий вид:
с . |
. Ja |
И + |
2,0a[ l - ( P o T |
l d |
|
Л |
|
b + |
|||
' И - |
fe* (£oP [1 - |
( < Л f — 'У |
cosec Ѳ |
||
у ___________________________________ \ Р |
/а__________ |
||||
1 + |
*2 (Е0)п [1 - |
( и £ у п\ |
|
cosec О |
|
Относительно величин /, |
п и k 2 см. разд. |
I, |
Б, 3. |
Электронно-зондовый микроанализ |
137 |
Г. Некоторые специфические особенности анализа легких элементов
Количественный анализ легких элементов (Z < 12) представляет во многих отношениях экстремальный случай, и допущения и при ближения, использованные при выводе поправок для средних и тя желых элементов, могут оказаться недействительными. С другой стороны, возможны некоторые упрощения. Варианты существую щей теории были рассмотрены Данкамбом и Мелфордом [135].
Минимальная энергия пучка, обеспечиваемая в типичном микро анализаторе, составляет ~ 5 кэВ. Это означает, что при анализе легких элементов отношение перенапряжения (U0) всегда высокое.
Например, для углерода это отношение составляет ~20 при 5 кэВ, тогда как для элемента со средним атомным номером, таким, как медь, отношение перенапряжения составляет всего лишь 5 при 45 кэВ, т. е. при напряжении, близком к максимальному для мик роанализатора. Работа при высоком перенапряжении предъявляет более жесткие требования к выражениям для сечения ионизации, которые должны быть точными в значительно более широком диапа зоне энергий. Предположение о постоянном значении ф, принятое, например, Филибером [69—71], в этих условиях будет, по-видимо му, менее приемлемым.
По мере торможения электрона и понижения отношения пере напряжения эффективность порождаемого им рентгеновского излу чения возрастает по крайней мере до тех пор, пока U не понизится
до ~ 3 (см. рис. 5). Высокое перенапряжение сильнее влияет на распределение по глубине порождаемого рентгеновского излучения, увеличивая его долю в результате генерации от более глубоких слоев, чем в случае низкого перенапряжения. Данкамб и Мелфорд сравнивают метод Монте-Карло в работе Бишопа [77] для распреде ления по глубине С/(«-излучения в Fe3C при Е 0 = 20 кэВ, которо
му присуща отмеченная выше тенденция, с экспоненциальной мо делью, предложенной Филибером, где такая тенденция отсутствует. Данкамб и Мелфорд придают известное значение результатам рас четов Бишопа, так как полученные Бишопом значения Ф (рг) у поверхности ф0 для C/C« в Fe3C при 5, 10, 20 и 30 кэВ согласу ются со значениями, полученными ими экспериментально, с точ ностью до ±10% . Эти авторы приходят к выводу, что теоретическая модель Филибера, находящаяся в хорошем согласии с эксперимен тальными данными Кастена [3, 79] при отношении перенапряжения ~ 3 , не может дать надежные значения / (х) при более высоких пе ренапряжениях.
Возможно, что поправка Боровского и Рыдника [31] приемлема для легких элементов. Как показали эти авторы, для высокого отношения перенапряжения (UQ > 1) уравнение (46) упрощается
138 П . Мартин, Д . Пул
(если рассматривать только изменение с U0) и принимает вид
I А = const • |
( A^n~ l ln U A |
(66) |
|
l + |
fe2 (£ 0) « ^ ^ ( ^ ) " c o s e c 6 |
Эта зависимость хорошо согласуется с измерениями Кэмпбелла
1136] для /(-линий С, В и Be [31].
Помимо поправки на поглощение Боровского и Рыдника, дан ные, которые могут быть использованы при количественном анализе легких элементов, содержатся в экспериментальных работах Грина [28, 73], в работе [94], продолжающей расчеты Арчарда и Мулви, а также в расчетах методом Монте-Карло по Бишопу; однако они, за малым исключением, применимы только в случае углерода.
Р и с . |
14. |
Расчетные |
кри |
|
вые |
поверхностной |
иони |
||
зации |
фо |
как |
функции |
|
атомного номера |
7. [137]. |
|||
------------ расчет обратного |
рассея |
|||
ния (по Данкамбу); ^ |
эксперимен |
тальные значения для пленки угле рода (по Мелфорду); % расчеты
методом Монте-Карло (по Бишо пу); Н эксперименты с мечеными слоями (по Кастену).
Поверхностная ионизация <p{Uü,Z)
Отсутствие достаточно надежных данных по / (у) Для легких
элементов обусловливает интерес, проявляемый к так называемой модели тонкой пленки, которая была предложена Данкамбом и Мелфордом [135] и развита Данкамбом [137]. Это упрощение, в частности, дает возможность учитывать влияние атомного номера и поглощения при помощи одного простого уравнения. Модель основана на том, что при анализе легких элементов и образец и эталон сильно поглощают. При очень сильном поглощении из образца^выходит только то рентгеновское излучение, которое порож дено в тонком поверхностном слое. В этом случае эффективная толщина эмитирующего слоя Д (pz) обратно пропорциональна у.
Учитывая, что интенсивность испускаемого излучения составляет Ф0Д (pz), имеем
л |
_ _ |
,„л |
„л |
|
Га |
?0а |
У-А |
(67) |
|
ІА |
~ °А |
тА |
' ѵл |
|
1А |
|
'РОЛ |
У-а |
|