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

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

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

где Лі интенсивность, излучаемая элементарным сло­ ем алюминия в идеальном случае — в отсутствие рассея­ ния электронов, при поверхностной плотности слоя, рав­ ной единице, и при энергии электронов Е0. В кольцевой

Рис. 2. Схема к расчету поправки на флуоресценцию. 1 — алюминий, 2 — магний.

объем, вырезаемый в слое магния массовой толщиной р х £ , расположенном на глубине p z b двумя конусами с углами полураствора г|з и і|з + dip, из источника S доходит интен­ сивность

а2!^ = —-— /А іср (pz) sin гр • ехр І

dipd(pz).

 

C OS V|)

Интенсивность излучения, прошедшего в глубь сквозь этот

объем,

равна

 

 

 

 

d2l\ =

Hi

(P2 i -

Pz )

X

 

-^-ІА\Ц> (pz) sin ip • ехр

COSVJ)

 

 

 

 

 

 

 

 

X

ехр

со 5 г|з

dip d (pz),

 

 

 

 

где (ІА! И u^f массовые коэффициенты поглощения Ка- излучения алюминия соответственно в алюминии и магюга. Интенсивность же излучения, поглощенного в слое

магния,

равна

, А1 ,

 

 

 

Ліф (pz) sin ip • ехр

 

 

d4\ =

cos \p

X

 

 

 

 

 

 

X 1 — exp

СОЗ ip

dip d(pz).

 

 

 

 

Таким образом, полная поглощенная интенсивность от

очи ика во всем слое

магния

составляет

 

dh =

/ А 1 Ф

(pz) d (pz)

\ {ехр у

 

 

 

 

о

 

 

 

 

— ехр '

cos i|>

sin l])tttj\

Обозначим

 

 

 

 

 

 

 

 

а

- (ХАГ(pz! — pz),

р (J-AFPIS,

 

 

я 2

 

 

 

 

О

Тогда интенсивность вторичного излучения в телесном

угле Q

равна

1 - 1 / А І Ф

(pz) d (pz) { / (a) - / (a + p)},

dl) =

±-v>K (Mg)

где ojjf (Mg) — выход флуоресценции для /іГ-излучения магния, а г — отношение массовых коэффициентов погло­ щения по обе стороны от ІГ-края поглощения магния

( г > 1 ) .

Область образца с К ^ дает вклад в излучение

по

о

Аналогично, вклад для области z1 < z < оо равен

ПО

где

 

 

Т |АА1 (('2 — f>3

 

Интегралы

// и 7/

были рассчитаны на

электронной вы­

числительной машине.

 

 

Вклад флуоресцентного излучения магния по отноше­

нию к его

первичному

характеристическому излучению

равен

 

 

 

 

 

 

1

г _ і / ° ,

+

 

AcP / =

- l ;

o H M g ) V - # - ^ '

Б принятой для расчета области ускоряющих напри* жений для определения отношения /лі//ме мы^выбрали сечение ионизации в форме, предложенной Росселандом. Полученные при этом результаты отличаются от найден­ ных при выборе сечения ионизации в виде (l/U) In U не более, чем на 5 %. Вычисление указанного отношения, кро­ ме того, требует априорного знания функций cp (pz)j в первом приближении можно взять экспериментальные значения этих функций и затем улучшать результат мето­ дом последовательных приближений. Таким образом, мы не можем учесть вклад вторичного возбуждения при нор­ мировке получаемых кривых <p (pz) с мечеными слоями различной толщины, и остается принять, что распределе­ ние по глубине вторичной флуоресценции для тонкого слоя мишени не зависит от толщины этого слоя. Это сооб­ ражение имеет то обоснование, что функция ф (pz) для вто­ ричного возбуждения меняется значительно более медлен­ но, чем соответствующая функция для первичного воз­ буждения излучения.

Чтобы получить возможно более точное значение ф (0), мы измерили эту величину отдельно. Она представляет собой сумму вкладов от нерассеянного первичного пучка и от потока обратно рассеянных электронов из мишени. Соответствующее измерение интенсивности проводилось от слоя магния, расположенного на поверхности мишени. Экспериментальные значения этой интенсивности следует исправить на возможное рассеяние электронов в изоли­ рованном меченом слое, служившем для эталонирования, и на вклад вторичного возбуждения в излучение магние­ вого слоя на алюминиевой мишени. Соответствующие зна­ чения этой интенсивности приведены в таблице.

Ускоряющее

29

25

20

15

І0

напряжение,

К1

Ф(0) 1,52 1,51 1,45 1,38 1,26

Точность полученных результатов мы оцениваем в 5% для 29, 25 и 20 кви 10% для 15 и 10 ке; при низких напря­ жениях поправки труднее оценить. Величина ф (0) в дей­ ствительности представляет только теоретический интерес: для построения кривых поглощения излучения }(%)

Рис. 5. Кривые

lg / (х)> при 29 кв. 1 — А1 с меченым слоем Си И]; 2 - А *

 

с меченым слоем Mg (настоящая работа).

значения ф (pz) необходимо знать лишь с точностью до постоянного коэффициента пропорциональности.

На рис. 3 представлены экспериментальные кривые для различных ускоряющих напряжений на электронном зонде. На этом графике показаны также кривые распреде­ лений ф (pz), полученных Бишопом [2] методом МонтеКарло для алюминиевой мишени.

Применения. Определенные экспериментально функ­ ции ф (pz) мы использовали для построения кривых по­ правки на поглощение для алюминиевого излучения.

Пусть

j"

cp(pz)exp(— %pz)d(pz)

fix) - F(0) M

-

 

І ф (pz)rf(pz)

 

0

представляет собой отношение интенсивностей испущен­ ного излучения при наличии и отсутствии его поглощения. Оба интеграла здесь мы вычислили графически.

На рис. 4 приведено семейство кривых lg / (%), получен­ ных также и другими авторами: Бишопом [2] на осно­ ве вычисленной им функции ф (pz), Грином [3] экспери­

ментально (отметим, что он использовал значение и.д{

=

= 410 вместо более точного ц.д| = 380), Данкамбом

и

Шилдс [4], модифицировавшими формулу Филибера; так­ же приведены результаты расчета по формуле Филибера при 15 кв.

На рис. 5 представлены кривые lg / (%) для алюминия при 29 кв, найденные с помощью графиков распределения интенсивности по глубине с мечеными слоями меди [1] и магния. Первую из этих кривых следует использовать при анализе следов меди (или другого элемента со средним атомным номером) в алюминии, а вторую, например, для анализа сплавов A l — Mg .

ЛИТЕРАТУРА

 

 

 

 

 

 

 

 

1.

C a s t a i n g

R.,

D e s c a m p s

J . ,

Journal

de

Physique,

2.

p.

304 (1955).

 

Частное

сообщение.

 

 

 

B i s h o p

H . E.,

 

 

 

3.

G r e e M . , These,

Universale de Cambridge (1962).

 

см. также

4.

D u n c u m b

P.,

S h i e l d

P., Частное сообщение,

 

Meeting

of

the Electrochemical Society,

October

1964,

Washing­

5.

ton

D. C.

 

J . , Metaux,

Corros.

40,

216, 466

(1964).

P l i i l i

b e r t

ВЫЧИСЛЕНИЕ ПРОНИКНОВЕНИЯ ЭЛЕКТРОНОВ И ВОЗБУЖДЕНИЯ РЕНТГЕНОВСКОГО ИЗЛУЧЕНИЯ В ТОЛСТЫХ МИШЕНЯХ

Г. В и won

Введение. Точность количественного локального рент­ геноспектрального анализа нередко ограничивается нели­ нейностью соотношения между весовой концентрацией элеметна в образце и эффективностью возбуждения его характеристического рентгеновского излучения. Для раз­ работки метода учета этого «эффекта атомного номера» необходимо располагать детальным знанием рассеяния электронов в образце. Однако это рассеяние является слож­ ным процессом, и использование транспортного уравне­ ния для его описания имеет ограниченную ценность, по­ скольку решение этого уравнения возможно лишь в специ­ альном случае бесконечной мишени [1]. Ранее с ограни­ ченным успехом использовались теории, базирующиеся на упрощенных моделях проникновения электронов в ве­ щество. В последнее время развитие цифровых вычисли­ тельных машин позволило обратиться к другим, более перспективным методам решения этой задачи. В методе Монте-Карло ЦВМ моделирует траектории отдельных элек­ тронов в мишени. Грин [2] выполнил такой расчет для меди, основываясь на результатах измерения рассеяния в тонких пленках. Однако его подход по-прежнему не свободен от некоторых недостатков; мы разработали более гибкий метод, использующий только теоретические данные по рассеянию электронов. Краткий очерк метода расчета приводится в настоящей работе; более детальное его описание находится в печати * ) .

Описание расчета. Каждая электронная траектория разбивалась на 25 шагов. Рассеяние на каждом этапе вы­ числялось на основе теории многократного рассеяния Гаудсмита — Саундерсона [4]. Сечение однократного рас­ сеяния о* (•&) бралось в виде экранированного резерфор-

*) Н. Е. B i s h o p , Proc. Phys. Soc. 85, 855 (1965).. (Прим.'перев.)

довского сечения [5], умноженного на (Z + 1)

Z для учета

процессов неупругих столкновений [6]:

 

 

Z(Z

+ l)

1

 

(1)

а (в) = 5,2-10"2 1

(1 +

2а — cos fl)2 '

где а == 3 , 4 - Ю - 8 Zl'>IE, Е

— энергия электронов в

кэв,

Z — атомный номер рассеивателя,

•& — угол

рассеяния.

Потери энергии описывается уравнением Бете для непре­

рывных потерь

[7]:

dE

= 7,85 -10і ~ 4 - In

d(ps)

где s — траекторный пробег электрона, А — атомный вес, е — основание натуральных логарифмов. Значения средне­ го потенциала ионизации / взяты согласно Нелмс [8], что­ бы иметь возможность сравнить наши результаты с по­ лученными Спенсером [9].

Число траекторий, входивших в каждый расчет, со­ ставляло 5000. Это дает статистическую ошибку около + 0 , 0 1 в значениях коэффициента обратного рассеяния. Об­ щая точность расчета проверялась сравнением средней и среднеквадратичной глубин проникновения электронов, полученных методом Монте-Карло и вычисленных из урав­ нений Льюиса [1].

Достигнутое хорошее согласие говорит о том, что любое расхождение между расчетом и опытом обусловлено скорее выбором исходных данных по рассеянию, нежели ошибками в самих расчетах.

Сравнение с экспериментом. Описываемым методом был выполнен ряд расчетов для чистых элементов и некоторых бинарных систем при условии нормального падения элект­ ронов на поверхность мишеней. Полученные значения ко­ эффициентов обратного рассеяния электронов для ряда элементов сравниваются с последними экспериментальны­ ми результатами [12] в табл. 1. Абсолютная статистическая ошибка в этих коэффициентах составляет около + 0,01; однако, поскольку для задания каждой совокупности траекторий использовалась одна и та же система псевдослу­ чайных чисел, относительные ошибки должны быть зна­ чительно меньше указанной [11]. При 30 кэв значения коэффициента обратного рассеяния г) систематически на 8% выше экспериментальных значений, за исключением углерода. При 10 кэв положение иное. Все значения т) при

Т а б л и ц а 1

Сравнение вычисленных

методом Монте-Карло

(М.-К.)

и экспериментальных значений коэффициентов обратного

рассеяния

 

 

 

 

 

 

 

Е0 =• 10 кэз

 

Ее = 30

кэв

Z

Опыт

м.-к.

Отноше­

Опыт

м.-к.

Отноше­

 

 

 

 

ние

 

 

ние

6

0,072

0,0694

1,04

0,060

0,0598

1,00

13

0,171

0,1818

0,97

0,155

0,1664

0,93

22

0,268

0,2982

0,90

0,254

0,2722

0,93

29

0,339

0,3658

0,93

0,319

0,3390

0,94

47

0,420

0,4834

0,87

0,420

0,4470

0,94

79

0,501

0,6141

0,82

0,521

0,5678

0,92

этой энергии выше, чем при 30 кэв. Это согласуется с опы­

том

в области малых

атомных номеров,

но

для высоких

 

 

 

 

 

 

атомных

номеров

 

коэффици­

 

 

 

 

 

 

енты оказываются

преувели­

 

 

 

 

 

 

ченными.

Наиболее

вероят­

 

 

 

 

 

 

ным объяснением этого факта

 

 

 

 

 

 

может быть несправедливость

 

 

 

 

 

 

борцовского

 

приближения,

 

 

 

 

 

 

использованного

при

выводе

 

 

 

 

 

 

сечения

рассеяния

(1).

 

 

 

 

 

 

 

Описываемый

расчет

по­

 

 

 

 

 

 

зволяет также

найти

энерге­

 

 

 

 

 

 

тическое

распределение

об­

 

 

 

 

 

 

ратно

рассеянных

электро­

 

 

 

 

 

 

нов.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На

рис.

1

сравниваются

 

 

 

 

 

 

результаты для меди

и золо­

 

 

 

 

 

 

та, полученные для

рассея­

 

 

 

 

 

 

ния во всю заднюю полусфе­

Рис.

1. Энергетическое распределе­

ру,

с

экспериментальными

ние

обратно рассеянных

электро­

данными, определенными при

нов. Сплошные кривые —

экспери­

ментальные

результаты для

угла

угле

выхода

электронов

45°

рассеяния 135°. Черные и светлые

[12].

Согласие

в

низкоэнер­

точки — результаты

расчета

мето­

дом Монте-Карло для

меди и

золо­

гетической части распределе­

 

та

соответственно.

 

 

 

 

 

 

 

 

ния

хорошее, но

 

оно

суще­

ственно

нарушается в высокоэнергетической части. Это —

результат использования уравнения

(2),

которое

не учи-

тывает статистической природы потерь энергии электро­ нами в мишени и тем самым смещает пик наиболее вероят­ ной энергии. Можно сделать поправку на этот эффект, однако она вызовет существенное увеличение объема вы­ числений.

Электронные траектории можно также использовать для получения данных об эффективности и пространственном распределении актов возбуждения рентгеновского излуче­ ния в мишени. На рис. 2 показан пример результатов,

Y—\

с» з/ш ал

А1 :аш

0.5\

 

 

О

Т,5\

2,0

as

1,0

1,5

0,5

1,0

6)

fiz

 

Рис. 2. Вычисленное распределение по глубине актов ионизации для иониза­ ционного потенциала E J J = EJU. Светлые точки — результаты Кастена и

Декампа для меди с цинком в качестве меченого слоя. С ними следует сравни­ вать кривую для 1/U = 0,3.

вычисленных для меди (а) и алюминия (б) при начальной энергии электронного пучка 30 кэв; на графиках также при­ ведены экспериментальные результаты Кастена и Декам­ па [13]. При расчете использовалось сечение ионизации в форме, предложенной Уортингтоном и Томлином [14]:

Ц = Щ^Ъ1л-^,

6 = 0,35-1,37.

(3)

За исключением слишком быстрого спадания рассчитан­ ной кривой на больших глубинах, общая форма экспери­ ментальной и теоретических кривых одинакова. Различие в абсолютных величинах несущественно, будучи опреде­ лено нормировкой ф (pz). Эта нормировка чувствительна

кточной форме сечения ионизации, принятого в расчетах,

атакже зависит от конечной толщины пленки, используе­ мой для нормировки экспериментального распределения. Результаты Кастена и Энока [15] для алюминия с меченым слоем магния очень хорошо согласуются с вычисленными значениями при 30 кэв, но при 10 кэв согласие почти столь

же неудовлетворительное, как и отмеченное выше для меди.

Поправки на поглощение и атомный номер. Расчеты эффективности возбуждения и пространственного распреде­ ления актов возбуждения рентгеновского излучения в об­ разце могут быть использованы для определения поправок

Аюолишц номер Z

 

а)

в) " '

Рис. 3. Изменение средней глубины ионизации в зависимости от атомного номера при различных избытках напряжения V для Е0 = 10 (а) и 30 кэв (б).

на поглощение и на атомный номер. Здесь нет возможности воспроизвести полученные кривые / (%), однако, поскольку в первом приближении поглощение в мишени пропорцио­ нально средней глубине ионизации, из кривых рис. 3 можно получить некоторые указания на вид зависимости поправки на поглощение от атомного номера мишени. При больших избытках напряжения U доля поглощенного излу­ чения уменьшается с ростом Z при данном U и, наоборот, возрастает при малых U. Росту средней глубины возбуж­ дения для золота при 10 кэв можно не придавать значения ввиду неточности наших данных по рассеянию в этой области энергий.

Результаты расчетов методом Монте-Карло позволили найти возбуждение характеристического рентгеновского излучения элемента А, имеющего потенциал возбуждения

E J U .

Определим эффективность возбуждения РА следующим образом:

ПА = РАС А-

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