Файл: Кудрин, Л. П. Статистическая физика плазмы.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 230

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

сравним энергию двух состояний свободного электронного газа с одинаковым числом спинов, направленных «вверх» и «вниз» (немагнитное состояние), и с одинаковой ориентацией всех спинов (ферромагнитное состояние). Согласно принципу Паули, во втором случае энергия Ферми выше, чем в случае неферро­ магнитного состояния, так как на упорядочение спинов нужно затратить определенную энергию. С другой стороны, получается выигрыш в обменной энергии, поскольку в ферромагнитном со­ стоянии все электроны (а не половина их) дают вклад в обмен­ ную энергию. Численно оказывается, что если pF— импульс Ферми в немагнитном состоянии, то соответствующий импульс в ферромагнитном состоянии Рмакс= 21/3/?^ при условии, что пол­ ное число электронов в обоих случаях одинаково.

Следовательно, энергия, необходимая для упорядочения

спинов,

 

£ к и н (2Ч'Рр) £ к„н (Pf) =

(2,21/rs) -0,588 Ry.

Выигрыш в обменной энергии составляет

£ о б м (2V ) - £ о6м (Pf) = -

(0,916/0 (2V*- l) =

= - (0,916/0-0,26 Ry.

Таким образом, полная разность энергий

ОО,

Это выражение

обращается в нуль при rs= rs= 5,46, поэтому

при rs> r s Д<0,

и система электронов должна быть ферромаг­

нитной. Однако столь серьезного вывода нельзя делать на основе модели, не учитывающей кулоновскую корреляцию электронов. Дело в том, что учет кулоновского взаимодействия приводит к дальнейшему понижению энергии системы.

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

Итак, свойства квазичастиц, полученные теоретически в при­ ближении Хартри — Фока, противоречат известным эксперимен­ тальным данным. Особенно грубо эта модель описывает плот­ ность состояний на поверхности Ферми. Чтобы ликвидировать логарифмическую расходимость, содержащуюся в формуле (21.26), нужно учесть кулоновскую корреляцию и экранирование дальнодействующих кулоновских сил.

240


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

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

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

Ео

 

£

<0 I и ш 10>

 

 

(21.27)

 

 

ЕпЕа

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

Я1п1 =

V

Як а+_

, а

р,

(21.28)

Ал

о Чi о

q, о о

 

р.

q, ка,, а'

 

 

 

 

 

Какие состояния п дают вклад в сумму (21.27) и как ее вычислить? Попробуем сначала проинтерпретировать физически записанные выше выражения. Акт рассеяния, описываемый гамильтонианом (21.28), можно представить следующим обра­ зом: электроны, находящиеся в состояниях р и —q, переходят соответственно в состоянии р + q и —(q + k). Вследствие принципа Паули состояния р и —q должны находиться внутри фермиповерхности SF, а состояния р + к и —(q + k) — вне ее. Матрич­ ный элемент такого перехода, если он возможен, £/к= 4яе2/£2. Соответствующая энергия возбуждения

£

(k) =

(р + к)2

,

 

(— q — к)2

_

р2 _

jE _

_ к (р +

q +

к)

Р' 4

2т

 

 

 

 

т

 

Перепишем теперь

выражение

(21.27)

в

следующем

виде:

 

 

Е,

£

<0 I Нint 1 П> < n

I

flint I 0>

(21.29)

 

 

 

Е п

Е 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из этой записи следует, что, попав

в

возбужденное

состоя­

ние п,

характеризуемое наличием электронов в состояниях, р + к

и — (q + k)

и дырок

в состояниях

р,

—q,

необходимо еще раз

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

241


по тому же «пути», по которому система была возбуждена. Соответствующий матричный элемент равен Uк-

Схема (а) и диаграмма (б)

прямого процесса изображены

на рис.

28.

Электроны, находящиеся первоначально в состоя­

ниях р

и

—q,

в

результате

взаимодействия

друг с другом

«выбиваются»

из

сферы Ферми в состояния

р+ k и (q+ k).

Рис. 28. Иллюстрация прямого взаимодействия элект­ рон— дырка во втором порядке теории возмущений.

Этот процесс изображен сплошными стрелками на рис. 28, а. Затем вследствие повторного взаимодействия друг с другом они возвращаются обратно в свои первоначальные состояния, что показано на рисунке пунктирными стрелками. Графически этот процесс представлен на рис. 28, б. Сплошные линии соответ­ ствуют функциям распространения электронов и дырок, а пунк­ тирные линии описывают взаимодействие с передачей импульсов. Следовательно, вклад прямых процессов в Е2

q. р, к

X h2к (р

i ( l

И.р-|-к) tlq (1

^q4-k)-

(21.30)

к)

 

 

 

Ограничения, налагаемые принципом Паули на промежуточ­ ные состояния п в формуле (21.29), учтены в функциях рас­ пределения

 

п =

(0

при

| р |

>

pF;

 

 

 

 

Р

h

при

I р I

<

pF,

 

 

 

а множитель 4 в выражении (21.30)

учитывает

суммирование

по спинам.

возможный

процесс — обменного

типа. При

этом

Второй

электроны

из возбужденных состояний

р+ k

и

(q+ k)

пере­

ходят в состояния —q и р. Соответствующий матричный эле­ мент равен

— t/p+q_k = — 4ле2/(р + Q+ к)2-

242


Знак минус учитывает изменение порядка следования опера­ торов рождения и уничтожения. Отметим, что обменный процесс возможен, если спины электронов в состояниях р и q парал­ лельны, тогда как для прямого процесса такого условия нет.

Рис. 29. Иллюстрация обменного взаимодействия элект­ рон— дырка во втором порядке теории возмущений.

Типичный пример обменного процесса изображен на рис. 29. Вклад процессов обменного типа в Е2 имеет вид:

Eib) = 2

4ле2

4пе2

 

 

k2

(р f q+ к)2X

 

 

 

 

 

 

Р, q, к

 

 

 

X /г2к (р + q + к) ■Яр (1

Яр-fk ) (1

^q

(21.31)

к)-

Покажем теперь, что вклад от прямых процессов логариф­

мически расходится. Для этого рассмотрим выражение

(21.30)

в предельном случае малых передач импульса.

При этом абсо­

лютные значения |р| и |q| должны быть близки к pF, так как условия p <.Pf и |p + k| > p F ограничивают область возможных значений р (и аналогично q) тонким слоем в ближайшей окрест­ ности сферы Ферми. Толщина этого слоя порядка |к|. Ограни­ чения, обусловленные принципом Паули, в случае малых пере­ дач импульса описываются выражением

Яр (1 — Яр+к) — Яр (1

Яр кV рЯр) — Яр

(кр)

Ь(Р—РгУ>

(21.32)

 

 

 

L Р

 

(кр) > 0.

 

 

 

 

 

Здесь использовано равенство

 

 

 

 

 

VрЯр = (р/р) б (р — pF),

 

 

а также

условие

|p + k |> /jF,

которое и привело к не­

равенству

(кр) 5г0.

Следовательно,

для

оценки выражения

243


(21.30) в пределе малых k можно заменить р и q на pF и положить

^ tip (1 " tip . j,) Nkjkp

/Zq (1

Hq к)

р

q

 

Переходя в выражении (21.30) от суммирования к интегриро­ ванию, получаем

lim Eia) ~ —

т

С

dk

— In k.

h2kkp

J

k

k-+0

 

Легко видеть, что эта логарифмическая особенность непосред­ ственно связана с дальнодействующим характером кулоновских сил. Если бы взаимодействие достаточно быстро убывало с рас­ стоянием, то матричный элемент Uк при малых k стремился бы

кнекоторой константе и никакой расходимости не возникало. Эта расходимость задержала систематическое исследование

электронного газа в рамках теории возмущений примерно на 20 лет. Расходимость поправки к энергии во втором порядке теории возмущений представляет лишь одну из трудностей, воз­ никающих при вычислении энергии системы многих частиц по теории возмущений. Другая принципиальная трудность возни­ кает при рассмотрении зависимости от объема каждого из чле­ нов теории возмущений. Поскольку нас интересуют объемные свойства системы, в частности энергия основного состояния элек­ тронного газа, приходящаяся на одну частицу, то необходимо, чтобы получаемые выражения имели смысл в асимптотическом пределе, т. е. когда число частиц и объем системы достаточно велики:

N -»■ оо, V -»• оо, n — N/V — const,

где п — средняя плотность электронов в системе.

Пусть некоторое возмущенное состояние электронного газа |\F ) отличается от исходного, определенного волновой функ­ цией «вакуума» |0 ) , возбуждением нескольких пар электрон — дырка. Примером простой функции, описывающей такое воз­ мущенное состояние, может служить волновая функция

 

V I

> = I 0 > +

- “г

5 ] ф1 (Рь Р *)4 (Pi) 2-q (Р*)» (21-33)

 

 

 

 

q. Pi, Ра

где

-S

 

 

/ ч

— оператор рождения пары; cLq— оператор уничтожения

пары электрон — дырка.

скалярное произведение < ОIT ) в

 

Легко

показать,

что

случае макроскопически больших систем оказывается экспонен­

циально

малой величиной.

Так

что условно можно

записать

< 0 |XF )

~ ех р (—V). Это

легко

понять, поскольку

функция

(21.33)

отвечает возбуждению

всего

лишь двух пар

на всю

бесконечную систему. Любая волновая

функция, описывающая

244