ными по ЭПР, а также двухфотонная спектроскоппя [2,8,10, И , 28, 30]. Главным инструментом теоретической интерпретации спектральных измерений является аппарат теории кристаллического поля [31—35], активно использующий методы теории групп [36].
2.1.Основные концепции современной теории кристаллического поля
В настоящей главе собран экспериментальный материал по энергии штарковских уровней активаторных ионов в некоторых лазерных диэлектрических кри сталлах. Но прежде чем рассматрпвать эти данные, коротко остановимся на основных положениях теории кристаллического поля. Читателям, желающим более подробно ознакомиться с техникой расчетов спектров активированных кристаллов, можно рекомендовать ряд обстоятельных монографий, например [28, 36, 44].
Основное приближение теории кристаллического поля связано с предполо жением о возможности введения эффективного гамильтониана активаторного иона с незаполненной электронной (d пли /) оболочкой в кристалле
Я = Я 0 + ЯИР, |
(2.1) |
где Я„ — гамильтониан свободного иона, Укр — энергия взаимодействия иона с решеткой, усредненная по квантово-механическим состояниям последней, так что FKp можно представить и как энергию электронов, локализованных на ионе-актпваторе, в эффективном внутрнкрпсталлическом электрическом поле. После выделения в Я 0 центрально-симметричного самосогласованного потен циала, определяющего одноэлектронные энергии, для определения уровней энергии каждой электронной конфигурации, задаваемой распределением элек тронов по состояниям с фиксированными величинами главного (п) и орбиталь ного (Z) квантовых чисел, требуется найти собственные значения оператора
V = Vee + Vs0 + УКр, |
(2.2) |
рассматриваемого как возмущение, где Vee — энергия |
электростатического |
взаимодействия электронов иона-активатора, У50 — оператор спип-орбиталь- ного взаимодействия.
В общем случае гамильтониан взаимодействия активаторного иона с кри сталлической решеткой записывается через одно- н двухэлектроииые опера
торы: |
|
|
= |
J v . t a . r ) , |
(2.») |
J
где гi — радиус-вектор S-го электрона активатора. В пространстве состояний электронной конфигурации (nl)N операторы Уг и У2 можно разложить по не
приводимым тензорным операторам Uqf) ранга |
к = 2т! (О |
тп/ |
1) [37]: |
|
VI = ^ B ^ )U ^ \ |
|
|
(2.4) |
|
У2 == S |
|
|
(2.5) |
|
|
PP'hg |
|
|
| р — р' | |
фигурные скобки в (2.5) означают кронекеровское произведение, |
||||
А |
I Р + р' |- Параметры кристаллического |
поля В ^ |
и ВдС) (рр1) являют |
|
ся функционалами от электронной плотности в кристалле с параметрической зависимостью от структуры решетки. В приближении парных взаимодействий
ионов (модель «суперпозиции») параметры кристаллического поля В ^ пред ставляются суммой по решетке:
= s Вт (Я„) С?>(0„, ф„), |
(2.6) |
где Rm, 0m и <pm — полярные координаты и Сдк) — функции направляющих
косинусов иона т (в системе координат с центром на активаторе), имеющие те же свойства преобразования, что и тензорные операторы TJq^ Различные па
раметры 2?дК) связываются уравнениями, вытекающими из инвариантности оператора (2.3) при преобразованиях точечной группы симметрии активатора.
В диэлектрических кристаллах волновые функции локализованных d- и /-электронов активаторных попов заметно перекрываются только с волновыми функциями внешних электронных оболочек ближайших соседей-лигандов, поэ тому координаты лигандов полностью определяют величины параметров
Bq') (рр1) в операторе (2.5), обусловленном прямым обменным взаимодействием электронов активатора и лигандов [37].
Многочисленные микроскопические расчеты энергетических спектров акти ваторных центров, содержащих исследуемый ион и его ближайшее окружение, позволили выяснить многие вопросы, касающиеся природы и относительной роли различных механизмов взаимодействия активатора и решетки (см., на пример, [38, 39]). В частности, вклады двухэлектронных операторов в наблю даемые штарковскпе расщепления составляют, по-видимому, не более несколь ких процентов [37].
При интерпретации спектральных данных параметры кристаллического поля определяют путем подгонки спектров, рассчитанных с гамильтонианом (2.4), к результатам экспериментальных измерений. Так как апроксимацпя Укр одноэлектронными операторами заведомо не допускает точного описанпя структуры спектра (отметим также, что в наблюдаемые штарковскпе расщепле ния заметный вклад (до 10%) впосят смещения уровней энергии, индуциру емые электроп-фоноиным взаимодействием), подобная процедура подгонки часто приводит к неоднозначным результатам, лишенным физического содер жания. Анализ экспериментальных данных должен основываться на предва рительном расчете, исходя из микроскопической модели актпваторного центра, затравочных параметров кристаллического поля, уточняемых в последующем с сохранением качественных соотношений между ними сравнением рассчитан ной и экспериментальной штарковской структуры, а также спектров ЭПР. Та кой подход не только позволяет получить физически обоснованные параметры кристаллического поля, но и в ряде случаев помогает расшифровке спектра и правильному выбору модели актпваторного центра [40, 41].
Первые расчеты параметров кристаллического поля базировались на точеч ной модели ионной решетки: активаторный ион рассматривался в электрическом
поле точечных ионов с эффективными зарядами zme. Параметры |
в этом слу |
чае пропорциональны |
|
^ a„)e‘Zmirky R ^ \ |
(2.7) |
где <УС> — среднее значение к-й степени радиуса оптического электрона; щ. — константы линейного экранирования, введение которых позволяет учесть для Ln- и Ас-ионов влияние внешних 5s2, 5р6 или 6s2, 6рб оболочек, деформирован ных в кристаллическом поле [42]. Вычисленные в рамках точечной модели па раметры кристаллического поля при использовании величин <гк)>, полученных с одноэлектронпыми хартри-фоковскими радиальными функциями свободных ионов [28], обычно совпадают с экспериментальными данными по знаку, но существенно различаются по величине. Учет электрической дипольной (для ионов матрицы в узлах без ииверспд) п квадрупольной поляризации ионов в кристалле не улучшает согласия с экспериментом; просуммированные по
решетке вклады в Вд^ от полей точечпых зарядов, дипольных и квадрупольных моментов (найденных.с поляризуемостью свободных ионов) близки по величине и могут иметь разные злаки [43]. Пренебрежение эффектами блпзкодействня — пространственным распределением электронной плотности на лигандах —
не позволяет в рамках точечной модели корректно описать ни соотношений ме жду различными параметрами кристаллического поля, ни характера их изме нения при всесторонней и одноосной деформации решетки. Поскольку квантовомеханические расчеты энергетического спектра активаторных ионов весьма трудоемки, а получаемые численные оценки имеют небольшую точность, интер претация штарковской структуры оптических спектров активаторов в последние годы проводится в рамках различных моделей с полуфеиоменалогической па раметризацией энергии оптических электронов. В частности, широкое распро странение получил упрощенный вариант модели «суперпозиции» (см. обзор [39]). Он основывается на предположении о полной взаимной компенсации вкла дов в параметры кристаллического поля от электрических мультиполей, лока лизованных на ионах матрицы вне активаторного центра [45]. Соответствен© в суммах по решетке (2.6) учитываются только вклады лигандов, и при заданной геометрии ближайшего окружения величины штарковских расщеплений пол
ностью определяются одноименными параметрами |
(Вт). В предположении |
||||
степенной зависимости параметров |
от |
расстояния R m между лигандом |
|||
и активатором В ()!) (R m) = (i?0/i?m)t,i‘5 (J,') (R 0) |
(взаимодействие |
фиксированной |
|||
пары ионов |
активатор-лиганд характеризуется соответственно |
четырьмя (f2, |
|||
i4, В(2\ В(4)) |
в случае d-электронов и шестью |
(t2, t4, t6, В^2\ BW, |
ВW) в случае |
||
/-электронов параметрами). Используемые в рамках модели «суперпозиции»
параметры кристаллического поля В ^ и Bq^ хорошо согласуются с результа тами микроскопических расчетов; в ряде случаев рассмотрение штарковской структуры спектров активированных кристаллов на основе этой модели поз волило выполнить анализ локальной деформации решетки вблизи актпваторных ионов [46, 47], однако искусственное завышение роли близкодействия часто не позволяет получить физически обоснованных степенных показателей f2 для квадрупольвых параметров кристаллического поля [39].
Дальнейшее развитие задачи параметризации эффективного гамильтониана взаимодействия активаторных ионов с решеткой получила в работах [48, 49], где была сформулирована модель обменных зарядов, использованная в даль нейшем в расчетах статических и динамических спектральных характеристик как высокосимметричных (кубических), так и низкосимметричных активатор ных центров Ln-ионов [40,41, 50]. Модель обменных зарядов в теории кристал лического поля основывается на тех же представлениях о распределении элек тронной плотностп в кристалле, что и аналогичная модель в теории динамики кристаллических решеток — на связях между ионами с перекрывающимися электронными оболочками помещаются фиктивные положительные заряды, пропорциональные квадратам интегралов перекрытия волновых функций [51]. Электростатические взаимодействия валентных электронов с точечными зарядами и дипольными моментами на ионах матрицы рассматриваются в явпом виде, а параметризуется только та часть потенциала кристаллического поля, которая создается обменными зарядами (при этом нужно учитывать отличие поляри зуемости ионов в кристалле от поляризуемости свободного иона, а также по ляризацию, обусловленную деформацией электронных оболочек соседних ионов
сперекрывающимися волновыми функциями).
Влюбой теории кристаллического поля, основанной на рассмотрении парных взаимодействий, зависимость вкладов отдельных ионов в параметры кристал лического поля от угловых координат 0т и (pw остается той же, что и в точечной модели (см. выражение (2.6)); поэтому учет перекрывания и неортогональности волновых функций валентных электронов активатора и лигандов сводится к
к перенормировке средних Чтобы перейти от точечной модели к модели
обменных зарядов, достаточно в параметрах B{1i) (Rm), |
записанных в точечной |
|||
модели (см. выражение (2.7)), заменить |
на |
|
||
<^> |
2 (2*+,) |
О — «,) ■Й'’(Ят ), |
|
(2.8) |
(21 + 1, |
|
|
||
где Zm — валентность лг-го лиганда, I — орбитальный момент электронов ак
тиватора, G,. — безразмерные параметры модели, 5^ (Rm) — линейные ком бинации квадратов интегралов перекрытия, которые в случае лигандов с внеш ними заполненными s2pe оболочками равны
sil) = s* -|- s l + P p s l
p f = - p<3) = » /„ p f = l / a
для /-электронов,
P f ] = 1, P?> = - y a, Рц2) =0
для fZ-элсктронов. Здесь Ss, S0 и Sn — двухцентровые интегралы перекрытия волновых функций /- и d-электронов с s-, р,- и /^-функциями электронов лиганда соответственно.
В рамках простейшего варианта модели с использованием только одного подгоночного параметра G —(?|г, характеризующего величины обменных зарядов н постоянного для конкретной пары активатор-лнганд, возможно удовлетвори тельное описание как штарковской структуры спектра активатора в одной матрице, так и ее изменений в пределах гомологического ряда кристаллов. Следует отметить, что реализация расчета параметров кристаллического поля обязательно требует предварительного рассмотрения локальной деформации решетки вблизи активатора.
Различные внутрнкристаллпческие поля в зависимости от соотношения
между Vee, |
Fso и У1(), принято подразделять, следуя |
[31], на три типа: слабое |
|
поле (У^, |
У?0^ У |ф ), средпее {Vee >> Укр !>> Vs0) |
п сильное (У1ф |
Vee |
У80). Случай слабого поля реализуется для Ln- и Ас-попов, у которых внут ренние незаполненные оболочки 4/ н 5/ сжаты п заэкранированы внешними за полненными оболочками 5sap6 п 6s2p° соответственно. Если связь между элек
тронами имеет в основном нормальный характер (Уее У*0, Ln-поны), электро статическое взаимодействие электронов приводит к образованию термов, ха рактеризуемых определенными значениями полного орбитального п спинового моментов количества движения (Ln S). Сшш-орбитальное взаимодействие У^0 = = XLS часто снимет (2L + 1)(25 1)-кратное вырождение термов и формирует мультиплетпую структуру спектров свободного иона (здесь X — константа спин-орбитального взаимодействия). При рассмотрешш штарковского расщепле ния в слабом поле в качестве функций нулевого приближения выбираются соб ственные функции | L, S, J, /- ) полного момента J = L + S. Этот случай показан на рис. 2.1, где приведена схема расщепления энергетических уровней (на примере основного состояния 41»/,) нона Nds+ в кристалле YA103.
В случае среднего поля в качестве волновых функций нулевого приближения используются собственные функции | L, L:, S, S:y операторов L п S. Штарковское расщепление в среднем поле меньше энергетических зазоров между разными термами одппаковой мультиплетностн (с одинаковыми значениями полпого спина), но больше расщеплеппя топкой структуры в свободном ноне. Сшш-орбптальное взаимодействие наряду с низкоснмметрпчными составляю щими кристаллического поля обычно рассматриваются как возмущение после учета кубической компоненты кристаллического поля. Следует отметить, что в чистом виде случай среднего поля не встречается. Обычно на ионы группы же леза действует промежуточное кристаллическое поле (Fcf ^ FHp), ц следует рассматривать в качестве возмущения сумму Vee -f- FKp, либо учитывать вза имодействие термов при последовательной дпагонализации матриц F«., FKp.
Сильное кристаллическое поле разрывает связь между орбитальными и спиновы м и моментами электропов па незаполненной оболочке; электронная конфигурация задается в этом случае числами заполнения одпоэлектрониых
|
|
Nd: |
ф д а Л we |
|
/— |
Л |
foau/wdtucmfae |
|
|
— ч |
|
I |
|
/2^си 1 |
|
|
- Г |
|
|
|
|
|
Сявн-ербвтал&яее
8еавмаеевсяг0ие
/ ^ Г
/- 1
Л \к
/0 си
rfj/rif/ffffo ircm lJ zi* 0 faympaxpi/cmaswffYfwsrtf
лелем
288СМ1 |
/ // |
Рис. 2.1. Схема расщеп ления уровнен нона Nd3+ в кристалле YA10g
I |
— термы, |
характеризуе |
|||
|
мые квантовыми |
чис |
|||
|
лами L и S ] |
|
|
|
|
I I |
— мультиплсты, |
харак |
|||
|
теризуемые |
квантовы" |
|||
|
ми числами |
L, S и J \ |
|||
I I I |
— пггарковские |
уровни, |
|||
|
соответствующие |
не |
|||
|
приводимым представ |
||||
|
лениям |
группы |
ло“ |
||
|
калькой |
симметрии |
|||
1
штарковских уровней, энергия которых находится при диагонализации опера тора ГКр в базисе одноэлектронных функций ( llzy. В сильном кристалличе ском поле, действующем на ионы элементов групп палладия и платины, особен но существенны эффекты, которые обусловлены образованием ковалентной связи между актпваторными ионами и лигандами.
Количество штарковских уровней в кристаллических полях разной симме трии в зависимости от значений квантовых чисел I (сильное поле), L (среднее) и / (слабое) определяется методами теории групп. В табл. 2.1 приведены числа штарковских компонент для активаториых центров с Ьп3+-ионами с различ ной симметрией прп значениях / от 0 до 8 (целым значениям J характеризуются Ьп8+-ионы с четным числом электронов на незаполненной 4/-оболочке) и от 1/2 до 17/2. Эта таблица позволяет сравнивать количество штарковских уровней при различной симметрии (структуре) активаторного центра. Например, для иона Nd3+, самого распространенного лазерного активатора в полях любой симметрии, кроме кубической, количество штарковских уровней у данного мультиплета всегда будет одинаковым. Это обстоятельство в значительной степени облегчает расшифровку оптических спектров для построения схемы энергети ческих уровней. Для Ьп3+-ионов с четным числом электронов, например для Рга+ или Но3+, в полях различной симметрии количество штарковских уровней будет различным. В этом случае для интерпретации оптических спектров знание точной симметрии активаторного центра является необходимым.
2*2. Энергия уровней генерирующих активаторов
влазерных кристаллах (экспериментальные данные)
Кнастоящему времени СИ возбуждено на волнах различных межштарковских переходов генерационных каналов (см. табл. 1.10) активаториых ионов у бо лее 400 лазерных кристаллов (матрица-основа -J- один тип генерирующего