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

книги / Принципы динамической теории решетки

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
19.18 Mб
Скачать

282 Глава 4

1

М

4 sinh2 (tXj + аз) —^ (Pi + Рч)2

А лАо

(4.3.25)

(Pi - Pi)2 ~

4 sinh2 у (*! - <x2)

Знаки корней (4.3.24) можно выбрать четырьмя различными спо­ собами, но лишь два из них дают независимые решения (Р iP2 > 0 и PiP2 < 0)» Здесь мы ограничимся лишь случаем PjP2 > 0* Из (4.3,24) и (4.3.25) находим

А = 2 | / £

8in h £ ,

А =

2]

/ £ s i n h £ ,

]

 

 

 

Jlf

 

2

 

 

|/ Jtf

2

 

 

 

 

/ sinh -4- («1 +

<*г)\

 

 

 

 

 

(4.3.26)

 

 

 

 

 

 

 

А\Аг =

 

 

 

 

 

 

 

 

 

 

\ sinh

(«! — “ г) /

 

 

 

 

 

 

Рассмотрим поведение функции xpt

в области, где <p j = осxl - р xt =0.

При этом

 

 

 

 

 

 

 

 

 

(р2 =

(Х2£ — p2t =

— 9l +

 

 

 

 

 

(4.3.27)

где

 

 

«I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.3.28)

 

*1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поскольку рг/р 2 = sh(a1/2 )/sh (a a/2 ) > a^/ajj при ocj >

a2 > 0. Таким

образом, выражение (4.3.23) принимает вид

 

 

 

y>i

1 +

exp (<*jZ — pxt)

при t ->■ —с»,

1

 

(4.3.29)

y>i cz А2 -f exp (ocxl — pxt)

при t -> +oo. J

 

 

 

Аналогично в области, где <р2 = a2l - Р2* 550, имеем

 

У/ ^

i4j + exp (a2l p2t)

при t -> — оо, |

 

(4.3.30)

У/ ^

1 -^2exp («2Я— р2*)

при *

+оо. j

 

 

 

В соответствии с (4.3.17), (4.3.29) и (4.3.30) (ср. (4.3.18) и (4.3.20))

можем записать, что асимптотически

 

 

 

 

 

4аЬ

 

2

 

 

 

< =

(1)

при

« - - « о ,

 

 

 

 

 

 

w

 

(4.3.31)

e-^i — 1 =

рл aech2-i- («,г —

 

Т й+) ,

г =

при

I '-*■+ оо,

Ангармонические кристаллы и структурные фазовые переходы

283

где

и

<5+ =1пЛ2.

(4.3.32)

д~ = In А г

Таким образом, при *-»-«> это решение описывает два солитона, расположенных так, что медленный солитон (2) находится впереди быстрого солитона (1), а при t -> + *> они расположены в обратном по­ рядке (рис. 4.13). Подразумевается при этом, что солитоны устой­ чивы по отношению к столкновениям15.

В случае р,р2 < 0 солитоны движутся в противоположных направ­ лениях. За более детальным рассмотрением этого случая и /V солитонной проблемы мы отсылаем читателя к работам Тоды [390, 3911 и цитированной там литературе.

Получим теперь периодические решения для решетки Тоды. При этом мы используем следующую теорему для эллиптических функций

Якоби:

d

sn и сп. и dn и sn2v

 

 

 

8 П 2 +

v) — sn2(и — v) = 2 dv

1 к2sn2и sn2v

(4.3.33)

где эллиптические функции Якоби sn, сп и dn определяются соотно­ шениями (об используемых ниже свойствах этих функций см., напри-

Рис. 4.13.

Двухсолитонное решение для решетки Тоды при Р1Р2 > 0 (см.

текст). а)

при t - ° р ; б) t г* + °°;. 1 - быстрый солитон, 2 — медленный

солитон.

 

15 Интересно отметить, что численные методы, использованные для ис­ следования столкновений уединенных волн в одномерной цепочке с потен­ циалом Леннард-Джонса (см. (1.4.46)) привели к результатам, которые ана­ логичны полученным для решетки Тоды [ 398].

284

Глава 4

 

мер, приложение 1 в [ 389])

9

г. г

 

= Г de

= ( __

d#

J Vl — к2 sin2 &

J У(1 ж2) (1 - к2х2)

о0

snw = sin 97 = sn (и, к) ,

спи = cos= сп (м, А:),

dn и = y i — Р sn2г* = dn (u, A;).

i

Интегрируя (4.3,33) при учете, что dn2u = 1 — k 2s n 2u

= - &2sn u

cn и, получаем

«"(«)

€(u +

v) + t(u v) — 2€(u) =

(1 /sn2 v) — 1 -f-

где мы ввели функцию

и

(4.3.34)

и d(dn и ) / du =

(4.3.35)

€(и) =

J dи' dn2 и' .

(4.3.36)

 

о

 

Дзета-функция Якоби Z (и), или zn и, определяется как

 

Z(u) =

€(и) - Еи/К

(4.3.37)

и является периодической функцией с периодом 2К* Здесь К и Е - пол­ ные эллиптические интегралы

я/2

d6>

 

= Г-==

 

к = Я(*) = Г

0

(4.3.38)

J yi -

fc2 sin2

J

У(Г^ х2) (1 - ifc2s 2)

О

 

 

о

 

 

я/2

 

 

 

1

_______

Е = Е(к) = Jd0 yi — sin2 0 =

JАх

(4.3.39)

О

 

 

 

о

 

С помощью (4.3.35) и (4.3.37) находим, что

 

(l/sn*«)------------- ------------------------1 + EIK

+ Z'(u)

= Z(u + v) + Zlu — v) — 2Z(u).

 

^

+

' 840)

Чтобы сравнить уравнение (4.3.40) с уравнением движения для ре* шегки Тоды, запишем

3,(t) = j di's,(t')

(4.3.41)

Ангармонические кристаллы и структурные фазовые переходы

285

и продифференцируем (4.3.13) по времени. В результате получим

^ In (1 +

s{/а) =

(*/-1+ */м — &«)•

(4.3.42)

Далее, с помощью (4.3.41) запишем (4.3.12) в виде

 

е"Ьг*— 1 = *,/«.

 

(4.3.43)

Сравнивая (4.3.40) и (4.3.42), получаем

 

*I =

2 J M ( U),

 

(4.3.44)

и =

2(11). ±

vt) К ,

 

(4.3.45)

о =

2КIX,

 

 

(4.3.46)

(l/(sn2v) -

1 + EjK) (2K vf

= 1.

(4.3.47)

Используя (4.3.43), (4.3. 44) -

(4.3.47) и (4.3.37), окончательно

находим решение (так называемую кноидальную волну)

 

е-ьи _ 1 =

м (2K vf {dn2[2(//Я

rf) X] - В Д ,

(4.3.48)

 

 

ab

 

 

где частота v и длина волны Л удовлетворяют дисперсионному соот­ ношению (4.3.47)

2Kv =

<ab

{

1

 

я ^-1'2

(4.3.49)

Л/ \sn2(2KjX)

+ к )

 

 

Функция dn2(2xK) - Е/К является периодической функцией от х

с периодом, равным 1. Ее фурье-разложение имеет вид

 

, 2/0

 

 

Е

л2

~

п cos 2лпх

(4.3.50)

dn (2хК)

 

 

к -

Ka ^

ainh (nnK'jK )’

К' = К { р

-

к2).

 

 

 

(4.3.51)

Для малых k

 

величина К'

очень большая, и приближенно можно за­

писать

 

 

 

 

 

 

 

(4.3.52)

8П X C=Lsin X.

В этом случае кноидальная волна (4.3.48) сводится к косинусоидальной

286

Глава 4

волне в гармонической решетке:

 

 

(4.3.53)

со = 2 YabjM sin (я/А).

(4.3.54)

Таким образом, параметр k определяет амплитуду волны. Если же мы используем тождество

Ш к " <4-3-55>

то кноидальную волну можно представить в виде

где

(4.3.57)

<х = пК!ХК,1 p=nKvlK\

Это означает, что кноидальная волна представляет собой беско­ нечную последовательность одинаковых s e c h 2 - импульсов (солито-

нов); рис. см. 4.14. Солитоны при этом не являются независимыми друг от друга. Их скорость определяется дисперсионным соот­ ношением (4*3.49), а не формулой (4.3.21). В пределе Л -> k -» 1 и а бесконечно волновая последовательность сводится к одиночному солитону на бесконечной решетке.

С помощью численных экспериментов Накамуре [ 276] удалось за­ метить, что примесный атом в решетке Тоды может вызывать появ­ ление довольно хорошо определенной (незатухающей ) локализованной моды. Локализация оказывается очень сильной в случае сильной свя­ зи с примесью; в то же время она ослабляется при малой примесной массе [ 430]. Кроме того, при некоторых условиях падающий на при­ месь солитон может служить источником отраженной и / или прошед-

dnz(2xK) - £ (kz-0,99Z1)

х

Рис. 4.14. Форма кноидальной волны.

Ангармонические кристаллы и структурные фазовые переходы

287

Рис. 4.15. Рассеяние солитонов на примеси в решетке Тоды в случае сильной

связи с примесью. V^ (е ^ * - 1 ).а^ {1 + (а' - 1) (5^ + 5 ^ } ,b^nb {1+

+ (Ь'- 1 ) (6/о + 8 ^ )} ,а' = Ь' = 5 (см. (4.3.7)); а) Д= 4, б) а= 8 (см. (4.3.7)).

Амплитуды нормированы на амплитуду падающего солитона. Время измеря­ ется в единицах {M /a b )1/ 2 . Для очень крутого фронта падающего солитона (большие (X) возбуждение локализованных мод оказывается менее эффек­

тивным. (Согласно [276] .)

шей уединенных волн, которые сопровождаются волнами "ряби" (рис. 4.15). Отметим также другие актуальные проблемы, это - двухатомная ре­ шетка Тоды [266] и статистическая механика решетки Тоды (114,263].

4 .3 .3 . Солитоны и фононы в модели двух потенциальных

ям для структурны х фазовых переходов

Для описания структурных (а также несобственных сегнетоэлектрических) переходов широко используется однокомпонентная d-мерная модель, определяемая следующим классическим гамильтонианом:

288 Глава 4

н ~ ?

Ш + f [ т “ ,г + Т М'4] + Т £

{щ ~ Щ,)2 > (4.3.58)

с л < 0; В,

С > 0 (см ., например, [ 32, 72, 73, 340]). Здесь иг и pt -

смещение и импульс атома, расположенного в /-м узле регулярной

(простой кубической для d = 3) решетки, М -

атомная масса и третье

суммирование в (4.3.58) осуществляется лишь по ближайшим сосе­ дям. Этот гамильтониан соответствует набору связанных гармони­ ческих осцилляторов с идентичными одночастичными потенциалами, представляющими собой две потенциальные ямы (рис. 4.16). Он мо­ делирует кристалл с двумя подрешетками, где атомы одной подре­ шетки (предполагается, что она жестко зафиксирована) создают двухъямный потенциал для подвижных атомов другой подрешетки.

Как можно показать, модель (4.3.58) описывает систему, в ко­ торой (в случае d > 1) при ненулевой температуре Тс происходит фазовый переход в упорядоченное состояние с конечным средним смещением (< иг > ф 0). Здесь возможны два предельных случая - это переход порядок — беспорядок и переход типа смещения. Разни­ цу между этими предельными случаями легко понять с помощью изо­ браженной на рис. 4.16 модели. При низких температурах все подвиж­ ные атомы расположены на дне, скажем, левой потенциальной ямы. С ростом температуры возможна, однако, реализация двух различ­ ных ситуаций. В одном случае наиболее вероятное положение ато­ мов соответствует вершине потенциального барьера (переход типа смещения), а в другом случае оно отвечает дну потенциальной ямы, в результате чего ямы оказываются заполненными равновероятно (переход порядок - беспорядок).

Деформационные переходы можно классифицировать на основе от­ ношения f = EQ/ EC где Е0 = А2/4В и Ес ЛС\ А\/ В обозначают соот-

Рис. 4.16. Схематическое изображение модели (4.3.58) дл" d " 1* (Согласно

[73] .)

Ангармонические кристаллы и структурные фазовые переходы

289

ветственно глубину ямы (высоту потенциального барьера) и энергию взаимодействия ближайших соседей в разных потенциальных ямах (рис, 4Л7 и рис. 4.18) (минимум двухьямного потенциала соответству­ ет смещению ±и0 = ±у/\ Л |/£). Отметим, что в случае d = 3 отноше­ ние kTc /Ec порядка единицы [330]. При условии £ >> 1 (предельный случай перехода порядок - беспорядок) каждый (подвижный) атом ока­ зывается локализован вблизи дна потенциальной ямы при всех темпе­ ратурах, кроме Т » Тс . Таким образом, в гармоническом приближе­ нии все колебания вблизи высокотемпературного положения равнове-. сия (т.е. около вершины потенциального барьера) неустойчивы. В дан­ ном предельном случае основные динамические процессы - это прыж­ ковые процессы между соседними ямами. Такая ситуация может быть описана псевдоспиновой моделью, а при достаточно высоких темпера­ турах — моделью, соответствующей невзаимодействующим двухъямным осцилляторам. В пределе £ « 1, т.е. перехода типа смещения, неустойчива небольшая группа длинноволновых колебаний вблизи

Рис. 4.17. Распределение плотности вероятности для координат частицы Р (и)

в случае перехода порядок - беспорядок (Е0/£ с = 6,25) в модели (4.3.58) при d = 1 и различных темпаретурах. (Согласно [20] .)

290

Глава 4

Рис. 4.18. Распределение плотности вероятности Р (и) для координат чатицы в случае перехода типа смещения в модели (4.3.58) при d = 1 и различных тем­ пературах; Сплошная кривая: к Т /Е 0 - 0,5; штрихпунктирная: к Т /Е 0 = 2 ; штриховая: к Т /Е 0 - 10. (Согласно [20] .)

высокотемпературного положения равновесия. Для описания этого предельного случая оказывается продуктивной концепция мягкой фо­ нонной моды.

Здесь мы ограничимся рассмотрением предельного случая £ « 1 перехода типа смещения в одномерной системе [ 20, 224, 399]. Посколь­ ку при £ « 1 разность |их + J относительно мала, мы предпола­ гаем, что гамильтониан (4.3.58) можно записать в континуальном представлении:

'“ / x f

=j дхЩр(х), и(х), ди(х)1дх)у

где L - постоянная кристаллической решетки, xt = IL == х - ная пространственная переменная и с% = 2L2C/M.

Из канонических уравнений

(4.3.69)

непрерыв­

ээе

д

гэе

ди

дх д(ди/дх) 1

(4.3.60)

и плотности гамильтониана К (4.3.59) получаем следующее уравнение движения для поля смещений и(х):

Мй + Аи + ВигМс02и" = 0.

(4.3.61)

Ангармонические кристаллы и структурные фазовые переходы

291

Если нас интересуют решения (4.3.61) с постоянным профилем, т.е. решения, зависящие лишь от одной переменной х - vt с постоян­ ной скоростью v, то (4.3.61) можно записать в виде

§

+ п -

= 0;

 

(4.3.62)

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

ulu0 = rj,

 

 

(4.3.63)

vt)/£ = в,

 

 

где

М(с02 — v2)l\A\

(квадрат длины).

(4.3.04)

f 2 =

Рассмотрим вначале решения уравнения (4.3.62), описывающие

колебания с малой амплитудой. При г|а << r\«

1 (4.3.62) сводится

к уравнению

 

 

 

1?" + V =

0,

 

(4.3.66)

которое имеет решения вида

 

г) =

осsin (s + 0),

 

(4.3.66)

где а -

амплитуда и © -

фаза. Используя (4.3.66) и (4.3.63), полу­

чаем

<хщsin [(а? — vt) f-i + 0].

 

и =

(4.3.67)

Это — фонон с волновым вектором k = g-1, частотой v/ § и фазовой скоростью v. Из (4.3.67) и (4.3.64) находим дисперсионное соотно­ шение

о.>2(к) = v2k2 = с02к2 + А/М.

(4.3.68)

Поскольку А < 0, частота со(А) действительна лишь при к > (| А\/ЪЛс*)Ъ Решение (4.3.67) описывает колебания с малой амплитудой около положения и = 0.

Колебания с малой амплитудой около ± и0

можно получить, под­

ставляя

 

*7 = ± 1 + У

(4.3.69)

в (4.3.62) и сохраняя лишь линейные по у члены. В результате имеем уравнение

V" ~ = 0»

(4.3.70)

с решениями вида

 

у = sin }/28 -f 0).

(4.3.71)

Соседние файлы в папке книги