Учебники / Физика конденсированных сред
.pdf
§ 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αβ (kω) , |
gαβ (kω) , |
χαβ (kω) |
зависят от волнового вектора 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.
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 + τA−1 .
Предположим теперь, что несохраняющаяся величина имеет статические корреляции бесконечно большого радиуса и
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 – равновесный магнитный момент образца. Поскольку поворот вектора спонтанной намагниченности не связан с
§ 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 – расходимости Боголюбова. В следующем параграфе на качественном
