Влияние эффектов дальнодействия на мультикритическое поведение неупорядоченных сжимаемых систем Текст научной статьи по специальности «Физика»

Математические структуры и моделирование 2005, вып. 15, с. 55-65

УДК 536.763/764

ВЛИЯНИЕ ЭФФЕКТОВ ДАЛЬНОДЕЙСТВИЯ НА МУЛЬТИКРИТИЧЕСКОЕ ПОВЕДЕНИЕ НЕУПОРЯДОЧЕННЫХ СЖИМАЕМЫХ СИСТЕМ

С.В. Белим

В рамках теоретико-полевого подхода рассмотрено мультикритическое поведение систем, описываемых двумя флуктуирующими параметрами порядка, с замороженными дефектами структуры при различных значениях параметра дальнодействия непосредственно в трехмерном пространстве. Выявлены устойчивые фиксированные точки ренормгруппового преобразования. Определены типы мультикритического поведения при различных значениях параметра дальнодействия.

Как хорошо известно, критические свойства систем задаются малым количеством параметров. Таких, как размерность, симметрия параметра порядка и скорость убывания взаимодействия с расстоянием. Особый интерес представляют системы, в которых кроме обычного блнзкодействня присутствуют эффекты дальнодействия. В классической изингоподобной системе взаимодействие между флуктуациями убывает экспоненциально с расстоянием по закону ехр(—г/го), в связи с чем рассматривается взаимодействие только между ближайшими соседями, в связи с чем данные системы можно охарактеризовать как близкодействующие. При убывании взаимодействия с расстоянием г по закону r~D~a) где D - размерность пространства, уже нельзя ограничиваться взаимодействием между ближайшими соседями, и возникают эффекты дальнодействия. Как показано в работе [1] в рамках ^-разложения (є = 2а — D) для значений параметра дальнодействия в интервале 1,5 < а < 2 наблюдается негауссово критическое поведение, отличное от критического поведения близкодействующих систем и существенно зависящее от параметра дальнодействия а. В интервале значений а < 1,5 в сисеме наблюдается гауссово критическое поведение. Как показали расчеты непосредственно в трехмерном пространстве для критического поведения однородных систем с эффектами дальнодействия подтверждается предсказание ^-разложения [2]. Введение в систему замороженных примесей существенно изменяет картину критического поведения [3]. Так, для неупорядоченных систем так же, как и для однородных систем при значениях параметра дальнодействия а < 1,5, наблюдается гауссово критическое

Copyright © 2005 С.В. Белим.

Омский государственный университет.

E-mail: belim@univer. omsk. su

Работа поддержана грантом РФФИ N 04-02-16002.

56

С. В. Белим Влияние эффектов дальнодействия...

поведение. Однако интервал 1,5 < а < 2 разбивается на два. При значениях 1,8 < а < 2 критическое поведение носит негауссовый характер, существенно зависящий от параметра дальнодействия а. Однако при 1,5 < а < 1,8 в системе происходит срыв на фазовый переход первого рода. Упругие деформации приводят как к смене режима критического поведения, так и к появлению на фазовой диаграмме вещества трикритических линий [4]. Однако интервалы значений параметра дальнодействия, в рамках которых критическое поведение качественно носит одинаковый характер, остаются неизменными.

Как было показано в работе [5], в однородных системах, описываемых двумя флуктуирующими параметрами порядка, эффекты дальнодействия могут приводить к изменению режима мультикритического поведения. Для значений параметра дальнодействия 2 > а > 1,6 наблюдается бикритическое поведение, тогда как для 1,5 < а < 1,6 поведение становится тетракритическим. В интервале значений параметра порядка а < 1,5 устойчивой становится гауссова критическая точка.

Присутствие замороженных точечных дефектов структуры приводит к изменению режима поведения близкодействующих систем как в бикритической, так и в тетракритической области [6]. В указанной работе показано, что влияние 5-кореллированных примесей приводит к развязыванию параметров порядка в мультикритических точках. В работе [7] выявлено, что упругие деформации для однородных систем приводят к смене бикритического поведения тетракритическим. Для неупорядоченных систем [2] деформационные степени свободы, не меняя типа мультикритического поведения, изменяют режим тетракритического поведения.

Предметом данной статьи является исследование влияния эффектов дальнодействия на неупорядоченные сжимаемые системы, описываемые двумя параметрами порядка.

При структурных фазовых переходах с отсутствием пьезоэффекта в парафазе упругие деформации играют роль вторичного параметра порядка, флуктуации которого в большинстве случаев не являются критическими [8]. В связи с тем, что в критической области основной вклад в стрикционные эффекты дает зависимость обменного интеграла от расстояния, рассматриваются лишь упруго-изотропные системы.

Репличный гамильтониан неупорядоченной сжимаемой системы с эффектами дальнодействия имеет вид:

пп

m

Математические структуры и моделирование. 2005. Вып. 15.

57

+2-U03 / dDqidDq2dDq3 ^ (Фд1Фд2)(Ф5зФ-«г1-«г2-?з)

** а, 6=1

/III

dDq3dDq3dDq3 £(Ф^Фу (Ф^Ф*,,.,^) -

a=1

/ТП

АЛЛ - (1)

а=1

/га

dDqidDq3dDq3 ^Ф^Н^Ф^.,,) +

а=1

/ш « m

dDq\dDq2yq\ ^52^-91-92 + ^2 / dD q\dD qyyq\ Фд2^-ді-д2 +

а=1 ^ а=1

о /* m n° Г m С Я

+ ПУ0 / ^5^Ф2Ф-9 + #ї/о / <Рд^Ф“Ф“, + 2/? / dDqyqy-q + 2^УІ

^ а=1 ^ а=1 ^

где Ф и Ф — m-мерные флуктуирующие параметры порядка, г^оі и Щ)2 — положительные константы,

ті ~ Т - ТіІ/Ть

7-2 ~ Т - Тс2І/Тс2,

ГС1 и Тс2 — температуры фазового перехода для первого и второго параметра порядка соответственно,

з

Щ) =

ск=1

где иар — тензор деформаций, щ и д2 — параметры квадратичной стрикции, (3 — постоянная, характеризующая упругие свойства кристалла, D — размерность пространства. Свойства исходной системы могут быть получены в пределе m —> 0. Неотрицательные константы 50i, 5q2, ^оз описывают взаимодействие критических флуктуаций через поле примесей. Взаимодействие примесей с упругими деформациями носит линейный характер и при усреднении по примесям приводит к переопределению констант #01, $021 $оз- В данном гамильтониане уже проведено интегрирование по слагаемым, зависящим от нефлуктуирующих переменных, не взаимодействующих с параметром порядка. В (1) выделены слагаемые у^) описывающие однородные деформации. Как показано в работе [9], такое разделение необходимо, так как неоднородные деформации yq отвечают за обмен акустическими фононами и приводят к эффектам дальнодействия, которые отсутствуют при однородных деформациях.

Определим эффективный гамильтониан системы, зависящий только от сильно флуктуирующих параметров порядка Ф и Ф, следующим образом:

ехр{-Я[Ф,Ф]} = Я J ехр{—Я0[Ф, Ф, 2/]} JJ dyq,

(2)

58

С. В. Белим Влияние эффектов дальнодействия...

Если эксперимент осуществляется при постоянном объеме, ТО у о является константой, интегрирование в (2) проводится только по неоднородным деформациям, и однородные деформации вклада в эффективный гамильтониан не вносят. При постоянном давлении в гамильтониан добавляется слагаемое РП, объем представляется в терминах компонент тензора деформации в виде

П = П0[1 + uaa + uaauyy + 0(u3)\ (3)

OL=l ОіфУ

и интегрирование в (2) осуществляется также и по однородным деформациям. Как отмечено в [10], учет в (3) квадратичных слагаемых может оказаться важным в случае высоких давлений и кристаллов с большими стрикционными эффектами. В результате:

-І Г m Л Г m

Н=2 d°q^ + & Е ф.ф% + 2 /dDq{j2 + q2) Е ф.ф% +

** а=1 ^ а=1

/га

dDqidDq3dDq3 V (Ф^ХФ^Ф-,!-*!-*.) +

а,6=1

/га

dDqidDq3dDq3 Y, (ФХФ^)(Ф»зФ-»1-»2-»з) +

~ U_1

(2,6—1 171

+2^03 / dDqidDq2dDq3 (Ф“іФ“2)(Ф5зФ-?і-?2-?з) ^ (2,6=1

С /> ™

--S- dDqldDq.1dDq3Y1{%3%3)i^-,

- —ql—qc2—q?y)

(2=1

m

—? / ^?1^22^9зЕ(Ф?1Ф92)(ФдЗФ%1-д2-д3) -

J (2=1

/га

dDqldDq2dDq3 Х^ФУО^Ф %-*!-«>) +

(2=1

2 2/* 171

+£Чг1 / <*с<гЛ ^(Ф^Ф-зіМФ^Ф-й) +

J a,6=1

m

2 , іС«і<іі’«2^(Ф:іФ“„1)(Ф“2Ф-»2) +

4 - tt>2

+(^2-вд) / 5](Ф“1ФД)(ф;2Ф

b ) -q2J>

(2,6=1

l>01 = «01 - z\/2, v02 = u02-zl/2, V03 = U03 - Z1Z2/2, _ 91 _ 92 _ 9i _ g2

Zl VT 22 VT

(4)

Данный гамильтониан приводит к широкому разнообразию мультикритиче-

Математические структуры и моделирование. 2005. Вып. 15.

59

ских точек. Возможно как бикритическое

(«з + CiZ2 - wxw2 - <5з)/2)2 > Д + (zf -w\- 5і)/2)(г>2 + ~w\~ 82)/2),

так и тетракритическое

(V3 + (Z1Z2 - Wiw2 - (53)/2)2 < (vi + (z\ -w\- 8i)/2){v2 + (z2 -w\- 52)/2)

поведение. В первом случае в мультикритической точке пересекаются две линии фазовых переходов второго рода и одна линия фазовых переходов первого рода, во втором — четыре линии фазовых переходов второго рода. В непосредственной окрестности мультикритической точки система демонстрирует специфическое критическое поведение, характеризующееся конкуренцией типов упорядочения. При этом в бикритической точке происходит вытеснение одного критического параметра другим, тетракритическая же точка допускает существование смешанной фазы с сосуществующими типами упорядочения. Кроме того, стрикционные эффекты могут приводить к мультикритическим точкам более высокого порядка.

В рамках теоретико-полевого подхода [11] асимптотическое критическое поведение и структура фазовых диаграмм во флуктуационной области определяется ренормгрупповым уравнением Каллана-Симанчика для вершинных частей неприводимых функций Грина. Для вычисления /3- и д-функций как функций, входящих в уравнение Каллана-Симанчика, перенормированных вершин взаимодействия щ, гд2, Щ, Ф, 52, 53, 9ъ 92, 9^\ 9^ или более удобных для определения мультикритического поведения модели комплексных вершин щ, z2, гщ, гс2, щ, r>2, г>3, ф, 52, 53 был применен стандартный метод, основанный на диаграммной технике Фейнмана и процедуре перенормировки [12]. В результате в рамках двухпетлевого приближения были получены следующие выражения для /5-функций:

Дл = —щ + 36щ + Фщ — 24щф — 1728 — 1 — - GJ vf —

- 192 ^2J- 1 - g'J V!vl - 64(2 J- l)dj +

+ 96 ^2J — 1 — — G^j V1V3S2, + 32(2J — +

+ 2304(2J - 1 - 1 G)v\81 - 672 ^2 J- 1 - G/j vx Sf + 16(2 J- 1)ЩЬ Pv2 = —v2 + 36v% + 4г| — 24^2 — 1728 ^2 J — 1 — ^ G^j v\ —

- 192 ^2 J — 1 — ^ G^j v2v\ — 64(2 J — 1)г>| + 96 ^2 J — 1 — G^j v2v353 +

+ 32(2J - 1)Щ3 + 2304 (2J - 1 - 1 G^j v\82 ~

- 672 \ 2J — 1 — G j v28\ + 16(2J — 1)Щ2,

(5)

—Z\ + 24^1^1 -|- 2— \Q8\Z\ — 4$$Z2 2^і^2 —

576 (2J- 1 - -бА v\zx - 32(2J - \)v\zi - 16(2J- 1 - G)vjz2+

Xz

I +

+ 1 Oo Ov co II 1 + + Jo Oo to

to 00 00 О 4o 1 On co 1—1 CO to or 1—1 to 00 1 On to

4o ^) —i- N N 4o +

1 1 00 to О to О 1 1—1 C7)

1—1 1—^ co to + 1 1—^ 1 1—^ 1—1 On to to

to to On со со со + к 4 1—1 to C2 1 CO 1 Ю 1 CO 1 Ю C2 co to On co 1 to

On

to to

+

jo

J*

+

I + +

+

05

О

С. В. Белим Влияние эффектов дальнодействия..

Математические структуры и моделирование. 2005. Вып. 15.

61

- 32(2 J - l)z26% + 128(2 J - 1 - bG)v3zx83 + 128(2J - 1 - 5G)v3z283, f3z2 = ~z2 + 2Av2z2 + 2z\ — 1682z2 — A83z\ + 2z2^i 4” ~

- 32(2J - l)Zl6% + 128(2J - 1 - 5G)v3z283 + 128(2J - 1 - 5G>3^3, f3w\ — —W\ 2Av3w3 З- 2z\w\ — 2— 108\W\ — A83w2 2w\Z<2

Известно, что ряды теории возмущений являются асимптотическими, а вершины взаимодействия флуктуаций параметров порядка во флуктуационной области достаточно велики, чтобы можно было непосредственно применять выражения (5). Поэтому с целью извлечения из полученных выражений нужной физической информации был применен обобщенный на многопараметрический случай метод Паде-Бореля. При этом прямое и обратное преобразования Боре-ля имеют вид

/Д, 1>2,1>з, 8г,82,83l Zi, z2l Wi,w2) =

+ 128(2J - 1 - bG)v3wx83 + 128(2J - 1 - 5G)v3w283,

(3W2 = ~w2 + 2Av2w2 + 2z%w2 - 2w\ - 16<52w2 - A83wi + 2w2z\ +

+ 128(2J - 1 - 5G)v3w283 + 128(2J - 1 - bG)v3wx83.

...no^i1 vl2Д8\4Д8г36z\7zi8w[9wl10 =

h v Лю

(X)

0

F(vi, v2l v3, 8i,82, 83, zuz2, wuw2)

(6)

62

С. В. Белим Влияние эффектов дальнодействия...

- £

Чт--МО

І1.-.І7--- у*2 Дз ДА £*6 Ц > Ц ЦЮ _

Г (Іі +... + гіо)! 1 2 3 1 2 3 1 2 1 2

Для аналитического продолжения борелевского образа функции вводится ряд по вспомогательной переменной 9:

F(vі, v2, v3, 5Ъ 62, S3, z1} z2, w1} w2, 9) = (7)

(X)

/c=0 , • • • ,^io

к которому применяется аппроксимация Паде [L/M] в точке 9 = 1.

В двухпетлевом приближении для вычисления /5-функций был использован аппроксимант [2/1]. Природа критического поведения определяется существованием устойчивой фиксированной точки, удовлетворяющей системе уравнений:

рг(в*г^*2,в*,5*г,5*2,5*,

Zi,z2,w1,w2

= О

(* = !,.., 10). (8)

Требование устойчивости фиксированной точки сводится к условию, чтобы собственные значения bi матрицы

—

ад(щ*, vl 61 6%, 61 щ*, zZ, wl wl)

dvj

(S)

/ _ * * * c* c* c* * * * * \

[Vi.Vj = v1,v2,v3,d1,d2,d3,z1,z2,w1,w2)

лежали в правой комплексной полуплоскости.

Полученная система просуммированных /5-функций содержит широкое разнообразие фиксированных точек, лежащих в физической области значений вершин (vi > 0, Si > 0, і = 1,2,3). Полный анализ фиксированных точек, соответствующих критическому поведению только одного параметра порядка, приведен в работе [3]. Фиксированные точки однородной системы, описываемой двумя параметрами порядка, с эффектами дальнодействия были получены в работе [5]. Фиксированные точки, отвечающие за мультикритическое поведение сжимаемых однородных, приведены в Таблице 1, неупорядоченных систем — в Таблице 2. В обеих таблицах введен дополнительный параметр

p=(v3 + (ziZ2 - wxw2 - <53)/2)2 -

-Д + (z\ -w\- ф)/2)(у2 + (z2 -w\- 82)/2),

определяющий тип мультикритического поведения. При значениях р > 0 в системе наблюдается бикритическое поведение, в обратном случае р < 0 - тетра-критическое.

Анализ значений фиксированных точек однородных сжимаемых систем с эффектами дальнодействия и их устойчивости позволяет сделать ряд выводов. Для всех значений параметра дальнодействия 1, 5 < a < 2 на фазовой диаграмме вещества может наблюдаться тетракритическое поведение (р < 0), для этого

Математические структуры и моделирование. 2005. Вып. 15.

63

Таблица 1. Значения фиксированных точек и собственных значений матрицы устойчивости

однородных систем.

N Ц, «2> Vl Щ 5 ^2 «У W2 Й) &2, h h, b5 ^6) У V

а = 1,6

1Д 0,027427 0,224319 0,111301 0,157 1,118 1,798 -0,002092

0,027427 -0,224319 -0,111301 0,738 0,138 0,256

0,026699 0,919

1,2 0,027427 0,224319 0,224319 0,157 1,118 -1,050 -0,000039

0,027427 -0,224319 -0,224319 0,738 0,138 -1,188

0,026699 0,919

а = L,7

2,1 0,031287 0,248013 0 0,113 1,675 2,546 -0,003849

0,031287 -0,248013 0 0,629 0,095 0,041

0,031334 0,809

2,2 0,031287 0,248013 0,248013 0,113 1,675 -1,580 0,000039

0,031287 -0,248013 -0,248013 0,629 0,095 -1,675

0,031334 0,809

а = L,8

3,1 0,033682 0,266919 0 0,090 1,831 2,980 -0,004802

0,033682 -0,266919 0 0,571 0,104 0,115

0,034575 0,753

3,2 0,033682 0,266919 0,266919 0,090 1,831 -1,954 0,000061

0,033682 -0,266919 -0,266919 0,571 0,104 -2,079

0,034575 0,753

а = L,9

4Д 0,035842 0,297071 0 0,069 2,079 3,765 -0,079943

0,035842 -0,297071 0 0,505 0,125 0,049

0,039202 0,702

4,2 0,035842 0,297071 0,297071 0,069 2,079 -1,954 0,000252

0,035842 -0,297071 -0,297071 0,505 0,125 -2,079

0,039202 0,702

константы, характеризующие стрикционное взаимодействие флуктуирующих параметров порядка с деформационными степенями свободы, должны быть разного знака. Для систем, характеризующихся стрикционными константами одного знака, устойчивых фиксированных точек не существует, что свидетельствует о срыве на фазовый переход первого рода. Устойчивые фиксированные точки (ФТ), описывающие тетракритическое поведение системы в интервале значений параметра порядка 1,6 < а < 2 (ФТ 2.1, 3.1, 4.1), характеризуются нулевыми значениями эффективных зарядов и свидетельствующих о том, что данная тетракритическая точка должна наблюдаться при постоянном объеме системы. При значении параметра дальнодействия а = 1,6 эффективные заряды и тої, отличны от нуля (ФТ 1.1). По-видимому, данное отличие

64

С. В. Белим Влияние эффектов дальнодействия...

Таблица 2. Значения фиксированных точек и собственных значений матрицы устойчивости

неупорядоченных систем.

N Ц, «2> Vl а* л* л* ? °2> °3 b1, Ьз &4) ^5, Ьб

а = 1, 6

1,3 0 -0,630515 11,225 22,854

0 -0,630515 15,069 22,854

0 -3,713534 15,069 17,439

а = 1,7

2,3 0 -0,320512 2,239 11,106

0 -0,320512 2,239 11,116

0 -1,932190 5,197 9,250

а = 1, 8

3,3 0 -0,256437 2,068 8,642

0 -0,256437 2,068 8,642

0 -1,545749 3,890 7,441

а = 1, 9

4,3 0 -0,223270 1,224 3,149

0 -0,223270 1,338 7,349

0 -0,790114 1,338 7,349

обусловлено тем, что в отсутствие упругих деформаций система демонстрирует тетракритическое поведение при а = 1,6 и бикритическое при остальных значениях параметра дальнодействия 1, 6 < а < 2. Кроме тетракритических точек упругие деформации приводят к возможности появления на фазовой диаграмме вещества критических точек более высокого порядка. Так, фиксированные точки 1,2, 2,2, 3,2, 4,2 характеризуются равенством эффективных зарядов щ = w\ и Д2 = ^2, то есть в них пересекаются по две линии трикритического поведения, вследствие чего они соответствуют критическим точкам четвертого порядка. Данные фиксированные точки не являются устойчивыми в рамках выбранного приближения, однако для полного разрешения вопроса об их устойчивости необходим учет в исходном гамильтониане кубических слагаемых по деформационным степеням свободы.

Как хорошо видно из таблицы 2, фиксированные точки неупорядоченной «жесткой» системы лежат в нефизической области значений (фф2 < 0). Более того, все устойчивые фиксированные точки имеют достаточно большие значения по модулю, для которых неприменимы положения используемой теории. Отсюда можно сделать вывод о неустойчивости всех фиксированных точек, лежащих в физической области значений эффективных зарядов, и, как следствие, об отсутствии на фазовой диаграмме вещества мультикритических точек. Нахождение устойчивых фиксированных точек сжимаемых неупорядоченных систем лишено смысла, так как они характеризуются теми же значениями эффективных зарядов Vi, Si, что и несжимаемые системы, а потому также будут

65

лежать в нефизической области. Таким образом, введение в систему замороженных дефектов структуры приводит к размытию мультикритического поведения.

Литература

1. Fisher М. Е., Ма S.-k., Nickel В. G. Critical exponents for long-range intraction // Phys. Rev. Lett. 1972. V.29. P.917.

2. Белим C.B. Влияние эффектов дальнодействия на критическое поведение трехмерных систем// Письма в ЖЭТФ. 2003. № 77. В. 2. С. 118-120.

3. Белим С.В. Влияние упругих деформаций на критическое поведение неупорядоченных систем с эффектами дальнодействия//ЖЭТФ. 2003. Т.124. Вып.6. С.507-513.

4. Белим С.В. Прудников В.В. Трикритическое поведение сжимаемых систем с замороженными дефектами структуры // ФТТ. 2001. Т.43, Вып.7. С.1299.

5. Белим С.В. Влияние эффектов дальнодействия на мультикритическое поведение однородных систем // ЖЭТФ. 2004. Т.125. Вып.1. С.152.

6. Прудников В.В., Прудников П.В., Федоренко А.А.// ФТТ. 2000. Т.42, N.l. С.158.

7. Белим С.В. Влияние стрикционных эффектов на мультикритическое поведение однородных систем // Письма в ЖЭТФ. 2002. N.75. Вып.9. С.547-550.

8. Imry Y. Tricritical Points in Compressible Magnetic Systems // Phys. Rev. Lett. B. 1974. V.33. P.1304.

9. Ларкин А.И., Пикин С.A. // ЖЭТФ. 1969. Т.56. С.1664.

10. Bergman D.J., Halperin В.I. Critical behavior of an Ising model on a cubic compressible lattice // Phys.Rev.B. 1976. V.13, N.4. P.2145.

11. Amit D. Field theory the renormalization group and critical phenomena. New York: McGraw-Hill, 1976.

12. Zinn-Justin J. Quantum field theory and critical phenomena. Oxford: Clarendon Press, 1989.