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

книги из ГПНТБ / Немошкаленко В.В. Теоретические основы рентгеновской эмиссионной спектроскопии

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

330 Глава 4. Одноэлектронное и многоэлектронное приближения

наблюдается в значениях nd для хрома и ванадия, что, возможно, объясняется тем, что с уменьшением атомного номера элемента d-волновые функции становятся менее локализованными и возрас­ тает число электронов, находящихся между атомными сферами. Действительно, согласно значениям, полученным в работе [314], эта величина равна для никеля 0,49, а для хрома — 1,23. Кроме

того, в работе [289]

для большинства элементов принималось произ­

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

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

веса s-электронов

в валентной полосе (0,5):

 

 

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

интенсивности в рентгеновских эмиссионных

полосах меди вычислялось не

только методом

функции Грина.

Рис. 145. Экспериментальная

(---- ) н рассчитанная

(-------- ) /(-полоса

меди [327].

Первый расчет распределения интенсивности в эмиссионной /СРг,5-по­ лосе меди был выполнен методом ячеек [327] (рис. 145). Однако согласие с экспериментом при этом было несколько хуже, хотя рассчитанная эмиссионная /Г-полоса, как и экспериментальная, имела два максимума. Это объясняется непригодностью метода ячеек для описания зонной структуры переходных металлов.

Эмиссионные Lin- и Мщ-полосы меди рассчитаны также мето­ дом ППВ [311 ]. Найдены парциальные плотности состояний Ns (Е), Np (Е ) и Nd (Е ) и полная плотность состояний N (Е ). Расчет вы­ полнен для 7000 точек в 1/48 части зоны Бриллюэна. Детали расчета

и

использовавшаяся

интерполяционная схема

подробно описаны

в

работе [328]. Как

отмечают авторы работы

[311], парциальная

плотность состояний //-симметрии Np (Е ) мала по всей изученной области энергий. Это, собственно говоря, и послужило основанием для расчета распределения интенсивности в Бщ- и Мщ-полосах, которые в основном определяются парциальной плотностью Nd(E).

На рис. 146, 147 приведены результаты теоретического расчета эмиссионных Бщ- и Л4ш-полос меди. Полосы имеют хорошо выра­

женную

тонкую структуру. Применение процедуры, описанной

в

работе

[294], для

учета Оже-эффекта в валентной полосе при

Г0

= 2

эв в основном

не повлияло на результат. Полоса не расши­

рилась,

несколько лишь сгладилась ее тонкая структура. Вместе

с тем, как отмечалось выше, остовные Lm- и Мщ-уровни меди очень широкие и связанный с этим эффект «замазывания» тонкой

Одноэлектронное приближение

331

структуры N (Д)-кривой должен быть большим. Действительно, экспериментально наблюдающиеся формы эмиссионных Дц- и Мш-полос меди гладкие и не имеют тонкой структуры. Экспери­ ментальная Lin-полоса [329] может быть сопоставлена с теоретиче­ ски рассчитанной, так как из экспериментальных исследований известно, что энергетическое разделение Дьш-дублета составляет около 20 эв. Для Миди-дублета эта величина равна всего 2,5 эв. Поэтому на рис. 147 кроме суммарной Мц,ш-полосы [330] приве-

Рис.

146.

Рассчитанная [311]

Рис. 147. Сопоставление теорети­

(—)

и

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

ческой Л41П-полосы меди [311] (---- )

[329]

(-------- ) Д п -полоса

с экспериментальной Мц

[330]

 

 

меди.

(------) и «экспериментальной» Л4щ-

 

 

 

полосой, полученной при разложе­

 

 

 

нии Л4П (ц-полосы (-------- ).

дена выделенная из нее УИцрполоса в предположении указанного

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

разделения Мц.пгДублета и отношения интен-

Ми

= 2, т. е. отношения статистических весов рассмат-

сивностеи —г-:—

риваемых уровней. Согласие рассчитанных и экспериментальных кривых, особенно в области главного максимума, неудовлетво­ рительное. На экспериментальных кривых наблюдается лишь низкоэнергетический наплыв, имеющий аналог на рассчитанной кривой интенсивности этих полос. Основная причина расхождения экспериментальной и теоретической форм Lu- и Мш-полос, по-ви­ димому, заключается в значительной ширине остовного уровня. Этот эффект может существенно сгладить тонкую структуру в об­ ласти главного максимума, где Оже-расширение играет значитель­ но меньшую роль. Нельзя согласиться с предположением авторов работы [311] о том, что основная причина отсутствия в спектре Д - и М 1-полос состоит в малой плотности состояний р-симметрии в ва­ лентной полосе меди, как это следует из проведенного ими расчета парциальной плотности состояний Np (Е). Выше отмечалось, что даже при такой небольшой плотности состояний хорошо фикси­ руется /(-полоса меди.

Рассчитано распределение интенсивности в рентгеновских эмисси­ онных полосах германия. Германий — полупроводник, электронная

Рис. 148. Эмиссионная /^-полоса германия:

р а с с ч и т а н н а я / — к - р — О П д - м е т о д о м м 2 — с у ч е т о м ш и р и н ы остов*- н о г о у р о в н я и О ж е - р а с ш н р е н и я у р о в н е й в в а л е н т н о й п о л о с е , 3 —» э к с п е р и м е н т а л ь н а я [ 2 7 4 ] .

р а с с ч и т а н н а я / — k - р

— О П В м е т о д о м и 2 — с у ч е т о м ш и р и н ы

о с т о в н о г о у р о в н я н

О ж е р а с ш и р е н и я у р о в н е й в в а л е н т н о й

 

п о л о с е .

р а с с ч и т а н н а я / — к • р — О П В - м е т о д о м н 2 — с у ч е т о м ш и р и н ы о с т о в н о г о у р о в н я и О ж е - р а с ш и р е н н я у р о в н е й в в а л е н т н о й п о л о с е .

 

Одноэлектронное приближение

 

333

структура

которого сходна с

электронной

структурой

крем­

ния. В рассматривавшейся работе

Клима

[299]

наряду с

К- и

L-полосами кремния рентгеновские эмиссионные полосы германия

рассчитаны

к • р — ОПВ-методом

(рис.

148— 150). Рассчитанная

форма /(-полосы германия хорошо согласуется с результатами эксперимента [274].

Элементы второго большого периода. Структура валентной по­ лосы и распределение интенсивности в рентгеновских эмиссионных

L-, М- и /V-полосах ниобия, мо­

 

 

 

либдена, родия и палладия с

 

 

 

ОЦК

и

ГЦК

решетками

рас­

 

 

 

считывались

методом

функции

 

 

 

Грина соответственно

в 55

и 85

 

 

 

неэквивалентных точках непри­

 

 

 

водимой части зоны Бриллюэ-

 

 

 

на [331]. При помощи метода

 

 

 

интерполирования

число

то­

 

 

 

чек

в

полной

зоне

Бриллю-

 

 

 

эна

было

увеличено

соответ­

 

 

 

ственно

до

91 520 и

178 870. К

 

 

 

сожалению,

расчет

проводился

 

 

 

в нерелятивистском

приближе­

Рис. 151. Энергетическая зависимость

нии, а при расчете

зонной

матричных

элементов

радиальных

структуры

тяжелых

 

металлов

факторов

вероятности перехода Р2р d:

(см. гл. 3) учет релятивистских

/ — д л я м о л и б д е н а , 2 — н и оби я ,

3 — р о д н я и

эффектов приводит к появлению

 

4 — п а л л а д и я .

 

дополнительных

расщеплений в

 

 

 

центре валентной зоны. Однако полученные даже в таком прибли­ жении результаты [331] важны потому, что, как и в случае эле­ ментов первого большого периода [314], получена не только об­ щая плотность состояний, но и парциальные плотности состояний Ns (Е ), Np (Е ), Nd (Е ) и соответствующие числа заполнения. Значения ns, пр и nd хорошо согласуются с аналогичными значе­ ниями для элементов первого большого периода, полученными на основе измерений относительных интегральных интенсивностей

[289].

При расчете распределения интенсивности в рентгеновских эмиссионных полосах учитывались следующие матричные элементы

радиального фактора вероятности

перехода:

для Бщ-полосы —

P ip,d (Е), для Му-полосыР3а,р (Е)

и для /7ш-полосы — P4pd (Е)

(соображения те же, что и при расчете К-, Е-

и М-полос элементов

первого большого периода). В образовании Му-полосы принимают участие также / 3d переходы, однако они во внимание не прини­ мались, так как в зонном расчете /-электроны не рассматривались.

Результаты расчета энергетической зависимости матричных эле­ ментов P2p,d (рис. 151), P3d'p (рис. 152) и P4p,d (рис. 153) подтверждают

334 Глава 4. Одноэлектронное и многоэлектронное приближения

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

для

железа [310]. Радиальные р-функ-

ции валентных электронов

слабо зависят от энергии, вследствие

чего слабо зависит

от энергии и матричный элемент Р^,р (Е ).

 

 

2

Рис. 152. Энергетическая зави­

Рис. 153. Энергетическая зависи­

симость матричных

элементов

мость матричных элементов ради­

радиальных

факторов

вероят­

альных факторов вероятности пе­

ности

перехода P3d р:

рехода P4pd:

/ — д л я н и о б и я , 2 — м о л и б д е н а .

/ — д л я н и о б и я , 2 — м о л и б д е н а .

С другой стороны, заметное повышение степени локализации d-функций с ростом энергии приводит к тому, что матричный эле­ мент P 2p,d возрастает в два раза при переходе от дна полосы к

Рис. 154. Сопоставление экспериментальных [332, 333] (------) и рассчитанных [331] без учета (—О —) и с учетом (-------- ) Оже-расширения /Иу - и Nn ш -полос ниобия и мо­

либдена.

уровню Ферми у ниобия и молибдена и примерно в полтора раза у родия и палладия.

Рассчитанные спектры размывались на ширину остовного уровня и Оже-расширение в валентной полосе [294]. Результаты

 

Эмиссионные

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

спектры

неупорядоченных сплавов

335

расчета

Му- и

Л^п.ш-полос

ниобия

и молибдена приведены

на

рис.

154,

а Z-ш-полос ниобия, молибдена, родия и палладия— на

рис.

155.

При учете Оже-расширения спектр сильно сглаживается

и теоретически рассчитанное распределение интенсивности хорошо согласуется с экспериментально определенным. Расхождение со­ стоит в том, что в области максимума Л у Му-полос ниобия и молиб­ дена на теоретически рассчитанных кривых наблюдается расщепле-

- 5

0

5

- 5

0

5 Е,рид

Рис. 155. Сопоставление экспериментальных ( ----

) и рассчитанных [331 ]

без учета (—О —) и с учетом (--------

 

) Оже-расширения Дц-полос

ниобия,

молибдена, родия и палладия.

ние, которого нет на экспериментальной кривой. В остальных слу­ чаях согласие удовлетворительное как по форме полос, так и по положению максимумов. Рассчитанные полосы во всех случаях несколько уже, однако, как отмечалось в гл. 3, ширина валентной зоны существенно зависит от выбора потенциала. В данной работе он определялся несамосогласованным методом. Вместе с тем, как показали расчеты, одноэлектронное приближение может весьма успешно применяться для количественной интерпретации рентге­ новских эмиссионных полос. Дальнейшее уточнение расчетов как в рамках одноэлектронного подхода, так и с учетом многоэлектрон­ ных эффектов должно улучшить согласие с экспериментом.

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

Пусть система электронов сплава описывается гамильтониа­ ном Н0. Введем в рассмотрение оператор возмущения учиты­ вающий взаимодействие между системой и внешним полем — линейно-поляризованной плоской волной с волновым вектором q и

336 Глава 4. Одноэлектронпое и многоэлектронное приближения

комплексным единичным вектором поляризации п. Опреатор Н' можно записать в виде

Н' = Н ехр (— Ш ) +

Н+ ехр Ш ,

 

где

 

 

Н = 2 (X \ W \ l')a ta v ,

(358)

XX

 

 

W = ехр i<\ ■rn •р,

 

|Я ) — произвольный полный набор

ортонормированных

одно­

электронных функций, о?}- и ах — операторы рождения и уничто­ жения электрона в состоянии X. В качестве полного набора ортонормпрованных одноэлектронных функций удобно использовать волновые функции остовных |R , п) и внешних валентных |R, I) электронов (R — узел решетки, п и I — номера соответственно остовных и валентных состояний). Будем считать, что электроны, находящиеся на остовных уровнях, не взаимодействуют друг с другом и с валентными электронами. Таким образом, гамильтониан Н0 в этом приближении можно записать в виде

Н0 = Н 0 -{- V eR„aitiORji,

Rn

где eR„ — энергия остовного уровня, Н0 содержит только операторы рождения и уничтожения валентных электронов. В таком прибли­ жении гамильтониан Н0 не учитывает влияние дырки остовного уровня на валентные электроны. Как показано в работе [334], этот эффект существен лишь в области энергий вблизи энергии Ферми. В случае эмиссии система электронов находится в возбуж­ денном состоянии |i) с дыркой на остовном уровне. Вследствие взаимодействия с электромагнитным полем система переходит из начального в конечное стационарное состояние. В случае эмиссии возможны также спонтанные переходы. При учете спонтанных переходов изменяется лишь коэффициент в формуле для интенсив­ ности излучения, а это не отражается на форме получаемых спект­ ров. Интенсивность рентгеновской эмиссии в линейном приближе­ нии по внешнему полю описывается соотношением

 

/ N

- 2 6 (£ .« ~

E i +

<*){i\H |m) (m |Я + 1i),

(359)

 

 

m

 

 

 

где

Em и |m) — соответственно

собственные значения

энергии

и

собственные

функции

гальмильтониана # 0. Учитывая, что

рентгеновский переход осуществляется между остовными и валент­

ными состояниями,

оператор (358) можно представить в виде

^ = 2 2

O ^ R t/t .R jn sflR ^ O R ^ + ®?R2n!,R|/tO^,n,flR1/I)-

Ri l\ Rj/ij

 

Эмиссионные

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

спектры неупорядоченных сплавов 337

Считая, что в начальном состоянии дырка

находится

в зоне J R , п),

а в конечном состоянии — в валентной зоне, валентные электроны

описываем соответственно

волновыми

функциями

|0') и |т')

гамильтониана Н0. Все остовные состояния заполнены. Для

матричного элемента </|#|т> с точностью

до знака

получаем

выражение

 

 

 

 

 

 

 

 

(i IН |т) =

2

(0' |at,!,

|m').

 

(360)

 

 

 

R,/,

 

 

 

 

 

Как показано в работе [335], подстановка матричного элемента

(360) в выражение

(359) приводит к соотношению

 

 

 

 

/ (со) — 1т 0 ((.I — со — eR„) х

 

 

 

X

^

У ]

U ^ R

n . R ,

 

( с о

- f -

8 R n ) ^(361)R „ / : , R n ,

К.г.

 

 

 

 

 

 

 

где GRtiuR2u (со)— опережающая функция Грина,

 

 

 

 

0 Н

-

‘ ■ если “ > 0 '

 

 

 

 

 

 

 

(0, если со < 0 .

 

 

 

Из определения функции 0 (со) вытекает неравенство

со < (.1 — eR„,

свидетельствующее о том, что испускание фотонов происходит в результате перехода электрона из заполненных состояний валент­ ной зоны. Интенсивность рентгеновской эмиссии I (со) определя­ лась в предположении, что в кристалле есть одна вакансия п на узле R, которая может быть заполнена. Если считать, что вакан­ сия может быть на любом узле, то выражение (361) следует просум­ мировать по всем узлам с вакансией.

Рассмотрим анализ [335] соотношения (361). Считая, что вакан­ сия расположена на узле 0, с точностью до численного множителя можно записать

I (со— Eon) —

2 1^0n,R1/,(JRl/„R2/!ll7R2/2,0/i.

 

Ri/t R2/2

В случае многокомпонентного сплава необходимо произвести усред­ нение по всем атомам, окружающим исследуемый атом:

Лх (® — Е0/>) — ' т 2 2 ^^.Ri/i^Ri/t.Rs/г (ш) ^ R./o,0n,

(362)

r , ; , r j 2

 

где а — сорт атома с вакансией на узле 0. Таким образом, интен­ сивность рентгеновского излучения можно представить в виде бес­ конечного ряда одночастичных функций Грина. Коэффициентами перед функциями Грина являются произведения матричных эле­ ментов перехода W. Матричные элементы функций Грина и опера­ тора W определены на одноэлектронных узельных функциях

338

Глава 4. Одноэлектропное и многоэлектронное приближения

|R, />

и |R, п),

составляющих

полный

ортонормированный

набор.

Функции

|R , /> и (R, п)

могут

быть вычислены по

атомным функциям методом, предложенным в работе [336]. Для остовных состояний нет необходимости строить такие функции, так как атомные остовные функции перекрываются незначительно,

а в этом

случае базисные

функции |R, п)

практически

будут

совпадать

с

|R, /г>а.

Базисная функция

валентных

электронов

]R, 1)а является

линейной комбинацией атомных функций

] R, />а,

IR,

0

=

I R, 0 a

- 4 - 2

l R + Р - O aa( R + P . n R ,

Da,

(363)

 

 

 

 

р

 

 

 

 

 

 

где р — радиус-вектор

ближайшего

соседнего

узла.

 

 

Рассмотрим

двухкомпонентный

сплав,

состоящий

из

атомов

сорта А и В, беспорядочно расположенных в узлах простой куби­ ческой решетки. Пусть Сд и Св — концентрация соответственно компонентов Л и В в сплаве. Для простоты будем полагать, что каждый атом имеет только одно валентное s-состояние. Собствен­ ные функции этого состояния одинаковы для атомов обоих сортов и настолько локализованы, что перекрываются только с волновыми функциями соседних узлов. Пусть гамильтониан Н0 в одноэлект­ ронном приближении имеет вид

 

 

 

Но = 2

LRj^rfaijcin.i,

 

 

 

 

 

 

R,R.

 

 

 

 

где L =

Т + V (Т — оператор

кинетической энергии

электрона,

V — потенциал решетки). С точностью до первого порядка

по пе­

рекрытию

атомных функций

[R, 1)а, используя разложение ба­

зисной функции

(363),

матричный

элемент оператора

L

можно

записать

в

виде

 

= <j( R i . l\T +

V | R 2, l)a

 

 

 

 

 

^ r,/,r2/

 

 

 

 

------ g-

„ ( R i. 1 1R 2. D a (e (R l .

/ I T + V |R lf D a +

 

 

 

 

+ a ( R 2> ^ I T + V |R 2, / ) о) (1 — S r jRj).

 

 

Выделим в этом выражении диагональную и недиагональную части:

£ r,/,R2/= a(Rl> Ц Т + У I ^1> Oa^RjRj = Е Rj^RjR,.

Если Ra и R2 — соседние узлы, недиагональные члены отличаются от нуля. Тогда

^RH,R2/ =

a ( R . ЦТ + ^ I R +

P> D a

------ 5- ( E r

+ ^R+p) a(R> l I R

+ P> Da-

Предположим, что матричные элементы не зависят от сорта атомов,

находящихся на узлах R и R + р. Следовательно, гамильтониан

Эмиссионные рентгеновские спектры неупорядоченных сплавов 339

Н0, для которого должна быть вычислена функция Грина, имеет вид

Я 0= 2

+ Г V 6 Rl,R1+ p ) a R 1;aR „/.

Ri Rs

р

Найдем матричный элемент Ц?0и,дг с точностью до второго порядка по перекрытию атомных функций, используя разложение (363):

^On.R/ = а(0> I ^ I R> О а

- 4

- S e(0,ft|^|R + p, Ofle(R + p,Z|R, 1)а.

(364)

z

р

 

Так как атомные функции остовных электронов не перекрываются, то, оставляя в этом выражении только члены первого порядка по перекрытию их, получаем

^ ( K R / = a(0> n \ W\ 0 , / ) q6 ro

-----Y а(0/г |И7 |0/)a a(0/1 R/)a (1 — 6r0).

Введем обозначения

 

a(0,n\W\0,l) = W,

---- i- (0, 11p, l)a = 5.

В этих обозначениях формула (364) имеет вид

^on,R/ =

2 Srp-

 

р

Подставляя это выражение в формулу для интенсивности рентге­ новского излучения (362), находим

 

 

(СО -

eon) ~

I w f

Im t e ( со) + S 2[G ^(co) +

 

 

 

 

 

+ Gop> (со)] + S2 2

GPlPo (со)l .

(365)

 

 

 

 

 

 

Pi Pj

j

 

Следует помнить,

что в этом выражении узел 0 занят атомом сорта

а,

а р,

рх,

р2 — радиусы-векторы

ближайших атомов-соседей

(усреднение

выполнено

по

конфигурациям атомов, занимающих

все

узлы,

кроме

узла

0).

Численное

значение функции

Грина в

(362) зависит лишь от того, какого сорта атом, А или В, занимает

узел 0, и не зависит от положения этого узла в кристалле.

В по­

следнем члене оставляем только основной вклад 2 Gpp, т. е.

члены

р

 

С Pi = Р'2 = Р-

Таким образом, формула (365) связывает интенсивность излу­ чения со степенью перекрытия атомных функций. Мнимая часть

функции Gоо (со), как известно (см. стр. 234), называется локальной

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