Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Статистическая физика. Часть 2. Теория конденси...doc
Скачиваний:
0
Добавлен:
01.07.2025
Размер:
3.85 Mб
Скачать

2 YnVTe3

fyaJV, Т

3/2

(85,23)

(А. А. Веденов, 1959). Согласно общей теореме о малых добав­ках, эта же формула, выраженная через другие термодинами­ческие переменные, дает поправку к другим термодинамическим потенциалам.

Для невырожденной плазмы все производные дпа/дца = па/Т, и тогда (85,23) переходит в формулу

3/2]

КОр "

(85,24)

для поправки к свободной энергии, совпадающую с V (78,12).

В случае сильного вырождения электронов в плазме <^ре) производная дпе/дце ~ ne/[ie<^.ne/T. В сумме по а в (85,23) можно тогда вообще пренебречь электронным членом, и мы снова по­лучаем формулу (85,24) с той лишь разницей, что сумма в ней берется лишь по сортам ионов в плазме. Таким образом, при сильном вырождении электроны вообще не влияют на радиус экранирования и на корреляционную часть термодинамических величин плазмы.

Глава IX

ГИДРОДИНАМИЧЕСКИЕ ФЛУКТУАЦИИ

§ 86. Динамический формфактор жидкости

Рассмотренная в V § 116 корреляционная функция флукту­ации плотности является частным случаем более общей функции, связывающей флуктуации плотности не только в различных точках пространства, но и в различные моменты времени. В классической теории эта функция определяется как среднее значение

na(t; l-j, г2) = <8я(г£, ri)8n(tif r2)>, (86,1)

где t = tt t2\ из определения а вынесен множитель n = NlV средняя плотность числа частиц. Для однородной и изотропной среды (жидкость, газ) функция (86,1) зависит от rj и г2 только через расстояние r = \rt — г2| между двумя точками, что и будет предполагаться ниже.

В квантовой теории аналогичная функция определяется с помощью симметризованного произведения зависящих от вре­мени (гейзенберговских) операторов плотности как

no(t, r) = ^<8n(ti, rj)6n(i2, r2) + 8n(t2, г2)8я(^, rj) (86,2)

(в соответствии с общим способом определения согласно V (118,4)). Некоторые преимущества, однако, имеет в данном случае несимметричное определение

па (г, r) = <8n(tit г^ЬпЦ,, г2)>, (86,3)

для которого сохраним обозначение а (/, г)1). В противополож­ность функции a(t, г) функция o(t, г) не является четной по переменной t; очевидно, что

o(t,r) = ^[a(t,r)+a(-t,r)]. (86,4)

J) Именно эта функция является непосредственно наблюдаемой величиной, например, при неупругом рассеянии нейтронов в жидкости —см. задачу.

Фурье-образ функции a(t, г) по времени и координатам

СО

о-(ю,к)= а(©, k) = l S о (t, r)dtd3x (86,5)

00

называют динамическим формфактором среды. Ввиду изотропии функции a(t, г) он зависит только от абсолютной величины волнового вектора. Из (86,4) следует, что фурье-образ функции a(t, г)

о (со, k) = 1 [а (со, k) + а (- со, *)]. (86,6)

Чисто пространственная корреляция флуктуации плотности жидкости определяется функцией (86,1) при £ = 0: о (г) = = a (t = 0, г) = о (t = 0, г). Эта функция связана с введенной в V § 116 (и использованной в § 83) функцией v(r) согласно а (г) = v (г) + 8 (г); их фурье-образы: а (/г) = v (k) -f- 1. Функцию а (k) или v(&) называют статическим формфактором жидкости. Функции а (со, k) и a (k) связаны друг с другом интегральным соотношением

СО 00

o(fe)= j а (со, fe)e-to*g|<=o= j а (со, fe)g. (86,7)

со со

Шредингеровский (не зависящий от времени) оператор плот­ности дается суммой

я(г) = (г-гв), (86,8)

а

взятой по всем частицам среды; координаты частиц га играют роль параметров (ср. (24,4)). Нам понадобятся ниже компоненты фурье-разложения этого оператора по координатам

пк = J п (г) e-(krd»x=2e"(kr«. (86,9)

а

Переход к зависящему от времени (гейзенберговскому) опера­тору происходит по общему правилу

п (t, г) = ехр (iHt/fbfn (г) ехр (— iHt/A), (86,10)

где Н—гамильтониан системы. Этот оператор может быть пред­ставлен выражениями (86,8—9) с заменой в них га на ra (t) гейзенберговские операторы координат частиц.

Согласно основным принципам статистики, усреднение <. ..> можно понимать по-разному, в зависимости от того, через какие термодинамические переменные должен быть выражен результат.

TaKj если функция ст определяется при заданных полной энергии и числе частиц системы, то усреднение производится, по опре­деленному (m-му) стационарному состоянию, т. е. взятием соответствующего диагонального матричного элемента. Для однородной системы (жидкость) зависимость матричных элемен­тов оператора bn(t, г) от времени и координат дается формулой

<m\bn(t, r)|Z> = <m|6u(0)|/>exp[f(<de,f-k»,r)], (86,11)

вполне аналогичной (8,4) (в правой части стоит матричный элемент шредингеровского оператора бп(г), взятого в точке г = 0). С учетом этой формулы пишем

na(t, г) = 2 | бп (tit rf) | /></ ] бп (t2, г2) | т> = = 2 \<т | бп (0) | /> |« ехр [I (co^-k^r)].

Фурье-образ этой функции по (со, к) = (2я)< 2 |<m I (0) I />|2 б (со—со^) б (k—k,J. (86,12)

Суммирование в этих формулах производится по всем состояниям системы с заданным (Nm) числом частиц (поскольку оператор б/г не меняет этого числа).

Если же мы хотим выразить формфактор через температуру и химический потенциал жидкости, то выражение (86,12) должно еще быть усреднено по распределению Гиббса:

по (со, k)

= (2я)« £ ехр Га-Е»-Р"»\{<т | бп (0) | />|' б (со-со, J6 (к-кш)

I, т

(86,13)

(причем во всех членах суммы N[ = Nm). Выписав такую же формулу для а(—со, —к) = о(—со, к), взаимно переобозначив в ней индексы суммирования / и т и заменив в экспоненци­альном множителе Е1 Ет + Ет-\-1ш (последнее равен­ство—следствие наличия б-функции), получим

а(— со, £) = а(со, к)е~^Т .(86,14)

и затем, согласно (86,6),

а(со, k) = \{\+e-^lT)o{«>, k). (86,15)

Отметим, что из (86,13) (или (86,12)) следует, что функция о (со, fe)^0 при всех значениях ее аргументов. Из соотношения

же (86,14) следует, что при нулевой температуре

a(co,fe) = 0 при со<0, 7 = 0. (86,16)

В макроскопическом пределе (N,V—^оо при заданном отно­шении N(V) «частокол» б-функций в (86,13) размазывается в непрерывную функцию, но б-функционные пики в сг (со, k) остаются при значениях со=со(&), отвечающих незатухающим элементарным возбуждениям (как это следует из рассуждений, подобных изложенным в § 8). Такие пики возникают, однако, лишь для возбуждений без изменения числа частиц1).

Покажем, каким образом формфактор жидкости может быть связан с величинами, фигурирующими в общей формулировке флуктуационно-диссипационной теоремы (D. Pines, Ph. Nozie-res, 1958).

Пусть на каждую частицу жидкости действует некоторое внешнее поле, сообщающее частице потенциальную энертию U (/, г). Тогда оператор возмущения, действующий на жидкость в целом, будет

9{t) = lh(t,r)U(t,t)d»x. (86,17)

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

8Й(со, tt) = $сс(со,|г—tt\)U(co, r2)d*x2, (86,18)

где функция а (со, г) играет роль обобщенной восприимчивости. Фурье-компонента по времени от корреляционной функции o(t, г) есть, в обозначениях флуктуационно-диссипационной теоремы:

по (со, г) = (бп(г1)вя(г|))в, r = ri—г,.

Согласно этой теореме, эта функция выражается через обобщен­ную восприимчивость формулой

по (со, г) = flcth Ima(co, г). (86,19)

Такой же формулой выражается фурье-компонента по коорди­натам а (со, k) через а (со, k), после чего, согласно (86,15),

г) Так, в ферми-жидкости а (со, k) имеет б-функционную особенность при со= ku00скорость нулевого звука), ноне имеет таких особенностей, отве­чающих фермионной ветви спектра — см. § 91.

находим для динамического формфактора

по (со, k) = Ima(co,fe). (86,20)

1 _е-лю

Важность этих формул связана прежде всего с тем, что ими устанавливается связь динамического формфактора с функцией с известными общими аналитическими свойствами (по перемен­ной со); для функции а (со, k) эти свойства описаны в V § 123. Они позволяют также применить к вычислению формфактора общую формулу (ср. (75,11)), согласно которой

а (со, k) =

со

= тИе''"кГ)<"^' ^ "(°' 0) °)"(*» r)>dtd3x. (86,21) о

Выразив операторы плотности через гр-операторы (n = W+W), можно привести это выражение к виду двухчастичной функции Грина, для вычисления которой применима диаграммная тех­ника.

Задача

Выразить через динамический формфактор вероятность неупругого рас­сеяния медленных нейтронов в жидкости, состоящей из одинаковых атомов (G. Placzek, 1952).

Решение. Согласно методу псевдопотенциала (см. III § 1.51), рассея­ние медленных нейтронов может быть описано как результат взаимодействия с потенциальной энергией

0)

М

2п%!

- an (г),

где п (г) — оператор плотности (86,8); М — приведенная масса атома и нейтрона; а—длина рассеяния медленного нейтрона на отдельном атоме (т. е. взятое с обратным знаком предельное значение амплитуды рассеяния). Вероятность перехода из некоторого начального (i) состояния системы жидкость + нейтрон в конечное (/) состояние в некотором интервале d\f есть

СО

(2)

(см. III (40,5)); для недиагональных матричных элементов U ц в (1) можно писать Ьп вместо п. Волновую функцию начального состояния нейтрона (с импульсом р и энергией е) нормируем на одну частицу в объеме V, а вол­новую функцию конечного состояния (импульс р' и энергия е') нормируем на б-функцию от р/2и. Тогда d\f = d*p'/(2rih)3, а матричный элемент возму­щениягде fik = p— р*, Асо = е— г', a bnp(t, г)—матричный элемент по отношению к волновым функциям жидкости. Подставив это выражение в dwj;, просумми­руем вероятность перехода по всем возможным конечным состояниям жидко­сти. При этом квадрат модуля интеграла записываем в виде двойного интег­рала (по dt dt' d3x d3x') и замечаем, что

2 bnfi (t, г) brift (tr, г')* = 2 О б"/« С г) =

f f

= <i\6n(t\ r')bn{t, t)\iy = no(t' — t, r'—t)

(причем о выражено в функции от полной энергии жидкости в состоянии i). Интегрирование по d(t' t)ds(x'—х) дает а (со, k), а еще одно интегрирова­ние (скажем, по dt d3x) дает просто объем V и полный интервал времени /. Опустив множитель t, получим в результате вероятность рассеяния в единицу времени

4яФ-2 , ., d3p"

w = по?а (со, k) j— i (3)

Мг (2я£)3

Это выражение остается, конечно, справедливым и после усреднения по рас­пределению Гиббса, т. е. при формфакторе, выраженном через температуру.

Отметим, что свойство формфактора (86,16) в применении к рассеянию нейтронов выражает собой тот факт, что при Т = 0 жидкость может только приобретать, но не отдавать энергию. Соотношение же (86,14) выражает собой принцип детального равновесия, так как процессы рассеяния с передачей энергии и импульса (со, к) и (—со, — к) являются взаимно обратными.