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

книги из ГПНТБ / Кудрин, Л. П. Статистическая физика плазмы

.pdf
Скачиваний:
6
Добавлен:
22.10.2023
Размер:
17.6 Mб
Скачать

При этом многоэлектронная волновая функция в модели Зоммерфельда есть произведение одноэлектронных волновых функций

^ = П

(20.2)

i-=1 1

где Хг — пространственная и спиновая координата электрона. Поскольку электроны не взаимодействуют, то энергия всей системы

Значения вектора к представляют дискретный ряд и опре­ деляются периодическими граничными условиями:

к = (2я/1) п,

где компоненты вектора п — целые числа. Следовательно, плот­ ность узлов для одного направления спина равна

dnxdnydnz/d3k = (L/2л)3 = р (к).

Тогда переход от суммирования по к к интегрированию в к-про- страпстве для любой функции /(к) совершается по правилу*

V/(k)к = (L/2jt)»j7 (k)dk.

Какие волновые векторы к должны фигурировать в волно­ вой функции (20.2)? Согласно принципу Паули, каждому зна­ чению к можно поставить в соответствие не больше двух электронов (с противоположными спинами). При пулевой тем­ пературе (р = оо) энергия системы минимальна, если электроны заполняют в соответствии с принципом Паули N наинизших состояний. Таким образом, приходим к распределению по им­ пульсам, представляющему собой две заполненные сферы Ферми— по одной для каждого направления спина, в которых заняты все состояния вплоть до некоторого максимального им­ пульса hk р.

Величина kF определяется из условия, чтобы в занятой части пространства импульсов содержалось N/2 векторов, т. е.

N = 2 n&dk (L/2n)3= L ^ /Зл2.

о

Таким образом, максимально возможное значение kF опреде­ ляется плотностью электронов n = N(V:

kF = (Зл2п ) и .

(20.3)

* Отметим, что рассматриваемый объем системы V=L3 предполагается достаточно большим, так что спектр дискретных значений k можно считать квззинепрерывным.

220

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

Ve =i / n = (4я/3) г3 .

Удобно ввести безразмерный параметр, который будет часто использоваться в дальнейшем: гя = Го,'ао, где flo — боровский радиус. Тогда

 

 

kFr0 — (9л/4)'/:1яй 1,92,

 

а максимальная энергия (энергия Ферми)

 

е

h2k2p

_

/ 9Я V/:

те4 « 3,68r~2 Ry,

(20.4>

р

~biT

~

\ ~ Г )

2

 

где Ry = me'lj2h2. Поскольку концентрация электронов в боль­ шинстве металлов соответствует значениям параметра rs —

= 2-=-5,5, то

1,82Л-1 kF> 0,662Л~1; 12,5 зв > ер > 1,66 эв.

(20.5>

Среднюю энергию, приходящуюся на частицу, условно можно> записать в виде (она, конечно, совпадает со средней кинетиче­ ской энергией, приходящейся на частицу в газе Зоммерфельда):

««„и = 2 в* = (3/5) ef » (2,2 1/г2) Ry.

(20.6>

к

 

Отметим, что eF и еКШ1 пропорциональны п2/3. Энергия Фер­ ми' eF имеет простой термодинамический смысл. Действительно,, функция распределения Ферми по квантовым состояниям

/к (Р) = [ехр <(ек — р) Р) -г I ] - 1

(20.7)

в пределе низких температур (|}—>-оо) обращается в единицу при

всех значениях ек<Ср и в нуль

при

ек>р. Отсюда следует,

что при абсолютном

нуле

p = eJ(-. Температуру, определяемую;

соотношением T0=\/$qX ef ,

называют

т е м п е р а т у р о й в ы ­

р о ж д е н и я .

 

 

 

 

 

 

Полная энергия газа получается умножением ек на числа

состояний 2V-4nh3k2dk/(2лИ )3

и интегрированием

по всем воз­

можным импульсам:

 

 

 

 

 

 

Е =

ГFk*dk -

— (Зл2/г/з — V n u .

(20.8)

2тл2

,!

 

10 4

т

 

 

Поскольку

О

 

 

 

 

 

 

 

 

 

 

 

Е = Nix + (S/P) -

PV = р

+ р

+

Q,

 

 

 

д\1

 

ор

 

22 Е

тде S — энтропия, легко получить выражение для термодинами­ ческого потенциала вырожденного электронного газа

Q = — — (Зл2)2/а — п и .

(20.9)

5 т

Поделив это выражение па объем системы, получим выражение для давления с обратным знаком. Таким образом, уравнение

•состояния зоммерфельдовского электронного газа имеет вид

Р = (1/5) (Зл2)2/*(ПУт) п и .

(20.10)

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

При конечных температурах функция распределения Ферми плавно меняется вблизи границы Ферми ер. Говорят, что сту­ пенчатая функция распределения, соответствующая вырожден­ ному электронному газу, «размазывается». Если температура системы невелика, то размазывается лишь малая область вблизи поверхности Ферми, т. е. бkfkpta ((Зе^)-1. Физически это ■означает, что только доля электронов порядка (ре^)-1 принимает участие в тепловом возбуждении. Остальные 1—(Ре^)"1 элек­ тронов «выморожены» согласно принципу Паули. Поэтому априори ясно, что в этих условиях теплоемкость электронного таза, близкого к состоянию полного вырождения, должна быть мала, ибо малая доля электронов способна участвовать в явле­

ниях

переноса.

 

 

 

 

 

 

 

Выражения (20.9) и (20.10) являются первыми членами

разложения более общих выражений при (З^оо,

так как число

•состояний в общем

случае необходимо записывать в виде

 

 

2 у

h3k-dk

 

1

 

 

 

 

 

(2nh)3

е.хр[(ей — р)(3] + Г

 

Тогда, как нетрудно показать:

 

 

 

 

 

 

 

 

 

со

 

 

 

0

_____2 _

2 1/’ т

, / »

 

е 3^ d z

(20. 11)

 

 

 

3

л2Д3

J ехр [(ё — р) Р] + 1

 

 

 

 

*

Если речь идет об абстрактной

задаче

с электронным газом, то это

* верное

рассуждение.

В реальном

кристалле

температура,

рассматриваемая

как возмущение, должна быть меньше температуры Дебая, которая порядка сотен градусов.

222

Интеграл можно представить в виде

 

 

 

 

 

 

 

 

 

 

f (е) de

 

 

 

оо

 

 

 

 

 

 

/

 

 

 

 

 

1_ (

 

.+ (г' РД dz « 1 / (е) de - f

J

 

exp [(е — ц) Р] +

1

 

 

 

р

f)

exp г +

1

,)

 

 

 

о

 

 

оо

 

 

—М-Р

 

 

о

 

 

 

 

 

 

 

f[p +

(z/P)]-/[u-( г/P)] dz

 

 

 

 

 

 

 

, J _

(•

 

 

 

 

 

 

 

Р

Jо

 

 

ехр г +

1

 

 

 

 

 

Интегралы

вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J - ^ L V =

(1~

2У~1)Г (У) £(£/)

 

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

выражаются через Г-функцшо

и ^-функцию

Римана

(Цх) —

ОО

j

 

\ *

 

 

 

 

 

 

 

 

 

 

 

= 2

 

)

. В результате получим

 

 

 

 

 

 

/=i

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в

 

 

 

 

 

 

 

 

 

 

 

 

/

=

 

f /

(е) de f - f -

Г 2 / '

(P) +

- | £

- Г

4 Г (P)

+

• •

(2 0 .12>

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

Используя это разложение

для

/(е )= е 3/2,

получаем для

термодинамического

потенциала

с

точностью

до

двух

членов

 

 

 

 

Q =

й„ +

V'p- 2

{2^

. f

U ,

 

 

 

(20.13)

где fio — термодинамический

потенциал

полностью

вырожден­

ного ферми-газа. Рассматривая второй член как малую до­ бавку, легко получить выражение для любой термодинамической функции. Например, выражая ц через Т и V, имеем для сво­ бодной энергии

 

 

F = F0

 

/3

 

 

 

 

2

 

 

 

 

 

Поскольку S =

(<3Q/<3p) (. Р2, то из выражения

(20.13)

получим**

 

 

s = ( y J *377 v $ ~'

n 'u = Cv ■

 

 

Отметим,

что

CV — CP с точностью

до

членов, пропорциональ­

ных р-3.

Действительно, из третьего

начала

термодинамики

(теорема

Нернста) следует, что энтропия

S—>-0

при Р-*-оо

как р-"

[6].

Следовательно, разность

СРCv

обращается в

нуль как

p_2n_I. В данном случае СР—С у~р~3.

 

 

* £(3/2) =2,612; £(5/2) = 1,341; Г (3/2) = 1/2;

Г (5/2) =3/4.

 

 

** Предлагаем читателю убедиться в том, что выражение для теплоемко­ сти совпадает с выражением для энтропии S.

223

Теплоемкость удобно выразить через плотность одночастич­ ных состояний на поверхности Ферми pF:

о

( dn

N

\ f

d n \

1 1

 

mkp

P f

\Лг

j e—

A

dk J

 

 

n2ft2

Тогда

 

f

лгк

 

л2m

k

 

 

Cv

 

(20.14)

 

Ч

зр )

N

h%

' P ’

 

 

 

где Тс— постоянная Больцмана.

 

 

 

Интересно,

что измерение теплоемкости при достаточно низ­

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

точно хорошее

количественное совпадение

с формулой (20.14)

(табл. 6). Отношение m*/m равно

отношению эксперименталь-

 

 

 

 

 

Таблица 6

Сравнение экспериментальных и теоретических значений теплоемкости [8],

 

 

см. формулу (20.14)

 

 

 

 

 

!

 

 

 

Металл

Na

m*jm

J

Металл

Na

ni*Ini

C u

1

1,5

!

A1

3

1,6

 

Ag

1

1,o

 

Ga

3

0,4

Be

2

0,46

 

In

3

1,3

Mg

2

1,33

 

T1

3

1,15

Ca

2

0,8

 

La

3

4,3

Zn

2

0,9

 

Sn

4

1 ,2

Cd

2

0,75

 

P b

4

2 ,1

Mg

2

2

1

 

 

 

ного значения

Су

к теоретическому. Величина пг* называется

э ф ф е к т и в н о й

м а с с о й .

Na — число валентных

электронов

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

лов, которые иногда

называют в е щ е с т в а м и т и п а

э л е к ­

т р о н н о г о

газа .

Это означает, что ни периодический потен­

циал ионов,

ни взаимодействие электронов друг с другом

(кото­

рое не является слабым), ни взаимодействие электронов с фоно­ нами не оказывает существенного качественного влияния на электронную часть теплоемкости. Не исключено, что влияние решетки в какой-то мере компенсирует коллективные эффекты, связанные с дальнодействием кулоновских сил. Столь хорошего соответствия простейшей теории с опытом не наблюдается ни для полуметаллов типа Ш или Sb, ни для переходных метал­ лов, для которых влияние периодического потенциала играет большую роль.

.224

Если модель свободных электронов приводит к разумному результату для теплоемкости при наличии малого параметра (Рер)-1, характеризующего слабое возбуждение системы, то не

исключено, что и для

с п и н о в о й

или

п а р а м а г н и т н о й

в о с п р и и м ч и в о с т и

%s можно получить удовлетворительное

выражение даже в столь грубом приближении.

 

В присутствии магнитного поля

напряженностью Н элек­

троны с разной ориентацией спина

(по

полю и против

поля)

обладают разной энергией:

 

 

 

Ekt = (E2k2l2m) +

и*#,

(20.15)

 

 

 

 

Ek i = (h2k2/2m) \ieH,

 

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

спино­

вой поляризации электронного газа запишем «парциальные плот­

ности» электронов с двумя

возможными ориентациями спина:

Ь3

 

КрЬл

KF

n ^ f ( l + P ) ^

(20.16)

п , = ^ ( \ - Р ) = 6л2

6я2

Считая, что магнитное поле слабо поляризует «замороженные»

электроны в сфере Ферми, что имеет место при реЯ<Се^, раз­ ложим энергию основного состояния электронного газа, прихо­ дящуюся на одну частицу, в ряд по степеням отношения Р/п:

Е(Р) = Е (0) + аР2 — Рр,Я,

где

Е(0 ) — энергия

(на

одну

частицу)

основного

состояния

системы в отсутствие поля. Тогда число поляризованных спинов

 

 

 

 

Р = \хеН/2а.

 

 

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

спиновая восприимчивость единицы объема

 

 

 

 

Х* =

пр«/а.

 

(20.17)

то

Поскольку

в модели

Зоммерфельда

электроны

свободны,

а определяется

лишь

кинетической

энергией

электронов:

а=(2/3)е^. Следовательно, парамагнитная восприимчивость идеального электронного газа

3_

(20.18)

2 '

Спиновая восприимчивость электронов проводимости в щелоч­ ных металлах (Li, Na) была получена экспериментально с по­ мощью комбинации эффектов электронного спинового и ядерного магнитного резонансов [8]. Сравнение данных опыта с результатом (20.18) приведено в табл. 7.

8 Зак. 635

225

 

Т а б л и ц а 7

Можно видеть, что модель

Значение xit 10—6 CGS

Зоммерфельда приводит к не­

 

 

Эксперимен­

ожиданно

удовлетворитель­

Металл

Расчетное

ному согласию с эксперимен­

тальное

 

 

 

том. Несколько позже увидим,

Li

0,81

2,08+0,1

что это

согласие

с

экспери­

ментом

 

именно

для

теплоем­

Na

0,64

0,89+0,04

кости

и

парамагнитной вос­

даже на

основе

 

приимчивости можно

понять

феноменологической

теории

электронной

жидкости.

 

 

 

 

 

 

 

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

 

dgk

 

h-kr

 

3k

k ~ k f

nr

где fit — одноэлектроиная

энергия в периодическом поле ре-

шетки. Если

поверхность

Ферми

несферична (что имеет место

в реальных

металлах),

то

эффективная масса определяется

усреднением градиента одиоэлектронпой энергии в к-простран- стве:

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

§ 21. ПОПЫТКА УЛУЧШЕНИЯ МОДЕЛИ ВЫРОЖДЕННОГО ЭЛЕКТРОННОГО ГАЗА

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

226

выражение (20.2)], то, минимизируя энергию системы, получаем уравнение для одночастичных волновых функций фа :

 

 

1

Фа (г) = е«фв (г).

(21. 1)

р+ ос

I

r - Г ' I

 

 

 

части этого уравнения, которое

назы­

Второй член в левой

 

вается у р а в н е н и е м

 

X а р т р и, описывает среднее

поле,

действующее на данный электрон. Метод Хартри является

самосогласованным,

если первоначально выбранные функции

— собственные функции уравнения

(21.1). Если это не так, то

применяется метод

последовательных

итераций. В случае газа

свободных электронов функции фа— плоские волны, удовлетво­ ряющие уравнению

-£ -ф а ( г ) - е а фа (г).

(21.2)

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

С учетом условия нормировки Ф выражается через детер­ минант Слэтера:

 

Фх (ri)

Фа (П) •

■Фд, (гх)

Ф -

Фх (С)

Ф2 (Га) •

% (Га)

У N

 

 

 

 

 

 

 

 

Ф хМ

Ф2(» )

. % ( r N )

Тогда одноэлектронные волновые функции, согласно вариацион­ ному принципу, удовлетворяют уравнению Хартри — Фока:

+

(Г)

 

Р

— е2

Фа (Г')^> Фр (г) = еафа (г). (21.3

Это уравнение отличается от уравнения Хартри обменным членом. Задачу можно по-прежнему решать с помощью итера-

8* 227

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

/ W

ikr = - еа V f dr'e~'k r' — - Ц - eikr eik'r =

 

 

-Й—JJ

I г— г' I

 

 

 

ч *

 

 

 

= еа " V Г dr' expt ~ ‘ (k~ k ') (r-г'И eikr = _ y i

4ле* eikr

(2j 4)

Z a )

I r - Г ' i

Z A I k - k ' |2

k'

 

k'

 

 

Отметим, что, как и в последующих формулах, необходимо сум­ мирование по к', соответствующим электронам с параллельными спинами. Кроме того, сумма по к' не содержит члена с к= к'.

Из выражения (21.4) видно, что плоские волны действи­ тельно являются собственными функциями обменной части га­

мильтониана НобмЭнергия электрона в приближении Хартри—■ Фока

£& ,а

1i2k~

4ле2

hW

(21.5)

у

- - V . UaN,q'vq+k, а ■

 

jAJ I k —k' |2

2m

 

 

 

k'

q

 

Следовательно, энергия основного состояния в расчете на один электрон есть

Ех- ф

 

_____1_ 4 ^ 1 Что?2 .1 _1__ / 2,21

0 , 9 1 6 \

Ашк 2т

■к' |:i} N

Ry-

 

к

к-к'

 

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

удобно выразить через флуктуации плотности пк, являющиеся фурье-компонентами плотности п:

п( г) = Z s (r — r/) =

S nkexp(ikr).

(21.7)

i

k

 

Соотношение (21.7) написано для единичного объема системы.

Поскольку Пк = Хехр(—ikr), а фурье-компонепта

кулоновской

энергии взаимодействия

Uh = iae2jk2, то

 

Н =

S

и к( пк + пи

(21.8)

 

V 2

 

 

 

 

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

228

Впредставлении

Я= £ 6ра£с<

р, о

вторичного квантования (см. приложение)

а +

S

^k<r+.

а+ . ,а

,(21.9)

 

к p+k, a

q—к, а'

q, а

р, а’ '

'

P.q.k

где ер=р2/2т — энергия свободного электрона, а операторы рождения и уничтожения частиц удовлетворяют правилам анти­ коммутации для фермионов:

 

ар,о’ V .a 'l = [ < « :

a?,a'] = °;

 

/N

/S

"1

 

 

a t . O’

ap'. o' I = V

p' 6o, o'-

Рассмотрим

теперь

в

рамках

приближений Хартри и

Хартри— Фока

энергию

основного

состояния и элементарные

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

Если полностью пренебречь потенциальной энергией, то волновая функция основного состояния системы, которую обо­ значим |0> , будет описывать состояние с полностью заполнен­ ной сферой Ферми. По определению волновой функции основ­ ного состояния к заполненной сфере Ферми нельзя добавить

частицу, т. е. ap,a|0> = 0 при |р |>/?*■, и нельзя уничтожить электрон вне сферы Ферми, поскольку их там нет, т. е. dp,a|0> = = 0 при |p |> p F. Поэтому плотность

п

 

< 0

at o ap,G

1 при р < pF;

(21. 10)

Р, о

О при р > рр,

 

\

 

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

Е к ин = ом о

2

sDa f

a

р, a

0\ = (2,21/ $ Ry. (21. 11)

 

р, a

 

 

Р, О

 

 

 

 

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

Е г = -L / о

U

а+ , , а

, а

 

О4), (21.12)

a+,,

 

N

Р+к, a

q—к, ст'

q, а'

р , а

 

 

р, q, к, а , а '

 

 

 

 

где, как и в формуле (21.11), записано усреднение по волновым функциям «вакуума», которым в задаче многих фермионов яв­ ляется заполненная сфера Ферми.

229

Соседние файлы в папке книги из ГПНТБ