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

Учебники / Физика конденсированных сред

.pdf
Скачиваний:
72
Добавлен:
30.03.2022
Размер:
3.97 Mб
Скачать
fαβ (ω) = 2
gαβ (ω).

§ 12. Флуктуационно-диссипационная теорема

411

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

Подставляя результат (6.188) в определение функций fαβ (ω) и gαβ (ω) , получаем

 

 

 

 

 

 

)

 

*

 

 

 

−∞

 

fαβ (ω) =

i

 

 

1

 

e−β ω

 

 

 

 

dt1 e−iωt1

M α M β (t1)

0; (6.189)

 

 

 

 

 

 

 

2

 

 

 

 

 

)

 

*

 

 

 

−∞

 

gαβ (ω) =

1

 

1 + e−β ω

 

 

 

dt1 e−iωt1

M α M β (t1)

0. (6.190)

 

 

 

 

Объединяя эти два результата, получаем искомую взаимосвязь функций fαβ (ω) и gαβ (ω) :

i 1 − e−β ω

1 + e−β ω

(6.191)

Вспоминая взаимосвязь функции fαβ (ω) с мнимой частью магнитной восприимчивости (6.186), можно выразить спектральную интенсивность симметризованной корреляционной функции gαβ (ω) через мнимую часть тензора магнитной восприим-

чивости

1

+ e−β ω

 

 

 

 

 

 

 

gαβ (ω) =

 

 

 

 

 

·

χαβ (ω) − χβα(ω)

. (6.192)

2i

1

− e−β ω

Если умножить числитель и знаменатель последнего выра-

жения на eβ ω/2

и ввести обозначение

 

1

 

 

 

 

 

 

 

 

2i

 

 

χαβ (ω) − χβα(ω) = Imχαβs ,

 

которое имеет смысл мнимой части симметричной составляющей тензора магнитной восприимчивости, то выражение (6.192) можно представить в компактной форме

gαβ (ω) = · cth

β ω

· Im χαβs (ω).

(6.193)

2

412 Глава 6. Метод НСО

Выражение (6.193) и является формулировкой ФДТ Каллена – Велтона. Результаты (6.193) и (6.186) остаются в силе и в пространственно неоднородном случае, когда функции

 

 

 

 

fαβ () ,

gαβ () ,

χαβ ()

зависят от волнового вектора k и

частоты ω .

 

 

В предыдущем параграфе была введена еще одна функ-

 

 

 

 

ция Sαβ (k, ω) . Установим

ее связь с функциями fαβ (k, ω)

и χαβ (k, ω) . Для простоты рассмотрим вначале пространственно однородный случай. Тогда, используя определения

(6.171) и (6.170) и полагая k = 0 , получаем

 

 

 

 

 

 

 

 

 

 

 

 

Sαβ (ω) =

 

dt1 eiωt1

M α(t1) M β

0 =

 

 

 

)

−∞

 

*

)

 

*

 

*

 

−∞

 

 

−∞

 

)

 

=

 

dt1 eiωt1 M α M β (

 

t1)

0 =

 

dt1 e−iωt1

M

α M β (t1)

0.

При записи последнего равенства в этой формуле сделана замена переменных t1 → −t1 . С учетом этого результата выражения (6.189), (6.190) перепишем следующим образом:

 

1

 

 

β ω

 

fαβ (ω) =

i

 

1

 

e−β ω

Sαβ (ω);

(6.194)

 

 

 

 

 

 

 

 

 

 

 

gαβ (ω) =

 

1 + e

 

Sαβ (ω).

(6.195)

2

 

Поскольку соотношения (6.194), (6.195) непосредственно обобщаются на пространственно неоднородный случай, учитывая (6.193), получаем

gαβ (k, ω) =

1

1 + e−β ω Sαβ (k, ω) = · cth

β ω

· Im χαβs (k, ω).

2

2

Отсюда в пределе малых частот β ω 1 следует простое равенство

s

 

βω

 

 

 

 

Im χαβ

(k, ω) =

2

Sαβ (k, ω).

(6.196)

 

 

 

 

§ 12. Флуктуационно-диссипационная теорема

413

Наконец, подставляя в это выражение значение функции

Sαβ (k, ω) (6.179), справедливое в пределе малых k , получаем представление для мнимой части магнитной восприимчивости

s

 

 

Dk2 ω

 

 

Im χαβ

(k, ω) = χ

ω2

+ (Dk2)2

,

(6.197)

 

 

 

 

справедливое в длинноволновом приближении k 0 . В связи с полученным результатом важно пояснить, что структура мнимой части магнитной восприимчивости (6.197) «навязана» законами сохранения и свойствами симметрии рассматриваемой системы и не зависит от конкретного вида гамильтониана системы.

Задача 6.1

Доказать, что

· lim S(k, t = 0) = 1/β χ,

0 k

где χ – статическая восприимчивость в пространственно однородном случае.

Решение

Будем исходить из определения

lim S(k, t = 0) = lim dr Sp M (r)M (0)ρ0

}

e−ik

r ,

 

 

 

 

{

 

 

 

 

 

 

 

 

 

 

k→0

 

k→0

 

 

 

 

 

 

 

 

 

 

 

 

где M - намагниченность системы вдоль направления внешнего маг-

 

амплитудой h . Далее, M (0) =

 

 

 

Mi;

 

нитного поля с

 

 

 

 

 

 

 

 

 

i

 

 

lim drM (r)e−ikr = lim

dr

Miδ(r

ri)e−ikr =

Mi.

k→0

 

k→0

i

 

 

 

 

 

 

 

 

 

i

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

п

 

п

ρ

0

}

,

 

 

 

 

lim

S(k, t = 0) = Sp

M

M

 

 

 

 

 

k→0

где Mп – полный магнитный момент образца.

414

Глава 6. Метод НСО

Магнитная восприимчивость образца в феноменологической термодинамике определяется соотношением

χ = lim )M *h. h→0 ∂h

Если считать, что гамильтониан системы H во внешнем поле h можно представить в виде H = H0 − hMп , то среднюю намагниченность можно вычислить, используя усреднение по равновесному ансамблю

)M *h = Sp{Mп exp[−β(H0 − hMп]}. Sp{exp[−β(H0 − hMп)]}

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

h→0

∂h ) *

 

· )

п *0 )

 

п

*0

 

· { п

п

0

}

 

lim

M

h = β

 

M 2

M

 

2

= β

Sp M

M ρ

 

 

.

 

 

 

 

 

 

Последнее равенство в этом выражении следует из того факта, что в

отсутствие спонтанной намагниченности

M

п*0

= 0 . Таким образом,

мы доказали, что действительно

)

 

 

 

 

 

·

χ.

 

 

 

 

lim S(k, t = 0) = 1

 

 

k→0

§13. Дальние корреляции и медленные моды

В§ 11 этой главы на примере явления спиновой диффузии рассмотрены условия возникновения гидродинамических мод,

т.е. слабозатухающих в пределе k → 0 коллективных возбуждений. Показано, в частности, что если динамическая переменная удовлетворяет некоторому закону сохранения и является квазиинтегралом движения, то соответствующая автокорреляционная функция в комплексной плоскости z будет иметь гидродинамический полюс. Обобщим эти результаты, используя метод операторов проектирования Мори, развитый в § 6 – 9 главы 6.

§ 13.

Дальние корреляции и медленные моды

415

 

 

 

 

 

Определим автокорреляционную функцию операторов A(k, t)

A(k , t ) соотношением

 

 

1

CAA(k, t) = A(k, t), A+(k, 0) =

 

 

= Sp {Δ(A(k, t)ρ0τ A+(k, 0)ρ01−τ }.

(6.198)

0

 

 

 

В силу однородности пространства, как следует из формулы (5.168) в § 8 главы 5, отличными от нуля будут лишь такие

средние, для которых k = k . Будем считать оператор A само-

сопряженным. Тогда

+

 

 

 

A (k) =

A(−k) .

 

Произведем преобразование Лапласа корреляционной функ-

ции

 

 

 

 

t , определив корреляционную

CAA(k, t) по переменной

функцию

 

 

 

 

 

CAA(k, z) :

 

 

 

 

 

 

 

 

 

 

 

 

CAA(k, z) =

dt CAA(k, t) e−zt.

(6.199)

 

 

 

0

 

 

 

Как следует из § 6 настоящей главы, для функции CAA(k, z) , которая отличается от функции Θ(z) (6.106) только множи-

+ 1

телем вида ( A, A ) и дополнительной зависимостью от k , повторяя выкладки, приведшие нас от (6.94) к (6.106), можно получить представление

 

 

 

β1χAA(k)

 

 

 

 

CAA(k, z) =

 

 

 

 

 

,

 

 

 

 

 

 

 

 

z + SA(k, z) · [χAA(k)]1

 

где

 

 

 

 

 

 

SA(k, z) = β (1 P)A˙

1

 

P)A˙

+(k)

(k),

 

(1

z − (1 P)iL

(6.200)

. (6.201)

В формулах (6.200), (6.201) использовано определение

) + * 1 CAA(k, t = 0) = A(k)A (k) = β χAA(k).

416 Глава 6. Метод НСО

Кроме того, для простоты мы положили, что

˙ + + 1

iΩ = A(k), A (k) A(k), A (k) = 0.

Если величина A(k, t) была бы сохраняющейся величиной,

S

то функция A(k, z) , как это показано в § 11, была бы пропорциональна k2 и мы имели бы гидродинамический полюс.

В этом легко убедиться, если считать, что A(k) является единственным базисным оператором, удовлетворяющим уравнению

движения

 

 

 

 

 

 

 

 

 

 

˙

 

 

 

 

 

 

 

 

 

A(k) = a · A(k) + ikJA(k),

 

 

где a – некоторый с-числовой коэффициент, а

 

 

JA(k) – вектор

потока, связанный с физической величиной A . Тогда, посколь-

 

 

 

 

 

 

 

 

содержит конструкцию

ку коррелятор случайных сил SA(k, z)

˙

 

 

 

 

 

 

вклада не даст. Поэтому

(1 P)A(k) , составляющая a

· A(k)

 

˙

 

 

 

 

 

 

 

 

 

 

(1 P)A(k) k и

 

 

 

 

 

2

.

 

 

 

Если

 

 

SA(k, z) k

 

 

 

 

 

 

(k, z)

 

[χ (k)]1 = 0

 

 

lim

S

 

·

 

 

 

 

A

 

AA

 

 

 

 

 

k→0

 

 

 

 

 

 

 

 

 

 

и остается конечной величиной при k 0 , то это будет означать, что корреляционная функция CAA(0, t) будет удовлетворять уравнению

 

 

 

CAA(0, t) =

CAA(0, t)

,

(6.202)

 

∂t

τA

τ

1

= lim SA(k)[χAA(k)]1,

 

 

 

A

 

 

 

 

 

 

 

 

k→0

 

 

 

демонстрируя обычное релаксационное поведение. Действительно, если выполнить трансформу Лапласа уравнения (6.202) с использованием соотношения (6.105), получаем

1

−CAA(0, 0) + z CAA(0, z) = τA CAA(0, z),

§ 13. Дальние корреляции и медленные моды

417

или

CAA(0, 0)

CAA(0, z) = z + τA1 .

Предположим теперь, что несохраняющаяся величина имеет статические корреляции бесконечно большого радиуса и

lim

 

M0

 

 

 

 

χAA(k) RAk2

,

k→0

 

 

 

где M0 и RA – некоторые константы. Очевидно, что в этом случае

S · 1 2

A(k, z) [χAA(k)] k

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

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

Если теперь вдоль оси X приложить внешнее поле hx(r) , то под действием этого поля возникнет отличная от нуля намагниченность Mx(r) . Производя фурье-преобразование материального уравнения, получаем

) * Mx(k) = χxx(k)hx(k).

Рассмотрим в этом уравнении переход k → 0 . Очевидно, что в этом предельном случае

lim

 

*

= M

 

= χ

h ,

(k)

 

k→0)Mx

 

 

0

 

xx x

где M0 – равновесный магнитный момент образца. Поскольку поворот вектора спонтанной намагниченности не связан с

j

418 Глава 6. Метод НСО

какой-либо работой, то такой поворот будет происходить и в бесконечно малом поле. Этот факт дает основание заключить,

→ ∞

что lim 0 χxx(k) .

k

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

Таким образом, можно считать, что расходимость χxx(k) связана со статическими корреляциями большого радиуса, которые, в свою очередь, обусловлены спонтанным нарушением симметрии в основном состоянии. Обсудим более подробно, что имеют в виду, когда говорят о нарушении симметрии в основном состоянии. Рассмотрим систему спинов, взаимодействие которых описывается гамильтонианом Гейзенберга

H =

1

 

J | ri − rj |

 

 

(6.203)

2

i=j

SiSj .

 

 

 

 

 

 

Хорошо известно, что в этой системе сохраняется полный спин,

поэтому

 

 

S, H = 0.

 

 

i

Полагая, что результирующий магнитный момент направлен вдоль оси Z , вычислим среднее значение z -компоненты полного спина

i

Siz

*

0

1

Sp i

Six,

j

Sjy ρ0

=

= i

)

 

 

 

 

 

 

 

 

 

 

1

 

Sp Six, ρ0

Sjy = 0.

 

=

 

 

(6.204)

i

i

§ 13. Дальние корреляции и медленные моды

419

Из этого результата следует, что равновесный статистиче-

ский оператор

 

 

 

1

 

 

ρ

=

 

e−βH .

 

 

 

0

 

 

Z

 

Более правильно считать, что

 

 

 

ρ0 =

 

1

 

PZ e−βH ,

 

 

Z

 

 

 

 

 

 

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

том, что

ρ0 Z1 exp −β(H − Siz h)!,

i

где h – бесконечно малый параметр. В этом случае z -компонента полного спина коммутирует с ρ0 , а x - и y -компоненты не коммутируют. Внешнее поле h , каким бы малым оно ни было, придает несколько больший статистический вес состояниям, в которых суммарный спин вдоль оси Z отличен от нуля. В ферромагнитном состоянии этого вполне достаточно, чтобы выстроить все спины параллельно оси Z .

Сингулярность χxx(k) связана именно с нарушенной симметрией (магнитные моменты выстроены вдоль оси Z , хотя гамильтониан системы инвариантен относительно вращений). Именно поэтому поворот результирующего магнитного момента происходит в бесконечно слабом поле h , приложенном вдоль

оси X , что и приводит к сингулярности χxx(k) . При этом, по-

скольку χxx(k) является четной функцией k ( χxx(r) зависит только от | r |), эта особенность имеет вид 1/k2 .

Интересно отметить, что χzz (k) не имеет какой-либо особенности в своем поведении, поскольку увеличение магнитного момента вдоль оси Z при приложении бесконечно малого поля вдоль этой оси также будет бесконечно малым.

Обсуждаемый здесь результат квадратичной сингулярности статической восприимчивости в системах со спонтанно нарушенной симметрией носит название теоремы об 1/k2 – расходимости Боголюбова. В следующем параграфе на качественном

420

Глава 6. Метод НСО

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

вмонографиях [44, 49].

§14. Неравенство Боголюбова и теорема

об 1/k2 - расходимости

Приведем несколько упрощенный вывод неравенства Боголюбова и теоремы о сингулярности статических компонент обобщенной восприимчивости в системах с нарушенной симметрией основного состояния на примере статической магнитной восприимчивости. Установим сначала дисперсионные соотношения Крамерса – Кронига и правило сумм, связывающее действительную и мнимую части обобщенной восприимчивости.

Обобщенную восприимчивость определим как отклик значения оператора физической величины B на внешнее воздействие, определяемое возмущением −A+F :

)

 

*

 

 

 

δ B(k, ω) = χBA(k, ω) F (k, ω),

 

χBA(k, ω) = χBA

(k, ω) + BA

(k, ω),

(6.205)

где χBA(k, ω) и χBA(k, ω) – действительная и мнимая части тензора обобщенной восприимчивости. Предположим теперь, что

 

(6.206)

ω lim χBA(k, ω) 0.

→±∞

 

Если это не так, то следует произвести перенормировку обоб-

щенной восприимчивости, например, вычитая из χBA(k, ω) зна-

чение χBA(k, ) , чтобы предельный переход (6.206) был справедлив и для этого случая. В этом случае есть все основания

считать, что χBA(k, ω) является аналитической функцией в комплексной плоскости z . Тогда на основании теоремы Коши для аналитических функций

dz

χBA(k, z)

= 0,

(6.207)

z − ω

 

 

 

 

C