Материал: Многоуровневые функциональные схемы кристаллических лазеров

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

ными по ЭПР, а также двухфотонная спектроскоппя [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- один тип генерирующего