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

Москвин А.С. Панов Ю.Д. Атомы в кристаллах. - Изд-во УрГУ, Екатеринбург, 2009

.pdf
Скачиваний:
54
Добавлен:
13.03.2016
Размер:
610.15 Кб
Скачать

Энергетическая щель между ветвями адиабатического потенциала (янтеллеровское расщепление) в точках минимумов равна

V 2

ε+0) − ε0) = 2 E = 4EJT . (5.39)

KE

Электронные функции, соответствующие двум ветвям ε± , находятся

ˆ

как собственные векторы матрицы H1 (5.34).

ψ

= sin ϕ ϕ1 + sin ϕ ϕ2

)

при VE > 0

(5.40)

ψ+

= cos ϕ2

ϕ1 + sin ϕ2

ϕ2

 

 

 

 

2

2

 

 

 

 

(при VE < 0 вектора меняются местами).

ˆ

 

 

 

 

 

 

 

, желоб "мексиканской

С учетом квадратичной вибронной связи H2

шляпы" гофрируется, на нем появляются три эквивалентных минимума при ϕ = 0, 23π , 43π , (или ϕ = π3 , π, 53π ), разделенных потенциальными барьерами с высотой, пропорциональной константе квадратичной вибронной связи. Эти минимумы определяют три эквивалентные тетрагональные искажения октаэдра (вдоль осей x, y или z).

Решение системы вибронных уравнений (5.20) даже в случае простой линейной E − e-задачи сложно. Задача может быть упрощена в предельных случаях слабой или сильной вибронной связи, определяемых величиной отношения энергии EJT ян-теллеровской стаби-

лизации к энергии кванта колебаний в желобе (вдоль координаты ρ) p

~ωE = ~ KE/ME

EJT 1 − слабая вибронная связь,

~ωE

(5.41)

EJT 1 − сильная вибронная связь.

~ωE

Впределе слабой вибронной связи для решения системы (5.20) обычно пользуются теорией возмущений.

Впределе сильной вибронной связи для линейной E − e-задачи становится применимым рассмотрение нижних вибронных состояний

врамках адиабатического приближения.

Вэтом приближении возможно разделение радиальных (ρ) и угловых (ϕ) переменных в вибронной функции для каждой ветви адиабатического потенциала, так что

ψ±(r, ρ, ϕ) = ϕ±(r, ϕ) Φ(ϕ)

χ(±)(ρ)

,

(5.42)

 

 

ρ

81

где Φ(ϕ) является собственной функцией оператора "z-компоненты мо-

ˆ

 

 

 

 

 

 

 

 

 

 

мента" Lz = −i~

 

∂ϕ

 

 

 

 

 

 

 

 

 

−i~

 

1

eimϕ

 

 

 

 

Φ = mΦ ,

где Φ =

 

(5.43)

 

∂ϕ

 

 

с собственным значением, принимающим полуцелые значения m =

±12 , ±32 , . . ., что вытекает из требования 2π-периодичности полной функции ψ±(r, ρ, ϕ) (электронные функции ϕ±(r, ϕ) имеют период 4π). "Радиальные" функции χ(±)(ρ) являются решением уравнения ангармонического осциллятора

2ME ∂ρ2 + Vm±(ρ) χ(ρ) = E χ(ρ) ,

~2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~2m2

(5.44)

V ±(ρ) =

1

K

 

ρ2

V ρ +

,

 

 

2MEρ2

m

2

 

 

E

 

± | E |

 

причем потенциальная энергия такого осциллятора отличается от адиабатического потенциала ε± наличием центробежной энергии вращения.

Зависимость Vm±(ρ) представлена на Рис.5.3. Отметим, что наличие минимума Vm+(ρ) при ρ+0 обусловлено только "центробежной силой". В гармоническом приближении для нижних вибронных состояний уравнение (5.44) переходит в обычное уравнение линейного гармонического осциллятора. В этом случае вибронные состояния и вибронный спектр примут вид

 

ψ(±)(r, ρ, ϕ) ϕ±(r, ϕ)

eimϕ χn(ρ − ρ0±)

,

 

(5.45)

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

 

p

 

 

 

 

 

 

2

 

3

 

 

 

 

 

 

 

 

2πρ0±

 

 

 

Emn = ~ωm

n + 2

+

 

 

~ |m| ωm , ωm = 3 ~ m

1

2

,

+

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

VE

 

 

Emn= ~ωE

n + 2

+

2ME0)2 − EJT ,

 

 

 

 

 

 

| |

 

(5.46)

 

1

 

 

 

~2m2

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

(ρ−ρ )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

χn(ρ − ρ0) = Nn Hn(ρ − ρ0) e

0

 

 

 

 

(5.47)

 

2l2

 

 

 

(Hn полином Эрмита, l =

~/M ω

 

) обычная осцилляторная

функция.

 

 

 

 

 

 

p

 

E

E

 

 

 

 

 

 

 

 

 

 

 

Вибронная функция (5.45) имеет наглядный физический смысл. Сомножитель χn(ρ − ρ0) описывает гармоническое колебание величины

82

деформации октаэдра, сомножитель eimϕ вращение волны деформации вокруг центра октаэдра. Электронная волновая функция ϕ±(r, ϕ) адиабатически следует за этой волной деформации. Двукратное вырождение энергетических уровней (Emn ≡ E|m|n ) связано с эквивалентностью вращений по и против часовой стрелки в желобе.

Наиболее прост вибронный спектр нижней ветви адиабатического

потенциала это суперпозиция электронной (−EJT ), "классических"

колебательной (~ωE n + 12 ) и вращательной ([~2m2]/[2ME0 )2]) энергий.

В рассматриваемом пределе сильной вибронной связи для нижних вибронных состояний ρ0 ≈ |VE |/(2KE) и квант вращательной энергии

много меньше кванта колебательной энергии:

 

 

~2

=

2~2KE2

=

~2KE

=

~ωE

~ωE ~ωE .

(5.48)

 

2ME 0)2

ME VE2

ME EJT

EJT

 

Рассмотрим особенности решения вибронной E −e-задачи при наличии внешних возмущений электронных или механических (упругих).

В качестве простейшего электронного возмущения, снимающего вырождение E -уровня, рассмотрим действие низкосимметричного КП,

описываемого эффективным гамильтонианом

 

 

ˆ

σˆz = b (sin ϕ0 σˆx + cos ϕ0 σˆz) ,

 

H1 = b2 σˆx + b3

 

b = q

 

 

 

 

 

 

 

 

(5.49)

b22 + b32

,

b2 = b sin ϕ0 , b3 = b cos ϕ0 ,

 

где b2,3 (b, ϕ0 ) параметры НКП. Оператор типа

ˆ

легко полу-

H1

 

ˆ

 

 

 

 

, Q3

некоторыми

чить из гамильтониана H1 (5.32) полагая в нем Q2

статическими смещениями, тогда b2 = VEQ2(0) , b3 = VE Q3(0) .

Собственные значения гамильтониана

 

 

 

 

 

Hˆ1= Hˆ1 + Hˆ1

 

(5.50)

легко находятся:

 

 

 

 

 

 

ε±(ρ, ϕ) = ± "1 + ρVE

 

+ 2

ρVE cos(ϕ − ϕ0)#

1

 

ρ |VE| . (5.51)

 

 

 

b

 

2

b

 

2

 

 

 

 

 

 

 

 

Так что с учетом линейного вибронного взаимодействия адиабатический потенциал

 

1

 

ε± =

2 KEρ2 + ε±(ρ, ϕ) = ε±(ρ, ϕ)

(5.52)

83

теперь теряет аксиальную симметрию и становится функцией как ρ, так и ϕ. Элементом симметрии потенциала остается только плоскость

ϕ= ϕ0 , ϕ = ϕ0 ± π. "Мексиканская шляпа" деформируется: дно желоба на ветви εподнимается и опускается, достигая экстремальных положений −(EJT ± b) как раз в точках пересечения с плоскостью симметрии "шляпы". Сечение адиабатического потенциала плоскостью

ϕ= ϕ0 , ϕ = ϕ0 ± π, когда

ε±(ρ) =

1

KEρ2 ± (ρ |VE| ± b)

(5.53)

2

(с независимыми вариантами знаков перед b) приведено на рис.5.4 При ρ = 0 расстояние между ветвями потенциала равно 2 |b| вели-

ˆ

ˆ

чине расщепления электронного E -уровня под действием Vнкп =

H1 .

Эта же величина (2 |b|) определяет и расстояние между максимальным и минимальным положениями дна желоба.

При низких температурах система будет находиться в единственном абсолютном минимуме при ρ = ρ0 и ϕ = ϕ0 (или ϕ = ϕ0 ± π, в зависимости от соотношения знаков VE и b). Другими словами, если в отсутствие низкосимметричного электронного возмущения среднее значение hQ2i = hQ3i = 0 и октаэдр ML6 в "среднем" не искажен, то наличие сколь угодно слабого НКП будет приводить к стабилизации определенного искажения комплекса с величиной ρ = ρ0 , вообще не зависящей от НКП (!) и величиной ϕ, определяемой не величиной, а "ориентацией" НКП. Параметр b, характеризующий величину НКП будет обуславливать всего лишь интервал температур kT |b|, в котором эффект стаблилизации искажений октаэдра проявляется наиболее эффективно. Таким образом, мы имеем пример системы с гигантским "электромеханическим" эффектом, когда малое электронное воздействие может приводить к аномально большим упругим деформациям системы.

С другой стороны рассматриваемая нами задача классический пример т.н. псевдоэффекта Яна-Теллера, связанного с проявлением вибронных эффектов для систем с близкими электронными уровнями.

Малое механическое возмущение, действующее на упругую подси-

стему октакомплекса ML6 представим как

 

 

ε = p2Q2 + p3Q3 = p ρ cos(ϕ − ϕ0) ,

(5.54)

где

 

p = q

 

, p2 = p sin ϕ0 , p3 = p cos ϕ0 ,

 

p22 + p32

(5.55)

84

а p2 , p3 величины, определяемые различными компонентами тензора механических напряжений. С учетом возмущения ε адиабатический потенциал линейной E − e-задачи примет вид

 

1

 

ε± =

2 KEρ2 ± |VE| ρ + p ρ cos(ϕ − ϕ0)

(5.56)

искаженной "мексиканской шляпы" с плоскостью симметрии ϕ = ϕ0 , (ϕ0 ± π). В этой плоскости (рис.5.5)

ε±(ρ) =

1

KEρ2 ± (|VE| ± p) ρ

(5.57)

2

(с независимыми вариантами знаков перед p).

Для нижней ветви адиабатического потенциала при этом имеем два

минимума в точках

 

 

 

 

 

 

 

 

 

 

ρ =

|VE | , ±p

= ρ

 

ρ

,

ρ =

p

,

(5.58)

 

 

 

±

KE

0 ±

0

 

0

KE

 

 

ε±) = −

1

KE 0

±

ρ0)2 .

 

 

(5.59)

 

 

 

 

 

2

 

 

Эти минимумы наибольший и наименьший во всем деформированном желобе различаются по энергии на величину

ε

)

ε

(ρ )

= 2K

ρ

 

ρ

= 2 V

ρ

=

2 |VE|

p = 2 p ρ . (5.60)

|

+

 

 

|

 

E

0

0

| E |

0

 

KE

0

Таким образом, как в случае слабого электронного возмущения (|b| EJT ), так и в случае слабого механического возмущения (|p| |VE|), вибронная (ян-теллеровская) E − e-система будет стабилизироваться при низких температурах (kT |b|, или kT |VE | p/KE ) в единственном минимуме адиабатического потенциала при ρ = ρ0 (или ρ =

ρ0 + ρ0 ≈ ρ0 ) и ϕ = ϕ0 , (ϕ0 ± π).

Малые механические и электронные возмущения могут приводить к аномально сильным деформациям октаэдра ML6 .

5.4Линейная вибронная E − (b1 + b2)-задача в квадратных комплексах.

Квадратные комплексы типа ML4 (рис. 5.6) имеют симметрию точечной группы D4h , которой соответствует десять неприводимых представлений a1 , a2 , b1 , b2 , e (соответственно четные и нечетные).

85

Qγν

Симметризованные комбинации

Трансформационные

 

декартовых смещений атомов 1 ÷ 4

свойства

Qa1G

21 (y1 + x2 − y3 − x4)

dz2

Qa2G

21 (x1 − y2 − x3 + y4)

sz, Lz

QeU

21 (−x1 + x2 − x3 + x4)

x, px

 

21 (−y1 + y2 − y3 + y4)

y, py

Qb1G

21 (y1 − x2 − y3 + x4)

dx2−y2

Qb2G

21 (x1 + y2 − x3 − y4)

dxy

Таблица 5.2: Симметризованные смещения лигандов квадратного комплекса ML4

Симметризованные смещения атомов L квадратного комплекса ML4 представлены в табл.5.2.

Рассмотрим вибронные эффекты для комплексов ML4 с двукратно вырожденным электронным состоянием типа E (Eu или Eg ). Учитывая, что для симметричного квадрата

[D(E)]2 = A1 + B1 + B2 ,

(5.61)

мы видим, что активными в вибронном взаимодействии могут быть только Qb1g - и Qb2g -моды смещений (рис.5.6).

Обозначая электронные E -функции, преобразующиеся как px (dxz)

и py (dyz) соответственно ϕ1 и ϕ2

и учитывая, что для электронных

 

 

 

 

ˆ

 

 

 

 

 

ˆ

, преобразующихся как Qb1g ≡ Q1 и Qb2g ≡ Q2 ,

операторов Vb1g и Vb2g

справедливы следующие соотношения

 

 

 

 

Dϕ1

Vˆb1g

 

ϕ1E = − Dϕ2

Vˆb1g

ϕ2E = V1

Vˆb1g = V1σˆz ,

(5.62)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

1

 

Vˆ

 

ϕ

2

 

= ϕ

2

 

Vˆ

 

ϕ

= 0

 

 

D

 

b1g

 

 

 

 

D

 

b1g

 

 

1

E

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ1

 

ˆ

 

 

ϕ1

= ϕ2

 

ˆ

 

ϕ2

= 0

 

 

 

 

 

Vb2g

 

Vb2g

 

Vˆb2g = V2σˆx ,

(5.63)

D

 

 

 

 

 

 

 

E

D

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

1

 

Vˆ

 

 

ϕ

2

= ϕ

2

 

Vˆ

 

ϕ

1

= V

 

 

D

 

b2g

 

 

E

D

b2g

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

представим гамильтониан линейного вибронного взаимодействия для E − (b1 + b2)-задачи в виде

ˆ

= V1Q1 σˆz + V2Q2 σˆx

(5.64)

H1

86

с двумя константами линейной вибронной связи V1 и V2 .

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

 

V1Q1 = ρ VE cos ϕ , V2Q2 = ρ VE sin ϕ ,

 

ρ VE = q(V1Q1)2 + (V2Q2)2, (ρ > 0, VE > 0) ,

(5.65)

ˆ

представим H1 в виде, полностью аналогичном вибронному гамильтониану для линейной E −e-задачи. Собственные функции и собственные

ˆ

значения H1 имеют вид

ψ+ = cos ϕ2 ϕ1 + sin ϕ2 ϕ2 ,

ε+ = +ρVE ,

ψ= − sin ϕ2 ϕ1 + cos ϕ2 ϕ2 ,

(5.66)

ε= −ρVE ,

Адиабатический потенциал имеет две ветви (см. Рис.5.7)

 

 

 

1

 

 

 

 

 

±

(5.67)

ε

 

=

1

K

Q2

+ K

Q2

 

 

ρV =

 

±

 

2

1

1

2

2

 

 

 

E

 

 

 

 

K1Q12 + K2Q22 ± q

 

.

 

 

=

 

(V1Q1)2 + (V2Q2)2

 

 

2

Экстремальные точки поверхности ε(Q1, Q2) расположены симметрично относительно нулевой точки:

(0)

 

V

 

(0)

 

 

 

(0)

(1)

 

V 2

 

1

 

 

 

 

 

1

 

Q1

= ±

 

, Q2

= 0 , ε(Q

 

) = −EJT

= −

 

,

K1

 

2K1

 

 

 

 

 

 

 

 

 

 

 

V 2

(5.68)

(0)

 

 

(0)

 

 

V

 

(0)

(2)

 

 

 

 

 

2

 

= −

2

 

Q1

= 0 , Q2

= ±

 

, ε(Q

 

) = −EJT

 

,

K2

 

2K2

причем если EJ(1)T > EJ(2)T , то экстремумы вдоль оси Q1 являются абсолютными минимумами, а экстремумы вдоль оси Q2 седловыми точками, если же EJ(1)T < EJ(2)T , то наоборот, первые экстремумы седловые точки, а вторые минимумы.

На рис.5.7 изображен адиабатический потенциал E −(b1 +b2)-задачи при EJ(1)T > EJ(2)T ; на рис.5.8 сечение вдоль координаты Q1 .

87

Электронные волновые

Q(0)1 =6 0, Q(0)2 = 0 : ψ

Q(0)1 = 0, Q(0)2 =6 0 : ψ

функции в экстремальных точках имеют вид:

(r, Q(0)) =

( ϕ2, V1Q1 (0)

0

 

 

 

 

 

 

 

 

(0)

>

 

,

 

 

 

−ϕ1, V1Q1 < 0,

 

(r, Q ) =

1

 

 

 

 

 

(0)

 

 

 

1

 

 

 

 

(0)

 

(0)

 

2

1

 

ϕ2) , V2Q2

> 0,

 

2 1 + ϕ2) , V2Q2 < 0,

 

 

 

 

 

 

 

 

 

 

 

(5.69) Оператор квадратичного вибронного взаимодействия в квадратном комплексе с электронным E -уровнем может содержать только слагаемые типа Qa1g Qb1g и Qa1g Qb2g . Их учет при определенных условиях может приводить к появлению четырех минимумов адиабатического

потенциала ε(Q1, Q2).

88