
книги из ГПНТБ / Немошкаленко В.В. Теоретические основы рентгеновской эмиссионной спектроскопии
.pdf330 Глава 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 (Е). Выше отмечалось, что даже при такой небольшой плотности состояний хорошо фикси руется /(-полоса меди.
Рассчитано распределение интенсивности в рентгеновских эмисси онных полосах германия. Германий — полупроводник, электронная
|
Одноэлектронное приближение |
|
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), находим
|
|
1а (СО - |
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), называется локальной